A Provable Stable and Accurate Davies-like Relaxation Procedure Using Multiple Penalty Terms for Lateral Boundaries in Weather Prediction
|
|
- Willa Newman
- 6 years ago
- Views:
Transcription
1 Department of Mathematics A Provable Stable and Accurate Davies-like Relaxation Procedure Using Multiple Penalty Terms for Lateral Boundaries in Weather Prediction Hannes Frenander and Jan Nordström LiTH-MAT-R--24/9--SE
2 Department of Mathematics Linköping University S Linköping, Sweden.
3 A provable stable and accurate Davies-like relaxation procedure using multiple penalty terms for lateral boundaries in weather prediction. Hannes Frenander a, Jan Nordström b a Department of Mathematics, Computational Mathematics, Linköping University, SE Linköping, Sweden (hannes.frenander@liu.se). b Department of Mathematics, Computational Mathematics, Linköping University, SE Linköping, Sweden (jan.nordstrom@liu.se). Abstract A lateral boundary treatment using summation-by-parts operators and simultaneous approximation terms is introduced. The method, that we refer to as the multiple penalty technique, is similar to Davies relaxation and have similar areas of application. The method is proven, by energy methods, to be stable. We show how to apply this technique on the linearized Euler equations in two space dimensions, and that it reduces the errors in the computational domain. Keywords: Davies relaxation, lateral boundary conditions, summation-by-parts, weak boundary conditions, penalty terms. Introduction Accurate numerical calculations on large domains are often unfeasible. To circumvent this issue, a coarse global mesh is used, and inter-spaced local domains are inserted where accurate results are wanted. These local meshes are typically employed where something interesting happens, and where fine meshes are needed to capture the phenomena with a given accuracy. This can, for example, be used to model local weather phenomena or turbulence around airplane wings. However, this model requires that the local domains are matched correctly to the global one. H.C Davies [] introduced in 976 a lateral boundary procedure for weather prediction models where the coarse grid numerical results are interpolated into the local fine grid domain, a Preprint submitted to Elsevier September 26, 24
4 method later referred to as Davies relaxation. Other types of interpolation methods are also used within the weather prediction community, see [5],[2] [3]. In most applications, one aims for accurate results in the local domains, while the accuracy in the global domain is considered given. Therefore, the coupling between the domains are in most cases neglected. There are methods, see [4], where the coupling between the meshes is considered, such that the results from the local domain influence the global one. The results from the global area can be transferred to the local one by using so-called window functions, see [6],[7], which are used to blend the results from the domains. These window functions can, for example, be based on Fourier extensions close to the interface. However, this coupling often introduces stability issues unless optimally done [9],[],[],[2]. In other cases, the global data is mixed into global domain by interpolation [],[8], which is the technique we will adapt in this work. Moreover, we will neglect any coupling between the domains. We demonstrate our technique on the linearized Euler equations, discretized using operators with the summation-by parts (SBP) property [4],[5] and augmented with simultaneous approximation terms (SAT) [6]. This is often referred to as the SBP-SAT technique and a comprehensive review is given in [3]. If well-posed boundary conditions are available, this technique yields stable numerical schemes [7]. In [8], it was shown that if data is available inside the computational domain, additional penalty terms can be applied at these points without ruining stability. This method, which we call the multiple penalty technique (MPT), can be used to assimilate global data into local area models, constructing non-reflecting boundaries, raising accuracy of the scheme and increasing the rate of convergence to steady state, see [8],[2]. More importantly, it can be done in a provable stable way. The MPT is similar to the Davies relaxation mentioned above, but can also be proven stable. The rest of the paper will proceed as follows. In section 2, we derive well-posed boundary conditions for the linearized Euler equations och discretesize it by using the SBP-SAT technique. Stability conditions are derived och multiple penalties are applied such that stability is preserved. In section 3, numerical experiments are performed to illustrate the increased rate of convergence to steady-state when using multiple penalties. Numerical experiments are also performed on time-dependent problems, illustrating an increased convergence rate when MPT is used. The results of this paper are concluded in section 4. 2
5 Figure : Illustration of the problem set-up. The calculations are performed in the central domain where boundary data has been interpolated from the coarse, global domain. 2. Well-posedness and stability We consider the Euler equations in the dimensionless quantities ρ = ˆρ/ˆρ, u = û/û, v = ˆv/ˆv and p = ˆp/(ˆρ c 2 ) where ˆρ is the density, û, ˆv the horizontal and vertical component of the velocity, ˆp the pressure and c the speed of sound. The subscript denotes a reference state. For simplicity, and without restriction, a square domain x, y [, ] is considered and t. This domain is assumed to be contained within a global mesh, from which the boundary data is obtained. See Figure for an illustration. The Euler equations with boundary and initial conditions are U t + A(U)U x + B(U)U y = L w U(, y, t) = g w (y, t) L s U(x,, t) = g s (x, t) L e U(, y, t) = g e (y, t) L n U(x,, t) = g n (x, t) u(x, y, ) = f(x, y) where U = [ρ, u, v, p] T, L w,s,e,n are the boundary operators to be determined and g w,s,e,n are the boundary data. The matrices are u ρ v ρ A(U) = u /ρ u, B(U) = v v /ρ. ρ c 2 u ρ c 2 v 3 ()
6 The problem () is linearized by considering small deviations from the reference state; that is, we insert U = U + U where U is the reference state and U is small. The linearized problem is then transformed according to w = T s U = [p /(ρ c), u, v, p ρ c 2 ] T, with T s = /(ρ c) c 2 which leads to the symmetric and linearized form of the Euler equations,, w t + A s (U )w x + B s (U )w y = L s ww (, y, t) = g w(y, t) L s sw (x,, t) = g s(x, t) L s ew (, y, t) = g e(y, t) L s nw (x,, t) = g n(x, t) w (x, y, ) = f (x, y) (2) where f (x, y) is the perpetuated initial condition, g w,s,e,n are the perpetuated boundary data and ū c v c A s = T s ATs = c ū ū, B s = T s BTs = v c v. ū v Note that A s and B s are now constant symmetric matrices and ū, v are the normalized reference state velocities. A general technique to symmetrize the Euler and Navies-Stokes equations can be found in [2]. 2.. Well-posedness We now determine the boundary operators such that (2) becomes wellposed. If one multiply (2) with w T from the left, integrate over the domain and apply Greens theorem, one arrives at w 2 t = [w T A s w, w T B s w ] ˆnds (3) ω 4
7 where ˆn = [dy, dx] T / dx 2 + dy 2 is the outward pointing normal to the domain and ds = dx 2 + dy 2 is the differential along the boundary. In (3), the norm is defined as w 2 = Ω (w ) 2 dxdy. In our simplified model, the domain Ω is a square. The relation (3) can under these circumstances be written out explicitly as w 2 t = [w T A s w, w T B s w ] ˆnds = x= y= ω w (x,, t) T B s w (x,, t)dx w (, y, t) T A s w (, y, t)dy x= y= w (x,, t) T B s w (x,, t)dx w (, y, t) T A s w (, y, t)dy. All boundary operators that bounds the right-hand side of (4) results in an energy estimate, and well-posedness follows if the correct number of boundary conditions is used. By considering the eigenvalue decomposition of A s = S A Λ A S T A and B s = S B Λ B S T B where Λ A,B are the eigenvalue matrices of A s and B s and S A,B the corresponding similarity transformations, the most direct way to construct well-posed boundary conditions is then to choose (4) L s w = S A Λ + A ST A = A + s L s s = S B Λ + B ST B = B + s L s e = S A Λ A ST A = A s L s n = S B Λ B ST B = B s. (5) In (5), Λ + A,B and Λ A,B refers to the positive and negative parts of the eigenvalue matrices, respectively. Problem (2) can with the operators in (5) be written w t + A s (U )w x + B s (U )w y = A + s w (, y, t) = g w(y, t) B + s w (x,, t) = g s(x, t) A s w (, y, t) = g e(y, t) B s w (x,, t) = g n(x, t) w (x, y, ) = f (x, y). (6) 5
8 By using the eigenvalue decomposition as in (5) to describe the boundary conditions in (6), one obtains w t + A s (U )w x + B s (U )w y = Λ + Aŵ (, y, t) = ĝ w(y, t) Λ + B w (x,, t) = g s(x, t) Λ Aŵ (, y, t) = ĝ e(y, t) Λ B w (x,, t) = g n(x, t) w (x, y, ) = f (x, y). where ŵ = S A w, w = S B w, ĝ w,e = S A g w,e and g s,n = S B g s,n. Here, we study only subsonic velocities, such that ū 2 + v 2 < c 2. Using the boundary operators (7) in (4), one obtains v w 2 t = + x= x= g T s g T n ĝ w T y= y= ĝ T e v g s + v+ c g n v c ū ū ĝ w + ū+ c ĝ e ū c x= x= y= y= w (x,, t)b s w (x,, t)dx w (x,, t)b + s w (x,, t)dx w (, y, t)a s w (, y, t)dx w (, y, t)a + s w (, y, t)dx where we have used Λ A = diag(ū, ū, ū + c, ū c), Λ B = diag( v, v, v + c, v c) and ū 2 + v 2 < c 2. Equation (8) states that the norm w t is bounded by data. Since we also have used the correct number of boundary conditions, (2) is strongly well-posed. We state this result in the following proposition (7) (8) 6
9 Proposition. The problem (2) is strongly well-posed if L s w = A + s, L s s = B + s, L s e = A s, L s n = B s 2.2. Stability and error estimates The semi-discrete approximation of (6) using the SBP-SAT technique is v t + Px Q x I y A s v + I x Py Q y B s v = (Px E Py I 4 )Σ w (Ā+ s v g w)+ (Px (Px (P x E N P I 4 )Σ e (Ā s v g e)+ P y y E I 4 )Σ s ( B + s v g s)+ P y E N I 4 )Σ n ( B s v g n) (9) where P x,y are symmetric and positive definite and Q x,y satisfy the SBPproperty, Q x,y + Q T x,y = B x,y = E + E N = diag(,,...,, ). In (9), Ā ± S = I x I y A ± S and B s ± = I x I y B s ±. The matrices E,N are projection matrices, such that (E ) = (E N ) NN = and zero otherwise. The penalty matrices Σ w,s,e,n will be determined such that stability is obtained and g w,e,s,n are the boundary data. By multiplying (9) with v T P x P y I 4 from the left, one obtains ( v 2 P x P y I 4 ) t + v T (B x P y A s )v + v T (P x B y B s )v = 2α v T (E P y I 4 )(Ā+ s v g w) +2α 2 v T (E N P y I 4 )(Ā s v g e) +2α 3 v T (P x E I 4 )( B s + v g s) +2α 4 v T (P x E N I 4 )( B s v g s) () where we have used Σ w,e = α,2 I x P y I 4, Σ s,n = α 3,4 P x I y I 4. The constants α,2,3,4 are to be determined such that stability is obtained. In (), we have used the SBP-property of Q x,y,: Q x,y + Q T x,y = B x,y = E + E N. By 7
10 putting α,3 = and α 2,4 =, () becomes ( v 2 P x P y I 4 ) t = v T (E P y I 4 )Ā+ s v + 2v T (E P y I 4 )g w + v T (E P y A s )v + v T (E N P y I 4 )Ā s v 2v T (E N P y I 4 )g e v T (E N P y A + s )v v T (P x E I 4 ) B s + v + 2v T (P x E I 4 )g s + v T (P x E Bs )v + v T (P x E N I 4 ) B s v 2v T (P x E I 4 )g n v T (P x E N B s + )v. () By completing squares, () finally becomes ( v 2 P x P y I 4 ) t = (Ā+ s v g w) T (E P y (A s ) + )(Ā+ S v g w) + g w T (E P y (A s ) + )g w + v T (E P y A s )v + (Ā s v g e) T (E N P y (A s ) )(Ā s v g e) g w T (E N P y (A s ) )g w v T (E N P y A + s )v ( B s + v g s) T (P x E (Bs ) + )( B s + v g s) + g s T (P x E (Bs ) + )g s + v T (P x E Bs )v ( B s v g n) T (P x E N (Bs ) )( B s v g n) g n T (P x E N (Bs ) )g n v T (P x E N B s + )v. (2) where (A s ) ± and (Bs ) ± refers to the negative and positive parts of the inverses of A s and B s, respectively. Since (2) is bounded by the data on the right-hand side, the semi-discrete formulation (9) is stable. Note that the terms g w,e(e T,N P y (A s ) ± )g w,e T and g s,n(p T x (A s ) ± )g w,e T mimics the integrals in (8) containing the boundary data and the terms v T (E,N P y A ± s )v and v T (P x E,N B s ± )v mimics the remaining integrals. Hence, (2) mimics the continuous energy estimate, given by (8), if the small damping terms are neglected. 8
11 2.3. Error estimates By inserting the analytic solution, u, into (9) and subtracting it from (9) with the numerical solution v, one obtains e t + Px Q x I y A s e + I x Py Q y B s e = E P (P x (P x (P x (P x P y y I 4 )Σ w Ā + s e+ E N Py I 4 )Σ e Ās e+ Py E I 4 )Σ s B+ s e+ E N I 4 )Σ n B s e + T e (3) where e = v u and T e is the truncation error. By multiplying (3) with e T (P x P y I 4 ) from the left, using the SBPproperty of Q x,y and the same penalty matrices as in (), one obtains ( e 2 ) t e T (E P y A + s )e + e T (E N P y A s )e e T (P x E B + s )e + e T (P x E N B s )e + 2e T (P x P y I 4 )T e. (4) In [2], it was shown that equation (4) can be rewritten as e t et (P en)e 2 e 2 e + T e (5) where the matrix P en = P en T denotes the sum of all penalty matrices and boundary terms in (4). By assuming that η = et (P en)e < η 2 e 2 < for some constant η, (5) leads to e e η t t e ητ T e dτ T e max. (6) η The error e is thus bounded by the truncation error according to (6), and the error (3) is therefore bounded in time. It can be shown, see [2], that even the case where η is occasionally zero leads to an error bound. We conclude these results in the following proposition Proposition 2. The semi-discrete formulation (9) with the boundary operators given by (5) is error bounded in time if Σ e,w = α,2 I x P y I 4, Σ s,n = α 3,2 P x I y I 4 where α,3 = and α 2,4 =. 9
12 3. The multiple penalty technique Next, we move on to the MPT. We consider a local domain where the boundary data has been incorporated from global results. The MPT is applied at several grid points close to the boundaries where the applied data has been interpolated from the global mesh. This is similar to the method introduced Davies []. Here, the data is applied weakly with the SAT technique such that stability is preserved. This technique has been studied earlier in [8]. The technique will be demonstrated close to the eastern boundary; the analysis is similar at the other boundaries, and therefore neglected. Equation (3) with the MPT on the eastern boundary is e t + (Px Q x I y A s )e + (I x Py Q y B s )e = (P en)e + P x P y I 4 N p i= ϕ N i ((E N i Σ N i e L s e)e (7) where ϕ N i > are arbitrary chosen coefficients, P en is the standard penalty matrices and E N i a projection matrix with the entry N i, N i equal to one and zero otherwise and Σ N i e are the additional penalty matrices, where the superscript indicate at which grid points the penalties are applied. Note that the truncation error is neglected in (7). By multiplying (7) with e T (P x P y I 4 ), one arrives at e 2 t i e T ϕ i (P x E N i (Σ N i e + (Σ N i ) T e ) L s e)e (8) since the rest of the terms only contributes with a decay of energy, as concluded in the previous section. According to (5), the boundary operator L s e is symmetric and negative semi-definite, so the norm in (8) will remain bounded if the penalty matrices have a positive semi-definite symmetric part, Σ N i e + (Σ N i ) T e. The general form of the stability condition is (Σ N i e,w,s,n + (Σ N i ) T e,w,s,n) L s e,w,s,n. We summarize the results in the following proposition Proposition 3. The MPT applied close to the boundaries as in (7) with penalty matrices Σ N i e,w,s,n preserves stability if (Σ N i e,w,s,n + (Σ N i ) T e,w,s,n) L s e,w,s,n.
13 4. Numerical results for the steady-state problem In our test-cases we have chosen the reference state u = v =, c = 2, ρ = and introduced a disturbance in form of a Gaussian pressure pulse centered in the middle of the local domain, that is p (x, y, ) = e (x /2)2 (y /2) 2 ; see Figure 2 for an illustration. This local domain is contained in a coarse, global mesh, as illustrated in Figure. Since the pulse illustrate a disturbance or an error, we want it to disappear as fast as possible. The multiple penalties are applied according to the procedure described in section 3, such that stability is preserved. A square, uniform mesh is used and the number of multiple penalties close to each boundary is N MP = rn, where r is a number between zero and one and N is the number of grid points. The solution, u, is assumed to be known in these domains. For the numerical approximation, a SBP scheme with third order overall accuracy is used and a fourth order Runge-Kutta method is used to integrate in time. Figure 2 shows the pressure in the local domain at T =, which is advected out of the domain when integrated in time. The pulse will eventually be reflected at the boundaries, and parts of the error remain in the domain of interest. Figure 3 shows the error at T = 2. when standard penalties are used, and Figure 4 how the amplitude of the error is reduced by using multiple penalties. In Figures 5 and 6, we observe how the error decreases as one integrate in time. This is clarified by displaying the time-dependence of the error in Figures 7-9. The increased rate of convergence to steady-state when applying multiple penalties is consistent with previous work in [8]. 4.. Time-dependent solutions We will in this section demonstrate how the numerical errors are reduced also for time-dependent boundary data when multiple penalties are applied.
14 Figure 2: Initial condition in the local domain for the linearized pressure variable inside the local domain. In our example, we consider a uniform mesh with N = Figure 3: The error of the linearized pressure variable at T = 2. when standard penalties are used. N = 2 2
15 Figure 4: The error for the linearized pressure variable at T = 2. when multiple penalties (r =.2) are used. N = 4 8. x Figure 5: The error of the linearized pressure variable at T = 8. when standard penalties are used. Note that the z-axis has been re-scaled when comparing with Figure 3. N = 2 3
16 8. x Figure 6: The error for the linearized pressure variable at T = 8. when multiple penalties (r =.2) are used. Note that the z-axis has been re-scaled when comparing with Figure 4. N = 4 r= r=.2 r=.4 Norm of error t Figure 7: Error for the quantity p as a function of time. The number of grid points is N = 2 in each space direction. r = N MP /N 4
17 r= r=.2 r=.4 2 Norm of error t Figure 8: Error for the quantity p as a function of time. The number of grid points is N = 4 in each space direction. r = N MP /N r= r=.2 r=.4 2 Norm of error t Figure 9: Error for the quantity p as a function of time. The number of grid points is N = 8 in each space direction. r = N Mp /N 5
18 We consider a slightly modified version of (2) w t + A s (U )w x + B s (U )w y = F (x, y, t) L s ww (, y, t) = g w(y, t) L s sw (x,, t) = g e(x, t) L s ew (, y, t) = g s(y, t) L s nw (x,, t) = g n(x, t) w (x, y, ) = f (x, y) (9) where F (x, y, t) is a known forcing function. We create a manufactured analytic solution given by u(x, y, t) = exp( r 2 )sin(2π(r 2 t)), (2) by choosing F = u t + A s u x + B s u y and the boundary data accordingly. In (2), r 2 = (x.5) 2 + (y.5) 2 and the function describes an oscillating pulse centered at (.5,.5), as illustrated in Figure. We then start from the initial condition f (x, y) = and measure how fast the numerical solution converges to (2). In Figure -3, the difference between the numerical solution and the manufactured one is measured as a function of time. Since a (2) cannot, unlike the trivial solution in the previous section, be approximated exactly by the SBP-operators, a small error will always remain. However, by adding multiple penalties, the error will drop to the minimal level faster, as illustrated in Figure Conclusion Using the SBP-SAT technique, a lateral boundary procedure similar to the one introduced by H.C Davies [] has been applied to the linearized Euler equations in two space dimensions. The technique is straight-forward to implement and stability is easily shown by energy methods. These additional penalty terms contributes with a decay of energy and increases the rate of convergence to steady-state. The lateral boundaries also works as buffer zones, where spurious reflections are reduced. 6
19 Figure : The time-dependent solution of p at T =.7. N = 2 in each space-direction. r= r=.2 r=.4 Norm of error t Figure : Error for the quantity p as a function of time for time-dependent boundary data. N = 2 and r = N MP /N. 7
20 r= r=.2 r=.4 Norm of error t Figure 2: Error for the quantity p as a function of time for time-dependent boundary data. N = 4 and r = N MP /N. r= r=.2 r=.4 Norm of error t Figure 3: Error for the quantity p as a function of time for time-dependent boundary data. N = 8 and r = N MP /N. 8
21 The effect is in our examples large; the errors can in some cases be reduced by several orders of magnitude, both for the steady-state problem and with time-dependent boundary data. In the cases studied in this work, all data from the global mesh is considered given with arbitrary accuracy, and the domains are assumed to be uncoupled. References [] H.C Davies. A lateral boundary formulation for multi-level prediction models, Quart.J.Roy. Meteor. Soc., 2: 45-48, 976. [2] G.E Birchfield. Numerical prediction of hurricane movement with the use of a fine grid, Journal of Meteorology, 7: 46-44, 96. [3] G.E Hill. Grid telescoping in Numerical weather prediction, Journal of Applied Meteorology, 29-38, 967. [4] E.J Harrison, R.L Elsberry. A method for incorporating nested finite grids in the solution of systems of geophysical equations, Journal of the Atmospheric Sciences, 29: , 972 [5] D.L Williamson, G.L Browning. Formulation of the lateral boundary conditions for the NCAR limited-area model, Journal of Applied Meteorology, 3: 8-6, 973 [6] J.P Boyd. Limited-area Fourier spectral methods and data analysis schemes: windows, Fourier extension, Davies relaxation, and all that. Monthly Weather Review, 33: , 25 [7] B.Denis, J.Côté, R.Laprise. Spectral Decomposition of Two-Dimensional Atmospheric Fields on Limited-Area Domains Using Discrete Cosine Transform (DCT). Monthly Weather Review,3: , 22 [8] D.J Perkey, C.W Kreitzberg. A time-dependent lateral boundary scheme for limited area primitive equation models. Monthly Weather Review, 4: , 975. [9] J. Nordström, J. Gong. A stable hybrid method for hyperbolic problems. Journal of Computational Physics 22: , 26. 9
22 [] J. Gong, J. Nordström. A stable end efficient hybrid scheme for for viscous problems in complex geometries. Journal of Computational Physics 226: 29-39, 27. [] J. Nordström, F. Ham, M. Shoeybi, E. Weide, M. Svärd, K. Mattsson, G. Iaccarino, J. Gong. A hybrid method for unsteady inviscid fluid flow. Computers and Fluids 38: , 29. [2] J. Nordström, J. Gong, E. Weide, M. Svärd. A stable and conservative high order multiblock method for the compressible Navier- Stokes equations. Journal of Computational Physics 228: , 29. [3] M. Svärd, J. Nordström. Review of summation-by-parts schemes for initial-boundary-value problems. Journal of Computational Physics, 48: , 24 [4] B. Strand. Summation by parts for finite difference approximations for, Journal of Computational Physics : (994) d dx [5] K. Mattsson, J. Nordström. Summation by parts for finite difference approximations of second order derivatives, Journal of Computational Physics, 48: (999) [6] M. Carpenter, D. Gottlieb, S. Aberbanel. Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes, Journal of Computational Physics : , 994. [7] M. Carpenter,J. Nordström, D. Gottlieb. A stable and conservative interface treatment of arbitrary spatial accuracy, Journal of Computational Physics 48: , 999. [8] J. Nordström, Q. Abbas, B.A Erickson and H. Frenander. A flexible boundary procedure for hyperbolic problems: multiple penalty terms. Communications in Computational Physics 6(2): 54-57, 24. [9] C.W Rowley, T. Colonius. Discretely non-reflecting boundary conditions for linear hyperbolic systems, Journal of Computational Physics 57: 5-538, 2. 2
23 [2] J. Nordström, M. Svärd. Well-posed boundary conditions for the Navier- Stokes equations. SIAM Journal on Numerical Analysis 43: , 25, [2] B. Erickson, J. Nordström. Stable, high order accurate adaptive schemes for long time, highly intermittent geophysics problems. Journal of Computational and Applied Mathematics 27: , 24, 2
A Stable and Accurate Davies-like Relaxation Procedure using Multiple Penalty Terms for Lateral Boundary Conditions
A Stable and Accurate Davies-like Relaxation Procedure using Multiple Penalty Terms for Lateral Boundary Conditions Hannes Frenander Division of Computational Mathematics, Department of Mathematics, Linköping
More informationWell-posedness, stability and conservation for a discontinuous interface problem: an initial investigation.
Well-posedness, stability and conservation for a discontinuous interface problem: an initial investigation. Cristina La Cognata and Jan Nordström Abstract A robust interface treatment for the discontinuous
More informationSpectral analysis of the incompressible Navier-Stokes equations with different boundary conditions
Spectral analysis of the incompressible Navier-Stokes equations with different boundary conditions Cristina La Cognata, Jan Nordström Department of Mathematics, Computational Mathematics, Linköping University,
More informationNew Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5-94 New Diagonal-Norm Summation-by-Parts Operators for the First Derivative with Increased Order of Accuracy David
More informationEfficient wave propagation on complex domains
Center for Turbulence Research Annual Research Briefs 2006 223 Efficient wave propagation on complex domains By K. Mattsson, F. Ham AND G. Iaccarino 1. Motivation and objectives In many applications, such
More informationWell-posedness, Stability and Conservation for a Discontinuous Interface Problem
Department of Mathematics Well-posedness, Stability and Conservation for a Discontinuous Interface Problem Cristina La Cognata and Jan Nordström LiTH-MAT-R--214/16--SE Department of Mathematics Linköping
More informationREVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS
REVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS MAGNUS SVÄRD AND JAN NORDSTRÖM Abstract. High-order finite difference methods are efficient, easy to program, scale well in multiple
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
More informationStable and high-order accurate finite difference schemes on singular grids
Center for Turbulence Research Annual Research Briefs 006 197 Stable and high-order accurate finite difference schemes on singular grids By M. Svärd AND E. van der Weide 1. Motivation and objectives The
More informationConjugate Heat Transfer for the Unsteady Compressible Navier-Stokes Equations Using a Multi-block Coupling
Conjugate Heat Transfer for the Unsteady Compressible Navier-Stokes Equations Using a Multi-block Coupling Jan Nordström Department of Mathematics, Linköping University, SE-58 83 Linköping, Sweden Jens
More informationCharacterizing the Accuracy of Summation-by-Parts Operators for Second-Derivatives with Variable-Coefficients
1 Characterizing the Accuracy of Summation-by-Parts Operators for Second-Derivatives with Variable-Coefficients by Guang Wei Yu Supervisor: D. W. Zingg April 13, 2013 Abstract This paper presents the
More informationSUMMATION-BY-PARTS IN TIME: THE SECOND DERIVATIVE
Department of Mathematics SUMMATION-BY-PARTS IN TIME: THE SECOND DERIVATIVE Jan Nordström and Tomas Lundquist LiTH-MAT-R--2014/11--SE Department of Mathematics Linköping University S-581 83 Linköping SUMMATION-BY-PARTS
More informationNUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY
NUMERICAL SOLUTION OF THE LINEARIZED EULER EQUATIONS USING HIGH ORDER FINITE DIFFERENCE OPERATORS WITH THE SUMMATION BY PARTS PROPERTY Stefan Johansson Department of Information Technology Scientific Computing
More informationDiagonal-norm upwind SBP operators
Diagonal-norm upwind SBP operators Ken Mattsson June 8, 16 Abstract High-order accurate first derivative finite difference operators are derived that naturally introduce artificial dissipation. The boundary
More informationBoundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations
Boundary Conditions for a Divergence Free Velocity-Pressure Formulation of the Incompressible Navier-Stokes Equations Jan Nordström, Ken Mattsson and Charles Swanson May 5, 6 Abstract New sets of boundary
More informationSingularity of the discrete Laplacian operator
Singularity of the discrete Laplacian operator Andrea Alessandro Ruggiu Jan Nordström Singularity of the Laplacian A.A.Ruggiu, J.Nordström August 24, 2015 1 Motivation In the SBP SAT framework, for unsteady
More informationNew developments for increased performance of the SBP-SAT finite difference technique
New developments for increased performance of the SBP-SAT finite difference technique Jan Nordström 1 and Peter Eliasson 2 1 Department of Mathematics, Computational Mathematics, University of Linköping,
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 informationEnergy decay in vortices
UPTEC F11 031 Examensarbete 30 hp 07 Juni 2011 Energy decay in vortices Björn Lönn Abstract Energy decay in vortices Björn Lönn Teknisk- naturvetenskaplig fakultet UTH-enheten Besöksadress: Ångströmlaboratoriet
More informationWell Posed, Stable and Weakly Coupled Fluid Structure Interaction Problems
Well Posed, Stable and Weakly Coupled Fluid Structure Interaction Problems Jan Nordström and Sofia Eriksson April 26, 29 Abstract We investigate problems of fluid structure interaction type and aim for
More informationAccurate and stable finite volume operators for unstructured flow solvers
Center for Turbulence Research Annual Research Briefs 2006 243 Accurate and stable finite volume operators for unstructured flow solvers By F. Ham, K. Mattsson AND G. Iaccarino 1. Motivation and objectives
More informationEULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS
EULERIAN DERIVATIONS OF NON-INERTIAL NAVIER-STOKES EQUATIONS ML Combrinck, LN Dala Flamengro, a div of Armscor SOC Ltd & University of Pretoria, Council of Scientific and Industrial Research & University
More informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More informationOptimal diagonal-norm SBP operators
Optimal diagonal-norm SBP operators Ken Mattsson 1, Martin Almquist 1 and Mark H. Carpenter 2 1 Department of Information Technology, Uppsala University 2 Computational Aerosciences Branch, NASA Langley
More informationUNIVERSITY OF EAST ANGLIA
UNIVERSITY OF EAST ANGLIA School of Mathematics May/June UG Examination 2007 2008 FLUIDS DYNAMICS WITH ADVANCED TOPICS Time allowed: 3 hours Attempt question ONE and FOUR other questions. Candidates must
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationEnergy Stable Boundary Conditions for the Nonlinear Incompressible Navier-Stokes Equations
Energy Stable Boundary Conditions for the Nonlinear Incompressible Navier-Stokes Equations Jan Nordström, Cristina La Cognata Department of Mathematics, Computational Mathematics, Linköping University,
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationAbstract Simulations of sound pressure levels are often used as a decision basis when planning location for e.g wind turbines. In previous work, an e
Abstract Simulations of sound pressure levels are often used as a decision basis when planning location for e.g wind turbines. In previous work, an e cient model for computing sound propagation over irregular
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationStable Numerical Scheme for the Magnetic Induction Equation with Hall Effect
Stable Numerical Scheme for the Magnetic Induction Equation with Hall Effect Paolo Corti joint work with Siddhartha Mishra ETH Zurich, Seminar for Applied Mathematics 17-19th August 2011, Pro*Doc Retreat,
More informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationWell-Balanced Schemes for the Euler Equations with Gravity
Well-Balanced Schemes for the Euler Equations with Gravity Alina Chertock North Carolina State University chertock@math.ncsu.edu joint work with S. Cui, A. Kurganov, S.N. Özcan and E. Tadmor supported
More informationHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
High Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation Faheem Ahmed, Fareed Ahmed, Yongheng Guo, Yong Yang Abstract This paper deals with
More informationMultiscale method and pseudospectral simulations for linear viscoelastic incompressible flows
Interaction and Multiscale Mechanics, Vol. 5, No. 1 (2012) 27-40 27 Multiscale method and pseudospectral simulations for linear viscoelastic incompressible flows Ling Zhang and Jie Ouyang* Department of
More informationMAE210C: Fluid Mechanics III Spring Quarter sgls/mae210c 2013/ Solution II
MAE210C: Fluid Mechanics III Spring Quarter 2013 http://web.eng.ucsd.edu/ sgls/mae210c 2013/ Solution II D 4.1 The equations are exactly the same as before, with the difference that the pressure in the
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationA new well-posed vorticity divergence formulation of the shallow water equations
See discussions, stats, and author profiles for this publication at: http://www.researchgate.net/publication/280312608 A new well-posed vorticity divergence formulation of the shallow water equations ARTICLE
More informationCFD with OpenSource software
CFD with OpenSource software The Harmonic Balance Method in foam-extend Gregor Cvijetić Supervisor: Hrvoje Jasak Faculty of Mechanical Engineering and Naval Architecture University of Zagreb gregor.cvijetic@fsb.hr
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationSubsonic characteristic outflow boundary conditions with transverse terms: a comparative study
Center for Turbulence Research Annual Research Briefs 200 35 Subsonic characteristic outflow boundary conditions with transverse terms: a comparative study By G. Lodato, F. Ham AND H. Pitsch. Motivation
More informationMULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER EQUATION. Thomas Y. Hou. Danping Yang. Hongyu Ran
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 13, Number 5, December 2005 pp. 1153 1186 MULTISCALE ANALYSIS IN LAGRANGIAN FORMULATION FOR THE 2-D INCOMPRESSIBLE EULER
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationOn the Convergence Rates of Energy- Stable Finite-Difference Schemes
Department of Mathematics On the Convergence Rates of Energy- Stable Finite-Difference Schemes Magnus Svärd and Jan Nordström LiTH-MAT-R--217/14--SE Department of Mathematics Linköping University S-581
More informationCoupled High-Order Finite Difference and Unstructured Finite Volume Methods for Earthquake Rupture Dynamics in Complex Geometries.
UPTC F11040 xamensarbete 30 hp Juni 011 Coupled High-Order Finite Difference and Unstructured Finite Volume Methods for arthquake Rupture Dynamics in Complex Geometries Ossian OReilly Abstract Coupled
More informationOpportunities for efficient high-order methods based on the summation-by-parts property
AIAA Aviation -6 June 5, Dallas, TX nd AIAA Computational Fluid Dynamics Conference AIAA 5- Opportunities for efficient high-order methods based on the summation-by-parts property Jason E. Hicken, a David
More informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY ON INTERIOR NODAL SETS
6th European Conference on Computational Mechanics (ECCM 6 7th European Conference on Computational Fluid Dynamics (ECFD 7 11 15 June 018, Glasgow, UK ON THE BENEFIT OF THE SUMMATION-BY-PARTS PROPERTY
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationDecay profiles of a linear artificial viscosity system
Decay profiles of a linear artificial viscosity system Gene Wayne, Ryan Goh and Roland Welter Boston University rwelter@bu.edu July 2, 2018 This research was supported by funding from the NSF. Roland Welter
More informationarxiv: v1 [physics.comp-ph] 28 Jan 2008
A new method for studying the vibration of non-homogeneous membranes arxiv:0801.4369v1 [physics.comp-ph] 28 Jan 2008 Abstract Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo
More informationarxiv: v1 [math.ap] 2 Feb 2011
HIGHER ORDER FINITE DIFFERENCE SCHEMES FOR THE MAGNETIC INDUCTION EQUATIONS WITH RESISTIVITY U KOLEY, S MISHRA, N H RISEBRO, AND M SVÄRD arxiv:48v [mathap] Feb Abstract In this paper, we design high order
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationDavid Abrecht. February 17, 2012
Addendum to An analytical solution method for the unsteady, unbounded, incompressible three-dimensional Navier-Stokes equations in Cartesian coordinates using coordinate axis symmetry degeneracy David
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
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 informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationStable 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 informationCompressible Navier-Stokes (Euler) Solver based on Deal.II Library
Compressible Navier-Stokes (Euler) Solver based on Deal.II Library Lei Qiao Northwestern Polytechnical University Xi an, China Texas A&M University College Station, Texas Fifth deal.ii Users and Developers
More informationarxiv: v1 [math.na] 20 Nov 2013
REVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS MAGNUS SVÄRD AND JAN NORDSTRÖM arxiv:1311.4984v1 [math.na] 2 Nov 213 Abstract. High-order finite difference methods are efficient,
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationThe Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing
The Dynamics of Detonation with WENO and Navier-Stokes Shock-Capturing Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and
More informationVerified Calculation of Nonlinear Dynamics of Viscous Detonation
Verified Calculation of Nonlinear Dynamics of Viscous Detonation Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and Joseph
More informationρ Du i Dt = p x i together with the continuity equation = 0, x i
1 DIMENSIONAL ANALYSIS AND SCALING Observation 1: Consider the flow past a sphere: U a y x ρ, µ Figure 1: Flow past a sphere. Far away from the sphere of radius a, the fluid has a uniform velocity, u =
More informationFinite Volume Method
Finite Volume Method An Introduction Praveen. C CTFD Division National Aerospace Laboratories Bangalore 560 037 email: praveen@cfdlab.net April 7, 2006 Praveen. C (CTFD, NAL) FVM CMMACS 1 / 65 Outline
More informationMultiscale Computation of Isotropic Homogeneous Turbulent Flow
Multiscale Computation of Isotropic Homogeneous Turbulent Flow Tom Hou, Danping Yang, and Hongyu Ran Abstract. In this article we perform a systematic multi-scale analysis and computation for incompressible
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationImprovement of convergence to steady state solutions of Euler equations with. the WENO schemes. Abstract
Improvement of convergence to steady state solutions of Euler equations with the WENO schemes Shuhai Zhang, Shufen Jiang and Chi-Wang Shu 3 Abstract The convergence to steady state solutions of the Euler
More informationJACOBI S ITERATION METHOD
ITERATION METHODS These are methods which compute a sequence of progressively accurate iterates to approximate the solution of Ax = b. We need such methods for solving many large linear systems. Sometimes
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationPublications by Jan Nordström updated:
Publications by Jan Nordström updated: 2018-03-01 5 most cited publications (Google Scholar, Scopus, Web of Science) 1. M. H. Carpenter, J. Nordström & D. Gottlieb, A Stable and Conservative Interface
More informationExtension to moving grids
Extension to moving grids P. Lafon 1, F. Crouzet 2 & F. Daude 1 1 LaMSID - UMR EDF/CNRS 2832 2 EDF R&D, AMA April 3, 2008 1 Governing equations Physical coordinates Generalized coordinates Geometrical
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationEfficient outdoor sound propagation modelling in time-domain
Efficient outdoor sound propagation modelling in time-domain Timothy Van Renterghem Acoustics Group Department of Information Technology Ghent University/iMinds CEAS, X-Noise : Atm. and ground effects
More information3. FORMS OF GOVERNING EQUATIONS IN CFD
3. FORMS OF GOVERNING EQUATIONS IN CFD 3.1. Governing and model equations in CFD Fluid flows are governed by the Navier-Stokes equations (N-S), which simpler, inviscid, form is the Euler equations. For
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationSpace-time Discontinuous Galerkin Methods for Compressible Flows
Space-time Discontinuous Galerkin Methods for Compressible Flows Jaap van der Vegt Numerical Analysis and Computational Mechanics Group Department of Applied Mathematics University of Twente Joint Work
More informationSPECTRAL ELEMENT STABILITY ANALYSIS OF VORTICAL FLOWS
SPECTRAL ELEMENT STABILITY ANALYSIS OF VORTICAL FLOWS Michael S. Broadhurst 1, Vassilios Theofilis 2 and Spencer J. Sherwin 1 1 Department of Aeronautics, Imperial College London, UK; 2 School of Aeronautics,
More informationNumerical Solutions of the Burgers System in Two Dimensions under Varied Initial and Boundary Conditions
Applied Mathematical Sciences, Vol. 6, 22, no. 3, 563-565 Numerical Solutions of the Burgers System in Two Dimensions under Varied Initial and Boundary Conditions M. C. Kweyu, W. A. Manyonge 2, A. Koross
More informationImplementation of a symmetry-preserving discretization in Gerris
Implementation of a symmetry-preserving discretization in Gerris Daniel Fuster Cols: Pierre Sagaut, Stephane Popinet Université Pierre et Marie Curie, Institut Jean Le Rond D Alembert Introduction 10/11:
More informationSIMPLE Algorithm for Two-Dimensional Channel Flow. Fluid Flow and Heat Transfer
SIMPLE Algorithm for Two-Dimensional Channel Flow Fluid Flow and Heat Transfer by Professor Jung-Yang San Mechanical Engineering Department National Chung Hsing University Two-dimensional, transient, incompressible
More informationStable boundary treatment for the wave equation on second-order form
Stable boundary treatment for the wave equation on second-order form Ken Mattsson Frank Ham Gianluca Iaccarino June 24, 2008 Abstract A stable and accurate boundary treatment is derived for the secondorder
More informationFinal Exam May 4, 2016
1 Math 425 / AMCS 525 Dr. DeTurck Final Exam May 4, 2016 You may use your book and notes on this exam. Show your work in the exam book. Work only the problems that correspond to the section that you prepared.
More information