Downloaded 01/04/14 to Redistribution subject to SIAM license or copyright; see
|
|
- Dominic Short
- 6 years ago
- Views:
Transcription
1 SIAM J. NUMER. ANAL. Vol. 46 No. pp c 008 Society for Indstrial and Applied Mathematics Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see ON STABILITY OF STAGGERED SCHEMES AMY L. BAUER RAPHAËL LOUBÈRE AND BURTON WENDROFF Abstract. This paper investigates the theoretical stability bond of a Lagrangian staggered scheme sed to solve hydrodynamics eqations. We present the two-dimensional D wave eqation as a possible model for this stdy and by sing the nmerical radis of the amplification matrix we prove that the family of schemes defined with two time-centering parameters is limited by a nonclassical stability bond limit defined with an analytical crve. We frther show that D nmerical experiments agree with this theoretical reslt. Key words. staggered schemes stability wave eqations Eler eqations AMS sbject classifications. 65N1 65N 65M1 65J10 65J15 DOI / Introdction. The concept of stability was first introdced in the seminal paper [1] of Corant Friedrichs and Lewy in 198 where they discssed finite-difference methods for solving partial differential eqations. Amazingly in the very same paper they stdied one of the first staggered nmerical schemes for the one-dimensional 1D wave eqation v t = cw x w t = cv x which was frther developed in 1967 by Richtmyer and Morton in [6] and cast in its modern form as: v j v n j Δt = c wn j+1/ wn j 1/ Δx w j 1/ wn j 1/ Δt = c v j v j 1. Δx By considering these works as the genesis for nmerical analysis and comptational flid dynamics one wold expect the stability of staggered nmerical schemes to be at present well nderstood. However this is far from being the case. Althogh staggered schemes have since been widely employed comptationally not mch is known analytically abot the stability of staggered schemes. The se of staggered schemes also dates back to the inception of Los Alamos National Laboratory where the calclation of certain time-dependent flid flows played an important part in the wartime work of the laboratory preface to the first edition of [6]. Staggered schemes have freqently been sed to solve the compressible flid dynamics eqations in their Lagrangian formlation. In this case a staggered spatial placement of variables is sed where the position and velocity are defined at grid points and density internal energy and the pressre are defined at cell centers. Since this time other Lagrangian hydrodynamics staggered nmerical codes have been de- Received by the editors May ; accepted for pblication in revised form October ; pblished electronically March This work was performed nder the aspices of the U.S. Department of Energy at Los Alamos National Laboratory nder contract W-7405-ENG-36. The athors acknowledge the partial spport of the DOE/ASCR Program in the Applied Mathematical Sciences and the Laboratory Directed Research and Development program LDRD. The athors also acknowledge the partial spport of DOE s Accelerated Strategic Compting Initiative ASCI. Theoretical Division T-14 Los Alamos National Laboratory MS-B6 Los Alamos NM albaer@lanl.gov. Institt de Mathematiqes et CNRS UMR-519 Université Pal-Sabatier Tolose Tolose Cedex 9 France lobere@mip.ps-tlse.fr. Theoretical Division T-7 Los Alamos National Laboratory MS-B84 Los Alamos NM bbw@lanl.gov. 996
2 ON STABILITY OF STAGGERED SCHEMES 997 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see veloped arond the globe bt again little advancement has been made towards a deeper nderstanding and analytical proof of the stability of sch schemes. Some years ago Shashkov proposed a simple problem for the prpose of testing the stability of Lagrangian hydrodynamics codes. The initial data are given as zero velocity constant pressre and constant density. One then comptes the total kinetic energy which shold be zero for all time. Caramana recently tried this problem with his two-dimensional D staggered grid Lagrangian code and discovered that it was nstable for certain mesh ratios that folklore indicated shold be otherwise the instability manifested itself as an explosive growth of kinetic energy. The srprise here is not that a 1D symmetry-preserving D code cold be nstable for 1D data; the srprise lies in the contradiction of the folklore. The precise nmerical scheme we are referring to and which is the focs of or investigation is the staggered-grid predictor-corrector compatible Lagrangian hydrodynamics scheme. We refer the reader to [3] [4] for the details of this scheme which we omit here. Briefly there are two parameters α and β sch that if α = β =1/ the scheme resembles a time-centered Crank Nicolson scheme while if α = β = 1 it looks like backward Eler. For the Corant Friedrichs Lewy CFL condition CFL = 1 the scheme is stable for α = β =1/. However for α =1β=1/ the scheme is nstable for CFL > 0.71 indicating that as α increases the CFL limit decreases. Baer et al. then created a 1D problem by sing the same data; they seed their 1D staggered-grid predictor-corrector compatible hydrodynamics Lagrangian code with a small random pertrbation of the pressre. Moreover they created the mltidimension version of the problem and ran their D and Caramana s 3D code as well. They observed the same phenomenon in 1D D and 3D [3]. The stability limit was conjectred to be CFL = 1/ αβ and nmerical tests sing the fll nonlinear eqations show that this is tre provided α 1/ and β 1/ [3]. We prove the above conjectre in D which contains 1D as a special case. Of corse there is no possibility of doing this for the fll nonlinear problem the Eler eqations. Ultimately one applies a von Nemann analysis to the linearized system. Baer et al. [3] almost scceeded in doing this theoretical analysis in 1D resorting to nmerical sampling only at the very end of the proof. We have scceeded in D by sing the nmerical radis of the amplification matrix as a tool an idea apparently first applied in []. This proof is the goal of this paper which is articlated as follows: First we present the Lagrangian coordinates and jstify the se of the wave eqation as a possible linear model for the Eler eqations in 1D and D; second we introdce the notation; third we show that the staggered implicit scheme is nconditionally stable for α 1/ β 1/; forth we show that the predictorcorrector staggered scheme is stable for α 1/ β 1/ and CFL 1/ αβ. In a paragraph we show on a specific 1D example that the schemes for α<1/ or β<1/ cannot be stable as some Forier components are amplified. Finally we show D nmerical reslts for the wave and Eler eqations sing a compatible Lagrangian hydrodynamics code.. Lagrangian coordinates and wave eqation models..1. 1D. The hydrodynamic eqations in 1D Lagrangian coordinates expressed in mass variables m are d dt = p m dτ dt = m de = p dt m
3 998 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see where is the velocity τ is the specific volme e is the specific internal energy and p = pτe is the pressre. Total energy E and entropy S are conserved: and de dt = 1 d dt + de dt = p m p m = p m T ds dt = de dt + pdτ dt =0. Therefore for smooth flows the exact differential eqations for p are and taking the coefficients p dp dτ d dt = p m dp dt = dp dτ dτ dt + dp de dp de dt = dτ pdp de m and dp de constant coefficients. Now assme that dp dτ to be constant yields a linear system with dp < 0 and de > 0. As dp c = dτ + pdp de is the sond speed the linearization of the system with constant coefficients is obtained by freezing the sond speed; that is c is assmed to be constant in space and time. Conseqently the reslting system is jst the wave eqation d dt = p m dp = c dt m. By sing the new variables = /c and dm/c = da and dropping the prime notation this linearized system can be written as.1 d dt = p a. dp dt = a... D. The same linearized redction can be performed for Lagrangian coordinates in D bt it is not qite as straightforward. In D the Eler eqation system
4 ON STABILITY OF STAGGERED SCHEMES 999 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see is ρ d dt = p x ρ dv dt = p y ρ dτ dt = x + v y ρ de dt = p x + v y dx dt = dy dt = v. By following [7] [8] we make the change of variables: with the choice x = xa b t y = ya b t a = xa b 0 b = ya b 0. The Jacobian matrix of the transformation x y a b is given by x x a J = y a y the determinant of which is the Jacobian J = x y a y x a. For the two velocity eqations we now have Jρ d dt = y p a y p a Jρ dv dt = x p a x p a. Bt then by starting the linearization by freezing the entries of the Jacobian we obtain Jρ d x dt a + v y = J p a a Jρ d x dt + v y = J p. Ths by letting = x a + v y a v = x + v y
5 1000 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see we obtain We next need to evalate x + v y ρ d dt = p a ρ dv dt = p. in the energy eqation. First J = y y a v Jv = x x a v. Then taking the derivative of and v yields J x = y y a y a J v y = x x a x a Smming the previos eqations gives x J x + J v y = + x x a + y a y y a x a y a y v a + v a y a x v a x a v v. x x a + y y a x y + a a and in order to cancel the cross terms we are forced to assme here that.3 x x a + y y a =0 so that x + v y = d 1 a + d v x d 1 = 1 y J + x d = 1 y J +. a a So ρ d dt = p a ρ dv dt = p ρ dτ dt = d 1 a + d v ρ de d dt = p 1 a + d v. v
6 ON STABILITY OF STAGGERED SCHEMES 1001 Now ρ can be scaled into a b and as before a simple sbstittion yields Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see Finally retrning to v a b set d dt = p a dv dt = p dp dt = d c 1 a + d v. = c d 1 v = c d v a = c d 1 ã b = c d b to obtain the D wave eqation as a model:.4 d dt = p ã dv dt = p b dp dt = ã + v b. However this formlation critically depends on.3 so ltimately we present this as only a potential model Notation. This section defines the notation that will be sed throghot the rest of this paper. The standard inner prodct on a complex vector space of dimension N is f g = N i=1 f i ḡ i for two complex vectors f =f i i=1...n and g =g i i=1...n. Time. We will assme eqal time steps; that is the temporal interval [0T > 0] is discretized into eqal intervals [t n t ] with t = t n +Δt. Space. In D we define niform rectangles with vertices x y where x i+1j = x +Δx and y +1 = y +Δy. Nodal qantities are indexed by i j while cell qantities are indexed by i + 1 j+ 1. Discretization. The ratio of time step to space step in each dimension is given by λ x = Δt Δx λ y = Δt Δy. Hence any qantity P at some point i j at time t n is represented as P n. Similarly any cell-based qantity C at time t is represented as C. i+ 1 j+ 1 For any variable w defined at two time levels t >t n on a point or in a cell we define its interpolated vale at an intermediate time n + κ as: w n+κ = κw +1 κ w n 0 κ 1. We frther define a vector w = v p t and w αβ = n+α v n+α p n+β t where { } = { : i j} v = {v : i j} and p = p i+ 1 j+ 1 : i j. 1 The Jacobian matrix J contains information relating to the volme shape and orientation of an element after the transformation a b x y. Λ = J T J is the associated metric tensor: x λ11 λ Λ = 1 = a + x x x a + y y a λ 1 λ x x a + y y y a a + y Λ is a symmetric matrix and moreover the constraint.3 states that λ 1 = 0: The metric tensor Λ is diagonal meaning that the implied system of coordinates is orthogonal.
7 100 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see Bondary conditions. We assme throghot that all sms exist and that Forier transforms can be taken. Also in the sal way it is spposed that effectively there are no bondaries and therefore that indices can be shifted in a sm withot changing its vale. For example for all integers k p i+ 1 j = i+kj p i+ 1 +kj. i i Operators. There are two sets of data: nodal qantities and cell qantities. Consider the operator Q x that maps cell data to nodal data given by Q x p = 1 p i+ 1 j+ 1 + p i+ 1 j 1 p i 1 j+ 1 p i 1 j 1. By sing the sal scalar prodct in both nodal and cell spaces the adjoint Q x maps nodal data to cell data and is defined by Q x p = Q x p = 1 p i+ 1 j+ 1 + p i+ 1 j 1 p i 1 j+ 1 p i 1 j 1 = 1 pi+ 1 j i+1j i+1j+1 ; that is Q x i+ 1 j+ = i+1j i+1j+1. Similar definitions exist for Q y and Q y. The notation for a staggered scheme on a rectanglar grid for the D wave eqation is depicted in Figre 3.1 where p is defined at the center of each cell and and v are defined at nodes. For this scheme p has fractional indices whereas and v have integer indices. Moreover one defines interpolated vales at each midedge point; for example for a midedge point defined with indices i + 1 j +1 one has p i+ 1 j+1 = 1 p i+ 1 j+ 3 + p i+ 1 j 1 and i+ 1 j+1 = i+1j+1 whereas for a midedge point i j + 1 one has p + 1 = 1 p i+ 1 j+ 1 + p i+ 3 j+ 1 and + 1 = i+1j. By sing this notation a flly implicit staggered scheme applied to system.4 can be written as 3.1 = n + λ x p n+α i+ 1 pn+α j i 1 j v = v n + λ y p n+α p n+α p = p n i+ 1 j+ 1 i+ 1 j+ + λ 1 x n+β n+β i+1j λ y The difference operator is given by the matrix 0 0 Q x M = 0 0 Q y Q x Q y 0 together with Λ = λ x λ y v n+β i+ 1 vn+β. j+1 i+ 1 j
8 ON STABILITY OF STAGGERED SCHEMES 1003 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see p i j 1/ i j 1/ v i j 1/ p p p p p i+1/ j+1 i+1/ j+1 v i+1/ j+1 i 1/ j+3/ i 1/ j+1/ i j+1 i j i+1/ j+3/ i+1 j+1 i+1 j p p p i 1/ j 1/ * p i+1/ j+1/ * * * i+1/ j 1/ i+1/ i i+1 i+3/ j+3/ p i+3/ j+1/ i+3/ j 1/ j+1 j+1/ j p i+1 j+1/ i+1 j+1/ v i+1 j+1/ Fig Rectanglar staggered scheme. p is defined at cell centers with fractional indices: p i+ 1 j+ 1 whereas v are defined at nodes with integer indices:. Moreover v p are interpolated at each midedge point. For example p i+ 1 j+1 = 1 p i+ 1 j+ 3 + p i+ 1 j 1 and i+ 1 j+1 = i+1j+1 and similarly for v. Ths the implicit difference scheme 3.1 has the form w = w n + ΛMΛw αβ. Theorem 3.1. The staggered implicit scheme is stable for any λ x λ y if α 1 and β 1. Proof. Throghot the proofs all sms are taken over both indices i and j. Applying the energy method yields Hα β = This follows from the fact that and therefore for real w n n+β + v v n v n+β + p p n i+ 1 j+ 1 i+ 1 j+ 1 M = M p n+α =0. i+ 1 j+ 1 Mw w =0.
9 1004 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Now let s define UβVβPα sch that Hα β =Uβ+V β+p α and Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see U β = = β V β = v = β v P α = p i+ 1 j+ 1 = α p i+ 1 j+ 1 n n+β 1 β n β 1 v n v n+β 1 β v n β 1 v p n+α i+ 1 j+ 1 p n i+ 1 j+ 1 1 α α 1 p i+ 1 j+ 1 p n i+ 1 j+ 1 p n i+ 1 j+ 1. n v n We detail the steps in the proof by sing U and comment that the steps in the proof for V and P are similarly obtained. Differentiating U with respect to β gives U β = + n n. Then by the Schwarz ineqality U β > 0 U = 0 only in the trivial case that the initial data are constant. Since 1 U = 1 n it follows that for β 1 while for β< 1 U β 1 U β < 1 n n. By applying the same reasoning to V and P we see that for β 1 and α 1 n + v v n + p i+ p n 1 j+ 1 i+ 1 j+ 0; 1 that is we have stability becase + v + p i+ 1 j+ 1 n + v n + p n i+ 1 j+ 1.
10 ON STABILITY OF STAGGERED SCHEMES 1005 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see 4. Rectanglar predictor-corrector scheme. The predictor-corrector scheme in D [4] for the wave eqation can now be described. Predictor step. ũ = n + λ x Q x p n ṽ = v n + λ y Q y p n p i+ 1 j+ = p n 1 i+ 1 j+ λ 1 x Q x n+β λ i+ 1 j+ 1 y Q y v n+β ; i+ 1 j+ 1 recall that w n+β = β w +1 βw n for w = v. Corrector step. = n + λ x Qx p n+α v = v n + λ y Qy p n+α p = p n i+ 1 j+ 1 i+ 1 j+ λ 1 x Q x n+β i+ 1 j+ λ 1 y Q y v n+β ; i+ 1 j+ 1 recall that p n+α = α p +1 αp n and now w n+β = βw +1 βw n for w = v. By sbstitting the eqation reslting from the predictor step into the corrector step we obtain the following difference scheme: w = Sw n. Rather than write this ot in terms of the difference operators we immediately move to the Forier transforms of the variables and operators. The Forier transform employs a sbstittion of variables; for instance for p we have p n i+ 1 j+ 1 p 0 e θnδt+iδi+ 1 Δx+γj+ 1 Δy where θ is complex and δ γ are real. After factoring we obtain: Qx p = p 0 e θnδt+iδiδx+γjδy 1 i sinδδx e iγδy + e iγδy = p 0 e θnδt+iδiδx+γjδy i sinδδx cosγδy with i = 1. Then if one denotes ξ = δδx and η = γδy the dimensionless wave nmbers in the x and y directions respectively on an niform mesh the Forier transforms of the operator Q x Q y are 4.1 Qx p = p 0 e θnδt+iiξ+jη i sin ξ cos η 4. Qy p = p 0 e θnδt+iiξ+jη i cos ξ sin η. The same formlas can be similarly obtained for Q. By setting Φ x =λ x sin ξ cos η and Φ y =λ y sin η cos ξ
11 1006 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see and dropping the hat notation we obtain 1 αφ x αφ xφ y iφ x 1 αβφ x +Φ y S = αφ xφ y 1 αφ y iφ y 1 αβ Φ x +Φ y iφ x 1 αβφ x +Φ y iφ y 1 αβφ x +Φ y 1+αβ Φ x +Φ y β. Φ x +Φ y The Lax Richtmyer stability theory [6] tells s that we mst show that the amplification factor S has niformly bonded powers. First observe that SS S Ssoitis not sfficient to show that the eigenvales of S are less than 1. We shall show instead that the nmerical radis of S is bonded by 1 since this implies that the powers of S are niformly bonded by. A srvey of reslts abot the nmerical radis is given in [5] along with a direct proof that if the nmerical radis is bonded by one the powers are bonded by. We remark that the Kreiss matrix theorem implies a weaker reslt namely that the powers are bonded by a constant that depends on the size of the matrix [6]. The nmerical radis was apparently first sed as a stability analysis tool in [] and here we se the basic ideas of that work. We henceforth se the same symbols for grid fnctions and their Forier transforms realizing that this is an abse of the notation bt hoping that no confsion will occr. The nmerical radis of S denoted RS is RS = sp Sw w with w w =1. w The matrix S can be split into real and imaginary parts as where Now let A = B = S = A +ib 1 αφ x αφ x Φ y 0 αφ x Φ y 1 αφ y αβ Φ x +Φ y β Φ x +Φ y 1 αβφ x +Φ y 0 0 Φ x 0 0 Φ y Φ x Φ y 0. r = Aw w j = Bw w. Since A and B are real and symmetric r and j are real. Then Sw w = Aw w +i Bw w = r +ij = r + j. Theorem 4.1. The D staggered rectanglar scheme is stable if α 1 β 1 and 4αβ maxλ xλ y 1. Proof. Sppose that α 1 β 1 and 4αβ maxλ xλ y 1. Let First note that x= αβ Φ x y= αβ Φ y. z = x + y 1
12 ON STABILITY OF STAGGERED SCHEMES 1007 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see since z = x + y = αβ Φ x +Φ y = αβ 4λ x sin ξ cos η +4λ y sin η cos ξ 4αβ max λ xλ [ y cos η sin ξ + cos ξ sin η ] }{{} 1 cos η sin ξ + cos ξ sin η 1. Recall that + v + p = 1 and note that 4.3 since 4.4 Also by setting x + yv 1 p x + y x + yv =x + yvx + yv =x + y vv + xyv + v x + y vv + x v + y a=1/α = + v x + y 1 p x + y. b=1/β we have that a b are positive and bonded by. Then r =1 b x + yv a x + y p + a x + y p j = ab R p x + yv 1 x + y where Rq is the real part of q and Now set r = 1 b x + yv az 1 z p j 4ab 1 z p x + yv. By γ 1. We ths have that Also note that for all 0 p 1 x + yv = z 1 p γ. r = 1 γbz 1 p az 1 z p j 4a γb z 1 z p 1 p. p 1 p 1 4. Since γb has the same range as b we can replace γb byb. What we need to show then is that g zpa b [ = 1 bz 1 p az 1 z p ] + abz 1 z 1 for 0 a b and 0 zp 1.
13 1008 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see Then Lemma 4.. The fnction gzpa b 1 for 0 a b and 0 zp 1. Proof: Let q = p t = z. g = [ 1 bt 1 q at 1 t q ] + abt 1 t. We first show that g 1 on the bondary of the t q domain. bondary vales: and since b1 q g t =0=1 g t =1 1. Consider the t For the q bondary vales: [ g q =0 1=bt +bt + a1 t ] [ ] bt +bt + a1 t [ ] bt +t +1 t =0 while g q =1 1=at 1 t [ ] [ ] +at 1 t+b1 t +at + b1 t 0. To finish the proof we need only to show that for each fixed t g a qadratic fnction in q does not have a maximm for 0 <q<1. This is achieved by showing that either g q q =1 0 in which case g is either monotone nondecreasing or has a minimm in 0 <q<1 or that g g q q =0 0 and q q =1 0 in which case g is monotone nonincreasing in 0 <q<1. We proceed as follows: g q =κ[ 1 bt + κq ] where Note that 4.5 κ = bt at + at. 1 bt + κ =1 at1 t 1 a Case a bα β: In this case κ 0 so by 4.5 g q q =1=κ 1 bt + κ κ 1 0. Case a bα β: Then κ 0ift t 0 =1 b a and κ 0ift t 0. If t t 0 k 0. Then as in the case above g q q =1> 0. If t t 0 k 0. Then by 4.5 both g g q q =1< 0 and q q =0< 0. Ths ends the proof of the lemma and completes the proof of the theorem. This case is the most important since for the compatible Lagrangian hydrodynamics scheme applied to the Eler eqations one mst necessarily have β = 1 for the scheme to conserve total energy.
14 ON STABILITY OF STAGGERED SCHEMES 1009 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see 5. Instability for α < 1/ or β < 1/. Both the flly implicit and the predictor-corrector schemes are nstable if either α<1/ orβ<1/ for any mesh ratios. We jst need to show this in 1D for the special case of p = for then the schemes are being applied to d dt = x with one set of difference eqations sed to advance at half-integer indexes i +1/ and a different set at integer indexes i. Instability occrs if some Forier component is amplified. So by setting Φ y = 0 in the amplification factor S for the predictor-corrector factor we need only to show first that 1 αφ x +iφ x 1 αβφ x > 1ifα<1/. This is clearly the case if Φ x is positive and sfficiently small since 1 αφ x +iφ x 1 αβφ x 1 =1+Φ x 1 α +Φ 4 x α β }{{} β +Φ x αβ >0 } {{ } oφ 4 x A similar argment works for the second row of S if β < 1/. The flly implicit scheme yields to the same kind of analysis; we omit the details. 6. Nmerical reslts. We solve the simple problem already described in [3] where on a nit domain sing a qadranglar mesh of velocity is set to zero density and sond speed are set to nity by sing a perfect eqation of state with γ the ratio of specific heats eqal to 5/3 and reflective bondary conditions. Since for this problem no velocity shold develop simlations are rn for a very large nmber of time cycles typically 10 5 and with varying vales for the CFL nmber λ = λ x = λ y and parameters α β. A sensitive gage of instability is to keep track of the total kinetic energy K λ t n = 1 [ n +v n ] for a given CFL nmber λ at a given time t n. Since the density and domain size are scaled to nity this nmber shold remain at the sqare of machine precision abot in or case. For an nstable scheme it is observed that K λ t n grows by several orders of magnitde long before the 10 5 cycle limit is reached. In this way one can accrately identify the stability bondary in the CFL nmber and α space for a given β. D wave eqations. We solve the D wave eqations with the rectanglar predictor-corrector scheme. From Theorem 4.1 we have that 4αβ max λ xλ 1 y 1 = λ αβ provided that λ = λ x = λ y. Ths the maximm CFL nmber achievable denoted CFL max therefore depends on α and β and is described by the following fnction: 6.1 CFL max = CFL max α β = 1 αβ. In Figre 6.1 we plot the maximal CFL as a fnction of α for a sample of five β vales: β = { }. The maximm CFL nmbers obtained nmerically with the code for a pair of α β are sperimposed; + for β =0.5 for β =0.6 for β =0.75 for β =0.9 and for β =1.0. As expected they perfectly match the theoretical CFL limit 6.1 given by Theorem 4.1. D Eler eqations. We present the reslts obtained with or D Lagrangian predictor-corrector code for the Eler eqations by sing Lagrangian coordinates. In this context see [3] for example the conservation of total energy is ensred provided that β = 1. In Figre 6. we plot seven data points symbol + obtained by the code which represent the maximm CFL nmber that can be sed for a given α for which.
15 1010 A. L. BAUER R. LOUBÈRE AND B. WENDROFF Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see CFL Exact β= max CFL = Fig D wave eqations maximm CFL nmber λ as a fnction of α β fixed for the kinetic energy to remain on the order of machine precision for the test case. The vertical dashed line is the α =1/ limit and the horizontal one is the CFL =1limit. The theorem predicts the continos thick lines and the code prodces the data: + for β =0.5 for β =0.6 for β =0.75 for β =0.9 and for β =1.0. Any scheme defined by a vale α 1/ and β 1/ is stable with the following CFL nmber: λ 1/ αβ. CFL Exact D Lagrang code α Fig. 6.. Hydrodynamics eqations maximm CFL nmber λ as a fnction of α β =1/ for the kinetic energy to remain on the order of machine precision for the test case. The theorem predicts the continos line; the D code prodces the dashed line. Any scheme defined by a vale α 1/ is stable with the CFL nmber λ 1/ α and is nstable otherwise. α
16 ON STABILITY OF STAGGERED SCHEMES 1011 Downloaded 01/04/14 to Redistribtion sbject to SIAM license or copyright; see nmerical stability is preserved. The exact fnction CFL dashed line predicted by 6.1 is sperimposed. 3 We can easily see that any scheme defined by a vale α 1/ is stable with the CFL nmber λ 1/ α and is nstable otherwise. We comment that for α =1/ one exactly reaches the maximm CFL nmber λ =1. The scheme with α = β = 1/ is the only staggered predictor-corrector scheme that reaches the maximm CFL nmber. Therefore the choice of sing this scheme in several Lagrangian simlation codes at least for this optimal stability reason seems jstified. The reslts of or simlations perfectly match the predicted crve; this is also tre of or 1D and 3D codes. 7. Conclsion. We have proposed the D wave eqation as a linear constant coefficient model for D Lagrangian hydrodynamics and we have established a sfficient stability condition for the standard staggered-grid nmerical scheme. The 1D case is concrrently obtained as a special case of the D analysis. This stability bondary has been tested with the D compatible Lagrangian hydrodynamics scheme for the wave eqations and the Eler eqations in Lagrangian coordinates. The nmerical reslts fit the theoretical crve showing that the general belief that for sch nmerical schemes the more implicit the more stable is not always tre the classical choice made by several generations of Lagrangian code developers to se the α = β =1/ scheme is reasonable as it leads to optimal stability reslts at least for the artificial test case presented. For the Eler eqations these nmerical reslts have been obtained in 1D D and 3D. Acknowledgment. The athors thank Misha Shashkov for many fritfl discssions. REFERENCES [1] R. Corant K. O. Friedrichs and H. Lewy Uber die partiellen Differenzengleichngen der mathematischen Physik Math. Ann pp [] P. D. Lax and B. Wendroff Difference schemes for hyperbolic eqations with high order of accracy Comm. Pre Appl. Math. XVII 1964 pp [3] A. L. Baer D. E. Brton E. J. Caramana R. Lobère M. J. Shashkov and P. P. Whalen The internal consistency stability and accracy of the discrete compatible formlation of Lagrangian hydrodynamics J. Compt. Phys pp [4] E. J. Caramana D. E. Brton M. J. Shashkov and P. P. Whalen The constrction of compatible hydrodynamics algorithms tilizing conservation of total energy J. Compt. Phys pp [5] M. Goldberg and E. Tadmor On the nmerical radis and its applications Linear Algebra Appl pp [6] R. D. Richtmyer and K. W. Morton Difference Methods for Initial Vale Problems nd ed. Interscience [7] W. H. Hi P. Y. Li and Z. W. Li A nified coordinate system for solving the two-dimensional Eler eqations J. Compt. Phys pp [8] W. H. Hi The nified coordinate system in comptational flid dynamics Commn. Compt. Phys. 007 pp In this case CFL = 1/ αβ =1/ α becase β =1/.
Discontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationA Survey of the Implementation of Numerical Schemes for Linear Advection Equation
Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear
More informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationL 1 -smoothing for the Ornstein-Uhlenbeck semigroup
L -smoothing for the Ornstein-Uhlenbeck semigrop K. Ball, F. Barthe, W. Bednorz, K. Oleszkiewicz and P. Wolff September, 00 Abstract Given a probability density, we estimate the rate of decay of the measre
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationChem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics. Fall Semester Homework Problem Set Number 10 Solutions
Chem 4501 Introdction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Nmber 10 Soltions 1. McQarrie and Simon, 10-4. Paraphrase: Apply Eler s theorem
More informationActive Flux Schemes for Advection Diffusion
AIAA Aviation - Jne, Dallas, TX nd AIAA Comptational Flid Dynamics Conference AIAA - Active Fl Schemes for Advection Diffsion Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3, USA Downloaded
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationVectors in Rn un. This definition of norm is an extension of the Pythagorean Theorem. Consider the vector u = (5, 8) in R 2
MATH 307 Vectors in Rn Dr. Neal, WKU Matrices of dimension 1 n can be thoght of as coordinates, or ectors, in n- dimensional space R n. We can perform special calclations on these ectors. In particlar,
More information08.06 Shooting Method for Ordinary Differential Equations
8.6 Shooting Method for Ordinary Differential Eqations After reading this chapter, yo shold be able to 1. learn the shooting method algorithm to solve bondary vale problems, and. apply shooting method
More informationρ u = u. (1) w z will become certain time, and at a certain point in space, the value of
THE CONDITIONS NECESSARY FOR DISCONTINUOUS MOTION IN GASES G I Taylor Proceedings of the Royal Society A vol LXXXIV (90) pp 37-377 The possibility of the propagation of a srface of discontinity in a gas
More informationConditions for Approaching the Origin without Intersecting the x-axis in the Liénard Plane
Filomat 3:2 (27), 376 377 https://doi.org/.2298/fil7276a Pblished by Faclty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Conditions for Approaching
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationOPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1
OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationMath 116 First Midterm October 14, 2009
Math 116 First Midterm October 14, 9 Name: EXAM SOLUTIONS Instrctor: Section: 1. Do not open this exam ntil yo are told to do so.. This exam has 1 pages inclding this cover. There are 9 problems. Note
More informationStudy of the diffusion operator by the SPH method
IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait
More informationClassify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.
Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports
More informationMAXIMUM AND ANTI-MAXIMUM PRINCIPLES FOR THE P-LAPLACIAN WITH A NONLINEAR BOUNDARY CONDITION. 1. Introduction. ν = λ u p 2 u.
2005-Ojda International Conference on Nonlinear Analysis. Electronic Jornal of Differential Eqations, Conference 14, 2006, pp. 95 107. ISSN: 1072-6691. URL: http://ejde.math.txstate.ed or http://ejde.math.nt.ed
More informationA Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations
Applied Mathematics, 05, 6, 04-4 Pblished Online November 05 in SciRes. http://www.scirp.org/jornal/am http://d.doi.org/0.46/am.05.685 A Comptational Stdy with Finite Element Method and Finite Difference
More informationMomentum Equation. Necessary because body is not made up of a fixed assembly of particles Its volume is the same however Imaginary
Momentm Eqation Interest in the momentm eqation: Qantification of proplsion rates esign strctres for power generation esign of pipeline systems to withstand forces at bends and other places where the flow
More informationThe Heat Equation and the Li-Yau Harnack Inequality
The Heat Eqation and the Li-Ya Harnack Ineqality Blake Hartley VIGRE Research Paper Abstract In this paper, we develop the necessary mathematics for nderstanding the Li-Ya Harnack ineqality. We begin with
More informationPhysicsAndMathsTutor.com
C Integration - By sbstittion PhysicsAndMathsTtor.com. Using the sbstittion cos +, or otherwise, show that e cos + sin d e(e ) (Total marks). (a) Using the sbstittion cos, or otherwise, find the eact vale
More informationAn Investigation into Estimating Type B Degrees of Freedom
An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information
More informationPrandl established a universal velocity profile for flow parallel to the bed given by
EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More informationB-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Quadratic Optimization Problems in Continuous and Binary Variables
B-469 Simplified Copositive and Lagrangian Relaxations for Linearly Constrained Qadratic Optimization Problems in Continos and Binary Variables Naohiko Arima, Snyong Kim and Masakaz Kojima October 2012,
More informationEfficiency Increase and Input Power Decrease of Converted Prototype Pump Performance
International Jornal of Flid Machinery and Systems DOI: http://dx.doi.org/10.593/ijfms.016.9.3.05 Vol. 9, No. 3, Jly-September 016 ISSN (Online): 188-9554 Original Paper Efficiency Increase and Inpt Power
More informationFREQUENCY DOMAIN FLUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS
7 TH INTERNATIONAL CONGRESS O THE AERONAUTICAL SCIENCES REQUENCY DOMAIN LUTTER SOLUTION TECHNIQUE USING COMPLEX MU-ANALYSIS Yingsong G, Zhichn Yang Northwestern Polytechnical University, Xi an, P. R. China,
More informationLewis number and curvature effects on sound generation by premixed flame annihilation
Center for Trblence Research Proceedings of the Smmer Program 2 28 Lewis nmber and crvatre effects on sond generation by premixed flame annihilation By M. Talei, M. J. Brear AND E. R. Hawkes A nmerical
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationOn the Optimization of Numerical Dispersion and Dissipation of Finite Difference Scheme for Linear Advection Equation
Applied Mathematical Sciences, Vol. 0, 206, no. 48, 238-2389 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ams.206.6463 On the Optimization of Nmerical Dispersion and Dissipation of Finite Difference
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationThe spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam
J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College
More informationm = Average Rate of Change (Secant Slope) Example:
Average Rate o Change Secant Slope Deinition: The average change secant slope o a nction over a particlar interval [a, b] or [a, ]. Eample: What is the average rate o change o the nction over the interval
More informationQUANTILE ESTIMATION IN SUCCESSIVE SAMPLING
Jornal of the Korean Statistical Society 2007, 36: 4, pp 543 556 QUANTILE ESTIMATION IN SUCCESSIVE SAMPLING Hosila P. Singh 1, Ritesh Tailor 2, Sarjinder Singh 3 and Jong-Min Kim 4 Abstract In sccessive
More informationApplying Laminar and Turbulent Flow and measuring Velocity Profile Using MATLAB
IOS Jornal of Mathematics (IOS-JM) e-issn: 78-578, p-issn: 319-765X. Volme 13, Isse 6 Ver. II (Nov. - Dec. 17), PP 5-59 www.iosrjornals.org Applying Laminar and Trblent Flow and measring Velocity Profile
More informationComputational Geosciences 2 (1998) 1, 23-36
A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), 23-36 Abstract.
More informationKrauskopf, B., Lee, CM., & Osinga, HM. (2008). Codimension-one tangency bifurcations of global Poincaré maps of four-dimensional vector fields.
Kraskopf, B, Lee,, & Osinga, H (28) odimension-one tangency bifrcations of global Poincaré maps of for-dimensional vector fields Early version, also known as pre-print Link to pblication record in Explore
More informationGravitational Instability of a Nonrotating Galaxy *
SLAC-PUB-536 October 25 Gravitational Instability of a Nonrotating Galaxy * Alexander W. Chao ;) Stanford Linear Accelerator Center Abstract Gravitational instability of the distribtion of stars in a galaxy
More informationElements of Coordinate System Transformations
B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationChapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS
Chapter 3 MATHEMATICAL MODELING OF DYNAMIC SYSTEMS 3. System Modeling Mathematical Modeling In designing control systems we mst be able to model engineered system dynamics. The model of a dynamic system
More informationChapter 2 Difficulties associated with corners
Chapter Difficlties associated with corners This chapter is aimed at resolving the problems revealed in Chapter, which are cased b corners and/or discontinos bondar conditions. The first section introdces
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationOptimal Control, Statistics and Path Planning
PERGAMON Mathematical and Compter Modelling 33 (21) 237 253 www.elsevier.nl/locate/mcm Optimal Control, Statistics and Path Planning C. F. Martin and Shan Sn Department of Mathematics and Statistics Texas
More information5. The Bernoulli Equation
5. The Bernolli Eqation [This material relates predominantly to modles ELP034, ELP035] 5. Work and Energy 5. Bernolli s Eqation 5.3 An example of the se of Bernolli s eqation 5.4 Pressre head, velocity
More informationA Decomposition Method for Volume Flux. and Average Velocity of Thin Film Flow. of a Third Grade Fluid Down an Inclined Plane
Adv. Theor. Appl. Mech., Vol. 1, 8, no. 1, 9 A Decomposition Method for Volme Flx and Average Velocit of Thin Film Flow of a Third Grade Flid Down an Inclined Plane A. Sadighi, D.D. Ganji,. Sabzehmeidani
More informationEXERCISES WAVE EQUATION. In Problems 1 and 2 solve the heat equation (1) subject to the given conditions. Assume a rod of length L.
.4 WAVE EQUATION 445 EXERCISES.3 In Problems and solve the heat eqation () sbject to the given conditions. Assme a rod of length.. (, t), (, t) (, ),, > >. (, t), (, t) (, ) ( ) 3. Find the temperatre
More informationSTEP Support Programme. STEP III Hyperbolic Functions: Solutions
STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh
More information1 The space of linear transformations from R n to R m :
Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations
More informationTIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS
International Jornal of Wavelets, Mltiresoltion and Information Processing Vol. 4, No. (26) 65 79 c World Scientific Pblishing Company TIME ACCURATE FAST THREE-STEP WAVELET-GALERKIN METHOD FOR PARTIAL
More informationEssentials of optimal control theory in ECON 4140
Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationCRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE. Jingbo Xia
CRITERIA FOR TOEPLITZ OPERATORS ON THE SPHERE Jingbo Xia Abstract. Let H 2 (S) be the Hardy space on the nit sphere S in C n. We show that a set of inner fnctions Λ is sfficient for the prpose of determining
More information4 Primitive Equations
4 Primitive Eqations 4.1 Spherical coordinates 4.1.1 Usefl identities We now introdce the special case of spherical coordinates: (,, r) (longitde, latitde, radial distance from Earth s center), with 0
More informationEVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE
13 th World Conference on Earthqake Engineering Vancover, B.C., Canada Agst 1-6, 2004 Paper No. 3099 EVALUATION OF GROUND STRAIN FROM IN SITU DYNAMIC RESPONSE Ellen M. RATHJE 1, Wen-Jong CHANG 2, Kenneth
More informationAdmissibility under the LINEX loss function in non-regular case. Hidekazu Tanaka. Received November 5, 2009; revised September 2, 2010
Scientiae Mathematicae Japonicae Online, e-2012, 427 434 427 Admissibility nder the LINEX loss fnction in non-reglar case Hidekaz Tanaka Received November 5, 2009; revised September 2, 2010 Abstract. In
More informationDetermining of temperature field in a L-shaped domain
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 0, (:-8 Determining of temperatre field in a L-shaped domain Oigo M. Zongo, Sié Kam, Kalifa Palm, and Alione Oedraogo
More informationWhen are Two Numerical Polynomials Relatively Prime?
J Symbolic Comptation (1998) 26, 677 689 Article No sy980234 When are Two Nmerical Polynomials Relatively Prime? BERNHARD BECKERMANN AND GEORGE LABAHN Laboratoire d Analyse Nmériqe et d Optimisation, Université
More informationRemarks on strongly convex stochastic processes
Aeqat. Math. 86 (01), 91 98 c The Athor(s) 01. This article is pblished with open access at Springerlink.com 0001-9054/1/010091-8 pblished online November 7, 01 DOI 10.1007/s00010-01-016-9 Aeqationes Mathematicae
More informationApproach to a Proof of the Riemann Hypothesis by the Second Mean-Value Theorem of Calculus
Advances in Pre Mathematics, 6, 6, 97- http://www.scirp.org/jornal/apm ISSN Online: 6-384 ISSN Print: 6-368 Approach to a Proof of the Riemann Hypothesis by the Second Mean-Vale Theorem of Calcls Alfred
More informationPhysicsAndMathsTutor.com
. Two smooth niform spheres S and T have eqal radii. The mass of S is 0. kg and the mass of T is 0.6 kg. The spheres are moving on a smooth horizontal plane and collide obliqely. Immediately before the
More informationInertial Instability of Arbitrarily Meandering Currents Governed by the Eccentrically Cyclogeostrophic Equation
Jornal of Oceanography, Vol. 59, pp. 163 to 17, 3 Inertial Instability of Arbitrarily Meandering Crrents Governed by the Eccentrically Cyclogeostrophic Eqation HIDEO KAWAI* 131-81 Shibagahara, Kse, Joyo,
More informationFormal Methods for Deriving Element Equations
Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact
More informationAffine Invariant Total Variation Models
Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,
More informationInternational Journal of Physical and Mathematical Sciences journal homepage:
64 International Jornal of Physical and Mathematical Sciences Vol 2, No 1 (2011) ISSN: 2010-1791 International Jornal of Physical and Mathematical Sciences jornal homepage: http://icoci.org/ijpms PRELIMINARY
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationBertrand s Theorem. October 8, µr 2 + V (r) 0 = dv eff dr. 3 + dv. f (r 0 )
Bertrand s Theorem October 8, Circlar orbits The eective potential, V e = has a minimm or maximm at r i and only i so we mst have = dv e L µr + V r = L µ 3 + dv = L µ 3 r r = L µ 3 At this radis, there
More informationOn Generalizations of Asymptotically AdS 3 Spaces and Geometry of SL(N)
EJTP 9, No. 7 01) 49 56 Electronic Jornal of Theoretical Physics On Generalizations of Asymptotically AdS 3 Spaces and Geometry of SLN) Heikki Arponen Espoontie 1 A, 0770 Espoo, Finland Received 30 October
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationThe Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks
The Cryptanalysis of a New Pblic-Key Cryptosystem based on Modlar Knapsacks Yeow Meng Chee Antoine Jox National Compter Systems DMI-GRECC Center for Information Technology 45 re d Ulm 73 Science Park Drive,
More information1 Differential Equations for Solid Mechanics
1 Differential Eqations for Solid Mechanics Simple problems involving homogeneos stress states have been considered so far, wherein the stress is the same throghot the component nder std. An eception to
More informationVariability sustained pattern formation in subexcitable media
Variability sstained pattern formation in sbexcitable media Erik Glatt, Martin Gassel, and Friedemann Kaiser Institte of Applied Physics, Darmstadt University of Technology, 64289 Darmstadt, Germany (Dated:
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
More informationFigure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study).
Wedge Copla This docment explains the constrction and properties o a particlar geometrical copla sed to it dependency data rom the edram case stdy done at Portland State University. The probability density
More informationarxiv: v1 [physics.flu-dyn] 4 Sep 2013
THE THREE-DIMENSIONAL JUMP CONDITIONS FOR THE STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY, SINGULAR FORCES, AND AN INCOMPRESSIBLE INTERFACE PRERNA GERA AND DAVID SALAC arxiv:1309.1728v1 physics.fl-dyn]
More informationIntrodction In this paper we develop a local discontinos Galerkin method for solving KdV type eqations containing third derivative terms in one and ml
A Local Discontinos Galerkin Method for KdV Type Eqations Je Yan and Chi-Wang Sh 3 Division of Applied Mathematics Brown University Providence, Rhode Island 09 ABSTRACT In this paper we develop a local
More informationNumerical Study on Bouncing and Separation Collision Between Two Droplets Considering the Collision-Induced Breakup
Jornal of Mechanical Science and Technology (007) 585~59 Jornal of Mechanical Science and Technology Nmerical Stdy on Boncing and Separation Collision Between Two Droplets Considering the Collision-Indced
More informationHome Range Formation in Wolves Due to Scent Marking
Blletin of Mathematical Biology () 64, 61 84 doi:1.16/blm.1.73 Available online at http://www.idealibrary.com on Home Range Formation in Wolves De to Scent Marking BRIAN K. BRISCOE, MARK A. LEWIS AND STEPHEN
More informationDownloaded 07/06/18 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 4, No., pp. A4 A7 c 8 Society for Indstrial and Applied Mathematics Downloaded 7/6/8 to 8.83.63.. Redistribtion sbject to SIAM license or copyright; see http://www.siam.org/jornals/ojsa.php
More informationDYNAMICAL LOWER BOUNDS FOR 1D DIRAC OPERATORS. 1. Introduction We consider discrete, resp. continuous, Dirac operators
DYNAMICAL LOWER BOUNDS FOR D DIRAC OPERATORS ROBERTO A. PRADO AND CÉSAR R. DE OLIVEIRA Abstract. Qantm dynamical lower bonds for continos and discrete one-dimensional Dirac operators are established in
More informationMECHANICS OF SOLIDS COMPRESSION MEMBERS TUTORIAL 2 INTERMEDIATE AND SHORT COMPRESSION MEMBERS
MECHANICS O SOIDS COMPRESSION MEMBERS TUTORIA INTERMEDIATE AND SHORT COMPRESSION MEMBERS Yo shold jdge yor progress by completing the self assessment exercises. On completion of this ttorial yo shold be
More information9. Tensor product and Hom
9. Tensor prodct and Hom Starting from two R-modles we can define two other R-modles, namely M R N and Hom R (M, N), that are very mch related. The defining properties of these modles are simple, bt those
More informationWe automate the bivariate change-of-variables technique for bivariate continuous random variables with
INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,
More informationReduction of over-determined systems of differential equations
Redction of oer-determined systems of differential eqations Maim Zaytse 1) 1, ) and Vyachesla Akkerman 1) Nclear Safety Institte, Rssian Academy of Sciences, Moscow, 115191 Rssia ) Department of Mechanical
More informationExperimental Study of an Impinging Round Jet
Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a
More informationConvex Hypersurfaces of Constant Curvature in Hyperbolic Space
Srvey in Geometric Analysis and Relativity ALM 0, pp. 41 57 c Higher Edcation Press and International Press Beijing-Boston Convex Hypersrfaces of Constant Crvatre in Hyperbolic Space Bo Gan and Joel Sprck
More information4.2 First-Order Logic
64 First-Order Logic and Type Theory The problem can be seen in the two qestionable rles In the existential introdction, the term a has not yet been introdced into the derivation and its se can therefore
More informationLecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018
Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:
More informationDecoder Error Probability of MRD Codes
Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehighed Zhiyan Yan Department of Electrical
More informationGraphs and Their. Applications (6) K.M. Koh* F.M. Dong and E.G. Tay. 17 The Number of Spanning Trees
Graphs and Their Applications (6) by K.M. Koh* Department of Mathematics National University of Singapore, Singapore 1 ~ 7543 F.M. Dong and E.G. Tay Mathematics and Mathematics EdOOation National Institte
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More informationDecision Oriented Bayesian Design of Experiments
Decision Oriented Bayesian Design of Experiments Farminder S. Anand*, Jay H. Lee**, Matthew J. Realff*** *School of Chemical & Biomoleclar Engineering Georgia Institte of echnology, Atlanta, GA 3332 USA
More informationMEAN VALUE ESTIMATES OF z Ω(n) WHEN z 2.
MEAN VALUE ESTIMATES OF z Ωn WHEN z 2 KIM, SUNGJIN 1 Introdction Let n i m pei i be the prime factorization of n We denote Ωn by i m e i Then, for any fixed complex nmber z, we obtain a completely mltiplicative
More information