arxiv: v1 [math.na] 2 Feb 2018
|
|
- Barrie Edwards
- 5 years ago
- Views:
Transcription
1 A Kernel Based High Order Eplicit Unconditionally Stable Scheme for Time Dependent Hamilton-Jacobi Equations Andrew Christlieb, Wei Guo, and Yan Jiang 3 Abstract. arxiv:8.536v [math.na] Feb 8 In this paper, a class of high order numerical schemes is proposed for solving Hamilton-Jacobi (H-J) equations. This work is regarded as an etension of our previous work for nonlinear degenerate parabolic equations, see Christlieb et al. [3], which relies on a special kernel-based formulation of the solutions and successive convolution. When applied to the H-J equations, the newly proposed scheme attains genuinely high order accuracy in both space and time, and more importantly, it is unconditionally stable, hence allowing for much larger time step evolution compared with other eplicit schemes and saving computational cost. A high order weighted essentially non-oscillatory methodology and a novel nonlinear filter are further incorporated to capture the correct viscosity solution. Furthermore, by coupling the recently proposed inverse La-Wendroff boundary treatment technique, this method is very fleible in handing comple geometry as well as general boundary conditions. We perform numerical eperiments on a collection of numerical eamples, including H-J equations with linear, nonlinear, conve or non-conve Hamiltonians. The efficacy and efficiency of the proposed scheme in approimating the viscosity solution of general H-J equations is verified. Key Words: Hamilton-Jacobi equation; Kernel based scheme; Unconditionally stable; High order accuracy; Weighted essentially non-oscillatory methodology; Viscosity solution. Department of Computational Mathematics, Science and Engineering, Department of Mathematics and Department of Electrical Engineering, Michigan State University, East Lansing, MI, christli@.msu.edu Department of Mathematics and Statistics, Teas Tech University, Lubbock, TX, 749. Research is supported by NSF grant NSF-DMS weimath.guo@ttu.edu 3 Department of Mathematics, Michigan State University, East Lansing, MI, jiangyan@math.msu.edu
2 Introduction In this paper, we propose a class of high order non-oscillatory numerical schemes for approimating the viscosity solution to the Hamilton-Jacobi (H-J) equation { φt + H( φ,, t) =, Ω R d, φ(, ) = φ (), (.) with suitable boundary conditions. The H-J equations find applications in diverse fields, including optimal control, seismic waves, crystal growth, robotic navigation, image processing, and calculus of variations. The well-known property of the H-J equations is that the weak solutions may not be unique, and the concept of viscosity solutions was introduced to single out the unique physically relevant solution [5, 4]. The viscosity solution to the H-J equation (.) is known to be Lipschitz continuous but may develop discontinuous derivatives in finite time, even with a C initial condition. A variety of numerical schemes have been proposed to solve the H-J equation (.) in the literature, including for eample, the monotone schemes [6,, 6], the essentially non-oscillatory (ENO) schemes [8, 9] or the weighted ENO (WENO) schemes [4, 4], the Hermite WENO (HWENO) schemes [3, 3, 39, 43, 44], the discontinuous Galerkin methods [, 7, 9, 4, ]. More details about recent development of high-order numerical methods for solving the H-J equations can be found in the review paper [35]. Most of these methods are in the Method of Lines (MOL) framework, which means that the discretization is first carried out for the spatial variable, then the numerical solution is updated in time by coupling a suitable time integrator. The most commonly used time evolution methods are the strong-stability-preserving Runge-Kutta (SSP RK) schemes and SSP multi-step schemes [9, 34, 8], which can preserve the strong stability in some desired norm and prevent spurious oscillations near spatial discontinuities. However, it is well known that an eplicit time discretization does have a restriction on the time step in order to maintain linear stability. Recently, optimization algorithms based on the Hopf formula were developed for solving a class of H-J equations from optimal control and differential games [7, ]. Such algorithms are able to overcome the curse of dimensionality when solving (.) with large d and have attracted a lot of attention. However, the applicability of this type of algorithms depends on the Hopf formula, which is not available for general H-J equations, such as when the Hamiltonian H and the initial condition φ are both non-conve.
3 Another framework named the Method of Lines Transpose (MOL T ) has been eploited in the literature for solving time-dependent partial differential equations (PDEs). It is also known as Rothes method or transverse method of lines [3, 33, 5]. In such a framework for solving linear PDEs, the temporal variable is first discretized, resulting in a set of linear boundary value problems (BVPs) at discrete time levels. Each BVP can be inverted analytically in an integral formulation based on a Green s/kernel function and then the numerical solution is updated accordingly. As a notable advantage, the MOL T approach is able to use an implicit method but avoid solving linear systems at each time step, see [5]. Moreover, a fast convolution algorithm is developed to reduce the computational compleity of the scheme from O(N ) to O(N) [7,, ], where N is the number of discrete mesh points. Over the past several years, the MOL T methods have been developed for solving the heat equation [4, 6, 5, 3], Mawell s equations [8], the split Vlasov equation [], among many others. This methodology can be generalized to solving some nonlinear problems, such as the Cahn-Hilliard equation [4]. However, it rarely applied to general nonlinear problems, mainly because efficient fast algorithms of inverting nonlinear BVPs are lacking and hence the advantage of the MOL T is compromised. More recently, the authors proposed a novel numerical scheme for solving the nonlinear degenerate parabolic equations [3]. Following the MOL T philosophy, we epress the spatial derivatives in terms of a kernel-based representation. By construction, the solution is evolved by an eplicit SSP RK method, similar to the MOL approach. However, by carefully choosing a parameter β introduced in the formulation, the scheme is proven to be unconditionally stable, and hence allowing for large time step evolution and saving computational cost. This is considered the most remarkable property of the proposed method. The robust WENO methodology and a nonlinear filter are further employed to capture sharp gradients of the solution without producing spurious oscillations, and at the same time, the high order accuracy is attained in smooth regions. The scheme can be easily implemented in a nonuniform but orthogonal mesh and hence very effective in handling comple geometry. In this paper, we etend the work in [3] to solve H-J equations (.). In particular, under such a framework, the spatial derivatives are still reconstructed by the kernel-based approach, which is genuinely high order accurate in both space and time and essentially non-oscillatory due to the use of the WENO methodology as well as the nonlinear filter. This is very desired since the viscosity 3
4 solution of the H-J equations may develop discontinuous derivatives. More importantly, the scheme is unconditionally stable if the parameter β is appropriately chosen. By further incorporating the recently proposed high order inverse La-Wendroff (ILW) boundary treatment technique, the scheme is able to deal with many types of boundary conditions effectively such as non-homogeneous Dirichlet and Neumann boundary conditions. In summary, the proposed scheme for solving the H-J equations is robust, high order accurate, unconditionally stable, fleible and efficient. The rest of this paper is organized as follows. In Section, we present our numerical scheme for one-dimensional (D) H-J equations. The two-dimensional (D) case is considered in Section 3. In Section 4, a collection of numerical eamples is presented to demonstrate the performance of the proposed method. In Section 5, we conclude the paper with some remarks and future work. One-dimensional case In the D case, (.) becomes φ t + H(φ ) =, φ(, ) = φ (). (.) Assume the domain is a closed interval [a, b] which is partitioned with N + points a = < < < N < N = b, (.) with i = i i. Note that the mesh can be nonuniform. Let φ i (t) denote the solution φ( i, t) at mesh point i. As with many H-J solvers, we will construct the following semi-discrete scheme d dt φ i(t) + Ĥ(φ,i, φ+,i ) =, i =,..., N, (.3) where Ĥ is a numerical Hamiltonian. In this work, we use the local La-Friedrichs flu Ĥ(u, u + ) = H( u + u + ) α(u, u + ) u+ u, (.4) where α(u, u + ) = ma u H (u) with the maimum taken over the range bounded by u and u +. φ,i and φ+,i in (.3) are the approimations to φ at i obtained by left-biased and right-biased methods, respectively, to take into the account the direction of characteristics propagation of the H-J equation. This is the main part of this paper and will be detailed in the following subsections. 4
5 . Approimation of the first order derivative We start with a brief review on construction of the kernel-based representation of proposed in [3]. First, we introduce two operators L L and L R L L = I +, L R = I, [a, b], (.5) where I is the identity operator and > is a constant. Then we can invert L L and L R analytically as follows: L L [v, ]() = IL [v, ]() + A L e ( a), L R [v, ]() = IR [v, ]() + B R e (b ), (.6a) (.6b) where I L [v, ]() = I R [v, ]() = a b e ( y) v(y)dy, e (y ) v(y)dy, (.7a) (.7b) with constant A L and B R being determined by the boundary condition imposed for the operators. For eample, if assume L L and L R are periodic, i.e., L L [v, ](a) = L L [v, ](b), and L R [v, ](a) = L [v, ](b), R then we have A L = IL [v, ](b) µ, and B R = IR [v, ](a), (.8) µ with µ = e (b a). In addition, it is straightforward to show that D p L. p= I) = D R/(I D R ) = D p R, p= = L L I = L L (I L L ) = D L/(I D L ) = = I L R = L R (L R (.9a) (.9b) where the two operators D L and D R are defined by D L = I L L, D R = I L R, [a, b]. (.) 5
6 Hence, following the idea in [3], when φ is a periodic function, we can approimate the first derivative φ ± with (modified) partial sums in (.9), φ () Pk L [φ, ]() = k D p L [φ, ](), k =,, p= k (.a) D p L [φ, ]() D DL [φ, ](), k = 3, p= φ + () Pk R [φ, ]() = k D p R [φ, ](), k =,, p= k (.b) D p R [φ, ]() + D DR [φ, ](), k = 3. p= Note that by construction, D p L [φ, ]() only depends on the values of φ on [a, ], while Dp R [φ, ]() only depends on the values of φ on [, b]. This design in fact takes into account the direction of propagating characteristics of the H-J equation. Also note that there is an etra term for k = 3. As remarked in [3], such a term is needed for attaining unconditional stability of the scheme. D is defined as D [v, α]() = v() b a e y v(y)dy A e ( a) B e (b ). (.) The coefficients A and B are also obtained from the boundary condition. For instance, if we require D [v, α]() to be a periodic function, i.e., D [v, ](a) = D [v, ](b), and D [v, ](a) = D [v, ](b), then we have A = I [v, ](b) µ, B = I [v, ](a) µ (.3) with I [v, ]() = / b a e y v(y)dy. It is natural to require the following boundary conditions for the operators D p L (a) = Dp L (b), Dp L (a) = Dp L (b), D (a) = D (b), (.4) where p, if periodic boundary conditions of the solution are imposed. Moreover, an error estimate for the partial sums approimation (.) carried out in [3] shows that ã k φ() Pk L [φ, ]() C k+ φ(), (.5a) 6
7 ã k φ() Pk R [φ, ]() C k+ φ(). (.5b) In numerical simulations, we will take = β/(α t) in (.), with α being the maimum wave propagation speed. Here, t denotes the time step and β is a constant independent of t. Hence, the partial sums approimate φ with accuracy O( t k ). For time integration, we propose to use the classical eplicit SSP RK methods [9] to advance the solution from time t n to t n+. We denote φ n as the semi-discrete solution at time t n. In this paper, we use the following SSP RK methods, including the first order forward Euler scheme the second order SSP RK scheme φ n+ = φ n tĥ(φn,, φ n,+ ); (.6) φ () = φ n tĥ(φn, φ n+ = φn + and the third order SSP RK scheme φ () = u n tĥ(φn, φ () = 3 4 φn + 4 φ n+ = 3 φn + 3, φ n,+ ), Ä φ () tĥ(φ(),, φ (),+ ) ä ; (.7), φ n,+ ), Ä φ () tĥ(φ(),, φ (),+ ) ä, Ä φ () tĥ(φ(),, φ (),+ ) ä. (.8) In addition, linear stability of the proposed kernel-based schemes has been established in [3]. In particular, we proved the following theorem. Theorem.. For the linear equation φ t + cφ =, (i.e. the Hamiltonian is linear) with periodic boundary conditions, we consider the k th order SSP RK method as well as the k th partial sum in (.), with = β/( c t). Then there eists a constant β k,ma > for k =,, 3, such that the scheme is unconditionally stable provided < β β k,ma. The constants β k,ma for k =,, 3 are summarized in Table.. Table.: β k,ma in Theorem. for k =,, 3. k 3 β k,ma.43 7
8 . Non-periodic boundary conditions In this subsection, we will focus on the D H-J equation (.) with two commonly considered non-periodic boundary conditions, including Dirichlet and Neumann boundary conditions. Unlike periodic boundary conditions, it is not straightforward to specify the boundary conditions for the operators D p, where could be L, R or. The reason is that boundary conditions imposed on the operators D must be consistent with the boundary condition specified on φ (i.e. Dirichlet boundary condition) or its derivative (i.e. Neumann boundary condition). To address the difficulty, we propose to modify the partial sum (.). We will first investigate the operators D and prove a lemma in Subsection.., which is crucial for the scheme development as well as the error analysis. In Subsection.., under the assumption that we already have all boundary values as needed, e.g., the derivatives of φ at the boundary, we show that how to specify proper boundary conditions for the operators D, i.e., the coefficients A, B, A L and B R, such that the partial sum can approimate φ with high order accuracy and the Dirichlet or Neumann boundary condition is satisfied. In the Subsection..3, we will provide a systematic approach to reconstructing those needed boundary values... Operators with non-periodic boundary treatment We first study the operator D with a non-periodic boundary condition, where can be, L and R. Suppose C a and C b are given numbers. If we require D L [v, ](a) = C a, and D R [v, ](b) = C b, then, the coefficients are given as A L = v(a) C a, and B R = v(b) C b. (.9) If we require then, the coefficients are given as D [v, ](a) = C a, and D [v, ](b) = C b, A = Ä Ä µ I ä Ä µ [v, ](b) v(b) + C b I ää [v, ](a) v(a) + C a, (.a) B = Ä Ä µ I ä Ä µ [v, ](a) v(a) + C a I ää [v, ](b) v(b) + C b. (.b) 8
9 We start with proving the following lemma, which connects the function D [v, ]() with the derivatives of the function v() as well as their boundary values. Lemma.. Consider the boundary treatments (.9) and (.). Suppose v C k+ [a, b]. Then, we can obtain that D L [v, ]() = D R [v, ]() = p= p= ã p Ä p v() v(a)e p ( a)ä ã k+ I L [ k+ v, ]() + C a e ( a), (.a) ã p Ä p v() v(b)e p (b )ä ã k+ I R [ k+ v, ]() + C b e (b ). Suppose v C k+ [a, b]. Then, we have µc b C a µ e ( a) µc a C b µ e (b ) + µi [ k+ v, ](b) I [ k+ v, ](a) µ e ( a) + µi [ k+ (.b) ã p Ç D [v, ]() = p v() + µ p v(b) p v(a) µ p= e ( a) + µ p v(a) p å v(b) µ e (b ) ã k+ ÄI [ k+ v, ]() v, ](a) I [ k+ å v, ](b) µ e (b ) (.) Proof. For brevity, we only show the details of the proof for operator D L. Following a similar argument, one can easily prove the other cases. Using integration by parts repeatedly, we have I L [v, ]() =v() v(a)e ( a) IL [v (y), ]() = = p= ã p Ä p v() v(a)e p ( a)ä + ã k+ I L [ k+ v(y), ](). Therefore, with the coefficient A L = v(a) C a, D L [v, ]() can be rewritten as D L [v, ]() = p= ã p Ä p v() v(a)e p ( a)ä + ã k+ I L [ k+ v(y), ]() + C a e ( a). In light of Lemma., it is easy to check that given the boundary values φ (a) and φ (b) from the boundary condition, φ ± can be approimated by φ D L [φ, ](), φ + D R [φ, ](), (.3) 9
10 with order O(/) (i.e., first order accuracy since we choose = O(/ t)), by requiring D L (a) = φ (a), and D R (a) = φ (b). (.4) It seems that for the higher order scheme given in (.), the boundary condition for φ can be satisfied based on (.4) together with D p L [φ, ](a) = Dp R [φ, ](b) =, for p and D DL [φ, ](a) = D DR [φ, ](b) = if k = 3. However, due to the etra boundary terms appear in Lemma., one can check that the scheme (.) is not high order accurate. In fact, it is only first order in the case of non-periodic boundary conditions. To circumvent this problem, we will use a modified scheme instead of (.)... The modified partial sum Suppose we have obtained by some means the derivative values at boundary, i.e., m φ(a) and m φ(b), m. We propose the following modified partial sums for k 3 to handle non-periodic boundary conditions φ () P k L [φ, ]() = k D L [φ,p, ](), k =,, p= k D L [φ,p, ]() D [φ,3, ](), k = 3, p= (.5a) φ + () P k R [φ, ]() = And φ,p and φ,p are given as φ, = φ, φ, = D L [φ,, ] φ,3 = D L [φ,, ] + φ, = φ, φ, = D R [φ,, ] φ,3 = D R [φ,, ] + k D R [φ,p, ](), k =,, p= k D R [φ,p, ]() + D [φ,3, ](), k = 3. p= m= m= ã m m φ(a)e ( a), (m ) m= m= with the boundary conditions for the operators ã m m φ(a)e ( a), ã m m φ(b)e (b ), (m ) ã m m φ(b)e (b ), (.5b) (.6a) (.6b) D L [φ,, ](a) = φ (a), D R [φ,, ](b) = φ (b),
11 D L [φ,p, ](a) = D R [φ,p, ](b) =, for p, D [φ,3, ](a) = D [φ,3, ](b) =, could be or. One can easily check that the modified partial sum (.5) agrees with the derivative values at the boundary, i.e., P L k [φ, ](a) = φ (a), P R k [φ, ](b) = φ (b), meaning that the modified partial sum approimation (.5) is consistent with the boundary condition imposed on φ. Moreover, we have the following theorem: Theorem.3. Suppose φ C k+ [a, b], k =,, 3. Then, the modified partial sums (.5) satisfy Proof. With the help of Lemma., we have D L [φ,, ]() = φ ã φ() P k k L [φ, ]() C k+ φ(), (.7a) ã φ() P k k R [φ, ]() C k+ φ(). (.7b) p= ã p Ä p φ() φ(a)e p ( a)ä ã k+ I L [ k+ φ, ](). Thus, φ, () = p= ã p p φ() ã k+ I L [ k+ φ, ](). Using the fact the I L [v, ]() = v() I L [v, ](), we have the general form for q k Thus, we have Ñ ã q q φ, () = ã q q = = = q I q L [v, ]() = ( ) q p v() p + ( ) q I L [v, ](). q+k p=q+ q+k p=q+ p=q+ p= é ã p p φ() ã k+ I L [ k+ φ, ]() p= Ñ ã p p () ã k++q ã p p () + ã p p () k+q p=k+ é q ( ) q p p k+ φ() + ( ) q I L [ k+ φ, ]() p= ã p p φ() ã k+ I L [ k+ φ, ](). ã k+ I L [ k+ φ, ]()
12 Furthermore, we apply D L on φ, with D L (a) =, and obtain k D L [φ,, ]() = Therefore, = + q= Ñ k q= ã q Ä q φ, () φ q, (a)e ( a)ä ã k I L [ k φ,, ]() = (p ) p= + p=q+ ã p p φ() ã k+ I L [ k+ φ, ]() + ã k+ ÄI L ä [ k+ φ, ]() ã p Ä p φ() φ(a)e p ( a)ä + (k ) ã k+ ÄI L ä [ k+ φ, ](). D L [φ,, ]() + D L [φ,, ]() = φ () + (p ) p=3 + (k 3) which proves the case of k = for the approimation of φ. p=q+ é ã p p φ(a)e ( a) ã k+ I L [ k+ φ, ]() ã p Ä p φ() φ(a)e p ( a)ä ã k+ I L [ k+ φ, ]() + ã k+ ÄI L ä [ k+ φ, ](), Following a similar argument, one can easily prove the case of k = 3. With little modification, the proof holds for φ + as well...3 Boundary values reconstruction We have shown that the derivatives m φ(a) and m φ(b), m, are needed in (.6) when dealing with the non-periodic boundary conditions. For an outflow boundary condition in the sense of characteristics propagation, it is natural to employ high order etrapolation to obtain the boundary values. For an inflow boundary condition, it becomes more complicated. In [, 4], the ILW boundary treatment methodology was designed to reconstruct accurate solution values at the ghost points near the boundary for a class of fast sweeping finite difference WENO schemes that solves the static H-J equations. Over the past several years, ILW has eperienced systematic development in various settings, see [37, 38]. Here, we will apply this methodology to obtain the boundary values of solution derivatives at an inflow boundary. The main procedure of ILW is to convert the spatial derivatives into the time derivatives at an inflow boundary by repeatedly utilizing the underlying PDE. To illustrate the main idea, we consider the left boundary = a for eample:
13 If = a is an outflow boundary of the domain, where no physical boundary condition should be given, then we obtain all the derivatives m φ(a), m, through etrapolation with a suitable order of accuracy. If = a is an inflow boundary of the domain, with a Dirichlet boundary condition φ(a, t) = f(t). Then, the following ILW procedure is applied:. To obtain φ (a), we first note that φ t (a, t) = f (t). We evaluate (.) at = a and obtain f (t) + H(φ (a)) =. (.8) Then φ (a) can be solved from (.8). Note that there might be more than one root. In this case we should choose the root φ (a) so that H (φ (a)) >, (.9) which guarantees that the boundary = a is indeed an inflow boundary. Moreover, if the condition (.9) still cannot pin down a root, then we would choose the root which is closest to the value obtained through etrapolation.. To obtain φ (a), we first take the derivative with respect to t in the original H-J equation (.), yielding φ tt = H(φ ) t = H (φ )φ t = H (φ ) ( H(φ )) = ( H (φ ) ) φ. (.3) Note that at = a, φ tt = f (t) and φ (a) is available from the last step, we then solve this equation for φ (a); that is φ (a) = f (t) H (φ (a)). (.3) 3. To obtain φ (a), we further differentiate (.3) with respect to t and make use of the original H-J equation (.) again, we then have φ ttt = 3H (φ ) ( H (φ ) ) φ ( H (φ ) ) 3 φ. (.3) 3
14 At = a, note that φ ttt = f (t), then φ (a) is the only unknown quantity and can be easily solved φ (a) = f (t) + 3H (φ (a)) (H (φ (a))) φ (a) H (φ (a)) 3. (.33) 4. Repeating this procedure, we can compute any order derivative at = a as needed. If = a is an inflow boundary of the domain with a Neumann boundary condition φ (a, t) = g(t), the same ILW procedure can be applied, ecept that we start from solving the second derivative φ (a). Similarly, the ILW procedure for the right boundary = b can be developed. Note that in this case, the inflow condition (.9) should be changed into.3 Space Discretization H (φ (b)) <. (.34) In this subsection, we present the details about the spatial discretization of D L and D R. We denote L [v, ]( i ) as L i on each grid point i, where can be L and R. Note that I L i = e i I L i + J L i, i =,..., N, I L =, (.35a) I R i = e i+ I R i+ + J R i, i =,..., N, I R N =, (.35b) where, J L i i i+ = v(y)e (i y) dy, Ji R = v(y)e (y i) dy. (.36) i i Therefore, once J L i and J R i are computed for all i, we then can obtain I L i and I R i via the recursive relation (.35). In [3], we proposed a high order and robust methodology to compute J L i and J R i. In what follows we briefly review the underlying procedure for completeness of the paper. Note that, to approimate J L i with k th order accuracy, we may choose the interpolation stencil S(i) = { i r,..., i r+k }, which contains i and i. There is a unique polynomial p() of degree at most k that interpolates v() at the nodes in S(i). Then J L i is approimated by J L i i p(y)e (i y) dy. i (.37) 4
15 Note that, the integral on the right hand side can be evaluated eactly. Similarly, we can approimate J R i by J R i i+ p(y)e (y i) dy, (.38) i with polynomial p() interpolating v() on stencil S(i) = { i+r k,..., i+r } which includes i and i+. However, in some cases, the linear quadrature formula (.37) and (.38) with a fied stencil may generate the entropy-violating solutions, which will be demonstrated in Section 4. The main reason is that the viscosity solution of the H-J equation may involve discontinuous derivative. Hence, an approimation with a fied stencil such as (.37) and (.38) will develop spurious oscillations and eventually lead to failure of the scheme. We found that the WENO quadrature and the nonlinear filter, discussed in [3], are effective to control oscillations and capture the correct solution for solving nonlinear degenerate parabolic equations. In this work, we adapt this methodology to solve H-J equations as both types of equations may develop discontinuous derivatives for the solutions. The numerical evidence provided in Section 4 justifies the effectiveness of the WENO method and nonlinear filter in controlling oscillation as well as in capturing the viscosity solution to the H-J equation. In summary, when φ is periodic for eample, we will use the following modified formulation instead of (.) to approimate φ ± : φ,i = D L[φ, ]( i ) + φ +,i = D R[φ, ]( i ) p= p= σ p i,l Dp L [φ, ]( i), σ p i,r Dp R [φ, ]( i); (.39a) (.39b) while for non-periodic boundary conditions, (.5) are changed into φ,i = D L[φ,, ]( i ) + φ +,i = D R[φ,, ]( i ) p= p= σ p i,l D L[φ,p, ]( i ), σ p i,r D R[φ,, ]( i ). (.4a) (.4b) Note that the WENO quadrature is only applied for approimating the operators with p = in (.39), and cheaper high order linear formula are used for those with p >. The filter σ i,l and σ i,r are generated based on the smoothness indicators from the WENO methodology. 5
16 Figure.: The structure of the stencils in WENO integration. Here, the fifth order WENO-based quadrature for approimating J L i is provided as an eample. The corresponding stencil is presented in Figure.. We choose the big stencil as S(i) = { i 3,..., i+ } and the three small stencils as S r (i) = { i 3+r,..., i+r }, r =,,. For simplicity, we only provide the formulas for the case of a uniform mesh, i.e., i = for all i. Note that the WENO methodology is still applicable to the case of a nonuniform mesh, see [36].. On each small stencil S r (i), we obtain the approimation as i Ji,r L = e (i y) p r (y)d, i (.4) where p r () is the polynomial interpolating v() on nodes S r (i).. Similarly, on the big stencil S(i), we have J L i i = e (i y) p(y)d = i d r Ji,r, L (.4) r= with the linear weights d r satisfying r= d r =. 3. We change the linear weights d r into the following nonlinear weights ω r, which are defined as with ω r = ω r / ω s, r =,,, (.43) s= ω r = d r + τ ã 5. ɛ + β r Here, ɛ > is a small number to avoid the denominator becoming zero. We take ɛ = 6 in our numerical tests. The smoothness indicator β r is defined by β r = 3 l= i l 3 i i 6 Ç l å p r () d, (.44) l
17 which is used to measure the relative smoothness of the function v() in the stencil S r (i). In particular, we have the epressions as β = 3 ( v i 3 + 3v i 3v i + v i ) + 4 (v i 3 5v i + 7v i 3v i ), β = 3 ( v i + 3v i 3v i + v i+ ) + 4 (v i v i v i + v i+ ), β = 3 ( v i + 3v i 3v i+ + v i+ ) + 4 ( 3v i + 7v i 5v i+ + v i+ ). τ 5 is simply defined as the absolute difference between β and β τ 5 = β β. Moreover, we define a parameter ξ i as which will be use to construct the nonlinear filter. Here, ξ i = β min β ma, (.46) Ç å Ç å τ 5 τ 5 β ma = +, β min = +. ɛ + min(β, β ) ɛ + ma(β, β ) Note that we have employed the nonlinear weights based on the idea of the WENO-Z scheme proposed in [3], which enjoys less dissipation and higher resolution compared with the original WENO scheme. 4. Lastly, we define the approimation The filter σ i,l is defined as Ji L = ω r Ji,r. L (.48) r= σ i,l = min(ξ i, ξ i ). (.49) The process to obtain J R i and σ i,r is mirror symmetric to that of J L i and σ i,l with respect to point i. 3 Two-dimensional Case Consider the D H-J equation φ t + H(φ, φ y ) =, φ(, y, ) = φ (, y). (3.) 7
18 Let ( i, y j ) be the (i, j) node of a D orthogonal grid, which can be nonuniform. Denote by i = i i and y j = y j y j and denote φ i,j (t) as the solution φ( i, y j, t). Then, as with the D case, the semi-discrete scheme is d dt φ i,j + Ĥ(φ,i,j, φ+,i,j ; φ y,i,j, φ+ y,i,j ) =, (3.) where Ĥ(u, u; v, v) is a Lipschitz continuous monotone flu that is consistent with the Hamiltonian H. In this work, we employ the following local La-Friedrichs flu given as Ĥ(u, u + ; v, v + ) = H( u + u + where,, v + v + ) α (u, u + ) u+ u α (u, u + ) = α y (v, v + ) = ma H (u, v), u I(u,u + ),C v D ma H (u, v). A v B,v I(v,v + ) α y (v, v + ) v+ v, (3.3) Here, H i (u, v) denotes the partial derivative of H(u, v) with respect to the i th argument. [A, B] is the value range of u ± and [C, D] is the value range of v ±. In (3.), φ ±,i,j and φ± y,i,j are approimations to φ ( i, y j ) and φ y ( i, y j ), respectively, which can be computed directly by the proposed D formulation. In particular, when constructing φ ±, we fi y j and apply the D operators in -direction. For instance, if φ is periodic in -direction, then φ,i,j D p L [φ(, y j), ]( i ), p= φ +,i,j D p R [φ(, y j), ]( i ). p= Here, we choose = β/(α t) with α = ma A u B,C v D H (u, v). Similarly, to approimate φ ± y, we fi i and obtain φ y,i,j y D p L [φ( i, ), y ](y j ), p= φ + y,i,j y D p R [φ( i, ), y ](y j ), p= with α y = ma A u B,C v D H (u, v) and y = β/(α y t). Note that in the D case, we need to choose β ma as half of that for the D case to ensure the unconditional stability of the scheme. For a non-periodic boundary condition, the derivative values at an inflow boundary are still calculated through ILW. To illustrate the idea, we consider a rectangular domain [a, b ] [a y, b y ] for simplicity. Without loss of generality, assume the left boundary = a is an inflow Dirichlet boundary that satisfies φ(a, y, t) = f(y, t), y [a y, b y ] 8
19 and f(y, t) is a given function. For a fied time t, denote φ (k,k ), j = k +k k y k φ(a, y j, t), where (a, y j ) is a boundary grid point. For a k th order approimation, we need to compute the derivative values φ (k,), j for k k as in the D case. Note that φ t and φ y at boundary = a can be easily obtained, which are the functions of y. Hence, φ,, j can be computed through relation f t (y j ) + H(φ (,), j, f y(y j )) =. (3.4) In the case of more than one roots, we should use H (φ (,), j, f y(y j )), to single out the correct solution of φ (,), j that ensures the left boundary = a is an inflow boundary. Similar to the D case, if this condition still cannot pin down a root, then we would choose the root that is closest to the value obtained from etrapolation. If we take the derivative with respect to y for the original H-J equation (3.), then we obtain φ ty = H (φ, φ y )φ y H (φ, φ y )φ yy. (3.5) Notice that at the point (a, y j ), the equation becomes f ty (y j ) = H (φ (,), j, f y(y j ))φ (,), j H (φ (,), j, f y(y j ))f yy (y j ). (3.6) The only unknown quantity in (3.5) is φ (,), j, which can be easily solved. Furthermore, we take the derivative with respect to t for (3.) and get φ tt = (H ) φ + H H φ y + (H ) φ yy. (3.7) When evaluating equation (3.7) at point (a, y j ), the only unknown quantity is φ (,), j, which can be easily obtained by solving (3.5). Repeatedly applying this ILW procedure, we are able to obtain the derivative values φ (k,), j for k k. We now consider an inflow Neumann boundary condition imposed at boundary = a, or equivalently φ (a, y, t) = g(y, t) is specified, where g(y, t) is a given function. In this case, at a fied time t, we directly get φ (,), j = g(y j ). To obtain φ (,), j, we first differentiate the H-J equation with respect to, yielding φ t = H (φ, φ y )φ H (φ, φ y )φ y. 9
20 When evaluating at (a, y j ), the equation becomes f t (y j ) = H (φ (,), j, φ(,) j )φ (,), j H (φ (,), j, φ(,) j )g y(y j ), (3.8) where we have used the identity φ y (a, y) = g y (y). Note that, unlike the Dirichlet case, there are two unknowns in (3.8) including φ (,), j and φ (,), j. As suggested in [37], we can use the numerical differentiation to approimate the tangential derivatives at the boundary points, i.e., φ (,), j in this case. Once φ (,), j is obtained, we can solve for the only unknown quantity φ (,), j from (3.8). We then continue applying ILW to generate φ (k,), j for all k k and again the needed tangential derivatives are obtained by numerical differentiation. As in the D case, for an outflow boundary, we compute derivative values through etrapolation with suitable orders of accuracy. Lastly, we remark that, as a remarkable property, ILW can be conveniently applied to handle comple geometry. In fact, we only need to rotate the coordinate system locally at the underlying boundary point so that the new aes point to the normal and tangential directions of the boundary. By doing so, we are able to apply the above ILW procedure and generate the required derivative values at the inflow boundary. 4 Numerical results In this section, we present the results of the proposed schemes to demonstrate their efficiency and efficacy. For the D cases, we choose the time step as t = CFL α, (4.) where α is the maimum wave propagation speed. For the D cases, the time step is chosen as t = CFL ma(α /, α y / y), (4.) where α and α y are the maimum wave propagation speeds in and y directions, respectively. We use the fifth order linear or WENO quadrature formula to compute J L i and J R i. Unless additional declaration is made, linear quadrature formula will be adopted. Note that the lower order temporal error will always dominate the total numerical error in our eamples. Furthermore, ecept we make additional declaration, we only report the numerical results computed by the scheme with k = 3 for brevity and two CFLs including.5 and are used to demonstrate the performance.
21 Eample 4.. We first test the accuracy of the schemes for solving the linear problem φt + φ =, π π φ () = sin(). (4.3) with a Dirichlet boundary condition φ( π, t) = sin( π t). The eact solution is φ(, t) = sin( t). In Table 4., we test the problem on a uniform mesh and report the errors and the associated order of accuracy at time T =. Note that when k =, the second order accuracy is observed, which is unepected since the error estimate only implies first order accuracy. Such superconvergence is observed because when applied to the linear problem, the first order scheme with β = is equivalent to the second order Crank-Nicolson scheme in the MOL T framework []. For general nonlinear problems, the superconvergence cannot be observed anymore, demonstrating the sharpness of the error estimate. Table 4.: Eample 4.: L Errors and orders of accuracy for linear equation on the uniform meshes with a Dirichlet boundary condition. T =. CFL.5 N k =. β =. k =. β =. k = 3. β =.. error order error order error order 8.878E-3.3E-.989E E E E E E E E E E E E E E-6..37E E E E-.335E E E E E E E E E E E E E E E E E E-.69E E E E E E E E E E E E E E E E Eample 4.. Net, we solve the D Burgers equation: φt + (φ + ) =, φ () = cos(π) (4.4) with a -periodic boundary condition.
22 When T =.5/π, the solution is still smooth. The L errors and orders of accuracy are presented with uniform and nonuniform meshes in Table 4.3 and Table 4.4, respectively. The nonuniform meshes are obtained from uniform meshes with % random perturbation. The designed order of accuracy is observed for both types of meshes. Note that the schemes allow for large CFL numbers due to the notable unconditional stability property. When T = 3.5/π, the solution of the Burgers equation has developed discontinuous derivative. We benchmark our method by comparing the numerical results with the eact solution in Figure 4.. Here, a uniform mesh with N = 4 is used. It is observed that our schemes are able to approimate the non-smooth solution very well without producing noticeable oscillations. Large CFLs can be used to save computational cost as with the linear case. Furthermore, we note that, compared with low order schemes, the high order one is able to capture the solution structure more accurately. Table 4.3: Eample 4.: L Errors and orders of accuracy for D Burgers equation on the uniform meshes. T =.5/π. CFL.5 N k =. β =. k =. β =. k = 3. β =.. error order error order error order.593e-.757e- 6.35E E E E E E E E E E E E E E E E E- 5.89E E E E E E E E E E E E E E E E E E-.753E- 3.9E E E E E E E E E E E E E E E E Eample 4.3. Consider the D H-J equation with a non-conve Hamiltonian: φt cos(φ + ) =, φ () = cos(π) (4.5)
23 . CFL=.5 CFL= CFL=.4 Eact. CFL=.5 CFL= CFL=.4 Eact. CFL=.5 CFL= CFL=.4 Eact (a) k =. β =. Uniform mesh (b) k =. β =. Uniform mesh (c) k = 3. β =.. Uniform mesh.. CFL=.5 CFL= CFL=.4 Eact. CFL=.5 CFL= CFL=.4 Eact. CFL=.5 CFL= CFL=.4 Eact (d) k =. β =. Nonuniform mesh. (e) k =. β =. Nonuniform mesh. (f) k = 3. β =.. Nonuniform mesh. Figure 4.: Eample 4.: D Burgers equation. T = 3.5/π. N = 4 grid points. 3
24 Table 4.4: Eample 4.: L Errors and orders of accuracy for D Burgers equation on the nonuniform meshes. T =.5/π. CFL.5 N k =. β =. k =. β =. k = 3. β =.. error order error order error order.56e-.79e E E E E E E E E E E E E E E E E E- 6.3E- 4.75E E E E E E E E E E E E E E E E E-.5E- 3.79E E E E E E E E E E E E E E E E with a -periodic boundary condition. When T =.5/π, the solution is still smooth. We report L errors and orders of accuracy in Table 4.5 to numerically demonstrate that our schemes can achieve the designed order accuracy even if large CFLs are used. When T =.5/π, the solution has developed discontinuous derivative. In Figure 4.3, we compare the numerical results with the eact solutions on a uniform mesh with N = 4. Similar to the previous eample, our schemes are able to approimate the solution accurately without producing noticeable oscillations. Eample 4.4. In this eample, we solve the D Riemann problem with a non-conve Hamiltonian: φt + ( 4 φ ) ( φ 4 ) =,, φ (4.6) () =. Inflow Dirichlet boundary conditions are imposed with φ(, t) = and φ(, t) =. Numerical results computed by the proposed schemes with different orders are shown in Figure 4.4. In the simulation, we compute the solution up to T = and use a mesh with N = 8 grids. We remark that, to capture the correct viscosity solution, we need to incorporate the WENO quadrature as well as nonlinear filter, which is justified by the numerical evidence (Figure 4(d)). This eample and Eample 4.7 are the only eamples in which the nonlinear WENO quadrature and filter are required. 4
25 Table 4.5: Eample 4.3: Errors and orders of accuracy for D equation with a non-conve Hamiltonian. T =.5/π. CFL.5 N k =. β =. k =. β =. k = 3. β =.. error order error order error order 4.5E-3.66E-3 4.5E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E- 3.69E- 7.6E E E E E E E E E E E E E E E E In addition, we plot the numerical solutions with CFL=, to demonstrate the unconditional stability property of our schemes. We also would like to remark that one can save computational cost by choosing large CFLs, e.g. CFL=, but the approimation quality deteriorates due to the corresponding large temporal error. If small CFLs are used, e.g, CFL=, our numerical schemes are able to generate high quality numerical results, but the cost increases. The optimal choice of CFL is problem-dependent. Eample 4.5. Now, we consider the D Burgers equation φt + (φ + φ y + ) =,, y, φ (, y) = cos(π( + y)/), (4.7) with a 4-periodic boundary condition in each direction. When T =.5/π, the solution is smooth. We report the L error and orders of accuracy in Table 4.6 on uniform meshes. It is observed that our schemes can achieve the designed order accuracy. In Figure 4.5, we plot the numerical solutions at time T =.5/π, when discontinuous derivative has developed. Again, our schemes are able to capture the non-smooth structures very well and generate high quality numerical results. 5
26 .5 CFL=.5 CFL= CFL= Eact.5 CFL=.5 CFL= CFL= Eact.5 CFL=.5 CFL= CFL= Eact (a) k =. β =..5.5 (b) k =. β =..5.5 (c) k = 3. β =.. Figure 4.3: Eample 4.3: D HJ equation with a non-conve Hamiltonian. T =.5/π. N = 4 grid points CFL=.5 CFL= CFL= CFL= Ref.6.8 CFL=.5 CFL= CFL= CFL= Ref.5.5 (a) k =. β =..5.5 (b) k =. β = CFL=.5 CFL= CFL= CFL= Ref k=3,cfl=.5 Without filter Without WENO.5.5 (c) k = 3. β = (d) Comparison. Figure 4.4: Eample 4.4: D Riemann problem with a non-conve Hamiltonian. T =. N = 8 grid points. 6
27 Table 4.6: Eample 4.5: L errors and orders of accuracy for D Burgers equation. T =.5/π. k =. β =. k =. β =.5. k = 3. β =.6. CFL N N y error order error order error order.573e-.77e E E E E E E E E E E E E E E- 6.3E- 4.E E-..83E E E E E E E E E E E E-.96E E E E E E-.35.8E E E E E E E E y y (a) Surface. CFL=.5. (b) Surface. CFL=. Figure 4.5: Eample 4.5: D Burgers equation. T =.5/π. 4 4 grid points. k = 3 with β =.6. Eample 4.6. In this eample, we solve the following D HJ equation with a non-conve Hamiltonian φt cos(φ + φ y + ) =,, y, φ (, y) = cos(π( + y)/), (4.8) with a 4-periodic boundary condition in each direction. In table 4.7, we report the L error and the orders of accuracy to demonstrate that the proposed schemes can achieve the designed order accuracy if the solution is smooth. In Figure 4.6, the numerical solutions at T =.5/π are plotted. As with the former eample, the schemes are observed to be able to capture the non-smooth solution structures without producing noticeable oscillations. 7
28 Table 4.7: Eample 4.6: Errors and orders of accuracy for D equation with a non-conve Hamiltonian. T =.5/π. k =. β =. k =. β =.5. k = 3. β =.6. CFL N N y error order error order error order E-3.9E E E E E E E E E E E E-3 4.3E-3.84E E E E E E E E E E E E-3.978E E E E E E E E E E y y (a) Surface. CFL=.5. (b) Surface. CFL=. Figure 4.6: Eample 4.6: D equation with a non-conve Hamiltonian. T =.5/π. 4 4 grid points. k = 3 with β =.6. Eample 4.7. We consider the D Riemann problem with a non-conve Hamiltonian φt + sin(φ + φ y ) =,, y, φ(, y, ) = π( y ) (4.9) and an outflow boundary condition is imposed. We compute the solutions up to T =. A uniform mesh with N = N y = 8 is used. We plot the numerical solutions in Figure 4.7. It is observed that the proposed schemes are able to capture the viscosity solution accurately. In particular, the numerical result agrees well with those computed by other methods, see e.g., [, 4]. We remark that in this eample, we should incorporate the robust WENO quadrature as well as the nonlinear filter to ensure the convergence towards the viscosity solution. Otherwise, the rarefaction wave is missing, and the numerical solution ehibits spurious oscillations, see Figure 7(e)-7(f). 8
29 .5 y y.5.5 (a) Contour. CFL=.5. (b) Surface. CFL=.5..5 y y.5.5 (c) Contour. CFL=. (d) Surface. CFL=..5 3 y y.5.5 (e) Contour. Without WENO quadrature. (f) Surface. Without WENO quadrature. Figure 4.7: Eample 4.7: D Riemann problem with a non-conve Hamiltonian. T =. 8 8 grid points. k = 3 with β =.6. 9
30 .5.5 z y y (a) Solution. CFL=.5. (b) Optimal control sign(φ y). CFL= z y y (c) Solution. CFL=. (d) Optimal control sign(φ y). CFL=. Figure 4.8: Eample 4.4: D Control problem. T =. 6 6 grid points. k = 3 with β =.6. Eample 4.8. We consider solving the following problem from optimal control: φt + sin(y)φ + (sin() + sign(φ y ))φ y sin(y) ( cos()) =, π, y, π φ (, y) =, (4.) with periodic boundary conditions. We plot the numerical solutions at T = in Figure 4.8. A uniform mesh with N = N y = 6 is used. The optimal control sign(φ y ) which contains the most interesting information in the optimal control problem is shown as well. Again, our numerical schemes are able to resolve the non-smooth solution structures very nicely. Eample 4.9. In this eample, we simulate a propagating surface on a rectangular domain by solving the following equation φt ( ɛk)» φ + φ y + =, φ (, y) = 4 (cos(π) )(cos(πy) ), y (4.) 3
31 .8 t=.9.8 t=.9.6 t=.6.6 t= t=.3. t= t= t= y y.5 (a) Propagating surface. CFL=.5. (b) Propagating surface. CFL=. Figure 4.9: Eample 4.9: propagating surface. ɛ =. 6 6 grid points. k = 3 with β =.6. where K is the mean curvature defined by K = φ ( + φ y) φ y φ φ y + φ yy ( + φ ) (φ + φ y + ) 3/, (4.) and ɛ is a small constant. Periodic boundary conditions are imposed. Note that if ɛ, we need to discretize the second derivatives u, u y and u yy as well. The approach proposed in our previous work [3] is employed for the second derivative discretization. In the simulation, a uniform mesh with N = N y = 6 is used. In Figures , we plot the snapshots of numerical solutions for ɛ = (pure convection) at t =,.3,.6,.9 and for ɛ =. at t =,.,.3,.6, respectively. Eample 4.. As our last eample, we consider the same problem of a propagating surface (4.) with ɛ = but on the unit disk + y. The initial condition is ã π φ(, y, ) = sin ( + y ). A Dirichlet boundary condition φ(, y, t) = + t, for all + y = 3
32 .6.4 t= t=.6. t=.3. t= t=...6 t= t=.85 t= y y.5 (a) Propagating surface. CFL=.5. (b) Propagating surface. CFL=. Figure 4.: Eample 4.9: propagating surface. ɛ = grid points. k = 3 with β =.4. 3
33 .5 Y X Figure 4.: Eample 4.: The discretization of the domain which is embedded in a Cartesian grid with N = N y = 6. is imposed. We would like to use this eample to demonstrate the ability of the proposed methods to handle comple geometry and general boundary conditions. The domain is discretized by embedding the boundary in a regular Cartesian mesh. The mesh points include all the interior points as well as the intersections of all possible Cartesian grid lines = i and y = y j with the boundary curve, as illustrated in Figure 4.. In Figure 4., we plot several snapshots of the evolution of the surface. 5 Conclusion In this paper, we developed a class of high order non-oscillatory numerical schemes with unconditional stability for solving the nonlinear Hamilton-Jacobi (H-J) equations. The development of the schemes was based on our previous work [3], in which the spatial derivatives of a function were represented as a special kernel-based formulation. By incorporating a robust WENO methodology as well as a novel nonlinear filter, the proposed methods are able to suppress spurious oscillations and capture the correct viscosity solution. The high order eplicit strong-stability-preserving Runge-Kutta methods are employed for time integration, which do not require matri inversion or solving nonlinear algebraic equations. As the most notable advantage, the proposed schemes attain unconditional stability, high order accuracy (up to third order accuracy), and the essentially non-oscillatory property at the same time. In addition, the methods are able to handle complicated geometry and general boundary conditions if we further incorporate the inverse La-Wendroff 33
34 t=.8 t= t=. t=..5 t=.6.5 t=.6.5 t=..5 t=..5.5 Y X.5.5 Y X (a) Propagating surface. CFL=.5. (b) Propagating surface. CFL=. Figure 4.: Eample 4.: propagating surface on the unit disk. k = 3 with β =.6. boundary treatment technique. In the future, we plan to etend our schemes to solve other types of time-dependent problems, such as hyperbolic conservation laws. References [] R. Abgrall. Numerical discretization of the first-order Hamilton Jacobi equation on triangular meshes. Communications on pure and applied mathematics, 49(): , 996. [] J. Barnes and P. Hut. A hierarchical O(N log N) force-calculation algorithm. Nature, 34: , 986. [3] R. Borges, M. Carmona, B. Costa, and W. S. Don. An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws. Journal of Computational Physics, 7(6):39 3, 8. [4] M. Causley, H. Cho, and A. Christlieb. Method of lines transpose: Energy gradient flows using direct operator inversion for phase-field models. SIAM Journal on Scientific Computing, 39(5):B968 B99, 7. 34
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations 1
Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations Feng Zheng, Chi-Wang Shu 3 and Jianian Qiu 4 Abstract In this paper, a new type of finite difference Hermite weighted essentially
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationSemi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations
Semi-Discrete Central-Upwind Schemes with Reduced Dissipation for Hamilton-Jacobi Equations Steve Bryson, Aleander Kurganov, Doron Levy and Guergana Petrova Abstract We introduce a new family of Godunov-type
More informationReceived 6 August 2005; Accepted (in revised version) 22 September 2005
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol., No., pp. -34 Commun. Comput. Phys. February 6 A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a
More informationHigh order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes
High order finite difference Hermite WENO schemes for the Hamilton-Jacobi equations on unstructured meshes Feng Zheng, Chi-Wang Shu and Jianxian Qiu 3 Abstract In this paper, a new type of high order Hermite
More informationALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS
ALTERNATING EVOLUTION DISCONTINUOUS GALERKIN METHODS FOR HAMILTON-JACOBI EQUATIONS HAILIANG LIU AND MICHAEL POLLACK Abstract. In this work, we propose a high resolution Alternating Evolution Discontinuous
More informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationHybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation
Numerical Analysis and Scientific Computing Preprint Seria Hybrid semi-lagrangian finite element-finite difference methods for the Vlasov equation W. Guo J. Qiu Preprint #21 Department of Mathematics University
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationHIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS
June 6, 7 :7 WSPC/Lecture Notes Series: 9in x 6in chapter HIGH ORDER NUMERICAL METHODS FOR TIME DEPENDENT HAMILTON-JACOBI EQUATIONS Chi-Wang Shu Division of Applied Mathematics, Brown University Providence,
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationPositivity-preserving high order schemes for convection dominated equations
Positivity-preserving high order schemes for convection dominated equations Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Xiangxiong Zhang; Yinhua Xia; Yulong Xing; Cheng
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationA central discontinuous Galerkin method for Hamilton-Jacobi equations
A central discontinuous Galerkin method for Hamilton-Jacobi equations Fengyan Li and Sergey Yakovlev Abstract In this paper, a central discontinuous Galerkin method is proposed to solve for the viscosity
More informationMATH 220 solution to homework 1
MATH solution to homework Problem. Define z(s = u( + s, y + s, then z (s = u u ( + s, y + s + y ( + s, y + s = e y, z( y = u( y, = f( y, u(, y = z( = z( y + y If we prescribe the data u(, = f(, then z
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationMulti-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws
Multi-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws Bruno Costa a, Wai Sun Don b, a Departamento de Matemática Aplicada, IM-UFRJ, Caia Postal 6853, Rio de Janeiro, RJ, C.E.P. 945-97,
More informationHFVS: An Arbitrary High Order Flux Vector Splitting Method
HFVS: An Arbitrary High Order Flu Vector Splitting Method Yibing Chen, Song Jiang and Na Liu Institute of Applied Physics and Computational Mathematics, P.O. Bo 8009, Beijing 00088, P.R. China E-mail:
More informationA Third Order Fast Sweeping Method with Linear Computational Complexity for Eikonal Equations
DOI.7/s95-4-9856-7 A Third Order Fast Sweeping Method with Linear Computational Complexity for Eikonal Equations Liang Wu Yong-Tao Zhang Received: 7 September 23 / Revised: 23 February 24 / Accepted: 5
More informationA new type of Hermite WENO schemes for Hamilton-Jacobi equations
A new type of Hermite WENO schemes for Hamilton-Jacobi equations Jun Zhu and Jianxian Qiu Abstract In this paper, a new type of Hermite weighted essentially non-oscillatory (HWENO) schemes is designed
More informationWeighted ENO Schemes
Xiaolei Chen Advisor: Prof. Xiaolin Li Department of Applied Mathematics and Statistics Stony Brook University, The State University of New York February 7, 014 1 3 Mapped WENO-Z Scheme 1D Scalar Hyperbolic
More informationHigh order finite difference WENO schemes with positivity-preserving. limiter for correlated random walk with density-dependent turning rates
High order finite difference WENO schemes with positivity-preserving limiter for correlated random walk with density-dependent turning rates Yan Jiang, Chi-Wang Shu and Mengping Zhang 3 Abstract In this
More informationSECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS. Min Tang. (Communicated by the associate editor name)
Manuscript submitted to AIMS Journals Volume X, Number X, XX 2X Website: http://aimsciences.org pp. X XX SECOND ORDER ALL SPEED METHOD FOR THE ISENTROPIC EULER EQUATIONS Min Tang Department of Mathematics
More informationA parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows.
A parametrized maximum principle preserving flux limiter for finite difference RK-WENO schemes with applications in incompressible flows Tao Xiong Jing-ei Qiu Zhengfu Xu 3 Abstract In Xu [] a class of
More informationExplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code
Eplicit Jump Immersed Interface Method: Documentation for 2D Poisson Code V. Rutka A. Wiegmann November 25, 2005 Abstract The Eplicit Jump Immersed Interface method is a powerful tool to solve elliptic
More informationHigh Order Semi-Lagrangian WENO scheme for Vlasov Equations
High Order WENO scheme for Equations Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Andrew Christlieb Supported by AFOSR. Computational Mathematics Seminar, UC Boulder
More informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationA discontinuous Galerkin finite element method for directly solving the Hamilton Jacobi equations
Journal of Computational Physics 3 (7) 39 415 www.elsevier.com/locate/jcp A discontinuous Galerkin finite element method for directly solving the Hamilton Jacobi equations Yingda Cheng, Chi-Wang Shu *
More informationIn particular, if the initial condition is positive, then the entropy solution should satisfy the following positivity-preserving property
1 3 4 5 6 7 8 9 1 11 1 13 14 15 16 17 18 19 1 3 4 IMPLICIT POSITIVITY-PRESERVING HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR CONSERVATION LAWS TONG QIN AND CHI-WANG SHU Abstract. Positivity-preserving
More informationBound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu
Bound-preserving high order schemes in computational fluid dynamics Chi-Wang Shu Division of Applied Mathematics Brown University Outline Introduction Maximum-principle-preserving for scalar conservation
More information30 crete maximum principle, which all imply the bound-preserving property. But most
3 4 7 8 9 3 4 7 A HIGH ORDER ACCURATE BOUND-PRESERVING COMPACT FINITE DIFFERENCE SCHEME FOR SCALAR CONVECTION DIFFUSION EQUATIONS HAO LI, SHUSEN XIE, AND XIANGXIONG ZHANG Abstract We show that the classical
More informationA class of the fourth order finite volume Hermite weighted essentially non-oscillatory schemes
Science in China Series A: Mathematics Aug., 008, Vol. 51, No. 8, 1549 1560 www.scichina.com math.scichina.com www.springerlink.com A class of the fourth order finite volume Hermite weighted essentially
More informationFUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC PROBLEMS WITH DISCONTINUITIES
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK FUNDAMENTALS OF LAX-WENDROFF TYPE APPROACH TO HYPERBOLIC
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationFLUX DIFFERENCE SPLITTING AND THE BALANCING OF SOURCE TERMS AND FLUX GRADIENTS
FLUX DIFFERENCE SPLITTING AND THE BALANCING OF SOURCE TERMS AND FLUX GRADIENTS M.E.Hubbard The University of Cambridge, DAMTP, Silver Street, Cambridge, CB3 9EW, United Kingdom and P.Garcia-Navarro Area
More informationMulti-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws
Multi-Domain Hybrid Spectral-WENO Methods for Hyperbolic Conservation Laws Bruno Costa Wai Sun Don February 5, 6 Abstract In this article we introduce the multi-domain hybrid Spectral-WENO method aimed
More informationAn adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions
An adaptive mesh redistribution method for nonlinear Hamilton-Jacobi equations in two- and three-dimensions H.-Z. Tang, Tao Tang and Pingwen Zhang School of Mathematical Sciences, Peking University, Beijing
More informationAn Improved Non-linear Weights for Seventh-Order WENO Scheme
An Improved Non-linear Weights for Seventh-Order WENO Scheme arxiv:6.06755v [math.na] Nov 06 Samala Rathan, G Naga Raju Department of Mathematics, Visvesvaraya National Institute of Technology, Nagpur,
More informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
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 informationFunctions with orthogonal Hessian
Functions with orthogonal Hessian B. Dacorogna P. Marcellini E. Paolini Abstract A Dirichlet problem for orthogonal Hessians in two dimensions is eplicitly solved, by characterizing all piecewise C 2 functions
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More informationAN EXACT RESCALING VELOCITY METHOD FOR SOME KINETIC FLOCKING MODELS
AN EXACT RESCALING VELOCITY METHOD FOR SOME KINETIC FLOCKING MODELS THOMAS REY AND CHANGHUI TAN Abstract. In this work, we discuss kinetic descriptions of flocking models, of the so-called Cucker- Smale
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationImproved Seventh-Order WENO Scheme
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Eposition 4-7 January 2, Orlando, Florida AIAA 2-45 Improved Seventh-Order WENO Scheme Yiqing Shen Gecheng Zha Dept.
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationNUMERICAL CONSERVATION PROPERTIES OF H(DIV )-CONFORMING LEAST-SQUARES FINITE ELEMENT METHODS FOR THE BURGERS EQUATION
SIAM JOURNAL ON SCIENTIFIC COMPUTING, SUBMITTED JUNE 23, ACCEPTED JUNE 24 NUMERICAL CONSERVATION PROPERTIES OF H(DIV )-CONFORMING LEAST-SQUARES FINITE ELEMENT METHODS FOR THE BURGERS EQUATION H. DE STERCK,
More informationA weighted essentially non-oscillatory numerical scheme for a multi-class traffic flow model on an inhomogeneous highway
Manuscript A weighted essentially non-oscillatory numerical scheme for a multi-class traffic flow model on an inhomogeneous highway Peng Zhang, S.C. Wong and Chi-Wang Shu. Department of Civil Engineering,
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationHierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws
Hierarchical Reconstruction with up to Second Degree Remainder for Solving Nonlinear Conservation Laws Dedicated to Todd F. Dupont on the occasion of his 65th birthday Yingjie Liu, Chi-Wang Shu and Zhiliang
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationFifth order multi-moment WENO schemes for hyperbolic conservation laws
J. Sci. Comput. manuscript No. (will be inserted by the editor) Fifth order multi-moment WENO schemes for hyperbolic conservation laws Chieh-Sen Huang Feng Xiao Todd Arbogast Received: May 24 / Accepted:
More informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationDevelopment and stability analysis of the inverse Lax-Wendroff boundary. treatment for central compact schemes 1
Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes François Vilar 2 and Chi-Wang Shu 3 Division of Applied Mathematics, Brown University, Providence,
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationYINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG
CENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION YINGJIE LIU, CHI-WANG SHU, EITAN TADMOR, AND MENGPING ZHANG Abstract. The central scheme of
More informationConvex ENO Schemes for Hamilton-Jacobi Equations
Convex ENO Schemes for Hamilton-Jacobi Equations Chi-Tien Lin Dedicated to our friend, Xu-Dong Liu, notre Xu-Dong. Abstract. In one dimension, viscosit solutions of Hamilton-Jacobi (HJ equations can be
More informationSMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES
Proceedings of ALGORITMY 2016 pp. 155 164 SMOOTHNESS INDICATORS FOR WENO SCHEME USING UNDIVIDED DIFFERENCES TAMER H. M. A. KASEM AND FRANÇOIS G. SCHMITT Abstract. The weighted essentially non-oscillatory
More informationChapter 3 Burgers Equation
Chapter 3 Burgers Equation One of the major challenges in the field of comple systems is a thorough understanding of the phenomenon of turbulence. Direct numerical simulations (DNS) have substantially
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 informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Mathematics, Vol.4, No.3, 6, 39 5. ANTI-DIFFUSIVE FINITE DIFFERENCE WENO METHODS FOR SHALLOW WATER WITH TRANSPORT OF POLLUTANT ) Zhengfu Xu (Department of Mathematics, Pennsylvania
More informationICES REPORT A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws
ICES REPORT 7- August 7 A Multilevel-WENO Technique for Solving Nonlinear Conservation Laws by Todd Arbogast, Chieh-Sen Huang, and Xikai Zhao The Institute for Computational Engineering and Sciences The
More informationLecture 16. Wednesday, May 25, φ t + V φ = 0 (1) In upwinding, we choose the discretization of the spatial derivative based on the sign of u.
Lecture 16 Wednesday, May 25, 2005 Consider the linear advection equation φ t + V φ = 0 (1) or in 1D φ t + uφ x = 0. In upwinding, we choose the discretization of the spatial derivative based on the sign
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationHigh-order finite difference methods with subcell resolution for 2D detonation waves
Center for Turbulence Research Annual Research Briefs 7 High-order finite difference methods with subcell resolution for D detonation waves B W. Wang, C.-W. Shu, H. C. Yee AND B. Sjögreen. Motivation and
More information2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13
2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13 REVIEW Lecture 12: Classification of Partial Differential Equations (PDEs) and eamples with finite difference discretizations Parabolic PDEs Elliptic
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationImproving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on
Improving Dynamical Core Scalability, Accuracy, and Limi:ng Flexibility with the ADER- DT Time Discre:za:on Matthew R. Norman Scientific Computing Group National Center for Computational Sciences Oak Ridge
More informationAnalysis and Computation for a Fluid Mixture Model of Tissue Deformations. by Qunlei Jiang
ABSTRACT JIANG, QUNLEI. Analysis and Computation for a Fluid Miture Model of Tissue Deformations. (Under the direction of Professors Zhilin Li and Sharon R. Lubkin). A fluid miture model of tissue deformations
More informationarxiv: v2 [math.na] 4 Mar 2014
Finite Difference Weighted Essentially Non-Oscillatory Schemes with Constrained Transport for Ideal Magnetohydrodynamics arxiv:1309.3344v [math.na] 4 Mar 014 Andrew J. Christlieb a, James A. Rossmanith
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
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 informationA New Fifth Order Finite Difference WENO Scheme for Hamilton-Jacobi Equations
A New Fift Order Finite Difference WENO Sceme for Hamilton-Jacobi Equations Jun Zu, 1 Jianxian Qiu 1 College of Science, Nanjing University of Aeronautics and Astronautics, Nanjing, Jiangsu 10016, People
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 informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationAdvection algorithms II. Flux conservation, subgrid models and flux limiters
Chapter 4 Advection algorithms II. Flu conservation, subgrid models and flu limiters In this chapter we will focus on the flu-conserving formalism of advection algorithms, and we shall discuss techniues
More informationA Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem
A Fourth Order Accurate Discretization for the Laplace and Heat Equations on Arbitrary Domains, with Applications to the Stefan Problem Frédéric Gibou Ronald Fedkiw April 27, 2004 Abstract In this paper,
More informationNumerical Smoothness and PDE Error Analysis
Numerical Smoothness and PDE Error Analysis GLSIAM 2013 Central Michigan University Tong Sun Bowling Green State University April 20, 2013 Table of contents Why consider Numerical Smoothness? What is Numerical
More informationA Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations
J Sci Comput (00) 45: 404 48 DOI 0.007/s095-009-9340-y A Central Discontinuous Galerkin Method for Hamilton-Jacobi Equations Fengyan Li Sergey Yakovlev Received: 6 August 009 / Revised: November 009 /
More informationA recovery-assisted DG code for the compressible Navier-Stokes equations
A recovery-assisted DG code for the compressible Navier-Stokes equations January 6 th, 217 5 th International Workshop on High-Order CFD Methods Kissimmee, Florida Philip E. Johnson & Eric Johnsen Scientific
More informationA STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS *
Journal of Computational Mathematics Vol.8, No.4, 010, 55 568. http://www.global-sci.org/jcm doi:10.408/jcm.1003-m0016 A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationCharbel Farhat, Arthur Rallu, Sriram Shankaran
A Higher-Order Generalized Ghost Fluid Method for the Poor for the Three-Dimensional Two-Phase Flow Computation of Underwater Eplosions and Implosions Charbel Farhat, Arthur Rallu, Sriram Shankaran Department
More informationRunge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter
Runge-Kutta discontinuous Galerkin method with a simple and compact Hermite WENO limiter Jun Zhu, inghui Zhong, Chi-Wang Shu 3 and Jianxian Qiu 4 Abstract In this paper, we propose a new type of weighted
More informationc 2012 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 5, No., pp. 544 573 c Society for Industrial and Applied Mathematics ARBITRARILY HIGH-ORDER ACCURATE ENTROPY STABLE ESSENTIALLY NONOSCILLATORY SCHEMES FOR SYSTEMS OF CONSERVATION
More informationNUMERICAL STUDIES ON HAMILTON-JACOBI EQUATION IN NON-ORTHOGONAL COORDINATES
J. Appl. Math. & Computing Vol. 21(2006), No. 1-2, pp. 603-622 NUMERICAL STUDIES ON HAMILTON-JACOBI EQUATION IN NON-ORTHOGONAL COORDINATES MINYOUNG YUN Abstract. The level-set method is quite efficient
More informationJournal of Computational Physics
Journal of Computational Physics 3 () 3 44 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp A local discontinuous Galerkin method
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationA Strong Stability Preserving Analysis for Explicit Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions.
A Strong Stability Preserving Analysis for Explicit Multistage Two-Derivative Time-Stepping Schemes Based on Taylor Series Conditions. Zachary Grant 1, Sigal Gottlieb 2, David C. Seal 3 1 Department of
More informationLevel Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows
Level Set and Phase Field Methods: Application to Moving Interfaces and Two-Phase Fluid Flows Abstract Maged Ismail Claremont Graduate University Level Set and Phase Field methods are well-known interface-capturing
More informationFAST HUYGENS SWEEPING METHODS FOR SCHRÖDINGER EQUATIONS IN THE SEMI-CLASSICAL REGIME
FAST HUYGENS SWEEPING METHODS FOR SCHRÖDINGER EQUATIONS IN THE SEMI-CLASSICAL REGIME SHINGYU LEUNG, JIANLIANG QIAN, AND SUSANA SERNA Abstract. We propose fast Huygens sweeping methods for Schrödinger equations
More informationarxiv: v1 [math.na] 15 Apr 2014
High order maimum principle preserving discontinuous Galerkin method for convection-diffusion equations Tao Xiong Jing-Mei Qiu Zhengfu Xu 3 arxiv:44.46v [math.na] 5 Apr 4 Abstract In this paper, we propose
More informationGeometric Modeling Summer Semester 2010 Mathematical Tools (1)
Geometric Modeling Summer Semester 2010 Mathematical Tools (1) Recap: Linear Algebra Today... Topics: Mathematical Background Linear algebra Analysis & differential geometry Numerical techniques Geometric
More information