Lagrangian cell-centered conservative scheme
|
|
- Agatha Lewis
- 6 years ago
- Views:
Transcription
1 Appl. Math. Mech. -Engl. E., 33(10), (2012) DOI /s c Shanghai University an Springer-Verlag Berlin Heielberg 2012 Applie Mathematics an Mechanics (English Eition) Lagrangian cell-centere conservative scheme uan-wen GE ( ) (Institute of Applie Physics an omputational Mathematics, Beijing , P. R. hina) Abstract This paper presents a Lagrangian cell-centere conservative gas ynamics scheme. The piecewise constant pressures of cells arising from the current time sub-cell ensities an the current time isentropic spee of soun are introuce. Multipling the initial cell ensity by the initial sub-cell volumes obtains the sub-cell Lagrangian masses, an iviing the masses by the current time sub-cell volumes gets the current time subcell ensities. By the current time piecewise constant pressures of cells, a scheme that conserves the momentum an total energy is constructe. The vertex velocities an the numerical fluxes through the cell interfaces are compute in a consistent manner ue to an original solver locate at the noes. The numerical tests are presente, which are representative for compressible flows an emonstrate the robustness an accuracy of the Lagrangian cell-centere conservative scheme. Key wors sub-cell force, Lagrange cell-centere scheme, Lagrangian cell-centere conservative gas ynamics scheme, piecewise constant pressure of cell hinese Library lassification O Mathematics Subject lassification 76N15, 65M06 1 Introuction We are intereste in solving the gas ynamics equation written in the Lagrangian form. In the present paper, we aim at presenting a Lagrangian cell-centere conservative gas ynamics scheme. This scheme is the extension of the cell-centere Lagrangian scheme for twoimensional compressible flow problem presente in Ref. [1]. Before escribing the metho, we briefly give a historical overview on the Lagrangian schemes. In Lagrangian hyroynamics methos, a computational cell moves at the flow velocity. In practice, this means that the cell vertices move at a compute velocity, an the cell faces are uniquely specifie by the vertex positions. This ensures that there is no mass flux crossing the bounary of the Lagrangian moving cell. Thus, Lagrangian methos can capture contact iscontinuity sharply in multimaterial flui flows. However, in the Lagrangian framework, one has to iscretize not only the gas ynamics equation but also the vertex motion in orer to move the mesh. Moreover, the numerical fluxes of the physical conservation laws must be etermine in a compatible way with the vertex velocity such that the geometric conservation law (GL) can be satisfie, namely, the rate of change of a Lagrangian volume has to be compute coherently with the Receive Aug. 29, 2011 / Revise May 2, 2012 Project supporte by the National Natural Science Founation of hina (No ) orresponing author uan-wen GE, Associate Professor, Ph.D., ge quanwen@iapcm.ac.cn
2 1330 uan-wen GE noe motion. This critical requirement is the cornerstone of any Lagrangian multi-imensional scheme. The most natural way to solve this problem is to use the staggere iscretization. In the staggere iscretization, the position, the velocity, an the kinetic energy are centere at the points, while the ensity, the pressure, an the internal energy are within the cells. The issipation of the kinetic energy into the internal energy through shock waves is ensure by an artificial viscosity term. Since the seminal works of von Neumann an Richtmyer [2] an Wilkins [3], a lot of evelopments have been mae to improve the accuracy an the robustness of the staggere hyroynamics [4 6]. More specifically, the construction of the compatible staggere iscretization leas to a scheme that conserves the total energy in a rigorous manner [7 8]. We also note the recent evelopment of a variational multi-scale stabilize approach in the finite element computation of Lagrangian hyroynamics piecewise linear approximation for the variables [9 10]. The case of 1 /P 0 finite element was stuie in Ref. [11], where the kinematic variables were represente using a piecewise linear continuous approximation, an the thermoynamic variables utilize a piecewise constant representation. An alternative option to the previous iscretizations is to erive a Lagrangian scheme base on the Gounov metho [12]. In comparison with staggere iscretizations, the Gounov-type methos exhibit the goo properties, i.e., they are naturally conservative, they o not nee an artificial viscosity, an they allow a straightforwar implementation of conservative remapping methos when they are use in the context of the arbitrary Lagrangian-Eulerian (ALE) strategy. In the Gounov-type methos, all conserve quantities incluing momentum an cell velocities are cell-centere. The cell-face quantities incluing a face-normal component of the velocity are available from the solution to an approximate Riemann problem at each cell face. However, the vertex velocity must be etermine to move the mesh. Aessio et al. [13] use a weighte least squares algorithm to compute the vertex velocity on the conition that the vertex velocity projecte in the normal irection of a face is equal to the Riemann velocity on that face. It turns out that this algorithm is capable of generating aitional spurious components in the vertex velocity fiel. Hence, it leas to the artificial gri motion which requires a very expensive treatment [14]. This rawback comes probably from the fact that the flux computation is not compatible with the noe isplacement, an hence the GL is not satisfie. An important achievement concerning the compatibility between the flux iscretization an the vertex velocity computation was introuce by Despŕes an Mazeran [15]. They presente a scheme in which the interface fluxes an the noe velocity were compute coherently with an approximate Riemann solver locate at the noes. This original approach leas to a first-orer conservative scheme which satisfies a local semiiscrete entropy inequality. The multi-imensional higher-orer extension of this scheme was evelope in Ref. [16]. A thorough stuy [1] of the properties of the Despŕes- Mazeran [15] noal solver showe a strong sensitivity to the cell aspect ratio, which coul lea to severe numerical instabilities. To overcome this ifficulty, Maire et al. [1] propose an alternative scheme that successfully solve the aspect ratio problem an kept the compatibility between flux iscretization an vertex computation. The main feature lay in the iscretization of the pressure graient, which meant that two pressures were set at each noe of a cell, an the pressure of each noal was associate with the irection of the outwar unit normal vectors relate to the eges originating from the noe. These noal pressures were linke to the noal velocity by the half Riemann problems. However, Maire et al. [1] faile to remove the numerical ifficulties of anomalous gri istortion that have plague the Lagrangian schemes because of their inception. The bane of Lagrangian hyroynamics calculations is the premature breakown of the gri topology that results in severe egraation of accuracy an run termination before the assumption of a Lagrangian zonal mass is prove to be vali. At short spatial gri scales, this is usually calle the hourglass moe or the keystone motion associate, in particular, with constraine gris such as quarilaterals an hexaherons in two an three imensions, respectively. At longer spatial lengths relative to the gri spacing, there is the spurious vorticity
3 Lagrangian cell-centere conservative scheme 1331 or the long-thin zone problem. In both cases, the result is anomalous gri istortion an tangling that has nothing to o with the actual solution. In this work, we show how such motion can be eliminate by the proper usage of piecewise constant pressures of cells. These piecewise constant pressures of cells arise from the current time isentropic spee of soun of cells an the current time sub-cell ensities. Multipling the initial cell ensity by the initial sub-cell volumes obtains the sub-cell Lagrangian masses, an iviing the sub-cell Lagrangian masses by the current time sub-cell volumes gets the current time sub-cell ensities. This paper is organize as follows. The cell-centere Lagrangian scheme for two-imensional compressible flows is introuce in Section 2. We also introuce the concept of sub-cell forces borrowe from the staggere iscretization framework [7] to erive a general form of the cell-centere iscretization. The Lagrangian cell-centere conservative gas ynamics scheme is introuce in Section 3. Extensive numerical experiments are reporte in Section 4. The experiments show not only the robustness an accuracy of the present metho but also the ability to successfully hanle complex two-imensional flows. onclusions an perspectives are given in Section 5. 2 Notation an previous work 2.1 Lagrangian hyroynamics Let D be an open subset of R 2, which is fille with an invisci ieal flui an equippe with the orthonormal frame (O, X, Y ) an the orthonormal basis (e X, e Y ). Define the unit vector e Z = e X e Y. Discretize the equations of Lagrangian hyroynamics as follows: V (t) V (t) V (t) V (t) ρ V = 0, (1) V U NS = 0, (2) S(t) ρ UV + PNS = 0, (3) S(t) ρev + PU NS = 0, (4) S(t) where enotes the Lagrangian time erivation. V (t) is the moving control volume, an S(t) is the bounary of V (t). ρ, U = (u, v) T,, an E are the specific volume, the velocity, the pressure, an the specific total energy of the flui, respectively. N enotes the outwar unit normal vector relate to the moving bounary S(t). Equations (1), (3), an (4) express the conservation of mass, momentum, an total energy, respectively. We note that Eq.(2) is also name the GL, an it is equivalent to the following local kinematic equation: X = U, X(0) = x, (5) where X stans for the coorinates efining the control volume surface at time t > 0, an x stans for the coorinates at time t = 0. Then, X = X(x, t) is implicitly efine by the local kinematic equation, which is also calle the trajectory equation. This enables us to efine the map M t : V (0) V (t), X X(x, t), where X is the unique solution to (5). With fixe t, this map transforms each flui particle from the position at time t = 0 to the position at time t. The thermoynamical closure of the
4 1332 uan-wen GE set of Eqs.(1), (3), an (4) is obtaine by the aitional equation of state in the form of P = P(ρ, ε), (6) where the specific internal energy ε is relate to the specific total energy by ε = E 1 2 U Notation an assumptions Let us consier the physical omain V (0) initially fille with the flui. We assume that we can map it by a set of quarilateral cells Ω (0) ( = 1, 2,, I) without gaps or overlaps, where I is the total number of the cells. Each cell is assigne a unique inex an enote by Ω (0). Using the map M t, we set Ω (t)=m t (Ω (0)). Here, we assume that Ω (t) is still a quarilateral, that is, the map M t is a continuous an linear function over each element of the mesh. Each point (vertex) of the mesh Ω (t) is assigne a unique inex, an we enote by () the counter clockwise orere list of points of cells. The vertices of cells Ω (t) are perioic, i.e., ()+1 = 1, 0 = (). () is the set of the cells aroun point (vertex), an the numbering is also perioic, i.e., Ω ()+1 (t) = Ω 1 (t). We use inex f to enote a generic face of cell Ω (t). L f is the length, Nf is the outwar unit normal vector relate to the face, an F is the set of faces of cells (see Fig. 1). We also introuce the face fluxes Uf, Π f, an (ΠU) f as follows: Uf = 1 L US, f f Π f = 1 L PS, f f (ΠU) f = 1 PUS. L f f (7) Fig. 1 Notation relate to quarilateral cell Ω (t) To get the iscrete evolution equations for the primary variables (τ, U, E), the control volume formulation is use to the quarilateral cell Ω (t). Let m enote the mass of cell Ω (t), which is a constant after Eq.(1). For a flow variable φ, we introuce its mass average value over the cell Ω (t), i.e., φ = 1 ρφv. m Ω (t)
5 Lagrangian cell-centere conservative scheme 1333 Then, Eqs. (1), (3), an (4) are τ m U m + E m + f F() f F() f F() L f U f N f = 0, L f Π f N f = 0, L f (ΠU) f N f = 0. (8) The local kinematic equations in a iscrete form at the point are X = U, X (0) = x, where X enotes the coorinates of point at time t > 0, x is the initial position, an U is the velocity of point. Equation (8) represents the face flux iscretization of the Lagrangian hyroynamics equations for the iscrete variables. In orer to compute the time evolution of the flow variables, we nee to calculate the face fluxes Uf, Π f, an (ΠU) f. Moreover, we also nee to compute the point velocity U to move the mesh. 2.3 ompatible iscretization of GL Since m /ρ = V, m, ρ, an V are the mass, the ensity, an the volume of cell Ω (t), respectively. Equation (2) can be rewritten as follows: V f F() L f U f N f = 0. (9) V is the volume of cell Ω (t), which is a function of coorinates X of point for (). We compute this volume by performing the triangular ecomposition of cell Ω (t) as follows: which is shown in Fig.2. V = 1 2 () (X X +) e Z, Fig. 2 Triangular ecomposition of quarilateral cell Ω (t) The time ifferentiation of this equation leas to V = () 1 2 (L N + L +N +) U, (10) where L an L + are the lengths, an N an N + are the outwar unit normal vectors relate to the eges an +. By shifting inices in the previous summation,
6 1334 uan-wen GE Eq.(10) becomes V = () 1 2 L +N + (U + U +). (11) Now, the comparison between Eqs.(9) an (11) shows that they are equivalent uner the conition U f = 1 2 (U + U +). (12) The only way to satisfy the compatibility conition (12) is to first compute the point velocity U an euce the face velocity Uf. Then, the compatibility of the face iscretization of the GL with the rate of change of the cell volume can be ensure. Introuce notation L = 1 2 L, N = N, L = 1 2 L +, N = N +. Equation (10) can be rewritten as V = () L N U, where N stans for the unit corner vector relate to point an cell Ω (t) efine by L N = L N + L N. 2.4 omputation of momentum flux Maire et al. [1] iscretize the momentum flux by introucing two pressures at each noe of cell Ω (t). These pressures are enote by Π an Π, as shown in Fig.3. They can be seen as noal pressures viewe from cell Ω (t) an relate to the two eges impinging at noe. Using these noal pressures, Maire et al. [1] propose the following efinition of the iscrete graient operator ( P) over cell Ω (t): ( P) = 1 V () (L Π N + L Π N ). Fig. 3 Localization of noal pressures given by half Riemann problems at point Using the iscrete graient operator, the momentum equation is rewritten as m U + () (L Π N + L Π N ) = 0. (13)
7 Lagrangian cell-centere conservative scheme 1335 Maire et al. [1] obtaine a noal flux iscretization for the momentum equation which is equivalent to its face flux iscretization (11) provie that the momentum face flux is written as Π f = 1 2 (Π + Π ). Introuce force F, which is a sub-cell force relate to point an cell Ω (t), Then, the momentum equation can be written as F = L Π N + L Π N. (14) m U + () F = 0. (15) Since the velocity of eges an + in the vicinity of point is equal to the noal velocity, the noal pressures Π an Π can be compute using the following half approximate Riemann problems: P Π = Z (U U ) N, P Π = Z(U U ) N. (16) Here, Z an Z are the mass fluxes swept by the waves. To etermine the coefficients, we follow the approach by setting Z = ρ a + ρ Γ (U U ) N, Z = ρ a + ρ Γ (U U ) N, (17) where a is the local isentropic spee of soun, an Γ is a material-epenent parameter that is given in terms of the ensity ratio in the limit of very strong shocks. In the case of Gamma law gas, one gets Γ = (γ +1)/2. Utilizing Eq.(14), the sub-cell force F can be rewritten as F = L P N M (U U ), (18) where N stans for the unit corner vector relate to point an cell Ω (t) efine by L N = L N + L N, an M = Z L (N N ) + Z L (N N ) is a 2 2 symmetric positive efinite matrix. 2.5 omputation of total energy flux The time rate of change of total energy is equal to the summation of the work performe by the sub-cell forces F over cell Ω (t). The substitution of the sub-cell force efinition of F in the previous equation leas to the following noe flux iscretization of the total energy equation: m E + () (L Π N + L Π N ) U = 0. (19) 2.6 onservation of momentum an total energy To examine the momentum conservation, let us write the global balance of momentum without taking the bounary conitions into account. The summation of the momentum equation (13) over all the cells leas to ( ) m U = F. (20) ()
8 1336 uan-wen GE Switching the summation over cells an the summation over noes in the right-han sie of Eq.(20) yiels ( ) m U = F, (21) () where () is the set of cells aroun point. The summation of the total energy equation (19) over all the cells yiels ( ) m E = F U. (22) () Switching again the summation over cells an the summation over noes in the right-han sie of Eq. (22) yiels ( ) m E = F U. (23) () Then, the total energy conservation is ensure provie that the sub-cell forces satisfy the following conition: F = 0. (24) () 2.7 onstruction of noal solver By Eq.(18), the sufficient conition (24) can be rewritten as Set () (L P N M (U U )) = 0. (25) M = () M. Then, the system satisfie by the point velocity U is written as M U = () (L P N + M U ). (26) 2.8 Problem statement V = (L N + L N) U, () U m + (L Π N + L Π N) = 0, () E m + (L Π N + L Π N) U = 0, M U = () (L P N + M U ), () P Π = Z (U U ) N, P Π = Z(U U ) N. (27)
9 Lagrangian cell-centere conservative scheme 1337 Maire et al. [1] claime that, when conition (24) is satisfie, the momentum an the total energy are conserve. However, we fin that conition (24) is incorrect, because two pressures are set at each noe of cell Ω (t), an the pressure iscontinuous line exists insie cell Ω (t), as shown in Fig.4. Similarly, the pressure iscontinuous line ( = 1, 2,, ()) exists insie each cell Ω (t) ( = 1, 2,, ()). Due to the presence of these pressure iscontinuous lines ( = 1, 2,, ()), conition (24) is incorrect. Our purpose is to present a Lagrangian cell-centere conservative gas ynamics scheme. Instea of following the proceure employe by Maire et al. [1], we fin it more transparent an physically motivate to utilize a piecewise constant pressure of cell Ω (t). We next efine the piecewise constant pressure of cell Ω (t). In Fig.5, we show a quarilateral cell labelle Ω (t), single out one of its efining points (vertices), an label it as. We efine the sie mipoints of a given cell Ω (t) an label them as B 1, B 2, B 3, an B 4. Using the cell points (vertices), +, k, an an the sie mipoints B 1, B 2, B 3, an B 4 of a given cell Ω (t), we subivie a quarilateral Ω (t) into eight triangle sub-cells B 1, B 2, + B 2, + B 3, k B 3, k B 4, B 4, an B 1. onsiering the triangle sub-cells B 1, B 2, + B 2, + B 3, k B 3, k B 4, B 4, an B 1, we enote the pressures by Π, Π, Π, Π, Π, Π, Π, an Π. Then, the pressure of + + k k cell Ω (t) is enote as the piecewise constant pressures Π, Π, Π +, Π +, Π k, Π k, Π, an Π in triangle sub-cells B 1, B 2, + B 2, + B 3, k B 3, k B 4, B 4, an B 1. onsiering the triangle sub-cell B 2, we enote its volume, pressure, mass, an ensity by Ṽ, Π, m, an ρ, respectively. onsiering the triangle sub-cell B 1, we enote its volume, pressure, mass, an ensity by Ṽ, Π, m, an ρ, respectively. The forces ue to these piecewise constant pressures can be calculate along the bounaries of these sub-cells. Then, we introuce force F, which is a sub-cell force relate to point an cell Ω (t), i.e., where F = L Π N + L Π N + D Π + D Π, (28) Π = a2 ( ρ ρ ) + P, Π = a2 ( ρ ρ ) + P, (29) in which a is the isentropic spee of soun of cell Ω (t), the length an the outwar unit normal vector D are relate to ege of triangle sub-cell B 1, an the length an the outwar unit normal vector D are relate to ege of the triangle sub-cell B 2. Fig. 4 Notation relate to noal solver at point of cell Ω (t) Fig. 5 Localization of piecewise pressure at point viewe from cell Ω (t)
10 1338 uan-wen GE 3 onservative spatial approximation 3.1 omputation of momentum flux Using the force, we propose the following efinition of the iscrete graient operator over cell Ω (t): ( P) = 1 V () Then, the momentum equation can be written as m U m U + () + () (L Π N + L Π N + D Π + D Π ). (30) (L Π N + L Π N + D Π + D Π ) = 0. (31) F = 0. (32) 3.2 omputation of total energy flux The time rate of change of total energy is equal to the summation of the work performe by the sub-cell forces F over cell Ω (t), i.e., m E + () F U = 0. (33) Substituting the sub-cell force F in Eq.(33) yiels the noe flux iscretization of the total energy equation, i.e., m E + () (L Π N + L Π N + D Π + D Π ) U = 0. (34) 3.3 Noe flux iscretization for cell Ω (t) Gathering the results from the previous sections, we write the semi-iscrete Lagrangian cell-centere conservative gas ynamics scheme of noe flux iscretization for the unknowns (τ, U, E ) as follows: τ m = (L N + L N) U, () U m + (L Π N + L Π N + D Π + D Π ) = 0, (35) () E m + (L Π N + L Π N + D Π + D Π ) U = 0. () 3.4 onservation of momentum an total energy To examine the momentum conservation, let us write the global balance of momentum without taking the bounary conitions into account. The summation of the momentum equation (32) over all the cells leas to ( ) m U = () F. (36)
11 Lagrangian cell-centere conservative scheme 1339 Switching the summation over cells an the summation over noes in the right-han sie of Eq.(36) yiels ( ) m U = () F. (37) Then, the momentum conservation is ensure provie that the sub-cell forces satisfy the following conition: () F = 0. (38) Since the sum of all forces acting on one point is equal to zero (see Fig. 6), conition (38) is satisfie, an the momentum is conserve. We claim that the total energy is also conserve. To emonstrate this property, we perform the summation of the total energy equation (33) over all the cells, ( ) m E = () F U. (39) Then, we switch again the summation over cells an the summation over noes in the right-han sie of Eq.(39) to obtain ( ) m E = () ( () F ) U. (40) F = 0. (41) Due to Eq.(38), the term between parentheses in the right-han sie of Eq.(40) is null, an the total energy is conserve. Fig. 6 Notation relate to noal solver at point 3.5 onstruction of solver at vertex velocity U We can introuce only one pressure at point. This pressure is etermine by the half Riemann problem efine in the irection of the unit corner vector N, i.e., P Π = Z (U U ) N, (42)
12 1340 uan-wen GE where Z is the acoustic impeance of cell Ω (t), Z = ρ a + ρ Γ (U U ) N. (43) This amounts to efine only one noal pressure Π for each cell Ω (t) that surrouns point, as shown in Fig. 7. The sub-cell force F corresponing to the single noal pressure Π is F = L P N M (U U ), (44) where N stans for the unit corner vector relate to point an cell Ω (t) efine by L N = L N + L N, an M = Z L N N is a 2 2 symmetric positive efinite matrix. Since the sum of all forces acting on point is equal to zero, we get () (L P N M (U U )) = 0. (45) Using conitions (28), (29), an (38), we write the linear system satisfie by the components (u, v ) of velocity U of an internal vertex. This generic vertex is not on the bounary of the omain so that it is surroune by cells Ω (t)( = 1, 2,, ()), as shown in Fig. 7. () (L a2 ( ρ ρ )N + L a2 ( ρ ρ )N + D a2 ( ρ ρ ) () + D a2 ( ρ ρ )) + (L P N + L P N ) = 0. (46) Fig. 7 Notation of one pressure Π at point of cell Ω (t) Substituting Eq.(45) to Eq.(46), we can obtain the following equation satisfying the point velocity U : () (L a 2 ( ρ ρ )N + L a2 ( ρ ρ )N + D a 2 ( ρ ρ ) () + D a2 ( ρ ρ )) + ( M (U U )) = 0. (47)
13 Set M = () () M U = Lagrangian cell-centere conservative scheme 1341 M. Then, the system satisfie by the point velocity U is written as () M U (L a2 ( ρ ρ )N + L a2 ( ρ ρ )N + D a2 ( ρ ρ ) + D a2 ( ρ ρ )). (48) Denote (u, v ) by the components of velocity U of an internal vertex satisfying { A u + v = S X, u + B v = S Y, (49) where A, B, an are efine by A = B = = () () () Z (L (N,X )2 + L (N,X )2 ), Z (L (N,Y )2 + L (N,Y )2 ), Z (L (N,X )(N,Y ) + L (N,X )(N,Y )). The right-han sie of Eq.(49) is the components (S X, S Y ) of vector Š efine by Š = () () M U (L a2 ( ρ ρ )N + L a2 ( ρ ρ )N + D a2 ( ρ ρ ) + D a2 ( ρ ρ )). The eterminant of Eq.(49) is = A B 2. We show that it is always positive. In orer to simplify the notation, we efine ( W X = Z1 L 1 N1, ) Z 1 L 1,X N1,, Z () L () N (), Z () L (),X,X N (),,X ( W Y = Z1 L 1 N1, ) Z 1 L 1,Y N1,, Z () L () N (), Z () L (),Y,Y N ().,Y We can immeiately obtain A = W X 2, B = W Y 2, an = W X, W Y, where, is the inner scalar prouct of R 2, an is its associate norm. From the auchy-schwarz inequality, we know that A B 2 0. In fact, A B 2 = 0 if an only if one of the two vectors is null or W X an W Y are colinear. This situation is generally impossible unless the eges aroun noe merge into a single line. The system (49) has always a unique solution which etermines the velocity (u, v ). 3.6 Bounary conitions In this paragraph, we explain our implementation of bounary conitions which is consistent with our internal solver. In the Lagrangian formalism, we have to consier two types of bounary conitions, the pressure or the normal component of the velocity. We use the same
14 1342 uan-wen GE type of notation as in Subsection 3.5. is on the bounary. It is surroune by () cells containe in the omain. There are () + 1 eges impinging on. They are numbere counter clockwise, as shown in Fig. 8. The first Ω 1 an the last cells Ω () have eges M 1 an M ()+1 on the bounary. The outwar normal vectors to the two bounary eges M 1 an M ()+1 coming out of are enote by N 1 respectively. an N() coherently with our notations, Fig. 8 Notation for bounary conitions ase1 A prescribe pressure We enote by Π 1 an Π ()+1 the pressures that are impose on the bounary eges M 1 an M ()+1 (see Fig. 8), respectively. For the conservation relation, we make a balance aroun the vertex ue to the bounary conitions an the conitions (45). We get () Set M = (L P N M (U U )) = L 1 N1 Π 1 + L () () M U = N () Π ()+1. (50) M. Then, the system satisfie by the point velocity U is written as () The components of U satisfy () M U + L P N L 1 N1 Π 1 L() N () Π ()+1. (51) { A u + v = S X, u + B v = S Y. (52) The coefficients A, B, an have alreay been efine in the previous paragraph. The right-han sie of Eq. (52) is the components ( S X, S Y ) of vector S efine by S = () () M U + L P N L 1 N1 Π 1 L () N () Π ()+1. The eterminant of Eq. (52) is = A B 2. It is positive, an system (52) always has a unique solution which etermines the velocity (u, v ). ase2 A prescribe normal velocity Let W 1 an W ()+1 be the values of the prescribe normal velocities on the bounary eges M 1 an M ()+1 coming out of. We istinguish the following two cases. (i) N 1 an N() are not colinear.
15 Lagrangian cell-centere conservative scheme 1343 In this case, the value of the point velocity U is given by the bounary conitions, an the components (u, v ) of the vertex velocity are the solution to the following linear system: N 1,X u + N1,Y v = W 1, N (),X u + N(),Y v = W ()+1. This linear system (53) always has a unique solution. Since the normal vectors are not colinear, the eterminant is not zero. (ii) N 1 an N() are colinear. In this case, the point velocity U is not given irectly, an we nee to know the balance of momentum aroun vertex taking the bounary conitions into account. We get () Set M = (L P N M (U U )) = (L 1 N1 + L() () M U = N () (53) )Π. (54) M. Then, the system satisfie by the point velocity U is written as () () M U + L P N (L 1 N1 + L() N () )Π, (55) where Π is an average pressure applie on the external sie of eges M 1 an M ()+1 coming out of. The pressure is unknown. However, we have an aitional equation corresponing to the bounary conition, i.e., (L 1 N1 + L() N () ) U = L 1 W 1 + L () W ()+1. (56) ombining Eq.(55) with Eq. (56), we get A u + v + GΠ = ŜX, u + B v + HΠ = ŜY, Gu + Hv = L1 W 1 + L () W ()+1, (57) where ŜX an ŜY are the components of vector Ŝ efine by Ŝ = () () M U + L P N. The coefficients A, B, an have alreay been efine in the previous paragraph, an we set G = L 1 N1 +,X L() N () an H =,X L1 N1 +,Y L() N (). The eterminant of (57),Y is = A H GH B G 2. Using the fact that A B 2 > 0, we can show that < 0, an Eq. (57) always has a unique solution (u, v ). 3.7 Time iscretization an mesh motion In this section, we iscretize the system with respect to time that escribes the evolution of the physical variables (τ, U, E ) in cell Ω (t). We compute the volume, momentum, an total energy fluxes Uf N f, Π f N f, an (ΠU) f N f on face + by the volume, momentum,
16 1344 uan-wen GE an total energy fluxes on points an + as follows: Uf N f = 1 2 (U + U +) Nf, L f Π f N f = 1 2 ((L Π N + L Π N + D Π + D Π ) +(L Π N + L Π N + D Π D Π )), + + L f (ΠU) f N f = 1 2 ((L Π N + L Π N + D Π + D Π ) U +(L Π N + L Π N + D Π D Π ) U + + +), where f enotes face + of cell Ω (t). Then, system (35) is written as τ m = 1 f F() 2 L f (U + U +) Nf, U m + 1 f F() 2 ((L Π N + L Π N + D Π + D Π ) +(L Π N + L Π N + D Π D Π )) = 0, + + E m + 1 f F() 2 ((L Π N + L Π N + D Π + D Π ) U +(L Π N + L Π N + D Π D Π ) U + + +) = 0. (58) For the time iscretization of Eq.(58), we use a classical forwar Euler scheme. We assume that the physical properties U,n, E,n, ρ,n, a,n, an P,n in Ω (t) an the geometrical characteristics X,n, Y,n ( = 1, 2,, ()) of cell at the beginning of the time step, i.e., t n, are known. We are going to compute their values at t n+1, an we set t = t n+1 t n. The noal solver (49) an system (52), (53), an (57) provie the vertex velocities (u, v ) from the physical variables an geometry characteristics evaluate at time t n. The explicit time integration of the trajectory equation provies the location of vertices for any time t [t n, t n+1 ], { X (t) = X,n + (t t n )u, Y (t) = Y,n + (t t n )v. (59) By this way, we get the location of the vertices at time t n+1, { X,n+1 = X,n + tu, Y,n+1 = Y,n + tv. (60) Then, we euce that L +N +(t) = ((Y +(t) Y (t)), (X +(t) X (t))) T is linear in time. Thus, tn+1 t n (L +N +)(t) = t 2 ((L +N +) n + (L +N +) n+1 ). The last results enable us to write an approximation of the volume equation which is coherent with the mesh motion, namely, m (τ,n+1 τ,n ) t 4 F() ((L f Nf ) n + (L f Nf ) n+1 )(U + U +) = 0. (61) f=1
17 Lagrangian cell-centere conservative scheme 1345 The approximation of the momentum an the total energy equation is fully explicit to conserve the momentum an the total energy exactly. The lengths of eges an the normal vectors have the same efinition as previously. We get m (U,n+1 U,n ) + t f F() 1 2 ((L Π N + L Π N + D Π + D Π ) +(L Π + +N + + L Π + +N + + +D Π D Π + +)) = 0, (62) an for the total energy, m (E,n+1 E,n ) + t f F() 1 2 ((L Π N + L Π N + D Π + D Π ) U +(L Π + +N + + L Π + +N + + +D Π D Π + +) U +) = 0. (63) The time step t is evaluate following criterion. At t n, for each cell Ω (t), we enote by λ,n the minimal value of the istance between two noes of cell Ω (t). We efine t = λ,n min, (64),2,,I ã,n where ã,n is the maximal value of the local isentropic spees of soun of triangle sub-cells an cell Ω (t). 3.8 Description of algorithm (i) Initialization We know ensity ρ 0 an coorinate of vertices (X0, Y 0 ) ( = 1, 2,, ()) of cells Ω (t) ( = 1, 2,, I) at the initial time, where I is the total number of cells. We compute the triangle sub-cell masses m, m, m, m, m, m, m an m from the initial + + k k ensity ρ 0 an the coorinate of vertices (X0, Y 0 ) ( = 1, 2,, ()) of cells Ω (t) ( = 1, 2,, I). At time t n, in each cell Ω (t) ( = 1, 2,, I), we know the flui variables U,n, E,n, ρ,n, a,n, an P,n an the geometrical characteristics (X,n, Y,n ) ( = 1, 2,, ()). We compute the triangle sub-cell volumes (Ṽ ) n, (Ṽ ) n, (Ṽ + ) n, (Ṽ + ) n, (Ṽ k ) n, (Ṽ k ) n, (Ṽ ) n, an (Ṽ ) n an ensities of triangle sub-cells ( ρ ) n, ( ρ ) n, ( ρ + ) n, ( ρ + ) n, ( ρ k ) n, ( ρ k ) n, ( ρ ) n, an ( ρ ) n. Using pressure P,n an ensities of triangle sub-cell ( ρ ) n, ( ρ ) n, ( ρ + ) n, ( ρ + ) n, ( ρ k ) n, ( ρ k ) n, ( ρ ) n, an ( ρ ) n, we get the isentropic spee of soun of triangle sub-cells (ã ) n, (ã ) n, (ã + ) n, (ã + ) n, (ã k ) n, (ã k ) n, (ã ) n, an (ã ) n. We compute the pressures of triangle sub-cells ( Π ) n, ( Π ) n an the sub-cell force ( F ) n by solving Eqs.(28) an (29). (ii) Noal solver For each internal vertex, we compute velocity (u, v ) by solving the linear system (49), n an for each bounary vertex, we compute velocity (u, v ) by solving the linear system (52), n (53), an (57). (iii) Time step limitations We compute t by solving Eq. (64). (iv) Upate the geometrical quantities We compute (X,n+1, Y,n+1 ) from Eq. (60). Then, we euce L +,n+1 an N +,n+1. (v) Upate the physical variables We get τ,n+1, U,n+1, an E,n+1 from Eqs.(61) (63).
18 1346 uan-wen GE (vi) Equation of state The internal energy ε,n+1 is given by ε,n+1 = E,n+1 U,n+1 2 /2. Then, we get the pressure P,n+1 an the isentropic soun spee a,n+1 from the equation of state. 4 Numerical results 4.1 Saltzman s shock tube We consier the movement of a planar shock wave on the artesian gris [1]. This is a well-known ifficult test case that can valiate the robustness of the present scheme when the mesh is not aligne with the flui flow. The computational omain is the rectangle (x, y) [0, 1] [0, 0.1]. The initial mesh is obtaine by transforming a uniform cells with the mapping { Xstr = x + (0.1 y)sin(πx), Y str = y. For the material, we use the equation of state of the monoatomic gas (γ = 5/3). The initial state is (ρ 0, P 0, U0 )=(1, 0, 0). The bounary conition at x = 0 is a normal velocity U = 1 (the inflow velocity). On all the other bounaries, we set the wall conitions. Figure 9 shows the initial gris. The exact solution is a planar shock wave that moves at the spee D = 4/3 from left to right. The correct answer for the ensity is 4.0 in the singly shocke region, 10.0 in the oubly shocke region an 20.0 in the triply shocke region. The locations of the shock front are 0.967, 0.927, an 0.96 at t =0.80, t =0.86, an t =0.93, respectively. Figures 10 an 11 show the gris at t =0.80 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme, respectively. Figures 12 an 13 show the contour plots of the ensity at t = 0.80 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme. The gris at t= 0.86 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme are shown in Figs. 14 an 15. The contour plots of the ensity at t= 0.86 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme are shown in Figs.16 an 17. The gris at t= 0.93 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme are shown in Figs. 18 an 19. The contour plots of the ensity at t= 0.93 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme are shown in Figs. 20 an 21. All gri istortion is spurious since this is a purely one-imensional problem. The comparison shows that the Lagrangian cell-centere conservative scheme preserves the Fig. 9 Initial gris for Saltzman s shock tube Fig. 10 Gris for Maire et al. [1] scheme at t= 0.80 Fig. 11 Gris for Lagrangian cell-centere conservative scheme at t=0.80
19 Lagrangian cell-centere conservative scheme 1347 Fig. 12 ontour plots of ensity for Maire et al. [1] scheme at t=0.80 Fig. 13 ontour plots of ensity for Lagrangian cell-centere conservative scheme at t= 0.80 Fig. 14 Gris for Maire et al. [1] scheme at t= 0.86 Fig. 15 Gris for Lagrangian cell-centere conservative scheme at t=0.86 Fig. 16 ontour plots of ensity for Maire et al. [1] scheme at t=0.86 Fig. 17 ontour plots of ensity for Lagrangian cell-centere conservative scheme at t= 0.86
20 1348 uan-wen GE one-imensional solution very well (all layers in the vertical irection) an eliminates the artificial gri istortion an hourglass type motion. The answers are extremely close to the true solution. The robustness of the Lagrangian cell-centere conservative scheme is clearly emonstrate by this test case. The Maire et al. [1] scheme oes not preserve the one-imensional solution, an the shock face is increasing tilte. This error comes probably from the fact that the Maire et al. [1] scheme is not the momentum conservation scheme. Fig. 18 Gris for Maire et al. [1] scheme at t= 0.93 Fig. 19 Gris for Lagrangian cell-centere conservative scheme at t=0.93 Fig. 20 ontour plots of ensity for Maire et al. [1] scheme at t=0.93 Fig. 21 ontour plots of ensity for Lagrangian cell-centere conservative scheme at t= Shock refraction problem We now consier a shock refraction problem in the artesian geometry using the same gris as in Ref. [14]. In this problem, a piston moves from the left sening a shock wave through an initially col material of a unit ensity with a gri that is tapere so that its right bounary is at an angle of 60 relative to the vertical irection. This shock then becomes incient on a secon gri that is slante but uniformly space at an angle of 60. Reflective bounary conitions are applie to the top an bottom of both gris. The region compose of the secon gri has an initial ensity of 1.5. The problem runs for t= 1.30 just before the shock starts to run off the secon gri. In Fig.22, the initial gris are shown. Figures 23 an 24 epict the gris at t=1.30 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme.
21 Lagrangian cell-centere conservative scheme 1349 Figures 25 an 26 show the contour plots of the ensity at t= 1.30 for the Maire et al. [1] scheme an the Lagrangian cell-centere conservative scheme. The importance of this problem is that we have both the physical vorticity an the spurious gri istortion [4,6,14]. The gri in Fig. 23 is noticeably istorte, particularly in the vicinity of the lower part of the interface, where the gri is on the verge of tangling. Figure 24 has no spurious gri istortion an hourglass type motion. The contour plots of the ensity in Fig. 25 o not present the physical vorticity, while the contour plots of the ensity in Fig. 26 present the physical vorticity. These results o not permit one to istinguish between the two computations. However, the Lagrangian cell-centere conservative scheme is clearly superior from the point of view of robustness of the computational physical vorticity an the quality. Fig. 22 Initial gris for shock refraction problem Fig. 23 Gris for Maire et al. [1] scheme at t =1.30 Fig. 24 Gris for Lagrangian cell-centere conservative scheme at t =1.30 Fig. 25 ontour plots of ensity for Maire et al. [1] scheme at t=1.30 Fig. 26 ontour plots of ensity for Lagrangian cell-centere conservative scheme at t =1.30
22 1350 uan-wen GE 5 onclusions In this paper, a Lagrangian cell-centere conservative scheme is evelope. This scheme relies on a genuinely two-imensional noal solver. The main feature of the algorithm is the introuction of a piecewise constant pressure of cells. Using a piecewise constant pressure of cells, we construct the scheme, for which the momentum an the total energy are conservative. The vertex velocities an the numerical fluxes through the cell interfaces are not compute inepenently contrary to stanar approaches, but they are evaluate in a consistent manner ue to an original solver locate at the noes. The Lagrangian cell-centere conservative scheme has goo robustness accoring to the numerical results obtaine from the test cases presente in this paper. The future research is to evelop the higher-orer extension of the present scheme. References [1] Maire, P. H., Abgrall, R., Breil, J., an Ovaia, J. A cell-centere Lagrangian scheme for compressible flow problems. SIAM J. Sci. omput., 29(4), (2007) [2] Von Neumann, J. an Richtmyer, R. D. A metho for the numerical calculations of hyroynamics shocks. J. Appl. Phys., 21, (1950) [3] Wilkins, M. L. alculation of elastic plastic flow. Methos in omputational Physics (e. Aler, B.), Vol.3, Acaemic Press, New York (1964) [4] aramana, E. J. an Shashkov, M. J. Elimination of artificial gri istorsion an hourglass-type motions by means of Lagrangian subzonal masses an pressures. J. omput. Phys., 142, (1998) [5] aramana, E. J., Shashkov, M. J., an Whalen, P. P. Formulations of artificial viscosity for multiimensional shock wave computations. J. omput. Phys., 144, (1998) [6] ampbell, J.. an Shashov, J.. A tensor artificial viscosity using a mimetic finite ifference algorithm. J. omput. Phys., 172, (2001) [7] aramana, E. J., Burton, D. E., Shashov, M. J., an Whalen, P. P. The construction of compatible hyroynamics algorithms utilizing conservation of total energy. J. omput. Phys., 146, (1998) [8] ampbell, J.. an Shashov, M. J. A compatible Lagrangian hyroynamics algorithm for unstructure gris. Selcuk J. Appl. Math., 4(2), (2003) [9] Scovazzi, G., hriston, M. A., Hughes, T. J. R., an Shai, J. N. Stabilize shock hyroynamics: I. a Lagrangian metho. omput. Methos Appl. Mech. Engrg., 196, (2007) [10] Scovazzi, G. Stabilize shock hyroynamics: II. esign an physical interpretation of the SUPG operator for Lagrangian computations. omput. Methos Appl. Mech. Engrg., 196, (2007) [11] Scovazzi, G., Love, E., an Shashkov, M. J. Multi-scale Lagrangian shock hyroynamics on 1/P 0 finite elements: theoretical framework an two-imensional computations. omput. Methos Appl. Mech. Engrg., 197, (2008) [12] Gounov, S. K., Zabroine, A., Ivanov, M., Kraiko, A., an Prokopov, G. Résolution Numérique es Problèmes Multiimensionnels e la Dynamique es Gaz, Eitions Mir, Moscow (1979) [13] Aessio, F. L., arroll, D. E., Dukowicz, J. K., Johnson, J. N., Kashiwa, B. A., Maltru, M. E., an Ruppel, H. M. aveat: a omputer oe for Flui Dynamics Problems with Large Distortion an Internal Slip, Technical Report LA MS, Los Alamos National Laboratory (1986) [14] Dukowicz, J. K. an Meltz, B. Vorticity errors in multiimensional Lagrangian coes. J. omput. Phys., 99, (1992) [15] Despŕes, B. an Mazeran,. Lagrangian gas ynamics in two imensions an Lagrangian systems. Arch. Rational Mech. Anal., 178, (2005) [16] arré, G., Delpino, S., Despŕes, B., an Labourasse, E. A cell-centere Lagrangian hyroynamics scheme on general unstructure meshes in arbitrary imension. J. omput. Phys., 228, (2009)
Stable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationOptimized Schwarz Methods with the Yin-Yang Grid for Shallow Water Equations
Optimize Schwarz Methos with the Yin-Yang Gri for Shallow Water Equations Abessama Qaouri Recherche en prévision numérique, Atmospheric Science an Technology Directorate, Environment Canaa, Dorval, Québec,
More informationarxiv: v1 [physics.flu-dyn] 8 May 2014
Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa
More informationLecture 2 Lagrangian formulation of classical mechanics Mechanics
Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationReALE: A Reconnection Arbitrary-Lagrangian-Eulerian method in cylindrical geometry
ReALE: A Reconnection Arbitrary-Lagrangian-Eulerian metho in cylinrical geometry Raphaël Loubère a, Pierre-Henri Maire b, Mikhail Shashkov c a Institut e Mathématiques e Toulouse, CNRS, Université e Toulouse,
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationSwitching Time Optimization in Discretized Hybrid Dynamical Systems
Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set
More informationLie symmetry and Mei conservation law of continuum system
Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive
More informationSchrödinger s equation.
Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of
More informationCode_Aster. Detection of the singularities and computation of a card of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size
More informationθ x = f ( x,t) could be written as
9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)
More informationII. First variation of functionals
II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent
More informationCode_Aster. Detection of the singularities and calculation of a map of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : DLMAS Josselin Clé : R4.0.04 Révision : Detection of the singularities an calculation of a map of size of
More informationMomentum and Energy. Chapter Conservation Principles
Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an
More informationMath 342 Partial Differential Equations «Viktor Grigoryan
Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More informationThe total derivative. Chapter Lagrangian and Eulerian approaches
Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function
More informationWUCHEN LI AND STANLEY OSHER
CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability
More informationLecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012
CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration
More informationEuler Equations: derivation, basic invariants and formulae
Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationSeparation of Variables
Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical
More informationEuler equations for multiple integrals
Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationIntroduction to the Vlasov-Poisson system
Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationHyperbolic Moment Equations Using Quadrature-Based Projection Methods
Hyperbolic Moment Equations Using Quarature-Base Projection Methos J. Koellermeier an M. Torrilhon Department of Mathematics, RWTH Aachen University, Aachen, Germany Abstract. Kinetic equations like the
More informationHyperbolic Systems of Equations Posed on Erroneous Curved Domains
Hyperbolic Systems of Equations Pose on Erroneous Curve Domains Jan Norström a, Samira Nikkar b a Department of Mathematics, Computational Mathematics, Linköping University, SE-58 83 Linköping, Sween (
More informationChapter 2 Lagrangian Modeling
Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie
More informationTechnion - Computer Science Department - M.Sc. Thesis MSC Constrained Codes for Two-Dimensional Channels.
Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Constraine Coes for Two-Dimensional Channels Keren Censor Technion - Computer Science Department - M.Sc. Thesis MSC-2006- - 2006 Technion
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationStructure preserving finite volume methods for the shallow water equations
Structure preserving finite volume methos for the shallow water equations by Ulrik Skre Fjorholm THESIS for the egree of MASTER OF SCIENCE Master i Anvent matematikk og mekanikk) Department of Mathematics
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationExamining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing
Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June
More informationThe derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)
Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)
More informationdiv(λ(x) u) = f in Ω, (1) u = 0 on Ω, (2) where we denote by Ω = Ω \ Ω the boundary of the domain Ω, under the following assumptions:
Discretisation of heterogeneous an anisotropic iffusion problems on general nonconforming meshes SUSHI: a scheme using stabilisation an hybri interfaces 1 R. Eymar 2, T. Gallouët 3 an R. Herbin 4 Abstract:
More informationVectors in two dimensions
Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication
More informationON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM
ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract
More informationCalculus of Variations
16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t
More informationarxiv: v1 [math-ph] 5 May 2014
DIFFERENTIAL-ALGEBRAIC SOLUTIONS OF THE HEAT EQUATION VICTOR M. BUCHSTABER, ELENA YU. NETAY arxiv:1405.0926v1 [math-ph] 5 May 2014 Abstract. In this work we introuce the notion of ifferential-algebraic
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationAgmon Kolmogorov Inequalities on l 2 (Z d )
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,
More informationSensors & Transducers 2015 by IFSA Publishing, S. L.
Sensors & Transucers, Vol. 184, Issue 1, January 15, pp. 53-59 Sensors & Transucers 15 by IFSA Publishing, S. L. http://www.sensorsportal.com Non-invasive an Locally Resolve Measurement of Soun Velocity
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationPolynomial Inclusion Functions
Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl
More informationNovel entropy stable schemes for 1D and 2D fluid equations
Novel entropy stable schemes for D an D flui equations Eitan Tamor an Weigang Zhong Department of Mathematics, Center for Scientific Computation an Mathematical Moeling (CSCAMM, Institute of Physical Science
More informationON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS
ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear
More informationTopological Sensitivity Analysis for Three-dimensional Linear Elasticity Problem
Topological Sensitivity Analysis for Three-imensional Linear Elasticity Problem A.A. Novotny, R.A. Feijóo, E. Taroco Laboratório Nacional e Computação Científica LNCC/MCT, Av. Getúlio Vargas 333, 25651-075
More informationLagrangian and Hamiltonian Mechanics
Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationSolution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010
NTNU Page of 6 Institutt for fysikk Fakultet for fysikk, informatikk og matematikk This solution consists of 6 pages. Solution to the exam in TFY423 STATISTICAL PHYSICS Wenesay ecember, 2 Problem. Particles
More informationNumerical Investigation of Non-Stationary Parameters on Effective Phenomena of a Pitching Airfoil at Low Reynolds Number
Journal of Applie Flui Mechanics, Vol. 9, No. 2, pp. 643-651, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Investigation of Non-Stationary Parameters on Effective
More information12.11 Laplace s Equation in Cylindrical and
SEC. 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential 593 2. Laplace s Equation in Cylinrical an Spherical Coorinates. Potential One of the most important PDEs in physics an engineering
More informationProblems Governed by PDE. Shlomo Ta'asan. Carnegie Mellon University. and. Abstract
Pseuo-Time Methos for Constraine Optimization Problems Governe by PDE Shlomo Ta'asan Carnegie Mellon University an Institute for Computer Applications in Science an Engineering Abstract In this paper we
More informationLinear First-Order Equations
5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)
More informationNOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,
NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which
More informationThe Ehrenfest Theorems
The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationConservation Laws. Chapter Conservation of Energy
20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action
More informationTime-of-Arrival Estimation in Non-Line-Of-Sight Environments
2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor
More information1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.
Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency
More informationAssignment 1. g i (x 1,..., x n ) dx i = 0. i=1
Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition
More informationQuantum Mechanics in Three Dimensions
Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.
More informationDiscrete Operators in Canonical Domains
Discrete Operators in Canonical Domains VLADIMIR VASILYEV Belgoro National Research University Chair of Differential Equations Stuencheskaya 14/1, 308007 Belgoro RUSSIA vlaimir.b.vasilyev@gmail.com Abstract:
More informationStudy on aero-acoustic structural interactions in fan-ducted system
Stuy on aero-acoustic structural interactions in fan-ucte system Yan-kei CHIANG 1 ; Yat-sze CHOY ; Li CHENG 3 ; Shiu-keung TANG 4 1,, 3 Department of Mechanical Engineering, The Hong Kong Polytechnic University,
More informationLower Bounds for the Smoothed Number of Pareto optimal Solutions
Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.
More informationand from it produce the action integral whose variation we set to zero:
Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine
More informationLecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.
b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference
More informationCOUPLING REQUIREMENTS FOR WELL POSED AND STABLE MULTI-PHYSICS PROBLEMS
VI International Conference on Computational Methos for Couple Problems in Science an Engineering COUPLED PROBLEMS 15 B. Schrefler, E. Oñate an M. Paparakakis(Es) COUPLING REQUIREMENTS FOR WELL POSED AND
More informationALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS
ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an
More informationCONSERVATION PROPERTIES OF SMOOTHED PARTICLE HYDRODYNAMICS APPLIED TO THE SHALLOW WATER EQUATIONS
BIT 0006-3835/00/4004-0001 $15.00 200?, Vol.??, No.??, pp.?????? c Swets & Zeitlinger CONSERVATION PROPERTIES OF SMOOTHE PARTICLE HYROYNAMICS APPLIE TO THE SHALLOW WATER EQUATIONS JASON FRANK 1 an SEBASTIAN
More informationFundamental Laws of Motion for Particles, Material Volumes, and Control Volumes
Funamental Laws of Motion for Particles, Material Volumes, an Control Volumes Ain A. Sonin Department of Mechanical Engineering Massachusetts Institute of Technology Cambrige, MA 02139, USA August 2001
More informationDomain Decomposition and Model Reduction of Systems with Local Nonlinearities
Domain Decomposition an Moel Reuction of Systems with Local Nonlinearities Kai Sun, Rolan Glowinsi, Matthias Heinenschloss, an Danny C. Sorensen Abstract The goal of this paper is to combine balance truncation
More informationTEST 2 (PHY 250) Figure Figure P26.21
TEST 2 (PHY 250) 1. a) Write the efinition (in a full sentence) of electric potential. b) What is a capacitor? c) Relate the electric torque, exerte on a molecule in a uniform electric fiel, with the ipole
More informationCalculus of Variations
Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,
More informationTMA 4195 Matematisk modellering Exam Tuesday December 16, :00 13:00 Problems and solution with additional comments
Problem F U L W D g m 3 2 s 2 0 0 0 0 2 kg 0 0 0 0 0 0 Table : Dimension matrix TMA 495 Matematisk moellering Exam Tuesay December 6, 2008 09:00 3:00 Problems an solution with aitional comments The necessary
More informationGeneralization of the persistent random walk to dimensions greater than 1
PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,
More information4. Important theorems in quantum mechanics
TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this
More informationSolutions to Math 41 Second Exam November 4, 2010
Solutions to Math 41 Secon Exam November 4, 2010 1. (13 points) Differentiate, using the metho of your choice. (a) p(t) = ln(sec t + tan t) + log 2 (2 + t) (4 points) Using the rule for the erivative of
More informationNoether s theorem applied to classical electrodynamics
Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationConvective heat transfer
CHAPTER VIII Convective heat transfer The previous two chapters on issipative fluis were evote to flows ominate either by viscous effects (Chap. VI) or by convective motion (Chap. VII). In either case,
More informationProblem set 2: Solutions Math 207B, Winter 2016
Problem set : Solutions Math 07B, Winter 016 1. A particle of mass m with position x(t) at time t has potential energy V ( x) an kinetic energy T = 1 m x t. The action of the particle over times t t 1
More informationContinuum Mechanics Lecture 4 Fluid dynamics
Continuum Mechanics Lecture 4 Flui ynamics Prof. http://www.itp.uzh.ch/~teyssier Outline - Flui kinematics - Mass an momentum conservation laws - The energy equation - Real fluis - Ieal fluis - Incompressible
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationImplicit Differentiation
Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,
More informationIntroduction to Markov Processes
Introuction to Markov Processes Connexions moule m44014 Zzis law Gustav) Meglicki, Jr Office of the VP for Information Technology Iniana University RCS: Section-2.tex,v 1.24 2012/12/21 18:03:08 gustav
More informationPartial Differential Equations
Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.
More informationThe Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis
The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into
More informationDelocalization of boundary states in disordered topological insulators
Journal of Physics A: Mathematical an Theoretical J. Phys. A: Math. Theor. 48 (05) FT0 (pp) oi:0.088/75-83/48//ft0 Fast Track Communication Delocalization of bounary states in isorere topological insulators
More informationChapter 4. Electrostatics of Macroscopic Media
Chapter 4. Electrostatics of Macroscopic Meia 4.1 Multipole Expansion Approximate potentials at large istances 3 x' x' (x') x x' x x Fig 4.1 We consier the potential in the far-fiel region (see Fig. 4.1
More informationControl Volume Derivations for Thermodynamics
Control olume Derivations for Thermoynamics J. M. Powers University of Notre Dame AME 327 Fall 2003 This ocument will give a summary of the necessary mathematical operations necessary to cast the conservation
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More information