Pressure-correction algorithm to solve Poisson system with constant coefficients for fast two-phase simulations
|
|
- Jordan Griffith
- 6 years ago
- Views:
Transcription
1 Center or Turbulence Research Annual Research Bries Pressure-correction algorithm to solve Poisson system with constant coeicients or ast two-phase simulations By D. Kim AND P. Moin 1. Motivation and objective It is a challenging problem to numerically simulate incompressible two-phase low owing to its complexity and high computational cost. The complexity o two-phase low are caused by the moving interace between the phases, where density and viscosity abruptly change. Various numerical methods such as the marker method, the volume-o-luid (VoF) method and the level set method, have been developed to handle arbitrarily shaped interaces in a ixed grid. Each o these phase tracking methods has its merits and demerits. In the present work, the level set method is used to track the phase interace. The advantages o the level set method is that the topology and movement o the interace are automatically captured. Thus, it is easy to implement numerically. The main drawback o the level set method is that the conservation o liquid volume is not guaranteed. Hybrid methods combining the level set method with VoF Sussman (2000) have been proposed to enhance the mass conservation property. However, the delocalization method has to be used to avoid local curvature error and it makes numerical scheme more complex. High computational cost is another big issue in turbulent simulations or two-phase low because it is diicult to apply ast solvers to the Poisson equations. For incompressible two-phase low, the continuity equation o divergence-ree velocity is used to orm the Poisson equation or the pressure-correction in the ractional step method. Since the phase interace is moving in time, the coeicients o the Poisson equation is also unction o time, while the coeicients is constant in case o single phase low where density is constant. Because the coeicients o the Poisson equation change in time, the gradientbased method such as the Biconjugate Gradient Stabilized Method is oten used instead o the much aster Multigrid method. This results in large extra computational cost since the Poisson solver is the main bottleneck in numerical simulations. In the present work, a new pressure-correction algorithm is developed to improve the computational speed o the Poisson solver in two-phase low simulations. In the pressurecorrection step, the mass conservation equation or VoF equation is used to obtain the Poisson equation with constant coeicients. Using this method, the mass conservation equation is ulilled at every time step instead o the divergence-ree continuity equation. Since the Poisson equation with constant coeicients is much easier to solve than variable coeicients, the computational speed is much improved. The canonical test problems are solved to demonstrate the accuracy and perormance o the method. 2. Governing equations The governing equations or incompressible, immiscible, two-phase low are as ollows:
2 368 D. Kim and P. Moin ψ t + (ψu) = 0, (2.1) u p + (uu) = t ρ + 1 ρ (µ( u + T u)) + g + T ρ, (2.2) where ψ, u, p, µ, and g are the volume o raction o luid 1, velocity, density, pressure, viscosity, and gravitational acceleration, respectively. T is the surace tension orce, which acts only on the phase interace (x ) and vanishes at other locations; it is described as ollows: T = σκδ(x x )(n), (2.3) where σ, κ, δ, and n are the surace tension coeicient, surace curvature, delta unction, and surace normal vector, respectively. The temporal evolution o the phase interace can be predicted by solving the level-set equation: G + u G = 0, (2.4) t where the iso-surace o G = 0 deines the location o the interace whereas G > 0 and G < 0 represent two dierent phase regions, respectively. In the computational domain, G is set to be a signed distance unction to the interace such that G = 1. (2.5) The interace normal vector n and the interace curvature κ can be calculated rom G as n = G G, (2.6) κ = n. (2.7) Assuming ρ and µ are constant within each luid, the density and viscosity are deined as a unction o G: ρ = ψρ 1 + (1 ψ)ρ 2, (2.8) µ = ψµ 1 + (1 ψ)µ 2, (2.9) ψ = H(G), (2.10) where the subscripts 1 and 2 denote properties in luid 1 and 2, respectively, and H is the Heaviside unction. 3. Numerical methodology In the present study, the Reined Level Set Grid (RLSG) method (Herrmann 2008) coupled with a Lagrangian particle-tracking method (Apte et al. 2003) is used to capture the phase interace. For the low solver, an unstructured-grid inite-volume DNS/LES solver named CDP is employed (Ham & Iaccarino 2004; Mahesh et al. 2004). The RLSG uses an auxiliary Cartesian grid (G-grid) whereas the CDP uses a base structured or unstructured
3 Pressure-correction algorithm or ast two-phase simulations 369 grid or luid low. Because separate grids are employed in the RLSG and CDP computations, an eective parallel coupling strategy is required or inormation exchange. The present numerical technique utilizes a coupling inrastructure, called CHIMPS (Alonso et al. 2006), to perorm the parallel data exchange between CDP and RLSG. In the present study, the velocities u n+1 i 3.1. Level set transport equation is used to advance Eq. 3.1 as G n+3/2 G n+1/2 + u n+1 G n+1/2 i = 0. (3.1) The velocities are obtained rom CDP and an iterative procedure is needed to solve both G n+3/2 and u n+1. This procedure will be explained in section 3.5. The numerical algorithm or the level set transport equation is based on the RLSG method (Herrmann 2008). We use a ith-order WENO scheme or Hamilton-Jacobi equations (Jiang & Peng 2000) in conjunction with a local Lax-Friedrichs entropy correction (Shu and Osher, 1989). Time integration is perormed using a third-order TVD Runge-Kutta method (Shu & Osher 1989) Reinitialization Due to numerical errors, the solution o Eq. 3.1 will not preserve G as a signed distance unction, i.e., G = 1. Thus, a reinitialization procedure has to be applied to reset G such that G = 1. There are several dierent methods to reinitialize G. Here, we solve the perturbed Hamilton-Jacobi PDE to a steady state (Chang et al. 1995). G τ + (V o V (τ))( P + κ)g = 0 (3.2) where V 0 is the total liquid volume at the initial state, and V (τ) is the liquid mass at pseudo time τ. κ is the local curvature and P is a positive stabilizing constant. Unortunately, this correction moves the G = 0 iso-surace, which is inevitable in a PDE-based reinitialization. However, we can conserve the global mass up to any accuracy depending on the criterion used to achieve the steady state solution. The global mass conservation should be always guaranteed by solving the Poisson equation deduced rom the continuity-constraint RLSG-CDP coupling The volume raction ψ cv o the cell volume V cv is calculated using a step unction developed by Van Der Pill et al. (2005) as, ψ cv = 1 H(G)dv. (3.3) V cv V cv The integration is done through CHIMPS. The density and viscosity at the cell center o the cell volume V cv are deined as ρ cv = ψ cv ρ 1 + (1 ψ cv )ρ 2, (3.4) µ cv = ψ cv µ 1 + (1 ψ cv )µ 2. (3.5)
4 370 D. Kim and P. Moin The properties at the ace o the cell volume are computed by averaging two neighboring properties at the cell center. The velocities, u i, in the level set equation (Eq. 3.1) are interpolated rom CDP grid using tri-linear interpolation Computation o low ields CDP is based on the ractional-step method or collocated variables on unstructured grids, which is especially suitable or large-scale parallel simulations o turbulent lows in complex geometries (Ham & Iaccarino 2004; Mahesh et al. 2004). The spatial discretization is perormed using a low dissipation second-order central-dierence scheme while time integration is perormed using a ully implicit Crank-Nicholson method. An upwind-biased scheme is employed to discretize nonlinear terms only at interace cells (Kim et al. 2007). The pressure gradient and surace tension orces are predicted by a balanced orce algorithm, which guarantees a discrete balance between them (Francois et al. 2006; Herrmann 2008). In the ollowing, the new pressure-correction algorithm is presented, which satisies the mass conservation equation (Eq. 2.1). The numerical algorithm is ψ n+3/2 cv ψcv n+1/2 + ψn+1/2 cv u n+1 = 0, (3.6) V cv u un + u n+1/2 u n+1/2 + u n+1/2 i,nbr A = V cv g + F viscous, (3.7) 2 ρ n+1/2 cv u n+1 ρn+1/2 cv u = pn+1/2 + T n+1/2. (3.8) where A is the ace area, u the ace normal velocity, g the gravity, T n+1/2 the surace tension orce. The viscous term F viscous is ully implemented and solved in the lux orm with the viscosity at the cell ace calculated by the harmonic mean o two neighboring viscosities at the cell volumes. Ater solving Eq. 3.7 to obtain u, the cell ace normal velocities u are calculated, u = 1 2 (u + u i,nbr)n i,. (3.9) Taking the divergence o Eq. 3.8 results in the ollowing Poisson equation, 2 p n+1/2 x 2 = 1 i ( ρn+1/2 cv u n+1 ρn+1/2 cv u ) + T n+1/2. (3.10) By using Eq. 3.4 and Eq. 3.6, the mass lux terms in the right hand side become ρ n+1/2 cv u n+1 = (ρ 1 ρ 2 ) ψ n+1/2 cv u n+1 = (ρ 1 ρ 2 ) ψn+3/2 cv ψcv n+1/2, (3.11) ρ n+1/2 cv u = (ρ 1 ρ 2 )u ψ n+1/2 cv + (ρ 2 + (ρ 1 ρ 2 )ψ n+1/2 cv ) u
5 Pressure-correction algorithm or ast two-phase simulations 371 = (ρ 1 ρ 2 ) ψ cv ψn+1/2 cv + ρ n+1/2 u cv, (3.12) where ψcv is ψcv = 1 H(G )dv. (3.13) V cv V cv G is computed rom the level set equation as G G n+1/2 + u i Then, Eq is written as and 2 p n+1/2 cv x 2 i p n+1/2 n = 1 A = 1 {(ρ 1 ρ 2 ) ψn+3/2 cv ψcv {(ρ 1 ρ 2 ) ψn+3/2 cv ψ cv G n+1/2 = 0, (3.14) + ρ n+1/2 cv V cv + ρ n+1/2 cv in terms o the inite volume ormulation, where T n+1/2 The surace tension T n+1/2 at the ace is deined, T n+1/2 = σκ n+1/2 ( ) ψ n+1/2 u } n+1/2 T +, (3.15) u A } + T n+1/2 A (3.16) is the ace normal surace tension.. (3.17) Pressure can be solved rom Eq by using Multigrid method because it has constant coeicients. Then, the ace normal velocity components to satisy the mass conservation equation (Eq. 3.6) is obtained, where P is deined as u n+1 = u P, (3.18) P = 1 ρ n+1/2 ( p n+1/2 ) 1 ρ n+1/2 T n+1/2. (3.19) The pressure gradient at the cell center, P is calculated rom the pressure gradient at the cell ace, P, using the ace-area weighted least-squares method (Mahesh et al. 2004) by minimizing ǫ = (P n i, P ) 2 A. (3.20) The velocities at the cell center are inally corrected, u n+1 = u P. (3.21) In the present method, the velocities at the ace do not necessarily satisy the divergence ree equation ( u = 0). Instead, the mass conservation equation (3.6) is guaranteed by solving Poisson solver.
6 372 D. Kim and P. Moin Figure 1. Inviscid column at equilibrium Iteration procedure at each time step The iteration process at each time step is ollowed: 1. Intermediate velocities u i are obtained rom the momentum equation. 2. Estimate G n+3/2 by solving the level set transport equation using u i. 3. Integrate ψcv n+3/2, κ n+3/2 cv and transer to CDP using CHIMPS. 4. Calculate ρ n+3/2 cv, µ n+3/2 cv rom ψcv n+3/2. 5. The corrected velocities u n+1 i and pressure p n+1/2 are obtained by solving the Poisson equation. 6. Interpolate u n+1 i on the CDP-grid to G-grid. 7. Repeat the process 2 using updated u n+1 i until convergence. 4. Results 4.1. Equilibrium inviscid column with exact curvature The equilibrium inviscid column is considered in order to demonstrate the balanced orce algorithm (Francois et al. 2006). As explained in Francois et al. (2006), no spurious currents should be produced around a column at equilibrium when the balanced orce algorithm is used. Here, we applied the parameters used in Francois et al. (2006). A column o radius R = 2 is placed at the center o 8 8 domain (igure 1). The surace tension coeicient is σ = 73 and the theoretical pressure jump across the interace is p ex = The densities inside and outside the column are ρ 1 and ρ 2, respectively. The same equidistance Cartesian grid is used or the CDP-grid and G-grid, and the grid size is set to 0.2. Table 1. shows the errors in velocities and pressure ater a single time step o = 10 6 when the exact curvature is used. The errors are machine precision zero even or large density ratios. In the practical simulations, however, the non-zero spurious current is produced because o errors in the curvature calculation Long-time evolution o the viscous sphere In the case o section 4.1, the exact curvature produces machine-zero spurious currents. However, the numerical errors in curvature evaluation result in spurious currents which are small, but non-zero. The long-time behavior o the equilibrium sphere, thus, has been simulated to check errors in velocities ater long time steps. A sphere o diameter D = 0.25 is placed at the center o a unit sized Cartesian box with Neumann boundary
7 Pressure-correction algorithm or ast two-phase simulations 373 ρ 1/ρ 2 L (u) L (p max p min p ex) Table 1. Errors in velocity and pressure ater single time step = 10 6 or dierent density ratios using the exact curvature. h, h G 1/128 1/160 1/192 Ca Ca (Renardy & Renardy) Table 2. Spurious current capillary numbers o the viscous sphere with those o Renardy & Renardy (2002). conditions at the sides. The surace tension coeicient is set to σ = 0.357, density and viscosity in both luids are ρ = 4 and µ = 1. The time step size is = Table 2 shows the maximum spurious current capillary number at t = compared with Renardy & Renardy (2002). Both methods exhibit good results Oscillation o column The oscillation period or the column without gravity has been perormed to veriy our method. The theoretical period in the linear regime is given by (Lamb 1945) as w 2 = n(n2 1)σ (ρ 1 + ρ 2 )R0 3. (4.1) A column o radius R 0 = 2 is placed at the center o domain and the parameters are ρ 1 = 1, ρ 2 = 1000, µ 1 = 0.01, µ 2 = , and σ = 1. The column is initially perturbed by a mode n = 2 perturbation with initial amplitude o A 0 = 0.01R 0. The same equidistance Cartesian grid is used or CDP-grid and G-grid. Figure 2 shows the error o the oscillation period with respect to the number o grids. The second order grid convergence is achieved as shown Broken dam problem The well-known broken problem is computed to compare with the experimental result o Martin & Moyce (1952). A water column with a = 0.06m wide and 2a = 0.12m high is initially located at rest. The column suddenly breaks down under the gravity g = 9.81m/s 2. No slip conditions are applied or all side boundaries. The Cartesian grid o = m is used here. Figure 3 shows the time history o the phase interace. In
8 374 D. Kim and P. Moin Figure 2. Error in oscillation period with respect to the grid size. Figure 3. Time history o the phase interace in the broken dam problem at the dierent nondimensional times t 2g/a=0, 2.0, 3.0, 5.0. order to compare with the experimental data (Martin & Moyce 1952), the ront position and the height o the column are presented in igure 4. It shows good agreement with the experimental data Computational perormance Since the algebraic multigrid method (AMG) can be used or the Poisson solver having constant coeicients, the present method shows much aster perormance than the previous method which uses the biconjugate gradient stabilized method (BCGSTAB) or the Poisson solver. Table 3 shows the computational costs measured in CPU seconds per time step on a tetrahedral mesh having 1 million control volumes. In the previous method using BCGSTAB, the level set advection takes only 0.06% o the total time. About 98.8%
9 Pressure-correction algorithm or ast two-phase simulations 375 Figure 4. (a) time history o the height o the column (b) time history o the ront position o the column. : present result; : experimental data (Martin & Moyce 1952) Previous method (BCGSTAB) Present method (AMG) G advection Reinitialization Momentum solver Poisson solver Total CPU time Table 3. Comparison o computational cost measured in CPU seconds per time step on a tetrahedral mesh o 1 millions. o the total time step is spent on solving the Poisson equation. However, in case o the present method using AMG, it is more than 10 times aster than the previous method even with several iterations. 5. Conclusions and uture work A new pressure correction algorithm is presented or ast and accurate two-phase simulations. A ractional step method is used to solve the Navier-Stokes equations, and the pressure is used to correct the velocity ield to satisy the mass conservation equation instead o divergence-ree velocity. This procedure yields the Poisson system having constant coeicients, which saves computational cost greatly. In the present method, the mass conservation equation is ulilled at every time step, but the incompressibility is not ully guaranteed due to the numerical errors. These errors can be eiciently reduced by reining the grid. The canonical test problems are solved to demonstrate the accuracy and much improved perormance o the method. The main disadvantage o the present method is that the incompressibility is not ulilled owing to the numerical errors. In the uture, we will ocus on reducing or eliminating the errors in incompressibility. Finally, the present method will be used to study air layer drag phenomenon in turbulent low using more than 500 millions grid points.
10 376 D. Kim and P. Moin Acknowledgments The work presented in this paper was supported by the Oice o Naval Research. REFERENCES Alonso, J. J., Hahn, S., Ham, F., Herrmann, M., Iaccarino, G., Kalitzin, G., LeGresley, P., Mattsson, K., Medic, G., Moin, P., Pitsch, H., Schlüter, J., Svärd, M., van der Weide, E., You, D. & Wu, X Chimps: A highperormance scalable module or multi-physics simulations. AIAA Paper Apte, S. V., Gorokhovski, M. & Moin, P Les o atomizing spray with stochastic modeling o secondary breakup. Int. J. Mult. Flow 29, Chang, Y. C., Hou, T. Y., Merriman, B. & Osher, S A level set ormulation o eulerian interace capturing methods or incompressible luid lows. Journal o Computational Physics 124, Francois, M. M., Cummins, S. J., Dendy, E. D., Kothe, D. B., Sicilian, J. M. & Williams, M. W A balanced-orce algorithm or continuous and sharp interacial surace tension models within a volume tracking ramework. Journal o Computational Physics 213, Ham, F. & Iaccarino, G Energy conservation in collocated discretization schemes on unstructured meshes. Annual Research Brie 2004, Herrmann, M A balanced orce reined level set grid method or two-phase lows on unstructured low solver grids. Journal o Computational Physics 227, Jiang, G.-S. & Peng, D Weighted eno schems or hamilton-jacobi equations. SIAM J. Sci. Comp. 21, Kim, D., Herrmann, M. & Moin, P Numerical simulation o high shear rate two-phase low with high density ratio. Bulletin o the American Physical Society. Lamb, H Hydrodynamics. Dover, New York. Mahesh, K., Constantinescu, G. & Moin, P A numerical method or large eddy simulation in complex geometries. Journal o Computational Physics 197, Martin, J. & Moyce, W An experimental study o the collapse o liquid columns on a rigid horizontal plane. Philos. Trans. R. Soc. A224, Renardy, Y. & Renardy, M Prost: a parabolic reconstruction o surace tension or the volume-o-luid method. Journal o Computational Physics 183, Shu, C.-W. & Osher, S Eicient implementation o essentially non-oscillatory shock-capturing schemes. Journal o Computational Physics 77, Sussman, M. & Puckett, E A coupled level set and volume-o-luid method or computing 3d and axisymmetric incompressible two-phase lows. Journal o Computational Physics 162, Van Der Pill, S. P., Segal, A. & Vuik, C A mass-conserving level-set method or modelling o multi-phase lows. Int. J. Numer. Meth. Fluids 47,
Toward two-phase simulation of the primary breakup of a round liquid jet by a coaxial flow of gas
Center for Turbulence Research Annual Research Briefs 2006 185 Toward two-phase simulation of the primary breakup of a round liquid jet by a coaxial flow of gas By D. Kim, O. Desjardins, M. Herrmann AND
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationA Cartesian adaptive level set method for two-phase flows
Center or Turbulence Research Annual Research Bries 2003 227 A Cartesian adaptive level set method or two-phase lows By F. Ham AND Y.-N. Young 1. Motivation and objectives Simulations o lows involving
More informationAn Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains
An Immersed Interace Method or the Incompressible Navier-Stokes Equations in Irregular Domains Duc-Vinh Le, Boo Cheong Khoo, Jaime Peraire Singapore-MIT Alliance Department o Mechanical Engineering, National
More informationReduction of parasitic currents in the DNS VOF code FS3D
M. Boger a J. Schlottke b C.-D. Munz a B. Weigand b Reduction of parasitic currents in the DNS VOF code FS3D Stuttgart, March 2010 a Institut für Aerodynamik und Gasdynamik, Universität Stuttgart, Pfaffenwaldring
More informationDEVELOPMENT AND VALIDATION OF HIGH-PRECISION CFD METHOD WITH VOLUME-TRACKING ALGORITHM FOR GAS-LIQUID TWO-PHASE FLOW SIMULATION ON UNSTRUCTURED MESH
DEVELOPMENT AND VALIDATION OF HIGH-PRECISION CFD METHOD WITH VOLUME-TRACKING ALGORITHM FOR GAS-LIQUID TWO-PHASE FLOW SIMULATION ON UNSTRUCTURED MESH I1. Kei Ito 1, 2, I2. Tomoaki Kunugi 1, I3. Hiroyuki
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationA Discontinuous Galerkin Conservative Level Set Scheme for Simulating Turbulent Primary Atomization
ILASS-Americas 3rd Annual Conference on Liquid Atomization and Spray Systems, Ventura, CA, May 011 A Discontinuous Galerkin Conservative Level Set Scheme for Simulating Turbulent Primary Atomization M.
More informationA semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows
A semi-implicit finite volume implementation of the CSF method for treating surface tension in interfacial flows M. Raessi, M. Bussmann*, and J. Mostaghimi Department of Mechanical and Industrial Engineering,
More informationDevelopment and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver
More informationA Level Set/Vortex Sheet Method for Modeling Phase Interface Dynamics During Primary Breakup
ILASS Americas 17th Annual Conference on Liquid Atomization and Spray Systems, Arlington, VA, May 2004 A Level Set/Vortex Sheet Method for Modeling Phase Interface Dynamics During Primary Breakup M. Herrmann
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationDetailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and Boundary Condition
ILASS-Americas 25th Annual Conference on Liquid Atomization and Spray Systems, Pittsburgh, PA, May 23 Detailed Numerical Simulation of Liquid Jet in Cross Flow Atomization: Impact of Nozzle Geometry and
More informationFEM-Level Set Techniques for Multiphase Flow --- Some recent results
FEM-Level Set Techniques for Multiphase Flow --- Some recent results ENUMATH09, Uppsala Stefan Turek, Otto Mierka, Dmitri Kuzmin, Shuren Hysing Institut für Angewandte Mathematik, TU Dortmund http://www.mathematik.tu-dortmund.de/ls3
More informationControlling the Heat Flux Distribution by Changing the Thickness of Heated Wall
J. Basic. Appl. Sci. Res., 2(7)7270-7275, 2012 2012, TextRoad Publication ISSN 2090-4304 Journal o Basic and Applied Scientiic Research www.textroad.com Controlling the Heat Flux Distribution by Changing
More informationA PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS
A PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS ROBERTO CROCE, MICHAEL GRIEBEL, AND MARC ALEXANDER SCHWEITZER Abstract. In this paper we present
More informationDevelopment and Testing of a Robust Free-Surface Finite Volume Method. Marwan Darwish, Ph.D.
Development and Testing o a Robust Free-Surace Finite Volume Method Marwan Darwish, Ph.D. Faculty o Engineering and Architecture, American University o Beirut July 2003 Summary A robust numerical technique
More informationBOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE WITH MAGNETIC FIELD IN A NANOFLUID
Proceedings o the International Conerence on Mechanical Engineering and Reneable Energy 7 (ICMERE7) 8 December, 7, Chittagong, Bangladesh ICMERE7-PI- BOUNDARY LAYER ANALYSIS ALONG A STRETCHING WEDGE SURFACE
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationAnalysis of Non-Thermal Equilibrium in Porous Media
Analysis o Non-Thermal Equilibrium in Porous Media A. Nouri-Borujerdi, M. Nazari 1 School o Mechanical Engineering, Shari University o Technology P.O Box 11365-9567, Tehran, Iran E-mail: anouri@shari.edu
More informationA dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries
Center for Turbulence Research Annual Research Briefs 2006 41 A dynamic global-coefficient subgrid-scale eddy-viscosity model for large-eddy simulation in complex geometries By D. You AND P. Moin 1. Motivation
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationA Critical Investigation of High-Order Flux Limiters In Multiphase Flow Problems
A Critical Investigation o High-Order Flux Limiters In Multiphase Flow Problems Chris Guenther Fluent In., 3647 Collins Ferry Rd., Morgantown, WV 26505, USA cpg@luent.com ABSTRACT. In recent years inite
More informationPressure corrected SPH for fluid animation
Pressure corrected SPH for fluid animation Kai Bao, Hui Zhang, Lili Zheng and Enhua Wu Analyzed by Po-Ram Kim 2 March 2010 Abstract We present pressure scheme for the SPH for fluid animation In conventional
More information.u= 0 ρ( u t +(u. )u)= ρ g p+.[µ( u+ t u)]
THETIS is a numerical simulation tool developed by University of Bordeaux. It is a versatile code to solve different problems: fluid flows, heat transfers, scalar transports or porous mediums. The potential
More informationFlow Field and Oscillation Frequency of a Rotating Liquid Droplet
Flow Field and Oscillation Frequency of a Rotating Liquid Droplet TADASHI WATANABE Center for Computational Science and e-systems Japan Atomic Energy Agency (JAEA) Tokai-mura, Naka-gun, Ibaraki-ken, 319-1195
More informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
More informationIterative Methods for Stokes/Darcy Coupling
Iterative Methods or Stokes/Darcy Coupling Marco Discacciati Ecole Polytechnique Fédérale de Lausanne Institut d Analyse et Calcul Scientiique Chair o Modelling and Scientiic Computing marco.discacciati@epl.ch
More informationAn Overview of Fluid Animation. Christopher Batty March 11, 2014
An Overview of Fluid Animation Christopher Batty March 11, 2014 What distinguishes fluids? What distinguishes fluids? No preferred shape. Always flows when force is applied. Deforms to fit its container.
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
More informationNumerical Studies of Droplet Deformation and Break-up
ILASS Americas 14th Annual Conference on Liquid Atomization and Spray Systems, Dearborn, MI, May 2001 Numerical Studies of Droplet Deformation and Break-up B. T. Helenbrook Department of Mechanical and
More informationBuoyancy Driven Heat Transfer of Water-Based CuO Nanofluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center
July 4-6 2012 London U.K. Buoyancy Driven Heat Transer o Water-Based CuO Nanoluids in a Tilted Enclosure with a Heat Conducting Solid Cylinder on its Center Ahmet Cihan Kamil Kahveci and Çiğdem Susantez
More informationSurface phase separation and flow in a simple model of drops and vesicles
Surface phase separation and flow in a simple model of drops and vesicles Tutorial Lecture 4 John Lowengrub Department of Mathematics University of California at Irvine Joint with J.-J. Xu (UCI), S. Li
More informationOn the Development of Implicit Solvers for Time-Dependent Systems
School o something FACULTY School OF OTHER o Computing On the Development o Implicit Solvers or Time-Dependent Systems Peter Jimack School o Computing, University o Leeds In collaboration with: P.H. Gaskell,
More informationA surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant
A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant Ashley J. James Department of Aerospace Engineering and Mechanics, University of Minnesota John Lowengrub
More informationOPTIMALLY STAGGERED FINNED CIRCULAR AND ELLIPTIC TUBES IN FORCED CONVECTION
OPTIMALLY STAGGERED FINNED CIRCULAR AND ELLIPTIC TUBES IN FORCED CONVECTION R. S. Matos a, T. A. Laursen b, J. V. C. Vargas a, and A. Bejan c, a Universidade Federal do Paraná Departamento de Engenharia
More informationAccurate and stable finite volume operators for unstructured flow solvers
Center for Turbulence Research Annual Research Briefs 2006 243 Accurate and stable finite volume operators for unstructured flow solvers By F. Ham, K. Mattsson AND G. Iaccarino 1. Motivation and objectives
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
More information2015 American Journal of Engineering Research (AJER)
American Journal o Engineering Research (AJER) 2015 American Journal o Engineering Research (AJER) e-issn: 2320-0847 p-issn : 2320-0936 Volume-4, Issue-7, pp-33-40.ajer.org Research Paper Open Access The
More informationNonlinear oscillations and rotations of a liquid droplet
Nonlinear oscillations and rotations of a liquid droplet Tadashi Watanabe watanabe.tadashi66@jaea.go.jp Abstract Nonlinear oscillations and rotations of a liquid droplet are simulated numerically by solving
More informationLecture 8 Optimization
4/9/015 Lecture 8 Optimization EE 4386/5301 Computational Methods in EE Spring 015 Optimization 1 Outline Introduction 1D Optimization Parabolic interpolation Golden section search Newton s method Multidimensional
More informationLarge Scale Simulations of Turbulent Flows for Industrial Applications. Lakhdar Remaki BCAM- Basque Centre for Applied Mathematics
Large Scale Simulations o Turbulent Flows or Industrial Alications Lakhdar Remaki BCAM- Basque Centre or Alied Mathematics Outline Flow Motion Simulation Physical Model Finite-Volume Numerical Method CFD
More informationCONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS. Convective heat transfer analysis of nanofluid flowing inside a
Chapter 4 CONVECTIVE HEAT TRANSFER CHARACTERISTICS OF NANOFLUIDS Convective heat transer analysis o nanoluid lowing inside a straight tube o circular cross-section under laminar and turbulent conditions
More informationA Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes
A Scalable, Parallel Implementation of Weighted, Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory SIAM Annual Meeting
More informationEquipment Safety in Renewable Energies Exploitation
European Association or the Development o Renewable Energies, Environment and Power Quality (EAEPQ International Conerence on Renewable Energies and Power Quality (ICREPQ 11 Las Palmas de Gran Canaria
More informationOblique Drop Impact on Deep and Shallow Liquid
1 2 3 4 5 6 7 8 9 11 10 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Commun. Comput. Phys. doi: 10.4208/cicp.XXX.XXX Oblique Drop Impact on Deep and Shallow Liquid B. Ray 1, G. Biswas 1,2, and A. Sharma 3
More informationA numerical analysis of drop impact on liquid film by using a level set method
Journal of Mechanical Science and Technology 25 (10) (2011) 2567~2572 wwwspringerlinkcom/content/1738-494x DOI 101007/s12206-011-0613-7 A numerical analysis of drop impact on liquid film by using a level
More informationarxiv: v1 [math.na] 29 Oct 2018
A Multirate Approach or Fluid-Structure Interaction Computation with Decoupled Methods Lian Zhang 1, Mingchao Cai 2 Mo Mu 3, arxiv:1810.11933v1 [math.na] 29 Oct 2018 Abstract We investigate a multirate
More informationThe Relationship between Discrete Calculus Methods and other Mimetic Approaches
The Relationship between Discrete Calculus Methods and other Mimetic Approaches Blair Perot Theoretical and Computational Fluid Dynamics Laboratory Woudschoten Conerence Oct. 2011 October 7, 2011 Background
More informationEfficient simulation techniques for incompressible two-phase flow
3D-Surface Engineering für Werkzeugsysteme der Blechformteilefertigung - Erzeugung, Modellierung, Bearbeitung - Efficient simulation techniques for incompressible two-phase flow O. Mierka, O. Ouazzi, T.
More informationSecond Order Slip Flow of Cu-Water Nanofluid Over a Stretching Sheet With Heat Transfer
Second Order Slip Flow o Cu-Water Nanoluid Over a Stretching Sheet With Heat Transer RAJESH SHARMA AND ANUAR ISHAK School o Mathematical Sciences, Faculty o Science and Technology Universiti Kebangsaan
More informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
More information2 Introduction Standard conservative schemes will admit non-physical oscillations near material interfaces, which can be avoided by the application of
A Conservative Approach to the Multiphase Euler Equations without Spurious Pressure Oscillations Xu-Dong Liu Ronald P. Fedkiw Stanley Osher June 8, 998 Abstract In [4], the pressure evolution equation
More informationLarge-eddy simulation analysis of turbulent combustion in a gas turbine engine combustor
Center for Turbulence Research Annual Research Briefs 2008 219 Large-eddy simulation analysis of turbulent combustion in a gas turbine engine combustor By D. You, F. Ham AND P. Moin 1. Motivation and objectives
More informationInvestigation of Cavitation Bubble Cloud with discrete Lagrangian Tracking
Investigation of Cavitation Bubble Cloud with discrete Lagrangian Tracking Xiuxiu Lyu, Xiangyu Hu, Nikolaus A. Adams Chair of Aerodynamics and Fluid Mechanics, Department of Mechanical Engineering, Technical
More informationWave regimes and mass transfer in two-layer falling films
Wave regimes and mass transer in two-layer alling ilms by GÖKÇEN ÇEKİÇ A thesis submitted to The University o Birmingham or the degree o Doctor o Philosophy School o Mathematics The University o Birmingham
More informationHeight functions for applying contact angles to 2D VOF simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids (2007) Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 0.002/fld.65 Height functions for applying
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More information1 General description of the OpenFOAM suite
Santiago Márquez Damián-Final Work-Computational Fluid Dynamics Description and utilization o interfoam multiphase solver 1 General description o the OpenFOAM suite OpenFOAM Open Field Operation and Manipulation
More informationIMPINGEMENT OF A DROPLET ONTO A DRY WALL: A NUMERICAL INVESTIGATION
IMPINGEMENT OF A DROPLET ONTO A DRY WALL:. Introduction A NUMERICAL INVESTIGATION N. Nikolopoulos, and G. Bergeles Department Mechanical Engineering Nat. Technical University of Athens 57 Zografos, Greece
More informationChapter 3 Water Flow in Pipes
The Islamic University o Gaza Faculty o Engineering Civil Engineering Department Hydraulics - ECI 33 Chapter 3 Water Flow in Pipes 3. Description o A Pipe Flow Water pipes in our homes and the distribution
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
More informationNumerical study of turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2014 41 Numerical study of turbulent two-phase Couette flow By A.Y. Ovsyannikov, D. Kim, A. Mani AND P. Moin 1. Motivation and objectives The motivation
More informationReview of Numerical Methods for Multiphase Flow
Title Review o Numerical Methods or Multiphase Flow M. Sommereld Mechanische Verahrenstechnik Zentrum ür Ingenieurwissenschaten 06099 Halle (Saale), Germany www-mvt.iw.uni-halle.de Content o the Lecture
More informationInvestigation of an implicit solver for the simulation of bubble oscillations using Basilisk
Investigation of an implicit solver for the simulation of bubble oscillations using Basilisk D. Fuster, and S. Popinet Sorbonne Universités, UPMC Univ Paris 6, CNRS, UMR 79 Institut Jean Le Rond d Alembert,
More informationMAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS
THERMAL SCIENCE: Year 8, Vol., No. B, pp. 383-39 383 MAGNETOHYDRODYNAMIC GO-WATER NANOFLUID FLOW AND HEAT TRANSFER BETWEEN TWO PARALLEL MOVING DISKS Introduction by Mohammadreza AZIMI and Rouzbeh RIAZI
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationNon-newtonian Rabinowitsch Fluid Effects on the Lubrication Performances of Sine Film Thrust Bearings
International Journal o Mechanical Engineering and Applications 7; 5(): 6-67 http://www.sciencepublishinggroup.com/j/ijmea doi:.648/j.ijmea.75.4 ISSN: -X (Print); ISSN: -48 (Online) Non-newtonian Rabinowitsch
More informationA Volume of Fluid Dual Scale Approach for Modeling Turbulent Liquid/Gas Phase Interfaces
ILASS-Americas 29th Annual Conference on Liquid Atomization and Spray Systems, Atlanta, GA, May 217 A Volume of Fluid Dual Scale Approach for Modeling Turbulent Liquid/Gas Phase Interfaces D. Kedelty,
More informationDirect numerical simulation of compressible multi-phase ow with a pressure-based method
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2902 Direct numerical simulation of compressible multi-phase ow with a pressure-based
More informationNON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION
NON-SIMILAR SOLUTIONS FOR NATURAL CONVECTION FROM A MOVING VERTICAL PLATE WITH A CONVECTIVE THERMAL BOUNDARY CONDITION by Asterios Pantokratoras School o Engineering, Democritus University o Thrace, 67100
More information3D Numerical Modelling of Convective Heat Transfer through Two-sided Vertical Channel Symmetrically Filled with Metal Foams
P Periodica Polytechnica Mechanical Engineering P 60(4), pp. 193-202, 2016 DOI: 10.3311/PPme.8511 Creative Commons Attribution b 3D Numerical Modelling o Convective Heat Transer through Two-sided Vertical
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationDescription of a One-Dimensional Numerical Model of an Active Magnetic Regenerator Refrigerator
This is a 1D model o an active magnetic regenerative rerigerator (AMRR) that was developed in MATLAB. The model uses cycle inputs such as the luid mass low and magnetic ield proiles, luid and regenerator
More informationGRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE
27th International Conference on Parallel Computational Fluid Dynamics Parallel CFD2015 GRAVITY-DRIVEN MOTION OF A SWARM OF BUBBLES IN A VERTICAL PIPE Néstor Balcázar, Oriol Lehmkuhl, Jesús Castro, Joaquim
More informationSurface tension in incompressible Rayleigh Taylor mixing flow
Journal of Turbulence Volume 7, No. 71, 6 Surface tension in incompressible Rayleigh Taylor mixing flow Y.-N. YOUNG and F. E. HAM Department of Mathematical Sciences, New Jersey Institute of Technology,
More informationNonlinear shape evolution of immiscible two-phase interface
Nonlinear shape evolution of immiscible two-phase interface Francesco Capuano 1,2,*, Gennaro Coppola 1, Luigi de Luca 1 1 Dipartimento di Ingegneria Industriale (DII), Università di Napoli Federico II,
More informationarxiv: v1 [physics.flu-dyn] 5 Aug 2017
On a relation between the volume of fluid, level-set and phase field interface models arxiv:1708.01805v1 [physics.flu-dyn] 5 Aug 2017 Abstract Tomasz Wac lawczyk, Department of Aeronautics, Institute of
More informationMechanical Engineering Research Journal BUOYANT FLOW OF NANOFLUID FOR HEAT-MASS TRANSFER THROUGH A THIN LAYER
Dept. o Mech. Eng. CUET Published Online March 2015 (http://www.cuet.ac.bd/merj/index.html) Mechanical Engineering Research Journal Vol. 9, pp. 712, 2013 M E R J ISSN: 1990-5491 BUOYANT FLOW OF NANOFLUID
More informationNumerical Simulation of Particle Mixing in Dispersed Gas-Liquid-Solid Flows using a Combined Volume of Fluid and Discrete Particle Approach
6 th International Conerence on Multiphase Flow, Numerical Simulation o Particle Mixing in Dispersed Gas-Liquid-Solid Flows using a Combined Volume o Fluid and Discrete Particle Approach Niels G. Deen,
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationGlobal Weak Solution of Planetary Geostrophic Equations with Inviscid Geostrophic Balance
Global Weak Solution o Planetary Geostrophic Equations with Inviscid Geostrophic Balance Jian-Guo Liu 1, Roger Samelson 2, Cheng Wang 3 Communicated by R. Temam) Abstract. A reormulation o the planetary
More informationx j r i V i,j+1/2 r Ci,j Ui+1/2,j U i-1/2,j Vi,j-1/2
Merging of drops to form bamboo waves Yuriko Y. Renardy and Jie Li Department of Mathematics and ICAM Virginia Polytechnic Institute and State University Blacksburg, VA -, U.S.A. May, Abstract Topological
More informationComputers & Fluids. Improving lattice Boltzmann simulation of moving particles in a viscous flow using local grid refinement
Computers & Fluids 00 (204) 2 Computers & Fluids Improving lattice Boltzmann simulation o moving particles in a viscous low using local grid reinement Songying Chen a,, Cheng Peng b, Yihua Teng c, Lian-Ping
More informationA Momentum Exchange-based Immersed Boundary-Lattice. Boltzmann Method for Fluid Structure Interaction
APCOM & ISCM -4 th December, 03, Singapore A Momentum Exchange-based Immersed Boundary-Lattice Boltzmann Method for Fluid Structure Interaction Jianfei Yang,,3, Zhengdao Wang,,3, and *Yuehong Qian,,3,4
More informationRotating Flow of Magnetite-Water Nanofluid over a Stretching Surface Inspired By Non-Linear Thermal Radiation and Mass Transfer
International Journal o Mathematics Research. ISSN 0976-5840 Volume 9, Number (017), pp. 89-97 International Research Publication House http://www.irphouse.com Rotating Flow o Magnetite-Water Nanoluid
More informationTHE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION
THE MODEL OF DRYING SESSILE DROP OF COLLOIDAL SOLUTION I.V.Vodolazskaya, Yu.Yu.Tarasevic Astrakan State University, a Tatiscev St., Astrakan, 41456, Russia e-mail: tarasevic@aspu.ru e ave proposed and
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationSolving Partial Differential Equations Numerically. Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah
Solving Partial Dierential Equations Numerically Miklós Bergou with: Gary Miller, David Cardoze, Todd Phillips, Mark Olah Overview What are partial dierential equations? How do we solve them? (Example)
More informationExperiments and Numerical Simulations of the Flow Within a Simplified Model of a Surge Chamber
Experiments and Numerical Simulations o the Flow Within a Simpliied Model o a Surge Chamber by Alexandre Massé* 1,3, Maryse age2 and Laurent Mydlarski3 1 Groupe-Conseil LaSalle, LaSalle (Qc), H8R 1R8,
More informationSIMULATION OF A FLOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT
SIMULATION O A LOW THROUGH A RIGID POROUS MEDIUM SUBJECTED TO A CONSTRAINT R. M. Saldanha da Gama a, J. J. Pedrosa ilho a, and M. L. Martins-Costa b a Universidade do Estado do Rio de Janeiro Departamento
More informationFlow and Heat Transfer Analysis of Copper-water Nanofluid with Temperature Dependent Viscosity Past a Riga Plate
Journal o Magnetics (), 181-187 (017) ISSN (Print) 16-1750 ISSN (Online) 33-6656 https://doi.org/10.483/jmag.017...181 Flo and Heat Transer Analysis o Copper-ater Nanoluid ith Temperature Dependent Viscosity
More informationSimulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang,
Simulation of T-junction using LBM and VOF ENERGY 224 Final Project Yifan Wang, yfwang09@stanford.edu 1. Problem setting In this project, we present a benchmark simulation for segmented flows, which contain
More informationTransport Properties: Momentum Transport, Viscosity
Transport Properties: Momentum Transport, Viscosity 13th February 2011 1 Introduction Much as mass(material) is transported within luids (gases and liquids), linear momentum is also associated with transport,
More informationHYDROMAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID
Rita Choudhury et al. / International Journal o Engineering Science and Technology (IJEST) HYDROAGNETIC DIVERGENT CHANNEL FLOW OF A VISCO- ELASTIC ELECTRICALLY CONDUCTING FLUID RITA CHOUDHURY Department
More information