Finite difference Hermite WENO schemes for the Hamilton-Jacobi equations 1
|
|
- Lily Holmes
- 5 years ago
- Views:
Transcription
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 non-oscillatory (HWENO) schemes are constructed for solving Hamilton-Jacobi (HJ) equations. Point values of both the solution and its first derivatives are used in the HWENO reconstruction and evolved via time advancing. While the evolution of the solution is still through the classical numerical flues to ensure convergence to weak solutions, the evolution of the first derivatives of the solution is through a simple dimension-by-dimension non-conservative procedure to gain efficiency. The main advantages of this new scheme include its compactness in the spatial field and its simplicity in the reconstructions. Etensive numerical eperiments in one and two dimensional cases are performed to verify the accuracy, high resolution and efficiency of this new scheme. Keywords: HWENO method; Hamilton-Jacobi equation; finite difference method. AMS(MOS) subject classification: 65M6, 35L65 The research of F. Zheng and J. Qiu is partially supported by NSFC grants 579 and NSAF grant U6347. The research of C.-W. Shu is partially supported by ARO grant W9NF-5--6 and NSF grant DMS School of Mathematical Sciences and Fujian Provincial Key Laboratory of Mathematical Modeling and High- Performance Scientific Computing, Xiamen University, Xiamen, Fujian 365, P. R. China. fzbz88-3@63.com. 3 Division of Applied Mathematics, Brown University, Providence, RI 9, USA. shu@dam.brown.edu 4 School of Mathematical Sciences and Fujian Provincial Key Laboratory of Mathematical Modeling & High-Performance Scientific Computing, Xiamen University, Xiamen, Fujian 365, P.R. China. jqiu@mu.edu.cn.
2 Introduction In this paper, we design a new finite difference Hermite weighted essentially non-oscillatory (HWENO) method to solve the following Hamilton Jacobi (HJ) equation t +H( ) =, (,) = () Ω R d where = (,,..., n ) T. We consider only up to two dimensions for simplicity and hence use and y instead of and in the sequel. The HJ equations come from many applications, from control theory and geometric optics, to image processingand level set method andso on. It is well known that solutions of the nonlinear HJ equation are typically continuous. However, the derivatives of the solutions could be discontinuous even though the initial condition () is smooth enough. The (weak) solution may not be unique unless suitable assumptions (viscosity solutions) are made, see, e.g. [4]. It is well known that HJ equations are closely related to conservation laws, so many successful numerical methods for solving conservation laws can be easily adapted to solve HJ equations. Hu and Shu [6] proposed a discontinuous Galerkin (DG) method to solve the HJ equation. Li and Shu [] reinterpreted and simplified the two dimensional method of Hu and Shu. Zhu and Qiu [, 3] used HWENO method for the HJ equations on both structured and unstructured meshes. Tao and Qiu [7] made use of the central HWENO method to solve the HJ equation successfully. As most of the finite volume methods resort to solving the conservation laws for the derivatives of the solution that seems to be not direct, Cheng and Shu [] designed DG methods to directly solve the HJ equation. Later, Yan and Osher [8] proposed a local DG (LDG) method to solve the HJ equation directly. In [3] Cheng and Wang improved the work in [] by utilizing the Roe speed and entropy fi at the cell interface, and based on [3], Zheng and Qiu [] developed HWENO schemes to directly solve the HJ equations. Comparing to the finite volume or DG methods, which evolve cell averages or a complete polynomial, the finite difference methods, which evolve only point values, may be easier to implement
3 and more efficient. Now, let us briefly review the early works. Crandall and Lion [5] introduced first order monotone schemes which can converge to the viscosity solution. Osher and Shu [] presented higher order finite difference ENO schemes, and Jiang and Peng [7] developed higher order finite difference weighted ENO (WENO) schemes. Qiu et al. [4, ] put forward Hermite WENO (HWENO) schemes for HJ equations with Runge-Kutta or La-Wendroff time discretization. Finite difference schemes on unstructured meshes have also been designed. Abgrall [] proposed first order monotone schemes on triangular meshes. Lafon and Osher [9] designed second order schemes. Higher order WENO schemes were developed by Zhang and Shu [9], and higher order HWENO schemes were proposed by Zhu and Qiu []. For a detailed review of high order HJ equations about both the finite difference and finite volume method, we recommend [5]. In this paper, following the methods developed in [3, 4, ], both and its first derivative (e.g. in one dimension) at the grid points are used in the HWENO method to reconstruct point values of the derivatives, and both of them are evolved by time marching. Comparing with the classical WENO method developed by Jiang and Shu in [8], the HWENO method requires etra work and storage but it is much more compact with the same order of accuracy. Comparing with DG methods, the HWENO method could achieve high order accuracy by adaptive stencils and hence could maintain the essentially non-oscillatory property. Notice that the HJ equations cannot be written in a conservation form, hence it would seem more natural to use finite difference methods based on point values instead of finite volume methods based on cell averages. In this way, we could avoid the costly multi-dimensional reconstructions and could use dimension-by-dimension interpolations, thereby reducing the computational cost and improving efficiency. It is important to use numerical flues(monotone Hamiltonians) to evolve the solution itself for HJ equations to ensure convergence to weak solutions. Since the derivatives of the HJ solution satisfy conservation laws, previous HWENO schemes tend to use conservative approimations to evolve them as well, hence causing complications and etra computational cost in multi-dimensional reconstruction. In this paper, we design a new HWENO scheme which evolves the solution itself through classical numerical 3
4 flues to ensure convergence to weak solutions, and evolves the derivatives of the solution in a nonconservative fashion, thereby allowing efficient dimension-by-dimension interpolation resulting in a gain of simplicity and efficiency. Etensive numerical eperiments are performed to illustrate the good performance of our schemes. This paper is organized as follows. In Section, we describe the detailed steps about construction and implementation of the finite difference HWENO schemes in both one and two dimensions for the HJ equations. In Section 3, we present etensive numerical results to verify the accuracy, stability and resolution of our method. Finally, a conclusion is given in Section 4. The numerical method for the Hamilton-Jacobi equations In this section, we will give the framework of the schemes first and then the detailed steps of the HWENO reconstructions for both one and two dimensional Hamilton-Jacobi equations.. One dimensional Hamilton-Jacobi equations The Hamilton-Jacobi equations in the one dimensional case can be written as: t +H( ) = (,) = () [a,b]. (.) For simplicity, we consider a uniform mesh that is defined as a = < < < N < N = b. However, thisassumptionisnotneededforourschemes. Wedenote j = ( j,t)asthenumerical approimation to the viscosity solution and u j = ( j,t) as the numerical approimation to its first derivative. Then, by taking the spatial derivative on both sides of (.), we obtain the following system of equations: d j dt = H( ) =j du j dt = H (u j )u =j where the H (u) = H u. We replace H( ) =j with a suitable monotone numerical flu, denoted by Ĥ(u j,u+ j ), and we will use the simple local La-Friedrichs flu as an eample in the sequel. For stability, we also split H (u j ) into a positive part and a negative part locally, denoted by H + j and (.) 4
5 H j respectively, and approimate u by upwind methods. Then, we have the following scheme: where d j dt = Ĥ(u j,u+ j ) du j dt = H+ j u j H j u+ j (.3) Ĥ(u j,u+ j ) = H(u j +u+ j ) α(u+ j u j ) (.4) with α = ma H (u), D = [min(u u D j,u+ j ),ma(u j,u+ j )], and H j = (H ( u j +u+ j ) H ( u j +u+ j ) ), H + j = (H ( u j +u+ j )+ H ( u j +u+ j ) ). (.5) u ± j andu± j,theleftandrightlimitsofthepointvaluesofu( j,t)andu ( j,t), willbeapproimated by the HWENO method which we will describe in detail in the net subsection. Notice that, even though we already have u j as part of our numerical solution, we would still want to obtain upwindbiased approimations u ± j through the point values of and u nearby (of course without using u j itself) in order to use the monotone flu to ensure stability and convergence to weak solutions. We observe that either H j or H+ j is equal to zero. In practical, we need to compute u j only when H + j is not equal to zero, and compute u+ j only when H j is not equal to zero. This will reduce the computational cost because only one of the u j and u+ j needs to be computed. After the spatial discretization, we can rewrite the scheme as U t = L(U), where L denotes the operator of the spatial discretization, and then use the third-order total variation diminishing (TVD) Runge-Kutta time discretization [6] to solve the semi-discrete form (.3): U () = U n + tl(u n ) U () = 3 4 Un + 4 (U() + tl(u () )) (.6) U n+ = 3 Un + 3 (U() + tl(u () )). HWENO reconstruction in one dimension u ± j. In this subsection, we will describe the HWENO reconstruction procedure for u ± j = ± j and 5
6 Step. Reconstruction of j by the HWENO method from the point values { i,u i }.. GiventhesmallstencilsS = { j, j, j },S = { j, j, j+ }, S = { j, j, j, j+ }, and the big stencil T = {S,S,S }, we construct Hermite cubic polynomials p (), p (), p (), and a fifth-degree polynomial q() such that: p ( j+i ) = j+i, i =,,, p ( j ) = u j p ( j+i ) = j+i, i =,,, p ( j+ ) = u j+ p ( j+i ) = j+i, i =,,,, q( j+i ) = j+i, i =,,,, q ( j± ) = u j± In fact, we only need the derivative values of these polynomials at the point = j, which have the following epressions: p ( j ) = j +4 j 5 i +4u j p ( j ) = j +4 j 5 j+ +u j+ 4 p ( j ) = j 6 j +3 j + j+ 6 q ( j ) = j +8 j 9 j j+ +9u j +3u j+ 8.. For each small stencil S m,m =,,, we compute the smoothness indicator, which measures the smoothness of the polynomials in each stencil: the smaller the indicator is, the smoother the polynomial is in the stencil. We use the formula similar to [8] to figure out the indicators, denoted by β m,m =,,: β m = 3 l= l ( l lp m( j )) m =,, (.7).3. We compute the linear weights, denoted by γ m,m =,,, satisfying: q ( j ) = γ m p m ( j ) m= for all the point values of { i } and {u i } in the big stencil T, which leads to γ = 4, γ = 3, γ = 5 6
7 .4. We compute the nonlinear weights based on the linear weights and the smoothness indicators by: ω m = ω m k= ω, ω m = k γ m (β m +ε), m =,, (.8) where ε is a small positive number to avoid the denominator becoming zero. In our numerical tests, we use ε = 6. The final HWENO approimation epression is: j = ω m p m ( j ), The reconstruction of u + j is mirror symmetric with respect to j of the above procedure. m= Step. Reconstruction of u j by the HWENO method from the point values { i,u i }.. GiventhesmallstencilsS = { j, j, j },S = { j, j, j+ }, S = { j, j, j, j+ }, and the big stencil T = {S,S,S }, we construct Hermite quartic polynomials p (), p (), p (), and a sith-degree polynomial q() such that: p ( j+i ) = j+i, i =,,, p ( j+i) = u j+i, i =, p ( j+i ) = j+i, i =,,, p ( j+i ) = u j+i, i =, p ( j+i ) = j+i, i =,,,, q( j+i ) = j+i, i =,,,, p ( j ) = u j q ( j+i ) = u j+, i =,, Again, we only need the second order derivative values of these polynomials at the point j, which have the following epressions: p ( j ) = u j +8u j + j +6 j 7 j ( ) p ( j ) = 4u j u j+ + j 8 j +7 j+ ( ) p ( j ) = 6u j j + j 5 j +4 j+ 6( ) q ( j ) = 8u j +8u j 6u j+ + j +54 j 8 j +6 j+ 8( ).. For each small stencil S m,m =,,, we compute the smoothness indicators respectively: β m = 4 l=3 l ( l lp m( j )), m =,, (.9) 7
8 .3. We compute the linear weights, denoted by γ m,m =,,, satisfying: q = γ m p m ( j ) m= for all the point values of { i } and {u i } in the big stencil T, which leads to γ = 4, γ = 3, γ = 5.4. We compute the nonlinear weights ω m as in (.8), the final HWENO approimation epression is: u j = ω m p m ( j ), m= The reconstruction of u + j is mirror symmetric with respect to j of the above procedure..3 Two dimensional Hamilton-Jacobi equations The Hamilton-Jacobi equation in the two dimensional case is written as t +H(, y ) = (,y,) = (,y) (,y) [a,b] [c,d] (.) For simplicity, we also assume the computational domain has been uniformly meshed as a = < < < N < N = b and c = y < y < < y N < y N = d. Furthermore, we define ij = ( i,y j,t) as the numerical approimation to the viscosity solution, and u ij = ( i,y j,t) and v ij = y ( i,y j,t) as the numerical approimations to its first order partial derivatives with respect to the variables and y, respectively. Then, by taking spatial derivatives on both sides of (.), we have the following system of equations: d ij = H(, y ) dt du ij = H (u,v)u H (u,v)v dt dv ij = H (u,v)u y H (u,v)v y dt (.) 8
9 where the H (u) = H u and H (u) = H v. As u y is equal to v in the smooth case, we can also rewrite the system of equations as the following: d ij = H(, y ) dt du ij = H (u,v)u H (u,v)u y dt dv ij = H (u,v)v H (u,v)v y dt (.) In the same way as before, we replace H(, y ) with a monotone numerical flu, and split H (u,v) and H (u,v) into a positive part and a negative part respectively. Then we discretize (.) into the following d ij dt du ij dt dv ij dt = Ĥ(u ij,u+ ij,v ij,v+ ij ) = H + ij u ij H ij u+ ij H iju y ij = H ij v ij H + ij v y ij H ij v+ y ij (.3) where Ĥ(u ij,u+ ij,v ij,v+ ij ) refers to a two-dimensional monotone flu, such as the simple local La- Friedrichs flu defined as Ĥ(u ij,u+ ij,v ij,v+ ij ) = H(u ij +u+ ij, v ij +v+ ij ) α(u+ ij u ij ) β(v+ ij v ij ), where α = ma u D,v E H (u,v) and β = ma u D,v E H (u,v). With regard to α, we take D as a local region and E as a global region, namely D = [min(u ij,u+ ij ),ma(u ij,u+ ij )] and E = [min(v,v + ),ma(v,v + )] [c,d]. We compute the coefficient β similarly, ecept that we then take D globally and E locally. This way of computing the local La-Friedrichs flu is needed to ensure monotonicity, see []. H ij, H ij are the average values of H and H at ( i,y j ) defined as the following: H ij = H ( u ij +u+ ij, v ij +v+ ij ), H ij = H ( u ij +u+ ij, v ij +v+ ij ), and H, H+, H, H+ as the negative and positive parts of H and H respectively, defined as follows: H ij = (H ij H ij ), H + ij = (H ij + H ij ), 9
10 H ij = (H ij H ij ), H + ij = (H ij + H ij ). ± ij, u± ij and v± ij are the left and right limits of the point values u( i,y j,t), u ( i,y j,t) and v ( i,y j,t) with respect to the variable, and ± y ij, u ± y ij and v y ± ij are the left and right limits of the point values v( i,y j,t), u y ( i,y j,t) and v y ( i,y j,t) with respect to the variable y. The values of ± ij, ± y ij, u ± ij and v± y ij are reconstructed with the one dimensional method in each direction with the other direction fied. As to the mied derivatives v,u y, we simply use the following fourth order central approimations in the and y directions, respectively: v ij = v i+,j +8v i+,j 8v i,j +v i,j, u yij = u i,j+ +8u i,j+ 8u i,j +u i,j. y In practice, we only need to compute either u ij or u+ ij but not both. v y ij and v + y ij can be computed similarly. Then, we can rewrite the scheme (.3) as U t = L(U), where L denotes the operator of spatial discretization, and use the third-order total variation diminishing (TVD) Runge-Kutta time discretization (.6) to solve the semi-discrete form (.3). 3 Numerical results In this section, we provide numerical eperiments for the fifth order HWENO method in one and two dimensional cases. In all the accuracy tests, we set t =.6 5 3/α in the one dimensional case and t = α β.6 y 5 3 in the two dimensional case, in order to guarantee that spatial numerical errors dominate. In other tests, we simply take t =.6 /α and t = α.6 + β.6 y and two dimensional cases respectively, unless otherwise indicated. 3. One dimensional case: accuracy test in the one Eample 3.: We solve the following linear scalar equation: t + =,
11 with the initial condition (, ) = sin(π), and periodic boundary condition. We compute the solution upto t =, thenumerical errorsandnumerical ordersof accuracy for thehweno method are shown in Table. We can see that the scheme achieves or eceeds fifth order accuracy. Table : t + =,(,) = sin(π). Periodic boundary conditions. t = N L error L order L error L order L error L order.58e-.44e-.6e- 9.4E E E E E E E E E E E E E E E Eample 3.: We solve the Burgers equation: t + ( +) =, with the initial condition (, ) = cos(π), and periodic boundary condition. We compute the solution up to t =.5 π. At this time, the solution is still smooth. The numerical results are shown in Table. Again, we can see that the scheme can reach its designed order of accuracy. Table : t + ( +) =,(,) = cos(π). Periodic boundary conditions. t =.5/π N L error L order L error L order L error L order.e E-4 8.4E-4.64E E E E E E E E E E E E E E E Eample 3.3: We solve the nonlinear scalar equation: t cos( +) =, < < with the initial condition (,) = cos(π), and periodic boundary condition. When t =.5 π, the solution is still smooth. The numerical results are shown in Table 3. Again, we observe that the scheme can achieve its designed accuracy.
12 Table 3: t cos( +) =,(,) = cos(π). Periodic boundary conditions. t =.5/π N L error L order L error L order L error L order.84e E-4 7.7E E-4.5.E E E E E E E E E E E E E E One dimensional case: discontinuous derivatives Eample 3.4: We solve the linear equation: t + = with the initial condition (,) = (.5) together with the periodic boundary condition, where () = ( π 3 )(+)+ cos( 3π ) 3 < 3, 3 +3cos(π) <, 3 5 3cos(π) < 3, 8+4π +cos(3π) 3 +6π( ) <. 3 We plot the result at t =. and t = 8. in Figure (a) and Figure (b) respectively, and observe that the HWENO schemes have good resolution for the corner singularity. Eample 3.5: We solve one dimensional Burgers equation: t + ( +) = with (,) = cos(π) and periodic boundary conditions. We compute the solution up to t = 3.5/π. Atthistime, thediscontinuousderivativehasalreadyappearedinthesolution. Weshowthe numerical results with the meshes N = 4 and N = 8 in Figure (a) and Figure (b) respectively. Again, we can see that the scheme gives high resolution for this problem.
13 - - T= T=8 - HWENO Eact - HWENO Eact (a) (b) Figure : One dimensional linear equation. N = cells. (a) t =., (b) t = 8.. Solid line: the eact solution; Square symbol: the HWENO scheme. N=4 N=8 -. HWENO Eact -. HWENO Eact (a) (b) Figure : One dimensional Burgers equation. t = 3.5/π. (a) N = 4 and (b) N = 8. Solid line: the eact solution; Square symbol: the HWENO scheme. 3
14 Eample 3.6: We solve the nonlinear equation with a non-conve / non-concave flu t cos( +) = with the initial data (, ) = cos(π) and periodic boundary conditions. This time, we compute the solution up to t =.5/π. We observe the results with N = 4 and N = 8, which are shown in Figure 3, and find that the method gives high resolution in this case. N=4 N=8.5 HWENO Eact.5 HWENO Eact (a) (b) Figure 3: H(u) = cos(u + ). t =.5/π. (a) N = 4 and (b) N = 8. Solid line: the eact solution; Square symbol: the HWENO scheme. Eample 3.7: We solve the problem t + 4 ( )( 4) = < < with the initial data (,) =. As the derivative of (,) is discontinuous, the initial value of u is undefinedat = which is a grid point. We simply take u(,) = in our code, and observe that the value of u(,) makes little influence on the final result in our numerical eperiment. We plot the results at t = with N = 4 and N = 8 cells in Figure 4, and can see that the scheme gives good results for this problem. 4
15 N=4 -. N=8 -.4 HWENO Eact -.4 HWENO Eact (a) (b) Figure 4: H(u) = (/4)(u )(u 4). t =. (a) N = 4 and (b) N = 8. Solid line: the eact solution; Square Symbol: the HWENO scheme. 3.3 Two dimensional case: accuracy test Eample 3.8: We solve the two dimensional linear scalar equation t + + y =,y with the initial data (,y,) = sin(π( + y)/) and periodic condition. The result at t = is listed in Table 4. We see that the scheme can achieve the designed accuracy. Table 4: t + + y =,(,y,) = sin(π(+y)/). Periodic boundary conditions. t = N L error L order L error L order L error L order.85e-.75e-.45e-.6e E E E E E E E E E E E E E E Eample 3.9: We solve the two dimensional Burgers equation t + ( + y +) =,y 5
16 with the initial data (,y,) = cos( π ( + y)) and periodic boundary condition. We compute the result up to t =.5/π and the solution is still smooth at that time. Again, from Table 5, we can see that the scheme achieves its designed order. Table 5: t + ( + y + ) =, (,y,) = cos( π ( + y)). Periodic boundary conditions. t =.5/π N L error L order L error L order L error L order.78e-3 9.6E-4 8.3E E E E E E E E E E E E E E E E- 5.9 Eample 3.: We solve t cos( + y +) =,y with the initial data (,) = cos( π (+y)) and periodic boundary condition. We compute the result up to t =.5/π and the solution is still smooth at that time. Again, from Table 6, we can see the scheme achieves the designed order of accuracy. Table 6: t cos( + y + ) =, (,y,) = cos( π ( + y)). Periodic boundary conditions. t =.5/π N L error L order L error L order L error L order.35e-3.5e E-4 3.E E E E E E E E E E E E E E E Two Dimensional case: discontinuous derivatives Eample 3.: We solve the two dimensional Burgers equation t + ( + y +) =,y 6
17 with the initial data (,y,) = cos( π ( + y)) and periodic boundary condition. We compute the result up to t =.5/π and the derivative discontinuity has appeared in the solution. We plot the results with 4 4 cells in Figure 5 and observe high resolution in this eample. Z Y X y y (a) (b) Figure 5: Two dimensional Burgers equation. t =.5/π with the HWENO schemes with 4 4 cells. Contours of the solution in Figure 5(a) and the surface of the solution in Figure 5(b) Eample 3.: We solve a problem from optimal control: t +sin(y) +(sin()+sign( y )) y sin (y)+cos() =, π <,y < π with (,y,) = and periodic boundary conditions. In this case, the equation can be denoted as t +H(, y,,y) =, and our scheme can be obtained through the following system of equations as before: d ij = H(, y,,y) dt du ij = H (u,v,,y)u H (u,v,,y)u y H (u,v,,y) dt dv ij = H (u,v,,y)v H (u,v,,y)v y H y (u,v,,y) dt Following (.3), we can easily deal with the terms H(, y,,y),h (u,v,,y),h (u,v,,y). The additional terms H (u,v,,y) and H y (u,v,,y) can be approimated as H (u,v,,y) =i,y=y j H ( u ij +u+ ij 7, v ij +v+ ij, i,y j )
18 H y (u,v,,y) =i,y=y j H y ( u ij +u+ ij, v ij +v+ ij, i,y j ) The time step is still taken as t = α. The solution and optimal control ω = sign() at.6 + β.6 y t = are plotted in Figure 6(a) and Figure 6(b), respectively. Again, we can observe that high resolution is achieved by our method. Z Z Y Y X X ω y 3 y (a) (b) Figure 6: Theoptimal control problem. t = with the HWENO schemes with 6 6 cells. Surface of the solution in Figure 6(a) and of the optimal control ω = sign( y ) in Figure 6(b) Eample 3.3: We solve the problem t cos( + y +) =, <,y < with (,y,) = cos( π (+y)) and periodic boundary conditions. We compute the solution up to t =.5/π and the solution has already developed discontinuous derivatives. The results are shown in Figure 7. We observe high resolutions for this eample. Eample 3.4: We solve the problem with another neither conve nor concave Hamiltonian t +sin( + y ) =, <,y < with (,y,) = π( y ). In this case, u(,y,) is undefined along = and v(,y,) is undefined along y =, and we simply take u(,y,) = and v(,,) = respectively. We 8
19 Z Y X y y (a) (b) Figure 7: Two dimensional equation with a neither conve nor concave Hamiltonian. t =.5/π with the HWENO schemes with 4 4 cells. Contours of the solution in Figure 7(a) and surface of the solution in Figure 7(b) compute the solution up to t =. The solution is shown in Figure 8. Again, we observe our schemes can achieve high resolution. Z Y X.5 y y (a) (b) Figure 8: Two dimensional equation with a neither conve nor concave Hamiltonian. t = by the HWENO schemes with 4 4 cells. Contours of the solution in Figure 8(a) and surface of the solution in Figure 8(b) 9
20 Eample 3.5: We solve the two dimensional Eikonal equation t + + y + =,,y < with the initial data (,y,) = 4 (cos(π) )(cos(πy) ). We compute the solution up to t =.6. The solution is shown in Figure 9. High resolutions are observed with our scheme. Z X Y y (a) y (b) Figure 9: Eikonal equation. t =.6 by the HWENO schemes with 6 6 cells. Contours of the solution in Figure 9(a) and surface of the solution in Figure 9(b) Eample 3.6: We solve t ( εk) + y + =,,y < where K is the mean curvature defined by: (,y,) = 4 (cos(π) )(cos(πy) ) K = (+ y ) y y + yy (+ ) (+ + y )3/ andεisasmallconstant, withtheinitial data(,y,) = 4 (cos(π) )(cos(πy) ) andperiodicboundarycondition. Whenε, theequationcanbedenotedas t +H(, y,, y, yy ) =
21 . The system of equations to be approimated then becomes d ij = H(, y,, y, yy ) dt du ij = H u H y u y H u H y u y H yy u yy dt dv ij = H v H y v y H v H y v y H yy v yy dt The procedure to deal with the terms H, H and H y are the same as before. Here, we simply use fourth order central differences in each direction with the other direction fied to approimate the terms u,u yy : u ij = u i+,j 6u i+,j +3u i,j 6u i,j +u i,j u yyij = u i,j+ 6u i,j+ +3u i,j 6u i,j +u i,j y As to the term u y, we take the following fourth order approimation using the values {u kl,k = i,,i+,l = j,,j +}: u y ij = 6 y (74u i+,j +63u i+,j +63u i,j+ +74u i,j+ 74u i+,j+ 63u i+,j+ 44u i,j+ +63u i,j+ 63u i+,j+ +44u i+,j+ +44u i,j 63u i,j +63u i+,j 44u i+,j 63u i,j 74u i,j ). The approimation to the terms v,v y,v yy can be obtained in a similar way. The time step is taken as t = α, where γ.6 + β.6 y + γ.3 + γ.3 y + γ 3 = ma H, γ = ma H y,.3 y y γ 3 = ma H yy and α, β are defined the same as before. The results of ε = (pure convection) yy and ε =. by the HWENO method with 6 6 cells are presented in Figure (a) and Figure (b) respectively. The surfaces at t = for ε = and for ε =., and at t =. for ε =., are shifted downward in order to show the details of the solution at later time. 4 Conclusion In this paper, we present a high order scheme based on the finite difference framework for the Hamilton-Jacobi equations in one and two dimensions. The main advantage of this scheme is
22 Z Z X Y X Y t=.9 t=.6.5 t=.6 t=.3.5 t=.3 t= t= -.35 t= y y (a) (b) Figure : Propagating surface. 6 6 cells. Figure (a): ε = ; Figure (b): ε =. its compactness and efficiency. Both the solution values and its first derivatives are used in the HWENO reconstruction and evolved via time marching. Etensive numerical eperiments in one dimensional and two dimensional cases show that the scheme can achieve high order accuracy in the smooth region and can maintain high resolution when the derivative becomes discontinuous. References [] R. Abgrall. Numerical discretization of the first-order Hamilton-Jacobi equation on triangular meshes. Commum. Pur. Appl. Math., 49: , 996. [] Y. Cheng and C.-W. Shu. A Discontinuous Galerkin finite element method for directly solving the Hamilton-Jacobi equations. J. Comput. Phys., 3:398 45, 7. [3] Y. Cheng and Z. Wang. A new Discontinuous Galerkin finite element method for directly solving the Hamilton-Jacobi equations. J. Comput. Phys., 68:34 53, 4.
23 [4] M. Crandall and P. L. Lions. Viscosity solutions of Hamilton-Jacobi equations. Trans. Americ. Math. Soc., 77: 4, 983. [5] M. Crandall and P. L. Lions. Two approimations of solutions of Hamilton-Jacobi equations. Trans. Americ. Math. Soc., 43: 9, 984. [6] C. Hu and C.-W. Shu. A Discontinuous Galerkin finite element method for Hamilton-Jacobi equations. SIAM J. Sci. Comput., :666 69, 999. [7] G. Jiang and D. Peng. Weighted ENO schemes for Hamilton-Jacobi equations. SIAM J. Sci. Comput., :6 43, 999. [8] G. S. Jiang and C.-W. Shu. Efficient implementation of weighted ENO schemes. J. Comput. Phys., 6: 8, 996. [9] F. Lafon and S. Osher. High order two dimensional nonoscillatory methods for solving Hamilton-Jacobi scalar equations. J. Comput. Phys., 3:35 53, 996. [] F. Li and C.-W. Shu. Reinterpretation and simplied implementation of a Discontinuous Galerkin method for Hamilton-Jacobi equations. Appl. Math. Lett., 8:4 9, 5. [] S. Osher and C.-W. Shu. High-order essentially nonoscillatory schemes for Hamilton-Jacobi equations. SIAM J. Numer. Anal, 8:97 9, 99. [] J. Qiu. Hermite WENO schemes with La-Wendroff type time discretizations for Hamilton- Jacobi equations. J. Comput. Math., 5:3 44, 7. [3] J. Qiu and C.-W. Shu. Hermite WENO schemes and their application as limiters for Runge- Kutta Discontinuous Galerkin method: one-dimensional case. J. Comput. Phys., 93:5 35, 4. [4] J. Qiu and C.-W. Shu. Hermite WENO schemes for Hamilton-Jacobi equations. J. Comput. Phys., 4:8 99, 5. 3
24 [5] C.-W. Shu. High order numerical methods for time dependent Hamilton-Jacobi equations. in Mathematics and Computation in Imaging Science and Information Processing, S.S. Goh, A. Ron and Z. Shen, Editors, Lecture Notes Series, Institute for Mathematical Sciences, National University of Singapore, World Scientific Press, :47 9, 7. [6] C.-W. Shu and S. Osher. Efficient implementation of essentially non-oscillatory shock capturing schemes. J. Comput. Phys., 77:439 47, 988. [7] Z. Tao and J. Qiu. Dimension-by-dimension moment-based central Hermite WENO schemes for directly solving Hamilton-Jacobi equations. Adv. Comput. Math., submitted. [8] J. Yan and S. Osher. A local Discontinuous Galerkin method for directly solving Hamilton- Jacobi equations. J. Comput. Phys., 3:3 44,. [9] Y. T. Zhang and C.-W. Shu. High order WENO schemes for Hamilton-Jacobi equations on triangular meshes. SIAM J. Sci. Comput, 4:5 3, 3. [] F. Zheng and J. Qiu. Directly solving the Hamilton-Jacobi equations by Hermite WENO schemes. J. Comput. Phys., 37:43 445, 6. [] J. Zhu and J. Qiu. Hermite WENO schemes for Hamilton-Jacobi equations on unstructured meshes. J. Comput. Phys., 54:76 9, 3. [] J. Zhu and J. Qiu. Finite volume Hermite WENO schemes for solving the Hamilton-Jacobi equation. Commun. Comput. Phys., 5:959 98, 4. [3] J. Zhu and J. Qiu. Finite volume Hermite WENO schemes for solving the Hamilton-Jacobi equations II: unstructured meshes. Computers Math. Appl., 68:37 5, 4. 4
High 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationarxiv: v1 [math.na] 2 Feb 2018
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,
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 informationFourier analysis for discontinuous Galerkin and related methods. Abstract
Fourier analysis for discontinuous Galerkin and related methods Mengping Zhang and Chi-Wang Shu Abstract In this paper we review a series of recent work on using a Fourier analysis technique to study the
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 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 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 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 informationTOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES
MATHEMATICS OF COMPUTATION Volume 67 Number 221 January 1998 Pages 73 85 S 0025-5718(98)00913-2 TOTAL VARIATION DIMINISHING RUNGE-KUTTA SCHEMES SIGAL GOTTLIEB AND CHI-WANG SHU Abstract. In this paper we
More informationA High Order WENO Scheme for a Hierarchical Size-Structured Model. Abstract
A High Order WENO Scheme for a Hierarchical Size-Structured Model Jun Shen 1, Chi-Wang Shu 2 and Mengping Zhang 3 Abstract In this paper we develop a high order explicit finite difference weighted essentially
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 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 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 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 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 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 informationHermite WENO schemes and their application as limiters for Runge Kutta discontinuous Galerkin method: one-dimensional case
Journal of Computational Physics 93 (3) 5 35 www.elsevier.com/locate/jcp Hermite WENO schemes and their application as limiters for Runge Kutta discontinuous Galerkin method: one-dimensional case Jianian
More informationAnti-diffusive finite difference WENO methods for shallow water with. transport of pollutant
Anti-diffusive finite difference WENO methods for shallow water with transport of pollutant Zhengfu Xu 1 and Chi-Wang Shu 2 Dedicated to Professor Qun Lin on the occasion of his 70th birthday Abstract
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 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 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 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 informationAccepted Manuscript. A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
Accepted Manuscript A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang, Hongkai Zhao PII: S-()- DOI:./j.jcp... Reference: YJCPH To
More informationDISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS
DISCONTINUOUS GALERKIN METHOD FOR TIME DEPENDENT PROBLEMS: SURVEY AND RECENT DEVELOPMENTS CHI-WANG SHU Abstract. In these lectures we give a general survey on discontinuous Galerkin methods for solving
More informationDownloaded 12/03/12 to Redistribution subject to SIAM license or copyright; see
SIAM J. SCI. COMPUT. Vol. 33, No., pp. 93 938 c Society for Industrial and Applied Mathematics Downloaded /3/ to 35..9.73. Redistribution subject to SIAM license or copyright; see http://www.siam.org/journals/ojsa.php
More informationA Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations
A Second Order Discontinuous Galerkin Fast Sweeping Method for Eikonal Equations Fengyan Li, Chi-Wang Shu, Yong-Tao Zhang 3, and Hongkai Zhao 4 ABSTRACT In this paper, we construct a second order fast
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 informationSemi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations
IMA Journal of Numerical Analysis 005 5, 113 138 doi: 10.1093/imanum/drh015 Semi-discrete central-upwind schemes with reduced dissipation for Hamilton Jacobi equations STEVE BRYSON Programme in Scientific
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 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 informationSimulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes
J Sci Comput (2011) 47: 281 302 DOI 10.1007/s10915-010-9437-3 Simulations of Shallow Water Equations with Finite Difference Lax-Wendroff Weighted Essentially Non-oscillatory Schemes Changna Lu Jianxian
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 Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws
J Sci Comput (011) 46: 359 378 DOI 10.1007/s10915-010-9407-9 A Speed-Up Strategy for Finite Volume WENO Schemes for Hyperbolic Conservation Laws Fei Teng Li Yuan Tao Tang Received: 3 February 010 / Revised:
More informationAdaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear Degenerate Reaction-Diffusion Problems
EPJ Web of Conferences 108, 0019 (016) DOI: 10.1051/ epjconf/ 0161080019 C Owned by the authors, published by EDP Sciences, 016 Adaptive WENO Schemes for Singular in Space and Time Solutions of Nonlinear
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 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 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 informationEulerian Vortex Motion in Two and Three Dimensions
Eulerian Vortex Motion in Two and Three Dimensions Igor Yanovsky June 2006 1 Objective An Eulerian approach is used to solve the motion of an incompressible fluid, in two and three dimensions, in which
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, College of Science, Nanjing University of Aeronautics and
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 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 informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
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 informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More 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 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 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 informationAN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION
AN OPTIMALLY ACCURATE SPECTRAL VOLUME FORMULATION WITH SYMMETRY PRESERVATION Fareed Hussain Mangi*, Umair Ali Khan**, Intesab Hussain Sadhayo**, Rameez Akbar Talani***, Asif Ali Memon* ABSTRACT High order
More informationCENTRAL DISCONTINUOUS GALERKIN METHODS ON OVERLAPPING CELLS WITH A NON-OSCILLATORY HIERARCHICAL RECONSTRUCTION
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 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 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 informationEntropy stable high order discontinuous Galerkin methods. for hyperbolic conservation laws
Entropy stable high order discontinuous Galerkin methods for hyperbolic conservation laws Chi-Wang Shu Division of Applied Mathematics Brown University Joint work with Tianheng Chen, and with Yong Liu
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 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 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 informationEssentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws
NASA/CR-97-0653 ICASE Report No. 97-65 Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws Chi-Wang Shu Brown University Institute for Computer
More informationA Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws
A Fourth-Order Central Runge-Kutta Scheme for Hyperbolic Conservation Laws Mehdi Dehghan, Rooholah Jazlanian Department of Applied Mathematics, Faculty of Mathematics and Computer Science, Amirkabir University
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 informationSOLUTION OF STEADY-STATE HAMILTON-JACOBI EQUATION BASED ON ALTERNATING EVOLUTION METHOD
ITALIAN JOURNAL OF PURE AND APPLIED MATHEMATICS N. 38 2017 (184 193 184 SOLUTION OF STEADY-STATE HAMILTON-JACOBI EQUATION BASED ON ALTERNATING EVOLUTION METHOD A. Tongxia Li Department of Mathematics Hulunbeier
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 informationThe Fast Sweeping Method. Hongkai Zhao. Department of Mathematics University of California, Irvine
The Fast Sweeping Method Hongkai Zhao Department of Mathematics University of California, Irvine Highlights of fast sweeping mehtod Simple (A Gauss-Seidel type of iterative method). Optimal complexity.
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 informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
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 informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
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 informationA genuinely high order total variation diminishing scheme for one-dimensional. scalar conservation laws 1. Abstract
A genuinely high order total variation diminishing scheme for one-dimensional scalar conservation laws Xiangxiong Zhang and Chi-Wang Shu 3 Abstract It is well known that finite difference or finite volume
More informationIMA Preprint Series # 1993
AN ADAPTIVE METHOD WITH RIGOROUS ERROR CONTROL FOR THE HAMILTON-JACOBI EQUATIONS. PART I: THE ONE-DIMENSIONAL STEADY STATE CASE By Bernardo Cockburn and Bayram Yenikaya IMA Preprint Series # 993 ( September
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 informationSolution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers
Solution of Two-Dimensional Riemann Problems for Gas Dynamics without Riemann Problem Solvers Alexander Kurganov, 1, * Eitan Tadmor 2 1 Department of Mathematics, University of Michigan, Ann Arbor, Michigan
More informationStrong Stability-Preserving (SSP) High-Order Time Discretization Methods
Strong Stability-Preserving (SSP) High-Order Time Discretization Methods Xinghui Zhong 12/09/ 2009 Outline 1 Introduction Why SSP methods Idea History/main reference 2 Explicit SSP Runge-Kutta Methods
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 informationMath 660-Lecture 23: Gudonov s method and some theories for FVM schemes
Math 660-Lecture 3: Gudonov s method and some theories for FVM schemes 1 The idea of FVM (You can refer to Chapter 4 in the book Finite volume methods for hyperbolic problems ) Consider the box [x 1/,
More informationA LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS
A LEVEL SET METHOD FOR THE COMPUTATION OF MULTIVALUED SOLUTIONS TO QUASI-LINEAR HYPERBOLIC PDES AND HAMILTON-JACOBI EQUATIONS SHI JIN AND STANLEY J. OSHER Abstract. We develop a level set method for the
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 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 informationcomputations. Furthermore, they ehibit strong superconvergence of DG solutions and flues for hyperbolic [,, 3] and elliptic [] problems. Adjerid et al
Superconvergence of Discontinuous Finite Element Solutions for Transient Convection-diffusion Problems Slimane Adjerid Department of Mathematics Virginia Polytechnic Institute and State University Blacksburg,
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 informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
More informationA Fifth Order Flux Implicit WENO Method
A Fifth Order Flux Implicit WENO Method Sigal Gottlieb and Julia S. Mullen and Steven J. Ruuth April 3, 25 Keywords: implicit, weighted essentially non-oscillatory, time-discretizations. Abstract The weighted
More informationMaximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments
Maximum-principle-satisfying and positivity-preserving high order schemes for conservation laws: Survey and new developments By Xiangxiong Zhang 1 and Chi-Wang Shu 1 Department of Mathematics, Brown University,
More informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
More information