Stabilization and shock-capturing parameters in SUPG formulation of compressible flows
|
|
- Tabitha Potter
- 5 years ago
- Views:
Transcription
1 Comput. Methods Appl. Mech. Engrg. 195 (006) Stabilization and shock-capturing parameters in SUPG formulation of compressible flows Tayfun E. Tezduyar *, Masayoshi Senga Mechanical Engineering, Rice University MS 31, 6100 Main Street, Houston, TX , USA Received 7 October 004; received in revised form 10 January 005; accepted 10 May 005 Abstract The streamline-upwind/petrov Galerkin (SUPG) formulation is one of the most widely used stabilized methods in finite element computation of compressible flows. It includes a stabilization parameter that is known as s. Typically the SUPG formulation is used in combination with a shock-capturing term that provides additional stability near the shock fronts. The definition of the shock-capturing term includes a shock-capturing parameter. In this paper, we describe, for the finite element formulation of compressible flows based on conservation variables, new ways for determining the s and the shock-capturing parameter. The new definitions for the shock-capturing parameter are much simpler than the one based on the entropy variables, involve less operations in calculating the shock-capturing term, and yield better shock quality in the test computations. Ó 005 Elsevier B.V. All rights reserved. Keywords: Compressible flows; Finite element formulation; SUPG stabilization, Stabilization parameters; Shock-capturing parameter 1. Introduction In finite element computation of flow problems, the streamline-upwind/petrov Galerkin (SUPG) formulation for incompressible flows [1,], the SUPG formulation for compressible flows [3 5], and the pressurestabilizing/petrov Galerkin (PSPG) formulation for incompressible flows [6] are some of the most prevalent stabilized methods. Stabilized formulations such as the SUPG and PSPG formulations prevent numerical instabilities in solving problems with high Reynolds or Mach numbers and shocks or thin boundary layers, as well as when using equal-order interpolation functions for velocity and pressure. * Corresponding author. Tel.: ; fax: addresses: tezduyar@rice.edu (T.E. Tezduyar), msenga@rice.edu (M. Senga) /$ - see front matter Ó 005 Elsevier B.V. All rights reserved. doi: /j.cma
2 16 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) The SUPG formulation for incompressible flows was first introduced in [1], with further studies in []. The SUPG formulation for compressible flows was first introduced, in the context of conservation variables in [3]. A concise version of that was published as an AIAA paper [4], and a more thorough version with additional examples as a journal paper [5]. After that, several SUPG-like methods for compressible flows were developed. Taylor Galerkin method [7], for example, is very similar, and under certain conditions is identical, to one of the SUPG methods introduced in [3 5]. Later, following [3 5], the SUPG formulation for compressible flows was recast in entropy variables and supplemented with a shock-capturing term [8]. It was shown in [9] that the SUPG formulation introduced in [3 5], when supplemented with a similar shock-capturing term, is very comparable in accuracy to the one that was recast in entropy variables. The stabilized formulation introduced in [10] for advection diffusion reaction equations also included a shock-capturing (discontinuity-capturing) term, and precluded augmentation of the SUPG effect by the discontinuity-capturing effect when the advection and discontinuity directions coincide. A stabilization parameter, known as s, is embedded in the SUPG and PSPG formulations. It involves a measure of the local length scale (also known as element length ) and other parameters such as the element Reynolds and Courant numbers. Various element lengths and ss were proposed starting with those in [1 5], followed by the one introduced in [10], and those proposed in the subsequently reported SUPG-based methods. Here we will call the SUPG formulation introduced in [3 5] for compressible flows (SUPG) 8, and the set of ss introduced in conjunction with that s 8. The s used in [9] with (SUPG) 8 is a slightly modified version of s 8. A shock-capturing parameter, which we will call here d 91, was embedded in the shock-capturing term used in [9]. Subsequent minor modifications of s 8 took into account the interaction between the shock-capturing and the (SUPG) 8 terms in a fashion similar to how it was done in [10] for advection diffusion reaction equations. All these slightly modified versions of s 8 have always been used with the same d 91, and we will categorize them here all under the label s 8-MOD. To be used in conjunction with the SUPG/PSPG formulation of incompressible flows, the discontinuitycapturing directional dissipation (DCDD) stabilization was introduced in [11,1] for flow fields with sharp gradients. This involved a second element length scale, which was also introduced in [11,1] and is based on the solution gradient. This new element length scale is used together with the element length scales already defined in [10]. Recognizing this second element length as a diffusion length scale, new stabilization parameters for the diffusive limit were introduced in [1 14]. The DCDD stabilization was originally conceived in [11,1] as an alternative to the LSIC (least-squares on incompressibility constraint) stabilization. The DCDD takes effect where there is a sharp gradient in the velocity field and introduces dissipation in the direction of that gradient. The way the DCDD is added to the formulation precludes augmentation of the SUPG effect by the DCDD effect when the advection and discontinuity directions coincide. Partly based on the ideas underlying the new ss for incompressible flows and the DCDD, new ways of calculating the ss and shock-capturing parameters for compressible flows were introduced in [14 17]. Like the parameters developed earlier, these new parameters are intended for use with the SUPG formulation of compressible flows based on conservation variables. In this paper, we describe how the new parameters are defined.. Navier Stokes equations of compressible flows Let X R n sd be the spatial domain with boundary C, and (0, T) be the time domain. The symbols q, u, p and e will represent the density, velocity, pressure and the total energy, respectively. The Navier Stokes equations of compressible flows can be written on X and "t (0, T) as ou ot þ of i ox i oe i ox i R ¼ 0; ð1þ
3 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) where U =(q,qu 1,qu,qu 3,qe) is the vector of conservation variables, and F i and E i are, respectively, the Euler and viscous flux vectors: u i q 0 u i qu 1 þ d i1 p T i1 F i ¼ u i qu þ d i p ; E i ¼ T i. ðþ B C B u i qu 3 þ d i3 p T i3 A u i ðqe þ pþ q i þ T ik u k Here d ij are the components of the identity tensor I, q i are the components of the heat flux vector, and T ij are the components of the Newtonian viscous stress tensor: T ¼ kð$ uþi þ leðuþ; ð3þ where k and l (=qm) are the viscosity coefficients, m is the kinematic viscosity, and e(u) is the strain-rate tensor: eðuþ ¼ 1 ðð$uþþð$uþt Þ. ð4þ It is assumed that k = l/3. The equation of state used here corresponds to the ideal gas assumption. The term R represents all other components that might enter the equations, including the external forces. Eq. (1) can further be written in the following form: ou ot þ A ou i o ox i ox i where K ij ou ox j R ¼ 0; A i ¼ of i ou ; K ou ij ¼ E i. ox j Appropriate sets of boundary and initial conditions are assumed to accompany Eq. (5). ð5þ ð6þ 3. SUPG formulations 3.1. Semi-discrete Given Eq. (5), we form some suitably-defined finite-dimensional trial solution and test function spaces S h U and Vh U. The SUPG formulation of Eq. (5) can then be written as follows: find Uh S h U such that 8W h V h U : W h ou h þ A h ou h ow h i dx þ K h ou h ij dx W h H h dc X ot ox i X ox i ox j C H W h R h dx þ Xn el h ow s SUPG A h ou h k þ A h ou h i o K h ou h ij R h dx X e¼1 X e ox k ot ox i ox i ox j þ Xn el h ow ou h m SHOC dx ¼ 0. ð7þ e¼1 X e ox i ox i Here H h represents the natural boundary conditions associated with Eq. (5), and C H is the part of the boundary where such boundary conditions are specified. The SUPG stabilization and shock-capturing
4 164 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) parameters are denoted by s SUPG and m SHOC. They were discussed in Section 1 and will further be discussed in Section Space time The space time version of Eq. (7) can be written based on the deforming-spatial-domain/stabilized space time (DSD/SST) formulation introduced in [6,18,19]. The finite element formulation of the governing equations is written over a sequence of N space time slabs Q n, where Q n is the slice of the space time domain between the time levels t n and t n+1. At each time step, the integrations involved in the finite element formulation are performed over Q n. The finite element interpolation functions are discontinuous across the space time slabs. We use the notation ðþ n and ðþþ n to denote the values as t n is approached from below and above respectively. Each Q n is decomposed into space time elements Q e n, where e =1,,...,(n el) n. The subscript n used with n el is to account for the general case in which the number of space time elements may change from one space time slab to another. For each slab Q n, we form some suitably-defined finite-dimensional trial solution and test function spaces ðs h U Þ n and ðvh U Þ n. In the computations reported here, we use first-order polynomials as interpolation functions. The subscript n implies that corresponding to different space time slabs we might have different discretizations. The DSD/SST formulation of Eq. (5) can then be written as follows: given ðu h Þ n, find Uh ðs h U Þ n such that 8Wh ðv h U Þ n : W h ou h þ A h ou h h ow i dq þ K h ou h ij Q n ot ox i Q n ox i ox j W h R h dq þ ðw h Þ þ n ðuh Þ þ n ðuh Þ n dx Q n X þ Xðn elþ n ow h s SUPG A h ou h k þ A h ou h i o e¼1 Q e ox n k ot ox i ox i þ Xðn elþ n ow h ou h m SHOC dq ¼ 0. ox i ox i e¼1 Q e n dq K h ij ou h ox j P n W h H h dp R h Here P n is the lateral boundary of the space time slab. The solution to Eq. (8) is obtained sequentially for all space time slabs Q 0,Q 1,Q,..., Q N 1, and the computations start with ðu h Þ 0 ¼ Uh 0, where U 0 is the specified initial value of the vector U. dq ð8þ 4. Calculation of the stabilization parameters for compressible flows and shock-capturing Various options for calculating the stabilization parameters and defining the shock-capturing terms in the context of the (SUPG) 8 formulation were introduced in [14 17]. In this section we describe those options. For this purpose, we first define the acoustic speed as c, and define the unit vector j as j ¼ $qh k$q h k. ð9þ As the first option in computing s SUGN1 for each component of the test vector-function W, the stabilization parameters s q SUGN1, su SUGN1 and se SUGN1 (associated with q, qu and qe, respectively) are defined by the following expression: s q SUGN1 ¼ su SUGN1 ¼ se SUGN1 ¼ Xn en a¼1 ju h $N a j! 1. ð10þ
5 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) As the second option, they are defined as s q SUGN1 ¼ su SUGN1 ¼ se SUGN1 ¼ Xnen a¼1 ðcjj $N a jþju h $N a jþ! 1. ð11þ In computing s SUGN, the parameters s q SUGN, su SUGN and se SUGN are defined as follows: s q SUGN ¼ su SUGN ¼ se SUGN ¼ Dt ; ð1þ where Dt is the time step. In computing s SUGN3, the parameter s u SUGN3 is defined by using the expression s u SUGN3 ¼ h RGN 4m ; ð13þ where h RGN ¼ Xnen a¼1! 1 jr $N a j ; r ¼ $kuh k k$ku h kk. The parameter s e SUGN3 is defined as ð14þ s e SUGN3 ¼ ðhe RGN Þ 4m e ; ð15þ where m e is the kinematic viscosity for the energy equation, h e RGN ¼ Xnen a¼1 jr e $N a j! 1 ; r e ¼ $hh k$h h k ð16þ and h is the temperature. The parameters ðs q SUPG Þ UGN, ðsu SUPG Þ UGN and ðse SUPG Þ UGN are calculated from their components by using the r-switch : 1 ðs q SUPG Þ UGN ¼ 1 1 r; ðs q þ ð17þ SUGN1 Þr ðs q SUGN Þr 1 ðs u SUPG Þ UGN ¼ r; þ þ ð18þ ðs u SUGN1 Þr ðs u SUGN Þr ðs u SUGN3 Þr 1 ðs e SUPG Þ UGN ¼ r. þ þ ð19þ ðs e SUGN1 Þr ðs e SUGN Þr ðs e SUGN3 Þr This r-switch was first introduced in [0]. Typically, r =. As the first option in defining the shock-capturing term, first the shock-capturing viscosity m SHOC is defined: m SHOC ¼ s SHOC ðu int Þ ; where ð0þ s SHOC ¼ h SHOC u cha h SHOC ¼ h JGN ; k$q h b kh SHOC ; ð1þ q ref ðþ
6 166 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) h JGN ¼ Xnen a¼1 jj $N a j! 1. ð3þ Here q ref is a reference density (such as q h at the inflow, or the difference between the estimated maximum and minimum values of q h ), u cha is a characteristic velocity (such as u ref or ku h k or acoustic speed c), and u int is an intrinsic velocity (such as u cha or ku h kor acoustic speed c). Typically, u int = u cha = u ref. The parameter b influences the smoothness of the shock-front. It is set as b = 1 for smoother shocks and b = for sharper shocks (in return for tolerating possible overshoots and undershoots). The compromise between the b = 1 and selections is defined as the following averaged expression for s SHOC : s SHOC ¼ 1 ððs SHOCÞ b¼1 þðs SHOC Þ b¼ Þ. ð4þ As another option for calculating the shock-capturing parameter, m SHOC is defined as m SHOC ¼kY 1 k Xn sd i¼1 Y 1 ou h ox i! b= 1 b h SHOC ; ð5þ where Y is a diagonal scaling matrix constructed from the reference values of the components of U: 3 ðu 1 Þ ref ðu Þ ref Y ¼ 0 0 ðu 3 Þ ref 0 0 ; ð6þ ðu 4 Þ ref ðu 5 Þ ref ¼ ouh ot þ A h i ou h ox i ð7þ or ¼ A h ou h i ð8þ ox i and b =1orb =. In a variation of the expression given by Eq. (5), m SHOC is defined by the following expression: m SHOC ¼kY 1 k Xn sd i¼1 Y 1 ou h ox i! b= 1 b ky 1 U h k 1 b h SHOC. ð9þ The compromise between the b = 1 and selections is defined as the following averaged expression for m SHOC : m SHOC ¼ 1 ððm SHOCÞ b¼1 þðm SHOC Þ b¼ Þ. ð30þ Based on Eq. (5), a separate m SHOC can be calculated for each component of the test vector-function W: ðm SHOC Þ I ¼jðY 1 Þ I j X n sd i¼1 Y 1 ou h ox i I! b= 1 b h SHOC ; I ¼ 1; ;...; n sd þ. ð31þ
7 Similarly, a separate m SHOC for each component of W can be calculated based on Eq. (9): X n ðm SHOC Þ I ¼jðY 1 sd Þ I j Y 1 ou h! b= 1 ox i¼1 i I Given m SHOC, the shock-capturing term is defined as S SHOC ¼ Xn el $W h : ðj SHOC $U h ÞdX; X e e¼1 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) jðy 1 U h Þ I j 1 b b h SHOC ; I ¼ 1; ;...; n sd þ. ð3þ where j SHOC is defined as j SHOC = m SHOC I. As a possible alternative, it is defined as j SHOC = m SHOC jj. If the option given by Eq. (31) or Eq. (3) is exercised, then m SHOC becomes an (n sd +) (n sd + ) diagonal matrix, and the matrix j SHOC becomes augmented from an n sd n sd matrix to an (n sd (n sd + )) ((n sd +) n sd ) matrix. To preclude compounding, m SHOC can be modified as follows: m SHOC m SHOC switchðs SUPG ðj uþ ; s SUPG ðjj uj cþ ; m SHOC Þ; ð34þ where the switch function is defined as the min function or as the r-switch used earlier in this section. For viscous flows, the above modification would be made separately with each of s q SUPG, su SUPG and s e SUPG, and this would result in m SHOC becoming a diagonal matrix even if the option given by Eq. (31) or Eq. (3) is not exercised. ð33þ 5. Test computations The test computations were carried out by using the space time SUPG formulation described in Section 3.. We used two steady-state, inviscid test problems: oblique shock and reflected shock. These were y line plots M =.0 Shock 10 M = x Fig. 1. Oblique shock. Problem description.
8 168 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) used in many earlier publications, and here we compute each of them with two different options for the shock-capturing parameter. In the option denoted by CY1, we use Eq. (5) with Eq. (30), and in the option denoted by CYU1, Eq. (9) with Eq. (30). In both options, we use for the expression given by Eq. (8), and set j SHOC = m SHOC I. With both options, as stabilization parameters, we use Eq. (11), and in Eqs. (17) (19) we do not include s SUGN. In both problems, the time-step size is The num- 1.6 CY1 1.5 τ 8-MOD with δ 91 Exact solution 1.4 Density y CYU1 τ 8-MOD with δ 91 Exact solution 1.4 Density y Fig.. Oblique shock. Density along x = 0.9, obtained with CY1 (top) and CYU1 (bottom), compared with the solution obtained with the s 8-MOD and d 91 combination.
9 ber of time steps, nonlinear iterations, and the inner and outer GMRES iterations are 100, 3, 30, and, respectively. The results are compared to those obtained with the s 8-MOD and d 91 combination. The version of s 8-MOD used in this paper for comparison is similar to the one given in [1]: s 8-MOD ¼ maxð0; fðs a s d ÞÞ; ð35þ where s a ¼ h BGN ; u cc s d ¼ d 91 ðu cc Þ ; u cc ¼ c þ u h $ku h k k$ku h kk ; ð36þ f ¼ T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) u ccdt=h BGN 1 þ u cc Dt=h BGN ; h BGN ¼ Xnen a¼1 1 $ku h k k$ku h kk $N a!. ð37þ Oblique shock. Fig. 1 shows the problem description. This is a Mach uniform flow over a wedge at an angle of 10 with the horizontal wall. The solution involves an oblique shock at an angle of 9.3 emanating from the leading edge. The computational domain is a square with 0 6 x 6 1 and 0 6 y 6 1. The inflow conditions are given as M =.0, q = 1.0, u 1 = cos10, u = sin10, andp = This results in an exact solution with the following outflow data: M = 1.64, q = 1.46, u 1 = 0.887, u = 0.0, and p = All essential boundary conditions are imposed at the left and top boundaries, slip condition at the wall, and no boundary conditions at the right boundary. The mesh is uniform and consists of 0 0 elements. Fig. shows the density along x = 0.9, obtained with CY1 and CYU1, compared with the solution obtained with the s 8-MOD and d 91 combination. In addition to being much simpler, the new shock-capturing parameters yield shocks with better quality. Reflected shock. Fig. 3 shows the problem description. This problem consists of three flow regions (R1, R and R3) separated by an oblique shock and its reflection from the wall. The computational domain is a rectangle with 0 6 x and 0 6 y 6 1. The inflow conditions in R1 are given as M =.9, q = 1.0, u 1 =.9, u = 0.0, and p = Specifying these conditions and requiring the incident shock to be at an angle of 9 results in an exact solution with the following flow data: R: M =.378, q = 1.7, u 1 =.619, u = 0.506, and p = 1.58; R3: M = 1.94, q =.687, u 1 =.401, u = 0.0, and p =.934. All essential boundary conditions are imposed at the left and top boundaries, slip condition at the wall, and no boundary conditions at the right boundary. The mesh is uniform and consists of 60 0 elements. Fig. 4 shows the density along y = 0.5, obtained with CY1 and CYU1, compared with the solution obtained with the s 8-MOD and d 91 combination. Again, the new, much simpler shock-capturing parameters yield shocks with better quality. y M =.378 M =.9 line plots 1 9 Shock 3 3 M = 1.94 x Fig. 3. Reflected shock. Problem description.
10 1630 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) CY1 τ 8-MOD with δ 91 Exact solution Density x 3.5 CYU1 τ 8-MOD with δ 91 Exact solution Density x Fig. 4. Reflected shock. Density along y = 0.5, obtained with CY1 (top) and CYU1 (bottom), compared with the solution obtained with the s 8-MOD and d 91 combination. 6. Concluding remarks We described, for the streamline-upwind/petrov Galerkin (SUPG) formulation of compressible flows based on conservation variables, new ways for determining the stabilization and shock-capturing parameters. The stabilization parameter, which is typically known as s, plays an important role in determining the accuracy of the solutions. The shock-capturing term provides additional stabilization near the shocks,
11 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) and how the shock-capturing parameter it involves is defined influences the quality of the solution near the shocks. These new ways of calculating the ss and shock-capturing parameters are partly based on the ideas underlying the ss and and DCDD stabilization developed for incompressible flows. Compared to the earlier shock-capturing parameter that was derived based on the entropy variables, the new ones are much simpler, involve less operations in calculating the shock-capturing term, and gave better shock resolution in the test computations we carried out. Acknowledgements This work was supported by the US Army Natick Soldier Center (Contract No. DAAD16-03-C-0051), NSF (Grant No. EIA ), and NASA Johnson Space Center (Grant No. NAG9-1435). References [1] T.J.R. Hughes, A.N. Brooks, A multi-dimensional upwind scheme with no crosswind diffusion, in: T.J.R. Hughes (Ed.), Finite Element Methods for Convection Dominated Flows, AMD-vol. 34, ASME, New York, 1979, pp [] A.N. Brooks, T.J.R. Hughes, Streamline upwind/petrov Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier Stokes equations, Comput. Methods Appl. Mech. Engrg. 3 (198) [3] T.E. Tezduyar, T.J.R. Hughes, Development of time-accurate finite element techniques for first-order hyperbolic systems with particular emphasis on the compressible Euler equations, NASA Technical Report NASA-CR-0477, NASA, 198, Available from: < [4] T.E. Tezduyar, T.J.R. Hughes, Finite element formulations for convection dominated flows with particular emphasis on the compressible Euler equations, in: Proceedings of AIAA 1st Aerospace Sciences Meeting, AIAA Paper , Reno, Nevada, [5] T.J.R. Hughes, T.E. Tezduyar, Finite element methods for first-order hyperbolic systems with particular emphasis on the compressible Euler equations, Comput. Methods Appl. Mech. Engrg. 45 (1984) [6] T.E. Tezduyar, Stabilized finite element formulations for incompressible flow computations, Adv. Appl. Mech. 8 (199) [7] J. Donea, A Taylor Galerkin method for convective transport problems, Int. J. Numer. Methods Engrg. 0 (1984) [8] T.J.R. Hughes, L.P. Franca, M. Mallet, A new finite element formulation for computational fluid dynamics: VI. Convergence analysis of the generalized SUPG formulation for linear time-dependent multi-dimensional advective diffusive systems, Comput. Methods Appl. Mech. Engrg. 63 (1987) [9] G.J. Le Beau, T.E. Tezduyar, Finite element computation of compressible flows with the SUPG formulation, in: Advances in Finite Element Analysis in Fluid Dynamics, FED-vol. 13, ASME, New York, 1991, pp [10] T.E. Tezduyar, Y.J. Park, Discontinuity capturing finite element formulations for nonlinear convection diffusion reaction problems, Comput. Methods Appl. Mech. Engrg. 59 (1986) [11] T.E. Tezduyar, Adaptive determination of the finite element stabilization parameters, in: Proceedings of the ECCOMAS Computational Fluid Dynamics Conference 001 (CD-ROM), Swansea, Wales, United Kingdom, 001. [1] T.E. Tezduyar, Computation of moving boundaries and interfaces and stabilization parameters, Int. J. Numer. Methods Fluids 43 (003) [13] T. Tezduyar, S. Sathe, Stabilization parameters in SUPG and PSPG formulations, J. Computat. Appl. Mech. 4 (003) [14] T.E. Tezduyar, Finite element methods for fluid dynamics with moving boundaries and interfaces, in: E. Stein, R. De Borst, T.J.R. Hughes (Eds.), Encyclopedia of Computational Mechanics. vol. 3: Fluids, John Wiley & Sons, 004 (Chapter 17). [15] T.E. Tezduyar, Stabilized finite element methods for computation of flows with moving boundaries and interfaces, in: Lecture Notes on Finite Element Simulation of Flow Problems, Japan Society of Computational Engineering and Sciences, Tokyo, Japan, 003. [16] T.E. Tezduyar, Stabilized finite element methods for flows with moving boundaries and interfaces, HERMIS: Int. J. Comput. Math. Appl. 4 (003) [17] T.E. Tezduyar, Determination of the stabilization and shock-capturing parameters in SUPG formulation of compressible flows, in: Proceedings of the European Congress on Computational Methods in Applied Sciences and Engineering, ECCOMAS 004 (CD-ROM), Jyvaskyla, Finland, 004. [18] T.E. Tezduyar, M. Behr, J. Liou, A new strategy for finite element computations involving moving boundaries and interfaces the deforming-spatial-domain/space time procedure: I. The concept and the preliminary tests, Comput. Methods Appl. Mech. Engrg. 94 (199)
12 163 T.E. Tezduyar, M. Senga / Comput. Methods Appl. Mech. Engrg. 195 (006) [19] T.E. Tezduyar, M. Behr, S. Mittal, J. Liou, A new strategy for finite element computations involving moving boundaries and interfaces the deforming-spatial-domain/space time procedure: II. Computation of free-surface flows, two-liquid flows, and flows with drifting cylinders, Comput. Methods Appl. Mech. Engrg. 94 (199) [0] T.E. Tezduyar, Y. Osawa, Finite element stabilization parameters computed from element matrices and vectors, Comput. Methods Appl. Mech. Engrg. 190 (000) [1] S.K. Aliabadi, T.E. Tezduyar, Parallel fluid dynamics computations in aerospace applications, Int. J. Numer. Methods Fluids 1 (1995)
ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS
ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 c ECCOMAS ADAPTIVE DETERMINATION OF THE FINITE ELEMENT STABILIZATION PARAMETERS Tayfun E. Tezduyar Team for Advanced
More informationSUPG STABILIZATION PARAMETERS CALCULATED FROM THE QUADRATURE-POINT COMPONENTS OF THE ELEMENT-LEVEL MATRICES
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 24 P. Neittaanmäki, T. Rossi, K. Majava, and O. Pironneau (eds.) W. Rodi and P. Le Quéré (assoc. eds.) Jyväskylä,
More informationFINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS
FINITE ELEMENT SUPG PARAMETERS COMPUTED FROM LOCAL DOF-MATRICES FOR COMPRESSIBLE FLOWS Lucia Catabriga Department of Computer Science, Federal University of Espírito Santo (UFES) Av. Fernando Ferrari,
More informationFinite element stabilization parameters computed from element matrices and vectors
Comput. Methods Appl. Mech. Engrg. 190 (2000) 411±430 www.elsevier.com/locate/cma Finite element stabilization parameters computed from element matrices and vectors Tayfun E. Tezduyar *, Yasuo Osawa Mechanical
More informationStabilized Formulations and Smagorinsky Turbulence Model for Incompressible Flows
WCCM V Fifth World Congress on Computational Mechanics July 7 2, 22, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Stabilized Formulations and Smagorinsky Turbulence Model for
More informationSolution techniques for the fully discretized equations in computation of fluid structure interactions with the space time formulations
Comput. Methods Appl. Mech. Engrg. 95 (2006) 5743 5753 www.elsevier.com/locate/cma Solution techniques for the fully discretized equations in computation of fluid structure interactions with the space
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More informationSpace-time finite element computation of compressible flows involving moving boundaries and interfaces*
Computer Methods in Applied Mechanics and Engineering 107 (1993) 209-223 North-Holland CMA 409 Space-time finite element computation of compressible flows involving moving boundaries and interfaces* S.K.
More informationComputer modeling of cardiovascular fluid structure interactions with the deforming-spatial-domain/stabilized space time formulation
Comput. Methods Appl. Mech. Engrg. 195 (2006) 1885 1895 www.elsevier.com/locate/cma Computer modeling of cardiovascular fluid structure interactions with the deforming-spatial-domain/stabilized space time
More informationA unified finite element formulation for compressible and incompressible flows using augmented conservation variables
. 2% mm d BB ELSEVIER Comput. Methods Appl. Mech. Engrg. 161 (1998) 229-243 Computer methods in applied mechanics and englneerlng A unified finite element formulation for compressible and incompressible
More informationComputational Methods and Applications in Parachute Maneuvers
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM 04, Sept. 5-10, 2004, Beijing, China c 2004 Tsinghua University Press & Springer-Verlag Computational Methods and Applications in Parachute Maneuvers
More informationImproved discontinuity-capturing finite element techniques for reaction effects in turbulence computation
Comput. Mech. (26) 38: 356 364 DOI.7/s466-6-45- ORIGINAL PAPER A. Corsini F. Rispoli A. Santoriello T.E. Tezduyar Improved discontinuity-capturing finite element techniques for reaction effects in turbulence
More informationAerodynamic Interactions Between Parachute Canopies
K. Stein Department of Physics, Bethel College, St. Paul, MN 55112 T. Tezduyar V. Kumar S. Sathe Mechanical Engineering, Rice University, MS 321, Houston, TX 77005 R. Benney Natick Soldier Center, Natick,
More informationOn finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion
On finite element methods for 3D time dependent convection diffusion reaction equations with small diffusion Volker John and Ellen Schmeyer FR 6.1 Mathematik, Universität des Saarlandes, Postfach 15 11
More informationT.E. Tezduyar, S. Sathe, R. Keedy. K. Stein
III Congreso Internacional sobre Metodos Numericos en Ingenieria y Ciencias Aplicadas S. Gallegos, I. Herrera, S. Botello, F. Zarate, y G. Ayala Editores) c ITESM, Monterrey 004 CIMNE, Barcelona 004 SPACE
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 219 On finite element methods for 3D time dependent convection diffusion reaction equations
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationStabilized and Coupled FEM/EFG Approximations for Fluid Problems
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM4, Sept. 5-, 24, Beijing, China c 24 Tsinghua University Press & Springer-Verlag Stabilized and Coupled FEM/EFG Approximations for Fluid Problems
More informationOn the interaction between dynamic model dissipation and numerical dissipation due to streamline upwind/petrov Galerkin stabilization
Comput. Methods Appl. Mech. Engrg. 194 (25) 1225248 www.elsevier.com/locate/cma On the interaction between dynamic model dissipation and numerical dissipation due to streamline upwind/petrov Galerkin stabilization
More informationEvaluating the LCD algorithm for solving linear systems of equations arising from implicit SUPG formulation of compressible
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2003; 00: 22 [Version: 2002/09/8 v2.02] Evaluating the LCD algorithm for solving linear systems of equations arising
More informationDue Tuesday, November 23 nd, 12:00 midnight
Due Tuesday, November 23 nd, 12:00 midnight This challenging but very rewarding homework is considering the finite element analysis of advection-diffusion and incompressible fluid flow problems. Problem
More informationStabilized methods for compressible flows
Stabilized methods for compressible flows Thomas J. R. Hughes Institute for Computational Sciences and Engineering, The University of Texas at Austin, 1 University Station C0200, Austin, TX 78712-0027,
More informationThe two-dimensional streamline upwind scheme for the convection reaction equation
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 575 591 The two-dimensional streamline upwind scheme for the convection reaction equation Tony W. H. Sheu*,1
More informationParachute uid±structure interactions: 3-D computation
Comput. Methods Appl. Mech. Engrg. 190 (2000) 373±386 www.elsevier.com/locate/cma Parachute uid±structure interactions: 3-D computation Keith Stein a, *, Richard Benney a, Vinay Kalro b, Tayfun E. Tezduyar
More informationMeshfree Petrov-Galerkin Methods for the Incompressible Navier-Stokes Equations
Meshfree Petrov-Galerkin Methods for the Incompressible Navier-Stokes Equations Thomas-Peter Fries and Hermann Georg Matthies Institute of Scientific Computing, Technical University of Braunschweig, Brunswick,
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More information^Laboratorio Nacional de Computaqao Cientifica, Rua Lauro Muller 455, Rio de Janeiro, RJ , Brasil
A space-timefiniteelement formulation for the shallow water equations F.L.B. Ribeiro,' A.C. Galeao," L. Landau* "Programa de Engenharia Civil, COPPE / Universidade Federal do Rio de Janeiro, Caixa Postal
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationInexact Newton-type methods for the solution of steady incompressible viscoplastic flows with the SUPG/PSPG finite element formulation
Comput. Methods Appl. Mech. Engrg. 195 (2006) 3145 3167 www.elsevier.com/locate/cma Inexact Newton-type methods for the solution of steady incompressible viscoplastic flows with the SUPG/PSPG finite element
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
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 informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationFINITE ELEMENT FORMULATION FOR THE VORTICITY -STREAM FUNCTION FORM OF THE INCOMPRESSIBLE EULER EQUATIONS ON MULTIPLY-CONNECTED DOMAINS* T.E.
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING 73 (1989) 331-339 NORTH-HOLLAND FINITE ELEMENT FORMULATION FOR THE VORTICITY -STREAM FUNCTION FORM OF THE INCOMPRESSIBLE EULER EQUATIONS ON MULTIPLY-CONNECTED
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More informationComputation of Cardiovascular Fluid-Structure Interactions with the DSD/SST Method
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM 04, Sept. 5-10, 2004, Beijing, China c 2004 Tsinghua University Press & Springer-Verlag Computation of Cardiovascular Fluid-Structure Interactions
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationA variational multiscale stabilized formulation for the incompressible Navier Stokes equations
Comput Mech (2009) 44:145 160 DOI 10.1007/s00466-008-0362-3 ORIGINAL PAPER A variational multiscale stabilized formulation for the incompressible Navier Stokes equations Arif Masud Ramon Calderer Received:
More informationSimulation of free surface fluids in incompressible dynamique
Simulation of free surface fluids in incompressible dynamique Dena Kazerani INRIA Paris Supervised by Pascal Frey at Laboratoire Jacques-Louis Lions-UPMC Workshop on Numerical Modeling of Liquid-Vapor
More informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationSTABILIZED FINITE ELEMENT METHODS OF GLS TYPE FOR MAXWELL-B AND OLDROYD-B VISCOELASTIC FLUIDS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
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 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 informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationNumerical solution of hyperbolic heat conduction in thin surface layers
International Journal of Heat and Mass Transfer 50 (007) 9 www.elsevier.com/locate/ijhmt Numerical solution of hyperbolic heat conduction in thin surface layers Tzer-Ming Chen * Department of Vehicle Engineering,
More informationCompatible algorithms for coupled flow and transport q
Comput. Methods Appl. Mech. Engrg. 193 (24) 2565 258 www.elsevier.com/locate/cma Compatible algorithms for coupled flow and transport q Clint Dawson *, Shuyu Sun, Mary F. Wheeler Center for Subsurface
More informationY Zβ Discontinuity Capturing for Advection-Dominated Processes with Application to Arterial Drug Delivery
Y Zβ Discontinuity Capturing for Advection-Dominated Processes with Application to Arterial Drug Delivery Y. Bazilevs a,1, V.M. Calo a,1, T.E. Tezduyar b,2, and T.J.R. Hughes a,3 a Institute for Computational
More informationApplication of multi-scale finite element methods to the solution of the Fokker Planck equation
Comput. Methods Appl. Mech. Engrg. 194 (2005) 1513 1526 www.elsevier.com/locate/cma Application of multi-scale finite element methods to the solution of the Fokker Planck equation Arif Masud a,1, Lawrence
More informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4607 SEMESTER: Autumn 2012-13 MODULE TITLE: Introduction to Fluids DURATION OF
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationPerformance evaluation of different model mixers by numerical simulation
Journal of Food Engineering 71 (2005) 295 303 www.elsevier.com/locate/jfoodeng Performance evaluation of different model mixers by numerical simulation Chenxu Yu, Sundaram Gunasekaran * Food and Bioprocess
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 informationSlip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics
Slip flow boundary conditions in discontinuous Galerkin discretizations of the Euler equations of gas dynamics J.J.W. van der Vegt and H. van der Ven Nationaal Lucht- en Ruimtevaartlaboratorium National
More informationIntroduction. Statement of Problem. The governing equations for porous materials with Darcy s law can be written in dimensionless form as:
Symbolic Calculation of Free Convection for Porous Material of Quadratic Heat Generation in a Circular Cavity Kamyar Mansour Amirkabir University of technology, Tehran, Iran, 15875-4413 mansour@aut.ac.ir
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 informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationSub-Cell Shock Capturing for Discontinuous Galerkin Methods
Sub-Cell Shock Capturing for Discontinuous Galerkin Methods Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 39, U.S.A. A shock capturing strategy for higher order
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 informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6. Mathematik Preprint Nr. 228 A Variational Multiscale Method for Turbulent Flow Simulation with Adaptive Large Scale
More informationForchheimer law derived by homogenization of gas flow in turbomachines
Journal of Computational and Applied Mathematics 215 (2008) 467 476 www.elsevier.com/locate/cam Forchheimer law derived by homogenization of gas flow in turbomachines Gottfried Laschet ACCESS e.v., IntzestraYe
More informationA STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION
A STABILIZED FINITE ELEMENT METHOD FOR INCOMPRESSIBLE VISCOUS FLOWS USING A FINITE INCREMENT CALCULUS FORMULATION Eugenio Oñate International Centre for Numerical Methods in Engineering Universidad Politécnica
More informationA hierarchical multiscale framework for problems with multiscale source terms
Available online at www.sciencedirect.com Comput. Methods Appl. Mech. Engrg. 197 (2008) 2692 2700 www.elsevier.com/locate/cma A hierarchical multiscale framework for problems with multiscale source terms
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
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 informationA Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods
ScientifiComputing A Review of Petrov-Galerkin Stabilization Approaches and an Extension to Meshfree Methods Thomas-Peter Fries, Hermann G. Matthies Institute of Scientific Computing Technical University
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
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 informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationEuclidean invariance and weak-equilibrium condition for the algebraic Reynolds stress model
J. Fluid Mech. 2006), vol. 569, pp. 399 408. c 2006 Cambridge University Press doi:10.1017/s0022112006003041 Printed in the United Kingdom 399 Euclidean invariance and weak-equilibrium condition for the
More informationReview of fluid dynamics
Chapter 2 Review of fluid dynamics 2.1 Preliminaries ome basic concepts: A fluid is a substance that deforms continuously under stress. A Material olume is a tagged region that moves with the fluid. Hence
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationIMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1
Computational Methods in Applied Mathematics Vol. 1, No. 1(2001) 1 8 c Institute of Mathematics IMPROVED LEAST-SQUARES ERROR ESTIMATES FOR SCALAR HYPERBOLIC PROBLEMS, 1 P.B. BOCHEV E-mail: bochev@uta.edu
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 15 (21) 3965 3973 Contents lists available at ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns A similarity solution in
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationA generalized method for advective-diffusive computations in engineering
Fluid Structure Interaction and Moving Boundary Problems 563 A generalized method for advective-diffusive computations in engineering H. Gómez, I. Colominas, F. Navarrina & M. Casteleiro Department of
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationComputational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor
Computational and Experimental Studies of Fluid flow and Heat Transfer in a Calandria Based Reactor SD Ravi 1, NKS Rajan 2 and PS Kulkarni 3 1 Dept. of Aerospace Engg., IISc, Bangalore, India. ravi@cgpl.iisc.ernet.in
More informationNumerical simulation of compressible viscous flow
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 15, Issue 1 Ver. IV (Jan. - Feb. 2018), PP 44-56 www.iosrjournals.org Numerical simulation of compressible
More informationANALYSIS OF FLOW IN A CONCENTRIC ANNULUS USING FINITE ELEMENT METHOD
Nigerian Journal of Technology (NIJOTECH) Vol 35, No 2, April 2016, pp 344 348 Copyright Faculty of Engineering, University of Nigeria, Nsukka, Print ISSN: 0331-8443, Electronic ISSN: 2467-8821 wwwnijotechcom
More informationNUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD
NUMERICAL SIMULATION OF THE FLOW AROUND A SQUARE CYLINDER USING THE VORTEX METHOD V. G. Guedes a, G. C. R. Bodstein b, and M. H. Hirata c a Centro de Pesquisas de Energia Elétrica Departamento de Tecnologias
More informationCommun Nonlinear Sci Numer Simulat
Commun Nonlinear Sci Numer Simulat 7 (0) 3499 3507 Contents lists available at SciVerse ScienceDirect Commun Nonlinear Sci Numer Simulat journal homepage: www.elsevier.com/locate/cnsns Application of the
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
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 informationFLUID MECHANICS. Atmosphere, Ocean. Aerodynamics. Energy conversion. Transport of heat/other. Numerous industrial processes
SG2214 Anders Dahlkild Luca Brandt FLUID MECHANICS : SG2214 Course requirements (7.5 cr.) INL 1 (3 cr.) 3 sets of home work problems (for 10 p. on written exam) 1 laboration TEN1 (4.5 cr.) 1 written exam
More informationRevisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model
Revisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model R. S. Myong a and S. P. Nagdewe a a Dept. of Mechanical and Aerospace Engineering
More informationA NEW MODELING OF SHEET CAVITATION CONSIDERING THE THERMODYNAMIC EFFECTS. Yumiko Sekino Ebara Cooporation, Ltd., Ota-ku, Tokyo JAPAN
Cav3-GS-6-3 Fifth International Symposium on Cavitation (CAV3) Osaka, Japan, November -4, 3 A NEW MODELING OF SHEET CAVITATION CONSIDERING THE THERMODYNAMIC EFFECTS Takashi Tokumasu Institute of Fluid
More information- Marine Hydrodynamics. Lecture 4. Knowns Equations # Unknowns # (conservation of mass) (conservation of momentum)
2.20 - Marine Hydrodynamics, Spring 2005 Lecture 4 2.20 - Marine Hydrodynamics Lecture 4 Introduction Governing Equations so far: Knowns Equations # Unknowns # density ρ( x, t) Continuity 1 velocities
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationFinite Element Model of a complex Glass Forming Process as a Tool for Control Optimization
Excerpt from the Proceedings of the COMSOL Conference 29 Milan Finite Element Model of a complex Glass Forming Process as a Tool for Control Optimization Felix Sawo and Thomas Bernard Fraunhofer Institute
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 information