A Volume of Fluid Dual Scale Approach for Modeling Turbulent Liquid/Gas Phase Interfaces
|
|
- Rodney Powell
- 6 years ago
- Views:
Transcription
1 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, J. Uglietta, and M. Herrmann School for Engineering of Matter, Transport and Energy Arizona State University Tempe, AZ USA Abstract In atomizing flows, the properties of the resulting liquid spray are determined by the interplay of fluid and surface tension forces. The resulting flow and interface dynamics typically span 4-6 orders of magnitude in length scales, making direct numerical simulations exceedingly expensive. This motivates the need for modeling approaches based on spatial filtering or ensemble averaging. In this contribution, advances to a dual-scale modeling approach (Herrmann, 213) are presented to describe turbulent two-phase interface dynamics in a large-eddy-simulation-type spatial filtering context. Spatial filtering of the governing equations introduces several sub-filter terms that require modeling (Toutant et al., 28). These include terms associated with sub-filter acceleration, transport, viscous effects, and surface tension. Instead of developing individual closure-models for each of the terms associated with the phase interface, the proposed dual-scale approach uses an exact closure by explicitly filtering a fully resolved realization of the phase interface. This resolved realization is maintained on a high-resolution over-set mesh using a Refined Local Surface Grid approach (Herrmann, 28) employing an un-split, geometric, bounded, and conservative Volume-of-Fluid method (Owkes and Desjardins, 214). The advection equation for the phase interface on this DNS scale requires a model for the fully resolved interface advection velocity. This velocity is the sum of the LES filtered velocity, readily available from the LES approach solving the filtered Navier Stokes equation, and the sub-filter velocity fluctuation that has two contributions. The first is due to sub-filter turbulent eddies, the second is due to sub-filter surface tension forces. We propose four different methods to reconstruct the sub-filter turbulent eddy fluctuation velocity using the fractal interpolation approach of Scotti & Meneveau (1999). Results of the dual-scale model are compared to recent direct numerical simulations of a phase interface in homogeneous isotropic turbulence (McCaslin and Desjardins, 215). Corresponding Author: marcus.herrmann@asu.edu
2 Introduction Atomization in turbulent environments involves a vast range of length and time scales. Predictive simulations aiming to resolve all relevant scales thus require enormous computational resources, taxing even the most powerful computers available today [1]. Since primary atomization is governed by the dynamics of the interface, a need therefore exists for appropriate interface models that make the computational cost of predicting the atomization outcome more tractable. A wide range of phenomenological models aiming to represent statistically the essential features of atomization have been proposed in the past. Although these aim to introduce the dominant mechanisms for breakup, they often use round blobs injected from the nozzle exit and hence neglect all details of the interface dynamics. Other modeling approaches to atomization include stochastic models [2, 3] representing the interface by constituent stochastic particles and the mean interface density transport equation model for Reynolds-Averaged Navier-Stokes (RANS) approaches [4, 5]. The former treats the interface dynamics in a stochastic sense but requires the a priori knowledge of the breakup mechanism, whereas the latter is affected by the drawbacks of the RANS approach: the transport of the mean interface density is modeled by a diffusion-like hypothesis, thereby neglecting the spatial grouping effects of liquid elements [1]. In the context of Large Eddy Simulations (LES), [6, 7, 8, 9] have proposed models to close the unclosed terms arising from the introduction of spatial filtering into the governing equations. However, these models typically neglect the contribution of the sub-filter surface tension term and are based on a cascade process hypothesis that may be questionable in the context of surface tension-driven atomization. In [1, 11], a dual-scale approach for LES of interface dynamics was proposed and a model for the sub-filter surface tension induced motion of phase interfaces was developed. The purpose of this contribution is to develop a model for the sub-filter phase interface motion induced by sub-filter turbulent velocity fluctuations. Combining such a model with the surface tension model proposed in [1, 11] will result in a LES model applicable to atomizing flows. Governing equations The equations governing the fully resolved motion of an unsteady, incompressible, immiscible, twofluid system in the absence of surface tension are the Navier-Stokes equations, ρu t + (ρu u) = p + (µ ( u + T u )), (1) where u is the velocity, ρ the density, p the pressure, and µ the dynamic viscosity. Here, we neglect surface tension to solely focus on the turbulence induced dynamics of phase interfaces. Furthermore, the continuity equation results in a divergence-free constraint on the velocity field u =. (2) Assuming ρ and µ are constant within each fluid, density and viscosity can be calculated from ρ = ψρ l + (1 ψ)ρ g (3) µ = ψµ l + (1 ψ)µ g, (4) where indices l and g denote values in liquid and gas, respectively, and ψ is a volume-of-fluid scalar that is ψ = in the gas and ψ = 1 in the liquid with ψ t = u ψ = (uψ) + ψ u. (5) Here, the last term on the right-hand side is zero for incompressible flows, see Eq. (2). Filtered governing equations Introducing spatial filtering into Eqs. (1) and (2) and assuming that the filter commutes with both the time and spatial derivatives, the filtered governing equations can be derived [7], ρ u t + (ρū u) = p + (6) (µ( u + T u)) + τ 1 + (τ 2 + τ 3 ), ū =, (7) where indicates spatial filtering, and τ 1 = ρ u ρu (8) t t τ 2 = ρ u u ρu u (9) τ 3 = µ( u + T u) µ( u + T u), (1) where τ 1, τ 2, and τ 3 are associated, respectively, with acceleration, advection, and viscous effects [7]. Using Eqs. (3) and (4), the filtered density and viscosity in Eq. (6) are ρ = ρ l ψ + ρ g (1 ψ) (11) µ = µ l ψ + µ g (1 ψ), (12) 2
3 where ψ = G(x)ψdx, (13) and G is a normalized spatial filter function. The dual-scale approach to modeling subfilter interface dynamics Instead of relying on a cascade process by which dynamics on a sub-filter scale can be inferred from the dynamics on the resolved scale, the dual-scale approach proposed in [11] aims to maintain a fully resolved realization of the immiscible interface geometry at all times, expressed, for example, in terms of a volume-of-fluids scalar ψ. Then ψ can be calculated exactly by explicit filtering using Eq. (13). Although this is an exact closure, the problem of modeling is of course simply shifted to the problem of maintaining a fully resolved realization of the interface geometry, i.e., describing the fully resolved motion of the interface, Eq. (5). Since the fully resolved velocity is the sum of the filtered velocity and the sub-grid velocity, u = u + u sg, this results in ψ t = ((u + u sg) ψ) + ψ (u + u sg ), (14) where the only term requiring modeling is u sg. In [11], a model for u sg is proposed consisting of three contributions, u sg = u + δu + u σ, (15) where u is due to sub-filter turbulent eddies, δu is attributed to the interface velocity increment due to relative sub-filter motion between the two immiscible fluids, and u σ is due to sub-filter velocities induced by sub-filter surface tension forces. The focus of the current contribution is on the first term; for a modeling outline of the second term, the reader is referred to [11, 14], and for modeling of the last term, the reader is referred to [1]. Sub-filter turbulent fluctuation velocity models We propose to reconstruct the sub-filter turbulent fluctuation velocity u using fractal interpolation [15]. To demonstrate fractal interpolation in one dimension, consider 3 adjacent LES scale nodes x, x 1, and x 2 with velocities u, u 1, and u 2. Following [15, 16] the fractal interpolation operator W F I can be written as W F I (x) = u + u 1 u x 1 x (x x ) + d 1 (u(2x x ) u (16) ) u 2 u (2x x ) if x [x, x 1 ] x 2 x W F I (x) = u 1 + u 2 u 1 x 2 x 1 (x x 1 ) + d 2 (u(2x x ) u (17) ) u 2 u (2x x ) if x [x 1, x 2 ] x 2 x Here d 1 < 1 and d 2 < 1 are stretching factors making W F I a contractive mapping [15]. Successively applying the fractal interpolation operator W F I starting with the LES filter velocities, generates the fully resolved turbulent fluctuation velocity. Two different approaches are pursued to extend the above formulation to three dimensions. In the first, fractal interpolation is first performed in one spatial direction only, followed by separate 1D fractal interpolations in the other two directions [15]. The second approach replaces the one-dimensional linear interpolation operators in Eqs. (16) and (17) with three-dimensional trilinear interpolation operators and performs the fractal interpolation directly in three dimensions [17]. Furthermore, two different approaches are used to determine the values of the stretching factors d 1 and d 2. The first, follows the so-called ZE1 model of [15] by using d 1 = d 2 = ±2 1/3 with the sign chosen randomly with equal probability. This choice of d i generates a velocity signal that satisfies the 5/3 kinetic energy spectrum of turbulence. The second approach chooses d i with random sign and absolute value determined from a log- Poisson distribution [16, 17], P with ( d i = ( ) γ ) 1 β n = e 2 λ λn, n =, 1,... (18) n! λ = 1 3γ ln 2 (19) 1 β3 and γ = 1/9 and β = (2/3) 1/3 [17]. Finally, two different approaches are analyzed concerning the spatial location of the LES velocities u, u 1, and u 2. If the LES flow solver utilizes a staggered mesh layout, face normal velocities are not colocated and hence the fractal interpolation has to be performed for different locations depending on the spatial component of the velocity vector. However, if one first interpolates the staggered face velocities 3
4 to the cell corners (nodes), then the velocity components are co-located and fractal interpolation for all components can be performed at the same location. However, this interpolation step is in essence an additional spatial filter of the LES velocity before fractal interpolation is performed. Numerical methods Equation (14) is solved using an unsplit geometric transport scheme for volume-of-fluid scalars that ensures both discrete volume conservation of each fluid and boundedness of the volume-of-fluid scalar ψ [18]. Geometric reconstruction of the interface within each computational cell is done using PLIC reconstruction, employing analytical formulas [19] using ELVIRA estimated normals [2, 21]. To efficiently solve Eq. (14) for the fully resolved immiscible interface, the RLSG method [12] is employed. By design, it solves the interface capturing advection equation on a separate, highly resolved Cartesian overset grid of mesh spacing h G, independent of the underlying LES flow solver grid of mesh spacing h. In the dual scale LES approach, h G needs to be chosen sufficiently small to maintain a fully resolved realization of the phase interface. The unsplit, geometric advection scheme of [18] requires face-centered velocities that are discretely divergence-free to ensure both conservation and boundedness. While discretely divergence free filtered velocities u are available on the flow solver mesh due to the application of a projection step of the velocities in a standard fractional step method, such velocities u are not directly available on the fine overset mesh. Since u = u + u, both u and u need to be discretely divergence free on the fine overset mesh. To ensure hg u = if h u =, we employ the optimal constrained approach for divergence-free velocity interpolation [22], recursively applying the interpolation technique for nested staggered meshes up to the refinement level of the overset mesh. The sub-filter fluctuation velocity u is calculated from u = WF k Iu u, (2) where the superscript k indicates k-times application of the fractal interpolation operator W F I, with ( ) h log h G k =. (21) log(2) u is then projected into the subspace of solenoidal velocity fields using the projection/correction step of a standard fractional step method. Although the fractal interpolation can be applied to the temporally evolving LES velocity u to obtain a time de- Figure 1. PLIC geometry on LES mesh cell (left), fully resolved PLIC geometry (center), difference used to calculate C mix (right). pendent u, here, we use a single snapshot of the LES field only (the initial velocity field), and thus u is frozen in time. Finally, to calculate ψ, Eq. (13) is evaluated by setting the filter size to the local flow solver mesh spacing h and evaluating the integral by explicitly summing the VoF scalar ψ of those overset-mesh cells that are contained within a given LES flow solver cell. Comparison Metrics Assuming that the phase interface is initially planar with normal orientation in the y-direction, α(x, z, t) defines the liquid volume fraction that is contained within a square column normal to the planar interface and with cross sectional area equal to h 2 G. The corresponding quantity on the LES mesh, α(x, z, t), can be calculated using a column with cross sectional area h 2. Then, the sub-filter liquid column height RMS can be defined as α 1 (t) = L 2 ((α(x, z, t) α(x, z, t)) 2 dxdz. L L (22) To take into account that even on an LES mesh scale the PLIC reconstruction of the filtered phase interface geometry provides a level of sub-filter geometry resolution beyond the mere liquid column height α, we define a local quantity C mix that measures the sub-filter variation of the phase interface on the overset mesh compared to the LES mesh. C mix is calculated by first performing a PLIC reconstruction on the LES mesh using the filtered volume fraction ψ, see Fig. 1. The reconstructed planar representation of the phase interface is then used to calculate liquid volume fractions on the overset mesh ψ hg. The root mean square of the difference of these volume fractions and the fully resolved volume frac- 4
5 Method interpolation scaling factors velocity location 3D-ZE1-N 3D ZE1 node 1 8 3D-P-N 3D Poisson node 3D-P-F 3D Poisson face 1D-P-F 1D Poisson face 1 6 Table 1. Fractal interpolation methods used. E(k) k Figure 2. Comparison of LES filtered velocity (left) and fractal interpolated velocity using method 3D- ZE1-N. tions ψ, see Fig. 1, is C mix, 1 ( ) 2 C mix (t) = ψ L 3 ψhg dxdydz. L L L (23) It should be pointed out that without a dual scale model, both α and C mix remain zero for all time. Results An initially flat interface is placed inside a box of fully developed isotropic turbulence. Both density and viscosity ratio are unity, and no surface tension forces are present with a Reynolds number of Re λ = 313 and W e =. Direct numerical simulation results using a mesh for this case are reported in [23]. Here, we present LES results using the dual-scale approach employing a LES mesh resolution of 32 3 and an overset mesh resolution of To study the impact of the sub-filter mesh resolution, an overset mesh resolution of is also studied. The sub-filter velocity u is obtained from fractal interpolation using the combination of approaches summarized in Table 1. Figure 2 shows as an example a comparison of the LES filtered velocity field and the fractal interpolated velocity using method 3D-ZE1-N. Noticeably more small scale structure is visible in the fractal interpolated velocity compared to the filter scale LES velocity. This is corroborated by the kinetic energy Figure 3. Kinetic energy spectrum for overset mesh using only trilinear interpolation (blue), fractal interpolation 3D-ZE1-N (red), sub-filter velocity u only (green), and 5/3 reference line (dashed). spectrum shown in Fig. 3 that compares the fractal interpolated velocity spectrum to one obtained using only trilinear interpolation. Note that the fractal interpolated spectrum agrees well with the theoretical 5/3 spectrum, but does not show any viscous dissipation range. The velocity spectrum using trilinear interpolation only, shows significantly lower kinetic energy on the sub-filter scale with incorrect scaling. Figure 4 compares realizations of the phase interface geometry using fractal interpolation 3D-ZE1- N and the interface geometry that would be obtained without any dual scale model. Including the fractal interpolated sub-filter velocity in the advection velocity of the resolved realization of the interface clearly leads to significantly more surface corrugations, as expected. To enable a quantitative comparison of the different fractal interpolation methods, Figs. 5 to 8 show comparisons using the two comparison metrics α and C mix. Figure 5 shows the sub-filter liquid column height RMS α using the four methods listed in Table 1 compared to the result of the DNS [23]. No significant difference between the four fractal interpolation methods can be discerned. Compared to the DNS, the dual scale LES slightly over-predicts α. This is due to the fact that the dual scale model uses a frozen fluctuation velocity field, whereas the DNS is simulated in decaying turbulence. Interactive coupling of the dual scale approach to the time evolving LES field, to be done in future work, should eliminate this discrepancy. Figure 6 analyzes the impact of the overset 5
6 .15 ` t/t Figure 6. Impact of overset mesh resolution on α : 3D-P-N-256 (red), 3D-P-F-256 (green), 3D-P-N-512 (black), 3D-P-F-512 (blue), and DNS (symbol)..2 Cmix Figure 4. Comparison of phase interface geometry: no dual scale model (left) and fractal interpolation velocity method 3D-ZE1-N (right) t/t Figure 7. Comparison of sub-filter variation Cmix : 3D-ZE1-N-256 (black), 3D-P-N-256 (red), 3D-P-F256 (green), and 1D-P-F-256 (blue) Cmix ` t/t t/t Figure 5. Comparison of sub-filter liquid column height RMS α : 3D-ZE1-N-256 (black), 3D-P-N-256 (red), 3D-P-F-256 (green), 1D-P-F-256 (blue), and DNS-124 (symbol). Figure 8. Impact of overset mesh resolution on Cmix : 3D-P-N-256 (red), 3D-P-F-256 (green), 3DP-N-512 (black), and 3D-P-F-512 (blue). 6
7 mesh resolution, comparing results obtained using the overset mesh to the mesh. No significant differences are present, indicating that in this metric, the phase interface is well resolved even on the overset mesh. Figure 7 compares the sub-filter variation C mix for the four different fractal interpolation methods. Unlike α, small differences between the methods are discernible. The successive application of 1D interpolations results in slightly higher sub-filter variation than the 3D trilinear interpolation methods, indicating that the trilinear interpolation introduces some minor dissipation in the fractal interpolated velocities. Furthermore, using filter scale face velocities instead of interpolated node velocities also results in higher C mix. This can be explained by the fact that the interpolation to the nodes is in effect an additional filter operation reducing slightly the kinetic energy at the LES filter scale before fractal interpolation is performed. While this does not impact α, C mix seems to be sensitive enough to pick up this effect. Finally, using stretching factors obtained by either method ZE1 or the more complex log-poisson distribution seems to have no significant impact on C mix. Finally, Fig. 8 analyzes the impact of the overset mesh resolution on C mix, comparing results obtained using the overset mesh to the mesh. Unlike α, increasing the overset mesh resolution increases C mix noticeably, because a higher overset mesh resolution enables the resolution of smaller scales in the resolved realization, and hence potentially larger C mix. Unfortunately, DNS results for comparison are not yet available. However, proper model verification would require that the metric C mix converges under overset mesh refinement. This will be confirmed in future work. Conclusions A dual-scale modeling approach for phase interface dynamics in turbulent flows is presented that is based on a volume-of-fluid approach. The method uses overset high-resolution meshes to capture a resolved realization of the phase interface geometry that can be explicitly filtered to close the terms that require modeling in the filtered Navier-Stokes equations. This contribution focuses on different fractal interpolation techniques to generate sub-filter velocity fluctuations that, as shown, have a significant impact on interface metrics as compared to the case of no dual-scale model, where these metrics remain zero for all time. Results show that no differences between the different fractal interpolation methods can be discerned using the sub-filter column height RMS, however, small differences are present using the more sensitive sub-filter variation metric. Necessary future work will include a more detailed comparison to DNS data both for planar interfaces and parallel liquid sheets in homogenous isotropic turbulence. Acknowledgements The support of NASA TTT grant NNX16AB7A is gratefully acknowledged. References [1] M Gorokhovski and M Herrmann. Annu. Rev. Fluid Mech., 4(1): , 28. [2] M Gorokhovski, J Jouanguy, and A Chtab. Proc. Int. Conf. Liq. Atom. Spray Syst., 26. [3] M Gorokhovski, J Jouanguy, and A Chtab. Int. Conf. Multiph. Flow, Leipzig, Germany, 27. [4] A Vallet and R Borghi. Proc. Int. Conf. Multiph. Flow, Lyon, France, [5] A Vallet, A Burluka, and R Borghi. Atom. Sprays, 11: , 21. [6] E Labourasse, D Lacanette, A Toutant, P Lubin, S Vincent, O Lebaigue, J P Caltagirone, and P Sagaut. Int. J. Multiphase Flow, 33(1):1 39, January 27. [7] A Toutant, E Labourasse, O Lebaigue, and O Simonin. Comput. Fluids, 37(7): , August 28. [8] A Toutant, M Chandesris, D Jamet, and O Lebaigue. Int. J. Multiphase Flow, 35(12): , December 29. [9] A Toutant, M Chandesris, D Jamet, and O Lebaigue. Int. J. Multiphase Flow, 35(12): , December 29. [1] M Herrmann. Comput. Fluids, 87():92 11, October 213. [11] M Herrmann and M Gorokhovski. ICLASS 29, 11th Triennial International Annual Conference on Liquid Atomization and Spray Systems, Vail, CO, 29. [12] M Herrmann. J. Comput. Phys., 227(4): , 28. [13] M Herrmann. ICLASS-215, Tainan, Taiwan, August
8 [14] M Herrmann and M Gorokhovski. Proceedings of the 28 CTR Summer Program, pp , Stanford, CA, 28. Center for Turbulence Research, Stanford University. [15] A Scotti and C Meneveau. Physica D: Nonlinear Phenomena, 127(3-4): , March [16] Ke-Qi Ding, Zhi-Xiong Zhang, Yi-Peng Shi, and Zhen-Su She. Phys. Rev. E, 82(3):36311, September 21. [17] Z X Zhang, K Q Ding, Y P Shi, and Z S She. Physical Review E, 211. [18] Mark Owkes and Olivier Desjardins. J. Comput. Phys., 27: , August 214. [19] Ruben Scardovelli and Stéphane Zaleski. J. Comput. Phys., 164(1): , October 2. [2] J E Pilliod. An analysis of piecewise linear interface reconstruction algorithms for volume-offluid methods. Master s thesis, University of California Davis, [21] James Edward Pilliod Jr and Elbridge Gerry Puckett. J. Comput. Phys., 199(2):465 52, September 24. [22] Antonio Cervone, Sandro Manservisi, and Ruben Scardovelli. Comput. Fluids, 47(1):11 114, August 211. [23] J McCaslin and O Desjardins. Proceedings of the 214 CTR Summer Program, Stanford, CA, 214. Center for Turbulence Research, Stanford University. 8
Development 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 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 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 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 informationDirect Numerical Simulations of Gas-Liquid Flows
Direct Numerical Simulations of Gas-Liquid Flows 1 Gretar Tryggvason*; 1 Jiacai Lu; 2 Ming Ma 1 Johns Hopkins University, Baltimore, MD, USA; 2 University of Notre Dame, Notre Dame, IN, USA Introduction
More informationAn evaluation of a conservative fourth order DNS code in turbulent channel flow
Center for Turbulence Research Annual Research Briefs 2 2 An evaluation of a conservative fourth order DNS code in turbulent channel flow By Jessica Gullbrand. Motivation and objectives Direct numerical
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 informationImplementation of a symmetry-preserving discretization in Gerris
Implementation of a symmetry-preserving discretization in Gerris Daniel Fuster Cols: Pierre Sagaut, Stephane Popinet Université Pierre et Marie Curie, Institut Jean Le Rond D Alembert Introduction 10/11:
More informationDirect numerical simulation of surface tension effects on interface dynamics in turbulence
ILASS Americas 27th Annual Conference on Liquid Atomization and Spray Systems, Raleigh, NC, May 215 Direct numerical simulation of surface tension effects on interface dynamics in turbulence J.O. McCaslin
More informationModel Studies on Slag-Metal Entrainment in Gas Stirred Ladles
Model Studies on Slag-Metal Entrainment in Gas Stirred Ladles Anand Senguttuvan Supervisor Gordon A Irons 1 Approach to Simulate Slag Metal Entrainment using Computational Fluid Dynamics Introduction &
More informationModelling of turbulent flows: RANS and LES
Modelling of turbulent flows: RANS and LES Turbulenzmodelle in der Strömungsmechanik: RANS und LES Markus Uhlmann Institut für Hydromechanik Karlsruher Institut für Technologie www.ifh.kit.edu SS 2012
More informationThe effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets
ILASS Americas, 26 th Annual Conference on Liquid Atomization and Spray Systems, Portland, OR, May 204 The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization
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 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 informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationNumerical Simulation of Gas-Jet Wiping in Steel Strip Galvanizing Process
, pp. 214 220 Numerical Simulation of Gas-Jet Wiping in Steel Strip Galvanizing Process Delphine LACANETTE, Stéphane VINCENT, Eric ARQUIS and Pascal GARDIN 1) Transferts, Ecoulements, Fluides, Energétique
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
More informationImpact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions
Journal of Physics: Conference Series PAPER OPEN ACCESS Impact of numerical method on auto-ignition in a temporally evolving mixing layer at various initial conditions To cite this article: A Rosiak and
More informationBuoyancy Fluxes in a Stratified Fluid
27 Buoyancy Fluxes in a Stratified Fluid G. N. Ivey, J. Imberger and J. R. Koseff Abstract Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed
More informationA unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation
A unifying model for fluid flow and elastic solid deformation: a novel approach for fluid-structure interaction and wave propagation S. Bordère a and J.-P. Caltagirone b a. CNRS, Univ. Bordeaux, ICMCB,
More informationNumerical study of stochastic particle dispersion using One-Dimensional-Turbulence
ILSS-mericas 29th nnual Conference on Liquid tomization and Spray Systems, tlanta, G, May 217 Numerical study of stochastic particle dispersion using One-Dimensional-Turbulence Marco Fistler *1, David
More informationNumerical simulation of phase separation and a priori two-phase LES filtering
Available online at www.sciencedirect.com Computers & Fluids 37 (2008) 898 906 www.elsevier.com/locate/compfluid Numerical simulation of phase separation and a priori two-phase LES filtering S. Vincent
More informationEnsemble averaged dynamic modeling. By D. Carati 1,A.Wray 2 AND W. Cabot 3
Center for Turbulence Research Proceedings of the Summer Program 1996 237 Ensemble averaged dynamic modeling By D. Carati 1,A.Wray 2 AND W. Cabot 3 The possibility of using the information from simultaneous
More informationSpectral analysis of energy transfer in variable density, radiatively heated particle-laden flows
Center for Turbulence Research Proceedings of the Summer Program 24 27 Spectral analysis of energy transfer in variable density, radiatively heated particle-laden flows By H. Pouransari, H. Kolla, J. H.
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 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 informationLARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS
The 6th ASME-JSME Thermal Engineering Joint Conference March 6-, 3 TED-AJ3-3 LARGE EDDY SIMULATION OF MASS TRANSFER ACROSS AN AIR-WATER INTERFACE AT HIGH SCHMIDT NUMBERS Akihiko Mitsuishi, Yosuke Hasegawa,
More informationDynamics of the Coarse-Grained Vorticity
(B) Dynamics of the Coarse-Grained Vorticity See T & L, Section 3.3 Since our interest is mainly in inertial-range turbulence dynamics, we shall consider first the equations of the coarse-grained vorticity.
More informationBefore we consider two canonical turbulent flows we need a general description of turbulence.
Chapter 2 Canonical Turbulent Flows Before we consider two canonical turbulent flows we need a general description of turbulence. 2.1 A Brief Introduction to Turbulence One way of looking at turbulent
More informationToward 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 information+ = + t x x x x u. The standard Smagorinsky model has been used in the work to provide the closure for the subgridscale eddy viscosity in (2):
International Conference on Methods of Aerophysical Research, ICMAR 008 LARGE EDDY SIMULATION OF TURBULENT ROUND IMPINGING JET B.B. Ilyushin, D.V. Krasinsky Kutateladze Institute of Thermophysics SB RAS
More informationEngineering. Spring Department of Fluid Mechanics, Budapest University of Technology and Economics. Large-Eddy Simulation in Mechanical
Outline Geurts Book Department of Fluid Mechanics, Budapest University of Technology and Economics Spring 2013 Outline Outline Geurts Book 1 Geurts Book Origin This lecture is strongly based on the book:
More informationTurbulence: Basic Physics and Engineering Modeling
DEPARTMENT OF ENERGETICS Turbulence: Basic Physics and Engineering Modeling Numerical Heat Transfer Pietro Asinari, PhD Spring 2007, TOP UIC Program: The Master of Science Degree of the University of Illinois
More informationLecture 14. Turbulent Combustion. We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing.
Lecture 14 Turbulent Combustion 1 We know what a turbulent flow is, when we see it! it is characterized by disorder, vorticity and mixing. In a fluid flow, turbulence is characterized by fluctuations of
More informationModeling flame brush thickness in premixed turbulent combustion
Center for Turbulence Research Proceedings of the Summer Program 2006 299 Modeling flame brush thickness in premixed turbulent combustion By E. Knudsen, O. Kurenkov, S. Kim, M. Oberlack AND H. Pitsch Turbulent
More informationApplication of stochastic models to primary air-blast atomization and cavitation
Center for Turbulence Research Proceedings of the Summer Program 2012 325 Application of stochastic models to primary air-blast atomization and cavitation By M. Gorokhovski, R. Zamansky, C. Le Ribault
More informationIntroduction to Turbulence and Turbulence Modeling
Introduction to Turbulence and Turbulence Modeling Part I Venkat Raman The University of Texas at Austin Lecture notes based on the book Turbulent Flows by S. B. Pope Turbulent Flows Turbulent flows Commonly
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 informationMassimo GERMANO Politecnico di Torino
Hybrid Massimo GERMANO Politecnico di Torino Martín SÁNCHEZ-ROCHA Dassault Systèmes SIMULIA Corporation Suresh MENON Georgia Institute of Technology 64th Annual APS-DFD Meeting Baltimore, Maryland November
More informationCVS filtering to study turbulent mixing
CVS filtering to study turbulent mixing Marie Farge, LMD-CNRS, ENS, Paris Kai Schneider, CMI, Université de Provence, Marseille Carsten Beta, LMD-CNRS, ENS, Paris Jori Ruppert-Felsot, LMD-CNRS, ENS, Paris
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationTurbulence Models for Flows with Free Surfaces and Interfaces
AIAA JOURNAL Vol. 44, No. 7, July 2006 Turbulence Models for Flows with Free Surfaces and Interfaces Ebrahim Shirani, Ali Jafari, and Nasser Ashgriz University of Toronto, Toronto, Ontario M5S 3G8, Canada
More informationIMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE FLUCTUATIONS
Proceedings of the Combustion Institute, Volume 9, 00/pp. 1971 1978 IMPROVED POLLUTANT PREDICTIONS IN LARGE-EDDY SIMULATIONS OF TURBULENT NON-PREMIXED COMBUSTION BY CONSIDERING SCALAR DISSIPATION RATE
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 informationThree-dimensional wall filtering formulation for large-eddy simulation
Center for Turbulence Research Annual Research Briefs 6 55 Three-dimensional wall filtering formulation for large-eddy simulation By M. Shoeybi AND J. A. Templeton 1. Motivation and objectives Large-eddy
More informationComputational Fluid Dynamics 2
Seite 1 Introduction Computational Fluid Dynamics 11.07.2016 Computational Fluid Dynamics 2 Turbulence effects and Particle transport Martin Pietsch Computational Biomechanics Summer Term 2016 Seite 2
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More informationA Solution Method for the Reynolds-Averaged Navier-Stokes Equation
A Solution Method for the Reynolds-Averaged Navier-Stokes Equation T.-W. Lee Mechanical and Aerospace Engineering, Arizona State University, Tempe, AZ, 8587 Abstract- Navier-Stokes equations are difficult
More informationInteraction(s) fluide-structure & modélisation de la turbulence
Interaction(s) fluide-structure & modélisation de la turbulence Pierre Sagaut Institut Jean Le Rond d Alembert Université Pierre et Marie Curie- Paris 6, France http://www.ida.upmc.fr/~sagaut GDR Turbulence
More informationNumerical Methods in Aerodynamics. Turbulence Modeling. Lecture 5: Turbulence modeling
Turbulence Modeling Niels N. Sørensen Professor MSO, Ph.D. Department of Civil Engineering, Alborg University & Wind Energy Department, Risø National Laboratory Technical University of Denmark 1 Outline
More informationTurbulence Modeling I!
Outline! Turbulence Modeling I! Grétar Tryggvason! Spring 2010! Why turbulence modeling! Reynolds Averaged Numerical Simulations! Zero and One equation models! Two equations models! Model predictions!
More informationNumerical simulations of drop impacts
Numerical simulations of drop impacts Christophe Josserand Institut D Alembert, CNRS-UPMC L. Duchemin, Z. Jian, P. Ray and S. Zaleski Numerical simulations of drop impacts Christophe Josserand Institut
More informationOptimizing calculation costs of tubulent flows with RANS/LES methods
Optimizing calculation costs of tubulent flows with RANS/LES methods Investigation in separated flows C. Friess, R. Manceau Dpt. Fluid Flow, Heat Transfer, Combustion Institute PPrime, CNRS University
More informationThe Johns Hopkins Turbulence Databases (JHTDB)
The Johns Hopkins Turbulence Databases (JHTDB) HOMOGENEOUS BUOYANCY DRIVEN TURBULENCE DATA SET Data provenance: D. Livescu 1 Database Ingest and Web Services: C. Canada 1, K. Kalin 2, R. Burns 2 & IDIES
More informationOn the transient modelling of impinging jets heat transfer. A practical approach
Turbulence, Heat and Mass Transfer 7 2012 Begell House, Inc. On the transient modelling of impinging jets heat transfer. A practical approach M. Bovo 1,2 and L. Davidson 1 1 Dept. of Applied Mechanics,
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationDNS of the Taylor-Green vortex at Re=1600
DNS of the Taylor-Green vortex at Re=1600 Koen Hillewaert, Cenaero Corentin Carton de Wiart, NASA Ames koen.hillewaert@cenaero.be, corentin.carton@cenaero.be Introduction This problem is aimed at testing
More informationIntermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint
Intermittency of quasi-static magnetohydrodynamic turbulence: A wavelet viewpoint Naoya Okamoto 1, Katsunori Yoshimatsu 2, Kai Schneider 3 and Marie Farge 4 1 Center for Computational Science, Nagoya University,
More informationA G-equation formulation for large-eddy simulation of premixed turbulent combustion
Center for Turbulence Research Annual Research Briefs 2002 3 A G-equation formulation for large-eddy simulation of premixed turbulent combustion By H. Pitsch 1. Motivation and objectives Premixed turbulent
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
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 informationFragmentation under the scaling symmetry and turbulent cascade with intermittency
Center for Turbulence Research Annual Research Briefs 23 197 Fragmentation under the scaling symmetry and turbulent cascade with intermittency By M. Gorokhovski 1. Motivation and objectives Fragmentation
More informationRegularization modeling of turbulent mixing; sweeping the scales
Regularization modeling of turbulent mixing; sweeping the scales Bernard J. Geurts Multiscale Modeling and Simulation (Twente) Anisotropic Turbulence (Eindhoven) D 2 HFest, July 22-28, 2007 Turbulence
More informationDNS, LES, and wall-modeled LES of separating flow over periodic hills
Center for Turbulence Research Proceedings of the Summer Program 4 47 DNS, LES, and wall-modeled LES of separating flow over periodic hills By P. Balakumar, G. I. Park AND B. Pierce Separating flow in
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationProbability density function (PDF) methods 1,2 belong to the broader family of statistical approaches
Joint probability density function modeling of velocity and scalar in turbulence with unstructured grids arxiv:6.59v [physics.flu-dyn] Jun J. Bakosi, P. Franzese and Z. Boybeyi George Mason University,
More informationA DNS database to improve atomization modeling
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 A DNS database to improve atomization modeling B. Duret, T. Menard,
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 informationJournal of Fluid Science and Technology
Science and Technology Effect of Molecular Diffusivities on Countergradient Scalar Transfer in a Strong Stable Stratified Flow (Study on the Linear and Nonlinear Processes by using RDT) Kouji NAGATA, Takashi
More informationZonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach
Zonal hybrid RANS-LES modeling using a Low-Reynolds-Number k ω approach S. Arvidson 1,2, L. Davidson 1, S.-H. Peng 1,3 1 Chalmers University of Technology 2 SAAB AB, Aeronautics 3 FOI, Swedish Defence
More informationA high resolution collision algorithm for anisotropic particle populations
A high resolution collision algorithm for anisotropic particle populations Philipp Pischke 1 and Reinhold Kneer1 1 Institute of Heat and Mass Transfer, RWTH Aachen University July 2014 Collision algorithm,
More informationDatabase analysis of errors in large-eddy simulation
PHYSICS OF FLUIDS VOLUME 15, NUMBER 9 SEPTEMBER 2003 Johan Meyers a) Department of Mechanical Engineering, Katholieke Universiteit Leuven, Celestijnenlaan 300, B3001 Leuven, Belgium Bernard J. Geurts b)
More informationBasic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue
11 Basic Features of the Fluid Dynamics Simulation Software FrontFlow/Blue Yang GUO*, Chisachi KATO** and Yoshinobu YAMADE*** 1 FrontFlow/Blue 1) is a general-purpose finite element program that calculates
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 information3D Numerical Simulation of Supercritical Flow in Bends of Channel
3D Numerical Simulation of Supercritical Flow in Bends of Channel Masoud. Montazeri-Namin, Reyhaneh-Sadat. Ghazanfari-Hashemi, and Mahnaz. Ghaeini- Hessaroeyeh Abstract An attempt has been made to simulate
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 informationDirect Numerical Simulation of fractal-generated turbulence
Direct Numerical Simulation of fractal-generated turbulence S. Laizet and J.C. Vassilicos Turbulence, Mixing and Flow Control Group, Department of Aeronautics and Institute for Mathematical Sciences, Imperial
More informationOn the decay of two-dimensional homogeneous turbulence
On the decay of two-dimensional homogeneous turbulence J. R. Chasnov a) The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong Received 19 June 1996; accepted 3 September
More informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
More informationA NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS
June 3 - July 3, 25 Melbourne, Australia 9 B-3 A NEW AUTONOMIC CLOSURE FOR LARGE EDDY SIMULATIONS Ryan N. King Department of Mechanical Engineering University of Colorado Boulder, CO, 839, USA ryan.n.king@colorado.edu
More informationTurbulence models and excitation of solar oscillation modes
Center for Turbulence Research Annual Research Briefs Turbulence models and excitation of solar oscillation modes By L. Jacoutot, A. Wray, A. G. Kosovichev AND N. N. Mansour. Motivation and objectives
More informationDynamic k-equation Model for Large Eddy Simulation of Compressible Flows. Xiaochuan Chai and Krishnan Mahesh
40th Fluid Dynamics Conference and Exhibit 8 June - July 00, Chicago, Illinois AIAA 00-506 Dynamic k-equation Model for Large Eddy Simulation of Compressible Flows Xiaochuan Chai and Krishnan Mahesh University
More informationmeters, we can re-arrange this expression to give
Turbulence When the Reynolds number becomes sufficiently large, the non-linear term (u ) u in the momentum equation inevitably becomes comparable to other important terms and the flow becomes more complicated.
More informationJournal of Computational Physics 157, (2000) doi: /jcph , available online at
Journal of Computational Physics 157, 796 800 000) doi:10.1006/jcph.1999.639, available online at http://www.idealibrary.com on NOTE A Simple Method for Computing Far-Field Sound in Aeroacoustic Computations
More informationEulerian models. 2.1 Basic equations
2 Eulerian models In this chapter we give a short overview of the Eulerian techniques for modelling turbulent flows, transport and chemical reactions. We first present the basic Eulerian equations describing
More informationBeyond the Point Particle: LES-Style Filtering of Finite-Sized Particles
ILASS Americas th Annual Conference on Liquid Atomization and Spray Systems, Chicago, IL, May 7 Beyond the Point Particle: LES-Style Filtering of Finite-Sized Particles Brooks Moses and Chris Edwards Department
More informationLES of turbulent shear flow and pressure driven flow on shallow continental shelves.
LES of turbulent shear flow and pressure driven flow on shallow continental shelves. Guillaume Martinat,CCPO - Old Dominion University Chester Grosch, CCPO - Old Dominion University Ying Xu, Michigan State
More informationMulti-Scale Modeling of Turbulence and Microphysics in Clouds. Steven K. Krueger University of Utah
Multi-Scale Modeling of Turbulence and Microphysics in Clouds Steven K. Krueger University of Utah 10,000 km Scales of Atmospheric Motion 1000 km 100 km 10 km 1 km 100 m 10 m 1 m 100 mm 10 mm 1 mm Planetary
More informationComputational Analysis of an Imploding Gas:
1/ 31 Direct Numerical Simulation of Navier-Stokes Equations Stephen Voelkel University of Notre Dame October 19, 2011 2/ 31 Acknowledges Christopher M. Romick, Ph.D. Student, U. Notre Dame Dr. Joseph
More informationNEAR-WALL TURBULENCE-BUBBLES INTERACTIONS IN A CHANNEL FLOW AT Re =400: A DNS/LES INVESTIGATION
ABSTRACT NEAR-WALL TURBULENCE-BUBBLES INTERACTIONS IN A CHANNEL FLOW AT Re =400: A DNS/LES INVESTIGATION D. Métrailler, S. Reboux and D. Lakehal ASCOMP GmbH Zurich, Technoparkstr. 1, Switzerland Metrailler@ascomp.ch;
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 informationOn the (multi)scale nature of fluid turbulence
On the (multi)scale nature of fluid turbulence Kolmogorov axiomatics Laurent Chevillard Laboratoire de Physique, ENS Lyon, CNRS, France Laurent Chevillard, Laboratoire de Physique, ENS Lyon, CNRS, France
More informationTurbulence Modeling. Cuong Nguyen November 05, The incompressible Navier-Stokes equations in conservation form are u i x i
Turbulence Modeling Cuong Nguyen November 05, 2005 1 Incompressible Case 1.1 Reynolds-averaged Navier-Stokes equations The incompressible Navier-Stokes equations in conservation form are u i x i = 0 (1)
More informationWALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID
WALL RESOLUTION STUDY FOR DIRECT NUMERICAL SIMULATION OF TURBULENT CHANNEL FLOW USING A MULTIDOMAIN CHEBYSHEV GRID Zia Ghiasi sghias@uic.edu Dongru Li dli@uic.edu Jonathan Komperda jonk@uic.edu Farzad
More informationEffects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence
Effects of Forcing Scheme on the Flow and the Relative Motion of Inertial Particles in DNS of Isotropic Turbulence Rohit Dhariwal PI: Sarma L. Rani Department of Mechanical and Aerospace Engineering The
More informationFinal abstract for ONERA Taylor-Green DG participation
1st International Workshop On High-Order CFD Methods January 7-8, 2012 at the 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee Final abstract for ONERA Taylor-Green DG participation JB Chapelier,
More information