Antony Jameson. Stanford University Aerospace Computing Laboratory Report ACL
|
|
- Regina Parks
- 6 years ago
- Views:
Transcription
1 Formulation of Kinetic Energy Preserving Conservative Schemes for Gas Dynamics and Direct Numerical Simulation of One-dimensional Viscous Comressible Flow in a Shock Tube Using Entroy and Kinetic Energy Preserving Schemes Antony Jameson Stanford University Aerosace Comuting Laboratory Reort ACL 007 Moana Surfrider, Waikiki March 007
2 Abstract This aer follows u on the author s recent aer The Construction of Discretely Conservative Finite Volume Schemes that also Globally Conserve Energy or Enthaly. In the case of the gas dynamics equations the revious formulation leads to an entroy reserving EP) scheme. It is shown in the resent aer that it is also ossible to construct the flux of a conservative finite volume scheme to roduce a kinetic energy reserving KEP) scheme which exactly satisfies the global conservation law for kinetic energy. A roof is resented for three dimensional discretization on arbitrary grids. Both the EP and KEP schemes have been alied to the direct numerical simulation of one-dimensional viscous flow in a shock tube. The comutations verify that both schemes can be used to simulate flows with shock waves and contact discontinuities without the introduction of any artificial diffusion. The KEP scheme erformed better in the tests.
3 Introduction Stemming from the ioneering work of Godunov [], rocedures for the construction of nonoscillatory shock caturing schemes are by now well established [, 3, 4, 5, 6]. In general they add artificial diffusion either exlicitly ot imlicitly via uwind oerators in order to satisfy total variation diminishing TVD) or local extremum diminishing LED) roerties [7, 8]. However there is a risk that the artificial diffusion may comromise the accuracy of viscous flow simulations. If the discrete scheme can be constructed to satisfy a global energy estimate of some kind, then it should by stable, at least for smooth solutions, without the need for artificial viscosity. The use of energy estimates to establish stability has a long history, and it is discussed in the classical book of Morton and Richtmyer [9]. One route to achieving discrete energy estimates is to use difference oerators which have a skew-symmetric form. Tyically these oerators are derived by slitting the governing equations in a mixture of conservation and quasilinear form [0]. However, in the treatment of comressible flows with shock waves it is beneficial to write the discrete equations in conservation form. According to the theorem of Lax and Wendroff [], this is sufficient to guarantee that the discrete solution satisfies the correct shock jum conditions as long as it converges in the limit as the mesh interval is reduced. In inviscid comressible fluid flow the total energy, consisting of the sum of the internal and kinetic energy, is conserved, not just the kinetic energy. Corresondingly, it follows from the governing equations that entroy is conserved if no shock waves are formed. In a recent aer [] the author showed that a semi-discrete scheme in conservation form can be constructed so that it globally conserves a generalized entroy function in a smooth flow. Such an entroy reserving EP) scheme is obtained by an aroriate formulation of the numerical fluxes across the interfaces between the grid cells. It requires the use of entroy variables in the evaluation of the flux, although the standard conservative variables are udated at each time ste. There is some latitude in the definition of the generalized entroy function hu) of the state vector u. Harten [3] has given conditions such hu) is a convex function, so that the solution cannot become unbounded if hu) remains bounded. Multilying the governing equations for u by t wt = h then roduces the evolution equation u for h, while t wt reresents the entroy variables. Entroy variables have been used by Hughes, Mallet and Franca [4], and also by Gerritsen and Olsson [5], who roosed a non conservative entroy reserving scheme. 3
4 Although a bound on the kinetic energy does not assure a bound on the solution in a comressible flow, correct simulation of the evolution of kinetic energy is a crucial requirement for accurate simulations of turbulence, where there is an energy cascade between the different eddy scales. In the resent work it is shown that the interface fluxes of a semidiscrete conservative scheme can be constructed in an alternative way which assures that the global discrete kinetic energy evolves in a manner that exactly corresonds to the true equation for kinetic energy. There is some latitude in the definition of the fluxes for such a kinetic energy reserving KEP) scheme, rovided that the fluxes for the continuity and momentum equations satisfy a comatibility condition. Section resents the derivation of the entroy reserving EP) and kinetic energy reserving KEP) schemes for the one dimensional gas dynamics equations. In Section 3 the formulation of the KEP scheme is extended to multi-dimensional viscous comressible flow simulations on arbitrary grids. The KEP roerty is again assured by a comatibility condition between the fluxes for the continuity and momentum equations. The roof is comleted by showing that the finite volume scheme satisfies a discrete Gauss theorem for the ressure and stress terms, rovided that these are both evaluated at the cell interfaces by arithmetic averages of their values in the neighboring cells. In Section 4 both the EP and the KEP schemes are alied to the direct numerical simulation DNS) of one dimensional viscous flow in a shock tube. It is demonstrated in numerical exeriments that both schemes can successfully resolve the shock wave, contact discontinuity and exansion fan without adding any artificial diffusion, rovided that a fine enough mesh is used with a number of cells of the order of the Reynolds number. Reresentative numerical calculations include simulations at a Reynolds number in the range based on the seed of sound on meshes with cells. 4
5 Entroy and kinetic energy reserving schemes for the one-dimensional gas dynamics equations This section resents the formulation of fully conservative schemes for one dimensional gas dynamics which are either entroy reserving EP) or kinetic energy reserving KEP). In both cases the global conservation roerty is obtained by roer construction of the interface flux between each air of neighboring cells. The detailed roof for the EP scheme has been given in [], but in order to facilitate the comarison of the EP and KEP schemes it is outlined here Consider the gas dynamics equations in the conservation form Here the state and flux vectors are u = ρ ρv ρe u t + fu) = 0.) x, f = ρv ρv + ρvh.) where ρ is the density, v is the velocity and, E and H are the ressure, energy and enthaly. Also = γ )ρ where γ is the ratio of secific heats. ) E v, H = E + ρ.3) In the absence of shock waves the entroy ) s = log ρ γ.4) is constant, satisfying the advection equation Consider the generalized entroy function s t + v s x = 0.5) hs) = ρgs).6) 5
6 where it has been shown by Harten [3] that h is a convex function of u rovided that / d g dg ds ds < γ.7) Then h satisfies the entroy conservation law t hu) + F u) = 0.8) x where the entroy flux is F = ρvgs).9) Moreover, introducing the entroy variables w T = h u.0) it can be verified that h u f u = F u and hence on multilying.) by w T we recover the entroy conservation law.8) where now the Jacobian matrix f w = f uu w is symmetric. Accordingly f can be exressed as the gradient of a scalar function G, f = G w.) and the entroy flux can be exressed as F = f T w G.) Different choices of the entroy function gs) have been discussed by Harten [3], Hughes, Franca and Mallet [4], and Gerritsen and Olsson [5]. Suose now that.) is aroximated in semi-discrete form on a grid with cell intervals x j, j =, n as where the numerical flux f j+ du j x j dt + f j+ f j = 0.3) is a function of u i over a range of i bracketing j. In order 6
7 to construct an entroy reserving EP) scheme multily.3) by w T and sum by arts to obtain n x j wj T j= du j dt = n j= w T j f j+ f j n = w T f wn T f n+ + j= ) f T j+ w j+ w j ) At interior oints evaluate f T j+ as the mean value of G wj+ in the sense of Roe [4] such that G wj+ w j+ w j ) = G w j+ ) G w j ).4) Also evaluate the boundary fluxes as f = f w ), f n+ = f w n ).5) Then the interior fluxes cancel, and using.0) amd.), we obtain the entroy conservation law in the discrete form n j= x j dh j dt = F w ) F w n ).6) G wj+ can be constructed to satisfy.4) exactly by evaluating it as the integral G wj+ = 0 G w ŵθ)) dθ.7) where ŵθ) = w j + θ w j+ w j ).8) since then G w j+ ) G w j ) = = 0 0 G w ŵθ)) w θ dθ G w ŵθ)) dθ w j+ w j ) Thus we obtain: 7
8 Theorem. The semi-discrete conservation law.3) satisfies the semi-discrete entroy conservation law.6) is the numerical flux is calculated as f j+ = 0 fŵθ)dθ, j =, n where ŵθ) is defined by.8), and the boundary fluxes are defined by.5) The construction of a kinetic energy reserving KEP) scheme requires a different aroach in which the fluxes of the continuity and momentum equations are searately constructed in a comatible manner. Denoting the secific kinetic energy by k, Thus ] k = ρ v, k [ u = v, v, 0 k t = v v ρ ρv) t t = { )} v + ρ v x + v x.9) Suose that the semi-discrete conservation scheme.3) is written searately for the continuity and momentum equations as dρ j x j dt + ρv) j+ ρv) j d x j dt ρv) j + ρv ) j+ ρv ) j + j+ j = 0.0) = 0.) 8
9 Now multilying.0) by v j and.) by v j, adding them and summing by arts, n j= d x j v j dt ρv) j v j = n j= v j n j= = v ρv) n + j= j= ρv j ) j+ ) dρ j = dt v j j+ j n j= ρv j ) j ) d v ) j x j ρ j dt ) n j= + v ρv ) + v + v n { ρv) j+ v j+ v j v j ρv ) j+ ρv ) j ρv) n+ v n ρv ) n+ } ) ρv ) j+ v j+ v j ) ) v n n+ n + j+ v j+ v j ).) Each term in the first sum containing the convective terms can be exanded as { v j+ + v j ρv) j+ } ρv ) j+ v j+ v j ) and will vanish if ρv ) j+ Now evaluating the boundary fluxes as v j+ + v j = ρv) j+.3) ρv) ρv) n+ = ρ v, ρv ) = ρ n v n, ρv ) n+ = ρ v, = ρ n vn, n+ =.4) = n.) reduces to the semi-discrete kinetic energy conservation law n j= v ) j x j ρ j ) ) v vn = v + ρ v n n + ρ n + n j= j+ v j+ + v j ).5) 9
10 Denoting the arithmetic average of any quantity q between j + and j as q = q j+ + q j ) the interface ressure may be evaluated as Also if one sets j+ =.6) ρv) j+ ) ρv j+ = ρ v.7) = ρ v.8) condition.3) is satisfied. Consistently one may set The foregoing argument establishes ρvh) j+ = ρ v H.9) Theorem. The semi-discrete conservation law.3) satisfies the semi-discrete kinetic energy global conservation law.5) if the fluxes for the continuity and momentum equations satisfy condition.3) and the boundary fluxes are calculated by equations.4). Condition.3) allows some latitude in the construction of the fluxes. For examle it is also satisfied if one sets ρv) j+ ) ρv j+ ρvh) j+ = ρv = ρv v = ρv H instead of equations.7).9). 0
11 3 Kinetic energy reserving KEP) scheme for multidimensional viscous flow The extension of the entroy reserving EP) scheme to multi-dimensional gas dynamics on general grids has been given in []. In this section it is shown how to extend the kinetic energy reserving KEP) scheme to multi-dimensional viscous comressible flow. Denote the density, velocity comonents, ressure, energy and enthaly by, v i,, E and H, where the suerscrit i is used to denote the i th coordinate direction. This is convenient because subscrits will be needed to denote the grid location. Reeated suerscrits will be used to indicate summation over the coordinate directions. Also ) The viscous stress tensor is = γ ) ρe vi, H = E + ρ 3.) ) v σ ij i = µ x + vj ij vk + λδ 3.) j x i x k where δ ij is the Kronecker delta, and µ and λ are the viscosity coefficients. Usually λ = 3 µ. The heat flux is q j = κ T x j 3.3) where T is the temerature and κ is the coefficient of heat conduction. equations can now be written in conservation form as The governing u t + x f i u) = 0 3.4) i where the state and flux vectors are u = ρ ρv ρv ρv 3 ρe, f i = ρv i ρv i v σ i + δ i ρv i v σ i + δ i ρv i v 3 σ i3 + δ i3 ρv i H v j σ ij q j 3.5)
12 The kinetic energy is again denoted by k where [ ] k = ρ vi, k u = vi, v, v, v 3, 0 3.6) and the kinetic energy conservation law is k t + ) } {v j + ρ vi + v i σ ij x j = vj vi σij 3.7) xj x j Suose now that the domain is covered by a grid, and the equations are discretized in finite volume form. The discrete variables may be associated either with nodes, each of which is surrounded by a control volume, or with the control volumes themselves. In the first case the control volumes may be taken as dual cells where the nodes are the vertices of the rimary cells. In the second case the control volumes are simly the rimary cells, and the discrete variables may be regarded as cell averages. Each interior control volume, say o, is bounded by faces not necessarily lanar) with a directed face area S o for the face searating the control volumes o and. Each boundary control volume is closed by an outer face with directed area S o which is the negative of the sum S o of the face areas between o and its neighbors. The control volumes may be tetrahedral, hexahedral or of mixed olyhedral form. The semi-discrete finite volume scheme to be considered has the form vol o du o dt + S i of i o = 0 3.8) for every interior control volume, where S i o are the rojected areas of the face S o in the coordinate directions. At a boundary control volume b there is an additional contribution Sb if b i, where the boundary flux vector will be evaluated as f i b = f i u b ) 3.9) Multilying 3.8) by w T o = ) k u o 3.0)
13 and summing over the nodes o vol o dk o dt = o w T o So i fo i b w T b S i b f i b 3.) where the last sum reresents the contribution of the boundaries. Every interior face aears twice in the sums on the right hand side of 3.) with oosite signs for its discrete face area. Thus its contribution is ) w T wo T S i o fo i 3.) Consider now contributions of the convective terms to 3.). these are ) v j vo j S i o ρv i v j) v o j ) vo j S ) o i ρv i o = { ) v j vo j S i ρv i o v j) o Thus the convective contributions of the interior faces will vanish if ρv i ) o v j ) } vo j ρv i v j) o = ρv i ) o v j + v j o) 3.3) Defining the average of any quantity q between its values at o and as q = q + q o ) 3.4) condition 3.3) will be satisfied if one set ρv i ) o = ρ v i 3.5) ρv i v j) o = ρ v i v j 3.6) To comlete the derivation it is now necessary to examine the contributions of the ressure and stress tensor to 3.). These can be exressed as vo T o SoP i o i b v T b S i bp i b 3.7) 3
14 where v = v v v 3, P i = δ i σ i δ i σ i δ i3 σ i3 3.8) Let P i o be evaluated as The first term in 3.7) is now P i o = P i + P i o) = P i 3.9) vo T Po i So i o o v T o PS i o i At every interior control volume o the fluxes v T SoP i o i generated by its neighbors can be associated with o, so the contributions at o can be written as vt o Po i So i + vt S i o Po i Here Su o = 0 because it is the sum of the face areas of a closed volume. Thus one can add vo T Po i Si o to obtain P o it v + v o ) S i o = σ ij o This reresents the finite volume discretization of v j So i + o v i So i vi vj σij xi x i for the control volume o. A boundary control volume b receives the contributions P b it v Sb i from the neighbors while it returns the contributions P it b v b Sb i Pb it v b Sb i 4
15 giving the total P b it v + v b ) S i o P it b = P it b = σ ij b v b v + v b ) S i b + P it b Sb i Pb it v b Sb i ) v j Sb i + v j b Si b + P b v b Sb i Pb it v b Sb i v i S i b + v i bs i b ) σ ij b vj b Si b P v i bs i b Combining this result with the result for the interior control columes, and including the convective terms for the boundaries, we finally obtain the semi-discrete global kinetic energy conservation law o vol o dk o dt = b + o S i o { v i o o vo i o + ρ o ) v i S i o σ ij o v i oσ ij o } ) v j So i 3.0) where the first sum on the right hand side over the boundary control volumes reresent the flux through the boundaries and the second sum is the finite volume discretization of This establishes theorem 3. D ) vi vj σij dv xi x i Theorem 3. The semi-discrete conservation law 3.8) satisfies the semi-discrete global kinetic energy conservation law 3.0) if the interface fluxes satisfy condition 3.3) and the boundary fluxes are evaluated by equation 3.9). Note that the discrete viscous terms o σ ij v j S i o guarantee dissiation of the discrete energy because vj S i o is the consistent discretization of vj x i in the control volume. Then slitting vj x i into its symmetric and anti-symmetric arts 5
16 as v j x = ) v j i x + vi + ) v j i x j x vi i x j the contribution of the viscous terms amounts to { ) ) ) v j µ x + vi v j i x j x + vi v k } + λ i x j x k o Setting λ = µ this is 3 o v j µ x + vi i x ) vk j δij 3 x k The discretization of the viscous terms covered by the roof does not have the most comact ossible stencil. Consequently it might allow an udamed odd-even mode. 6
17 4 Direct numerical solution of one-dimensional viscous flow in a shock tube This section resents the results of numerical exeriments in which both the entroy reserving EP) and the kinetic energy reserving KEP) schemes have been alied to the direct numerical simulation DNS) of one dimensional viscous flow in a shock tube. In an analysis of discrete solution methods for the viscous Burgers equation [] it was established that the EP scheme will satisfy the conditions for a local extremum diminishing LED) scheme if the local cell Reynolds number Re c. Here Re c is defined as v x ν where v is the velocity and ν is the kinematic viscosity. This indicates the need to use a mesh with a number of cells roortional to the global Reynolds number vl, where L is the ν global length scale. The comressible Navier Stokes equations are not amenable to such a simle analysis, but it can still be exected that the number of mesh cells needed to fully resolve shock waves and contact discontinuities will be roortional to the Reynolds number, given that the shock thickness is roortional to the coefficient of viscosity, as has been shown by G.I. Taylor and W.D. Hayes [6, 7]. In the numerical exeriments the viscous stress was calculated at each cell interface by discretizing the velocity and temerature gradient on a comact stencil as ) v x ) T x j+ j+ = x v j+ v j ) = x T j+ T j ) The coefficient of viscosity was calculated as a function of the temerature by Sutherland s law µ = T 3/ T + 0.3) 7
18 taking λ = µ, the viscous stress and heat flux assume the form 3 σ j+ q j+ = 4 ) u 3 µ x ) T = κ x where the coefficient of heat conduction κ was obtained from the Prandtl number P r as The Prandtl number was taken to be.75. κ = µc P r Numerical exeriments were erformed using three different flux formulas j+ j+. Simle averaging: f j+ = f u j+) + f u j )). The entroy reserving EP) scheme: f j+ = where w denote the entroy variables and 0 f ŵθ)) dθ ŵθ) = w j + θ w j+ w j ) 3. The kinetic energy reserving KEP) scheme: ρv) j+ = ρ v ) ρv = ρ v j+ ρvh) j+ = ρ v H In the EP scheme the entroy variables were taken to be w T = h u 8
19 where ) h = ρe s γ+ γ+ = ρ ρ γ Accordingly the entroy variables assume the comaratively simle form w = u 3 u u, u = w 3 w w where = γ ) γ + e s γ γ+ γ+ = γ + γ It was remarked in [] that the energy or entroy reserving roerty could be imaired by the time discretization scheme. One solution to this difficulty is to use an imlicit timesteing scheme of Crank Nicolson tye in which the satial derivatives are evaluated using the average value of the state vectors between the beginning and the end of each time ste, ū j = u n+ j ) + u n j This requires the use of inner iterations in each time ste. In order to avoid this cost, Shu s total variation diminishing TVD) scheme [8] was used for the time integration in all the numerical exeriments. Writing the semi-discrete scheme in the form du dt + Ru) = 0 4.) where Ru) reresents the discretized satial derivative, this advances the solution during one time ste by the three stage scheme u ) = u 0) t Ru 0) ) u ) = 3 4 u0) + 4 u) 4 t Ru) ) u 3) = 3 u0) + 3 u) 3 t Ru) ) The test case for the numerical exeriments was the well known examle originally roosed by Sod [9]. The shock tube extends over the range 0 x, with a discontinuity 9
20 in the initial data at x =.5. The left and right states are L =.0, R =. ρ L =.0, ρ R =.5 v L = 0, v R = 0 The Reynolds number is based on the seed of sound of the left state Re = ρ Lc L µ, c γl L = ρ L The numerical exeriments confirm that the EP and KEP schemes both enable direct numerical simulation DNS) of the viscous flow in the shock tube, rovided that a fine enough mesh is used, while significant oscillations can be observed in the solution when simle flux averaging is used scheme ). It is interesting that the oscillations are rimary observed in the exansion region. As a reresentative examle, Figures show the solutions for a Reynolds number Re = 5000 which were obtained with the three schemes on a grid with 4096 cells, at the time t =.36. The figures also dislay the exact inviscid solution with a solid line.with simle flux averaging there are oscillations in the entroy measured by 0 of the order of.0. With the EP scheme these oscillations are reduced to the ρ γ ρ γ 0 order of.00, while with the KEP scheme they are further reduced to the order of.000. When the Reynolds number is increased to million, and the calculations are erformed on a mesh with cells, it remains true that the only observable oscillations aear in the exansion region, but with reduced amlitude, while the three schemes exhibit the same order of merit. On the other hand the oscillations in the exansion region are amlified as the Reynolds number and number of mesh cells are reduced, as is illustrated in Figures In all cases surrisingly to the author) the KEP scheme erforms better than the EP scheme for this model roblem, while simle flux averaging is inadequate. 0
21 a) Pressure b) Density c) Velocity d) Energy Figure 4.: Simle averaging of the flux: 4096 mesh cells, Reynolds number 5000, Comuted solution values +, Exact inviscid solution
22 a) Pressure b) Density c) Velocity d) Energy Figure 4.: Entroy reserving scheme: 4096 mesh cells, Reynolds number 5000, Comuted solution values +, Exact inviscid solution
23 a) Pressure b) Density c) Velocity d) Energy Figure 4.3: Kinetic energy reserving scheme: 4096 mesh cells, Reynolds number 5000, Comuted solution values +, Exact inviscid solution 3
24 a) Pressure b) Density c) Velocity d) Energy Figure 4.4: Simle averaging of the flux coarse mesh): 5 mesh cells, Reynolds number 500, Comuted solution values +, Exact inviscid solution 4
25 a) Pressure b) Density c) Velocity d) Energy Figure 4.5: Entroy reserving scheme coarse mesh): 5 mesh cells, Reynolds number 500, Comuted solution values +, Exact inviscid solution 5
26 a) Pressure b) Density c) Velocity d) Energy Figure 4.6: Kinetic energy reserving scheme coarse mesh): 5 mesh cells, Reynolds number 500, Comuted solution values +, Exact inviscid solution 6
27 5 Conclusion As a sequel to [], the derivations in this aer establish that it is ossible to construct semi-discrete aroximations to the comressible Navier Stokes equations in conservation form which also discretely reserve the conservation of either entroy the EP scheme) or kinetic energy the KEP scheme). Both these schemes enable the direct numerical simulation of one dimensional viscous flow in a shock tube, rovided that the number of cells in the comutational mesh is of the order of the Reynolds number. The erformance of both the EP and the KEP schemes imroves as the Reynolds number and the number of mesh cells are simultaneously increased. For the model roblem examined in this aer, one-dimensional viscous flow in a shock tube, the KEP scheme erforms better than the EP scheme. The Kolmogoroff scale for the small eddies that can ersist in a viscous turbulent flow is of the order of Re 3 4. Accordingly it aears that by using a mesh with the order of Re 3 cells, direct numerical simulation DNS) of viscous turbulent flow with shock waves will be feasible in the future for high Reynolds number flows. Current high-end comuters attain comuting seeds of the order of 00 teraflos 0 4 floating oint oerations/second). This is about million times faster than high-end comuters 5 years ago. A further increase by a factor of million to 0 0 flos could enable DNS of viscous comressible flow at a Reynolds number of million. This is still short of the flight Reynolds numbers of long range transort aircraft in the range of million, but the eventual use of DNS for comressible turbulent flows can clearly be anticiated 7
28 Acknowledgement The author has benefited tremendously from the continuing suort of the Air Force Office of Scientific Research over the last fifteen years, most recently through Grant Number AF- F , under the direction of Dr. Fariba Fahroo. He is also indebted to Nawee Butsuntorn for his hel in rearing this work in L A TEX. 8
29 References [] S. K. Godunov, A Difference Method for the Numerical Calculation of Discontinous Solutions of Hydrodynamic Equations, Math. Sbornik, 47, 959, [] J. P. Boris and D. L. Book, Flux Corrected Transort, SHASTA, A Fluid Transort Algorithm that Works, J. Com. Phys.,, 973, [3] Bram Van Leer, Towards the Ultimate Conservative Difference Scheme, II, Monotomicity and Conservation Combined in a Second Order Scheme, J. Com. Phys., 4, 974, [4] P. L. Roe, Aroximate Reimann Solvers, Parameter Vectors and Difference Scheme, J. Com. Phys., 43, 98, [5] Amiram Harten, High Resolution Schemes for Hyerbolic Conservation Laws, J. Com. Phys., 49, 983, [6] M. S. Liou and C. J. Steffen, A New Flux Slitting Scheme, J. Com. Phys., 07, 993, 39. [7] Antony Jameson Analysis and Design of Numerical Schemes for Gas Dynamics Artificial Diffusion, Uwind Biasing, Limiters and Their Effect on Accuracy and Multigrid Convergence, International Journal of Comutational Fluid Dynamics, 4, 995, 7 8. [8] Antony Jameson Analysis and Design of Numerical Schemes for Gas Dynamics Artificial Diffusion and Discrete Shock Structure, International Journal of Comutational Fluid Dynamics, 5, 995, 38. [9] R. D. Richtmyer and K. W. Morton, Difference Methods for Initial Value Problems, Interscience, 967. [0] A.E. Honein and P. Moin, Higher Entroy Conservation and Numerical Stability of Comressible Turbulence Simulations, J. Com. Phys., 0, 004, [] P. D. Lax and B. Wendroff, Systems of Conservation Laws, Comm. Pure. Al. Math., 3, 960,
30 [] A. Jameson, The Construction of Discretely Conservative Finite Volume Schemes that Also Globally Conserve Energy or Entroy, Reort ACL 007, January 007. [3] Amiram Harten, On the Symmetric Form of Systems of Conservation Laws with Entroy, J. Com. Phys., 49, 983, [4] T. J. Hughes, L. P. Franca, and M. Mallet, A New Finite Element Formulation for Comutational Fluid Dynamics: I. Symmetric Forms of the Comressible Euler and Navier Stokes Equations and the Second Law of Thermodynamics, Comuter Methods in Alied Mechanics and Engineering, 54, 986, [5] Margot Gerritsen and Pelle Olsson, Designing an Efficient Solution Strategy for Fluid Flows, J. Com. Phys., 9, 996, [6] G.I. Taylor, The Conditions Necessary for Discontinuous Motion in Gases, Proc. Roy. Soc. London, A84, 90, [7] W.D. Hayes, Gasdynamic Discontinuities, Section D, Fundamentals of Gas Dynamics, edited by Howard W. Emmons, Princeton University Press, ). [8] C. W. Shu, Total-Variation-Diminishing Time Discretizations, SIAM J. Sci. Statist. Comuting, 9, 988, [9] G.A. Sod, Survey of several finite difference methods for systems of nonlinear hyerbolic conservation laws, J. Com. Phys., 6, 978, 3 30
Antony Jameson. AIAA 18 th Computational Fluid Dynamics Conference
Energy Estimates for Nonlinear Conservation Laws with Applications to Solutions of the Burgers Equation and One-Dimensional Viscous Flow in a Shock Tube by Central Difference Schemes Antony Jameson AIAA
More informationThe construction of discretely conservative finite volume schemes that also globally conserve energy or entropy
The construction of discretely conservative finite volume schemes that also globally conserve energy or entropy Antony Jameson Stanford University Aerospace Computing Laboratory Report ACL 007 January
More informationThe Construction of Discretely Conservative Finite Volume Schemes that Also Globally Conserve Energy or Entropy
J Sci Comput 008 34: 5 87 DOI 0.007/s095-007-97-7 The Construction of Discretely Conservative Finite Volume Schemes that Also Globally Conserve Energy or Entropy Antony Jameson Received: 30 July 007 /
More informationDirect Numerical Simulations of a Two-dimensional Viscous Flow in a Shocktube using a Kinetic Energy Preserving Scheme
19th AIAA Computational Fluid Dynamics - 5 June 009, San Antonio, Texas AIAA 009-3797 19th AIAA Computational Fluid Dynamics Conference, - 5 June 009, San Antonio, Texas Direct Numerical Simulations of
More informationUnderstanding DPMFoam/MPPICFoam
Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented
More information5. PRESSURE AND VELOCITY SPRING Each component of momentum satisfies its own scalar-transport equation. For one cell:
5. PRESSURE AND VELOCITY SPRING 2019 5.1 The momentum equation 5.2 Pressure-velocity couling 5.3 Pressure-correction methods Summary References Examles 5.1 The Momentum Equation Each comonent of momentum
More informationPaper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation
Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional
More informationChapter 6. Thermodynamics and the Equations of Motion
Chater 6 hermodynamics and the Equations of Motion 6.1 he first law of thermodynamics for a fluid and the equation of state. We noted in chater 4 that the full formulation of the equations of motion required
More informationRole of Momentum Interpolation Mechanism of the Roe. Scheme in Shock Instability
Role of Momentum Interolation Mechanism of the Roe Scheme in Shock Instability Xiao-dong Ren 1,2 Chun-wei Gu 1 Xue-song Li 1,* 1. Key Laboratory for Thermal Science and Power Engineering of Ministry of
More informationNUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS
NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS Tariq D. Aslam and John B. Bdzil Los Alamos National Laboratory Los Alamos, NM 87545 hone: 1-55-667-1367, fax: 1-55-667-6372
More informationDirect Numerical Simulations of Plunging Airfoils
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 010, Orlando, Florida AIAA 010-78 Direct Numerical Simulations of Plunging Airfoils Yves Allaneau
More informationWeek 8 lectures. ρ t +u ρ+ρ u = 0. where µ and λ are viscosity and second viscosity coefficients, respectively and S is the strain tensor:
Week 8 lectures. Equations for motion of fluid without incomressible assumtions Recall from week notes, the equations for conservation of mass and momentum, derived generally without any incomressibility
More informationChapter 1 Fundamentals
Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors
More informationI have not proofread these notes; so please watch out for typos, anything misleading or just plain wrong.
hermodynamics I have not roofread these notes; so lease watch out for tyos, anything misleading or just lain wrong. Please read ages 227 246 in Chater 8 of Kittel and Kroemer and ay attention to the first
More informationRigorous bounds on scaling laws in fluid dynamics
Rigorous bounds on scaling laws in fluid dynamics Felix Otto Notes written by Steffen Pottel Contents Main Result on Thermal Convection Derivation of the Model. A Stability Criterion........................
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationHigh speed wind tunnels 2.0 Definition of high speed. 2.1 Types of high speed wind tunnels
Module Lectures 6 to 1 High Seed Wind Tunnels Keywords: Blow down wind tunnels, Indraft wind tunnels, suersonic wind tunnels, c-d nozzles, second throat diffuser, shocks, condensation in wind tunnels,
More information16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE
16. CHARACTERISTICS OF SHOCK-WAVE UNDER LORENTZ FORCE AND ENERGY EXCHANGE H. Yamasaki, M. Abe and Y. Okuno Graduate School at Nagatsuta, Tokyo Institute of Technology 459, Nagatsuta, Midori-ku, Yokohama,
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 informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationHEAT, WORK, AND THE FIRST LAW OF THERMODYNAMICS
HET, ORK, ND THE FIRST L OF THERMODYNMIS 8 EXERISES Section 8. The First Law of Thermodynamics 5. INTERPRET e identify the system as the water in the insulated container. The roblem involves calculating
More informationHomework #11. (Due December 5 at the beginning of the class.) Numerical Method Series #6: Engineering applications of Newton-Raphson Method to
Homework #11 (Due December 5 at the beginning of the class.) Numerical Method Series #6: Engineering alications of Newton-Rahson Method to solving systems of nonlinear equations Goals: 1. Understanding
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More informationrate~ If no additional source of holes were present, the excess
DIFFUSION OF CARRIERS Diffusion currents are resent in semiconductor devices which generate a satially non-uniform distribution of carriers. The most imortant examles are the -n junction and the biolar
More informationMaximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming
Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More informationThe Role of Momentum Interpolation Mechanism of the Roe. Scheme in the Shock Instability
The Role of Momentum Interolation Mechanism of the Roe Scheme in the Shock Instabilit Xue-song Li Ke Laborator for Thermal Science and Power Engineering of Ministr of Education, Deartment of Thermal Engineering,
More informationModule 4 : Lecture 1 COMPRESSIBLE FLOWS (Fundamental Aspects: Part - I)
Module 4 : Lecture COMPRESSIBLE FLOWS (Fundamental Asects: Part - I) Overview In general, the liquids and gases are the states of a matter that comes under the same category as fluids. The incomressible
More informationQuantitative estimates of propagation of chaos for stochastic systems with W 1, kernels
oname manuscrit o. will be inserted by the editor) Quantitative estimates of roagation of chaos for stochastic systems with W, kernels Pierre-Emmanuel Jabin Zhenfu Wang Received: date / Acceted: date Abstract
More information4. Score normalization technical details We now discuss the technical details of the score normalization method.
SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules
More informationSteady, 1-d, constant area, adiabatic flow with no external work but with friction Conserved quantities
School of Aerosace Engineering Stead, -d, constant area, adiabatic flow with no external work but with friction Conserved quantities since adiabatic, no work: h o constant since Aconst: mass fluxρvconstant
More informationPhase transition. Asaf Pe er Background
Phase transition Asaf Pe er 1 November 18, 2013 1. Background A hase is a region of sace, throughout which all hysical roerties (density, magnetization, etc.) of a material (or thermodynamic system) are
More informationLinear diophantine equations for discrete tomography
Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationConvex Optimization methods for Computing Channel Capacity
Convex Otimization methods for Comuting Channel Caacity Abhishek Sinha Laboratory for Information and Decision Systems (LIDS), MIT sinhaa@mit.edu May 15, 2014 We consider a classical comutational roblem
More informationCFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids
CFL Conditions for Runge-Kutta Discontinuous Galerkin Methods on Triangular Grids T. Toulorge a,, W. Desmet a a K.U. Leuven, Det. of Mechanical Engineering, Celestijnenlaan 3, B-31 Heverlee, Belgium Abstract
More informationANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE
ANALYTICAL MODEL FOR THE BYPASS ALE IN A LOOP HEAT PIPE Michel Seetjens & Camilo Rindt Laboratory for Energy Technology Mechanical Engineering Deartment Eindhoven University of Technology The Netherlands
More informationLiquid water static energy page 1/8
Liquid water static energy age 1/8 1) Thermodynamics It s a good idea to work with thermodynamic variables that are conserved under a known set of conditions, since they can act as assive tracers and rovide
More informationLocation of solutions for quasi-linear elliptic equations with general gradient dependence
Electronic Journal of Qualitative Theory of Differential Equations 217, No. 87, 1 1; htts://doi.org/1.14232/ejqtde.217.1.87 www.math.u-szeged.hu/ejqtde/ Location of solutions for quasi-linear ellitic equations
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationEfficiencies. Damian Vogt Course MJ2429. Nomenclature. Symbol Denotation Unit c Flow speed m/s c p. pressure c v. Specific heat at constant J/kgK
Turbomachinery Lecture Notes 1 7-9-1 Efficiencies Damian Vogt Course MJ49 Nomenclature Subscrits Symbol Denotation Unit c Flow seed m/s c Secific heat at constant J/kgK ressure c v Secific heat at constant
More informationarxiv: v1 [physics.data-an] 26 Oct 2012
Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationJournal of Computational and Applied Mathematics. Numerical modeling of unsteady flow in steam turbine stage
Journal of Comutational and Alied Mathematics 234 (2010) 2336 2341 Contents lists available at ScienceDirect Journal of Comutational and Alied Mathematics journal homeage: www.elsevier.com/locate/cam Numerical
More informationTheory of turbomachinery. Chapter 1
Theory of turbomachinery Chater Introduction: Basic Princiles Take your choice of those that can best aid your action. (Shakeseare, Coriolanus) Introduction Definition Turbomachinery describes machines
More informationDetermination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure
Paer received: 7.10.008 UDC 61.64 Paer acceted: 0.04.009 Determination of Pressure Losses in Hydraulic Pieline Systems by Considering Temerature and Pressure Vladimir Savi 1,* - Darko Kneževi - Darko Lovrec
More informationwhether a process will be spontaneous, it is necessary to know the entropy change in both the
93 Lecture 16 he entroy is a lovely function because it is all we need to know in order to redict whether a rocess will be sontaneous. However, it is often inconvenient to use, because to redict whether
More informationThe Numerical Simulation of Gas Turbine Inlet-Volute Flow Field
World Journal of Mechanics, 013, 3, 30-35 doi:10.436/wjm.013.3403 Published Online July 013 (htt://www.scir.org/journal/wjm) The Numerical Simulation of Gas Turbine Inlet-Volute Flow Field Tao Jiang 1,
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationCompressible Flow Introduction. Afshin J. Ghajar
36 Comressible Flow Afshin J. Ghajar Oklahoma State University 36. Introduction...36-36. he Mach Number and Flow Regimes...36-36.3 Ideal Gas Relations...36-36.4 Isentroic Flow Relations...36-4 36.5 Stagnation
More informationJohn Weatherwax. Analysis of Parallel Depth First Search Algorithms
Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel
More informationAn Introduction to Theories of Turbulence. James Glimm Stony Brook University
An Introduction to Theories of Turbulence James Glimm Stony Brook University Topics not included (recent papers/theses, open for discussion during this visit) 1. Turbulent combustion 2. Turbulent mixing
More informationGeneralized Coiflets: A New Family of Orthonormal Wavelets
Generalized Coiflets A New Family of Orthonormal Wavelets Dong Wei, Alan C Bovik, and Brian L Evans Laboratory for Image and Video Engineering Deartment of Electrical and Comuter Engineering The University
More informationSpeed of sound measurements in liquid Methane at cryogenic temperature and for pressure up to 10 MPa
LNGII - raining Day Delft, August 07 Seed of sound measurements in liquid Methane at cryogenic temerature and for ressure u to 0 MPa Simona Lago*, P. Alberto Giuliano Albo INRiM Istituto Nazionale di Ricerca
More informationLocal Discontinuous Galerkin Methods for the Khokhlov Zabolotskaya Kuznetzov Equation
J Sci Comut (7 73:593 66 DOI.7/s95-7-5-z Local Discontinuous Galerkin Methods for the Khokhlov Zabolotskaya Kuznetzov Equation Ching-Shan Chou Weizhou Sun Yulong Xing He Yang Received: 5 Setember 6 / Revised:
More informationStatics and dynamics: some elementary concepts
1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and
More informationLarge Scale Simulations of Turbulent Flows for Industrial Applications. Lakhdar Remaki BCAM- Basque Centre for Applied Mathematics
Large Scale Simulations o Turbulent Flows or Industrial Alications Lakhdar Remaki BCAM- Basque Centre or Alied Mathematics Outline Flow Motion Simulation Physical Model Finite-Volume Numerical Method CFD
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationOPTIMUM TRANSONIC WING DESIGN USING CONTROL THEORY
OPTIMUM TRANSONIC WING DESIGN USING CONTROL THEORY Thomas V. Jones Professor of Engineering, Deartment of Aeronautics and Astronautics Stanford University, Stanford, CA 9435-435 jamesonbaboon.stanford.edu
More informationPretest (Optional) Use as an additional pacing tool to guide instruction. August 21
Trimester 1 Pretest (Otional) Use as an additional acing tool to guide instruction. August 21 Beyond the Basic Facts In Trimester 1, Grade 8 focus on multilication. Daily Unit 1: Rational vs. Irrational
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationLecture 8, the outline
Lecture, the outline loose end: Debye theory of solids more remarks on the first order hase transition. Bose Einstein condensation as a first order hase transition 4He as Bose Einstein liquid Lecturer:
More informationEstimation of Separable Representations in Psychophysical Experiments
Estimation of Searable Reresentations in Psychohysical Exeriments Michele Bernasconi (mbernasconi@eco.uninsubria.it) Christine Choirat (cchoirat@eco.uninsubria.it) Raffaello Seri (rseri@eco.uninsubria.it)
More informationResearch Article Comparison of HPM and PEM for the Flow of a Non-newtonian Fluid between Heated Parallel Plates
Research Journal of Alied Sciences, Engineering and Technology 7(): 46-434, 4 DOI:.96/rjaset.7.793 ISSN: 4-7459; e-issn: 4-7467 4 Maxwell Scientific Publication Cor. Submitted: November, 3 Acceted: January
More informationClassical gas (molecules) Phonon gas Number fixed Population depends on frequency of mode and temperature: 1. For each particle. For an N-particle gas
Lecture 14: Thermal conductivity Review: honons as articles In chater 5, we have been considering quantized waves in solids to be articles and this becomes very imortant when we discuss thermal conductivity.
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More informationGT NUMERICAL STUDY OF A CASCADE UNSTEADY SEPARATION FLOW
Proceedings of ASME urbo Exo 2004 2004 Power for Land, Sea, and Air Vienna, Austria, June 4-7, 2004 G2004-5395 NUMERICAL SUDY OF A CASCADE UNSEADY SEPARAION FLOW Zongjun Hu and GeCheng Zha Deartment of
More informationREFINED STRAIN ENERGY OF THE SHELL
REFINED STRAIN ENERGY OF THE SHELL Ryszard A. Walentyński Deartment of Building Structures Theory, Silesian University of Technology, Gliwice, PL44-11, Poland ABSTRACT The aer rovides information on evaluation
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationTowards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK
Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationNotes on pressure coordinates Robert Lindsay Korty October 1, 2002
Notes on ressure coordinates Robert Lindsay Korty October 1, 2002 Obviously, it makes no difference whether the quasi-geostrohic equations are hrased in height coordinates (where x, y,, t are the indeendent
More informationInformation collection on a graph
Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements
More informationNumerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model
Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model Q. Sun, Y.J. Liu, R.H, Jin School of Aeronautics, Northwestern Polytechnical University, Xi an 710072,
More information#A47 INTEGERS 15 (2015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS
#A47 INTEGERS 15 (015) QUADRATIC DIOPHANTINE EQUATIONS WITH INFINITELY MANY SOLUTIONS IN POSITIVE INTEGERS Mihai Ciu Simion Stoilow Institute of Mathematics of the Romanian Academy, Research Unit No. 5,
More informationLECTURE 7 NOTES. x n. d x if. E [g(x n )] E [g(x)]
LECTURE 7 NOTES 1. Convergence of random variables. Before delving into the large samle roerties of the MLE, we review some concets from large samle theory. 1. Convergence in robability: x n x if, for
More informationAnalysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia
roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.
More informationEvaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models
Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models Ketan N. Patel, Igor L. Markov and John P. Hayes University of Michigan, Ann Arbor 48109-2122 {knatel,imarkov,jhayes}@eecs.umich.edu
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationFUGACITY. It is simply a measure of molar Gibbs energy of a real gas.
FUGACITY It is simly a measure of molar Gibbs energy of a real gas. Modifying the simle equation for the chemical otential of an ideal gas by introducing the concet of a fugacity (f). The fugacity is an
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 informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationNumerical Linear Algebra
Numerical Linear Algebra Numerous alications in statistics, articularly in the fitting of linear models. Notation and conventions: Elements of a matrix A are denoted by a ij, where i indexes the rows and
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationRecursive Estimation of the Preisach Density function for a Smart Actuator
Recursive Estimation of the Preisach Density function for a Smart Actuator Ram V. Iyer Deartment of Mathematics and Statistics, Texas Tech University, Lubbock, TX 7949-142. ABSTRACT The Preisach oerator
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 informationTHREE-DIMENSIONAL POROUS MEDIA MODEL OF A HORIZONTAL STEAM GENERATOR. T.J.H. Pättikangas, J. Niemi and V. Hovi
THREE-DIMENSIONAL POROUS MEDIA MODEL OF A HORIZONTAL STEAM GENERATOR T.J.H. Pättikangas, J. Niemi and V. Hovi VTT Technical Research Centre of Finland, P.O.B. 000, FI 0044 VTT, Finland T. Toila and T.
More informationSYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY
SYMPLECTIC STRUCTURES: AT THE INTERFACE OF ANALYSIS, GEOMETRY, AND TOPOLOGY FEDERICA PASQUOTTO 1. Descrition of the roosed research 1.1. Introduction. Symlectic structures made their first aearance in
More informationWave Drift Force in a Two-Layer Fluid of Finite Depth
Wave Drift Force in a Two-Layer Fluid of Finite Deth Masashi Kashiwagi Research Institute for Alied Mechanics, Kyushu University, Jaan Abstract Based on the momentum and energy conservation rinciles, a
More informationMATH 2710: NOTES FOR ANALYSIS
MATH 270: NOTES FOR ANALYSIS The main ideas we will learn from analysis center around the idea of a limit. Limits occurs in several settings. We will start with finite limits of sequences, then cover infinite
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationSome results of convex programming complexity
2012c12 $ Ê Æ Æ 116ò 14Ï Dec., 2012 Oerations Research Transactions Vol.16 No.4 Some results of convex rogramming comlexity LOU Ye 1,2 GAO Yuetian 1 Abstract Recently a number of aers were written that
More informationChapter 9 Practical cycles
Prof.. undararajan Chater 9 Practical cycles 9. Introduction In Chaters 7 and 8, it was shown that a reversible engine based on the Carnot cycle (two reversible isothermal heat transfers and two reversible
More informationPressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids
American Journal of Alied Sciences (9): 19-195, 5 ISSN 1546-939 5 Science Publications Pressure-sensitivity Effects on Toughness Measurements of Comact Tension Secimens for Strain-hardening Solids Abdulhamid
More informationCFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS
CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS M. BORGHI - M. MILANI Diartimento di Scienze dell Ingegneria Università degli Studi di Modena Via Cami, 213/b 41100 Modena E-mail: borghi@omero.dsi.unimo.it
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 informationCarbuncle Phenomena and Other Shock Anomalies in Three Dimensions
Carbuncle Phenomena and Other Shock Anomalies in Three Dimensions Keiichi Kitamura * and Eiji Shima Jaan Aerosace Exloration Agency (JAXA) Sagamihara Kanagawa 5-50 Jaan and Phili. oe University of ichigan
More information