Parallel lattice fluid simulations for transport phenomena in porous media
|
|
- Samson Pope
- 5 years ago
- Views:
Transcription
1 Parallel lattice fluid simulations for transport phenomena in porous media T. Watanabel, K. Kono2 & H, 0hashi3 ljapan Atomic Energy Research Institute, Japan. Fuji Research Institute Corporation, Japan. 3Univemity of Tokyo, Japan. Abstract Parallel simulation programs using the lattice fluid methods are developed for studying the transport phenomena in porous media. The domain decomposition technique is applied and the message passing interface library is used for parallelization. Sample problems are solved to validate the simulation programs and evaluate parallel efficiency. It is shown for the LBGK simulations, in which the amount of data transfer between two processors is constant, that the on-nodeparallel calculation using shared memory is effective for the calculation time and the memory usage in comparison with the purely distributed calculations. The efficiency of data transfer for the RLG simulations, in which the amount of data transfer varies every time step, is found to be better for sending the number of data and some amount of particle data at the same time. 1 Introduction It is of importance for safety assessment of geological disposal of high-level radioactive waste (HLW) to predict the transport phenomena in the natural barrier such as rock matrix. Numerical simulations of those phenomena are extremely difficult since the material transport mechanism of HLW is complicated due to geometrical complexity of rock matrix and multiple physical and chemical phenomena in groundwater. A joint research project is being carried out by University of Tokyo, Fuji Research Institute and Japan Atomic Energy Research Institute to establish a simulation system for geological disposal of HLW. The project consists of a development of a multi-scale simulation method using lattice fluid models and a
2 ~~6 Applications of High-Performance Computing in Engineering VII development of virtual experiment environment using high performance computing and virtual reality techniques [1]. The lattice fluid models, which are based on the dynamics of discrete fluid particles, have recently progressed a great deal. Motions of fluid particles are calculated instead of solving continuous fluid equations, and macroscopic flow variables are obtained from the particle motion. Among particle simulation methods, the lattice gas automata (LGA) is one of the simple techniques for simulating complicated flow phenomena, as well as macroscopic flow fields. In the LGA introduced by Frisch, Hasslacher and Pomeau [2], space and time are discrete, and identical particles of equal mass populate a triangular lattice. The particles travel to neighboring sites at each time step, and obey simple collision rules that conserve mass and momentum. Macroscopic flow fields are obtained by coarse-grain averaging in space and time. The LGA has, however, some inherent drawbacks such as the lack of Galilean invariance, and several lattice fluid models have been proposed to overcome the drawbacks. The lattice BGK (LBGK) method and the real-coded lattice gas (RLG) method are used in our project to simulate the complicated transport phenomena in the natural barrier. In the LBGK method, the particle occupation variables, which are Boolean variables in the LGA, are replaced with single-particle distribution functions [3], and the collision process of particles is simplified to become a relaxation process of the distribution fimction toward the local equilibrium [4], The distribution fimctions are real variables as the ensemble average of particle occupation, and the local equilibrium distribution is chosen to recover the Navier-Stokes equations [5]. The RLG method has been proposed by Malevanets and Kapral [6]. In this method, the velocities of particles are represented by real variables, while they are Boolean variables in the LGA. The motion of particles are thus not restricted on the lattice. A stochastic propagation scheme with a multi-particle collision rule that conserve mass, momentum and energy is applied to the particles. Since the algorithm and programming are simple, and complex boundary geometries are easy to represent, the LBGK and the RLG methods have been applied to numerical simulations of various kinds of hydrodpamic flows [7,8,9,10]. The complicated structure of rock matrix and various kinds of flows with physical and chemical phenomena are modeled at the microscopic level using the LBGK and the RLG methods in this project. The domain decomposition technique is applied and the message passing interface library is used for parallelization. In this study, characteristics of the LBGK and the RLG methods in parallelization are described and the parallel efficiency is evaluated. 2 Lattice fluid models 2.2 Lattice BGK method In the LBGK method, a particle distribution functions fi(x,t) at space x and time t are introduced to represent a fluid, where the subscript i indicates the direction of
3 Applications of High-Performance Computing in Engineering VII 117 moving particle on a lattice. In this study, the 3-D 15-direction model is used. The lattice Boltzrnann equation is given by fi(x+e~, t+l)-fi(x, t)= Q~(x, t), (1) where e is a particle velocity vector and the right hand side is the collision operator. The linearized collision operator with a single time relaxation parameter is used in the LBGK method: Qi.+.jfq. r Here f~q is the local equilibrium velocity, and ~ is a characteristic momentum must be satisfied: (2) state depending on the local density and relaxation time. Conservation of mass and (3) where p and u is the density and the local velocity of the fluid, respectively. Using the Chaprnan-Enskog multiscale expansion, the equilibrium distribution func;on is obtained by - for I ei 1=O, U) + ~(ei c24)2 :U2)I for ei 1=1, (4) 212 fieq =L[~+:(ei tl)+~(ei u) ;U )1 for I ei 1= J 3, 893 so that the macroscopic mass and momentum conservation equations for continuous fluid are recovered. 2.2 Real-coded lattice gas method In the RLG method, fluid particles have a continuous velocity in both spatial direction and magnitude. The trajectories of the particles are traced in a shorttime interval by decoupling interparticle collisions, In the translational process,
4 ~~8 Applications of High-Performance Computing in Engineering VII locations of the particles are obtained using the velocities calculated in the previous collision process, r(t + At) = r(t)+ V(t), (5) where r and v represent the location and the velocity of a particle. In the collision process, velocities of particles are modified while their locations remain unchanged, The velocities of the particles are obtained using the average velocity in a collision cell, V, and the rotation operator, o: v(t + At) = ~ + C7(V(t) ~). (6) The vector of velocity difference between the average velocity and the particle velocity is rotated by some angles defined by the rotation operator. The mass, momentum and energy are conserved during the collision process. 3 Parallel efficiency The lattice fluid methods consist of a streaming process and a collision process: the streaming and the relaxation of the distribution fi,mction in the LBGK method, and the streaming of particles and the multi-particle collision in the RLG method, The streaming and the collision processes are calculated locally, and high efficiency of parallel computation is expected in the lattice fluid methods. The efficiency of parallel computation for the lattice fluid methods is studied here, The domain decomposition technique is applied and the message passing interface (MPI) library is used. The parallel computer system, HITACHI SR8000, with 4 nodes is used, where one node consists of a shared memory and 8 processor elements (PEs). In this systeq purely distributed calculations are possible using 32 PEs. Variables and calculation procedures are divided into 32 parts, and the data transfer among PEs is performed using the MPI even in the same node, The on-node-parallel calculations are also possible using the shared memory. Variables and calculation procedures are divided into 4 parts in this case, and the data transfer among nodes is performed using the MPI. Calculation procedures in the node are, then, divided into 8 parts using the shared memory, and the data transfer is not necessary inside the node. The purely MPI calculations and the (MPI + on-node-parallel) calculations are compared and the efficiency of data transfer is discussed in the following. The workstation cluster with 4 COMPAQ XP 1000 workstations are also used for comparison. 3.1 LBGK simulations The sample problem for the LBGK method is the phase separation process from a homogeneous mixture of two fluids in 64x64x192 lattice nodes. The twocomponent two-phase model of the LBGK is used [7], The periodic boundary
5 Applications of High-Performance Computing in Engineering VII 119 conditions are applied to the top, botto~ and side boundaries, Examples of the calculated results are shown in Figure 1, where the formation of interface between two fluids is depicted every 200 time steps. It is shown that the phase separation occurs and the size of each fluid region increases. One of the separated phases shown in Figure 1 is used as the complicated structure of rock matrix. The simulation region consists of 64x64 lattice points in the horizontal direction, and 192 lattice points in the vertical direction. Variables and calculation procedures are divided into 4 small parts in the vertical direction using the domain decomposition method. Each small part is assigned to each node. The variables and calculation procedures in the node are further divided into 8 small parts for the purely MPI calculations, while the calculation procedures alone are divided for the (MPI + on-node-parallel) calculations. The amount of data transfer between two processors is constant for the purely MPI calculations and that between two nodes is constant for the (MPI + on-nodeparallel) calculations. Several specific commands are inserted in the simulation program in order to use the on-node-parallel functions, and thus some more modifications are necessary for the (MPI + on-node-parallel) calculations. The speedup is shown in Figure 2, where the (MPI + on-node-parallel) cases are compared with the purely MPI cases on SR8000 and on the workstation cluster. The average speedup of five calculations is shown. For the (MPI + onnode-parallel) calculations using 8 PEs in one node, the calculation time is about 693 sec and larger than 560 sec of the purely MPI calculations. The on-nodeparallel ilmction needs additional calculation time, though the data transfer is not necessary in the node since the shared memory is used. The speedup of the (MPI + on-node-parallel) calculations inone node is thus lower than that of the MPI calculations, The calculation time of the (MPI + on-node-parallel) calculations with 4 nodes is about221 sec and slightly smaller than 237 sec of the purely MPI calculations. The speedup ratio is about 17.3 for the (MPI + on-node-parallel) calculations with 4 nodes, and about 16.2 for the purely MPI calculations using 32 PEs, The time for data transfer is shown in Figure 3, where the rate of the data transfer time to the elapsed time is shown. Domain decomposition is performed in one direction, and the amount of data transfer for one PE is constant with increase in the number of PEs, though the amount of calculation is reduced, The time for data transfer using 32 PEs is about 45 A of the total calculation time for the purely MPI calculations, while about 35 % for the (MPI + on-node-parallel) calculations. The rate of the data transfer time to the elapsed time is smaller for the (MPI + on-node-parallel) calculations, but the additional calculations are necessary, and thus the speedup is not much different between the purely MPI and the (MPI + on-node-parallel) cases as shown in Figure 2. It is seen in Figure 2, however, that the on-node-parallel finction is effective as the number of nodes and PEs increases. The normalized memory size of simulation program is shown in Figure 4, The memory size for a single PE is reduced and the total memory size is increased for both cases. The memory size for the case using 8 PEs on one node is about 130
6 ~20 Applications of High-Performance Computing in Engineering VII MB for the (MPI + on-node-parallel) calculations and 480 MB for the purely MPI calculations. The difference of memory size is due to an additional buffer memory needed for the data transfer by the MPI. Variables are divided into 8 parts in one node for the purely MPI calculations, while they are not divided for the (MPI + on-node-parallel) calculations, The amount of buffer memory is thus about 8 times larger for the purely MPI calculations. It is found that the shared memorys ystem is effective from the viewpoint of the memory size. 3.2 RLG simulations In parallel RLG simulations, the domain decomposition method is applied and the MPI library is used. The fluid particles are used in the RLG simulations, and the particle data are exchanged among PEs during the calculations. The amount of data transfer between two PEs varies every time step, since the motion of the fluid particles is not restricted on the lattice. For the data transfer of the parallel RLG simulations, two send-receive procedures are generally required. The number of data is sent first and all the particle data are sent next from one PE to other PEs. In the following, the efficiency of data transfer is studied for sending the number of data and some amount of particle data fust and the remaining particle data next. The sample problem for the RLG method is the thermal convection using 1.28 million particles in 40x40x20 collision cells. The top and bottom boundaries are cooled and heated, respectively, while periodic boundary conditions are applied to the side boundaries. This sample problem known as the Rayleigh-Benard convection is used for a validation of the heat transfer calculations in the RLG method. The simulation region is divided two-dimensionally in the horizontal directions in this case. Examples of the calculated results are shown in Figure 5, where the formation of convection pattern at the midelevation of the simulation region is depicted. The particle velocities are averaged in the collision cells and the distribution of vertical flow velocities are shown. The speedup is shown in Figure 6, where the purely MPI calculations are performed, The number of particles sent from one PE to other PEs in one time step is about 250 for the case with 32 PEs. In Figure 6, some amount of particle data (Num Data) is sent f~st with the number of particles: Num Data Oindicates that the number of particles alone is sent fust and all the particle data are sent next, Num Data 5 indicates that 5 particle data are sent first with the number of particles and the remaining particle data next, Num Data 50 indicates that 50 particle data are sent first with the number of particles and the remaining particle data next, The speedup of these three cases is not much different as shown in Figure 6. The speedup is about 19.0 using 32 PEs. The time for data transfer in this case is about 11 YO of the total computation time. The number of particles assigned to one PE decreases as the number of PEs increases. The amount of data transfer, which is the number of particles sent from one PE to other PEs, decreases as well, and the time for data transfer is relatively smaller for the RLG method than for the LBGK method. The average data transfer time in the case with 32 PEs is shown in Figure 7
7 Applications of High-Performance Computing in Engineering VII 121 against the number of particle data sent fust, It is clearly seen that the data transfer time becomes small when the number of particle data is 42 and 170. In our RLG progra~ each particle has 3 velocities and 3 coordinates in 4 bite real numbers, and thus 42 and 170 particle data correspond to 1024 B and 4096 B, respectively, These two values are the buffer size for the send-receive process of the MPI. Although the buffer size depends on the computer systems, the efficiency of data transfer would be better for sending some amount of particle data first with the number of data. 4 Conclusion Parallel simulation programs using the LBGK and the RLG methods have been developed for studying the transport phenomena in porous media. The domain decomposition technique has been applied and the MPI library was used. Sample problems were solved to validate the simulation programs and evaluate parallel efficiency. It was shown for the LBGK method, in which the amount of data transfer is constant, that the on-node-parallel calculation using shared memory was effective for the calculation time and the memory usage in comparison with the purely distributed calculations. The efficiency of data transfer for the RLG method, in which the amount of data transfer varied every time step, was found to be better for sending some amount of particle data fust with the number of data. Acknowledgments This work was a part of the IAE s (Institute of Applied Energy) IVNET (Innovative and viable nuclear energy technology) projects funded by the Agency of Natural Resources and Energy, Ministry of Economy, Trade and Industry (METI), References [1] Itoh, G., Nakamura, J., Kono, K., Watanabe, T., Ohashi, H., Chen, Y,, & Nagasaki, S. Pore-scale simulation for predicting material transport through porous media. Proc. of the 10 hint. Con$ on Nuclear Engineering, 2002, [2] Frisch, U., Hasslacher, B. & Pomeau, Y. Lattice-gas automata for the Navier-Stokes equations. Phys. Rev. Lett., 56, pp , [3] McNamara, G.G. & Zanetti, G. Use of the Boltzmann equation to simulate lattice-gas automata. Phys. Rev, Lett., 61, pp , [4] Bhatnagar, P,L., Gross, E,P, & Krook, M. A model for collision processes in gases. I: small amplitude processes in charged and neutral one-component system. Phys. Rev., 94, pp , [5] Chen, H., Chen, S. & Matthaeus, W.H, Recovery of the Navier-Stokes equations using a lattice-gas Boltzmarm method. Phys. Rev. A, 45, pp. R5339-R5342, [6] Malevanets, A, & Kapral, R, Continuous-velocity lattice-gas model for fluid
8 ~22 Applications of High-Performance Computing in Engineering VII flow. Europhys. Lett., 44, pp , [7] Rothman, D.H. & Zaleski, S. Lattice-gas models of phase separation: interfaces, phase transitions, and multiphase flow. Rev. Mod. Phys., 66, pp , 1994, [8] Chen, H. & Doolen, G.D, Lattice Boltzmann method for fluid flows, Annu. Rev. Fluid Mech., 30, pp , [9] Hashimoto, Y., Chen, Y. & Ohashi, H. Immiscible real-coded lattice gas. Computer Phys. Comm., 129, pp ,2000. [10] Sakai, T., Chen, Y. & Ohashi, H. Formation of micelle in the real-coded lattice gas. Computer Phys. Comm., 129, pp ,2000. Figure 1: Simulation of phase separation using the LBGK method.? I 15.0 Cl ,0 0,0 0 EKE!EEb 1 I I.. Y. D Number d Processors Figure 2: Speedup of the LBGK method.
9 Applications of High-Performance Computing in Engineering VII 123 ~MPl --U MPI + (h Node Parallel, -d I o Number of Processors Figure 3: Data transfer time vs. number of PEs L I J 1.00 Total D. n PE CH3MPI - u MPI + On Node Parallel 0.01 \ Number of Processors Figure 4: Memory size for the LBGK method,
10 ~24 Applications of High-Performance Computing in Engineering VII Figure 5: Simulation of thermal convection using the RLG method EEizJ 1,, 0.0 o Number of Processors Figure 6: Speedup of the RLG method I 1 I 20.0 o Number of Data Figure 7: Data transfer time vs. number of data.
A Particle Cellular Automata Model for Fluid Simulations
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 36(2), 2009, Pages 35 41 ISSN: 1223-6934 A Particle Cellular Automata Model for Fluid Simulations Costin-Radu Boldea Abstract. A new cellular-automaton
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More information1. Mathematical Modelling
1. describe a given problem with some mathematical formalism in order to get a formal and precise description see fundamental properties due to the abstraction allow a systematic treatment and, thus, solution
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 informationThe gas-kinetic methods have become popular for the simulation of compressible fluid flows in the last
Parallel Implementation of Gas-Kinetic BGK Scheme on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer Middle East Technical
More informationA New Approach to Reduce Memory Consumption in Lattice Boltzmann Method on GPU
Journal of Applied Fluid Mechanics, Vol. 10, No. 1, pp. 55-67, 2017. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. DOI: 10.18869/acadpub.jafm.73.238.26431 A New Approach to Reduce
More informationA Lattice-Boltzmann Based Method Applied to Digital Rock Characterization of Perforation Tunnel Damage
SCA2016-058 1/6 A Lattice-Boltzmann Based Method Applied to Digital Rock Characterization of Perforation Tunnel Damage Bernd Crouse, David M Freed, Nils Koliha, Gana Balasubramanian EXA CORP Rajani Satti,
More informationImproved bounce-back methods for no-slip walls in lattice-boltzmann schemes: Theory and simulations
PHYSICAL REVIEW E 67, 066703 2003 Improved bounce-back methods for no-slip walls in lattice-boltzmann schemes: Theory and simulations M. Rohde,* D. Kandhai, J. J. Derksen, and H. E. A. Van den Akker Kramers
More informationThe Lattice Boltzmann Method used for fluid flow modeling in hydraulic components
The 15th Scandinavian International Conference on Fluid Power, SICFP 17, June 7-9, 2017, Linköping, Sweden The Lattice Boltzmann Method used for fluid flow modeling in hydraulic components Bernhard Manhartsgruber
More informationUse of HPC to Predict Porous Media Properties and Flow in Hydrocarbon Reservoirs
Use of HPC to Predict Porous Media Properties and Flow in Hydrocarbon Reservoirs Henry Neeman, 2 Dimitrios V. Papavassiliou 1 1 School of Chemical Engineering and Materials Science 2 School of Computer
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 informationSimulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method
Simulation of Liquid-Gas-Solid Flows with the Lattice Boltzmann Method June 21, 2011 Introduction Free Surface LBM Liquid-Gas-Solid Flows Parallel Computing Examples and More References Fig. Simulation
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 informationLattice Boltzmann Method: A New Alternative for Loss Estimation Juan C. Díaz Cruz* a, Sergio E. Rodríguez a, Sergio A. Ordóñez b, Felipe Muñoz a a
799 VOL. 31, 2013 CHEMICAL ENGINEERING TRANSACTIONS A publication of The Italian Association of Chemical Engineering Online at: www.aidic.it/cet Guest Editors: Eddy De Rademaeker, Bruno Fabiano, Simberto
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationWebinar #3 Lattice Boltzmann method for CompBioMed (incl. Palabos)
Webinar series A Centre of Excellence in Computational Biomedicine Webinar #3 Lattice Boltzmann method for CompBioMed (incl. Palabos) 19 March 2018 The webinar will start at 12pm CET / 11am GMT Dr Jonas
More informationMaterials Modelling MPhil
Materials Modelling MPhil COURSE MP5: MESOSCALE AND MULTISCALE MODELLING COMPUTING CLASS 1 31/1/2007 14:00-16:00 Dissipative particle dynamics using Materials Studio 4.0 1 Aims and objectives Dissipative
More informationModeling Evaporating Liquid Spray
Tutorial 16. Modeling Evaporating Liquid Spray Introduction In this tutorial, FLUENT s air-blast atomizer model is used to predict the behavior of an evaporating methanol spray. Initially, the air flow
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 informationLATTICE-BOLTZMANN AND FINITE ELEMENT SIMULATIONS OF FLUID FLOW IN A SMRX STATIC MIXER REACTOR
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 31: 1019 1033 (1999) LATTICE-BOLTZMANN AND FINITE ELEMENT SIMULATIONS OF FLUID FLOW IN A SMRX STATIC MIXER REACTOR D. KANDHAI
More informationRelative Permeability Simulation using the Two-phase Lattice-Boltzmann Method
5th Conference & Exposition on Petroleum Geophysics, Hyderabad-2004, India PP 696-703 Relative Permeability Simulation using the Two-phase Lattice-Boltzmann Method Youngseuk Keehm, Tapan Mukerji & Amos
More informationINNOVATIVE CFD FOR SUPER-COMPUTER RESULTS ON YOUR DESKTOP
INNOVATIVE CFD FOR SUPER-COMPUTER RESULTS ON YOUR DESKTOP XFlow is a next generation CFD software that uses a proprietary, particle-based, meshless approach which can easily handle traditionally complex
More informationThe Lattice-Boltzmann Method for Gaseous Phenomena
The Lattice-Boltzmann Method for Gaseous Phenomena Xiaoming Wei 1, Wei Li 1, Klaus Mueller 1 and Arie Kaufman 1 Center for Visual Computing (CVC) and Department of Computer Science State University of
More informationTowards real-time prediction of Tsunami impact effects on nearshore infrastructure
Towards real-time prediction of Tsunami impact effects on nearshore infrastructure Manfred Krafczyk & Jonas Tölke Inst. for Computational Modeling in Civil Engineering http://www.cab.bau.tu-bs.de 24.04.2007
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 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 informationPreliminary Spray Cooling Simulations Using a Full-Cone Water Spray
39th Dayton-Cincinnati Aerospace Sciences Symposium Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray Murat Dinc Prof. Donald D. Gray (advisor), Prof. John M. Kuhlman, Nicholas L. Hillen,
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 Dynamic Solver Based on Cellular Discrete-Event
Computational Fluid Dynamic Solver Based on Cellular Discrete-Event Simulation Michael Van Schyndel, Gabriel Wainer, Mohammad Moallemi Dept. of Systems & Computer Engineering, Carleton University, Ottawa,
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Christian Wollblad Copyright 2017 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
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 informationHigh Performance Computing
High Performance Computing ADVANCED SCIENTIFIC COMPUTING Dr. Ing. Morris Riedel Adjunct Associated Professor School of Engineering and Natural Sciences, University of Iceland Research Group Leader, Juelich
More informationSoftware and Performance Engineering for numerical codes on GPU clusters
Software and Performance Engineering for numerical codes on GPU clusters H. Köstler International Workshop of GPU Solutions to Multiscale Problems in Science and Engineering Harbin, China 28.7.2010 2 3
More informationReal-time Thermal Flow Predictions for Data Centers
Real-time Thermal Flow Predictions for Data Centers Using the Lattice Boltzmann Method on Graphics Processing Units for Predicting Thermal Flow in Data Centers Johannes Sjölund Computer Science and Engineering,
More informationAn Overview of Computational Fluid Dynamics
An Overview of Computational Fluid Dynamics Dr. Nor Azwadi bin Che Sidik Faculty of Mechanical Engineering Universiti Teknologi Malaysia INSPIRING CREATIVE AND INNOVATIVE MINDS 1 What is CFD? C computational
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 informationModeling Evaporating Liquid Spray
Tutorial 17. Modeling Evaporating Liquid Spray Introduction In this tutorial, the air-blast atomizer model in ANSYS FLUENT is used to predict the behavior of an evaporating methanol spray. Initially, the
More informationAnimation of Fluids. Animating Fluid is Hard
Animation of Fluids Animating Fluid is Hard Too complex to animate by hand Surface is changing very quickly Lots of small details In short, a nightmare! Need automatic simulations AdHoc Methods Some simple
More informationDynamic Mode Decomposition analysis of flow fields from CFD Simulations
Dynamic Mode Decomposition analysis of flow fields from CFD Simulations Technische Universität München Thomas Indinger Lukas Haag, Daiki Matsumoto, Christoph Niedermeier in collaboration with Agenda Motivation
More informationLATTICE-BOLTZMANN AND COMPUTATIONAL FLUID DYNAMICS
LATTICE-BOLTZMANN AND COMPUTATIONAL FLUID DYNAMICS NAVIER-STOKES EQUATIONS u t + u u + 1 ρ p = Ԧg + ν u u=0 WHAT IS COMPUTATIONAL FLUID DYNAMICS? Branch of Fluid Dynamics which uses computer power to approximate
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
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 informationIntroduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich
Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat
More informationVirtual EM Inc. Ann Arbor, Michigan, USA
Functional Description of the Architecture of a Special Purpose Processor for Orders of Magnitude Reduction in Run Time in Computational Electromagnetics Tayfun Özdemir Virtual EM Inc. Ann Arbor, Michigan,
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 informationMore Animation Techniques
CS 231 More Animation Techniques So much more Animation Procedural animation Particle systems Free-form deformation Natural Phenomena 1 Procedural Animation Rule based animation that changes/evolves over
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 informationLarge Eddy Simulations of a Stirred Tank Using the Lattice Boltzmann Method on a Nonuniform Grid
Journal of Computational Physics 181, 675 704 (2002) doi:10.1006/jcph.2002.7151 Large Eddy Simulations of a Stirred Tank Using the Lattice Boltzmann Method on a Nonuniform Grid Zhenyu Lu, Ying Liao, Dongying
More informationA 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 informationParametric Study of Sloshing Effects in the Primary System of an Isolated LFR Marti Jeltsov, Walter Villanueva, Pavel Kudinov
1 Parametric Study of Sloshing Effects in the Primary System of an Isolated LFR 19.06.2013 Marti Jeltsov, Walter Villanueva, Pavel Kudinov Division of Nuclear Power Safety Royal Institute of Technology
More informationSailfish: Lattice Boltzmann Fluid Simulations with GPUs and Python
Sailfish: Lattice Boltzmann Fluid Simulations with GPUs and Python Micha l Januszewski Institute of Physics University of Silesia in Katowice, Poland Google GTC 2012 M. Januszewski (IoP, US) Sailfish:
More informationThe use of gas-kinetic schemes for the simulation of compressible flows become widespread in the two last
A Gas-Kinetic BGK Scheme for Parallel Solution of 3-D Viscous Flows on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer
More informationLarge Scale Parallel Lattice Boltzmann Model of Dendritic Growth
Large Scale Parallel Lattice Boltzmann Model of Dendritic Growth Bohumir Jelinek Mohsen Eshraghi Sergio Felicelli CAVS, Mississippi State University March 3-7, 2013 San Antonio, Texas US Army Corps of
More informationLattice-Boltzmann simulations in reconstructed parametrized porous media
International Journal of Computational Fluid Dynamics, Vol. 00, No. 00, December 2005, 115 Lattice-Boltzmann simulations in reconstructed parametrized porous media BENJAMIN AHRENHOLZ, JONAS TÖLKE, MANFRED
More informationAndrew Carter. Vortex shedding off a back facing step in laminar flow.
Flow Visualization MCEN 5151, Spring 2011 Andrew Carter Team Project 2 4/6/11 Vortex shedding off a back facing step in laminar flow. Figure 1, Vortex shedding from a back facing step in a laminar fluid
More informationPhys 113 Final Project
Phys 113 Final Project LBM model for capillary fluid modelling Andrey Sushko Introduction Lattice Boltzmann Method (LMB) models have recently emerged as a highly effective means of modelling complex fluid
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 informationCloud-based simulation of plasma sources for surface treatment
Cloud-based simulation of plasma sources for surface treatment Using the PlasmaSolve Simulation Suite (P3S) Adam Obrusnik, Petr Zikan June 7, 2018 Outline 1. About PlasmaSolve 2. PlasmaSolve Simulation
More informationHigh Scalability of Lattice Boltzmann Simulations with Turbulence Models using Heterogeneous Clusters
SIAM PP 2014 High Scalability of Lattice Boltzmann Simulations with Turbulence Models using Heterogeneous Clusters C. Riesinger, A. Bakhtiari, M. Schreiber Technische Universität München February 20, 2014
More informationInvestigation of mixing chamber for experimental FGD reactor
Investigation of mixing chamber for experimental FGD reactor Jan Novosád 1,a, Petra Danová 1 and Tomáš Vít 1 1 Department of Power Engineering Equipment, Faculty of Mechanical Engineering, Technical University
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationEfficient Parallel Multi-Copy Simulation of the Hybrid Atomistic-Continuum Simulation Based on Geometric Coupling
, pp.384-389 http://dx.doi.org/10.14257/astl.2016.139.77 Efficient Parallel Multi-Copy Simulation of the Hybrid Atomistic-Continuum Simulation Based on Geometric Coupling Qian Wang 1,2, Xin-Hai Xu 1,2,*,
More informationLecture 1.1 Introduction to Fluid Dynamics
Lecture 1.1 Introduction to Fluid Dynamics 1 Introduction A thorough study of the laws of fluid mechanics is necessary to understand the fluid motion within the turbomachinery components. In this introductory
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 informationImplementation of a Lattice Gas Simulation
Implementation of a Lattice Gas Simulation Tony H. Kim Spring 2007 1 Introduction In this paper I discuss the implementation of a lattice-gas simulation on a hexagonal grid. The application is written
More informationApplication of Lattice Boltzmann and Cellular Automata in Erosion Modeling for Deep-Water Turbidity Current
Application of Lattice Boltzmann and Cellular Automata in Erosion Modeling for Deep-Water Turbidity Current Siyao Xu 1 Tapan Mukerji 1 Interdisciplinary Program of Earth, Energy and Environment Sciences
More informationParallel Direct Simulation Monte Carlo Computation Using CUDA on GPUs
Parallel Direct Simulation Monte Carlo Computation Using CUDA on GPUs C.-C. Su a, C.-W. Hsieh b, M. R. Smith b, M. C. Jermy c and J.-S. Wu a a Department of Mechanical Engineering, National Chiao Tung
More informationComputational Fluid Dynamics using OpenCL a Practical Introduction
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 2011 http://mssanz.org.au/modsim2011 Computational Fluid Dynamics using OpenCL a Practical Introduction T Bednarz
More informationLattice Boltzmann with CUDA
Lattice Boltzmann with CUDA Lan Shi, Li Yi & Liyuan Zhang Hauptseminar: Multicore Architectures and Programming Page 1 Outline Overview of LBM An usage of LBM Algorithm Implementation in CUDA and Optimization
More informationClick to edit Master title style
New features in Serpent 2 for fusion neutronics 5th International Serpent UGM, Knoxville, TN, Oct. 13-16, 2015 Jaakko Leppänen VTT Technical Research Center of Finland Click to edit Master title Outline
More informationHeat generation analysis of a rubber wheel using the steady-state transport analysis capability in Abaqus
Heat generation analysis of a rubber wheel using the steady-state transport analysis capability in Abaqus R. K. Luo 1, X. P. Wu 2 and A. Spinks 1 1 Trelleborg IAVS, 1 Hoods Close, Leicester, UK, LE4 2BN
More informationFree Surface Lattice-Boltzmann fluid simulations. with and without level sets.
ree Surface Lattice-Boltzmann fluid simulations with and without level sets Nils Thürey, Ulrich Rüde University of Erlangen-Nuremberg System Simulation roup Cauerstr. 6, 91054 Erlangen, ermany Email: Nils.Thuerey@cs.fau.de
More informationImplementation of a Lattice Gas Simulation
Implementation of a Lattice Gas Simulation Tony H. Kim Spring 2007 1 Introduction In this paper I discuss the implementation of a lattice-gas simulation on a hexagonal grid. The application is written
More informationAn interactive thermal fluid simulator for the design of HVAC systems
An interactive thermal fluid simulator for the design of HVAC systems Michael Pfaffinger, Christoph van Treeck, André Borrmann, Petra Wenisch, Ernst Rank Computational Civil and Environmental Engineering,
More informationPSG2 / Serpent a Monte Carlo Reactor Physics Burnup Calculation Code. Jaakko Leppänen
PSG2 / Serpent a Monte Carlo Reactor Physics Burnup Calculation Code Jaakko Leppänen Outline Background History The Serpent code: Neutron tracking Physics and interaction data Burnup calculation Output
More informationSIMULATION OF FLOW FIELD AROUND AND INSIDE SCOUR PROTECTION WITH PHYSICAL AND REALISTIC PARTICLE CONFIGURATIONS
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22, 2012 SIMULATION OF FLOW FIELD AROUND AND INSIDE SCOUR PROTECTION WITH PHYSICAL AND REALISTIC
More informationTopology optimization of heat conduction problems
Topology optimization of heat conduction problems Workshop on industrial design optimization for fluid flow 23 September 2010 Misha Marie Gregersen Anton Evgrafov Mads Peter Sørensen Technical University
More informationParallelization of Scientific Applications (II)
Parallelization of Scientific Applications (II) Parallelization of Particle Based Methods Russian-German School on High Performance Computer Systems, June, 27 th until July, 6 th 2005, Novosibirsk 4. Day,
More informationRealistic Animation of Fluids
Realistic Animation of Fluids p. 1/2 Realistic Animation of Fluids Nick Foster and Dimitri Metaxas Realistic Animation of Fluids p. 2/2 Overview Problem Statement Previous Work Navier-Stokes Equations
More informationMONTEREY, CALIFORNIA DISSERTATION LATTICE BOLTZMANN METHODS FOR FLUID STRUCTURE INTERACTION. Stuart R. Blair. September Dissertation Supervisor:
MONTEREY, CALIFORNIA DISSERTATION LATTICE BOLTZMANN METHODS FOR FLUID STRUCTURE INTERACTION by Stuart R. Blair September 2012 Dissertation Supervisor: Young Kwon Approved for public release; distribution
More informationLattice Boltzmann Method for Simulating Turbulent Flows
Lattice Boltzmann Method for Simulating Turbulent Flows by Yusuke Koda A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree of Master of Applied Science
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 informationc Fluent Inc. May 16,
Tutorial 1. Office Ventilation Introduction: This tutorial demonstrates how to model an office shared by two people working at computers, using Airpak. In this tutorial, you will learn how to: Open a new
More informationNavier-Stokes & Flow Simulation
Last Time? Navier-Stokes & Flow Simulation Implicit Surfaces Marching Cubes/Tetras Collision Detection & Response Conservative Bounding Regions backtracking fixing Today Flow Simulations in Graphics Flow
More informationPerformance Analysis of a Parallel Implementation of the Lattice Boltzmann Method for Computational Fluid Dynamics
Grand Valley State University ScholarWorks@GVSU Masters Theses Graduate Research and Creative Practice 4-2002 Performance Analysis of a Parallel Implementation of the Lattice Boltzmann Method for Computational
More informationPore Scale modeling (cont d)
Pore Scale modeling (cont d) Thomas Dewers Geomechanics Department Sandia National Laboratory Albuquerque NM tdewers@sandia.gov Outline Examples of modeling methods Variety and structure of caprock pore
More informationlecture 8 Groundwater Modelling -1
The Islamic University of Gaza Faculty of Engineering Civil Engineering Department Water Resources Msc. Groundwater Hydrology- ENGC 6301 lecture 8 Groundwater Modelling -1 Instructor: Dr. Yunes Mogheir
More informationDynamic Simulation of Flight Test Manoeuvres on the Diamond D-Jet
Dynamic Simulation of Flight Test Manoeuvres on the Diamond D-Jet The numerical simulation of the complex fluid-structure interaction taking place when manoeuvring an aircraft remains a challenge. A realistic
More informationDeformation and the Change of Filtration Properties of Weaves - A Computational Approach ABSTRACT
Deformation and the Change of Filtration Properties of Weaves - A Computational Approach S. Rief, E. Glatt, A. Wiegmann, M. Kabel, H. Andrä Fraunhofer-Institut Techno- und Wirtschaftsmathematik, Fraunhofer-Platz
More informationVerification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard
Verification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard Dimitri P. Tselepidakis & Lewis Collins ASME 2012 Verification and Validation Symposium May 3 rd, 2012 1 Outline
More information4200: Chemical Engineering Courses
4200: Chemical Engineering Courses 101 TOOLS FOR CHEMICAL ENGINEERING 3 credits Course Instructors: Dr George Chase Corequisites: 3450:149. Introduction to Chemical Engineering. Basic concepts of engineering
More informationOptimisation of LESsCOAL for largescale high-fidelity simulation of coal pyrolysis and combustion
Optimisation of LESsCOAL for largescale high-fidelity simulation of coal pyrolysis and combustion Kaidi Wan 1, Jun Xia 2, Neelofer Banglawala 3, Zhihua Wang 1, Kefa Cen 1 1. Zhejiang University, Hangzhou,
More informationPARALLEL PROGRAM COMPLEX FOR UNSTEADY FLOW SIMULATION
WEST-EAST HIGH SPEED FLOW FIELD CONFERENCE 19-, November 007 Moscow, Russia PARALLEL PROGRAM COMPLEX FOR UNSTEADY FLOW SIMULATION Eugene V. Shilnikov* * Institute for Mathematical Modeling RAS 15047, Miusskaya
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationLecture 4: Locality and parallelism in simulation I
Lecture 4: Locality and parallelism in simulation I David Bindel 6 Sep 2011 Logistics Distributed memory machines Each node has local memory... and no direct access to memory on other nodes Nodes communicate
More informationHIGH PERFORMANCE COMPUTATION (HPC) FOR THE
HIGH PERFORMANCE COMPUTATION (HPC) FOR THE DEVELOPMENT OF FLUIDIZED BED TECHNOLOGIES FOR BIOMASS GASIFICATION AND CO2 CAPTURE P. Fede, H. Neau, O. Simonin Université de Toulouse; INPT, UPS ; IMFT ; 31400
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationXFlow HIGH FIDELITY COMPUTATIONAL FLUID DYNAMICS
XFlow HIGH FIDELITY COMPUTATIONAL FLUID DYNAMICS XFlow OVERVIEW In the traditional mesh-based approach to solving Computational Fluid Dynamics (CFD) problems, reliability is highly dependent on the quality
More informationRendering Smoke & Clouds
Rendering Smoke & Clouds Game Design Seminar 2007 Jürgen Treml Talk Overview 1. Introduction to Clouds 2. Virtual Clouds based on physical Models 1. Generating Clouds 2. Rendering Clouds using Volume Rendering
More information