Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM3
|
|
- Ginger Craig
- 5 years ago
- Views:
Transcription
1 Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol., pp (11 ARTICLE Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM3 Akio IKEBATA 1,, Yoshihiko MURAOKA 1 and Feng XIAO 1 TOTO LTD., -1-1, Nakashima, Kokurakita-ku, Kitakyushu 8-861, Japan Tokyo Institute of Technology, 459 Nagatsuta, Midori-ku, Yokohama, 6-85, Japan This talk presents the implementation of simulations on multiphase fluid dynamics on hardware of multiple GPGPU architecture by using robust and efficient numerical methods. An unsplit formulation for the advection computation is proposed to take the place of the original split formulation in the so-called VSIAM3 method. The new formulation improves dimensional symmetry of numerical results without losing the advantages of the original formulation. We also devised a numerical scheme named STAA to capture the moving interface between different fluids. The volume of fluid function is transported by the above advection scheme, and then modified by an anti-diffusion step based on a well designed numerical flux that rigorously conserves the mass and effectively eliminate the numerical smearing around the interface. Our experiments show that it is able to maintain the VOF function in a well-regulated shape with compact transition layer between different fluids. We have developed a 3D numerical model for multi-fluid simulations with above methods incorporated, and implemented it on a multiple GPGPU computer. It is shown that the present model is able to adequately make use of the advantage of the GPGPU architecture. Through this research, we have built a numerical tool with adequate reliability for sanitary ware product design. Computational efficiency which is of particular importance for real-case design process is achieved by implementing the code on a multiple GPGPU hardware. KEYWORDS: multiphase flow, numerical method, finite volume method, CIP, VSIAM3, GPGPU, OpenCL I. Introduction Tremendous efforts for developing numerical simulation technology for gas and liquid multiphase fluid have been dedicated for a long time. In the so-called one-fluid model, the interfaces between different fluids are explicitly solved by the interface capturing (or tracking schemes. Among the representative ones of this category, VOF (volume of fluid and level set methods have been suggested to identify different phases. 1 4 These methods transport the identification functions through the advection equations based on the velocity field computed from the Navier-Stokes equations. In practical applications, interface capturing schemes of both numerical accuracy and computational efficiency are of much interests. Following the recent trend in the expanding use of the new type hardware architectures, such as the GPGPU, numerical schemes that are well suited for multi-thread processing are particularly demanded. Level set method is suitable for GPGPU because the calculation of the scheme can be constructed as the assembly of simple vector processing. However, the mass conservation is not inherently guaranteed in the level set method. As a consequence, small liquid drops and bubbles tend to be lost. At this point, the VOF method that exactly conserves the fluid volume is more attractive. However, the geometrical reconstruction procedure in the conventional VOF method has algorithmic complication which prevents it Corresponding author, akio.ikebata@jp.toto.com from being well oriented for the GPU processing. In this work, we propose a new advection scheme using multi-moment concept which is very suitable for GPU computing. Different from the conventional VOF method, an interface capturing scheme that doesn t require geometrical reconstruction has been also suggested. An improved VSIAM3 model for multiphase simulation has been developed and implemented on a multiple GPU computer system. An integrated simulation tool which is able to use the advantage of GPU acceleration has been established for practical applications in product design process. II. The Advection Transport Scheme We consider the following two-dimensional advection equation ϕ t + u ϕ x + v ϕ =. (1 It can be alternatively written as ϕ t + uϕ x + vϕ ( u = ϕ x + v. ( For incompressible flow, it is obvious that the conservation of scalar value ϕ is guaranteed if Eq. ( is solved by the finite volume method (FVM. The CIP-CSL (Constrained Interpolation Profile- Conservative Semi-Lagrangian schemes are a class of c 11 Atomic Energy Society of Japan, All Rights Reserved. 491
2 49 Akio IKEBATA et al. conservative advection schemes of high numerical accuracy and computational efficiency. 5 9 The most remarkable feature of a CIP-CSL scheme which distinguishes it from the existing methods is using multi-moments and treating them simultaneously as the computational variables. When implementing the CIP-CSL schemes in multi-dimensions, dimensional splitting is usually used to get around the complexities of integrating the piece-wisely constructed interpolation functions over distorted volume. An economical splitting is introduced in the study to simplify the algorithm that uses all moments of three dimensions. 7 9 The simplified scheme makes use of only two types of moments, i.e. the volume integrated average (VIA and surface integrated average (SIA for discretizations in three dimensions. The resultant numerical formulations are called VSIAM3 (Volume/Surface Integrated Average based Multi-Moment Method. In two dimensional context, we define two kinds of moments, VIA and SIA as follows, and V ϕ 1 x y xi+ 1 yj+ 1 x i 1 y j 1 ϕ(x, ydxdy (3 yj+ 1 S x ϕ i+ 1 j ϕ(x i+ 1, ydy, (4 y j 1 S y ϕ + 1 xi+ 1 ϕ(x, y j+ 1 dy. (5 x i 1 Instead of the splitting in the existing VSIAM3, we propose here an alternative to update the SIA moment and numerical fluxes in an unsplit manner. In the unsplit method, SIA is calculated by semi- Lagrangian scheme based on Eq. (1, while VIA is updated with Eq. ( by a finite volume formulation over the control volume (mesh element. As shown in Fig. 1 we approximate the numerical fluxes by the Gaussian quadrature. The flux on right face of the control volume, for example, is calculated by F xi+ 1 j =u t { S x ϕ n (x i+ 1 g 1u t, y j g 1 v t + Sx ϕ n (x i+ 1 g u t, y j g v t }, where g 1 (3 3/6 and g (3 + 3/6. In the same way, the flux on bottom face of the control volume can be calculated by F y+ 1 = v t { S y ϕ n (x i g 1 u t, y j+ 1 g 1v t + Sy ϕ n (x i g u t, y j+ 1 g v t The VIA is then updated by V ϕ n+1 ( = V ϕ n F Fxi+ 1 j x i 1 j x ( + t V ϕ n u x + v. }. + F y+ 1 F y 1 y (6 (7 (8 + 1 V φ S x φ i 1 j S x φ i+ 1 Fig. 1 1 (u t, v t j (u t, v t X axis direction flux On the other hand, SIA is updated by S x ϕ n+1 i+ 1 j = Sx ϕ n (x i+ 1 u t, y j v t (9 S y ϕ n = Sy ϕ n (x i u t, y j+ 1 v t. (1 The key issue to compute Eqs. (6, (7, (9 and (1 is to properly update the numerical solutions for SIA variables. Without losing the generality, we assume a locally frozen velocity and denote the effective velocity components by u and v. The solutions of Sx ϕ n (x i+ 1 u t, y j v t and S y ϕ n (x i u t, y j+ 1 v t can be computed by the following two-step procedure as shown in Fig. Step 1 In x direction, the intermediate values for SIA S x ϕ (x i+ 1 u t, y j and flux Fx(x i+ 1 u t, y j are calculated by one dimensional CIP-CSL scheme. In the same manner, those in y direction, S y ϕ (x i, y j+ 1 v t and F y (x i, y j+ 1 v t, are calculated. Step From the intermediate results obtained at step 1, we update the SIA value by, S x ϕ n (x i+ 1 u t, y j v t = Sx ϕ (x i+ 1 u t, y j t = Sx ϕ (x i+ 1 u t, y j t ( v ϕ { v ϕ ( } v ϕ (11 The flux gradient in y direction on the right hand side of (11 can be approximated by a linear interpolation in terms of the intermediate flux function at step 1 as, v ϕ = F + F y with ( 1 F + + ξ Fy (x i, y j+ 1 v t ( 1 + ξ Fy (x i+1, y j+ 1 v t (1 (13 PROGRESS IN NUCLEAR SCIENCE AND TECHNOLOGY
3 Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM S x φ V φ i 1j S x φ i+ 1 j + 1 S x φ V φ i 1j S x φ i+ 1 j (u t,v t 1 (u t,v t Fig. Two step SIA calculation method, STEP1(left and STEP(right. ( 1 F + ξ ( 1 + ξ Fy (x i, y j 1 v t F y (x i+1, y j 1 v t ( Initial state MUSCL (κ = where ξ = u t/ x. The velocity gradient in y direction on the right hand side of Eq. (11 is approximated by central difference method. In the same manner, SIA Sy ϕ n (x i u t, y j+ 1 v t can be updated Original VSIAM UTI-VSIAM3 By repeating above procedure, fluxes can be calculated by Eqs. (6 and (7 with the effective velocity (u = g 1, u and v = g 1, v, and thus the VIA is updated by Eq. (8. The SIA is updated by letting u = u and v = v. This scheme can be straightforwardly extended to three dimensions. To distinguish the present method from the existing VSIAM3 algorithm, we call this scheme "UTI-VSIAM3" (Unsplit Time Integration VSIAM3. Because of unsplit time integration, the dimensional symmetry is much improved by the present scheme. We use UTI-VSIAM3 as the transport scheme for both moving interface between different fluids and the fluid dynamic solver. III. Numerical Tests We tested the present scheme by solving D advection problem. The advection velocity u and v are specified to give a solid rotation around the center of computational domain with a mesh of We specified the initial advected profile being centered over the computational domain and having a shape of letter "T" (for TOTO as the name of our company with discontinuous jumps between and 1. The result of UTI-VSIAM3 under cfl=.3 (775 cycle per rotation is shown in Fig. 3. It is obvious that the present scheme is more accurate than the conventional MUSCL scheme, and is comparable to the original VSIAM3. Moreover, without the directional sweep, the present scheme guarantees the dimensional symmetry, which provides a great convenience in optimizing the partitioning for parallel and GPU processing. It should be noted that the UTI-VSIAM3 also suffers from numerical diffusion as other Eulerian type schemes. Thus, extra numerical manipulation will be needed to reduce the smearing of the transition layer of the moving interface. In Fig. 3 Results of D advection test with various schemes (after 1 rotation. the present model, we make use of an efficient anti-diffusion method which is to be presented in the next section. 11 IV. Correction of Numerical Diffusion After the calculation of F / t + (u F =, where F denotes VOF function, we make modifications to reduce the numerical diffusion at interface between liquid (F = 1 and gas (F =. Our approach consists of two steps. At first step, we construct a level set function ϕ, that is a signed distance function from the interface based on the VOF function F. As the second step, the cells where ϕ > ϵ are corrected to F = or F = 1 by an effective anti-diffusion formulation of flux form so that the conservation of F is exactly guaranteed. In level set function calculation, initial value of ϕ must be set from F. We use the following expression, ϕ = α (F.5. (15 The level set function ϕ is then computed iteratively by ϕ n+1(x, y = ϕ n (x Cδ x, y Cδ y (16 ϕ n+1 (x, y = ϕ n+1(x, y + Csgn(ϕ δ (17 δ sgn(ϕ ϕ (18 ϕ where C CFL mesh spacing. If ϕ n+1(x, y ϕ n (x, y <, CFL must be reduced locally until ϕ n+1(x, y ϕ n (x, y is met. ϕ n+1(x, y can be interpolated by bi-linear function. VOL., OCTOBER 11
4 494 Akio IKEBATA et al Fig. 4 Result of D advection test with UTI-VSIAM3 + STAA (after 1 rotation In the anti-diffusion flux calculation, we evaluate the amount of VOF to be redistributed by F F F where F = for ϕ< and F =1 for ϕ>. F is then transported toward the interface along the opposite direction of the signed normal vector given by Eq. (18. In practice, we distribute F i,j to its neighboring cells (iup, j, (i, jup, (iup, jup according to the orientation of the interface, i.e. F i,j is tranported to (iup, j, (iup, jup if δ x > δ y and to (i, jup, (iup, jup if δ x < δ y. The corresponding ratio for distributing F i,j between (iup, j and (iup, jup is calulated by 1 δ y / δ x in the former case. In the latter case, the ratio between (i, jup and (iup, jup is calulated by 1 δ x / δ y. The result of D advection test is shown in Fig. 4. The sharpness of interface between F = 1 and F = are well maintained while the transported profile of shape "T" is faithfully reproduced. As shown in the latter numerical tests, the present scheme works well as an interface capturing scheme that has adequate numerical accuracy and robustness. This scheme is called STAA (Surface Tracking by Artificial Antidiffusion scheme. As discussed above, because both UTI-VSIAM3 and STAA schemes exactly guarantee numerical conservation, the VOF function is conserved. Meanwhile the level set function ϕ which is generated from the VOF function F provides other conveniences to get the geometric quantities of the interface, like normal vector and curvature, needed in surface tension calculation. Different from the conventional VOF method, there is no geometrical reconstruction involved in the STAA method, so it is computational efficient and can be extended to 3D straightforwardly. V. Navier Stokes Solver with Multiple GPGPU For solving liquid and gas multiphase flows, the following Navier Stokes equation and VOF advection equation must be computed. u t + (u u = 1 ρ P + 1 ρ : T + g + F (19 F t (u F = ( We will calculate above equations with CIP-CUP method. 1 CIP-CUP method is based on a fractional step solution procedure. First, the left hand side of the above equations, i.e. the advection part, is solved by UTI-VSIAM3, and the resulted F is then corrected by STAA. Second, the nonadvection terms except for pressure term is calculated by a method equivalent to the standard central difference. F includes surface tension force. It is calculated by the CSF (Continuum Surface Force formulation. 14 In this formulation, the curvature can be easily calculated by central difference of a level set function ϕ calculated in the STAA scheme. Finally, a pressure projection step is calculated after the following Poisson equation for pressure is solved ( 1 ρ P = 1 t ( 1 ρc s P t + u (1 where u denotes the intermediate velocity calculated from the previous sub-steps. We have implemented the code on a multiple GPGPU system in order to save computational time in practical applications. The multiple GPGPU system used in the present study is constructed by installing multiple GPU boards (Geforce GTX 95 GPUs into a host PC (Intel Xeon 3. GHz CPU. To retain the numerical accuracy of the original code, which appears to be essential in real applications, we use double-precision data even in GPU hardware environment. It should also be notified that to make the code easy to be read and accessed by users, as well as to be maintained and continuously developed by a community, sophisticated techniques that lead to the full use of the GPU advantage, such as the optimization of the local memory is got around in this study. All these consideration result in a strategy which may compromise to some extent in computational efficiency. Each GPGPU calculates one of the partitions divided by a domain decomposition technique. At each step, data exchange is conducted among GPUs which shares the partition boundaries of neighboring domains. Based on the current state of data transfer in the GPU architecture, the data to be exchanged among GPUs are firstly sent to the CPU memory buffers, and then transfered to the targeted GPUs. In these data transfer procedures, bottle neck due to the limited bandwidth of data pass between CPU memory and GPU memory may cause degradation of the scalability performance against the number of GPUs. In GPGPU computations, each thread of GPU calculates its assigned element locally. The numerical algorithms used in the present model is overall suited for GPU computation. In advection phase of u and F, UTI-VSIAM3 can be implemented by assembly of simple vector calculations because each element value of vector array can be calculated independently. In the STAA method, level set function reconstruction can be implemented in the same way. In anti-diffusion flux calculation, however, the redistribution of F to its neighborhood may cause memory access conflicts. We solve this problem by decomposing the loop into = 7 sub-loops in which the redistribution of F at the cells having same number as shown in Fig. 5 can be calculated independently by GPU threads. Other parts except for the pressure-projection are all based on explicit algorithms, and thus can also be implemented by PROGRESS IN NUCLEAR SCIENCE AND TECHNOLOGY
5 Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM Table 1 Comparison of PCG elapse time between CPU and GPU with OpenCL (cells=1,117,866, cycles=677 Y X Devices Elapse Mem. BW Ideal BW Xeon 1 Core.3 s 9.4 GB/s 3. GB/s W558 Core 1.4 s 15.5 GB/s 3. GB/s 3. GHz 3 Core 9.7 s 19.7 GB/s 3. GB/s 4 Core 9. s 1.3 GB/s 3. GB/s 1 GPU.59 s 74.1 GB/s GB/s Geforce GPU 1.48 s 19.6 GB/s 3.8 GB/s GTX 95 3 GPU 1.1 s 17.8 GB/s GB/s 4 GPU.95 s 1.7 GB/s GB/s Fig. 5 Calculation order of anti-diffusion flux on STAA simple vector calculations. The solution of Poisson equation of pressure, on the other hand, requires extra manipulations. Among the popular solvers for Poisson equation, preconditioned ICCG or Bi-CGSTAB are widely used due to their robustness and efficiency. However, recursive dependencies are usually involved in the preconditioner calculation. In the present study, we used PCG solver with the diagonal scaling preconditioner, and found that owing to the reduction of total floating point operations per cycle its computational time is comparable to that of the ICCG. The number of cycles in our test with ICCG was roughly 1/4 of PCG while ICCG calculation time per cycle using single CPU core was two times of PCG because of the heavy precondition calculation. Since the PCG is more parallelism oriented, we have not explored further the implementation of the ICCG on the GPU. We implemented present code with OpenCL in order to retain compatibility among various hardwares. It is known that a code using OpenCL appears to be more complicated compared to those using CUDA Runtime API. We simplified our code by using macro definition or our own basic API package. For example, although each arguments of calling kernel function must be set with "clsetkernelarg" one by one and called with "clenqueuendrangekernel", we wrapped these instructions with our API "CL_Exec", which was coded using variable argument system calls of C, "va_arg", "va_start", "va_end" etc. In the same way, we defined "CL_Put" and "CL_Get" with only four arguments which substitute the OpenCL APIs, "clenqueuewritebuffer" and "clenqueuereadbuffer" that have nine arguments. We measured the elapse time of PCG calculation as shown in Table 1. Generally, solution of the Poisson equation consumes the major part of computational time of the whole code. The calculation speed of 4 GPUs was almost 1 times faster than that of 4 Core CPUs. The scalability performance of 4 GPUs is almost /3. The reason of the performance degradation is partly due to the communication overhead among GPUs and the vector length. Because we do not use the "overlapped communication", the performance degradation due to communication seems not negligible. Nevertheless, the present program uses multi-threads parallel computation, so the major part of total communication time among GPUs is host-device communication time. The device-to-host and host-to-device communication on the PCIExpress(x16 band- width in our tests is about GB/s although we use the "pinned memory" objects. Nevertheless, we still consider it a big t= Fig. 6 t= t= t=1. Result of dam break test (X-Z plane at Y=75 mm gain of practical significance in terms of cost-performance improvement since a 1 time leap of the GPGPU code largely shortens the elapse time for each real-case application which usually involves number of runs in product design and optimization. We computed the dam break test problem and show the result in Fig. 6. The number of cells is with x = y = z = 3. mm. The gravity in the opposite direction of Z axis is 9.8 m/s. At t =. second (initial condition, F of the region on. x 15. mm and. z 15. mm is set to 1. to identify the water region. The densities and viscosities of water and air are specified cor- VOL., OCTOBER 11
6 496 Akio IKEBATA et al. Fig. 7 Result of simulation of a sanitary ware (mesh number = , CFL=.48 respondingly. It is obvious that the transition layer between water and air has been kept with a compact thickness in the simulation, and small droplets and bubbles are well resolved. Moreover, we would like to emphasize that the calculation was quite stable from the beginning to the end. VI. Application as a Design Tool for Products In our company, various products like sanitary wares, kitchens and bath rooms among many others, are produced. Because many products involve water flows, simulations of water/gas/body (deformable or un-deformable multiphase flows are important for design and improvement of our products. A few examples of applications in our product designs are shown in Figs. 7 and 8 where the free surface and moving solid bodies are reasonably reproduced. Our codes can simulate real-case problems which involve multiphases including arbitrarily deformable or undeformable bodies. The calculation of moving bodies is also accelerated with GPGPU computing. The moving bodies shown in Fig. 7 are calculated by method of a finite difference calculation procedure in which bodies represented by color functions are accelerated by integral of forces in each 1, 13 cells. In the simulation shown in Fig. 8, the bucket motion is specified by users. The motion of the solid bodies are calculated by coordinate transformation. As a matter of fact, the elapse times of these simulations are almost less than one or two seconds per time step with only one GPGPU system described above. Our simulation codes are very practical in respect of both the reliability of numerical outputs and the computational efficiency. With the results of multi-fluid simulations for product design, we are able to evaluate and optimize the performances of those products before manufacturing. Using this powerful tool, we can test different designs visually, which effectively reduces the cost and speeds up the product design procedure. VII. Conclusion Some recent efforts are made toward the establishment of an integrated simulation tool for designing products that involve multiphase flows. We have proposed a more accurate Fig. 8 Result of simulation of a kitchen (mesh number = , CFL=.48 and efficient formulation for the advection computation by using the unsplit solution procedure, so-called UTI-VSIAM3, to improve the original VSIAM3. In order to compute the moving interface in multiphase flows, a numerical correction method named STAA has been incorporated with UTI- VSIAM3, which results in a conservative and accurate interface capturing scheme well suited for GPU architecture. Owing to the simplicity of the proposed schemes, we were able to implement the multiphase flow solver on a multiple GPGPU hardware with practical computing performance and scalability. We expect the code will be used to the various applications in real-case design processes. References 1 D. L. Youngs, "Time-dependent multi material flow with large fluid distortion," Numerical Methods for Fluid Dynamics, K. W. Morton, M. J. Baines (eds., Academic Press, New York, 4, 73 (198. W. J. Rider, D. B. Kothe, "Reconstructing volume tracking," J. Comput. Phys., 141, 11 ( J. A. Sethian, "Level Set Methods and Fast Marching Methods," Cambridge University Press ( M. Sussman, P. Smereka, S. Osher, "A level set approach for computing solutions to incompressible two-phase flow," J. Comput. Phys., 114, 146 ( T. Yabe, R. Tanaka, T. Nakamura, F. Xiao, "Exactly conservative semi-lagrangian scheme (CIP-CSL in one dimension," Mon. Weather Rev., 19, 33 (1. 6 F. Xiao, T. Yabe, "Completely conservative and oscillationless semi-lagrangian schemes for advection transportation," J. Comput. Phys., 17, 498 (1. 7 F. Xiao, T. Yabe, X. Peng, H. Kobayashi, "Conservative and oscillation-less atmospheric transport schemes based on rational functions," J. Geophys. Res., 17[D], 469 (. 8 T. Nakamura, R. Tanaka, T. Yabe, K. Takizawa, "Exactly conservative semi-lagrangian scheme for multi-dimensional hyperbolic equations with directional splitting technique," J. Comput. Phys., 174, 171 (1. PROGRESS IN NUCLEAR SCIENCE AND TECHNOLOGY
7 Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM K. Takizawa, T. Yabe, T. Nakamura, "Multi-dimensional semi- Lagrangian scheme that guarantees exact conservation," Comput. Phys. Commun., 148, 137 (. 1 T. Yabe, P. Y. Wang, "Unified numerical procedure for compressible and incompressible fluid," J. Phys. Soc. Japan, 6, 15 ( A. Ikebata, F. Xiao, "Development of conservative front capturing scheme and applications to multi-fluid simulations," Proc. Annual Conference of Japan Society of Mechanical Engineering, 3, 31 (, [in Japanese]. 1 F. Xiao, T. Yabe, T. Ito, M. Tajima, "An algorithm for simulating solid objects suspended in stratified flow," Comput. Phys. Commun., 1, 147 ( F. Xiao, "A computational model for suspended large rigid body in 3D unsteady viscous flow," J. Comput. Phys., 155, 348 ( J. U. Brackbill, D. B. Kothe, C. Zemach, "A continuum method for modeling surface tension," J. Comput. Phys., 1, 335 (199. VOL., OCTOBER 11
A 3D VOF model in cylindrical coordinates
A 3D VOF model in cylindrical coordinates Marmar Mehrabadi and Markus Bussmann Department of Mechanical and Industrial Engineering, University of Toronto Recently, volume of fluid (VOF) methods have improved
More informationTwo-Phase flows on massively parallel multi-gpu clusters
Two-Phase flows on massively parallel multi-gpu clusters Peter Zaspel Michael Griebel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität Bonn Workshop Programming of Heterogeneous
More informationOverview of Traditional Surface Tracking Methods
Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More informationOn the thickness of discontinuities computed by THINC and RK schemes
The 9th Computational Fluid Dynamics Symposium B7- On the thickness of discontinuities computed by THINC and RK schemes Taku Nonomura, ISAS, JAXA, Sagamihara, Kanagawa, Japan, E-mail:nonomura@flab.isas.jaxa.jp
More informationDroplet collisions using a Level Set method: comparisons between simulation and experiments
Computational Methods in Multiphase Flow III 63 Droplet collisions using a Level Set method: comparisons between simulation and experiments S. Tanguy, T. Ménard & A. Berlemont CNRS-UMR6614-CORIA, Rouen
More informationInteraction of Fluid Simulation Based on PhysX Physics Engine. Huibai Wang, Jianfei Wan, Fengquan Zhang
4th International Conference on Sensors, Measurement and Intelligent Materials (ICSMIM 2015) Interaction of Fluid Simulation Based on PhysX Physics Engine Huibai Wang, Jianfei Wan, Fengquan Zhang College
More informationCUDA. Fluid simulation Lattice Boltzmann Models Cellular Automata
CUDA Fluid simulation Lattice Boltzmann Models Cellular Automata Please excuse my layout of slides for the remaining part of the talk! Fluid Simulation Navier Stokes equations for incompressible fluids
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationDrop Impact Simulation with a Velocity-Dependent Contact Angle
ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Drop Impact Simulation with a Velocity-Dependent Contact Angle S. Afkhami and M. Bussmann Department
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationWater. Notes. Free surface. Boundary conditions. This week: extend our 3D flow solver to full 3D water We need to add two things:
Notes Added a 2D cross-section viewer for assignment 6 Not great, but an alternative if the full 3d viewer isn t working for you Warning about the formulas in Fedkiw, Stam, and Jensen - maybe not right
More informationTutorial 17. Using the Mixture and Eulerian Multiphase Models
Tutorial 17. Using the Mixture and Eulerian Multiphase Models Introduction: This tutorial examines the flow of water and air in a tee junction. First you will solve the problem using the less computationally-intensive
More informationThree Dimensional Numerical Simulation of Turbulent Flow Over Spillways
Three Dimensional Numerical Simulation of Turbulent Flow Over Spillways Latif Bouhadji ASL-AQFlow Inc., Sidney, British Columbia, Canada Email: lbouhadji@aslenv.com ABSTRACT Turbulent flows over a spillway
More informationINTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS
INTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS Yoshinobu Yamamoto Department of Nuclear Engineering, Kyoto University Yoshida Sakyo Kyoto, 66-85,
More informationPhase-field simulation of two-phase micro-flows in a Hele-Shaw cell
Computational Methods in Multiphase Flow III 7 Phase-field simulation of two-phase micro-flows in a Hele-Shaw cell Y. Sun & C. Beckermann Department of Mechanical and Industrial Engineering, University
More informationA Contact Angle Model for the Parallel Free Surface Lattice Boltzmann Method in walberla Stefan Donath (stefan.donath@informatik.uni-erlangen.de) Computer Science 10 (System Simulation) University of Erlangen-Nuremberg
More informationCS 231. Fluid simulation
CS 231 Fluid simulation Why Simulate Fluids? Feature film special effects Computer games Medicine (e.g. blood flow in heart) Because it s fun Fluid Simulation Called Computational Fluid Dynamics (CFD)
More informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationImplicit surface tension model for stimulation of interfacial flows
surface tension model for stimulation of interfacial flows Vinh The University of Massachusetts Dartmouth March 3rd, 2011 About my project Project Advisor Dr. Mehdi Raessi Department of Mechanical Engineering
More information(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon
More informationEfficient Tridiagonal Solvers for ADI methods and Fluid Simulation
Efficient Tridiagonal Solvers for ADI methods and Fluid Simulation Nikolai Sakharnykh - NVIDIA San Jose Convention Center, San Jose, CA September 21, 2010 Introduction Tridiagonal solvers very popular
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationFinal drive lubrication modeling
Final drive lubrication modeling E. Avdeev a,b 1, V. Ovchinnikov b a Samara University, b Laduga Automotive Engineering Abstract. In this paper we describe the method, which is the composition of finite
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationAcknowledgements. Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn. SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar
Philipp Hahn Acknowledgements Prof. Dan Negrut Prof. Darryl Thelen Prof. Michael Zinn SBEL Colleagues: Hammad Mazar, Toby Heyn, Manoj Kumar 2 Outline Motivation Lumped Mass Model Model properties Simulation
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationSPH: Why and what for?
SPH: Why and what for? 4 th SPHERIC training day David Le Touzé, Fluid Mechanics Laboratory, Ecole Centrale de Nantes / CNRS SPH What for and why? How it works? Why not for everything? Duality of SPH SPH
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationA mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates
A mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates Peter Hjort Lauritzen Atmospheric Modeling & Predictability Section National Center for Atmospheric
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationLagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion
Lagrangian and Eulerian Representations of Fluid Flow: Kinematics and the Equations of Motion James F. Price Woods Hole Oceanographic Institution Woods Hole, MA, 02543 July 31, 2006 Summary: This essay
More informationShape of Things to Come: Next-Gen Physics Deep Dive
Shape of Things to Come: Next-Gen Physics Deep Dive Jean Pierre Bordes NVIDIA Corporation Free PhysX on CUDA PhysX by NVIDIA since March 2008 PhysX on CUDA available: August 2008 GPU PhysX in Games Physical
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationRobust Simulation of Sparsely Sampled Thin Features in SPH-Based Free Surface Flows
Copyright of figures and other materials in the paper belong to original authors. Robust Simulation of Sparsely Sampled Thin Features in SPH-Based Free Surface Flows Xiaowei He et al. ACM SIGGRAPH 2015
More informationAb initio NMR Chemical Shift Calculations for Biomolecular Systems Using Fragment Molecular Orbital Method
4 Ab initio NMR Chemical Shift Calculations for Biomolecular Systems Using Fragment Molecular Orbital Method A Large-scale Two-phase Flow Simulation Evolutive Image/ Video Coding with Massively Parallel
More informationSimulation of Fuel Sloshing Comparative Study
3. LS-DYNA Anwenderforum, Bamberg 2004 Netfreie Verfahren Simulation of Fuel Sloshing Comparative Study Matej Vesenjak 1, Heiner Müllerschön 2, Alexander Hummel 3, Zoran Ren 1 1 University of Maribor,
More informationEVALUATION OF VOLUME TRACKING ALGORITHMS FOR GAS-LIQUID TWO-PHASE FLOWS
Proceedings of ASME FEDSM 03 4 th ASME_JSME Joint Fluids Engineering Conference Honolulu, Hawaii, USA, July 6-10, 2003 F E DSM200 3-4 5164 EVALUATION OF VOLUME TRACKING ALGORITHMS FOR GAS-LIQUID TWO-PHASE
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationENERGY-224 Reservoir Simulation Project Report. Ala Alzayer
ENERGY-224 Reservoir Simulation Project Report Ala Alzayer Autumn Quarter December 3, 2014 Contents 1 Objective 2 2 Governing Equations 2 3 Methodolgy 3 3.1 BlockMesh.........................................
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationFree Surface Flow Simulations
Free Surface Flow Simulations Hrvoje Jasak h.jasak@wikki.co.uk Wikki Ltd. United Kingdom 11/Jan/2005 Free Surface Flow Simulations p.1/26 Outline Objective Present two numerical modelling approaches for
More informationMoment of fluid method for multimaterial flows
Moment of fluid method for multimaterial flows M. Jemison, M. Sussman, M. Shashkov Department of Mathematics, Florida State University October 1, 2012 Recent Motivating Work Volume preserving Moment-of-Fluid
More informationVolume Tracking on Adaptively Refined Grids with Curvature Based Refinement
Volume Tracking on Adaptively Refined Grids with Curvature Based Refinement Mayank Malik, Markus Bussmann Department of Mechanical & Industrial Engineering, University of Toronto mayank.malik@utoronto.ca,
More informationMETHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS
XVIII International Conference on Water Resources CMWR 2010 J. Carrera (Ed) CIMNE, Barcelona 2010 METHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS Poorya A. Ferdowsi
More informationSurface Tension Approximation in Semi-Lagrangian Level Set Based Fluid Simulations for Computer Graphics
Surface Tension Approximation in Semi-Lagrangian Level Set Based Fluid Simulations for Computer Graphics Israel Pineda and Oubong Gwun Chonbuk National University israel_pineda_arias@yahoo.com, obgwun@jbnu.ac.kr
More informationMoving Interface Problems: Methods & Applications Tutorial Lecture II
Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,
More informationExperimental Validation of the Computation Method for Strongly Nonlinear Wave-Body Interactions
Experimental Validation of the Computation Method for Strongly Nonlinear Wave-Body Interactions by Changhong HU and Masashi KASHIWAGI Research Institute for Applied Mechanics, Kyushu University Kasuga
More informationA Novel Approach to High Speed Collision
A Novel Approach to High Speed Collision Avril Slone University of Greenwich Motivation High Speed Impact Currently a very active research area. Generic projectile- target collision 11 th September 2001.
More informationSIMPLIFIED VOLUME OF FLUID METHOD (SVOF) FOR TWO-PHASE FLOWS
TASK QUARTERLY 12 No 3, 255 265 SIMPLIFIED VOLUME OF FLUID METHOD (SVOF) FOR TWO-PHASE FLOWS MACIEJ MAREK, WOJCIECH ANISZEWSKI AND ANDRZEJ BOGUSŁAWSKI Institute of Thermal Machinery, Czestochowa University
More informationStudy on the Numerical Accuracy for the CFD
Study on the Numerical Accuracy for the CFD T.Yamanashi 1, H.Uchida, and M.Morita 1 Department of Mathematics, Master s Research Course of Faculty of Science, Tokyo University of Science,1-3 Kagurazaka,
More informationAdarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs)
OBJECTIVE FLUID SIMULATIONS Adarsh Krishnamurthy (cs184-bb) Bela Stepanova (cs184-bs) The basic objective of the project is the implementation of the paper Stable Fluids (Jos Stam, SIGGRAPH 99). The final
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
More informationRealistic Animation of Fluids
1 Realistic Animation of Fluids Nick Foster and Dimitris Metaxas Presented by Alex Liberman April 19, 2005 2 Previous Work Used non physics-based methods (mostly in 2D) Hard to simulate effects that rely
More informationA brief description of the particle finite element method (PFEM2). Extensions to free surface
A brief description of the particle finite element method (PFEM2). Extensions to free surface flows. Juan M. Gimenez, L.M. González, CIMEC Universidad Nacional del Litoral (UNL) Santa Fe, Argentina Universidad
More informationPARALLEL SIMULATION OF A FLUID FLOW BY MEANS OF THE SPH METHOD: OPENMP VS. MPI COMPARISON. Pawe l Wróblewski, Krzysztof Boryczko
Computing and Informatics, Vol. 28, 2009, 139 150 PARALLEL SIMULATION OF A FLUID FLOW BY MEANS OF THE SPH METHOD: OPENMP VS. MPI COMPARISON Pawe l Wróblewski, Krzysztof Boryczko Department of Computer
More informationComputational Fluid Dynamics - Incompressible Flows
Computational Fluid Dynamics - Incompressible Flows March 25, 2008 Incompressible Flows Basis Functions Discrete Equations CFD - Incompressible Flows CFD is a Huge field Numerical Techniques for solving
More informationMultiphysics simulations of nuclear reactors and more
Multiphysics simulations of nuclear reactors and more Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg klasjareteg@chalmersse Division of Nuclear Engineering Department of Applied Physics Chalmers
More informationUse 6DOF solver to calculate motion of the moving body. Create TIFF files for graphic visualization of the solution.
Introduction The purpose of this tutorial is to provide guidelines and recommendations for setting up and solving a moving deforming mesh (MDM) case along with the six degree of freedom (6DOF) solver and
More informationCS-184: Computer Graphics Lecture #21: Fluid Simulation II
CS-184: Computer Graphics Lecture #21: Fluid Simulation II Rahul Narain University of California, Berkeley Nov. 18 19, 2013 Grid-based fluid simulation Recap: Eulerian viewpoint Grid is fixed, fluid moves
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationThe Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection
The Development of a Volume-of-Fluid Interface Tracking Method for Modeling Problems in Mantle Convection Jonathan Robey 1 September 12, 2016 1 Funded in part by a GAANN Fellowship 1 Motivation Mantle
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationUsing the Eulerian Multiphase Model for Granular Flow
Tutorial 21. Using the Eulerian Multiphase Model for Granular Flow Introduction Mixing tanks are used to maintain solid particles or droplets of heavy fluids in suspension. Mixing may be required to enhance
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationParticle-Based Fluid Simulation. CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016
Particle-Based Fluid Simulation CSE169: Computer Animation Steve Rotenberg UCSD, Spring 2016 Del Operations Del: = x Gradient: s = s x y s y z s z Divergence: v = v x + v y + v z x y z Curl: v = v z v
More informationSPH: Towards the simulation of wave-body interactions in extreme seas
SPH: Towards the simulation of wave-body interactions in extreme seas Guillaume Oger, Mathieu Doring, Bertrand Alessandrini, and Pierre Ferrant Fluid Mechanics Laboratory (CNRS UMR6598) Ecole Centrale
More informationFluid Simulation. [Thürey 10] [Pfaff 10] [Chentanez 11]
Fluid Simulation [Thürey 10] [Pfaff 10] [Chentanez 11] 1 Computational Fluid Dynamics 3 Graphics Why don t we just take existing models from CFD for Computer Graphics applications? 4 Graphics Why don t
More informationSIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH
Fifth International Conference on CFD in the Process Industries CSIRO, Melbourne, Australia 3-5 December 006 SIMULATION OF INTERFACIAL FLOWS USING FRONT TRACKING APPROACH Manasa Ranjan Behera and Murali.
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationParallel Interpolation in FSI Problems Using Radial Basis Functions and Problem Size Reduction
Parallel Interpolation in FSI Problems Using Radial Basis Functions and Problem Size Reduction Sergey Kopysov, Igor Kuzmin, Alexander Novikov, Nikita Nedozhogin, and Leonid Tonkov Institute of Mechanics,
More informationTransfer and pouring processes of casting by smoothed particle. hydrodynamic method
Transfer and pouring processes of casting by smoothed particle hydrodynamic method M. Kazama¹, K. Ogasawara¹, *T. Suwa¹, H. Ito 2, and Y. Maeda 2 1 Application development div., Next generation technical
More informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
More informationAdvances in 3-D Forging Process Modeling
Advances in 3-D Forging Process Modeling W J Slagter, C J L Florie and A C J Venis MSC.Software Corporation, Groningenweg 6, 2803 PV Gouda, The Netherlands MSC/SuperForge is a new code developed for performing
More informationNumerical simulations of dam-break floods with MPM
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (206) 000 000 www.elsevier.com/locate/procedia st International Conference on the Material Point Method, MPM 207 Numerical
More informationSweeping Jet Flushing Method for EDM
International Conference on Multidisciplinary Research & Practice P a g e 192 Sweeping Jet Flushing Method for EDM Kishor K Kulkarni # Production Engineering, ShivajiUniversity/Ashokrao Mane Group of Institution
More informationAurélien Thinat Stéphane Cordier 1, François Cany
SimHydro 2012:New trends in simulation - Hydroinformatics and 3D modeling, 12-14 September 2012, Nice Aurélien Thinat, Stéphane Cordier, François Cany Application of OpenFOAM to the study of wave loads
More informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
More information2.7 Cloth Animation. Jacobs University Visualization and Computer Graphics Lab : Advanced Graphics - Chapter 2 123
2.7 Cloth Animation 320491: Advanced Graphics - Chapter 2 123 Example: Cloth draping Image Michael Kass 320491: Advanced Graphics - Chapter 2 124 Cloth using mass-spring model Network of masses and springs
More informationFlux Vector Splitting Methods for the Euler Equations on 3D Unstructured Meshes for CPU/GPU Clusters
Flux Vector Splitting Methods for the Euler Equations on 3D Unstructured Meshes for CPU/GPU Clusters Manfred Liebmann Technische Universität München Chair of Optimal Control Center for Mathematical Sciences,
More informationModeling Unsteady Compressible Flow
Tutorial 4. Modeling Unsteady Compressible Flow Introduction In this tutorial, FLUENT s density-based implicit solver is used to predict the timedependent flow through a two-dimensional nozzle. As an initial
More informationA High Quality, Eulerian 3D Fluid Solver in C++ A Senior Project. presented to. the Faculty of the Computer Science Department of
A High Quality, Eulerian 3D Fluid Solver in C++ A Senior Project presented to the Faculty of the Computer Science Department of California Polytechnic State University, San Luis Obispo In Partial Fulfillment
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Pop Worksheet! Teams of 2. Hand in to Jeramey after we discuss. Sketch the first few frames of a 2D explicit Euler mass-spring simulation for a 2x3 cloth network
More informationInteractive Fluid Simulation using Augmented Reality Interface
Interactive Fluid Simulation using Augmented Reality Interface Makoto Fuisawa 1, Hirokazu Kato 1 1 Graduate School of Information Science, Nara Institute of Science and Technology, 8916-5 Takayama, Ikoma,
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationA TSTT integrated FronTier code and its applications in computational fluid physics
Institute of Physics Publishing Journal of Physics: Conference Series 16 (2005) 471 475 doi:10.1088/1742-6596/16/1/064 SciDAC 2005 A TSTT integrated FronTier code and its applications in computational
More informationAMR Multi-Moment FVM Scheme
Chapter 4 AMR Multi-Moment FVM Scheme 4.1 Berger s AMR scheme An AMR grid with the hierarchy of Berger s AMR scheme proposed in [13] for CFD simulations is given in Fig.4.1 as a simple example for following
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationNew Conjugate-Heat Transfer Solvers in the Compressible CESE Solver in LS-DYNA
New Conjugate-Heat Transfer Solvers in the Compressible CESE Solver in LS-DYNA Grant O. Cook, Jr. and Zeng-Chan Zhang Livermore Software Technology Corp., Livermore, CA 94551 Abstract Standard coupling
More informationChapter 2 Water and Bubbles
Chapter 2 Water and Bubbles 2.1 Animation of Bubbles in Liquid Abstract This section introduces a new fluid animation technique in which liquid and gas interact with each other, using the example of bubbles
More information