Introduction to Riemann Solvers
|
|
- Loraine Howard
- 5 years ago
- Views:
Transcription
1 CO 5 BOLD WORKSHOP 2, 2 4 June Introduction to Riemann Sovers Oskar Steiner
2 . Three major advancements in the numerica treatment of the hydrodynamic equations Three major progresses in computationa fuid dynamics of the past 4 years incude: the conservative formuation of the computationa scheme in terms of finite voumes, the technique of approimate Riemann-sovers for the computation of numerica fues, the fu-imiter technique for maintaining stabiity and monotonicity of higher-order accurate scheme. A these techniques are part of the numerica scheme of CO 5 BOLD.
3 2. Conservation aws finite voumes Consider the continuity equation: ρ t + (ρu) =. () Integration over a finite voume, V, and time period, T, eads to the integra form of this equation: V ρ(t,)dv V ρ(,)dv = T V (ρu) ndsdt (2) Soutions to Eq. (2) are caed weak soutions to the partia differentia equation (). Additionay to the soutions of Eq. (), the set of soutions to Eq. (2) encompasses discontinuous soutions, because no derivatives appear in Eq. (2). Discontinuous soutions to the Euer equations represent shock fronts of the rea word.
4 Conservation aws finite voumes (cont.) Consider the mass conservation in a one-dimensiona tube: A ( ρv) m ( ρv) 2 m(t + t) = m(t) + ρv A t ρv 2 A t (3) ρ (t + t) = ρ (t) t ( ρv 2 ρv ) (4) in the imit of and t ρ t = (ρv) But Eq. (3) is identica to the integra form Z V ρ(t,)dv Z V ρ(,)dv = Z T I V (ρu) n ds dt
5 Conservation aws finite voumes (cont.) The conservative, finite voume formuation has three highy desirabe properties: Conserved quantities (mass, momentum, energy) remain accuratey conserved Discontinuous soutions are incude by soving the integra form of the partia differentia equation It fufis one of two requirements of the theorem of La and Wendroff (9) that says: The approimate soution that is computed with a consistent and conservative scheme converges to a weak soution of the conservation aw.
6 Conservation aws finite voumes (cont.) Euer s equation in one dimension is given by where q t + f(q) =, q n+ i q = ρ ρu E C A f(q) = = q n i + t [f i /2 f i+/2 ], ρu ρu 2 + p u(e + p) C A In 3-D we have q t + f(q) + g(q) y + h(q) z =, with q = ρ ρu ρv ρw E C A f(q) = ρu ρu 2 + p ρuv ρuw u(e + p) C A
7 3. Riemann sovers A conservative finite voume scheme is an eact representation of the integra form of the partia differentia form of the conservation aw. The probem consists in computing the correct fu function f(q), i.e., ρv in the case of the continuity equation. It turns out that these fues can be computed eacty.
8 Riemann sovers (cont.) Idea of S.K. Godunov (959): Piecewise constant reconstruction with discontinuities at ce interfaces q q q i i+
9 Riemann sovers (cont.) The shock-tube probem t ρ ρ ρ r p p p r v v = v r =
10 Riemann sovers (cont.) The shock-tube probem t t ρ ρ ρ ρ ρ* ρ* r ρr ρ r p p p p p* p r p r v v v * v = v r = v v r
11 Riemann sovers (cont.) The shock-tube probem t t ρ ρ ρ ρ ρ* ρ* r ρr ρ r t p p p p p* t q * q * r p r p r q q r v v v * v = v r = v v r
12 Riemann sovers (cont.) The shock-tube probem t t ρ ρ ρ ρ ρ* ρ* r ρr ρ r t p p p p p* t q * q * r p r p r q q r v v v * v = v r = v v r q q = q r q = q r * q = q * q = q rf q = q
13 3.. The Riemann sover of Harten, La, and van Leer (HLL) Consider the system of one-dimensiona conservation aws q t + f(q) =, q(, ) = ( q if <, q r if >. T s t u = u a s r = q * q * r u+a q q r T s T s r r The integra form in the contro voume [, r ] [, T] is given by: Z r q(, T)d = Z r q(, )d + Z T f(q(, t))dt Z T f(q( r, t))dt
14 The HLL sover (cont.) Z r Z r q(, T) d = q(, )d + Z T f(q(, t)) dt Z T f(q( r, t)) dt = r q r q + T(f f r ), f = f(q ), f r = f(q r ) Z r q(, T) d = Z Ts q(, T) d + Z Tsr Ts q(, T) d + Z r Ts r q(, T) d = Z Tsr Ts q(, T) d + (Ts )q + ( r Ts r )q r T(s r s ) Z Tsr Ts q(, T)d := q h = s rq r s q + f f r s r s
15 The HLL sover (cont.) T u s = u a s r = T s q q * t q * r q r T u+a s r r Appying the integra form to the contro voume [, ] [, T] we obtain: Z Ts q(, T)d = Ts q + T(f f ), where f is the fu f(q) aong the t-ais. Hence, f = f s q T Z Ts q(, T)d. Doing the same for the contro voume [, r ] [, T] eads to f r = f r s r q r T Z Tsr q(, T)d. It foows that f = f r.
16 The HLL sover (cont.) Harten, La, and van Leer put forward the foowing approimation: s q t sr q h q(, t) = q r 8 >< >: q if t s, q h if s t s r, q r if t s r. f h = f + s (q h q ) or f h = f r + s r (q h q r ) f h = s rf s f r + s s r (q r q ) s r s
17 The HLL sover (cont.) The corresponding interce fu for the approimate Godunov method is then given by: f h i+/2 = f if s, s r f s f r + s s r (q r q ) s r s if s s r, f r if s r that can be used in the epicit conservative formua q n+ i = q n i + t [f i /2 f i+/2 ].
18 3.2. The HLLC sover The HLLC scheme is a modification of the HLL scheme in which the missing contact and shear waves are restored. s q t q* s* q* q r s r q(, t) = 8 >< >: q if t s, q if s t s, q r if s t s r, q r if t s r. Integrating over appropriate contro voumes, or more directy, by appying the Rankine-Hugoniot Conditions across each wave, we obtain f = f + s (q q ), f r = f + s (q r q ), f r = f r + s r (q r q r ).
19 The HLLC sover (cont.) The corresponding interce fu for the approimate Godunov method is then given by: f hc i+/2 = f if s, f + s (q q ) if s s, f r + s r (q r q r ) if s s r, f r if s r that can be used in the epicit conservative formua q n+ i = q n i + t [f i /2 f i+/2 ].
20 The HLLC sover (cont.) The intermediate states q and q r can be derived from f = f + s (q q ), f r = f + s (q r q ), f r = f r + s r (q r q r ), u = u r = u, p = p r = p, v = v, vr = v r, w = w, wr = w r, s = u, q k = ρ k sk u k s k s 2 «4 s v k w k E k ρ k + (s u k )[s + p k ρ k (s k u k ) 3 5, k =, r
21 3.3. Wave-speed estimates In order to competey determine the numerica fues in the HLL Riemann sover we need estimates for the wave speeds s and s r, and, for the HLLC sover, s. Given a positive speed s +, a simpe choice woud consist in s = s +, s r = s +. It is interesting to note that if we set s + equa to the maima speed according to the CFL-condition, i.e., s + = t, we obtain La-Friederichs numerica fu f i+/2 = 2 (f f r ) 2 which brings us back to a cassica scheme. t (q r q ),
22 Wave-speed estimates (cont.) More ingenious choices are motivated by a Roe average of the eft and right states, e.g., s = ũ ã, s r = ũ + ã, where ũ = ρ u + ρ r u r ρ + ρ r, ã =» (γ )( H 2ũ2 ) 2, with the enthapy H = (E + p)/ρ approimated as H = ρ H + ρ r H r ρ + ρ r.
23 Wave-speed estimates (cont.) In a different approach we first suppose to have estimates for p and u for the pressure and veocity in the Star Region. Then we compute the foowing wave speeds: s = u a r, s = u, s r = u r + a r r r, where r k = 8 >< >: if p p k raaction head» + γ + 2 2γ (p ) if p > p k shock p k p and u can be found from a inearization of the Riemann probem yieding p = 2 (p + p r ) 2 (u r u ) ρā, u = 2 (u + u r ) 2 where ρ = 2 (ρ + ρ r ), ā = 2 (a + a r ). (p r p ) ρā,
24 3.4. The Riemann sover of Roe Consider again the Riemann probem q t + f(q) =, 8 < q(, ) = : q if <, q r if >, where for the -spit three-dimensiona Euer equation ρ ρu ρu ρu 2 + p q = ρv, f(q) = ρuv B C ρw ρuw E u(e + p) C A.
25 The Roe sover (cont.) Using the chain rue, the conservation aw may be written as q t + f(q) = q t + A(q)q =, A(q) = f q. Roe s approach consists in repacing the Jacobian matri A(q) by a constant Jacobian à = Ã(q,q r ) resuting in the Riemann probem for the inear system q t + Ãq =, 8 < q(, ) = : q if <, q r if >, which can be soved eacty.
26 The Roe sover (cont.) Once the matri Ã(q,q r ), its eigenvaues λ i (q,q r ), and corresponding right eigenvectors k (i) (q,q r ) are avaiabe, the difference between right and eft state can be epanded in terms of the eigenvectors: mx q = q r q = α i k (i), i= from which one finds the wave strengths α i (q,q r ). ρ u L L w v L L p L t ρ * L ρ * R w v L L p* u* w v R R ρ u w v p R R R R R q i+/2 () = q + X α i k (i), λ i q i+/2 () = q r X α i k (i), λ i f i+/2 = f + X λ i α i λi k (i), f i+/2 = f r X α i λi k (i). or or λ i
27 The Roe sover (cont.) The -direction Jacobian matri for the Euer equations, A(q), is A = 2 4 ˆγH u 2 a 2 (3 γ)u ˆγv ˆγw ˆγ uv v u uw w u 2 u[(γ 3)H a2 ] H ˆγu 2 ˆγuv ˆγuw γu where ˆγ = γ and a = p γp/ρ. The eigenvaues are λ = u a, λ 2 = λ 3 = λ 4 = u, λ 5 = u + a. 3 5,
28 The Roe sover (cont.) The corresponding right eigenvaues are k () = 2 4 u a v w H ua 3 5, k (2) = 2 4 u v w 2 V 2 3 5, k (3) = 2 4 v 3 5, k (4) = 2 4 w 3 5, k (5) = 2 4 u + a v w H + ua 3 5, where H = E + p ρ, E = 2 ρv 2 + ρe, V 2 = u 2 + v 2 + w 2.
29 The Roe sover (cont.) Roe requires the constant Jacobian matri à = Ã(q,q r ) to satisfy the agebraic properties of the Jacobian A(q), i.e., λ i = λ i (q,q r ) R i, Ã(q,q) = A(q), f(q r ) f(q ) = Ã(q r q ). These conditions may be fufied with the foowing Roe averagged quantities to be used in the formuae for λ i and k (i) shown on the brevious pages: ũ = ṽ = w = ρ u + ρ r u r ρ + ρ r, ρ v + ρ r v r ρ + ρ r, ρ w + ρ r w r ρ + ρ r, H = ρ H + ρ r H r ρ + ρ r, ã = (γ )[ H 2 Ṽ 2 ] 2, Ṽ 2 = ũ 2 + ṽ 2 + w 2.
30 4. Higher order accurate methods Piece-wise inear, oca reconstruction (variabe etrapoation method) as part of the MUSCL scheme (Monotone Upstream-centred Scheme for Conservation Laws) u () i u L i u () i u () i+ q i () = q n i + i i [, ], q L i = q n i 2 i, u R i q R i = q n i + 2 i, eads to the Generaized Riemann Probem q t + f(q) =, 8 < q(, ) = : q i () if <, q i+ () if >.
31 Higher order accurate methods (cont.) a) u () i i+/2 b) t u () i+ Initia data (a) and soution structure (b) of the generaized Riemann probem. In the MUSCL-Hancock method intermediate boundary etrapoated vaues q R i obtained by and q L i+ are q L i = q L i + t 2 [f(ql i ) f(q R i )], q R i = q R i + t 2 [f(ql i ) f(q R i )], which then form the piece-wise constant data for the conventiona Riemann probem q t + f(q) =, 8 < q(, ) = : q R i if <, q L i+ if >.
32 Higher order accurate methods (cont.) An important property of the genera scaar conservation aw is monotonicity: u t + f(u) = (5) If two initia data functions v () and u () for Eq. (5) satisfy v () u (), then the corresponding soutions v(, t) and u(, t) satisfy v(, t) u(, t) t >. Correspondingy, a monotone scheme has the foowing property: if v n i u n i i then v n+ i u n+ i i. The Theorem of Godunov states that: There are no monotone, inear schemes for the soution of Eq. (5) of second or higher order of accuracy. (A inear scheme has generay the form q n+ i = k R X k= k L b k q n i+k, k L, k R N + )
33 Higher order accurate methods (cont.) In order to circumvent Godunov s theorem, higher order non-inear schemes were invented. One way of doing so consists in finding a sope imiter ξ i such that i = ξ i i, i = 2 ( + ω)(qn i q n i ) + 2 ( + ω)(qn i+ q n i ), ω [,]. ξ i depends in a non-inear way on the ratio r, where r = qn i q n i q n i+ qn i 2 ξ( r ) TVD region for sope imiters. For negative r the TVD region is the singe ine ξ =, for positive r the TVD region corresponds to the pink region. 2 r SUPERBEE, van Leer, MINMOD, etc type of sope imiters are a subset of this region.
34 Tabe of content. Three major advancements in the numerica treatment of the hydrodynamic Equations 2. Conservation aws finite voumes 3. Riemann sovers 3.. The Riemann sover of Harten, La, and van Leer (HLL) 3.2. The HLLC sover 3.3. Wave-speed estimates 3.4. The Riemann sover of Roe 4. Higher order accurate schemes References
35 References Laney, C.B.: 998, Computationa Gasdynamics, Cambridge University Press, Cambridge LeVeque, R.J.: 22, Finite Voume Methods for Hyperboic Probems, Cambridge University Press, Cambridge LeVeque, R.J.: 992, Numerica Methods for Conservation Laws, Birkhäuser Verag, Base LeVeque, R.J., Mihaas, D., Dorfi, E.A., and Müer, E.: 998, Computationa Methods for Astrophysica Fuid Fow, O. Steiner & A. Gautschy (eds.), Springer-Verag, Berin Toro, E.F.: 999, Riemann Sovers and Numerica Methods for Fuid Dynamics, Springer-Verag, Berin
Copyright information to be inserted by the Publishers. Unsplitting BGK-type Schemes for the Shallow. Water Equations KUN XU
Copyright information to be inserted by the Pubishers Unspitting BGK-type Schemes for the Shaow Water Equations KUN XU Mathematics Department, Hong Kong University of Science and Technoogy, Cear Water
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 information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
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 informationDual Integral Equations and Singular Integral. Equations for Helmholtz Equation
Int.. Contemp. Math. Sciences, Vo. 4, 9, no. 34, 1695-1699 Dua Integra Equations and Singuar Integra Equations for Hemhotz Equation Naser A. Hoshan Department of Mathematics TafiaTechnica University P.O.
More informationInternational Engineering Research Journal
Special Edition PGCON-MECH-7 Development of high resolution methods for solving D Euler equation Ms.Dipti A. Bendale, Dr.Prof. Jayant H. Bhangale and Dr.Prof. Milind P. Ray ϯ Mechanical Department, SavitribaiPhule
More informationLecture Notes for Math 251: ODE and PDE. Lecture 34: 10.7 Wave Equation and Vibrations of an Elastic String
ecture Notes for Math 251: ODE and PDE. ecture 3: 1.7 Wave Equation and Vibrations of an Eastic String Shawn D. Ryan Spring 212 ast Time: We studied other Heat Equation probems with various other boundary
More informationHomework #04 Answers and Hints (MATH4052 Partial Differential Equations)
Homework #4 Answers and Hints (MATH452 Partia Differentia Equations) Probem 1 (Page 89, Q2) Consider a meta rod ( < x < ), insuated aong its sides but not at its ends, which is initiay at temperature =
More informationA Finite Volume Code for 1D Gas Dynamics
A Finite Volume Code for 1D Gas Dynamics Michael Lavell Department of Applied Mathematics and Statistics 1 Introduction A finite volume code is constructed to solve conservative systems, such as Euler
More informationAXISYMMETRIC SOLUTIONS OF A TWO-DIMENSIONAL NONLINEAR WAVE SYSTEM WITH A TWO-CONSTANT EQUATION OF STATE
Eectronic Journa of Differentia Equations, Vo. 07 (07), No. 56, pp. 8. ISSN: 07-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu AXISYMMETRIC SOLUTIONS OF A TWO-DIMENSIONAL NONLINEAR
More informationStrauss PDEs 2e: Section Exercise 2 Page 1 of 12. For problem (1), complete the calculation of the series in case j(t) = 0 and h(t) = e t.
Strauss PDEs e: Section 5.6 - Exercise Page 1 of 1 Exercise For probem (1, compete the cacuation of the series in case j(t = and h(t = e t. Soution With j(t = and h(t = e t, probem (1 on page 147 becomes
More informationSolution of Wave Equation by the Method of Separation of Variables Using the Foss Tools Maxima
Internationa Journa of Pure and Appied Mathematics Voume 117 No. 14 2017, 167-174 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-ine version) ur: http://www.ijpam.eu Specia Issue ijpam.eu Soution
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationCoupling of LWR and phase transition models at boundary
Couping of LW and phase transition modes at boundary Mauro Garaveo Dipartimento di Matematica e Appicazioni, Università di Miano Bicocca, via. Cozzi 53, 20125 Miano Itay. Benedetto Piccoi Department of
More informationA Method for Simulating Two-Phase Pipe Flow with Real Equations of State
A Method for Simuating Two-Phase Pipe Fow with Rea Equations of State M. Hammer, A. Morin SINTEF Energy Research, P.O. Box 4761 Suppen, NO-7465 Trondheim, Norway Abstract Common two-fuid modes for pipe
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More information1. Introduction. The motivating question of this paper was the design of numerical schemes for P.D.E. systems that can be written under the form
SIAM J. NUMER. ANAL. Vo. 44, No. 1, pp. 3 321 c 26 Society for Industria and Appied Mathematics NUMERICAL METHODS FOR NONCONSERVATIVE HYPERBOLIC SYSTEMS: A THEORETICAL FRAMEWORK. CARLOS PARÉS Abstract.
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationSection 6: Magnetostatics
agnetic fieds in matter Section 6: agnetostatics In the previous sections we assumed that the current density J is a known function of coordinates. In the presence of matter this is not aways true. The
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 2: The Moment Equations
.615, MHD Theory of Fusion ystems Prof. Freidberg Lecture : The Moment Equations Botzmann-Maxwe Equations 1. Reca that the genera couped Botzmann-Maxwe equations can be written as f q + v + E + v B f =
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationLecture Notes 4: Fourier Series and PDE s
Lecture Notes 4: Fourier Series and PDE s 1. Periodic Functions A function fx defined on R is caed a periodic function if there exists a number T > such that fx + T = fx, x R. 1.1 The smaest number T for
More informationCHAPTER XIII FLOW PAST FINITE BODIES
HAPTER XIII LOW PAST INITE BODIES. The formation of shock waves in supersonic fow past bodies Simpe arguments show that, in supersonic fow past an arbitrar bod, a shock wave must be formed in front of
More information12.2. Maxima and Minima. Introduction. Prerequisites. Learning Outcomes
Maima and Minima 1. Introduction In this Section we anayse curves in the oca neighbourhood of a stationary point and, from this anaysis, deduce necessary conditions satisfied by oca maima and oca minima.
More informationFourier Series. 10 (D3.9) Find the Cesàro sum of the series. 11 (D3.9) Let a and b be real numbers. Under what conditions does a series of the form
Exercises Fourier Anaysis MMG70, Autumn 007 The exercises are taken from: Version: Monday October, 007 DXY Section XY of H F Davis, Fourier Series and orthogona functions EÖ Some exercises from earier
More informationMULTI-PERIOD MODEL FOR PART FAMILY/MACHINE CELL FORMATION. Objectives included in the multi-period formulation
ationa Institute of Technoogy aicut Department of echanica Engineering ULTI-PERIOD ODEL FOR PART FAILY/AHIE ELL FORATIO Given a set of parts, processing requirements, and avaiabe resources The objective
More information1D Heat Propagation Problems
Chapter 1 1D Heat Propagation Probems If the ambient space of the heat conduction has ony one dimension, the Fourier equation reduces to the foowing for an homogeneous body cρ T t = T λ 2 + Q, 1.1) x2
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 7
Strauss PDEs 2e: Section 4.3 - Exercise 1 Page 1 of 7 Exercise 1 Find the eigenvaues graphicay for the boundary conditions X(0) = 0, X () + ax() = 0. Assume that a 0. Soution The aim here is to determine
More informationMA 201: Partial Differential Equations Lecture - 10
MA 201: Partia Differentia Equations Lecture - 10 Separation of Variabes, One dimensiona Wave Equation Initia Boundary Vaue Probem (IBVP) Reca: A physica probem governed by a PDE may contain both boundary
More informationC. Fourier Sine Series Overview
12 PHILIP D. LOEWEN C. Fourier Sine Series Overview Let some constant > be given. The symboic form of the FSS Eigenvaue probem combines an ordinary differentia equation (ODE) on the interva (, ) with a
More informationJackson 4.10 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell
Jackson 4.10 Homework Probem Soution Dr. Christopher S. Baird University of Massachusetts Lowe PROBLEM: Two concentric conducting spheres of inner and outer radii a and b, respectivey, carry charges ±.
More informationHomotopy Perturbation Method for Solving Partial Differential Equations of Fractional Order
Int Journa of Math Anaysis, Vo 6, 2012, no 49, 2431-2448 Homotopy Perturbation Method for Soving Partia Differentia Equations of Fractiona Order A A Hemeda Department of Mathematics, Facuty of Science
More informationMultiscale Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Problems UNCORRECTED PROOF
AQ Mutiscae Domain Decomposition Preconditioners for 2 Anisotropic High-Contrast Probems Yachin Efendiev, Juan Gavis, Raytcho Lazarov, Svetozar Margenov 2, 4 and Jun Ren 5 Department of Mathematics, TAMU,
More informationFFTs in Graphics and Vision. Spherical Convolution and Axial Symmetry Detection
FFTs in Graphics and Vision Spherica Convoution and Axia Symmetry Detection Outine Math Review Symmetry Genera Convoution Spherica Convoution Axia Symmetry Detection Math Review Symmetry: Given a unitary
More informationA CLASS OF COMPUTATIONALLY FAST FIRST ORDER FINITE VOLUME SOLVERS: PVM METHODS.
A CLASS OF COMPUTATIONALLY FAST FIRST ORDER FINITE VOLUME SOLVERS: PVM METHODS. M. J. CASTRO DíAZ AND E. D. FERNÁNDEZ-NIETO Abstract. In this work, we present a cass of fast first order finite voume sovers,
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 informationEstablishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
Georgia Southern University Digita Commons@Georgia Southern Mathematica Sciences Facuty Pubications Department of Mathematica Sciences 2009 Estabishment of Weak Conditions for Darboux- Goursat-Beudon Theorem
More informationReliability: Theory & Applications No.3, September 2006
REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate
More informationMethods for Ordinary Differential Equations. Jacob White
Introduction to Simuation - Lecture 12 for Ordinary Differentia Equations Jacob White Thanks to Deepak Ramaswamy, Jaime Peraire, Micha Rewienski, and Karen Veroy Outine Initia Vaue probem exampes Signa
More informationVolume 13, MAIN ARTICLES
Voume 13, 2009 1 MAIN ARTICLES THE BASIC BVPs OF THE THEORY OF ELASTIC BINARY MIXTURES FOR A HALF-PLANE WITH CURVILINEAR CUTS Bitsadze L. I. Vekua Institute of Appied Mathematics of Iv. Javakhishvii Tbiisi
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationHigh Accuracy Split-Step Finite Difference Method for Schrödinger-KdV Equations
Commun. Theor. Phys. 70 208 43 422 Vo. 70, No. 4, October, 208 High Accuracy Spit-Step Finite Difference Method for Schrödinger-KdV Equations Feng Liao 廖锋, and Lu-Ming Zhang 张鲁明 2 Schoo of Mathematics
More informationOn norm equivalence between the displacement and velocity vectors for free linear dynamical systems
Kohaupt Cogent Mathematics (5 : 95699 http://dxdoiorg/8/33835595699 COMPUTATIONAL SCIENCE RESEARCH ARTICLE On norm equivaence between the dispacement and veocity vectors for free inear dynamica systems
More informationHigher dimensional PDEs and multidimensional eigenvalue problems
Higher dimensiona PEs and mutidimensiona eigenvaue probems 1 Probems with three independent variabes Consider the prototypica equations u t = u (iffusion) u tt = u (W ave) u zz = u (Lapace) where u = u
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f(q ) = exp ( βv( q )) 1, which is obtained by summing over
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 informationIntroduction to Simulation - Lecture 14. Multistep Methods II. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simuation - Lecture 14 Mutistep Methods II Jacob White Thans to Deepa Ramaswamy, Micha Rewiensi, and Karen Veroy Outine Sma Timestep issues for Mutistep Methods Reminder about LTE minimization
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationarxiv:quant-ph/ v3 6 Jan 1995
arxiv:quant-ph/9501001v3 6 Jan 1995 Critique of proposed imit to space time measurement, based on Wigner s cocks and mirrors L. Diósi and B. Lukács KFKI Research Institute for Partice and Nucear Physics
More informationMultiphasic equations for a NWP system
Mutiphasic equations for a NWP system CNRM/GMME/Méso-NH 21 novembre 2005 1 / 37 Contents 1 2 3 4 5 6 Individua mass species equations Momentum equation Thermodynamics equations Mutiphasics systems in practice
More informationTheory and implementation behind: Universal surface creation - smallest unitcell
Teory and impementation beind: Universa surface creation - smaest unitce Bjare Brin Buus, Jaob Howat & Tomas Bigaard September 15, 218 1 Construction of surface sabs Te aim for tis part of te project is
More informationMATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES
MATH 172: MOTIVATION FOR FOURIER SERIES: SEPARATION OF VARIABLES Separation of variabes is a method to sove certain PDEs which have a warped product structure. First, on R n, a inear PDE of order m is
More informationLecture 6: Moderately Large Deflection Theory of Beams
Structura Mechanics 2.8 Lecture 6 Semester Yr Lecture 6: Moderatey Large Defection Theory of Beams 6.1 Genera Formuation Compare to the cassica theory of beams with infinitesima deformation, the moderatey
More informationRESEARCH HIGHLIGHTS. WAF: Weighted Average Flux Method
RESEARCH HIGHLIGHTS (Last update: 3 rd April 2013) Here I briefly describe my contributions to research on numerical methods for hyperbolic balance laws that, in my view, have made an impact in the scientific
More informationVTU-NPTEL-NMEICT Project
MODUE-X -CONTINUOUS SYSTEM : APPROXIMATE METHOD VIBRATION ENGINEERING 14 VTU-NPTE-NMEICT Project Progress Report The Project on Deveopment of Remaining Three Quadrants to NPTE Phase-I under grant in aid
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationResearch Article New Iterative Method: An Application for Solving Fractional Physical Differential Equations
Abstract and Appied Anaysis Voume 203, Artice ID 6700, 9 pages http://d.doi.org/0.55/203/6700 Research Artice New Iterative Method: An Appication for Soving Fractiona Physica Differentia Equations A. A.
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationarxiv:nlin/ v2 [nlin.cd] 30 Jan 2006
expansions in semicassica theories for systems with smooth potentias and discrete symmetries Hoger Cartarius, Jörg Main, and Günter Wunner arxiv:nin/0510051v [nin.cd] 30 Jan 006 1. Institut für Theoretische
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationCHAPTER 10 TRANSVERSE VIBRATIONS-VI: FINITE ELEMENT ANALYSIS OF ROTORS WITH GYROSCOPIC EFFECTS
CHAPTER TRANSVERSE VIBRATIONS-VI: FINITE ELEMENT ANALYSIS OF ROTORS WITH GYROSCOPIC EFFECTS Previous, groscopic effects on a rotor with a singe disc were discussed in great detai b using the quasi-static
More informationSEMINAR 2. PENDULUMS. V = mgl cos θ. (2) L = T V = 1 2 ml2 θ2 + mgl cos θ, (3) d dt ml2 θ2 + mgl sin θ = 0, (4) θ + g l
Probem 7. Simpe Penduum SEMINAR. PENDULUMS A simpe penduum means a mass m suspended by a string weightess rigid rod of ength so that it can swing in a pane. The y-axis is directed down, x-axis is directed
More informationV.B The Cluster Expansion
V.B The Custer Expansion For short range interactions, speciay with a hard core, it is much better to repace the expansion parameter V( q ) by f( q ) = exp ( βv( q )), which is obtained by summing over
More informationPhysicsAndMathsTutor.com
. Two points A and B ie on a smooth horizonta tabe with AB = a. One end of a ight eastic spring, of natura ength a and moduus of easticity mg, is attached to A. The other end of the spring is attached
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationThe PLUTO code for astrophysical gasdynamics
Mem. S.A.It. Suppl. Vol. 13, 67 c SAIt 009 Memorie della Supplementi The PLUTO code for astrophysical gasdynamics A. Mignone 1, 1 Dipartimento di Fisica Generale Amedeo Avogadro, Università degli Studi
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationIntroduction. Figure 1 W8LC Line Array, box and horn element. Highlighted section modelled.
imuation of the acoustic fied produced by cavities using the Boundary Eement Rayeigh Integra Method () and its appication to a horn oudspeaer. tephen Kirup East Lancashire Institute, Due treet, Bacburn,
More informationEffect of Oxygen Injection into Argon Induction Plasmas on Chemically Non-Equilibrium Conditions
Proceedings of 17th Internationa Symposium on Pasma Chemistry, Toronto, Canada, August 7-12, 25 Effect of Oxygen Injection into Argon Induction Pasmas on Chemicay Non-Equiibrium Conditions Nobuhiko Atsuchi
More informationarxiv: v4 [math.na] 25 Aug 2014
BIT manuscript No. (wi be inserted by the editor) A muti-eve spectra deferred correction method Robert Speck Danie Ruprecht Matthew Emmett Michae Minion Matthias Boten Rof Krause arxiv:1307.1312v4 [math.na]
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationLecture Note 3: Stationary Iterative Methods
MATH 5330: Computationa Methods of Linear Agebra Lecture Note 3: Stationary Iterative Methods Xianyi Zeng Department of Mathematica Sciences, UTEP Stationary Iterative Methods The Gaussian eimination (or
More informationTHE ROLE OF ENERGY IMBALANCE MANAGEMENT ON POWER MARKET STABILITY
Proceedings of HICSS-31, Big Isand of Hawaii, January 6-9, 1998, Voume III, pp. 4-8. THE ROLE OF ENERGY IMBALANCE MANAGEMENT ON POWER MARKET STABILITY Fernando L. Avarado Department of Eectrica and C.
More information14-6 The Equation of Continuity
14-6 The Equation of Continuity 14-6 The Equation of Continuity Motion of rea fuids is compicated and poory understood (e.g., turbuence) We discuss motion of an idea fuid 1. Steady fow: Laminar fow, the
More informationIntuitionistic Fuzzy Optimization Technique for Nash Equilibrium Solution of Multi-objective Bi-Matrix Games
Journa of Uncertain Systems Vo.5, No.4, pp.27-285, 20 Onine at: www.jus.org.u Intuitionistic Fuzzy Optimization Technique for Nash Equiibrium Soution of Muti-objective Bi-Matri Games Prasun Kumar Naya,,
More informationInverse scattering transform for the integrable nonlocal nonlinear Schrödinger equation
Home Search Coections ournas About Contact us My IOPscience Inverse scattering transform for the integrabe nonoca noninear Schrödinger equation This content has been downoaded from IOPscience Pease scro
More informationMath-Net.Ru All Russian mathematical portal
Math-Net.Ru A Russian mathematica porta D. Zaora, On properties of root eements in the probem on sma motions of viscous reaxing fuid, Zh. Mat. Fiz. Ana. Geom., 217, Voume 13, Number 4, 42 413 DOI: https://doi.org/1.1547/mag13.4.42
More informationOn a geometrical approach in contact mechanics
Institut für Mechanik On a geometrica approach in contact mechanics Aexander Konyukhov, Kar Schweizerhof Universität Karsruhe, Institut für Mechanik Institut für Mechanik Kaiserstr. 12, Geb. 20.30 76128
More informationMultigrid Method for Elliptic Control Problems
J OHANNES KEPLER UNIVERSITÄT LINZ Netzwerk f ür Forschung, L ehre und Praxis Mutigrid Method for Eiptic Contro Probems MASTERARBEIT zur Erangung des akademischen Grades MASTER OF SCIENCE in der Studienrichtung
More informationIntroduction to LMTO method
1 Introduction to MTO method 24 February 2011; V172 P.Ravindran, FME-course on Ab initio Modeing of soar ce Materias 24 February 2011 Introduction to MTO method Ab initio Eectronic Structure Cacuations
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationMA 201: Partial Differential Equations Lecture - 11
MA 201: Partia Differentia Equations Lecture - 11 Heat Equation Heat conduction in a thin rod The IBVP under consideration consists of: The governing equation: u t = αu xx, (1) where α is the therma diffusivity.
More informationSECTION A. Question 1
SECTION A Question 1 (a) In the usua notation derive the governing differentia equation of motion in free vibration for the singe degree of freedom system shown in Figure Q1(a) by using Newton's second
More informationAerodynamic Design of Turbine Blades Using an Adjoint Equation Method
43rd AIAA Aerospace Sciences Meeting and Exhibit, January 10-13, 2005, Reno, Nevada Aerodynamic Design of Turbine ades Using an Adjoint Equation Method Hsiao-Yuan Wu, and Feng Liu Department of Mechanica
More informationA MULTI-ORDER DISCONTINUOUS GALERKIN MONTE CARLO METHOD FOR HYPERBOLIC PROBLEMS WITH STOCHASTIC PARAMETERS
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 A MULTI-ORDER DISCONTINUOUS GALERKIN MONTE CARLO METHOD FOR HYPERBOLIC PROBLEMS WITH
More informationPublished in: Proceedings of the Twenty Second Nordic Seminar on Computational Mechanics
Aaborg Universitet An Efficient Formuation of the Easto-pastic Constitutive Matrix on Yied Surface Corners Causen, Johan Christian; Andersen, Lars Vabbersgaard; Damkide, Lars Pubished in: Proceedings of
More informationMode in Output Participation Factors for Linear Systems
2010 American ontro onference Marriott Waterfront, Batimore, MD, USA June 30-Juy 02, 2010 WeB05.5 Mode in Output Participation Factors for Linear Systems Li Sheng, yad H. Abed, Munther A. Hassouneh, Huizhong
More informationRadar/ESM Tracking of Constant Velocity Target : Comparison of Batch (MLE) and EKF Performance
adar/ racing of Constant Veocity arget : Comparison of Batch (LE) and EKF Performance I. Leibowicz homson-csf Deteis/IISA La cef de Saint-Pierre 1 Bd Jean ouin 7885 Eancourt Cede France Isabee.Leibowicz
More informationIdentification of macro and micro parameters in solidification model
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vo. 55, No. 1, 27 Identification of macro and micro parameters in soidification mode B. MOCHNACKI 1 and E. MAJCHRZAK 2,1 1 Czestochowa University
More informationBohr s atomic model. 1 Ze 2 = mv2. n 2 Z
Bohr s atomic mode Another interesting success of the so-caed od quantum theory is expaining atomic spectra of hydrogen and hydrogen-ike atoms. The eectromagnetic radiation emitted by free atoms is concentrated
More informationCluster modelling. Collisions. Stellar Dynamics & Structure of Galaxies handout #2. Just as a self-gravitating collection of objects.
Stear Dynamics & Structure of Gaaxies handout # Custer modeing Just as a sef-gravitating coection of objects. Coisions Do we have to worry about coisions? Gobuar custers ook densest, so obtain a rough
More informationASummaryofGaussianProcesses Coryn A.L. Bailer-Jones
ASummaryofGaussianProcesses Coryn A.L. Baier-Jones Cavendish Laboratory University of Cambridge caj@mrao.cam.ac.uk Introduction A genera prediction probem can be posed as foows. We consider that the variabe
More informationHILBERT? What is HILBERT? Matlab Implementation of Adaptive 2D BEM. Dirk Praetorius. Features of HILBERT
Söerhaus-Workshop 2009 October 16, 2009 What is HILBERT? HILBERT Matab Impementation of Adaptive 2D BEM joint work with M. Aurada, M. Ebner, S. Ferraz-Leite, P. Godenits, M. Karkuik, M. Mayr Hibert Is
More informationA nodal collocation approximation for the multidimensional P L equations. 3D applications.
XXI Congreso de Ecuaciones Diferenciaes y Apicaciones XI Congreso de Matemática Apicada Ciudad Rea, 1-5 septiembre 9 (pp. 1 8) A noda coocation approximation for the mutidimensiona P L equations. 3D appications.
More informationTracking Control of Multiple Mobile Robots
Proceedings of the 2001 IEEE Internationa Conference on Robotics & Automation Seou, Korea May 21-26, 2001 Tracking Contro of Mutipe Mobie Robots A Case Study of Inter-Robot Coision-Free Probem Jurachart
More informationAn approximate method for solving the inverse scattering problem with fixed-energy data
J. Inv. I-Posed Probems, Vo. 7, No. 6, pp. 561 571 (1999) c VSP 1999 An approximate method for soving the inverse scattering probem with fixed-energy data A. G. Ramm and W. Scheid Received May 12, 1999
More information