Redistancing the Augmented Level Set Method
|
|
- Corey Mitchell
- 5 years ago
- Views:
Transcription
1 ILASS-Americas 23rd Annual Conference on Liquid Atomization and Spray Systems, Ventura, CA, May 2011 Redistancing the Augmented Level Set Method L. Anumolu and M.F. Trujillo Department of Mechanical Engineering University of Wisconsin Madison, Madison, WI Abstract Re-initialization of the level set function, through the solution of a Hamilton-Jacobi equation, has been found to lead to erroneous displacements of the zero level set iso-surface or interface. In the present work, we adopt an augmented-gradient methodology for solving this Hamilton-Jacobi equation, which updates both the level set function and its gradient. Even with this relatively accurate method, artificial displacements of the interface occur, leading to mass errors. Objects having coarse grid resolution are particularly vulnerable to this type of error. To address this issue, a hybrid scheme is proposed, where the underlying Hermite interpolation polynomials in nodes neighboring the interface are used to anchor the interface; thereby preventing erroneous movements. Away from this region the gradient-augmented solution for the re-distancing is employed. The results for 1D and 2D test cases show significant improvement with this hybrid treatment. Moreover, the introduction of a re-distancing scheme also improves the interface capturing behavior of the augmented level set method presented by Nave et al. [JCP 229:2010]. Corresponding Author: trujillo@engr.wisc.edu
2 Introduction The level set approach [1] has gained considerable popularity as a means of implicitly capturing the kinematics of a gas-liquid interface undergoing topological changes [2, 3]. Under this method, the level set function, ϕ, is generally defined as a signed distance function and the following transport equation is solved ϕ + u ϕ = 0, (1) t where, u, is the underlying flow field. We will interpret the ϕ = 0 iso-surface as denoting the material interface. A common error with the level set approach is the loss of mass that occurs due to erroneous displacements of the zero level set iso-surface, or the production of structures that are well below numerical resolution. One avenue that has been taken to avoid this issue is to combine the Volume-of- Fluid [4] and level set methods in the work of Sussman and Fatemi [5]. Other options include the use of Lagrangian particles within the level set method as reported by Enright et al. [6] and Wang et al. [7], or the use of a similar hybrid front tracking and level set method suggested by Nourgaliev et al. [8]. Improvements to the numerical accuracy for solving the ϕ transport equation, as in the gradient-augmented level set presented by Nave et al. [9], also contribute to the reduction of mass errors. Regardless of the method utilized, at some point the level set function will show significant distortion from its initial signed distance function characteristic, i.e ϕ will deviate substantially from one. In such cases small deviations in the level set function may cause large changes in the interface location, making the interface identification particularly susceptible to errors. It is therefore useful to recover the signed distance condition by replacing the deformed level set function with a new function without disturbing the zero level set [10]. This process is called reinitialization. This can be done by solving the Eikonal equation, ϕ = 1 using a fast marching method [11] or by recovering the property ϕ = 1 through the solution of the following Hamilton-Jacobi equation [2]. ϕ τ = sign(ϕo )(1 ϕ ), (2) where τ is the pseudo time for reinitialization and ϕ o is the initial deformed level set function. It can be shown that the analytical solution of this equation does not perturb the zero level iso-surface (i.e. ϕ o (x) = 0) as the re-initialization proceeds. In contrast, its numerical solution does resulting in artificial mass loss or gain. Since for practical problems, the only solution method is numerical, a need exists for preventing these artificial displacements. With this objective in mind, Sussman and Fatemi [12] proposed a scheme that anchors the interface by imposing local mass constraint. Russo and Smereka [13] suggested a simple modification by making use of truly upwind information for solving the reinitialization PDE (Eq. 2). Min [14] extended Russo s truly upwind method to higher order numerical accuracy by improving the calculation of gradients. In the present work, we present a hybrid treatment for re-initializing the level set function, which is based on a direct solution of the distance function for nodes neighboring the interface and an implementation of a gradient-augmented Hamilton-Jacobi method for nodes outside this region. This gradientaugmented approach was first introduced by van Leer [15] and has appeared in the solution of level set problems in the work of Nave et al. [9]. Here we apply this approach to the solution of the Hamilton- Jacobi equation. In the following section we describe our reinitialization scheme in this gradientaugmented framework and illustrate the artificial displacement of the interface that results when Eq. (2) is solved in the entire domain. This is followed by the presentation of our hybrid treatment, which is able to successfully address these artificial displacements. Various examples are used to illustrate the behavior. Gradient Augmented Level set Reinitialization In the work of Nave et al. [9], the characteristic form of the Courant-Isaacson-Rees (CIR) [16] method was used to solve for the level set function and its gradient, ϕ = ψ in a coupled fashion, without involving a re-initialization phase. This combined solution resulted in noticeable improvement in accuracy when compared to methods such as the 5th order WENO scheme [17]. Furthermore, the method is local and numerically stable. In our solution of the redistancing PDE, Eq. (2), we will employ a similar gradient-augmented approach. First, we take the gradient of Eq. (2) to yield an evolution equation for ϕ = ψ, namely ϕ τ +S(ϕo ) ϕ ϕ ψ = S(ϕo ) S(ϕ o ) ψ, (3) and we rewrite the reinitialization PDE (Eq. 2) as ϕ τ + S(ϕo ) ϕ ϕ ϕ = S(ϕo ). (4) 2
3 Putting both of these expressions in characteristic form, we have dϕ dτ = S(ϕo ), (5a) and dψ dτ = S(ϕo ) S(ϕ o ) ψ, (5b) along dx c dτ = S(ϕo ) ψ ψ, (5c) where S(ϕ o ) = sign(ϕ o ). We now proceed to solve these equations using the CIR methodology [16]. As an illustration of the procedure, consider a two-dimensional domain, where each node is denoted by indices (i, j). The field is described by ϕ n i,j and ψn i,j corresponding to time level τ n. Beginning with the initial field, we solve the characteristic equation Eq. (5c) backwards in time, such that at time level τ n + τ, each node (i, j) has been intercepted by its corresponding characteristic, i.e. x i,j = x c (τ n + τ). This is done using the following 3rd order Runge-Kutta method (RK) j 1 x (j) c = α jk x (k) c k=0 β jk τs(ϕ o ) ψ(k) ψ (k) ; j = 1, 2, 3, (6) where x c (τ n + τ) = x (0) c = x i,j and x c (τ n ) = x (3). The RK coefficients are represented by α jk and β jk. Here the origin of the characteristic is denoted by x c (τ n ) (= x f ). This location generally lies somewhere inside a given computational cell. In the following step, we integrate Eqs. (5a) and (5b) for ϕ and ψ, respectively from τ n to τ n + τ using the RK method. This is accomplish by solving the following equations This is obtained through fourth order Hermite interpolation, namely P (x c ) = 1 1 α=0 β=0 ϕi+α,j+β f α+1 (η x )f β+1 (η y ) β=2 ϕi+α,j+β 2 f α+1 (η x )f β+1 (η y ) y + 3 α=0 y 1 α=2 β=0 ϕi+α 2,j+β x + 3 α=2 f α+1 (η x )f β+1 (η y ) x 3 β=2 ϕi+α 2,j+β 2 xy f α+1 (η x )f β+1 (η y ) x y, (9) where P (x c ) represents either ϕ or ψ. When P (x c ) represents ψ then basis polynomials in Eq. (9) have to be differentiated with respect to x and y, and η x = x c x i x, η y = y c y j. y In the case of P (x c (τ n )) = ϕ(x c (τ n )), ϕ x, ϕ y and ϕ xy are the partial derivatives of ϕ, x and y are the respective cell spacing in the x and y directions, and basis polynomials are given by f β (η) = 1 3η 2 + 2η 3 if β = 1 3η 2 2η 3 if β = 2 η 2η 2 + η 3 if β = 3 η 2 + η 3 if β = 4 (10) The node values for ϕ in Eq. (9) are indicated by a superscript and their respective positions in a computation cell are illustrated in Fig. 1. j 1 ϕ (j) = α jk ϕ (k) + β jk τs (ϕ o ) ; j = 1, 2, 3. (7) k=0 and j 1 ψ (j) = α jk ψ (k) + β jk τ( S (ϕ o ) k=0 S (ϕ o ) ψ (k) ); j = 1, 2, 3. (8) resulting in the automatic update of the nodal values ϕ n+1 i,j and ψ n+1 i,j. This step requires that we determine the corresponding values for ϕ and ψ at x c (τ n ). Figure 1 Solution of Eqs. (6), (7) and (8) represents a relatively high accurate approach for re-initializing the level set function and its gradient. The number of pseudo-time steps employed is directly related to the extent of the re-initialization domain with respect to the interface. This is by virtue of the characteristics which propagate sign distance function behavior at velocities given by S(ϕ o ) ψ from the interface. We ψ 3
4 illustrate its performance with selected 1D and 2D test cases. The first test case is a 1D example taken from Czajkowski and Desjardins [18]. The initial field is given by function ϕ o = 0.05 sinh( 18x 5 ) and we computed re-initialization with 23 grid points spanning the entire domain. The results after re-initialization are shown globally in Fig. 2a and locally in Fig. 2b. These are contrasted to the exact solution. It is clear that re-initialization has effectively produced ϕ = 1 behavior everywhere. Nevertheless, as a result of this process the interface has erroneously shifted. This error is not surprising. As explained in the introduction, it has previously been discussed in the work of Sussman and Fatemi [12], Russo and Smereka [13] and Min [14]. atomization calculations, where the resulting small droplets are generally not sufficiently resolved due to computational restrictions. As suggested in the results of Fig. 3, it is this small droplet population that is more prone to drastic mass errors stemming from re-initialization using only the Hamilton-Jacobi equation. (a) (b) (a) Grid: 8x8. (b) Grid: 50x50. Figure 3: Reinitialization of a circular disk employing (a) 4 nodes and (b) 25 nodes across the disk diameter (Initial-Red, After reinitialization-using 200 time steps τ = x 2 Blue). The results clearly indicated that structures with low level of grid resolution are particularly susceptible to the interface drift error study in the present work. Figure 2: (a) Farfield and (b) nearfield view of the signed distance function generated using the gradient augmented scheme. A second case test consist of the re-initialization of the level set function beginning with the following 2D field ϕ o = (x x o ) 2 + (y y o ) 2 r. (11) where (x o, y o ) is the center and r is the radius of the disc. This level set field already satisfies the ϕ = 1 condition. This test case was implemented in the work of Sussman and Fatemi. [12], where their objective was to develop a reinitialization scheme to control the erroneous drift of the interface. In our calculations shown in Fig. 3, we employ two levels of grid resolution for the re-initialization process. The coarse and fine staggered mesh resolution corresponds to 4 (8x8 domain) and 25 (50x50 domain) grid cells spanning the diameter of the disk, respectively. In both cases artificial displacement of the interface occurs. However, this is much more noticeable for the coarser case; hence, objects with insufficient level of resolution are more susceptible to this type of error. This is particularly relevant in spray Hybrid Treatment From the test cases illustrated in the previous section, it is clear that the gradient augmented redistancing scheme Eqs. (5a), (5b) and (5c) fails to prevent erroneous interface displacement. In fact, any numerical solution of the Hamilton-Jacobi equation tends to incur this type of error as elaborated in [13, 12]. The underlying cause will be elucidated by analysis in a future publication. For the moment, we restrict the current work towards an immediate solution to this problem. Our approach is to directly solve for the shortest distance between all nodes neighboring the interface and the interface itself using the already existing Hermite interpolation polynomial basis. Once this has been done, the calculated distance is assigned to the corresponding level set function for each interface-neighboring nodes. This essentially anchors the level set field in a band around the interface. Subsequently, the Hamilton-Jacobi equation is solved for a prescribed time for the remaining nodes, where this prescribed time can be used as a way of controlling the extent of the re-initialized domain. While in theory the direct method applied to the interface-neighboring nodes can be equally applied to a larger number of nodes, in practice it is computationally more efficient 4
5 to employ the Hamilton-Jacobi method. Therefore, only the neighboring nodes, which are responsibe for the erroneous shift in the interface, are re-initialized with the direct method. The two methods combined provide a solution to the artificial interface motion, while minimizing computational expense. The first step in anchoring the neighboring nodes is their identification. This is done by solving the characteristic equation (5c) backwards in time to determine the origin (x f ) of each characteristic, which after τ intercepts its corresponding node. As an illustration, consider the cell depicted in Fig. 1, where the characteristic is shown originating from some point x f in the direction S (ϕ o ) ψ Since ψ. ϕ (x i,j+1 ) ϕ (x f ) < 0, node (i, j + 1) is a neighboring node. We now employ the Hermite interpolation basis along this characteristic line to solve for the shortest distance to the interface. Note that since this characteristic propagates normal to the interface, it already provides the direction of this shortest path. For the 2D geometry shown in Fig. 1, the Hermite basis polynomials are given by where a α (η x ) = and b α (η y ) = 1 3ηx 2 + 2ηx 3 if α = 1 3ηx 2 2ηx 3 if α = 2 η x 2ηx 2 + ηx 3 if α = 3 ηx 2 + ηx 3 if α = 4 1 3ηy 2 + 2ηy 3 if α = 1 3ηy 2 2ηy 3 if α = 2 η y 2ηy 2 + ηy 3 if α = 3 ηy 2 + ηy 3 if α = 4 η x = x x i x, η y = y y j y., (12), (13) Let us denote the coefficients of the above basis polynomials a α (η x ) and b α (η y ), respectively by A β α and Bα, β where β corresponds to the power of the corresponding term. For instance, for a 1 (η x ) = 1 3ηx 2 +2ηx 3 (α = 1), the coefficients A β α are A 3 1 = 2, A 2 1 = 3, A 1 1 = 0 and A 0 1 = 1. Along the characteristic direction, η y is a function of η x (η y = mη x + C), which can be substituted into b α (η y ). With this substitution, a α (η x ) and b α (η y ) can be introduced into the general expression for P (x) (Eq. 9) resulting in a sixth order polynomial with respect to η x. The coefficients of this polynomial can be expressed as C ζ = c ζ 11 ϕi,j + c ζ 12 ϕi,j+1 + c ζ 13 ϕi,j y +c ζ 21 ϕi+1,j + c ζ 22 ϕi+1,j+1 + c ζ 23 ϕi+1,j y y y + c ζ 14 ϕi,j+1 y y +c ζ 24 ϕi+1,j+1 y y + c ζ 31 ϕi,j x x + c ζ 32 ϕi,j+1 x x +c ζ 33 ϕi,j xy x y + c ζ 34 ϕi,j+1 xy x y + c ζ 41 ϕi+1,j x x +c ζ 42 ϕi+1,j+1 x x + c ζ 43 ϕi+1,j xy x y + c ζ 44 ϕi+1,j+1 xy x y, (14) where ζ [1, 7] N denotes (ζ 1) order for the corresponding term. Here, the values for c ζ γδ are given by c 1 γδ = A1 γb 1 δ, c 2 γδ = A1 γb 2 δ + A2 γb 1 δ, c 3 γδ = A1 γb 3 δ + A3 γb 1 δ + A2 γb 2 δ, c 4 γδ = A1 γb 4 δ + A2 γb 3 δ + A3 γb 2 δ + A4 γb 1 δ,, c 5 γδ = A2 γb 4 δ + A3 γb 3 δ + A4 γb 2 δ, c 6 γδ = A3 γb 4 δ + A4 γb 3 δ, and c 7 γδ = A4 γb 4 δ, with γ = (1, 2, 3, 4) and δ = (1, 2, 3, 4). Roots of the 6 th order polynomial are readily obtained from intrinsic functions in MATLAB. We discard all complex or imaginary roots, or roots outside the computational cell in question. Of the remaining valid roots, we choose the minimum one and assigned this to ϕ i,j+1. This process is repeated for all nodes neighboring the interface. Nodes that are farther removed from the interface are updated with the gradient augmented reinitialization scheme presented previously. We return to our former 1D case (i.e. ϕ o = 0.05 sinh(18x/5)), and apply the hybrid reinitialization. The results are shown in Fig. 4 which clearly mark the substantial improvement incurred by anchoring the interface. A more quantitative comparison is afforded through the use of the L 2 error norm, defined as L 2 = N (ϕ (x i ) ϕ exact (x i )) 2 (15) i=1 for N grid cells. In Table 1, the comparison between hybrid and Hamilton-Jacobi (conventional) re-initialization are compared for N=7, 20 and 49, which confirm the previous graphical findings. Furthermore the values clearly indicate that grid resolution plays a significant role in exacerbating reinitialization errors. 5
6 Grid cells Hybrid Conventional E E E E E E-2 Table 1: L 2 error comparison between the hybrid and conventional re-initialization scheme. (a) Conventional reinitialization scheme. (b) Hybrid reinitialization scheme. Figure 5: Mass errors are significantly reduced with the implementation of the hybrid re-initialization procedure. (a) Figure 4: (a) Farfield view of the comparison between hybrid and non-hybrid schemes and (b) near field view. Similarly, we return to the 2D case, corresponding to the ϕ o field given in Eq. (11). As discovered previously, the coarse grid re-initialization yielded unacceptable results. The hybrid treatment is now contrasted to these previous findings using the conventional re-initialization. The results are shown in Fig. 5, which again show significant improvement. Quantitative findings are determined with the L 1 error defined as Ω (b) 1 L H ( ϕ numerical) H ( ϕ expect) dxdy (16) where a numerical integration is performed following the procedure outlined by Sussman and Fatemi [12]. Their re-initialization procedure is compared to the one proposed here, and the results are shown in Table 2 after 4 and 8 steps of the re-initialization process. We report an improvement with the current method. x Steps L 1 Error L 1 Error [12] E E E E-04 Table 2: L 1 Error calculated from Eq. (16) and compared to result from [12], ( τ = x). Advection with re-initialization: In the original presentation of the augmentedlevel-set method, Nave et al. [9] argued that the method was sufficiently accurate for the omission of a re-initialization phase. While this may be true for some cases, the standard advection of the Zalesak disk [19] show noticeable improvement when exercising the hybrid re-initialization scheme. This is shown in Fig. 6. During hybrid re-initialization, the field was reinitialized for 6 pseudo time steps ( τ = x) every 314 advection steps. The case was executed in a domain of size 100x100, where the disk was centered at (50, 75) with a radius equal to 15. The width and height of the notch was 5 and 25, and the underlying flow field was given by u = π (100 y) 314 and v = π 314 (x 100). (17) As displayed in Fig. 6, the notch was preserved after 5 rotations, where it was lost in the pure advection case. To display the degree by which the level set function departs from ϕ = 1 condition, maps of this quantity are shown in Fig. 7 for both respective cases. Without re-initialization errors incurred during the advection of the interface lead to severe distortion of the level set function. Furthermore, errors originating from the boundary also shown signs of disturbing the ϕ = 1 condition. With proper re-initialization, these erroneous effects are greatly diminished. As a last test, we study the advection of a har- 6
7 (a) Advected for 1 rotation. (a) Advected with hybrid (b) Advected without reinireinitialization scheme. tialization. Figure 8: A = 4; H = 2 (Initial-Red, Advection with reinitialization-blue, Advection with out reinitialization-green). Figure 6: Zalesak disk advected for 5 rotations (Initial-Red, After 1 rotation-blue, After 4 rotations-black, After 5 rotations-green). (a) ψ for Fig (6a). Conclusions and Future work A semi-lagrangian method was presented for level set re-distancing. where the level set function ϕ and its gradient ϕ are updated simultaneously. Unlike the truly upwind schemes proposed in [13, 14] our method benefits from a local stencil and enjoys 4th order accuracy by virtue of the Hermite interpolation polynomials. The scheme presented addresses effectively the erroneous drift in the interface that plagues the conventional re-initialization procedure. This is accomplished by implementing a hybrid scheme, where the level set function is obtained directly from the Hermite polynomials for nodes neighboring the interface. This essentially anchors the interface and prevents artificial motion. For nodes away from this region, the re-distancing procedure employs a gradient-augmented methodology inspired by the advection scheme presented by Nave et al. [9]. In cases, where grid resolution is relatively coarse, re-initialization errors have the potential to dominate advection errors. It is in these cases, where the re-initialization scheme presented is truly beneficial. (b) ψ for Fig (6b). Figure 7: ψ at the near interface for 5th rotation of the Zalesak problem. monic disc, described by 2 2 (x xo ) + (y yo ) (r + A cos (4Hθ)) (18) where A = 4 is the harmonic amplitude and H = 2 are the number of harmonics per quarter circle. For this problem the domain size is 100x100, and the disk is centered at (50, 75) with a radius equal to 15. The flow field is the same as Zalesak s problem given in Eq. (17). Fig. 8 displays the results after 1 and 2 rotations corresponding to pure advection and hybrid re-initialization. In the case of re-initialization, the procedure was applied every quarter rotation resulting in a noticeable improvement. However, the improvement is not as drastic as the previous 1D case or circular disk reinitialization under coarse grid resolution. We find that for certain cases the error due to advection can easily overwhelm errors due to artificial interface displacement during re-initialization. Hence, even though these re-initialization errors are drastically minimized most of the error, originating from advection, dominates the solution. ϕo = (b) Advected for 2 rotations. Acknowledments Authors would like to acknowledge the support of the Office of Naval Research through Mark Spector, code 331, and the Caterpillar corporation. References [1] S. Osher and J. Sethian. J. Comput. Phys., 79:12 49, [2] M. Sussman, P. Smereka, and S. Osher. J. Comput. Phys., 114: , [3] M. Sussman, A.S. Almgren, J.B. Bell, P. Colella, and L.H. Howell. J. Comput. Phys., 148:81 124,
8 [4] C.W. Hirt and B.D. Nichols. J. Comput. Phys., 39: , [5] M. Sussman and E.G. Puckett. J. Comput. Phys., 162: , [6] D. Enright, R. Fedkiw, J. Ferziger, and I. Mitchell. J. Comput. Phys., 183:83 116, [7] Z. Wang, J. Yang, and F. Stern. J. Comput. Phys., 228: , [8] R. Nourgaliev, S. Kadioglu, and V. Mousseau. SIAM J. SCI. COMPUT., 32: , [9] J.-C. Nave, R. R. Rosales, and B. Seibold. J. Comput. Phys., 229: , [10] D. L. Chopp. J. Comput. Phys., 106:77 91, [11] J.A. Sethian. Proc. Nat. Acad. Sci, 93: , [12] M. Sussman and E. Fatemi. SIAM J. SCI. COMPUT., 20: , [13] G. Russo and P. Smereka. J. Comput. Phys., 163:51 67, [14] Chohong Min. J. Comput. Phys., 229: , [15] B. van Leer. J. Comput. Phys., 32: , [16] R. Courant, E. Isaacson, and M. Rees. Commun. Pure Appl. Math., 5: , [17] X.-D. Liu, S. Osher, and T. Chan. J. Comput. Phys., 115: , [18] M. Czajkowski and O. Desjardins. 22nd Annual Conference on Liquid Atomization and Spray Systems, Cincinnati, Ohio, May [19] S. T. Zalesak. J. Comput. Phys., 31: ,
NUMERICAL 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 informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationNumerical simulation of wave breaking in turbulent two-phase Couette flow
Center for Turbulence Research Annual Research Briefs 2012 171 Numerical simulation of wave breaking in turbulent two-phase Couette flow By D. Kim, A. Mani AND P. Moin 1. Motivation and objectives When
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationSTRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD
Shape optimization 1 G. Allaire STRUCTURAL OPTIMIZATION BY THE LEVEL SET METHOD G. ALLAIRE CMAP, Ecole Polytechnique with F. JOUVE and A.-M. TOADER 1. Introduction 2. Setting of the problem 3. Shape differentiation
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 informationHeat Transfer Equations The starting point is the conservation of mass, momentum and energy:
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 On Computational Investigation of the Supercooled Stefan Problem
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 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 informationA Discontinuous Galerkin Conservative Level Set Scheme for Simulating Turbulent Primary Atomization
ILASS-Americas 3rd Annual Conference on Liquid Atomization and Spray Systems, Ventura, CA, May 011 A Discontinuous Galerkin Conservative Level Set Scheme for Simulating Turbulent Primary Atomization M.
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 informationA PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS
A PARALLEL LEVEL-SET APPROACH FOR TWO-PHASE FLOW PROBLEMS WITH SURFACE TENSION IN THREE SPACE DIMENSIONS ROBERTO CROCE, MICHAEL GRIEBEL, AND MARC ALEXANDER SCHWEITZER Abstract. In this paper we present
More informationApplication of the immersed boundary method to simulate flows inside and outside the nozzles
Application of the immersed boundary method to simulate flows inside and outside the nozzles E. Noël, A. Berlemont, J. Cousin 1, T. Ménard UMR 6614 - CORIA, Université et INSA de Rouen, France emeline.noel@coria.fr,
More informationA MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD
Ninth International Conference on CFD in the Minerals and Process Industries CSIRO, Melbourne, Australia 10-12 December 2012 A MULTIPHASE FLUID-SOLID MODEL BASED ON THE LEVEL SET METHOD James E. HILTON
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 informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More information(1:1) 1. The gauge formulation of the Navier-Stokes equation We start with the incompressible Navier-Stokes equation 8 >< >: u t +(u r)u + rp = 1 Re 4
Gauge Finite Element Method for Incompressible Flows Weinan E 1 Courant Institute of Mathematical Sciences New York, NY 10012 Jian-Guo Liu 2 Temple University Philadelphia, PA 19122 Abstract: We present
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 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 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 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 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 informationComputational High Frequency Wave Propagation
Computational High Frequency Wave Propagation Olof Runborg CSC, KTH Isaac Newton Institute Cambridge, February 2007 Olof Runborg (KTH) High-Frequency Waves INI, 2007 1 / 52 Outline 1 Introduction, background
More informationDirect Numerical Simulation of Two-Phase Flows with phase change
Direct Numerical Simulation of Two-Phase Flows with phase change Sébastien Tanguy, Lucia Rueda Villegas, Catherine Colin, Grégory Huber Institut de Mécanique des Fluides de Toulouse Maître de conférence
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 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 informationDepartment of Mathematics University of California Santa Barbara, Santa Barbara, California,
The Ghost Fluid Method for Viscous Flows 1 Ronald P. Fedkiw Computer Science Department Stanford University, Stanford, California 9435 Email:fedkiw@cs.stanford.edu Xu-Dong Liu Department of Mathematics
More informationA level set projection model of lipid vesicles in general flows
, A level set projection model of lipid vesicles in general flows D. Salac a, M. Miksis a a Northwestern University, Engineering Sciences and Applied Mathematics, Evanston, IL, 60208 Abstract A new numerical
More informationA High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations
Motivation Numerical methods Numerical tests Conclusions A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations Xiaofeng Cai Department of Mathematics
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 informationChapter 1 Direct Modeling for Computational Fluid Dynamics
Chapter 1 Direct Modeling for Computational Fluid Dynamics Computational fluid dynamics (CFD) is a scientific discipline, which aims to capture fluid motion in a discretized space. The description of the
More informationA finite difference Poisson solver for irregular geometries
ANZIAM J. 45 (E) ppc713 C728, 2004 C713 A finite difference Poisson solver for irregular geometries Z. Jomaa C. Macaskill (Received 8 August 2003, revised 21 January 2004) Abstract The motivation for this
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 informationFinite Volume Method for Scalar Advection in 2D
Chapter 9 Finite Volume Method for Scalar Advection in D 9. Introduction The purpose of this exercise is to code a program to integrate the scalar advection equation in two-dimensional flows. 9. The D
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 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 information2 Introduction Standard conservative schemes will admit non-physical oscillations near material interfaces, which can be avoided by the application of
A Conservative Approach to the Multiphase Euler Equations without Spurious Pressure Oscillations Xu-Dong Liu Ronald P. Fedkiw Stanley Osher June 8, 998 Abstract In [4], the pressure evolution equation
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 informationTHE TOPOLOGICAL DESIGN OF MATERIALS WITH SPECIFIED THERMAL EXPANSION USING A LEVEL SET-BASED PARAMETERIZATION METHOD
11th. World Congress on Computational Mechanics (WCCM XI) 5th. European Conference on Computational Mechanics (ECCM V) 6th. European Conference on Computational Fluid Dynamics (ECFD VI) July 20-25, 2014,
More informationarxiv: v2 [physics.comp-ph] 4 Feb 2014
Fast and accurate solution of the Poisson equation in an immersed setting arxiv:1401.8084v2 [physics.comp-ph] 4 Feb 2014 Alexandre Noll Marques a, Jean-Christophe Nave b, Rodolfo Ruben Rosales c Abstract
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 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 informationModified Fourier Spectral Method for Non-periodic CFD
Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu Chaoqun Liu Technical Report 2011-14 http://www.uta.edu/math/preprint/ Modified Fourier Spectral Method for Non-periodic CFD Huankun Fu
More informationA modified quarter point element for fracture analysis of cracks
ndian Journal of Engineering & Materials Sciences Vol. 14, February 007, pp. 31-38 A modified quarter point element for fracture analysis of cracks Sayantan Paul & B N Rao* Structural Engineering Division,
More informationAn interface treating technique for compressible multi-medium flow with. Runge-Kutta discontinuous Galerkin method. Abstract
An interface treating technique for compressible multi-medium flow with Runge-Kutta discontinuous Galerkin method Chunwu Wang 1 and Chi-Wang Shu 2 Abstract The high-order accurate Runge-Kutta discontinuous
More informationConstruction of Shapes Arising from the Minkowski Problem Using a Level Set Approach
Construction of Shapes Arising from the Minkowski Problem Using a Level Set Approach Li-Tien Cheng June 15, 2006 Department of Mathematics University of California, San Diego 9500 Gilman Dr. La Jolla,
More informationA general well-balanced finite volume scheme for Euler equations with gravity
A general well-balanced finite volume scheme for Euler equations with gravity Jonas P. Berberich, Praveen Chandrashekar, Christian Klingenberg Abstract We present a second order well-balanced Godunov-type
More 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 informationSimulation of the evolution of concentrated shear layers in a Maxwell fluid with a fast high-resolution finite-difference scheme
J. Non-Newtonian Fluid Mech. 84 (1999) 275±287 Simulation of the evolution of concentrated shear layers in a Maxwell fluid with a fast high-resolution finite-difference scheme Raz Kupferman 1,a,*, Morton
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 informationA Level Set/Vortex Sheet Method for Modeling Phase Interface Dynamics During Primary Breakup
ILASS Americas 17th Annual Conference on Liquid Atomization and Spray Systems, Arlington, VA, May 2004 A Level Set/Vortex Sheet Method for Modeling Phase Interface Dynamics During Primary Breakup M. Herrmann
More informationA MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION
1 A MULTISCALE APPROACH IN TOPOLOGY OPTIMIZATION Grégoire ALLAIRE CMAP, Ecole Polytechnique The most recent results were obtained in collaboration with F. de Gournay, F. Jouve, O. Pantz, A.-M. Toader.
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 informationThree-dimensional high order large-scale numerical investigation on the air explosion
Three-dimensional high order large-scale numerical investigation on the air explosion Cheng Wang a, JianXu Ding a, Chi-Wang Shu b,, Tao Li a a State Key Laboratory of Explosion Science and Technology,
More informationLevel-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation
Level-Set Minimization of Potential Controlled Hadwiger Valuations for Molecular Solvation Bo Li Dept. of Math & NSF Center for Theoretical Biological Physics UC San Diego, USA Collaborators: Li-Tien Cheng
More informationChapter 1. Introduction
Chapter 1 Introduction Many astrophysical scenarios are modeled using the field equations of fluid dynamics. Fluids are generally challenging systems to describe analytically, as they form a nonlinear
More informationTotal Variation Theory and Its Applications
Total Variation Theory and Its Applications 2nd UCC Annual Research Conference, Kingston, Jamaica Peter Ndajah University of the Commonwealth Caribbean, Kingston, Jamaica September 27, 2018 Peter Ndajah
More informationNumerical Methods in Quantum Field Theories
Numerical Methods in Quantum Field Theories Christopher Bell 2011 NSF/REU Program Physics Department, University of Notre Dame Advisors: Antonio Delgado, Christopher Kolda 1 Abstract In this paper, preliminary
More informationEINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science. CASA-Report March2008
EINDHOVEN UNIVERSITY OF TECHNOLOGY Department of Mathematics and Computer Science CASA-Report 08-08 March2008 The complexe flux scheme for spherically symmetrie conservation laws by J.H.M. ten Thije Boonkkamp,
More informationNumerical Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 9 Finite Volume method II 9.1. Outline of Lecture Conservation property of Finite Volume method Apply FVM to
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
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 Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationModeling, Simulating and Rendering Fluids. Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan
Modeling, Simulating and Rendering Fluids Thanks to Ron Fediw et al, Jos Stam, Henrik Jensen, Ryan Applications Mostly Hollywood Shrek Antz Terminator 3 Many others Games Engineering Animating Fluids is
More informationAn Efficient, Interface-Preserving Level Set Re-distancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow
An Efficient, Interface-Preserving Level Set Re-distancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow Mark Sussman, Emad Fatemi To cite this version: Mark Sussman, Emad Fatemi.
More informationAN AUGMENTED IIM-LEVEL SET METHOD FOR STOKES EQUATIONS WITH DISCONTINUOUS VISCOSITY
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 193 210. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationarxiv: v2 [math.na] 29 Sep 2011
A Correction Function Method for Poisson Problems with Interface Jump Conditions. Alexandre Noll Marques a, Jean-Christophe Nave b, Rodolfo Ruben Rosales c arxiv:1010.0652v2 [math.na] 29 Sep 2011 a Department
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationPressure-correction algorithm to solve Poisson system with constant coefficients for fast two-phase simulations
Center or Turbulence Research Annual Research Bries 2009 367 Pressure-correction algorithm to solve Poisson system with constant coeicients or ast two-phase simulations By D. Kim AND P. Moin 1. Motivation
More informationGauge finite element method for incompressible flows
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2000; 34: 701 710 Gauge finite element method for incompressible flows Weinan E a, *,1 and Jian-Guo Liu b,2 a Courant Institute
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
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 informationSolving the Euler Equations!
http://www.nd.edu/~gtryggva/cfd-course/! Solving the Euler Equations! Grétar Tryggvason! Spring 0! The Euler equations for D flow:! where! Define! Ideal Gas:! ρ ρu ρu + ρu + p = 0 t x ( / ) ρe ρu E + p
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
More informationBy drawing Mohr s circle, the stress transformation in 2-D can be done graphically. + σ x σ y. cos 2θ + τ xy sin 2θ, (1) sin 2θ + τ xy cos 2θ.
Mohr s Circle By drawing Mohr s circle, the stress transformation in -D can be done graphically. σ = σ x + σ y τ = σ x σ y + σ x σ y cos θ + τ xy sin θ, 1 sin θ + τ xy cos θ. Note that the angle of rotation,
More informationc 2003 International Press
COMM. MATH. SCI. Vol., No. 3, pp. 593 6 c 3 International Press COMPUTATIONAL HIGH-FREQUENCY WAVE PROPAGATION USING THE LEVEL SET METHOD, WITH APPLICATIONS TO THE SEMI-CLASSICAL LIMIT OF SCHRÖDINGER EQUATIONS
More informationarxiv: v1 [physics.flu-dyn] 5 Aug 2017
On a relation between the volume of fluid, level-set and phase field interface models arxiv:1708.01805v1 [physics.flu-dyn] 5 Aug 2017 Abstract Tomasz Wac lawczyk, Department of Aeronautics, Institute of
More informationDirect Modeling for Computational Fluid Dynamics
Direct Modeling for Computational Fluid Dynamics Kun Xu February 20, 2013 Computational fluid dynamics (CFD) is new emerging scientific discipline, and targets to simulate fluid motion in different scales.
More informationExperimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column
Experimental and numerical study of the initial stages in the interaction process between a planar shock wave and a water column Dan Igra and Kazuyoshi Takayama Shock Wave Research Center, Institute of
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 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 informationOblique Drop Impact on Deep and Shallow Liquid
1 2 3 4 5 6 7 8 9 11 10 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Commun. Comput. Phys. doi: 10.4208/cicp.XXX.XXX Oblique Drop Impact on Deep and Shallow Liquid B. Ray 1, G. Biswas 1,2, and A. Sharma 3
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationDescription of the fire scheme in WRF
Description of the fire scheme in WRF March 8, 2010 1 Introduction The wildland fire model in WRF is an implementation of the semi-empirical fire propagation model developed in Coen (2005) and Clark et
More informationDUE: WEDS MARCH 26TH 2018
HOMEWORK # 2: FINITE DIFFERENCES MAPPING AND TWO DIMENSIONAL PROBLEMS DUE: WEDS MARCH 26TH 2018 NOTE: In this homework, you will choose ONE of the following three questions to perform and hand-in. Each
More informationDevelopment and analysis of a Lagrange-Remap sharp interface solver for stable and accurate atomization computations
ICLASS 2012, 12 th Triennial International Conference on Liquid Atomization and Spray Systems, Heidelberg, Germany, September 2-6, 2012 Development and analysis of a Lagrange-Remap sharp interface solver
More informationOptimal Interface Conditions for an Arbitrary Decomposition into Subdomains
Optimal Interface Conditions for an Arbitrary Decomposition into Subdomains Martin J. Gander and Felix Kwok Section de mathématiques, Université de Genève, Geneva CH-1211, Switzerland, Martin.Gander@unige.ch;
More informationNonlinear oscillations and rotations of a liquid droplet
Nonlinear oscillations and rotations of a liquid droplet Tadashi Watanabe watanabe.tadashi66@jaea.go.jp Abstract Nonlinear oscillations and rotations of a liquid droplet are simulated numerically by solving
More informationComputing High Frequency Waves By the Level Set Method
Computing High Frequency Waves By the Level Set Method Hailiang Liu Department of Mathematics Iowa State University Collaborators: Li-Tien Cheng (UCSD), Stanley Osher (UCLA) Shi Jin (UW-Madison), Richard
More informationSimulating Interfacial Tension of a Falling. Drop in a Moving Mesh Framework
Simulating Interfacial Tension of a Falling Drop in a Moving Mesh Framework Anja R. Paschedag a,, Blair Perot b a TU Berlin, Institute of Chemical Engineering, 10623 Berlin, Germany b University of Massachusetts,
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 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 informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationConvergence of a first order scheme for a non local eikonal equation
Convergence of a first order scheme for a non local eikonal equation O. Alvarez, E. Carlini, R. Monneau, E. Rouy March 22, 2005 Abstract We prove the convergence of a first order finite difference scheme
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 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 informationEulerian interface-sharpening methods for hyperbolic problems
Eulerian interface-sharpening methods for hyperbolic problems Application to compressible multiphase flow Keh-Ming Shyue Department of Mathematics National Taiwan University Taiwan 11:55-12:25, March 05,
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
More informationFinite volume method for two-phase flows using level set formulation
Finite volume method for two-phase flows using level set formulation Peter Frolkovič,a,1, Dmitry Logashenko b, Christian Wehner c, Gabriel Wittum c a Department of Mathematics and Descriptive Geometry,
More information