H(div) conforming and discontinuous Galerkin (DG) methods are designed for incompressible Euler's equation in two and three dimension. Error estimates are proved for both the semi-discrete method and fully-discrete method using backward Euler time stepping. Numerical examples exhibiting the performance of the methods are given.

No abstract uploaded!

We propose a computational efficient yet simple numerical algorithm to solve the surface eikonal equation on general implicit surfaces. The method is developed based on the embedding idea and the fast sweeping methods. We first approximate the solution to the surface eikonal equation by the Euclidean weighted distance function defined in a tubular neighbourhood of the implicit surface, and then apply the efficient fast sweeping method to numerically compute the corresponding viscosity solution. Unlike some other embedding methods which require the radius of the computational tube satisfies h = O ( x ) for some h = O ( x ) , our approach allows h = O ( x ) . This implies that the total number of grid points in the computational tube is optimal and is given by h = O ( x ) for a co-dimensional one surface in h = O ( x ) . The method can be easily extended to general static HamiltonJacobi equation defined on implicit

In our earlier work, we constructed uniformly high order accurate discontinuous Galerkin (DG) and finite volume schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics. In this paper, we present an extension of this framework to construct positivity-preserving high order essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) finite difference schemes for compressible Euler equations. General equations of state and source terms are also discussed. Numerical tests of the fifth order finite difference WENO scheme are reported to demonstrate the good behavior of such schemes.

No abstract uploaded!

The goal of this paper is to develop a novel semi-implicit spectral deferred correction (SDC) time marching method. The method can be used in a large class of problems, especially for highly nonlinear ordinary differential equations (ODEs) without easily separating of stiff and non-stiff components, which is more general and efficient comparing with traditional semi-implicit SDC methods. The proposed semi-implicit SDC method is based on low order time integration methods and corrected iteratively. The order of accuracy is increased for each additional iteration. And we also explore its local truncation error analytically. This SDC method is intended to be combined with the method of lines, which provides a flexible framework to develop high order semi-implicit time marching methods for nonlinear partial differential equations (PDEs). In this paper we mainly focus on the applications of the nonlinear PDEs with higher order spatial derivatives, e.g. convection diffusion equation, the surface diffusion and Willmore flow of graphs, the Cahn– Hilliard equation, the Cahn–Hilliard–Brinkman system and the phase field crystal equation. Coupled with the local discontinuous Galerkin (LDG) spatial discretization, the fully discrete schemes are all high order accurate in both space and time, and stable numerically with the time step proportional to the spatial mesh size. Numerical experiments are carried out to illustrate the accuracy and capability of the proposed semi-implicit SDC method.

No abstract uploaded!

Finding the maximum eigenvalue of a symmetric tensor is an important topic in tensor computation and numerical multilinear algebra. In this paper, we introduce a new class of structured tensors called W-tensors, which not only extends thewell-studied nonnegative tensors by allowing negative entries but also covers several important tensors arising naturally from spectral hypergraph theory.We then show that finding the maximum H-eigenvalue of an even-order symmetric W-tensor is equivalent to solving a structured semidefinite program and hence can be validated in polynomial time. This yields a highly efficient semidefinite program algorithm for computing themaximumH-eigenvalue of W-tensors and is based on a new structured sums-of-squares decomposition result for a nonnegative polynomial induced by W-tensors. Numerical experiments illustrate that the proposed algorithm can successfully find the maximum H-eigenvalue of W-tensors with dimension up to 10,000, subject to machine precision. As applications, we provide a polynomial time algorithm for computing the maximum H-eigenvalues of large-size Laplacian tensors of hyperstars and hypertrees, where the algorithm can be up to 13 times faster than the state-of-the-art numerical method introduced by Ng, Qi, and Zhou in 2009. Finally, we also show that the proposed algorithm can be used to test the copositivity of a multivariate form associated with symmetric extended Z-tensors,whose order may be even or odd.

Solitary waves bifurcated from edges of Bloch bands in two-dimensional periodic media are determined both analytically and numerically in the context of a two-dimensional nonlinear Schrödinger equation with a periodic potential. Using multiscale perturbation methods, the envelope equations of solitary waves near Bloch bands are analytically derived. These envelope equations reveal that solitary waves can bifurcate from edges of Bloch bands under either focusing or defocusing nonlinearity, depending on the signs of the second-order dispersion coefficients at the edge points. Interestingly, at edge points with two linearly independent Bloch modes, the envelope equations lead to a host of solitary wave structures, including reduced-symmetry solitons, dipole-array solitons, vortex-cell solitons, and so on—many of which have not been reported before to our knowledge. It is also shown analytically that the centers of envelope solutions can be positioned at only four possible locations at or between potential peaks. Numerically, families of these solitary waves are directly computed both near and far away from the band edges. Near the band edges, the numerical solutions spread over many lattice sites, and they fully agree with the analytical solutions obtained from the envelope equations. Far away from the band edges, solitary waves are strongly localized, with intensity and phase profiles characteristic of individual families.

In this paper we generalize a new type of limiters based on the weighted essentially non-oscillatory (WENO) finite volume methodology for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving nonlinear hyperbolic conservation laws, which were recently developed for structured meshes, to two-dimensional unstructured triangular meshes. The key idea of such limiters is to use the entire polynomials of the DG solutions from the troubled cell and its immediate neighboring cells, and then apply the classical WENO procedure to form a convex combination of these polynomials based on smoothness indicators and nonlinear weights, with suitable adjustments to guarantee conservation. The main advantage of this new limiter is its simplicity in implementation, especially for the unstructured meshes considered in this paper, as only information from immediate neighbors is needed and the usage of complicated geometric information of the meshes is largely avoided. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good performance of this procedure.

No abstract uploaded!

In this work, we investigate wave transmission property through a zero index metamaterial (ZIM) waveguide embedded with rectangular dielectric defects. We show that total reflection and total transmission (cloaking) can be achieved by adjusting the geometric sizes and/or permittivities of the defects. Our work provides another possibility of manipulating wave propagation through ZIM in addition to the widely studied dielectric defects with cylindrical geometries.

The objective of this paper is to design novel multi-layer neural network architectures for multiscale simulations of flows taking into account the observed data and physical modeling concepts. Our approaches use deep learning concepts combined with local multiscale model reduction methodologies to predict flow dynamics. Using reduced-order model concepts is important for constructing robust deep learning architectures since the reduced-order models provide fewer degrees of freedom. Flow dynamics can be thought of as multi-layer networks. More precisely, the solution (eg, pressures and saturations) at the time instant n+ 1 depends on the solution at the time instant n+ 1 and input parameters, such as permeability fields, forcing terms, and initial conditions. One can regard the solution as a multi-layer network, where each layer, in general, is a nonlinear forward map and the number of layers relates to the internal time steps. We will rely on rigorous model reduction concepts to define unknowns and connections for each layer. In each layer, our reduced-order models will provide a forward map, which will be modified (" trained") using available data. It is critical to use reduced-order models for this purpose, which will identify the regions of influence and the appropriate number of variables. Because of the lack of available data, the training will be supplemented with computational data as needed and the interpolation between data-rich and data-deficient models. We will also use deep learning algorithms to train the elements of the reduced model discrete system. We will present main ingredients of our approach and numerical results. Numerical

In this paper, we develop a discontinuous Galerkin (DG) method to solve the ideal special relativistic hydrodynamics (RHD) and design a bound-preserving (BP) limiter for this scheme by extending the idea in (X. Zhang and C.-W. Shu, Journal of Computational Physics, 229 (2010), 8918-8934). For RHD, the density and pressure are positive and the velocity is bounded by the speed of light. One difficulty in numerically solving the RHD in its conservative form is that the failure of preserving these physical bounds will result in ill-posedness of the problem and blowup of the code, especially in extreme relativistic cases. The standard way in dealing with this difficulty is to add extra numerical dissipation, while in doing so there is no guarantee in maintaining the high order of accuracy. Our BP limiter has the following features. It can theoretically guarantee to preserve the physical bounds for the numerical solution and maintain its designed high order accuracy. The limiter is local to the cell and hence is very easy to implement. Moreover, it renders $L^1$-stability to the numerical scheme. Numerical experiments are performed to demonstrate the good performance of this bound-preserving DG scheme. Even though we only discuss the BP limiter for DG schemes, it can be applied to high order finite volume schemes, such as weighted essentially non-oscillatory (WENO) finite volume schemes as well.

The immersed boundary method is one of the most useful computational methods in studying fluid structure interaction. On the other hand, the Immersed Boundary method is also known to require small time steps to maintain stability when solved with an explicit method. Many implicit or approximately implicit methods have been proposed in the literature to remove this severe time step stability constraint, but none of them give satisfactory performance. In this paper, we propose an efficient semi-implicit scheme to remove this stiffness from the immersed boundary method for the Navier–Stokes equations. The construction of our semi-implicit scheme consists of two steps. First, we obtain a semiimplicit discretization which is proved to be unconditionally stable. This unconditionally stable semi-implicit scheme is still quite expensive to implement in practice. Next, we apply the small scale decomposition to the unconditionally stable semi-implicit scheme to construct our efficient semi-implicit scheme. Unlike other implicit or semi-implicit schemes proposed in the literature, our semi-implicit scheme can be solved explicitly in the spectral space. Thus the computational cost of our semi-implicit schemes is comparable to that of an explicit scheme. Our extensive numerical experiments show that our semiimplicit scheme has much better stability property than an explicit scheme. This offers a substantial computational saving in using the immersed boundary method.

No abstract uploaded!

In this paper we propose a time-domain finite element method for modeling of electromagnetic cloaks. The permittivity and permeability of the cloak model are described by the Drude dispersion model. The model to be solved is quite challenging in that we have to solve a coupled problem with different partial differential equations given in different regions. Our method is based on a mixed finite element method using edge elements with different types of meshes in different regions. Numerical results demonstrate that our algorithm is quite effective for simulating cloaks in time-domain. To our knowledge, this is the first cloak simulation carried out by the time-domain finite element method.

One of the main challenges in computational simulations of gas detonation propagation is that negative density or negative pressure may emerge during the time evolution, which will cause blow-ups. Therefore, schemes with provable positivity-preserving of density and pressure are desired. First order and second order positivity-preserving schemes were well studied. For high order discontinuous Galerkin (DG) method, even though the characteristicwise TVB limiter can kill oscillations, it is not sufficient to maintain the positivity. A simple solution for arbitrarily high order positivity-preserving schemes solving Euler equations was proposed recently. In this paper, we first discuss an extension of the technique to design arbitrarily high order positivity-preserving DG schemes for reactive Euler equations. We then present a simpler and more robust implementation of the positivity-preserving limiter than the one before. Numerical tests, including very demanding examples in gaseous detonations, indicate that the third order DG scheme with the new positivity-preserving limiter produces satisfying results even without the TVB limiter.

Previous work has shown that in a two-dimensional periodic medium under focusing or defocusing cubic nonlinearities, gap solitons in the form of low-amplitude and slowly modulated single-Bloch-wave packets can bifurcate out from the edges of Bloch bands. In this paper, linear stability properties of these gap solitons near band edges are determined both analytically and numerically. Through asymptotic analysis, it is shown that these gap solitons are linearly unstable if the slope of their power curve at the band edge has the opposite sign of nonlinearity
here focusing nonlinearity is said to have a positive sign, and defocusing nonlinearity to have a negative sign. An equivalent condition for linear instability is that the power of the gap solitons near the band edge is lower than the limit power value on the band edge. Through numerical computations of the power curves, it is found that this condition is always satisfied, thus two-dimensional gap solitons near band edges are linearly unstable. The analytical formula for the unstable eigenvalue of gap solitons near band edges is also asymptotically derived. It is shown that this unstable eigenvalue is proportional to the cubic power of the soliton’s amplitude, and it induces width instabilities of gap solitons. A comparison between this analytical eigenvalue formula and numerically computed eigenvalues shows excellent agreement.

For decades, the widely used finite difference method on staggered grids, also known as the marker and cell (MAC) method, has been one of the simplest and most effective numerical schemes for solving the Stokes equations and Navier–Stokes equations. Its superconvergence on uniform meshes has been observed by Nicolaides (SIAM J Numer Anal 29(6):1579–1591, 1992), but the rigorous proof is never given. Its behavior on non-uniform grids is not well studied, since most publications only consider uniform grids. In this work, we develop the MAC scheme on non-uniform rectangular meshes, and for the first time we theoretically prove that the superconvergence phenomenon (i.e., second order convergence in the L2 norm for both velocity and pressure) holds true for the MAC method on non-uniform rectangular meshes. With a careful and accurate analysis of various sources of errors, we observe that even though the local truncation errors are only first order in terms of mesh size, the global errors after summation are second order due to the amazing cancellation of local errors. This observation leads to the elegant superconvergence analysis even with non-uniform meshes. Numerical results are given to verify our theoretical analysis.

In this paper, we propose a structured bisection method with adaptive randomized sampling for finding selected or all of the eigenvalues of certain real symmetric matrices A. For A with a low-rank property, we construct a hierarchically semiseparable (HSS) approximation and show how to quickly evaluate and update its inertia in the bisection method. Unlike some existing randomized HSS constructions, the methods here do not require the knowledge of the off-diagonal (numerical) ranks in advance. Moreover, for A with a weak rank property or slowly decaying off-diagonal singular values, we show an idea for aggressive low-rank inertia evaluation, which means that a compact HSS approximation can preserve the inertia for certain shifts. This is analytically justified for a special case, and numerically shown for more general ones. A generalized LDL factorization of the HSS approximation is then designed for the fast evaluation of the inertia. A significant advantage over standard LDL factorizations is that the HSS LDL factorization (and thus the inertia) of A − sI can be quickly updated with multiple shifts s in bisection. The factorization with each new shift can reuse about 60% of the work. As an important application, the structured eigensolver can be applied to symmetric Toeplitz matrices, and the cost to find one eigenvalue is nearly linear in the order of the matrix. The numerical examples demonstrate the efficiency and the accuracy of our methods, especially the benefit of low-rank inertia evaluations. The ideas and methods can be potentially adapted to other HSS computations where shifts are involved and to more problems without a significant low-rank property.

A convex non-convex variational model is proposed for multiphase image segmentation. We consider a specially designed non-convex regularization term which adapts spatially to the image structures for better controlling of the segmentation and easy handling of the intensity inhomogeneities. The nonlinear optimization problem is effciently solved by an Alternating Directions Methods of Multipliers procedure. We provide a convergence analysis and perform numerical experiments on several images, showing the effectiveness of this procedure.

Wave propagation problems arise in a wide range of applications. The energy conserving property is one of the guiding principles for numerical algorithms, in order to minimize the phase or shape errors after long time integration. In this paper, we develop and analyze a local discontinuous Galerkin (LDG) method for solving the wave equation. We prove optimal error estimates, superconvergence toward a particular projection of the exact solution, and the energy conserving property for the semi-discrete formulation. The analysis is extended to the fully discrete LDG scheme, with the centered second-order time discretization (the leap-frog scheme). Our numerical experiments demonstrate optimal rates of convergence and superconvergence. We also show that the shape of the solution, after long time integration, is well preserved due to the energy conserving property.

We present an efficient null space free Jacobi–Davidson method to compute the positive eigenvalues of time harmonic Maxwell’s equations. We focus on a class of spatial discretizations that guarantee the existence of discrete vector potentials, such as Yee’s scheme and the edge elements. During the Jacobi–Davidson iteration, the correction process is applied to the vector potential instead. The correction equation is solved approximately as in the standard Jacobi–Davidson approach. The computational cost of the transformation from the vector potential to the corrector is negligible. As a consequence, the expanding subspace automatically stays out of the null space and no extra projection step is needed. Numerical evidence confirms that the proposed scheme indeed outperforms the standard and projection-based Jacobi–Davidson methods by a significant margin.

We consider the numerical solution of the time-fractional diffusion-wave equation on two-dimensional and three-dimensional unbounded spatial domains. Introduce an artificial boundary and find the exact and approximate artificial boundary conditions for the given problem, which lead to a bounded computational domain. Using the exact or approximating boundary conditions on the artificial boundary, the original problem is reduced to an initial-boundary-value problem on the bounded computational domain which is respectively equivalent to or approximates the original problem. Finite difference methods are used to solve the reduced problems on the bounded computational domain and the stability of these finite difference methods are proved. The numerical results demonstrate that the method given in this paper is effective and feasible.