( ) A i,j. Appendices. A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in
|
|
- Annis Stanley
- 5 years ago
- Views:
Transcription
1 Appendices A. Sensitivity of the Van Leer Fluxes The flux Jacobians of the inviscid flux vector in Eq.(3.2), and the Van Leer fluxes in Eq.(3.11), can be found in the literature [9,172,173] and are therefore omitted here. The sensitivity of the Van Leer fluxes, for higher-order spatially accurate discrete sensitivity analysis on unstructured grids, has been performed and presented for the first time in the current research. These sensitivities may be derived by noting that for unstructured grids there are two means by which the mesh dependence on the design variables may influence the flux vector; the metric terms, N { η x η y η z } T, the cell face area, A f, and the upwind interpolation scheme, Q f { ρ u v w p} T. These dependencies may be symbolically written for the residual vector as R jκ ( i) + E i,j N N + E i,j N + N + E i,j Q f Q f + E i,j + A i,j Q f Q f + ( ) A i,j + E + i,j + E i,j (A.1) Recall from section 3.1.2, that the determination of the appropriate fluxes are based upon the Mach number normal to the cell face. For supersonic flow through the cell face, the flux vector given by Eq.(3.2) is used with the primitive variables interpolated to the cell interfaces. The sensitivity of this flux vector with respect to the metric dependence on the design variables, i.e., N η x η y T η z, may be expressed as 117
2 ρ Θ β k Q f ρu Θ + p η x β Q f k E N ρv Θ + p η y N β Q f k ρw Θ + p η z β Q f k ( ρe o + p) Θ β k Q f (A.2a) and with respect to the dependence of the upwind interpolation on the design variables Q f ρ u v w p T as E Q f Q f ρ Θ + ρ Θ N ρ uθ + ρ u Θ + ρu Θ + p η x β N k ρ vθ + ρ v Θ + ρv Θ + p η β y N k ρ wθ + ρ w Θ + ρw Θ + p η N β y k ( ) ρe o + p Θ + ( ρe o + p) Θ N (A.2b) where the sensitivity of the normal velocity is the combination Θ Θ + Θ β Q f k N (A.3a) with 118
3 Θ Q f u η x + v η y + w η z ; Θ N u η x + v η y + w η z (A.3b) and the sensitivity of the total enthalpy is given by ( ρe o + p) γ p + ρ u 2 + v 2 + w 2 + ρ u u + v v + w w γ 1 2 (A.4) For subsonic flow through the cell face, the split Van Leer fluxes in Eq.(3.11) are utilized. The sensitivity of these fluxes may be expressed as E N N + E Q f Q f f mass f mass f mass f mass f mom1 f mom2 f mom3 f mass f + f mom1 mass f mom2 + f mass + f mass f mom3 f f e + f e mass (A.5) where the sensitivity of the mass flux contains the contributions f mass f mass + f mass β Q k f N (A.6a) with f mass Q f ρa 2 M n 1 ( ) M n (A.6b) Q f f mass N 1 ρ 4 a M β n 1 k ( ) 2 + ρ a ( ) 2 + 2ρa( M n 1) M n M n 1 N (A.6c) The sensitivity of the f mom1 term is the combination 119
4 f mom1 f mom1 + f mom1 β Q k f N (A.7a) where f mom1 η x β Q k f ( Θ 2a) γ η x γ Θ Q f (A.7b) f mom1 N u + η x γ Θ 2 a β N k (A.7c) Analogous relations may be written for f mom2 sensitivity for the energy term is and f mom3. Similarly the f e f e Q f + f e N (A.8a) where f e 2 [ Q γ 2 1 ( 1 γ )Θ ( γ 1)a ] Θ f Q f (A.8b) f e N 2 ( γ 2 1 γ )Θ ( γ 1)a 1 { } Θ N + u u + v v + w w + { 2a ( γ 1)Θ } a (A.8c) In the above equations the sensitivities of the normal Mach number are given by M n 1 Θ ; a β Q f k Q f M n N 1 Θ a N Θ a a 2 (A.9) where the sensitivity for the speed of sound is a γ 2a 1 ρ p p ρ 2 ρ (A.10) 120
5 As seen, the sensitivity of the inviscid split fluxes requires the determination of the sensitivity of the metric terms and the sensitivity of the upwind interpolation scheme. The evaluation of the metric terms on unstructured grids and the corresponding sensitivity of these terms are presented below in Appendix C. The sensitivity of the upwind interpolation scheme, using both the inverse-distance and the psuedo-laplacian weighting methods, are presented in Appendix D. B. Jacobians/Sensitivity of the Boundary Conditions The boundary-condition types utilized in the current work are inviscid solid boundary (flow tangency) and characteristic inflow/outflow. These boundary conditions are those most commonly used in inviscid CFD analysis and may be easily found in the literature; for example, see reference 128. Furthermore, a first-order implementation of these boundary conditions is utilized in the current research. This implementation simply uses the cell center values of the primitive variables adjacent to the boundary faces as opposed to the interpolation of these variables to the boundary faces. A more detailed description and comparison of both implementations has been reported Ref. 129 for unstructured grid schemes. The most common approach, however, is the first-order one. The boundary values of the primitive variables Q b for flow tangency are specified as u b u o Θ o η x ; ρ b ρ o v b v o Θ o η y ; w b w o Θ o η z p b p o (B.1) where the normal velocity at the boundary is computed from the interior velocity components as Θ o u o η x + v o η y + w o η z (B.2) The state equation for the boundary conditions may be express as 121
6 B( Q o,q b,x ) ρ o ρ b u o Θ o η x u b v o Θ o η y v b 0 (B.3) w o Θ o η z w b p o p b The Jacobian of this state equation, with respect to the interior state variables, is derived as B Q o η x η x η y η x η z 0 0 η x η y 2 1 η y η y η z η x η z η y η z 1 η z (B.4) and with respect to the state variables on the boundary B Q b I [ ] B Q b 1 (B.5) where I is the identity matrix. The inverse of this matrix is required in the sensitivity analysis using either the pre-eliminated form of the sensitivity equation (e.g., Eqs.(3.26a and b)) or the incremental form (e.g., Eq.(4.12)). The sensitivity of the boundary state equation with respect to the design variables for the flow tangency condition may be written as 0 Θ o η η x + Θ x o B Θ η o η y + Θ y o Θ o η η z + Θ z o 0 (B.6) where 122
7 Θ o u o η x + v o η y + w o η z (B.7) For characteristic inflow/outflow boundary conditions, the Riemann invariants corresponding to the incoming and outgoing waves traveling in the characteristic directions defined normal to the boundary are applied at the far-field. The two locally one-dimensional Riemann invariants are given by where Θ o is given above and R + Θ o + 2a o γ 1 ; R Θ 2a γ 1 (B.8) a o γ p o ρ o 1 2 ; a γ p ρ 1 2 ; Θ u η x + v η y + w η z (B.9) The invariants are used to determine the locally normal velocity and speed of sound ( ) ; a γ 1 ( R + R ) (B.10) Θ 1 2 R+ + R For outflow (Θ > 0) on the boundary, the velocities are defined as u b u o + η x ( Θ Θ o ) ; v b v o +η y ( Θ Θ o ) ; w b w o +η z ( Θ Θ o ) (B.11) whereas for inflow (Θ < 0) on the boundary u b u + η x ( Θ Θ ) ; v b v +η y ( Θ Θ ) ; w b w +η z ( Θ Θ ) (B.12) The density on the boundary is computed from the entropy relation 4 ρ b a 2 γ S 1 γ 1 (B.13) with the entropy function S evaluated as Θ > 0 S S o p o γ ; Θ < 0 S S p γ ρ o ρ (B.14) 123
8 Then the corresponding pressure on the boundary is computed from the equation of state p b ρ b a 2 γ (B.15) Denoting for outflow ς o and for inflow ς, these boundary conditions may be written in state equation form as B( Q o,q b,x ) a 2 γ S ς u ς + η x Θ Θ ς v ς + η y Θ Θ ς w ς + η z Θ Θ ς ρ b a 2 p b γ 1 γ 1 ρb ( ) u b ( ) v b ( ) w b (B.16) The Jacobians of the state equation with respect to the interior state vector, for characteristic inflow/outflow boundary conditions may be derived as B Q o ξ 1 φ 1 Θ Θ η x 1+ η ρ x φ 3 η x φ 4 η x φ 5 η x o p o Θ Θ η y η y φ 3 1+ η y φ 4 η y φ 5 η y ρ o p o Θ Θ η z η z φ 3 η z φ 4 1+ η z φ 5 η z ρ o p o a a a a a ξ 2 ξ 2 ξ 2 ξ 2 ξ 2 ρ o u o v o w o p o ξ 1 a u o ξ 1 a v o ξ 1 a w o ξ 1 φ 2 (B.17a) where ξ 1 2a γ γ 1 γ S ς a 2 ( ) S ς 2 γ γ 1 (B.17b) 124
9 φ 1 a a S ς δ ρ o 2S ς ; ς ρ o φ 2 a p o a 2S ς δ ς S ς p o (B.17c) and φ 3 Θ u o δ ς Θ ς u o ; φ 4 Θ v o δ ς Θ ς v o ξ 2 2 ρ b a γ ; φ 5 Θ w o δ ς Θ ς w o (B.17d) (B.17e) where δ ς 1 for outflow (ς o) and δ ς 0 for inflow (ς ). The derivatives of the local normal velocity and speed of sound are Θ Q o γ p o Θ ρ o 2 2( γ 1)a o ρ o Θ u o η x 2 Θ v o η y 2 ; Θ w o η z 2 Θ p γ o ( γ 1)a o ρ o a γ 1 Θ Q o 2 Q o (B.17f) The derivative of the normal velocity for inflow ( ς ) is zero, and for outflow For the state variables on the boundary Θ ς o 0 η Q x η y η z 0 o { } T (B.17g) B I Q b [ ] + Elem(5,1); Elem(5,1) a 2 γ (B.18) Once again, the inverse of this matrix is needed in the sensitivity analysis, and from the above form, the inverse is simply B Q b 1 [ I ] Elem(5,1) (B.19) That is, only a change in sign of the element in location (5,1) is required for the inverse. 125
10 The sensitivity of the boundary state equation with respect to the design variables is given by 2 γ 2a a 2 γ 1 a γ ( γ 1)S ς γ S ς η x Θ Θ Θ β ς k ( ) +η x Θ ς B η y Θ Θ Θ β ς k ( ) +η y Θ ς η z Θ Θ Θ β ς k ( ) +η z Θ ς 2 ρ b a a γ β k (B.20) with the local speed of sound and normal velocity sensitivity given by a γ 1 u 4 o u Θ 1 2 ( ) η x u o + u ( ) η x ( ) η y + v o v ( ) η y + v o + v ( ) η z + w o w ( ) η z + w o + w (B.21) (B.22) and where Θ ς 1 2 u ς η x η y η + v β ς + w z k β ς k (B.23) is computed based on either an inflow or outflow condition. The sensitivity of the boundary state equation, like the sensitivity of the flux vectors for the interior cell state equation discussed in Appendix A, require the metric sensitivity derivatives. The evaluation of the metric terms on an unstructured grid, and the corresponding sensitivity, are discussed in Appendix C. 126
11 C. Sensitivity of the Metric Terms Expressions for the metric terms, and the corresponding sensitivity of these expressions, may be derived for a typical tetrahedral cell (see Fig. C.1) by first defining the edge vectors L 2 l 2 x L 3 l 3 x L 4 l 4 x y z { l 2 l 2 } T x n2 x n1 y z { l 3 l 3 } T x n3 x n1 y z { l 4 l 4 } T x n4 x n1 { } T (C.1a) ( ) ( y n2 y n1 ) ( z n2 z n1 ) { } T (C.1b) ( ) ( y n3 y n1 ) ( z n3 z n1 ) { } T (C.1c) ( ) ( y n4 y n1 ) ( z n4 z n1 ) where the sensitivities of these vectors with respect to the shape design variables are L 2 l x y z T 2 l 2 l 2 L 3 l x y z T 3 l 3 l 3 L 4 l x y T z 4 l 4 l 4 (C.2a) (C.2b) (C.2c) where, for example, l 2 x x n2 x n1 ; l 2 y y n2 y n1 ; l 2 z z n2 z n1 (C.3) In Eq.(C.3) the quantities x n, y n, and z n are the grid sensitivity terms discussed in Chapter 5. Similar expressions may be written for components of L 3 and L 4. The area of face n1-n2-n3 may be computed as A f 1 2 L 2 L a 2 2 x + a y + 2 [ az ] 12 ( ) 2 + l x 3 l z 2 l x z ( 2 l 3 ) 2 + l x 2 l y 3 l x y ( 3 l 2 ) 2 1 l y 2 2 l z 3 l z y 2 l 3 with the sensitivity of the face area given by 12 (C.4) 127
12 A f A f 1 a a x a y a 4 A x + a f β y + a z k β z k (C.5) where a x l y 2 z l β 3 + y l z 3 l2 l z 2 y l k β 3 + z l y 3 l2 k a y l x 3 z l β 2 + x l z 2 l3 l x 2 z l k β 3 + x l z 3 l2 k a z l x 2 y l β 3 + x l y 3 l2 l x 3 y l k β 2 + x l y 2 l3 k (C.6a) (C.6b) (C.6c) The normal direction, referred to as the direction cosines, for this face may be evaluated as N { η x η y η z } T L 2 L 3 L 2 L 3 (C.7) L 2 L 3 2 A f with the corresponding sensitivity N N 1 2 L 2 L 3 + L 2 L 3 1 ( A f ) 2 A f L 2 L 3 A f (C.8) Once again, similar expressions may be easily written for the other faces of the tetrahedron. The centroid of the tetrahedral cell is determined as x c 1 ( 4 x n1 + x n2 + x n3 + x n4 ) (C.9a) y c 1 ( 4 y n1 + y n2 + y n3 + y n4 ) (C.9b) where z c 1 ( 4 z n1 + z n2 + z n3 + z n4 ) (C.9c) 128
13 x c 1 4 x n1 + x n2 + x n3 + x n4 (C.10a) y c 1 4 y n1 + y n2 + y n3 + y n 4 (C.10b) z c 1 4 z n1 + z n2 + z n3 + z n4 (C.10c) are the coordinate sensitivities of the centroid location. D. Sensitivity of the Spatial Differencing It was discussed in section 3.1.3, that the development of a higher-order spatially accurate scheme requires the interpolation of the state variables to the cell interfaces. For unstructured grid schemes, this interpolation is mesh dependent. Thus, the sensitivity of the interpolation given in Eq.(3.15) may be written as Q f [( Q + Q r v ) ] [( Q v r ) ] (D.1a) introducing the expression for the solution gradient in Eq.(3.16) yields Q f Q n1 + Q n2 + Q n3 Q n4 (D.1b) where the variables at the nodes are obtained from the multidimensional weighted averaging given in equation The sensitivity of this averaging with respect to the design variables is given by Q n nc w c,i i1 Q c,i nc w c,i i1 nc wc,i Q c,i i1 nc w c,i i1 nc w c,i i1 2 (D.2) 129
14 The sensitivity of the weighting factors, w c,i, depends on the algorithm used. In the current research, an inverse-distance and a psuedo-laplacian weighting procedure have been utilized. The sensitivity of each scheme is presented below. D.1 Inverse Distance Weighting The weighting factors for the inverse-distance procedure were given in equation The sensitivity of these factors may be written as w c,i x ( x) + y ( y ) + z ( z) ( x 2 + y 2 + z 2 ) 32 (D.3) with ( x) x c,i x n ; ( y ) y c,i y n ; ( z ) z c,i z n (D.4) where the sensitivity of the coordinates of the centroid location are given in equations C.10a through C.10c. D.2 Psuedo-Laplacian Weighting The weighting factors for this method have been given previously in equations 3.19a and 3.19b. Sensitivity of the weighting factors may be expressed as ( ) ( ) ( ) w c,i λ x x x + λ β x + λ y y y + λ k β y + λ z z z + λ k β z k (D.5) where ( x), ( y), and ( z) have been given above in equation D.4. The derivatives of the Lagrange multipliers may be expressed as λ x v Ω β I v y I v z k ( ) Ω v v I y I v z + I v y I v z (D.6a) 130
15 λ y v Ω β I v x I v z k ( ) + Ω v I v x I v z + I v x I v z λ z v Ω β I v x v k ( I y ) v I v x Ω I v y + v I x v I y (D.6b) (D.6c) On examining the components of v Ω, v I x, v I y, and v I z, which are given in Eqs.(3.20d-i), it can be observed that the derivatives of these components only include the derivatives of the centroid location and nodal coordinates. Once again, the derivatives of the centroid location are given above in equations C.10a to C.10c. E. Sensitivity of Common Output Functions For aerodynamic calculations, typical output functions from which objective functions and constraints may be defined are the lift coefficient, drag coefficient, and lift-to-drag ratio. Each will be discussed to follow. Lift Coefficient The lift coefficient is computed as C L F z cosα F x sinα q A ref (E.1) where α is the free-stream angle of attack, q is the dynamic pressure, A ref is the reference area, and the forces in the z- and x-directions are nbf nbf F z η z A f ( p j p ) ; F x η x A f ( p j p ) (E.2) j 1 where nbf is the number of boundary faces over which the pressure is integrated. To compute the sensitivity derivatives in Eqs.(2.6a and b), the derivatives with respect to the state vector and the design variables are required. The derivatives of the lift coefficient are j1 131
16 C L Q i nbf δ ij A f j ( η z cosα η x sinα) q A ref (E.3) where δ ij is the Kronecker delta, and C L nbf ( p j p ) η z A f cosα η x A f j1 q sinα A ref A ref (E.4) with η z A f A ref A ref η z A A f + η f A z η z A ref f 2 (E.5) A ref A similar expression may be written for η x A f. Note, if the reference area is fixed A ref throughout the design then the sensitivity of the reference area is zero. If the reference area used is the actual wetted surface area then the sensitivity of the reference area becomes A ref nbf A f j1 ; A ref A f nbf j1 (E.6) where the sensitivity of the face area has been given in Eq.(C.5) and the sensitivity of the metric terms are given in equation C.8. Drag Coefficient The drag coefficient is computed as C D F x cosα + F z sinα q A ref (E.7) 132
17 with the derivatives of the drag coefficient given by and C D C D 0 Q i nbf δ ij A f ( η x cosα + η z sinα) j1 q A ref nbf ( p j p ) η x A f cosα + η z A f j1 q sinα A ref A ref (E.8) (E.9) Note that the above assumes that the side slip angle is zero, but may be easily incorporated. Lift-to-Drag Ratio The derivatives of the lift-to-drag ratio may be expressed in terms of the derivatives of the lift coefficient and drag coefficient as follows ( C L C D ) Q C L Q C L C D C D Q C D (E.10) similarly ( ) C L C D C L C L C D C D C D (E.11) where all the terms in Eqs.(E.10 and E.11) have been given above. 133
18 n 4 L 4 n 2 n 3 L 2 L 3 n 1 Figure C.1: Typical unstructured grid tetrahedral cell. 134
Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationEntry Aerodynamics MARYLAND U N I V E R S I T Y O F. Entry Aerodynamics. ENAE Launch and Entry Vehicle Design
Atmospheric Regimes on Entry Basic fluid parameters Definition of Mean Free Path Rarified gas Newtonian flow Continuum Newtonian flow (hypersonics) 2014 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu
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 informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
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 informationEstimation of grid-induced errors in computational fluid dynamics solutions using a discrete error transport equation
Graduate Theses and Dissertations Graduate College 2009 Estimation of grid-induced errors in computational fluid dynamics solutions using a discrete error transport equation Brandon Williams Iowa State
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics and Astronautics AIAA Aerospace Sciences Meeting, Reno, NV AIAA Paper 2007-0037
More informationDeforming Composite Grids for Fluid Structure Interactions
Deforming Composite Grids for Fluid Structure Interactions Jeff Banks 1, Bill Henshaw 1, Don Schwendeman 2 1 Center for Applied Scientific Computing, Lawrence Livermore National Laboratory, Livermore,
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 informationLecture 5.7 Compressible Euler Equations
Lecture 5.7 Compressible Euler Equations Nomenclature Density u, v, w Velocity components p E t H u, v, w e S=c v ln p - c M Pressure Total energy/unit volume Total enthalpy Conserved variables Internal
More informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
More informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
More informationEntropy stable schemes for compressible flows on unstructured meshes
Entropy stable schemes for compressible flows on unstructured meshes Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore deep@math.tifrbng.res.in http://math.tifrbng.res.in/
More informationEffect of Applied Magnetic Field on Shock Boundary Layer Interaction
University of Kentucky UKnowledge Mechanical Engineering Faculty Publications Mechanical Engineering 1-2012 Effect of Applied Magnetic Field on Shock Boundary Layer Interaction Ovais U. Khan University
More informationVectors 1C. 6 k) = =0 = =17
Vectors C For each problem, calculate the vector product in the bracet first and then perform the scalar product on the answer. a b c= 3 0 4 =4i j 3 a.(b c=(5i+j.(4i j 3 =0 +3= b c a = 3 0 4 5 = 8i+3j+6
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
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 informationSimulation of unsteady muzzle flow of a small-caliber gun
Advances in Fluid Mechanics VI 165 Simulation of unsteady muzzle flow of a small-caliber gun Y. Dayan & D. Touati Department of Computational Mechanics & Ballistics, IMI, Ammunition Group, Israel Abstract
More informationTHERMODYNAMICALLY BASED MOISTURE PREDICTION USING ROE S SCHEME
THERMODYNAMICALLY BASED MOISTURE PREDICTION USING ROE S SCHEME M. J. Kermani 1, A. G. Gerber 2, J. M. Stockie 3 University of New Brunswick Fredericton, New Brunswick, Canada, E3B 5A3 M.Kermani@unb.ca
More informationMultigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations. Todd A.
Multigrid Solution for High-Order Discontinuous Galerkin Discretizations of the Compressible Navier-Stokes Equations by Todd A. Oliver B.S., Massachusetts Institute of Technology (22) Submitted to the
More informationFOR many computational fluid dynamics (CFD) applications, engineers require information about multiple solution
AIAA Aviation 13-17 June 2016, Washington, D.C. 46th AIAA Fluid Dynamics Conference AIAA 2016-3809 Improved Functional-Based Error Estimation and Adaptation without Adjoints William C. Tyson *, Katarzyna
More informationEffects of the Jacobian Evaluation on Direct Solutions of the Euler Equations
Middle East Technical University Aerospace Engineering Department Effects of the Jacobian Evaluation on Direct Solutions of the Euler Equations by Ömer Onur Supervisor: Assoc. Prof. Dr. Sinan Eyi Outline
More informationOn Numerical Simulation of Flow Problems by Discontinuous Galerkin and Finite Volume Techniques
On Numerical Simulation of Flow Problems by Discontinuous Galerkin and Finite Volume Techniques Ondřej Winter 1, Petr Sváček 1 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Technical
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationDetailed Numerical Study on the Aerodynamic Behavior of a Gurney Flapped Airfoil in the Ultra-low Reynolds Regime
doi: 10.5028/jatm.v9i2.631 Detailed Numerical Study on the Aerodynamic Behavior of a Gurney Flapped Airfoil in the Ultra-low Reynolds Regime Mahzad Khoshlessan 1, Seyed Mohammad Hossein Karimian 1 Abstract:
More information. D CR Nomenclature D 1
. D CR Nomenclature D 1 Appendix D: CR NOMENCLATURE D 2 The notation used by different investigators working in CR formulations has not coalesced, since the topic is in flux. This Appendix identifies the
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
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 informationNumerical Simulation of Jet Aerodynamics Using the Three-Dimensional Navier-Stokes Code PAB3D
NASA Technical Paper 3596 Numerical Simulation of Jet Aerodynamics Using the Three-Dimensional Navier-Stokes Code PAB3D S. Paul Pao Langley Research Center Hampton, Virginia Khaled S. Abdol-Hamid Analytical
More informationHow Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm
AIAA 2 2673 How Many Steps are Required to Solve the Euler Equations of Steady, Compressible Flow: In Search of a Fast Solution Algorithm Antony Jameson Stanford University Stanford, CA 9435 D. A. Caughey
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationGeneral elastic beam with an elastic foundation
General elastic beam with an elastic foundation Figure 1 shows a beam-column on an elastic foundation. The beam is connected to a continuous series of foundation springs. The other end of the foundation
More informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationAdjoint-based aerodynamic shape optimization
IT Licentiate theses 2003-012 Adjoint-based aerodynamic shape optimization OLIVIER AMOIGNON UPPSALA UNIVERSITY Department of Information Technology Adjoint-based aerodynamic shape optimization BY OLIVIER
More informationADJOINT AND DEFECT ERROR BOUNDING AND CORRECTION FOR FUNCTIONAL ESTIMATES
ADJOINT AND DEFECT ERROR BOUNDING AND CORRECTION FOR FUNCTIONAL ESTIMATES Niles A. Pierce and Michael B. Giles Applied & Computational Mathematics, California Institute of Technology Computing Laboratory,
More informationAn Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations. Thesis by Cheng Chi
An Improved Ghost-cell Immersed Boundary Method for Compressible Inviscid Flow Simulations Thesis by Cheng Chi In Partial Fulfillment of the Requirements For the Degree of Masters of Science King Abdullah
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 informationA high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows
A high-order discontinuous Galerkin solver for 3D aerodynamic turbulent flows F. Bassi, A. Crivellini, D. A. Di Pietro, S. Rebay Dipartimento di Ingegneria Industriale, Università di Bergamo CERMICS-ENPC
More informationOn limiting for higher order discontinuous Galerkin method for 2D Euler equations
On limiting for higher order discontinuous Galerkin method for 2D Euler equations Juan Pablo Gallego-Valencia, Christian Klingenberg, Praveen Chandrashekar October 6, 205 Abstract We present an implementation
More informationTotal Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation
Total Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation D.I. Papadimitriou, K.C. Giannakoglou Laboratory of Thermal Turbomachines, National Technical University
More informationAnisotropic mesh refinement for discontinuous Galerkin methods in two-dimensional aerodynamic flow simulations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2008; 56:2111 2138 Published online 4 September 2007 in Wiley InterScience (www.interscience.wiley.com)..1608 Anisotropic
More informationShock Capturing for Discontinuous Galerkin Methods using Finite Volume Sub-cells
Abstract We present a shock capturing procedure for high order Discontinuous Galerkin methods, by which shock regions are refined in sub-cells and treated by finite volume techniques Hence, our approach
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationMAE SUMMER 2015 HOMEWORK 1 SOLUTION
MAE 04 - SUMMER 205 HOMEWORK SOLUTION Problem :. Plot the drag coefficient, c D, the aerodynamic efficiency, AE, and the center of pressure, x cp, of the airplane as a function of the angle of attack.
More informationCFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing. Lakhdar Remaki
CFD in Industrial Applications and a Mesh Improvement Shock-Filter for Multiple Discontinuities Capturing Lakhdar Remaki Outline What we doing in CFD? CFD in Industry Shock-filter model for mesh adaptation
More informationRANS Equations in Curvilinear Coordinates
Appendix C RANS Equations in Curvilinear Coordinates To begin with, the Reynolds-averaged Navier-Stokes RANS equations are presented in the familiar vector and Cartesian tensor forms. Each term in the
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationComparison of outflow boundary conditions for subsonic aeroacoustic simulations
Comparison of outflow boundary conditions for subsonic aeroacoustic simulations A. Fosso P., H. Deniau, N. Lamarque and T. Poinsot CERFACS, CFD Team, Av. G. Coriolis, 3157 Toulouse Cedex France SUMMARY
More informationSubsonic characteristic outflow boundary conditions with transverse terms: a comparative study
Center for Turbulence Research Annual Research Briefs 200 35 Subsonic characteristic outflow boundary conditions with transverse terms: a comparative study By G. Lodato, F. Ham AND H. Pitsch. Motivation
More informationIntroduction to numerical simulation of fluid flows
Introduction to numerical simulation of fluid flows Mónica de Mier Torrecilla Technical University of Munich Winterschool April 2004, St. Petersburg (Russia) 1 Introduction The central task in natural
More informationWe have read this thesis and recommend its acceptance:
To the Graduate Council: I am submitting herewith a thesis written by Nicholas G. Currier entitled A Hybrid Method for Flows in Local Chemical Equilibrium and Nonequilibrium. I have examined the final
More informationExplicit Incompressible Euler Solver
Explicit Incompressible Euler Solver (Written In Matlab) R D Teja Aerospace Undergraduate IIT Madras Contents: 1. Abstract 2. Introduction 2.1 Euler equations 2.2 Artificial compressibility equations 2.3
More informationNewtonian Analysis of Rarified Flows
Atmospheric Regimes on Entry Basic fluid parameters Definition of Mean Free Path Rarified gas Newtonian flow Continuum Newtonian flow (hypersonics) SphereConeAero so ware 2012 David L. Akin - All rights
More informationParallel Computation of Forced Vibration for A Compressor Cascade
44th AIAA Aerospace Sciences Meeting and Exhibit 9-2 January 26, Reno, Nevada AIAA 26-628 AIAA Paper 26-628 Parallel Computation of Forced Vibration for A Compressor Cascade Zongjun Hu and Gecheng Zha
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 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 informationReduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis
Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper
More informationTRANSONIC AIRFOIL CALCULATIONS USING THE EULER EQUATIONS
TRANSONIC AIRFOIL CALCULATIONS USING THE EULER EQUATIONS Antony Jameson Princeton University Princeton, NJ 08544 I am going to tell you about some of the transonic aerofoil calculations I have been doing
More informationWave propagation methods for hyperbolic problems on mapped grids
Wave propagation methods for hyperbolic problems on mapped grids A France-Taiwan Orchid Project Progress Report 2008-2009 Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan ISCM
More informationSynchronous Machine Modeling
ECE 53 Session ; Page / Fall 07 Synchronous Machine Moeling Reference θ Quarature Axis B C Direct Axis Q G F D A F G Q A D C B Transient Moel for a Synchronous Machine Generator Convention ECE 53 Session
More informationIntroduction to Finite Volume projection methods. On Interfaces with non-zero mass flux
Introduction to Finite Volume projection methods On Interfaces with non-zero mass flux Rupert Klein Mathematik & Informatik, Freie Universität Berlin Summerschool SPP 1506 Darmstadt, July 09, 2010 Introduction
More informationHYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS
1 / 36 HYPERSONIC AERO-THERMO-DYNAMIC HEATING PREDICTION WITH HIGH-ORDER DISCONTINOUS GALERKIN SPECTRAL ELEMENT METHODS Jesús Garicano Mena, E. Valero Sánchez, G. Rubio Calzado, E. Ferrer Vaccarezza Universidad
More informationSIMULATION OF UNSTEADY TURBOMACHINERY FLOWS USING AN IMPLICITLY COUPLED NONLINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 2011 GT2011 June 6-10, 2011, ancouver, Canada GT2011-46367 SIMULATION OF UNSTEADY TURBOMACHINERY FLOWS USING AN IMPLICITLY COUPLED NONLINEAR HARMONIC BALANCE METHOD Jonathan
More informationMatrices and Determinants
Chapter1 Matrices and Determinants 11 INTRODUCTION Matrix means an arrangement or array Matrices (plural of matrix) were introduced by Cayley in 1860 A matrix A is rectangular array of m n numbers (or
More informationAeroacoustic and Aerodynamics of Swirling Flows*
Aeroacoustic and Aerodynamics of Swirling Flows* Hafiz M. Atassi University of Notre Dame * supported by ONR grant and OAIAC OVERVIEW OF PRESENTATION Disturbances in Swirling Flows Normal Mode Analysis
More informationForce analysis of underwater object with supercavitation evolution
Indian Journal of Geo-Marine Sciences Vol. 42(8), December 2013, pp. 957-963 Force analysis of underwater object with supercavitation evolution B C Khoo 1,2,3* & J G Zheng 1,3 1 Department of Mechanical
More informationTwo-Dimensional Riemann Solver for Euler Equations of Gas Dynamics
Journal of Computational Physics 167, 177 195 (2001) doi:10.1006/jcph.2000.6666, available online at http://www.idealibrary.com on Two-Dimensional Riemann Solver for Euler Equations of Gas Dynamics M.
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationComparison of nonreflecting outlet boundary. conditions for compressible solvers on unstructured. grids
Comparison of nonreflecting outlet boundary conditions for compressible solvers on unstructured grids Victor Granet 1 CERFACS, CFD Team, 42 avenue Gaspard Coriolis, 31057 Toulouse Cedex, France and Renault
More informationSpectral Difference Method for Unstructured Grids II: Extension to the Euler Equations
Journal of Scientific Computing, Vol. 32, No. 1, July 2007 ( 2006) DOI: 10.1007/s10915-006-9113-9 Spectral Difference Method for Unstructured Grids II: Extension to the Euler Equations Z. J. Wang, 1,4
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationAn Adjoint Method Augmented with Grid Sensitivities for Aerodynamic Optimization
An Adjoint Method Augmented with Grid Sensitivities for Aerodynamic Optimization by Chad Oldfield A thesis submitted in conformity with the requirements for the degree of M.A.Sc. Graduate Department of
More informationFLOW OVER A GIVEN PROFILE IN A CHANNEL WITH DYNAMICAL EFFECTS
Proceedings of Czech Japanese Seminar in Applied Mathematics 2004 August 4-7, 2004, Czech Technical University in Prague http://geraldine.fjfi.cvut.cz pp. 63 72 FLOW OVER A GIVEN PROFILE IN A CHANNEL WITH
More informationOn a class of numerical schemes. for compressible flows
On a class of numerical schemes for compressible flows R. Herbin, with T. Gallouët, J.-C. Latché L. Gastaldo, D. Grapsas, W. Kheriji, T.T. N Guyen, N. Therme, C. Zaza. Aix-Marseille Université I.R.S.N.
More informationApplication Of Faraday s Law
Application Of Faraday s Law Dr Miguel Cavero September 2, 2014 Application Of Faraday s Law September 2, 2014 1 / 23 The PHYS120 Exam will be divided into three sections as follows: Section A: Short Questions
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationY. Abe, N. Iizuka, T. Nonomura, K. Fujii Corresponding author:
Seventh International Conference on Computational Fluid Dynamics (ICCFD7), Big Island, Hawaii, July 9-13, 2012 ICCFD7-2012-2801 ICCFD7-2801 Symmetric-conservative metric evaluations for higher-order finite
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationAppendix B. Coordinate Transformation
B Appendix B. Coordinate Transformation B.1. Coordinate Transformation Relations In this appendix we derive and discuss the transformation relations between curvilinear and orthogonal coordinate systems,
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 informationA Matrix Stability Analysis of the Carbuncle Phenomenon
A Matrix Stability Analysis of the Carbuncle Phenomenon Michael Dumbser a Jean-Marc Moschetta b Jérémie Gressier c a Universität Stuttgart, Institut für Aerodynamik und Gasdynamik, Pfaffenwaldring 21,
More informationAPPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM
AIAA 2003-3693 APPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM Vladimir V. Golubev* and Axel Rohde Embry-Riddle Aeronautical University Daytona
More informationCapSel Roe Roe solver.
CapSel Roe - 01 Roe solver keppens@rijnh.nl modern high resolution, shock-capturing schemes for Euler capitalize on known solution of the Riemann problem originally developed by Godunov always use conservative
More informationClassical differential geometry of two-dimensional surfaces
Classical differential geometry of two-dimensional surfaces 1 Basic definitions This section gives an overview of the basic notions of differential geometry for twodimensional surfaces. It follows mainly
More informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
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 informationNACA 0006 A NACA 0012 CT5
2012 10 30 5 Journal of Northwestern Polytechnical University Oct. Vol. 30 2012 No. 5 710072 NACA 0006 NACA 0012 Basic Finner MissileBFM 1. 58 ~ 2. 5 V211. 3 A 1000-2758201205-0784-05 CFD CFD 1 1-3 2 CFD
More informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationThe Most Marvelous Theorem in Mathematics. Dan Kalman American University
The Most Marvelous Theorem in Mathematics Dan Kalman American University www.dankalman.net Outline Overview of the theorem Ellipses Background Facts Proof Real Polynomials 5 4 3 2 Familiar functions: 4x
More informationActive Flux for Advection Diffusion
Active Flux for Advection Diffusion A Miracle in CFD Hiroaki Nishikawa National Institute of Aerospace! NIA CFD Seminar! August 25, 2015 In collaboration with the University of Michigan Supported by NASA
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationComparison of Approximate Riemann Solvers
Comparison of Approximate Riemann Solvers Charlotte Kong May 0 Department of Mathematics University of Reading Supervisor: Dr P Sweby A dissertation submitted in partial fulfilment of the requirement for
More informationImmersed boundary technique for compressible flow simulations on semi-structured meshes
Center for Turbulence Research Annual Research Briefs 2005 71 Immersed boundary technique for compressible flow simulations on semi-structured meshes By M. de Tullio AND G. Iaccarino 1. Motivation and
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 informationStable and high-order accurate finite difference schemes on singular grids
Center for Turbulence Research Annual Research Briefs 006 197 Stable and high-order accurate finite difference schemes on singular grids By M. Svärd AND E. van der Weide 1. Motivation and objectives The
More informationUncertainty Analysis of Computational Fluid Dynamics Via Polynomial Chaos
Uncertainty Analysis of Computational Fluid Dynamics Via Polynomial Chaos Rafael A. Perez Dissertation submitted to the faculty of the Virginia Polytechnic Institute and State University in partial fulfillment
More informationPLAXIS. Scientific Manual
PLAXIS Scientific Manual 2016 Build 8122 TABLE OF CONTENTS TABLE OF CONTENTS 1 Introduction 5 2 Deformation theory 7 2.1 Basic equations of continuum deformation 7 2.2 Finite element discretisation 8 2.3
More information