A 3D VOF model in cylindrical coordinates
|
|
- Joanna Summers
- 5 years ago
- Views:
Transcription
1 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 greatly, and they have been widely used to model multiphase phenomena associated with droplets, bubbles and sprays. Nearly all such methods, however, have been developed in Cartesian coordinates. Yet, there are many interesting multiphase flow problems that must be modeled in three-dimensional cylindrical coordinates; examples include the instability of liquid films inside or outside of a rotating cylinder, and the instability of immiscible interfacial flows between rotating concentric cylinders. There are other problems that can be solved on a Cartesian mesh, but for which a cylindrical mesh is a more suitable choice which leads to better accuracy at the same resolution; for example, simulating normal droplet impact onto a solid surface on a Cartesian mesh can produce numerical errors such as the asymmetric spread of a droplet that should remain axisymmetric[1]. To our knowledge, the only VOF model that has been developed in cylindrical coordinates is that of Chen et. al.[2] in which the volume fractions are advected using an algebraic method, i.e. an upwind/downwind scheme is used to calculate flux volumes. Such a VOF model leads to interface diffusion, and to unphysical volume fractions that need to be corrected. Our objective was to develop an accurate two-phase flow model in cylindrical geometry, with a focus on the accurate implementation of three aspects : interface tracking, the modeling of flows with high density ratios, and the modeling of surface tension. Here, we will briefly describe the methodology while focusing on the three above aspects. The incompressible Navier Stokes equations are discretized according to a colocated finite volume scheme. An explicit projection method is employed for time discretization. A one-phase approximation of the Navier- Stokes equations is used to model the two phases. Volume fractions f, are advected using a piecewise linear interface calculation (PLIC) method. In the PLIC method in three-dimensional domains, the interface is approximated by planes in interfacial cells. To define a plane, which can be written in the form of n.x = ρ, the plane normal n, and the plane constant, ρ, need to be determined. The normal is calculated from a height function (HF) which we describe when
2 we discuss the curvature calculation. The plane constant ρ is calculated so that the truncated volume, V truncated, is equal to the volume of the fluid in each cell, V fluid. Instead of solving for ρ directly, an iterative method is used to to find ρ so that V truncated (ρ) = V fluid. A geometrical algorithm has been developed to calculate the V truncated for a given plane of n and ρ, based on the algorithm first developed by Rider and Kothe [3] for determining the linear interface segments in a two dimensional Cartesian mesh. The reconstructed interface together with the velocity field are then used to estimate the volume fractions in the next time step. The volume fraction advection equation is solved using an operator split method. To model flows with high density ratios accurately, the mass (volume fraction) advection and momentum advection should be tightly coupled [4, 5]. Therefore, the Navier-Stokes equations are solved in a fully conservative form, where mass fluxes are calculated from the volume fraction advection algorithm. To model surface tension accurately two aspects are important: the surface tension model and the accuracy of curvature calculation. Here, we focus on the curvature calculation. Recent studies on the various methods of curvature calculation using volume fractions, have shown that Height Function (HF) method and the PROST (parabolic reconstruction of surface tension) method [6] are superior to other methods in terms of the magnitude of error and the order of convergence rate [7, 8]. The PROST method, however, is computationally more expensive, and the implementation of the method is not straightforward [8]. Therefore, the HF method was implemented to calculate curvatures. For a surface of the form R = h(θ, z), curvature, κ, can be derived as κ = 1 R 2 θ + R 2 + R 2 zr 23 (R2 θ RR zz + R 3 R zz 2RR θ R θz R z +RR 2 z R θθ R 2 R 2 z 2R2 θ + RR θθ R 2 ). R in the cell (i,j,k) is estimated by adding the volumes of one of the fluids in a column of cells in the R direction. The derivatives of R are approximated using second order finite difference schemes with R values calculated in a S D 3 3 stencil around cell (i,j,k), where S D is the stencil depth which is determined so that all of the cells in the first and last row are inside or outside of the heavy fluid. Curvature for surfaces of the forms Θ = h(r, z) and Z = h(r, θ) can be calculated similarly. Selection of the direction of the stencil (between R = h(θ, z), Θ = h(r, z), or Z = h(r, θ)) is based on the stencil depth, SD: the direction with the minimum SD leads to the most accurate estimation of curvature. We have completed the development of the model, and are beginning to validate our model. Here, we present the results for one such test. A simple but stringent test case for validating a surface tension calculation is the static drop test, where the pressure jump across the interface, and the magnitude of the induced currents are indicative of the accuracy of the surface tension implementation. The test was performed to compare the curvature
3 calculation using two methods: the HF method and a method based on calculating κ = ˆn from a convolved volume fraction field (CF). A continuum surface force model was employed for both cases. The test parameters correspond to the case tested by Gerlach et. al. [8] at high density ratios. Figure 1 shows the computational mesh. Figure 2 shows the pressure and velocity fields after one time step ( t = 1 5 ). The results show that HF method leads to a much smoother pressure field, and much smaller parasitic velocities. We propose to present the model and the results of a series validation tests, with an emphasis on those aspects that are particular to the cylindrical coordinate system. References [1] S. B. Johnson and J. P. Delplanque, Analyzing the results of droplet impact simulation by volume of fluid methods: Physics or numerical errors?, in ILASS Americas 18th Annual Conference on Liquid Atomization and Spray Systems, Irvine, CA, May 25. [2] L. Chen, S. V. Garimella, J. A. Reizes, and E. Leonardie, The development of a bubble rising in a viscous liquid, Journal of Fluid Mechanics, vol. 387, p. 61, [3] W. J. Rider and D. B. Kothe, Reconstructing volume tracking, Journal of Computational Physics, vol. 141, no. 2, p. 112, [4] M. Rudman, A volume-tracking method for incompressible multifluid flows with large density variations, International Journal for Numerical Methods in Fluids, vol. 28, no. 2, p. 357, [5] M. Bussmann, D. B. Kothe, and J. M. Sicilian, Modeling high density ratio incompressible interfacial flows, in Proceedings of the Sixth ASME/JSME Joint Fluids Engineering Conference, July 14-18, Montreal, Canada, 22. [6] Y. Renardy and M. Renardy, Prost: A parabolic reconstruction of surface tension for the volume-of-fluid method, Journal of Computational Physics, vol. 183, pp , DEC [7] S. J. Cummins, M. M. Francois, and D. B. Kothe, Estimating curvature from volume fractions, Computers and Structures, vol. 83, no. 6-7, p. 425, 25. [8] D. Gerlach, G. Tomar, G. Biswas, and F. Durst, Comparison of volume-of-fluid methods for surface tension-dominant two-phase flows, International Journal of Heat and Mass Transfer, vol. 49, pp , FEB 26.
4 Figure 1: The mesh used to test the static drop. The number of cells per radius is 16.
5 (a) (b) Figure 2: Comparison of pressure contours and velocity vectors for an static drop test using (a) HF method, and (b) CF method, after one time step ( t = 1 5 ). The radius of the drop is 2, the surface tension coefficient is 23.61, and the density ratio is 1 3. The maximum velocity is when using CF method and when using HF method.
Drop Impact Simulation with a Velocity-Dependent Contact Angle
ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Drop Impact Simulation with a Velocity-Dependent Contact Angle S. Afkhami and M. Bussmann Department
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
More 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 informationMETHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS
XVIII International Conference on Water Resources CMWR 2010 J. Carrera (Ed) CIMNE, Barcelona 2010 METHOD OF LEAST SQUARES FOR COMPUTING NORMALS AND CURVATURES FROM 2D VOLUME FRACTIONS Poorya A. Ferdowsi
More informationVolume Tracking on Adaptively Refined Grids with Curvature Based Refinement
Volume Tracking on Adaptively Refined Grids with Curvature Based Refinement Mayank Malik, Markus Bussmann Department of Mechanical & Industrial Engineering, University of Toronto mayank.malik@utoronto.ca,
More informationImplicit surface tension model for stimulation of interfacial flows
surface tension model for stimulation of interfacial flows Vinh The University of Massachusetts Dartmouth March 3rd, 2011 About my project Project Advisor Dr. Mehdi Raessi Department of Mechanical Engineering
More informationLeast Square Curvature Calculation Method for VOF Schemes
ILASS Americas 28th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2016 Least Square Curvature Calculation Method for VOF Schemes E. Cauble, M. Owkes Department of Mechanical
More informationHeight functions for applying contact angles to 3D VOF simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2009; 61:827 847 Published online 8 December 2008 in Wiley InterScience (www.interscience.wiley.com)..1974 Height functions
More informationAvailable online at ScienceDirect. Procedia IUTAM 15 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 15 (2015 ) 270 277 IUTAM Symposium on Multiphase flows with phase change: challenges and opportunities, Hyderabad, India (December
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More informationMoving Interface Problems: Methods & Applications Tutorial Lecture II
Moving Interface Problems: Methods & Applications Tutorial Lecture II Grétar Tryggvason Worcester Polytechnic Institute Moving Interface Problems and Applications in Fluid Dynamics Singapore National University,
More informationBubble Dynamics using Free Surfaces in a VOF framework
Bubble Dynamics using Free Surfaces in a VOF framework Sevilla Meeting: Numerical Challenges in two-phase flows 27 October 2014 Léon MALAN Supervisor: Prof. Stéphane ZALESKI 1 Contents Introduction: Simulation
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 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 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 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 informationAn accurate adaptive solver for surface-tension-driven interfacial flows
An accurate adaptive solver for surface-tension-driven interfacial flows Stéphane Popinet To cite this version: Stéphane Popinet. An accurate adaptive solver for surface-tension-driven interfacial flows.
More informationRobust Simulation of Sparsely Sampled Thin Features in SPH-Based Free Surface Flows
Copyright of figures and other materials in the paper belong to original authors. Robust Simulation of Sparsely Sampled Thin Features in SPH-Based Free Surface Flows Xiaowei He et al. ACM SIGGRAPH 2015
More 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 informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More 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 informationPROST: A Parabolic Reconstruction of Surface Tension for the Volume-of-Fluid Method
Journal of Computational Physics 183, 400 421 (2002) doi:10.1006/jcph.2002.7190 PROST: A Parabolic Reconstruction of Surface Tension for the Volume-of-Fluid Method Yuriko Renardy 1 and Michael Renardy
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationIsogeometric Analysis of Fluid-Structure Interaction
Isogeometric Analysis of Fluid-Structure Interaction Y. Bazilevs, V.M. Calo, T.J.R. Hughes Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA e-mail: {bazily,victor,hughes}@ices.utexas.edu
More informationTutorial 17. Using the Mixture and Eulerian Multiphase Models
Tutorial 17. Using the Mixture and Eulerian Multiphase Models Introduction: This tutorial examines the flow of water and air in a tee junction. First you will solve the problem using the less computationally-intensive
More informationAn accurate numerical method for computing surface tension forces in CFD codes
1 An accurate numerical method for computing surface tension forces in CFD codes Numerical experiments with surface tension M. Coquerelle 1 S. Glockner 1 1 Numerical Fluid Mechanics (MFN), département
More informationVOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 67 69 (997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS MURRAY RUDMAN Commonwealth Scientific and Industrial Research Organisation
More informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
More informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationA Finite Particle Approach for the Simulation of Multiphase Flows
ILASS-Americas 29th Annual Conference on Liquid Atomization and Spray Systems, Atlanta, GA, May 2017 A Finite Particle Approach for the Simulation of Multiphase Flows E.A. Wenzel and S.C. Garrick Department
More informationAccurate interface normal and curvature estimates on three-dimensional unstructured non-convex polyhedral meshes
Center for Turbulence Research Annual Research Briefs 2014 53 Accurate interface rmal and curvature estimates on three-dimensional unstructured n-convex polyhedral meshes By C. B. Ivey AND P. Moin 1. Motivation
More informationThe Surface Tension Challenge in Air-Water Interfaces using the Volume-of-Fluid Method
The Surface Tension Challenge in Air-Water Interfaces using the Volume-of-Fluid Method Konstantinos Politis a,d, Patrick Queutey b,d, Michel Visonneau c,d a. konstantinos.politis@ec-nantes.fr b. patrick.queutey@ec-nantes.fr
More informationA Sharp Interface Cartesian Grid Method for Simulating Flows with Complex Moving Boundaries
Journal of Computational Physics 174, 345 380 (2001) doi:10.1006/jcph.2001.6916, available online at http://www.idealibrary.com on A Sharp Interface Cartesian Grid Method for Simulating Flows with Complex
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 informationOverview of Traditional Surface Tracking Methods
Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx
More informationDroplet collisions using a Level Set method: comparisons between simulation and experiments
Computational Methods in Multiphase Flow III 63 Droplet collisions using a Level Set method: comparisons between simulation and experiments S. Tanguy, T. Ménard & A. Berlemont CNRS-UMR6614-CORIA, Rouen
More 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 informationApplication of coupled solver technology to CFD modeling of multiphase flows with CFX
Application of coupled solver technology to CFD modeling of multiphase flows with CFX A. Bums, A. Splawski, S. Lo & C. Guetari AEA Technology, 2000 Oxford Drive, Suite 61 0, Bethel Park, PA 151 02, USA
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 informationNUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE
NUMERICAL MODELING STUDY FOR FLOW PATTERN CHANGES INDUCED BY SINGLE GROYNE Jungseok Ho 1, Hong Koo Yeo 2, Julie Coonrod 3, and Won-Sik Ahn 4 1 Research Assistant Professor, Dept. of Civil Engineering,
More informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
More informationMIKE 21 & MIKE 3 FLOW MODEL FM. Transport Module. Short Description
MIKE 21 & MIKE 3 FLOW MODEL FM Short Description MIKE213_TR_FM_Short_Description.docx/AJS/EBR/2011Short_Descriptions.lsm//2011-06-17 MIKE 21 & MIKE 3 FLOW MODEL FM Agern Allé 5 DK-2970 Hørsholm Denmark
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 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 informationGradient Free Design of Microfluidic Structures on a GPU Cluster
Gradient Free Design of Microfluidic Structures on a GPU Cluster Austen Duffy - Florida State University SIAM Conference on Computational Science and Engineering March 2, 2011 Acknowledgements This work
More informationAnimation of Bubbles in Liquid
EUROGRAPHICS 2003 / P. Brunet and D. Fellner (Guest Editors) Volume 22 (2003), Number 3 Animation of Bubbles in Liquid Jeong-Mo Hong and Chang-Hun Kim Department of Computer Science, Korea University Abstract
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 informationMIKE 21 & MIKE 3 Flow Model FM. Transport Module. Short Description
MIKE 21 & MIKE 3 Flow Model FM Transport Module Short Description DHI headquarters Agern Allé 5 DK-2970 Hørsholm Denmark +45 4516 9200 Telephone +45 4516 9333 Support +45 4516 9292 Telefax mike@dhigroup.com
More informationA higher-order finite volume method with collocated grid arrangement for incompressible flows
Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More informationA TSTT integrated FronTier code and its applications in computational fluid physics
Institute of Physics Publishing Journal of Physics: Conference Series 16 (2005) 471 475 doi:10.1088/1742-6596/16/1/064 SciDAC 2005 A TSTT integrated FronTier code and its applications in computational
More 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 informationReconstructing Volume Tracking 1
JOURNAL OF COMPUTATIONAL PHYSICS 141, 112 152 (1998) ARTICLE NO. CP985906 Reconstructing Volume Tracking 1 William J. Rider and Douglas B. Kothe Applied Theoretical and Computational Physics Division,
More informationFEM techniques for interfacial flows
FEM techniques for interfacial flows How to avoid the explicit reconstruction of interfaces Stefan Turek, Shu-Ren Hysing (ture@featflow.de) Institute for Applied Mathematics University of Dortmund Int.
More informationCFD modelling of thickened tailings Final project report
26.11.2018 RESEM Remote sensing supporting surveillance and operation of mines CFD modelling of thickened tailings Final project report Lic.Sc.(Tech.) Reeta Tolonen and Docent Esa Muurinen University of
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 informationModelling of Levitation Melting using a Fixed Mesh Method
International Scientific Colloquium Modelling for Electromagnetic Processing Hannover, October 27-29, 2008 Modelling of Levitation Melting using a Fixed Mesh Method D. Hectors, E. Toorman, K. Van Reusel
More informationmidas NFX 2017R1 Release Note
Total Solution for True Analysis-driven Design midas NFX 2017R1 Release Note 1 midas NFX R E L E A S E N O T E 2 0 1 7 R 1 Major Improvements Midas NFX is an integrated finite element analysis program
More informationSimulation of Droplet Impingement on a Solid Surface by the Level Set Method
Simulation of Droplet Impingement on a Solid Surface by the Level Set Method Junling Hu *1, Ruoxu Jia 1, Kai-tak Wan 2, Xingguo Xiong 3 1 Department of Mechanical Engineering, University of Bridgeport,
More informationMultiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM3
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol., pp.491-497 (11 ARTICLE Multiphase Fluid Simulations on a Multiple GPGPU PC Using Unsplit Time Integration VSIAM3 Akio IKEBATA 1,, Yoshihiko MURAOKA 1 and
More informationMoment of fluid method for multimaterial flows
Moment of fluid method for multimaterial flows M. Jemison, M. Sussman, M. Shashkov Department of Mathematics, Florida State University October 1, 2012 Recent Motivating Work Volume preserving Moment-of-Fluid
More informationChapter 2 Water and Bubbles
Chapter 2 Water and Bubbles 2.1 Animation of Bubbles in Liquid Abstract This section introduces a new fluid animation technique in which liquid and gas interact with each other, using the example of bubbles
More informationANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step
ANSYS AIM Tutorial Turbulent Flow Over a Backward Facing Step Author(s): Sebastian Vecchi, ANSYS Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Governing Equation Start-Up Geometry
More informationINTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS
INTERACTION BETWEEN TURBULENT DYNAMICAL PROCESSES AND STATISTICS IN DEFORMED AIR-LIQUID INTERFACES, VIA DNS Yoshinobu Yamamoto Department of Nuclear Engineering, Kyoto University Yoshida Sakyo Kyoto, 66-85,
More informationSolved with COMSOL Multiphysics 4.2
Laminar Static Mixer Introduction In static mixers, also called motionless or in-line mixers, a fluid is pumped through a pipe containing stationary blades. This mixing technique is particularly well suited
More informationA Survey of Interface Tracking Methods in Multi-phase Fluid Visualization
A Survey of Interface Tracking Methods in Multi-phase Fluid Visualization Fang Chen 1 and Hans Hagen 2 1 University of Kaiserslautern AG Computer Graphics and HCI Group, TU Kaiserslautern, Germany chen@cs.uni-kl.de
More informationTwo-dimensional Numerical Simulation of Multiphase Flows in MATLAB
Two-dimensional Numerical Simulation of Multiphase Flows in MATLAB Saadi Daftari A Thesis in The Department of Mechanical and Industrial Engineering Presented in Partial Fulfillment of the Requirements
More informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationAn Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries
Journal of Computational Physics 156, 209 240 (1999) Article ID jcph.1999.6356, available online at http://www.idealibrary.com on An Accurate Cartesian Grid Method for Viscous Incompressible Flows with
More informationUsing the Eulerian Multiphase Model for Granular Flow
Tutorial 21. Using the Eulerian Multiphase Model for Granular Flow Introduction Mixing tanks are used to maintain solid particles or droplets of heavy fluids in suspension. Mixing may be required to enhance
More 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 informationDeforming fluid domains within the finite element method
Deforming fluid domains within the finite element method Five mesh-based tracking methods in comparison S. Elgeti H. Sauerland Aachen Institute for Advanced Study in Computational Engineering Science Financial
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 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 informationStrömningslära Fluid Dynamics. Computer laboratories using COMSOL v4.4
UMEÅ UNIVERSITY Department of Physics Claude Dion Olexii Iukhymenko May 15, 2015 Strömningslära Fluid Dynamics (5FY144) Computer laboratories using COMSOL v4.4!! Report requirements Computer labs must
More informationConvective Schemes for Capturing Interfaces of Free-Surface Flows on Unstructured Grids
Abstract Convective Schemes for Capturing Interfaces of Free-Surface Flows on Unstructured Grids M. Darwish and F. Moukalled American University of Beirut, Faculty of Engineering & Architecture, Mechanical
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 informationTutorial: Hydrodynamics of Bubble Column Reactors
Tutorial: Introduction The purpose of this tutorial is to provide guidelines and recommendations for solving a gas-liquid bubble column problem using the multiphase mixture model, including advice on solver
More informationThe Level Set Method THE LEVEL SET METHOD THE LEVEL SET METHOD 203
The Level Set Method Fluid flow with moving interfaces or boundaries occur in a number of different applications, such as fluid-structure interaction, multiphase flows, and flexible membranes moving in
More informationCFD Post-Processing of Rampressor Rotor Compressor
Gas Turbine Industrial Fellowship Program 2006 CFD Post-Processing of Rampressor Rotor Compressor Curtis Memory, Brigham Young niversity Ramgen Power Systems Mentor: Rob Steele I. Introduction Recent movements
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationMultiphase Flow Developments in ANSYS CFX-12
Multiphase Flow Developments in ANSYS CFX-12 Thomas Svensson Medeso 2008 ANSYS, Inc. All rights reserved. 1 ANSYS, Inc. Proprietary Outline Euler-Euler Wall Boiling Model Non-Drag Forces Euler-Lagrange
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationOn the thickness of discontinuities computed by THINC and RK schemes
The 9th Computational Fluid Dynamics Symposium B7- On the thickness of discontinuities computed by THINC and RK schemes Taku Nonomura, ISAS, JAXA, Sagamihara, Kanagawa, Japan, E-mail:nonomura@flab.isas.jaxa.jp
More informationThis tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following:
Tutorial 22. Modeling Solidification Introduction This tutorial illustrates how to set up and solve a problem involving solidification. This tutorial will demonstrate how to do the following: Define a
More informationStreamlining Aircraft Icing Simulations. D. Snyder, M. Elmore
Streamlining Aircraft Icing Simulations D. Snyder, M. Elmore Industry Analysis Needs / Trends Fidelity Aircraft Ice Protection Systems-Level Modeling Optimization Background Ice accretion can critically
More informationSELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND
Student Submission for the 5 th OpenFOAM User Conference 2017, Wiesbaden - Germany: SELECTIVE ALGEBRAIC MULTIGRID IN FOAM-EXTEND TESSA UROIĆ Faculty of Mechanical Engineering and Naval Architecture, Ivana
More informationMarker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 1, pp. 1-41 Commun. Comput. Phys. January 2009 REVIEW ARTICLE Marker-Based, 3-D Adaptive Cartesian Grid Method for Multiphase Flow Around Irregular Geometries
More informationSimulation of droplet impact on a solid surface using the level contour reconstruction method
Journal of Mechanical Science and Technology 23 (2009) 2434~2443 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 10.1007/s12206-009-0621-z Simulation of droplet
More informationMicrowell Mixing with Surface Tension
Microwell Mixing with Surface Tension Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract For many applications in the pharmaceutical
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 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 informationA Comparative CFD Analysis of a Journal Bearing with a Microgroove on the Shaft & Journal
Proceedings of International Conference on Innovation & Research in Technology for Sustainable Development (ICIRT 2012), 01-03 November 2012 182 A Comparative CFD Analysis of a Journal Bearing with a Microgroove
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 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 informationProgram: Advanced Certificate Program
Program: Advanced Certificate Program Course: CFD-Vehicle Aerodynamics Directorate of Training and Lifelong Learning #470-P, Peenya Industrial Area, 4th Phase Peenya, Bengaluru 560 058 www.msruas.ac.in
More informationConvection Cooling of Circuit Boards 3D Natural Convection
Convection Cooling of Circuit Boards 3D Natural Convection Introduction This example models the air cooling of circuit boards populated with multiple integrated circuits (ICs), which act as heat sources.
More informationMultiphysics simulations of nuclear reactors and more
Multiphysics simulations of nuclear reactors and more Gothenburg Region OpenFOAM User Group Meeting Klas Jareteg klasjareteg@chalmersse Division of Nuclear Engineering Department of Applied Physics Chalmers
More information