A Generalized Adaptive Collision Mesh for Multiple Injector Orifices
|
|
- Mary Brown
- 5 years ago
- Views:
Transcription
1 A Generalized Adaptive Collision Mesh for Multiple Injector Orifices Shuhai Hou, Sasanka Are, David P. Schmidt University of Massachusetts-Amherst ABSTRACT An algorithm for creating a generalized adaptive collision mesh has been proposed. Contrary to the common practice of using the fixed gas phase mesh to calculate droplet collision incidence, the proposed algorithm generates a collision mesh independent of the gas phase mesh. The collision mesh adaptively encompasses all the droplets in the computational domain at each time step. Schmidt and Rutland used such a mesh, but it was limited to a single spray. A new, more general algorithm has been developed. Calculations of multi-hole injectors have been performed to show that the mesh is capable of producing higher spatial resolution than the fixed gas phase mesh without incurring significant CPU cost. At the same time, the mesh algorithm maintains an adequate statistical sample of parcels in the collision cells. The continual random rotation of the mesh prevents collision from causing grid dependency in the spray shape. In addition to common multi-hole diesel injectors, the mesh algorithm is especially useful in simulating MEMS nozzles that are characterized by many closely spaced orifices. Due to the extremely small spacing between the orifices, inter-spray collisions are expected to play a significant role on spray behavior, requiring accurate calculation of droplet collision. INTRODUCTION It is intuitive to use the fixed gas phase mesh to calculate droplet collision. However, this approach suffers from two important weaknesses. First, the gas phase mesh is usually too coarse to evaluate collisions in the vicinity of the injector, where the droplet number density and collision frequency are very high. The gas phase mesh tends to produce large spatial discretization errors and under-resolves this dense spray area, aggravating grid dependency. Secondly, grid dependency is particularly severe when using a Cartesian mesh, as shown in Figure 1. The computed spray shape can be turned into a four-leafed clover by the grid dependency [1]. Since only parcels located in the same cell are allowed to collide, and the incidence of collision is proportional to the relative velocity between the pair parcels, the parcels whose trajectories are 90 degrees apart are most likely to collide because they have the greatest relative velocity. Parcels very close to each other, but in different quadrants, do not collide because their trajectories are almost aligned with the cell boundaries and more likely not in the same cell. Figure 1 A solid-cone spray directed towards the viewer. Figure 2 A cylindrical collision mesh superimposed on a spray. A Calculation made on a Cartesian mesh. Cartesian mesh gas phase mesh is partially cut away for visibility. This type of mesh has been used in past modeling of a single spray [1]. GENERALIZED ADAPTIVE COLLISION MESH In addressing the above shortcomings, Schmidt and Rutland took advantage of the fact that droplet collisions need not be directly dependent on the gas phase mesh [1]. They proposed a method for creating a collision mesh independent of the gas
2 phase. An example is shown in Figure 2. The collision mesh was based on a polar coordinate system aligned with the injection axis. The mesh would automatically form around an injection with any orientation. The angular orientation was chosen randomly every time step to produce a stochastic mesh. This method significantly improved spatial resolution and reduced grid dependency. It was successfully applied to spray bomb and engine simulations. However, this method is only applicable for a single spray because for multiple orifice nozzle sprays there is no single injection axis. Inspired by Schmidt and Rutland s ideas, an algorithm for creating a generalized, adaptive collision mesh is proposed. The collision mesh is independent of the gas phase mesh and created automatically at each time step. The new method can be used in multiple sprays with any configuration of orifices. It is based on the following steps: 1. Rotate the x, y, z axes by three random angles to generate a new coordinate system. The collision mesh will be oriented with this random coordinate system. This step is important to suppress the artifacts caused by Cartesian mesh. The randomness ensures that the collision mesh will be entirely different every time step. 2. Create a regular, three-dimensional, hexahedron encompassing all the parcels. 3. Check if the number of parcels in the hexahedron is below the maximum allowed number. If so, the mesh is completed. If not, continue to the next step. The maximum allowed number is used to make sure that a statistically adequate sample of parcels is used for the collision cells, and at the same time to control the resolution. 4. Halve the hexahedron in the x, y and z directions to make eight smaller hexahedra. 5. Check each hexahedron to see if the number of parcels in it is below the maximum allowed number. If so for all the interior hexahedra, the mesh refinement is completed. If not, repeat steps 4 and 5 for all the hexahedra in which the number of parcels exceeds the maximum allowed number. Though the creation of the collision mesh requires significant programming effort, its CPU cost is negligible. Typical runs show the collision meshing takes less than 0.1% of the total CPU time. Together with the NTC method, droplet collision represents less than 7% of the total CPU cost of the code when using about fifty thousand parcels. Figure 3 shows a picture of the collision mesh. It can be seen that the collision cells are smaller in the denser areas. Figure 4 shows its effectiveness in avoiding the grid dependency caused by a Cartesian mesh and traditional collision algorithm. Due to the stochastic nature of the mesh, it is entirely different every time step. Grid artifacts, which appear when the parcels are artificially forbidden to collide with near-neighbors, are prevented. Instead, the stochastic grouping means the parcels that are close to each other are likely to end up in the same collision cell during most time steps. The value of this approach is that the collision mesh is not tied to any one geometric configuration. Figures 5a and 5b show how the collision mesh can also be applied to a more typical spray configuration from a multi-hole diesel fuel injector. The mesh forms around any number of sprays in any orientation. Figure 3 A picture of the generalized adaptive collision mesh. Figure 4 Same case as in Figure 1, but with the generalized adaptive collision mesh.
3 Figure 5a and 5b Six sprays from a traditional orifice configuration are shown at the left. On the right, the collision mesh at one moment in time is shown superimposed around the spray. Part of the collision mesh and spray are cut away for clarity. NUMERICAL TESTS OF THE GENERALIZED ADAPTIVE COLLISION MESH In order to verify the characteristics of the generalized adaptive collision mesh, calculations were applied to a 5-hole MEMS nozzle sprays under typical diesel conditions. A three-dimensional, transient, CFD code named MOSES (MOdifiable Source Engine Simulation) was used for the present calculations. MOSES is currently under development by Convergent Thinking LLC in collaboration with the University of Massachusetts-Amherst. Numerical Schemes and Models A finite difference approach with a regular block-structured mesh is employed. With minor variations, the PISO method of Issa is used for solving the gas phase conservation equations for compressible flow [2]. All spatial derivatives are second-order central-difference, with the exception of the convection terms. For convection, since the main goal is to investigate the collision mesh, first-order upwinding is used for simplicity and efficiency. Turbulence is represented using the RNG k-ε model. All boundary conditions are no-slip walls. A semi-implicit Euler method is used for solving the ODE governing parcel trajectory. This method is exact with a constant ambient gas field. A predictor-corrector scheme is used to account for the variable gas velocity and density, with distinct time steps for individual parcels to control the spatial accuracy of the scheme. Parcels are sub-stepped with a time step approximately equal to the time required to traverse a cell. Tri-linear interpolation is used for gas to liquid coupling. Sirignano s interpolation method is used for liquid to gas coupling [3]. This method is second order accurate in space. The Enhanced Taylor Analogy Breakup model is used to model jet and droplet breakup [4]. For droplet turbulent dispersion, a time-correlated stochastic droplet dispersion model based on the k-ε model is used [5]. The NTC collision algorithm is used to calculate collision incidence [1]. The NTC algorithm is a popular method in calculating molecular collisions in gas dynamics. This algorithm has a computational cost linearly proportional to the total number of parcels in the cell. However, it had only been applied to the cases where the number of molecules in each parcel was constant. Schmidt and Rutland extended it to spray droplet collisions where each parcel contains varying number of droplets. Like O Rourke s collision algorithm, the NTC algorithm is first order accurate in time and second order accurate in space. Unlike O Rourke s collision algorithm, in which all the possible collision pairs are considered, the NTC algorithm only considers a randomly chosen subset of candidate pairs after scaling up the collision probability by the estimated maximum probability, reducing the computational cost. Thus more computational parcels can be used, allocating statistically enough parcels for each orifice in the MEMS nozzles. Detailed derivation and testing of the NTC algorithm can be found in reference [1]. O Rourke s model for collision outcomes is used to calculate the outcome of collision [6]. Post and Abraham have demonstrated the need for maintaining an adequate number of parcels per cell in collision calculation [7]. In the extreme, if there are less than two parcels in a cell, collision is totally suppressed. In order to observe this effect with the NTC algorithm, a spatially homogenous test case was run using the three-dimensional demonstration code of Garcia [8]. This is a freely-available program for demonstrating collision calculations. The number of parcels per cell was varied to see how this changed the predicted number of collisions over a short time interval. It was verified that the NTC collision algorithm is sufficiently accurate when at least five parcels are in a cell.
4 Computational Conditions Calculations were made on the 5-orifice MEMS nozzles as used by S. Baik, et al.[9]. Table 1 shows the computational conditions. Computational domain Gas phase mesh resolution mm mm Parcels used 100,000 Table 1 Computational conditions. RESULTS AND DISCUSSION The results and features of the generalized adaptive collision mesh are presented below. For the sake of clarity, three types of plots regarding the mesh are shown. One is Fraction of Cells versus Average Cell Size, where the Fraction of Cells refers to the percentage of the cells with size less than the indexed values. One is Fraction of Parcels versus Average Cell Size, where the Fraction of Parcels refers to the percentage of the parcels located in the cells with size less than the indexed values. One is Fraction of Parcels versus Number of Parcels in Cells, where Fraction of Parcels refers to the percentage of the parcels located in the cells in which the number of parcels is lager than the indexed values. The average cell size is calculated from the cube root of the cell volume. The term ASI refers to After Start of Injection. Figures 6, 7 and 8 show the results at 0.4ms ASI with the maximum allowed number of parcels (as used in the algorithm description in Step 3) equal to 160. It can be seen from Figure 6 that over 80% of the cells have average size less than 1mm, and over 50% less than 0.46mm. Also it is shown in Figure 7 that over 90% of the parcels are located in the cells with size less than 1mm, and almost 60% in less than 0.46mm. It is clear that the collision mesh produces better resolution than the gas phase mesh does, and thus is expected to produce higher spatial accuracy and to reduce grid dependency. Figure 8 verifies that the present approach maintains an adequate statistical sample size. Almost all parcels are located in the cells with more than five parcels in them, providing a minimum acceptable sample size for the stochastic collision calculation. Figure 6 Distribution of cells about cell size at 0.4ms ASI Figure 7 Distribution of parcels about cell size at 0.4ms ASI. Figures 9, 10 and 11 show the same plots for time at 3.0ms ASI. As we may expect, the spray gets more spread out and dilute, requiring fewer small cells. Figure 9 shows that over 55% of the cells have average size less than 1mm. Figure 10 shows that over 45% of the parcels are located in the cells with size less than 1mm. It is interesting to note that the distribution in Figure 11 differs little from that in Figure 8. In both cases, the vast majority of cells contain many more than five parcels, the minimum acceptable sample size. This observation suggests that the maximum allowed number of parcels could be decreased to get better resolution, still yet maintaining the minimum acceptable sample size in cells. Further calculations have shown that 100 is an optimal value of the maximum number of parcels for Step 3 of the algorithm.
5 Figure 8 Distribution of parcels about sample size at 0.4ms ASI. Figure 9 Distribution of cells about cell size at 3.0ms ASI. Figure 10 Distribution of parcels about cell size at 3.0ms ASI. Figure 11 Distribution of parcels about sample size at 3.0ms ASI CONCLUSIONS The proposed algorithm of a generalized adaptive collision mesh is demonstrated to achieve good spatial resolution, and at the same time assure statistically adequate sample size. The optimum value of the maximum number of parcels in a cell required by the algorithm is approximately 100. The computational cost was very modest; the liquid phase calculation consumed less then seven percent of the CPU time. The generalized adaptive collision mesh, together with the NTC collision algorithm, provides a useful technique to study the droplet collisions with multiple injector orifices. ACKNOWLEDGEMENTS This work is funded by Caterpillar Inc. and General Motors. Inc. The authors thank Prof. Franz X. Tanner for generously sharing the ETAB source code. They also thank Eric Pomraning and Keith Richards of Convergent Thinking LLC for the use of their gas phase solver. REFERENCES 1. David P. Schmidt, C. J. Rutland, A New Droplet Collision Algorithm, Journal of Computational Physics, v. 164, R. I. Issa, Solution of the Implicitly Discretised Fluid Flow Equations by Operator-Splitting, Journal of Computational Physics, v. 62, William A. Sirignano, Fluid Dynamics and Transport of Droplets and Sprays, p.p , Cambridge University Press, F. X. Tanner, Liquid Jet Atomization and Droplet Breakup Modeling of Non-Evaporating Diesel Fuel Sprays, SAE 1997 Transactions: Journal of Engines, Vol. 106, Sec. 3, p.p , Q. Zhou, M. A. Leschziner, A Lagrangian Particle Dispersion Model Based On A Time-Correlated Stochastic Approach, FED-Vol. 121, Gas-Solid Flows, ASME 1991.
6 6. P. J. O Rourke, Collective Drop Effects on Vaporizing Liquid Sprays, Ph.D. thesis, Department of Mechanical and Aerospace Engineering, Princeton University, S. Post and J. Abraham, Modeling Collisions Among Drops with a Non-Uniform Spatial Distribution, ILASS America 2002, Madison, WI. 8. Alejandro L. Garcia, Numerical Methods for Physics, 2 nd ed. Prentice-Hall, Seunghyun Baik, James P. Blanchard, Michael L. Corradini, Development of Micro-Diesel Injector Nozzles via MEMS Technology and Effects on Spray Characteristics, Ph.D. thesis, Engine Research Center and Engineering Physics, University of Wisconsin Madison, 2001.
Preliminary 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 informationNumerical Spray Calibration Process. ANSYS Inc. Pune, India. Laz Foley. ANSYS Inc. Evanston, IL
ILASS Americas, 23 rd Annual Conference on Liquid Atomization and Spray Systems, Ventura, CA, May 2011 Numerical Spray Calibration Process Padmesh Mandloi *,Jayesh Mutyal, Pravin Rajeshirke ANSYS Inc.
More informationAir Assisted Atomization in Spiral Type Nozzles
ILASS Americas, 25 th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 2013 Air Assisted Atomization in Spiral Type Nozzles W. Kalata *, K. J. Brown, and R. J. Schick Spray
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 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 informationUrea Injection Simulation with Adaptive Mesh Refinement for Engine Aftertreatment
ILASS Americas 26th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 2014 Urea Injection Simulation with Adaptive Mesh Refinement for Engine Aftertreatment Scott A. Drennan
More informationClick to edit Master title style
Click to edit Master title style LES LES Applications for for Internal Internal Combustion Engines Engines David Gosman & Richard Johns CD-adapco, June 2011 Some Qs and As Why would we use LES calculations
More information1. TopMath-Workshop Iffeldorf/Osterseen. Liquid Sprays. Ayoub Hmaidi Zentrum Mathematik, TU MÜNCHEN
1. TopMath-Workshop Iffeldorf/Osterseen Liquid Sprays Ayoub Hmaidi Zentrum Mathematik, TU MÜNCHEN What are Liquid Sprays? Why are Sprays important? Sprays occur in a large number of applications: Engines
More informationDeveloping LES Models for IC Engine Simulations. June 14-15, 2017 Madison, WI
Developing LES Models for IC Engine Simulations June 14-15, 2017 Madison, WI 1 2 RANS vs LES Both approaches use the same equation: u i u i u j 1 P 1 u i t x x x x j i j T j The only difference is turbulent
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationApplying Solution-Adaptive Mesh Refinement in Engine Simulations
International Multidimensional Engine Modeling User's Group Meeting April 11, 2016, Detroit, Michigan Applying Solution-Adaptive Mesh Refinement in Engine Simulations Long Liang, Yue Wang, Anthony Shelburn,
More informationSIMULATION OF THE SPRAY COATING PROCESS USING A PNEUMATIC ATOMIZER
ILASS-Europe 22 Zaragoza 9 11 September 22 SIMULATION OF THE SPRAY COATING PROCESS USING A PNEUMATIC ATOMIZER Q. Ye, J. Domnick, E. Khalifa Fraunhofer-Institute for Manufacturing Engineering and Automation
More informationCFD Best Practice Guidelines: A process to understand CFD results and establish Simulation versus Reality
CFD Best Practice Guidelines: A process to understand CFD results and establish Simulation versus Reality Judd Kaiser ANSYS Inc. judd.kaiser@ansys.com 2005 ANSYS, Inc. 1 ANSYS, Inc. Proprietary Overview
More informationValidation of an Automatic Mesh Generation Technique in Engine Simulations
International Multidimensional Engine Modeling User's Group Meeting April,, Detroit, Michigan Validation of an Automatic Mesh Generation Technique in Engine s Abstract Long Liang, Anthony Shelburn, Cheng
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 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 informationNumerical study & validation of a complete
Numerical study & validation of a complete SCR system using 1D-3D (CFD) coupling Presenter: Ashish Joshi Manager, Indian Operations Convergent Science Presenting on behalf of: Scott Drennan Director of
More informationCOMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING
2015 WJTA-IMCA Conference and Expo November 2-4 New Orleans, Louisiana Paper COMPUTATIONAL FLUID DYNAMICS USED IN THE DESIGN OF WATERBLAST TOOLING J. Schneider StoneAge, Inc. Durango, Colorado, U.S.A.
More informationKEY STAR TECHNOLOGIES: DISPERSED MULTIPHASE FLOW AND LIQUID FILM MODELLING DAVID GOSMAN EXEC VP TECHNOLOGY, CD-adapco
KEY STAR TECHNOLOGIES: DISPERSED MULTIPHASE FLOW AND LIQUID FILM MODELLING DAVID GOSMAN EXEC VP TECHNOLOGY, CD-adapco INTRODUCTION KEY METHODOLOGIES AVAILABLE IN STAR-CCM+ AND STAR-CD 1. Lagrangian modelling
More informationMcNair Scholars Research Journal
McNair Scholars Research Journal Volume 2 Article 1 2015 Benchmarking of Computational Models against Experimental Data for Velocity Profile Effects on CFD Analysis of Adiabatic Film-Cooling Effectiveness
More informationPDF-based simulations of turbulent spray combustion in a constant-volume chamber under diesel-engine-like conditions
International Multidimensional Engine Modeling User s Group Meeting at the SAE Congress Detroit, MI 23 April 2012 PDF-based simulations of turbulent spray combustion in a constant-volume chamber under
More informationExpress Introductory Training in ANSYS Fluent Workshop 02 Using the Discrete Phase Model (DPM)
Express Introductory Training in ANSYS Fluent Workshop 02 Using the Discrete Phase Model (DPM) Dimitrios Sofialidis Technical Manager, SimTec Ltd. Mechanical Engineer, PhD PRACE Autumn School 2013 - Industry
More informationCoarse Mesh CFD: Trend Analysis In a Fraction of the Time
Coarse Mesh CFD: Trend Analysis In a Fraction of the Time Y. He, C. J. Rutland, Z. Nagel, R. P. Hessel, R. D. Reitz, D.E. Foster Engine Research Center, University of Wisconsin-Madison In an effort to
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
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 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 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 informationContents Contents Contents... 1 Abstract... 3 Nomenclature... 4 Index of Figures... 6 Index of Tables... 8 Introduction... 9 Theory...
Contents Contents Contents... 1 Abstract... 3 Nomenclature... 4 Index of Figures... 6 Index of Tables... 8 1. Introduction... 9 1.1 Overview... 9 1.2 Task... 10 2. Theory... 11 2.1 Continuity and Momentum
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
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 informationMulti-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech
Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated
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 informationDirections: 1) Delete this text box 2) Insert desired picture here
Directions: 1) Delete this text box 2) Insert desired picture here Multi-Disciplinary Applications using Overset Grid Technology in STAR-CCM+ CD-adapco Dmitry Pinaev, Frank Schäfer, Eberhard Schreck Outline
More informationMultigrid Solvers in CFD. David Emerson. Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK
Multigrid Solvers in CFD David Emerson Scientific Computing Department STFC Daresbury Laboratory Daresbury, Warrington, WA4 4AD, UK david.emerson@stfc.ac.uk 1 Outline Multigrid: general comments Incompressible
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationHydrocyclones and CFD
Design of Hydrocyclone Separation Equipment Using CFD Coupled with Optimization Tools David Schowalter 1, Rafiqul Khan 1, Therese Polito 2, and Tim Olson 3. (1) Fluent Inc., 10 Cavendish Court, Lebanon,
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 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 informationICE Roadmap Japanese STAR Conference. Richard Johns
ICE Roadmap Japanese STAR Conference Richard Johns Introduction Top-Level Roadmap STAR-CCM+ and Internal Combustion Engines Modeling Improvements and Research Support Sprays LES Chemistry Meshing Summary
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 informationDevelopment of a CFD methodology for fuel-air mixing and combustion modeling of GDI Engines
Development of a CFD methodology for fuel-air mixing and combustion modeling of GDI Engines T. Lucchini, G. D Errico, L. Cornolti, G. Montenegro, A. Onorati Politecnico di Milano, Dipartimento di Energia,
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
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 information*Sebastian Hensel, Kai Herrmann, Reiner Schulz and German Weisser
Paper templeate (for draft and final papers) Numerical analysis and statistical description of the primary breakup in fuel nozzles of large two stroke engines for the application in CFD engine simulations
More informationCold Flow Simulation Inside an SI Engine
Tutorial 12. Cold Flow Simulation Inside an SI Engine Introduction The purpose of this tutorial is to illustrate the case setup and solution of the two dimensional, four stroke spark ignition (SI) engine
More informationarxiv: v1 [physics.flu-dyn] 30 Jan 2018
An Analysis of the Convergence of Stochastic Lagrangian/Eulerian Spray Simulations arxiv:1801.10070v1 [physics.flu-dyn] 30 Jan 2018 Abstract David P. Schmidt a,, Frederick Bedford b a Mechanical and Industrial
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationWAVE PATTERNS, WAVE INDUCED FORCES AND MOMENTS FOR A GRAVITY BASED STRUCTURE PREDICTED USING CFD
Proceedings of the ASME 2011 30th International Conference on Ocean, Offshore and Arctic Engineering OMAE2011 June 19-24, 2011, Rotterdam, The Netherlands OMAE2011-49593 WAVE PATTERNS, WAVE INDUCED FORCES
More informationStratified Oil-Water Two-Phases Flow of Subsea Pipeline
Stratified Oil-Water Two-Phases Flow of Subsea Pipeline Adib Zulhilmi Mohd Alias, a, Jaswar Koto, a,b,*, Yasser Mohamed Ahmed, a and Abd Khair Junaidi, b a) Department of Aeronautics, Automotive and Ocean
More informationBest Practices Workshop: Overset Meshing
Best Practices Workshop: Overset Meshing Overview Introduction to Overset Meshes Range of Application Workflow Demonstrations and Best Practices What are Overset Meshes? Overset meshes are also known as
More informationErosion modeling to improve asset life prediction
Erosion modeling to improve asset life prediction Erosion caused by solid particles in process flows impacting on the surfaces of downhole equipment or pipe walls is a common cause of damage and wear in
More informationDevelopment of Adaptive Mesh Refinement Scheme and Conjugate Heat Transfer Model for Engine Simulations
Graduate Theses and Dissertations Graduate College 2009 Development of Adaptive Mesh Refinement Scheme and Conjugate Heat Transfer Model for Engine Simulations Qingluan Xue Iowa State University Follow
More informationMetafor FE Software. 2. Operator split. 4. Rezoning methods 5. Contact with friction
ALE simulations ua sus using Metafor eao 1. Introduction 2. Operator split 3. Convection schemes 4. Rezoning methods 5. Contact with friction 1 Introduction EULERIAN FORMALISM Undistorted mesh Ideal for
More informationA steady-state Eulerian-Lagrangian solver for non-reactive sprays
ICLASS 212, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 212 A steady-state Eulerian-Lagrangian solver for non-reactive sprays A.
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More informationINVESTIGATION OF HYDRAULIC PERFORMANCE OF A FLAP TYPE CHECK VALVE USING CFD AND EXPERIMENTAL TECHNIQUE
International Journal of Mechanical Engineering and Technology (IJMET) Volume 10, Issue 1, January 2019, pp. 409 413, Article ID: IJMET_10_01_042 Available online at http://www.ia aeme.com/ijmet/issues.asp?jtype=ijmet&vtype=
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 informationValidation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons
Validation of an Unstructured Overset Mesh Method for CFD Analysis of Store Separation D. Snyder presented by R. Fitzsimmons Stores Separation Introduction Flight Test Expensive, high-risk, sometimes catastrophic
More informationSimulation of In-Cylinder Flow Phenomena with ANSYS Piston Grid An Improved Meshing and Simulation Approach
Simulation of In-Cylinder Flow Phenomena with ANSYS Piston Grid An Improved Meshing and Simulation Approach Dipl.-Ing. (FH) Günther Lang, CFDnetwork Engineering Dipl.-Ing. Burkhard Lewerich, CFDnetwork
More informationreal world design & problems fluid flow engineering solving Computational Fluid Dynamics System
C F D 2 0 0 0 Computational Fluid Dynamics System solving real world engineering design & problems aerospace architecture automotive biomedical chemical processing electrical cooling environmental marine
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationThe Elimination of Correlation Errors in PIV Processing
9 th International Symposium on Applications of Laser Techniques to Fluid Mechanics, Lisbon, Portugal, July, 1998 The Elimination of Correlation Errors in PIV Processing Douglas P. Hart Massachusetts Institute
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationHealthy Buildings 2017 Europe July 2-5, 2017, Lublin, Poland
Healthy Buildings 2017 Europe July 2-5, 2017, Lublin, Poland Paper ID 0122 ISBN: 978-83-7947-232-1 Numerical Investigation of Transport and Deposition of Liquid Aerosol Particles in Indoor Environments
More informationRotating Moving Boundary Analysis Using ANSYS 5.7
Abstract Rotating Moving Boundary Analysis Using ANSYS 5.7 Qin Yin Fan CYBERNET SYSTEMS CO., LTD. Rich Lange ANSYS Inc. As subroutines in commercial software, APDL (ANSYS Parametric Design Language) provides
More informationModeling of a DaimlerChrysler Truck Engine using an Eulerian Spray Model
Modeling of a DaimlerChrysler Truck Engine using an Eulerian Spray Model C. Hasse, S. Vogel, N. Peters Institut für Technische Mechanik RWTH Aachen Templergraben 64 52056 Aachen Germany c.hasse@itm.rwth-aachen.de
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationCFD Simulation for Stratified Oil-Water Two-Phase Flow in a Horizontal Pipe
CFD Simulation for Stratified Oil-Water Two-Phase Flow in a Horizontal Pipe Adib Zulhilmi Mohd Alias, a, Jaswar Koto, a,b,* and Yasser Mohamed Ahmed, a a) Department of Aeronautics, Automotive and Ocean
More informationEstimation of Flow Field & Drag for Aerofoil Wing
Estimation of Flow Field & Drag for Aerofoil Wing Mahantesh. HM 1, Prof. Anand. SN 2 P.G. Student, Dept. of Mechanical Engineering, East Point College of Engineering, Bangalore, Karnataka, India 1 Associate
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 informationON THE NUMERICAL MODELING OF IMPINGING JET HEAT TRANSFER
ON THE NUMERICAL MODELING OF IMPINGING JET HEAT TRANSFER Mirko Bovo 1,2, Sassan Etemad 2 and Lars Davidson 1 1 Dept. of Applied Mechanics, Chalmers University of Technology, Gothenburg, Sweden 2 Powertrain
More informationSpray Modeling for Outwardly-Opening Hollow-Cone Injector
Spray Modeling for Outwardly-Opening Hollow-Cone Injector Jaeheon Sim KAUST 2016-01-0844 Published 04/05/2016 Jihad Badra Saudi Aramco Ahmed Elwardany and Hong Im KAUST CITATION: Sim, J., Badra, J., Elwardany,
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationForces on particles and bubbles Literature
Lecture 1: Forces on particles and bubbles In order to describe the motion of particles and bubbles in fluids it is necessary to consider all relevant forces, external and fluid dynamic forces. Starting
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 informationDRAGON SOLUTIONS FOR BENCHMARK BWR LATTICE CELL PROBLEMS
DRAGON SOLUTIONS FOR BENCHMARK BWR LATTICE CELL PROBLEMS R. Roy and G. Marleau Institut de Génie Nucléaire École Polytechnique de Montréal P.O.Box 6079, Station CV, Montreal, Canada roy@meca.polymtl.ca
More informationNumerical and experimental study of spray cooling of a heated metal surface
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Numerical and experimental study of spray cooling of a heated metal
More informationAnalysis of fluid-solid coupling vibration characteristics of probe based on ANSYS Workbench
Analysis of fluid-solid coupling vibration characteristics of probe based on ANSYS Workbench He Wang 1, a, Changzheng Zhao 1, b and Hongzhi Chen 1, c 1 Shandong University of Science and Technology, Qingdao
More informationLearn the various 3D interpolation methods available in GMS
v. 10.4 GMS 10.4 Tutorial Learn the various 3D interpolation methods available in GMS Objectives Explore the various 3D interpolation algorithms available in GMS, including IDW and kriging. Visualize the
More informationONE DIMENSIONAL (1D) SIMULATION TOOL: GT-POWER
CHAPTER 4 ONE DIMENSIONAL (1D) SIMULATION TOOL: GT-POWER 4.1 INTRODUCTION Combustion analysis and optimization of any reciprocating internal combustion engines is too complex and intricate activity. It
More informationSimulation of Turbulent Axisymmetric Waterjet Using Computational Fluid Dynamics (CFD)
Simulation of Turbulent Axisymmetric Waterjet Using Computational Fluid Dynamics (CFD) PhD. Eng. Nicolae MEDAN 1 1 Technical University Cluj-Napoca, North University Center Baia Mare, Nicolae.Medan@cunbm.utcluj.ro
More informationform are graphed in Cartesian coordinates, and are graphed in Cartesian coordinates.
Plot 3D Introduction Plot 3D graphs objects in three dimensions. It has five basic modes: 1. Cartesian mode, where surfaces defined by equations of the form are graphed in Cartesian coordinates, 2. cylindrical
More informationPerformance Metrics of a Parallel Three Dimensional Two-Phase DSMC Method for Particle-Laden Flows
Performance Metrics of a Parallel Three Dimensional Two-Phase DSMC Method for Particle-Laden Flows Benzi John* and M. Damodaran** Division of Thermal and Fluids Engineering, School of Mechanical and Aerospace
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 informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationLaser speckle based background oriented schlieren measurements in a fire backlayering front
Laser speckle based background oriented schlieren measurements in a fire backlayering front Philipp Bühlmann 1*, Alexander H. Meier 1, Martin Ehrensperger 1, Thomas Rösgen 1 1: ETH Zürich, Institute of
More informationParticle Tracing Module
Particle Tracing Module Particle Tracing Module Released with version 4.2a in October 2011 Add-on to COMSOL Multiphysics Combines with any COMSOL Multiphysics Module Particle Tracing Particle tracing
More informationA STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS Muthukumaran.C.K.
More informationUsing the Discrete Ordinates Radiation Model
Tutorial 6. Using the Discrete Ordinates Radiation Model Introduction This tutorial illustrates the set up and solution of flow and thermal modelling of a headlamp. The discrete ordinates (DO) radiation
More informationFlow and Heat Transfer in a Mixing Elbow
Flow and Heat Transfer in a Mixing Elbow Objectives The main objectives of the project are to learn (i) how to set up and perform flow simulations with heat transfer and mixing, (ii) post-processing and
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 informationNumerical Simulations of Granular Materials Flow around Obstacles: The role of the interstitial gas
Numerical Simulations of Granular Materials Flow around Obstacles: The role of the interstitial gas Avi Levy, Dept. Mech. Eng., Ben Gurion University, Beer Sheva, Israel. Mohamed Sayed, CHC, National Research
More informationFLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016
FLUENT Secondary flow in a teacup Author: John M. Cimbala, Penn State University Latest revision: 26 January 2016 Note: These instructions are based on an older version of FLUENT, and some of the instructions
More informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationA new Eulerian computational method for the propagation of short acoustic and electromagnetic pulses
A new Eulerian computational method for the propagation of short acoustic and electromagnetic pulses J. Steinhoff, M. Fan & L. Wang. Abstract A new method is described to compute short acoustic or electromagnetic
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 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 informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
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 information