A GAS-KINETIC SCHEME FOR TURBULENT FLOW
|
|
- Roy Robinson
- 6 years ago
- Views:
Transcription
1 A GAS-KINETIC SCHEME FOR TURBULENT FLOW M. Righi Zurich University of Applied Sciences Abstract A gas-kinetic scheme is used in this study to simulate turbulent flow and the interaction between a shock wave and a turbulent boundary layer in particular. Gas-kinetic schemes belong to a class of their own of numerical schemes for fluid mechanics. These schemes are neither particle methods, such as DSMC and Molecular Dynamics, nor Lattice Boltzmann methods: they calculate the fluxes between numerical cells on the basis of the Boltzmann equation instead of the Navier-Stokes or Euler equations. Advective and viscous fluxes are coupled in gas-kinetic schemes; in other words, the effect of unresolved fluctuations (represented in conventional schemes by the viscous fluxes) may affect the transport terms (advection). In simulations of turbulent flow where the unresolved turbulent fluctuations are modeled by eddy viscosity, a gas-kinetic scheme may exploit the information provided by the turbulence model in a different way. The analysis proposed in this paper shows that gas-kinetic schemes have an in-built multiscalar sensor which controls the correction terms generated by the underlying kinetic theory. The simulations shown in this paper have been obtained with the RANS approach and a standard two-equation turbulence model (k-ω). It is shown that the gas-kinetic scheme provides good quality predictions in a number of test cases, where conventional schemes with the same turbulence model are known to fail. Introduction Turbulence is clearly a multiscale problem. In numerical simulations, spatial and temporal resolution may be insufficient to resolve all turbulent scales of motion down to the dissipative range. The effects of the unresolved scales of motion on the resolved ones must then be modeled. This approach has been inspired by the handling of thermal fluctuations, which can virtually never be resolved in the continuum regime, which are invariably modeled by diffusion proportional to a given viscosity. However, unlike thermal fluctuations, unresolved turbulent fluctuations may not be small and spread over a range without a clear separation from the resolved scales of motion. In theory, this situation, which is even a pre-requisite in Large Eddy Simulation, may make conventional advection-diffusion schemes physically inconsistent. A similar boundary is reached when (laminar) rarefied flows are simulated with conventional Navier-Stokes schemes. In rarefied flows, it may happen that even the scale of thermal fluctuations become not negligible with respect to the resolved scales. The unresolved-to-resolved scales ratio assumes large values in special flow regions, such as shocklayers, due to presence of strong gradients which may lead to small resolved scales and a high turbulence intensity, which may be related to large (in space and time) turbulent scales of motion. It is well accepted that conventional schemes, with the RANS approach in particular, may fail to accurately capture details of the shock-boundary layer interaction such as strength and position of the shock, position of separation and re-attachment, turbulence intensity downstream of the separated flow region. A comprehensive review can be found in [Babinsky and Harvey, 0]. In many cases, such as the design of supersonic or hypersonic vehicles, inaccurate prediction of this area may lead to unsuitable geometries or high development costs. In this study a gas-kinetic scheme is used, following the RANS approach, to simulate turbulent flow. Gas-kinetic schemes are derived form the Boltzmann equation and not from the Navier-Stokes or Euler equations. They are suitable for finite volume or finite elements schemes but, being relatively little known, they are often confused with particle methods such as DSMC or Molecular Dynamics or with Lattice Boltzmann methods. Numerous gas-kinetic schemes have been developed over the latest twenty years [Xu, 00, May et al., 007, Mandal and Deshpande, 99, Chou and Baganoff, 997, Xu and Prendergast, 99] with the aim to achieve a physically more consistent mathematical model of fluid mechanics. Gas-kinetic schemes are more accurate than conventional schemes, and might be able to resolve shock-layers. Besides, they are more suitable to high-order reconstruction [Li et al., 00, Xu et al., 005, Xuan and Xu, 0] and may be used as a platform to investigate rarefied flow [Liao et al., 007, Xu and Huang, 00]. Since conservation laws such as the Navier-Stokes and Euler equations can be derived from the Boltzmann equation, as shown in [Cercignani, 988, Xu, 998], a gas-kinetic scheme is always consistent with a conventional one, within the validity boundaries of the latter. Simulation of rarefied flow implies some multiscalar mechanism, i.e. the ability to correct the
2 conventional advective and diffusive fluxes whenever the unresolved scales of motion are not negligibly small. This property might prove interesting in the simulation of turbulent flow. To build the turbulent gas-kinetic scheme shown in this paper, the scheme developed by [Xu, 00] has been the starting point. An allied two-equation turbulence model (k-ω) is solved alongside the equations for the conservative variables (with Boltzmann-BGK fluxes). The strength of the scheme lies in its capability to evaluate advective and diffusive fluxes in a single operation, including therefore the effect of collisions on transport, or, in other words, to take into account the scale of the unresolved fluctuations - introducing a multiscale effect which is missing in Navier-Stokes schemes. Note that the inclusion of kinetic effects also generates a misalignment between turbulent stress and strain rate acting in a similar way to high-order turbulence models, as remarked in [Chen et al., 00]. The structure of this draft paper includes the description of the GKS used and modified for turbulent flows in section, the presentation of results obtained in numerical experiments and conclusions. Gas-kinetic schemes The state of a gas can be described by means of a distribution function f(x, v, t) - defined in the phase space. The conservative variables w = [ρ ρv ρv ρv ρe] T are recovered from f by taking moments in the phase space: w = ψfdξ, () where the elementary volume in phase space is dξ = dv dv dv dξ and: ψ = [ ( v v v ui + ξ )] T. () The numerical fluxes F related to a unit interface length normal to direction n, and a time step t are obtained by time integration: F n = t 0 fψu n dξ dt. () The distribution function f is assumed to be a solution of the Boltzmann-BGK equation [Bhatnagar et al., 95]: f t + (u )f = f eq f, () τ where τ is a relaxation time corresponding to the average fluctuation period and taken to be proportional to molecular viscosity and inversely proportional to static pressure in laminar flow, τ = µ p. A distribution function consistent with Eq. is obtained by means of the Chapman-Enskog expansion [Cercignani, 988]: f = f eq ɛ τdf eq + ɛ τd ( τdf eq ) (5) where τ is a reference relaxation time and ɛ is a small number representing the deviations from thermal equilibrium. It can be demonstrated [Cercignani, 988, Xu, 998] that Eq. 5 truncated at the first order ignores thermal fluctuations and is consistent with Euler fluxes. Eq. 5 truncated at the second order would generate Navier-Stokes fluxes whereas the next two truncation orders would provide the Burnett and super Burnett equations. The scheme described in this paper is a modification of the one developed in [Xu, 00] for laminar flow. It uses the MUSCL approach and assumes as such a conventional discontinuous reconstruction. The gas evolution is integrated in time at each time step. An initial distribution f 0 is generated on the basis of the gas states and adopting a second-order Chapman-Enskog expansion. Further, a distribution function f BGK is obtained as the analytical solution to the BGK equation, Eq. : τ f BGK (x, x, x, t, u, u, u, ξ) = t o f eq (x, x, x, t, u, u, u, ξ)e (t t )/τ dt +e t/τ f 0 (x u t, x u t, x u t) (6), where x = x u (t t ), x = x u (t t ), x = x u (t t ). The initial distribution f 0 is rather complex as it includes the spatial and temporal terms, necessary to obtain a proper Chapman-Enskog expansion. Details of the derivation can be obtained in [Xu, 00, Righi, 0, Righi, 0]. In this paper, the resulting distribution function, solution to the BGKmodel and compatible with the gas initial states, f BGK, is expressed in a compact form f BGK = f NS + e t/τ f, (7) where f NS is a second-order Chapman-Enskog expansion built around an average gas state at the interface, whereas f is a more complex function, obtained as difference between f NS and a Chapman- Enskog expansion obtained combining the left and right gas states. The presence of the blending function e t/τ is important and represents the kernel of the multiscalarity of the gas-kinetic scheme. It is important to remark that f NS generates fluxes which might be compared to the ones provided by a central conventional scheme. At first sight, Eq. 7, has the typical structure of upwind schemes: central fluxes plus a correction term whose role becomes more important the stronger the gradients involved. The fluxes are finally obtained inserting f = f BGK into Eq. ; they can be expressed into the form: F n = F NS n + α(ɛ) F n, (8)
3 where Fn NS are the fluxes associated with f NS and Fn the ones associated with f, α(ɛ) = ɛ( e /ɛ )) and ɛ = τ/ t. In Eq. 8 the parameter ɛ plays a very important role, as it controls the correction kinetic terms. Since ɛ can be seen as the ratio of unresolved to resolved timescales, the mechanism which adjusts the scheme to the size of unresolved scales is clearly visible. It is sometimes referred to as intrinsic multiscalarity of this scheme. This property has been exploited in laminar flow; since this scheme has also been used in the transitional regime [Liao et al., 007] with a suitable definition of τ. As is well known, the BGK model is consistent with conservationlaws with an unity Prandtl number; the heath flux must therefore be corrected for realistic fluids. A number of publications ([Li et al., 00, Xu et al., 005, Xuan and Xu, 0] and references therein) demonstrate the good qualities of the scheme. Following the RANS approach, the effect of unresolved turbulence can be modeled through an eddy viscosity µ t, determined by the modeled turbulent quantities. In a gas-kinetic scheme, the relaxation time τ must be related to the eddy viscosity. Trivially, one could set τ = µ t /p in analogy with laminar flow. However, the considerations exposed in [Chen et al., 00, Chen et al., 00] lead to: k/ω τ = µ p +, (9) T ( + η / ) where T is temperature and η = Sk/ε, S is a measure of the local velocity gradient. The turbulent quantities k and ω are provided by a standard k-ω model solved alongside. In turbulent flow, the quantity ɛ = τ/ t may assume large values, leading to numerical problems. It is proposed to replace it with: ɛ =, (0) τ τ where τ is the assumed timescale of the resolved flow and can be calculated form the gradients of one of the resolved variables (i.e. τ = ρ/dρ). Details are explained in [Righi, 0]. The expression of ɛ used in Eq. 0 is close to the one used in rarefied gas dynamics to estimate the local Knudsen number, with the difference that in this case timescales are used instead of spatial quantities. In the section dedicated to the results obtained, it can be seen that this timescale ratio or degree of rarefaction may assume values up to a few hundredths inside shocklayers. In rarefied gas dynamics the validity of Navier-Stokes schemes may become questionable for Knudsen number as high as one thousandth. The use of a gas-kinetic scheme with eddy viscosity may therefore exploit its better ability to cope with these rarefaction levels. Numerical implementation The gas-kinetic scheme described above has been implemented into a two-dimensional finite-volume steady-state solver. The results shown in this paper have been obtained with a second-order reconstruction. The minmod limiter has been used for the reconstruction of conservative variables and their gradients. The gas-kinetic scheme has been used in this study with CFL between and 0, which is slightly higher than the value used with the same solver with conventional fluxes. This not withstanding, the overall computational cost requested by the gas-kinetic scheme is higher by a factor of two on average. The solver is sequential and no attempts have been done to optimize the code for speed. No wall functions have been used, the viscous sublayer has been accurately resolved in every flow case. Results and Discussion In this paper results from four test cases are shown: the flow around the RAE8 supercritical airfoil at the conditions of Case 0 in the experiments published in [Cook et al., 979], the supersonic compression ramp investigated in [Settles et al., 979] and the one at lower Reynolds number investigate in [Bookey et al., 005], the hypersonic compression ramp investigated in [Dolling and Erengil, 99]. Transonic flow around a RAE 8 airfoil in Case 0 In Case 0 the airfoil is immersed in a flow at Mach M = 0.75 at an angle of attack α =.9 and Reynolds number Re c = based on airfoil chord. The experimental investigation have highlighted an incipient separation of the boundary layer, which fails to be predicted by most conventional solvers with a linear twoequation turbulence model (refer to the discussion in [Wallin and Johansson, 000]). More sophisticated turbulence models such as algebraic stress models may capture the incipient separation correctly but the calibration of all coefficients is seldom suitable for many flow classes. Fig. shows the results obtained with the gas-kinetic scheme and the ones obtained with a conventional solver based on MUSCL reconstruction and Roe s approximate Riemann solver. Remarkably, the gas-kinetic scheme provides a rather accurate prediction, whereas the conventional Navier-Stokes scheme places the shock slightly downstream of the position indicted by the measurements. Supersonic compression corner at high and low Reynolds number The supersonic flow impinging on a compression corner has been studied by many investigators in the past. One of the most popular campaign is the one carried out at Princeton and published in [Settles et al., 976, Settles et al., 979]. The flow is characterized by a Mach M =.85 and a Reynolds number Re θ = 000, based on momentum thick-
4 Cp Cf (a) x/c (b) x/c Figure : RAE8 airfoil (Case 0 Re = , M = 0.75, angle of attack α =.9 ). ( ) Gas-kinetic scheme (GKS) on finest grid, ( ) GKS on medium grid, ( ) GKS on coarsest grid, ( ) Navier-Stokes (Roe s approximate Riemann solver) on finest grid, ( o ): experimental data from Cook [Cook et al., 979]. (a) pressure coefficient, (b) skin friction coefficient. ness. Four different corners have been used: 8, 6, 0 and. The 8 corner does not separate the flow and the 6 one generates only an incipient separation. The time-averaged separation at 0 spans about δ and the one at about δ (where δ is the incoming boundary layer thickness). Results from conventional Navier-Stokes scheme are not reported in this paper, as for these flow cases, they tend to become more dependent on numerical details such as reconstruction technique and factorization. In order to avoid proposing one particular implementation, the reader is referred to the literature, for instance on the comprehensive review in [Babinsky and Harvey, 0], [Goldberg et al., 998] and [Menter and Rumsey, 99]. As is well accepted, conventional schemes fail to predict the right shock position and separation length: the whole interaction region is often translated downstream, the extension of the separation is often underestimated. The turbulent gas-kinetic scheme developed in this study does not provide a perfect agreement with experiments but does better than conventional schemes. Figures and show the reasonably good agreement of predictions in all four cases in terms of wall pressure and skin friction coefficient. In order to generate data for a benchmark case for LES and DNS methods, the same flow, at much a lower Reynolds number of Re θ = 00, has been investigated many years later also in Princeton and published in [Bookey et al., 005]. The first experiment of this campaign concerns a supersonic compression corner with an angle of. Within the same investigation an shock impinging with an angle of onto a turbulent boundary layer has also been investigated. The two flow generates a similar shock-boundary layer interaction. Fig. shows the pressure distributions for both the compression corner and the impinging shock, which are in reasonably good agreement with the experiments. For both flow cases, the one investigated by Settles and the one investigated by Bookey, the sensitivity to the Reynolds number has been investigated with the gas-kinetic scheme. Fig. 5 shows the interaction length as well as the separation point as a function of the Reynolds number. Both quantities are compared to the values estimated with an empirical law devised by Settles, details can be found in [Délery et al., 986]. The prediction is very accurate at lower Reynolds number and barely acceptable in the other case. Moreover, the effect of wall temperature on the interaction has also been investigated with the gaskinetic scheme. In this case the experimental data published in [Spaid and Frishett, 97] has been used. Fig. 6 shows the separation length and the separation
5 α = 8.0 Cf 0 α = 8 α = 6.0 Cf 0 α = 6 α = 0.0 Cf 0 α = 0 α =.0 Cf 0 α = Figure : Pressure calculated for four different compression corner flows, characterized by angles values of 8, 6, 0 and (freestream conditions: M =.85, Re = per length unit, δ 0 = m). ( ) Gas-kinetic scheme (GKS) on finest grid, ( ) GKS on medium grid (only shown for α = ), ( ) GKS on coarsest grid (only shown for α = ), ( o ): experimental data from Settles [Settles et al., 979]. Figure : Skin friction calculated for four different compression corner flows, characterized by angles values of 8, 6, 0 and (freestream conditions: M =.85, Re = per length unit, δ 0 = m). ( ) Gas-kinetic scheme (GKS) on finest grid, ( ) GKS on medium grid (only shown for α = ), ( ) GKS on coarsest grid (only shown for α = ), ( o ): experimental data from Settles [Settles et al., 979].
6 5 (a) L0/δ e + 06 e + 06 e + 06 Re.00 (b).50 L0/δ e + 07 e + 08 e + 08 Re Figure : Distribution of static pressure calculated for the compression corner (a) with an angle of and for the reflected shock (b) originating form a compression corner at. Freestream conditions (both flows) are M =.90, Re θ = 00. ( ) GKS on grid, ( ) GKS on grid, ( o ): experimental data from Bookey [Bookey et al., 005]. Grid sizes compression corner: 8 9, 5 68 respectively, grid size impinging shock: 8 08, 96 0 respectively. Grid and have different resolution and have been generated with different algorithms. point as a function of the wall temperature. The prediction can be considered as acceptable. The distribution of the time scales ratio ɛ is shown in the Figures 7 and 8. In both cases, it is evident that ɛ assumes the larger values in the interaction area, that is, where the results obtained with the gas-kinetic scheme and conventional schemes differ most. Hypersonic compression corner at Mach.95 This flow case, investigated by Dolling et al. [Dolling and Erengil, 99], extends the Mach number range to the hypersonics. However, adiabatic wall conditions are assumed. Reynolds number is based on boundary layer thickness. Results from RANS and hybrid simulations can be found in Edwards et al. [Edwards, 008]. In Fig. 9 the pressure distribution predicted by the turbulent gaskinetic scheme is compared to the experimental values [Dolling and Erengil, 99], showing an acceptable agreement. In Fig. 9 the distribution of the time scales ratio ɛ is shown. It is evident that the highest Figure 5: Interaction length as a function of Reynolds number; the prediction obtained with the gas-kinetic scheme are compared to the empirical law by Settles from [Délery et al., 986]. values arise in the interaction area, that is where the results provided by the gas-kinetic scheme and conventional schemes differ most. 5 Conclusions In this paper a gas-kinetic scheme is used to simulate turbulent flow, following the RANS approach with a standard two-equation model. The test cases proposed contain strong interactions between a shock and the turbulent boundary layer, causing a substantial separation. It is shown that the results obtained with gaskinetic scheme - unlike the ones, available in the literature, obtained with conventional schemes - are in good or very good agreement with the experimental data. The rationale is put forward, that the higher accuracy of the scheme proposed originates form the multiscalar mechanisms, a legacy of the kinetic theory of gases, which corrects the numerical fluxes as a function of the unresolved-to-resolved scales ratio. The idea of eddy viscosity assumes that the unresolved scales of motion (turbulent fluctuations) are small and well separated from the smallest resolved scales. As is well known, this condition is not always met in and around shocklayers, where the resolved scales are small and turbulence level is high. In theory, this situation might make conventional schemes unsuitable; on the other hand, it might provide a clear advantage for gas-kinetic schemes.
7 xm/δ T w /T t xs/δ T w /T t Figure 6: Effects of wall temperature for a 0 angle ramp, experimental data from [Spaid and Frishett, 97]. T w is wall temperature, T t is stagnation temperature, x m is the starting point of interaction from the corner, x s is the starting point of separation from the corner. Figure 9: Compression corner M =.95, Experimental values from Dolling et al. [Dolling and Erengil, 99]. Results obtained on a grid with size: 8 60 The behavior of gas-kinetic schemes in turbulent simulations needs additional validation, in particular in three dimensional geometries. However, the results obtained so far are very encouraging; two-dimensional test cases are proposed ranging from transonic to hypersonic (M = 5) where the gas-kinetic scheme predicts position and size of the interaction in an acceptable way. It is worthwhile observing that the same approach might be suitable for Large Eddy Simulation too, as the cut-off scale in LES is by definition in the inertial subrange, i.e. between adjacent smallest-resolved and largest-unresolved scales. 6 References Figure 7: Compression corner at higher Reynolds number investigated in [Settles et al., 976], distribution of the time scales ratio ɛ. [Babinsky and Harvey, 0] Babinsky, H. and Harvey, J. K. (0). Shock wave-boundary-layer interactions. Cambridge University Press New York. [Bhatnagar et al., 95] Bhatnagar, P., Gross, E., and Krook, M. (95). A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Phys. Rev., 9():5 55. Figure 8: Compression corner at lower Reynolds number investigated in [Bookey et al., 005], distribution of the time scales ratio ɛ. [Bookey et al., 005] Bookey, P., Wyckham, C., and Smits, A. (005). Experimental investigations of Mach shock-wave turbulent boundary layer interactions. AIAA Paper No [Cercignani, 988] Cercignani, C. (988). The Boltzmann equation and its applications. Springer, New York.
8 [Chen et al., 00] Chen, H., Kandasamy, S., Orszag, S., Shock, R., Succi, S., and Yakhot, V. (00). Extended Boltzmann kinetic equation for turbulent flows. Science, 0(56):6 66. [Chen et al., 00] Chen, H., Orszag, S., Staroselsky, I., and Succi, S. (00). Expanded analogy between Boltzmann kinetic theory of fluids and turbulence. J. Fluid Mech., 59():0. [Chou and Baganoff, 997] Chou, S. and Baganoff, D. (997). Kinetic flux vector splitting for the Navier Stokes equations. J. Comput. Phys., 0():7 0. [Cook et al., 979] Cook, P., McDonald, M., and Firman, M. (979). Aerofoil RAE 8 pressure distributions, and boundary layer andwake measurements. Experimental data base for computer program assessment. AGARD Advisory. [Délery et al., 986] Délery, J., Marvin, J., and Reshotko, E. (986). Shock-wave boundary layer interactions. AGARDograph, 80. [Dolling and Erengil, 99] Dolling, D. S. and Erengil, M. E. (99). Unsteady wave structure near separation in a Mach 5 compression rampinteraction. AIAA J, 9(5): [Edwards, 008] Edwards, J. R. (008). Numerical simulations of shock/boundary layer interactions using time-dependent modeling techniques: A survey of recent results. Prog. Aerosp. Sci., (6):7 65. [Goldberg et al., 998] Goldberg, U., Peroomian, O., and Chakravarthy, S. (998). Application of the k- e-r Turbulence Model to Wall-Bounded Compressive Flows. AIAA Paper No [Li et al., 00] Li, Q., Xu, K., and Fu, S. (00). A high-order gas-kinetic Navier Stokes flow solver. J. Comput. Phys., 9(9): [Liao et al., 007] Liao, W., Luo, L., and Xu, K. (007). Gas-kinetic scheme for continuum and near-continuum hypersonic flows. J. Spacecraft Rockets, (6): 0. [Mandal and Deshpande, 99] Mandal, J. and Deshpande, S. (99). Kinetic flux vector splitting for Euler equations. Comput. Fluids, ():7 78. [May et al., 007] May, G., Srinivasan, B., and Jameson, A. (007). An improved gas-kinetic BGK finite-volume method for three-dimensional transonic flow. J. Comput. Phys., 0(): [Menter and Rumsey, 99] Menter, F. and Rumsey, C. (99). Assessment of two-equation turbulence models for transonic flows. AIAA Paper No. 9-, (Paper No. 9-). [Righi, 0] Righi, M. (0). A finite-volume gaskinetic method for the solution of the Navier-Stokes equations. Royal Aeronautical Society, Aeronaut. J., (7)(9). [Righi, 0] Righi, M. (0). A Modified Gas- Kinetic Scheme for Turbulent Flow. Commun. Comput. Phys., 6():9 6. [Settles et al., 979] Settles, G., Fitzpatrick, T., and Bogdonoff, S. (979). Detailed study of attached and separated compression corner flowfields in high Reynolds number supersonic flow. AIAA J., 7(6): [Settles et al., 976] Settles, G., Vas, I., and Bogdonoff, S. (976). Details of a shock-separated turbulent boundary layer at a compression corner. AIAA J., (): [Spaid and Frishett, 97] Spaid, F. W. and Frishett, J. C. (97). Incipient separation of a supersonic, turbulent boundary layer, including effects of heat transfer. AIAA Journal, 0(7):95 9. [Wallin and Johansson, 000] Wallin, S. and Johansson, A. (000). An explicit algebraic Reynolds stress model for incompressible and compressible turbulent flows. J. Fluid Mech., 0:89. [Xu, 998] Xu, K. (998). Gas-kinetic schemes for unsteady compressible flow simulations. Von Karman Institute, Computational Fluid Dynamics, Annual Lecture Series, 9th, Rhode-Saint-Genese, Belgium. [Xu, 00] Xu, K. (00). A gas-kinetic BGK scheme for the Navier Stokes equations and its connection with artificial dissipation and Godunov method. J. Comput. Phys., 7():89 5. [Xu and Huang, 00] Xu, K. and Huang, J. (00). A unified gas-kinetic scheme for continuum and rarefied flows. J. Comput. Phys., 9(0): [Xu et al., 005] Xu, K., Mao, M., and Tang, L. (005). A multidimensional gas-kinetic BGK scheme for hypersonic viscous flow. J. Comput. Phys., 0():05. [Xu and Prendergast, 99] Xu, K. and Prendergast, K. (99). Numerical Navier-Stokes solutions from gas kinetic theory. J. Comput. Phys., ():9 7. [Xuan and Xu, 0] Xuan, L. and Xu, K. (0). A new gas-kinetic scheme based on analytical solutions of the BGK equation. J. Comput. Phys.
A Modified Gas-Kinetic Scheme for Turbulent Flow
Commun. Comput. Phys. doi: 10.4208/cicp.140813.021213a Vol. 16, No. 1, pp. 239-263 July 2014 A Modified Gas-Kinetic Scheme for Turbulent Flow Marcello Righi School of Engineering, Zurich University of
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationKinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods
Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
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 informationAn improved unified gas-kinetic scheme and the study of shock structures
IMA Journal of Applied Mathematics (2011) 76, 698 711 doi:10.1093/imamat/hxr002 Advance Access publication on March 16, 2011 An improved unified gas-kinetic scheme and the study of shock structures KUN
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationcompression corner flows with high deflection angle, for example, the method cannot predict the location
4nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 4, Reno, Nevada Modeling the effect of shock unsteadiness in shock-wave/ turbulent boundary layer interactions AIAA 4-9 Krishnendu Sinha*, Krishnan
More informationThe Equivalent Differential Equation of the Gas-Kinetic BGK Scheme
The Equivalent Differential Equation of the Gas-Kinetic BGK Scheme Wei Gao Quanhua Sun Kun Xu Abstract The equivalent partial differential equation is derived for the gas-kinetic BGK scheme (GKS) [Xu J.
More informationConnection Between the Lattice Boltzmann Equation and the Beam Scheme
NASA/CR-1999-209001 ICASE Report No. 99-10 Connection Between the Lattice Boltzmann Equation and the Beam Scheme Kun Xu Hong Kong University of Science and Technology, Kowloon, Hong Kong Li-Shi Luo ICASE,
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
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 informationA Gas-Kinetic BGK Scheme for the Navier Stokes Equations and Its Connection with Artificial Dissipation and Godunov Method 1
Journal of Computational Physics 171, 289 335 (2001) doi:10.1006/jcph.2001.6790, available online at http://www.idealibrary.com on A Gas-Kinetic BGK Scheme for the Navier Stokes Equations and Its Connection
More information2 Three Dimensional Finite Volume Gas-Kinetic Schemes
Gas-Kinetic Finite Volume Methods K. Xu, L. Martinelli, A. Jameson Department of M.A.E, Princeton University, Prificeton N.J. 08544, USA 1 Introduction Gas-kinetic schemes developed from the BGK model
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows Kun Xu a,, Juan-Chen Huang b a Mathematics Department Hong Kong University of Science and Technology Clear Water Bay, Kowloon, Hong Kong b
More informationIntroduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)
Introduction to Aerodynamics Dr. Guven Aerospace Engineer (P.hD) Aerodynamic Forces All aerodynamic forces are generated wither through pressure distribution or a shear stress distribution on a body. The
More informationStability and consistency of kinetic upwinding for advection diffusion equations
IMA Journal of Numerical Analysis (006 6, 686 7 doi:0.093/imanum/drl005 Advance Access publication on March, 006 Stability and consistency of kinetic upwinding for advection diffusion equations MANUEL
More informationA gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions
Acta Mech. Sin. 2017) 334):733 741 DOI 10.1007/s10409-017-0695-2 RESEARCH PAPER A gas-kinetic theory based multidimensional high-order method for the compressible Navier Stokes solutions Xiaodong Ren 1
More informationRegularization of the Chapman-Enskog Expansion and Its Description of Shock Structure
NASA/CR-2001-211268 ICASE Report No. 2001-39 Regularization of the Chapman-Enskog Expansion and Its Description of Shock Structure Kun Xu Hong Kong University, Kowloon, Hong Kong ICASE NASA Langley Research
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationAssessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow
Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering - University
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More information1. Introduction Some Basic Concepts
1. Introduction Some Basic Concepts 1.What is a fluid? A substance that will go on deforming in the presence of a deforming force, however small 2. What Properties Do Fluids Have? Density ( ) Pressure
More informationDirect modeling for computational fluid dynamics
Acta Mech Sin (5) 3(3):33 38 DOI.7/s49-5-453- RESEARCH PAPER Direct modeling for computational fluid dynamics Kun u Received: 5 January 5 / Revised: 5 February 5 / Accepted: 3 March 5 / Published online:
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationNew DNS Results of Shockwave/Turbulent Boundary Layer Interaction
36th AIAA Fluid Dynamics Conference, June 5 8, 26/San Francisco,California New DNS Results of Shockwave/Turbulent Boundary Layer Interaction M. Wu and M.P. Martin Mechanical and Aerospace Engineering Department
More informationA hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models
A hybrid method for hydrodynamic-kinetic flow - Part II - Coupling of hydrodynamic and kinetic models Alessandro Alaia, Gabriella Puppo May 31, 2011 Abstract In this work we present a non stationary domain
More informationWall treatments and wall functions
Wall treatments and wall functions A wall treatment is the set of near-wall modelling assumptions for each turbulence model. Three types of wall treatment are provided in FLUENT, although all three might
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 informationA combined application of the integral wall model and the rough wall rescaling-recycling method
AIAA 25-299 A combined application of the integral wall model and the rough wall rescaling-recycling method X.I.A. Yang J. Sadique R. Mittal C. Meneveau Johns Hopkins University, Baltimore, MD, 228, USA
More informationTAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF
TAU Extensions for High Enthalpy Flows Sebastian Karl AS-RF Contents Motivation Extensions available in the current release: Numerical schemes for super- and hypersonic flow fields Models for gas mixtures,
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationContinuity Equation for Compressible Flow
Continuity Equation for Compressible Flow Velocity potential irrotational steady compressible Momentum (Euler) Equation for Compressible Flow Euler's equation isentropic velocity potential equation for
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationEquivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme
Equivalence between kinetic method for fluid-dynamic equation and macroscopic finite-difference scheme Pietro Asinari (1), Taku Ohwada (2) (1) Department of Energetics, Politecnico di Torino, Torino 10129,
More informationComputers and Mathematics with Applications
Computers and Mathematics with Applications 59 (2010 2194 2199 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Toward
More informationProjection Dynamics in Godunov-Type Schemes
JOURNAL OF COMPUTATIONAL PHYSICS 142, 412 427 (1998) ARTICLE NO. CP985923 Projection Dynamics in Godunov-Type Schemes Kun Xu and Jishan Hu Department of Mathematics, Hong Kong University of Science and
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationNUMERICAL SIMULATION OF HIGH-SPEED SEPARATION FLOW IN THE AEROSPACE PROPULSION SYSTEMS
NUMERICAL SIMULATION OF HIGH-SPEED SEPARATION FLOW IN THE AEROSPACE PROPULSION SYSTEMS E.V. Larina*, I.A. Kryukov**,I.E. Ivanov * * Moscow Aviation Institute (National Research University), Russia ** Institute
More informationScaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09
Scaling Parameters in Rarefied Flow and the Breakdown of the Navier-Stokes Equations Mechanical Engineering Research Report No: 2004/09 Michael Macrossan, Centre for Hypersonics, University of Queensland
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 informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationAER1310: TURBULENCE MODELLING 1. Introduction to Turbulent Flows C. P. T. Groth c Oxford Dictionary: disturbance, commotion, varying irregularly
1. Introduction to Turbulent Flows Coverage of this section: Definition of Turbulence Features of Turbulent Flows Numerical Modelling Challenges History of Turbulence Modelling 1 1.1 Definition of Turbulence
More informationMicrochannel flow in the slip regime: gas-kinetic BGK Burnett solutions
J. Fluid Mech. (2004), vol. 513, pp. 87 110. c 2004 Cambridge University Press DOI: 10.1017/S0022112004009826 Printed in the United Kingdom 87 Microchannel flow in the slip regime: gas-kinetic BGK Burnett
More informationGas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations
Gas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations Iizln Xu* and Antony ~arneson~ Department of Mechanical and Aerospace Engineering Princeton University Princeton, New Jersey
More informationTurbulent eddies in the RANS/LES transition region
Turbulent eddies in the RANS/LES transition region Ugo Piomelli Senthil Radhakrishnan Giuseppe De Prisco University of Maryland College Park, MD, USA Research sponsored by the ONR and AFOSR Outline Motivation
More informationNONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS
June - July, 5 Melbourne, Australia 9 7B- NONLINEAR FEATURES IN EXPLICIT ALGEBRAIC MODELS FOR TURBULENT FLOWS WITH ACTIVE SCALARS Werner M.J. Lazeroms () Linné FLOW Centre, Department of Mechanics SE-44
More informationApplication of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows
Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows Timothy R. Deschenes and Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor,
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
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 informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationThe Simulation of Wraparound Fins Aerodynamic Characteristics
The Simulation of Wraparound Fins Aerodynamic Characteristics Institute of Launch Dynamics Nanjing University of Science and Technology Nanjing Xiaolingwei 00 P. R. China laithabbass@yahoo.com Abstract:
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationGiven the water behaves as shown above, which direction will the cylinder rotate?
water stream fixed but free to rotate Given the water behaves as shown above, which direction will the cylinder rotate? ) Clockwise 2) Counter-clockwise 3) Not enough information F y U 0 U F x V=0 V=0
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationThere are no simple turbulent flows
Turbulence 1 There are no simple turbulent flows Turbulent boundary layer: Instantaneous velocity field (snapshot) Ref: Prof. M. Gad-el-Hak, University of Notre Dame Prediction of turbulent flows standard
More informationColloquium FLUID DYNAMICS 2012 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 2012 p.
Colloquium FLUID DYNAMICS 212 Institute of Thermomechanics AS CR, v.v.i., Prague, October 24-26, 212 p. ON A COMPARISON OF NUMERICAL SIMULATIONS OF ATMOSPHERIC FLOW OVER COMPLEX TERRAIN T. Bodnár, L. Beneš
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 informationAnalysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data
45th AIAA Aerospace Science Meeting and Exhibit, January 8 11, 25/Reno, Nevada Analysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data M. Wu and M.P. Martin Mechanical
More informationLattice Boltzmann Methods for Fluid Dynamics
Lattice Boltzmann Methods for Fluid Dynamics Steven Orszag Department of Mathematics Yale University In collaboration with Hudong Chen, Isaac Goldhirsch, and Rick Shock Transient flow around a car LBGK
More informationA New Class of Gas-Kinetic Relaxation Schemes for the Compressible Euler Equations
Journal of Statistical Physics, Vol. 81, Nos. 1/2, 1995 A New Class of Gas-Kinetic Relaxation Schemes for the Compressible Euler Equations Kun Xu I Received September 1. 1994; final December 23. 1994 Starting
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 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 informationMasters in Mechanical Engineering Aerodynamics 1 st Semester 2015/16
Masters in Mechanical Engineering Aerodynamics st Semester 05/6 Exam st season, 8 January 06 Name : Time : 8:30 Number: Duration : 3 hours st Part : No textbooks/notes allowed nd Part : Textbooks allowed
More informationSIMULATION OF GAS FLOW OVER MICRO-SCALE AIRFOILS USING A HYBRID CONTINUUM-PARTICLE APPROACH
33rd AIAA Fluid Dynamics Conference and Exhibit 3-6 June 3, Orlando, Florida AIAA 3-44 33 rd AIAA Fluid Dynamics Conference and Exhibit / Orlando, Florida / 3-6 Jun 3 SIMULATION OF GAS FLOW OVER MICRO-SCALE
More informationBulk equations and Knudsen layers for the regularized 13 moment equations
Continuum Mech. Thermon. 7 9: 77 89 DOI.7/s6-7-5- ORIGINAL ARTICLE Henning Struchtrup Toby Thatcher Bulk equations and Knudsen layers for the regularized 3 moment equations Received: 9 December 6 / Accepted:
More informationHigh Speed Aerodynamics. Copyright 2009 Narayanan Komerath
Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil
More informationGeneralized Local Equilibrium in the Cascaded Lattice Boltzmann Method. Abstract
Accepted for publication on Physical Review E (R), code ERR1034 Generalized Local Equilibrium in the Cascaded Lattice Boltzmann Method Pietro Asinari Department of Energetics, Politecnico di Torino, Corso
More informationCOMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL
COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the
More informationPredicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions
Predicting Breakdown of the Continuum Equations Under Rarefied Flow Conditions Iain D. Boyd Department of Aerospace Engineering, University of Michigan, Ann Arbor, MI 48109 Abstract. The breakdown of the
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationLES Study of Shock Wave and Turbulent Boundary Layer Interaction
51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 07-10 January 2013, Grapevine (Dallas/Ft. Worth Region), Texas AIAA 2013-0984 LES Study of Shock Wave and
More informationNumerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form
Numerical Simulation of Rarefied Gases using Hyperbolic Moment Equations in Partially-Conservative Form Julian Koellermeier, Manuel Torrilhon May 18th, 2017 FU Berlin J. Koellermeier 1 / 52 Partially-Conservative
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 informationThe lattice Boltzmann equation (LBE) has become an alternative method for solving various fluid dynamic
36th AIAA Fluid Dynamics Conference and Exhibit 5-8 June 2006, San Francisco, California AIAA 2006-3904 Direct and Large-Eddy Simulation of Decaying and Forced Isotropic Turbulence Using Lattice Boltzmann
More informationZonal modelling approach in aerodynamic simulation
Zonal modelling approach in aerodynamic simulation and Carlos Castro Barcelona Supercomputing Center Technical University of Madrid Outline 1 2 State of the art Proposed strategy 3 Consistency Stability
More informationEffective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body
Effective Boundary Conditions for Continuum Method of Investigation of Rarefied Gas Flow over Blunt Body I.G. Brykina a, B.V. Rogov b, I.L. Semenov c, and G.A. Tirskiy a a Institute of Mechanics, Moscow
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
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 informationJournal of Computational Physics
Journal of Computational Physics 229 (2010 7747 7764 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A unified gas-kinetic scheme
More informationREGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS
1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,
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 informationSEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS
The Seventh Nobeyama Workshop on CFD: To the Memory of Prof. Kuwahara Sanjo Kaikan, the University of Tokyo, Japan, September. 23-24, 2009 SEMICLASSICAL LATTICE BOLTZMANN EQUATION HYDRODYNAMICS Jaw-Yen
More informationComputation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations
5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 9-2 January 22, Nashville, Tennessee AIAA 22-699 Computation of NACA2 Airfoil Transonic Buffet Phenomenon with
More informationMass Transfer in Turbulent Flow
Mass Transfer in Turbulent Flow ChEn 6603 References: S.. Pope. Turbulent Flows. Cambridge University Press, New York, 2000. D. C. Wilcox. Turbulence Modeling for CFD. DCW Industries, La Caada CA, 2000.
More informationA Hybrid WENO Scheme for Simulation of Shock Wave-Boundary Layer Interaction
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-36 A Hybrid WENO Scheme for Simulation of Shock Wave-Boundary Layer
More informationFluid Mechanics Theory I
Fluid Mechanics Theory I Last Class: 1. Introduction 2. MicroTAS or Lab on a Chip 3. Microfluidics Length Scale 4. Fundamentals 5. Different Aspects of Microfluidcs Today s Contents: 1. Introduction to
More informationShock/boundary layer interactions
Shock/boundary layer interactions Turbulent compressible channel flows F.S. Godeferd Laboratoire de Mécanique des Fluides et d Acoustique Ecole Centrale de Lyon, France Journée Calcul Intensif en Rhône
More informationHypersonic Blunt Body Thermophysics Using a Unified Kinetic/Continuum Solver
41st AIAA Thermophysics Conference 22-25 June 2009, San Antonio, Texas AIAA 2009-3838 Hypersonic Blunt Body Thermophysics Using a Unified Kinetic/Continuum Solver Andrew J. Lofthouse U.S. Air Force Institute
More informationStudy on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow
Research Article Study on lattice Boltzmann method/ large eddy simulation and its application at high Reynolds number flow Advances in Mechanical Engineering 1 8 Ó The Author(s) 2015 DOI: 10.1177/1687814015573829
More informationA Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations
An Extended Abstract submitted for the 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada January 26 Preferred Session Topic: Uncertainty quantification and stochastic methods for CFD A Non-Intrusive
More informationDefense Technical Information Center Compilation Part Notice
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADP013688 TITLE: Large-Eddy Simulation of Mach 3.0 Flow Past a 24-Degree Compression Ramp DISTRIBUTION: Approved for public release,
More informationThe mean shear stress has both viscous and turbulent parts. In simple shear (i.e. U / y the only non-zero mean gradient):
8. TURBULENCE MODELLING 1 SPRING 2019 8.1 Eddy-viscosity models 8.2 Advanced turbulence models 8.3 Wall boundary conditions Summary References Appendix: Derivation of the turbulent kinetic energy equation
More informationSimulation of 2D non-isothermal flows in slits using lattice Boltzmann method
Simulation of 2D non-isothermal flows in slits using lattice Boltzmann method João Chambel Leitão Department of Mechanical Engineering, Instituto Superior Técnico, Lisboa, Portugal Abstract: The present
More information