An Eulerian hybrid WENO centered-difference solver for elastic–plastic solids

We present a finite-difference based solver for hyper-elastic and viscoplastic systems using a hybrid of the weighted essentially non-oscillatory (WENO) schemes combined with explicit centered difference to solve the equations of motion expressed in an Eulerian formulation. By construction our approach minimizes both numerical dissipation errors and the creation of curl-constraint violating errors away from discontinuities while avoiding the calculation of hyperbolic characteristics often needed in general finite-volume schemes. As a result of the latter feature, the formulation allows for a wide range of constitutive relations and only an upper-bound on the speed of sound at each time is required to ensure a stable timestep is chosen. Several one- and two-dimensional examples are presented using a range of constitutive laws with and without additional plastic modeling. In addition we extend the reflection technique combined with ghost-cells to enforce fixed boundaries with a zero tangential stress condition (i.e. free-slip).     

Unconditionally convergent nonlinear solver for hyperbolic conservation laws with S-shaped flux functions

This paper addresses the convergence properties of implicit numerical solution algorithms for nonlinear hyperbolic transport problems. It is shown that the Newton–Raphson (NR) method converges for any time step size, if the flux function is convex, concave, or linear, which is, in general, the case for CFD problems. In some problems, e.g., multiphase flow in porous media, the nonlinear flux function is S-shaped (not uniformly convex or concave); as a result, a standard NR iteration can diverge for large time steps, even if an implicit discretization scheme is used to solve the nonlinear system of equations. In practice, when such convergence difficulties are encountered, the current time step is cut, previous iterations are discarded, a smaller time step size is tried, and the NR process is repeated. The criteria for time step cutting and selection are usually based on heuristics that limit the allowable change in the solution over a time step and/or NR iteration. Here, we propose a simple modification to the NR iteration scheme for conservation laws with S-shaped flux functions that converges for any time step size. The new scheme allows one to choose the time step size based on accuracy consideration only without worrying about the convergence behavior of the nonlinear solver. The proposed method can be implemented in an existing simulator, e.g., for CO₂ sequestration or reservoir flow modeling, quite easily. The numerical analysis is confirmed with simulation studies using various test cases of nonlinear multiphase transport in porous media. The analysis and numerical experiments demonstrate that the modified scheme allows for the use of arbitrarily large time steps for this class of problems.     

Entropy viscosity method for nonlinear conservation laws

A new class of high-order numerical methods for approximating nonlinear conservation laws is described (entropy viscosity method). The novelty is that a nonlinear viscosity based on the local size of an entropy production is added to the numerical discretization at hand. This new approach does not use any flux or slope limiters, applies to equations or systems supplemented with one or more entropy inequalities and does not depend on the mesh type and polynomial approximation. Various benchmark problems are solved with finite elements, spectral elements and Fourier series to illustrate the capability of the proposed method.     

A high-order accurate hybrid scheme using a central flux scheme and a WENO scheme for compressible flowfield analysis

A high-order accurate hybrid central-WENO scheme is proposed. The fifth order WENO scheme [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] is divided into two parts, a central flux part and a numerical dissipation part, and is coupled with a central flux scheme. Two sub-schemes, the WENO scheme and the central flux scheme, are hybridized by means of a weighting function that indicates the local smoothness of the flowfields. The derived hybrid central-WENO scheme is written as a combination of the central flux scheme and the numerical dissipation of the fifth order WENO scheme, which is controlled adaptively by a weighting function. The structure of the proposed hybrid central-WENO scheme is similar to that of the YSD-type filter scheme [H.C. Yee, N.D. Sandham, M.J. Djomehri, Low-dissipative high-order shock-capturing methods using characteristic-based filters, J. Comput. Phys. 150 (1999) 199–238]. Therefore, the proposed hybrid scheme has also certain merits that the YSD-type filter scheme has. The accuracy and efficiency of the developed hybrid central-WENO scheme are investigated through numerical experiments on inviscid and viscous problems. Numerical results show that the proposed hybrid central-WENO scheme can resolve flow features extremely well.     

Compact third-order limiter functions for finite volume methods

We consider finite volume methods for the numerical solution of conservation laws. In order to achieve high-order accurate numerical approximation to non-linear smooth functions, we introduce a new class of limiter functions for the spatial reconstruction of hyperbolic equations. We therefore employ and generalize the idea of double-logarithmic reconstruction of Artebrant and Schroll [R. Artebrant, H.J. Schroll, Limiter-free third order logarithmic reconstruction, SIAM J. Sci. Comput. 28 (2006) 359-381].     

Chiral anomaly in Weyl systems: No violation of classical conservation laws

The anomalous term $～\mathbf{E}\mathbf{B}$ in the balance of the chiral density can be rewritten as quantum current in the classical balance of density. Therefore it does not violate conservation laws as sometimes claimed to be caused by quantum fluctuations.     

Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics

In a pair of earlier papers the author showed the importance of divergence-free reconstruction in adaptive mesh refinement problems for magnetohydrodynamics (MHD) and the importance of the same for designing robust second order schemes for MHD. Second order accurate divergence-free schemes for MHD have shown themselves to be very useful in several areas of science and engineering. However, certain computational MHD problems would be much benefited if the schemes had third and higher orders of accuracy. In this paper we show that the reconstruction of divergence-free vector fields can be carried out with better than second order accuracy. As a result, we design divergence-free weighted essentially non-oscillatory (WENO) schemes for MHD that have order of accuracy better than second. A multi-stage Runge–Kutta time integration is used to ensure that the temporal accuracy matches the spatial accuracy. While this is achieved quite simply up to third order in time, going beyond third order is most simply achieved by using the ADER-WENO schemes that are detailed in a companion paper. (ADER stands for Arbitrary Derivative Riemann Problem.) Accuracy analysis is carried out and it is shown that the schemes meet their design accuracy for smooth problems. Stringent tests are also presented showing that the schemes perform well on those tests.     

Shock waves in virus fitness evolution

We consider a nonlinear partial differential equation of conservation type to describe the dynamics of transmission with sampling of viral populations observed in aliquots of fixed particle number taken from an evolving clone at periodic intervals of time [I.S. Novella, E.A. Duarte, S.F. Elena, A. Moya, E. Domingo, Exponential increases of RNA virus fitness during large population transmissions, Proc. Natl. Acad. Sci. USA 92 (1995) 5841–5844]. The fast processes of mutation and natural selection that occur within each transmission are incorporated into the conservation equation as parameters. With this, the changes in time behavior of fitness function noticed in experimental data are related to a crossover exhibited by the solutions to this equation in the transient regime for pulse-like initial conditions. As a consequence, the average replication rate of the population is predicted to reach a plateau as a power t^-1/2$t^{-1/2}$.     

An efficient local time-stepping scheme for solution of nonlinear conservation laws

We develop an efficient local time-stepping algorithm for the method of lines approach to numerical solution of transient partial differential equations. The need for local time-stepping arises when adaptive mesh refinement results in a mesh containing cells of greatly different sizes. The global CFL number and, hence, the global time step, are defined by the smallest cell size. This can be inefficient as a few small cells may impose a restrictive time step on the whole mesh. A local time-stepping scheme allows us to use the local CFL number which reduces the total number of function evaluations. The algorithm is based on a second order Runge–Kutta time integration. Its important features are a small stencil and the second order accuracy in the L² and L^∞ norms.     

Dark optical solitons and conservation laws for parabolic and dual-power law nonlinearities in (2 + 1)-dimensions

This paper studies the dynamics of optical solitons with parabolic and dual-power law nonlinearities. The dark 1-soliton solution is first obtained by the ansatz method along with the necessary constraint conditions, for both of these nonlinearities. Subsequently, the invariance, conservation laws and double reductions of the governing nonlinear Schrödinger's equation are studied and the conserved densities are thus revealed.     

A class of large time step Godunov schemes for hyperbolic conservation laws and applications

A large time step (LTS) Godunov scheme firstly proposed by LeVeque is further developed in the present work and applied to Euler equations. Based on the analysis of the computational performances of LeVeque’s linear approximation on wave interactions, a multi-wave approximation on rarefaction fan is proposed to avoid the occurrences of rarefaction shocks in computations. The developed LTS scheme is validated using 1-D test cases, manifesting high resolution for discontinuities and the capability of maintaining computational stability when large CFL numbers are imposed. The scheme is then extended to multidimensional problems using dimensional splitting technique; the treatment of boundary condition for this multidimensional LTS scheme is also proposed. As for demonstration problems, inviscid flows over NACA0012 airfoil and ONERA M6 wing with given swept angle are simulated using the developed LTS scheme. The numerical results reveal the high resolution nature of the scheme, where the shock can be captured within 1–2 grid points. The resolution of the scheme would improve gradually along with the increasing of CFL number under an upper bound where the solution becomes severely oscillating across the shock. Computational efficiency comparisons show that the developed scheme is capable of reducing the computational time effectively with increasing the time step (CFL number).     

Bäcklund transformation,infinite number of conservation laws and fission properties of an integro-differential model for ocean internal solitary waves

This paper presents an analytical investigation of the propagation of internal solitary waves in the ocean of finite depth. Using the multi-scale analysis and reduced perturbation methods, the integrodifferential equation is derived, which is called the intermediate long wave(ILW) equation and can describe the amplitude of internal solitary waves. It can reduce to the Benjamin–Ono equation in the deep-water limit, and to the Kd V equation in the shallow-water limit. Little attention has been paid to the features of integro-differential equations, especially for their conservation laws. Here,based on Hirota bilinear method, B?cklund transformations in bilinear form of ILW equation are derived and infinite number of conservation laws are given. Finally, we analyze the fission phenomenon of internal solitary waves theoretically and verify it through numerical simulation. All of these have potential value for the further research on ocean internal solitary waves.     

Delta shock waves in conservation laws with impulsive moving source: some results obtained by multiplying distributions

The present paper concerns the study of a Riemann problem for the conservation law u_t + [? (u)]_x = kδ (x ? vt) where x, t, k, v and u = u(x,t) are real numbers. We consider ? an entire function taking real values on the real axis and δ stands for the Dirac measure. Within a convenient space of distributions we will explicitly see the possible emergence of waves with the shape of shock waves, delta waves and delta shock waves. For this purpose, we define a rigorous concept of a solution which extends both the classical solution concept and a weak solution concept. All this framework is developed in the setting of a distributional product that is not constructed by approximation. We include the main ideas of this product for the reader’s convenience. Recall that delta shock waves are relevant physical phenomena which may be interpreted as processes of concentration of mass or even as processes of formation of galaxies in the universe.     

A massively parallel multi-block hybrid compact–WENO scheme for compressible flows

A new multi-block hybrid compact–WENO finite-difference method for the massively parallel computation of compressible flows is presented. In contrast to earlier methods, our approach breaks the global dependence of compact methods by using explicit finite-difference methods at block interfaces and is fully conservative. The resulting method is fifth- and sixth-order accurate for the convective and diffusive fluxes, respectively. The impact of the explicit interface treatment on the stability and accuracy of the multi-block method is quantified for the advection and diffusion equations. Numerical errors increase slightly as the number of blocks is increased. It is also found that the maximum allowable time steps increase with the number of blocks. The method demonstrates excellent scalability on up to 1264 processors.     

A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids

Recently a new high-order formulation for 1D conservation laws was developed by Huynh using the idea of “flux reconstruction”. The formulation was capable of unifying several popular methods including the discontinuous Galerkin, staggered-grid multi-domain method, or the spectral difference/spectral volume methods into a single family. The extension of the method to quadrilateral and hexahedral elements is straightforward. In an attempt to extend the method to other element types such as triangular, tetrahedral or prismatic elements, the idea of “flux reconstruction” is generalized into a “lifting collocation penalty” approach. With a judicious selection of solution points and flux points, the approach can be made simple and efficient to implement for mixed grids. In addition, the formulation includes the discontinuous Galerkin, spectral volume and spectral difference methods as special cases. Several test problems are presented to demonstrate the capability of the method.     

A systematic methodology for constructing high-order energy stable WENO schemes

A third-order Energy Stable Weighted Essentially Non-Oscillatory (ESWENO) finite difference scheme developed by the authors of this paper [N.K. Yamaleev, M.H. Carpenter, Third-order energy stable WENO scheme, J. Comput. Phys. 228 (2009) 3025–3047] was proven to be stable in the energy norm for both continuous and discontinuous solutions of systems of linear hyperbolic equations. Herein, a systematic approach is presented that enables “energy stable” modifications for existing WENO schemes of any order. The technique is demonstrated by developing a one-parameter family of fifth-order upwind-biased ESWENO schemes including one sixth-order central scheme; ESWENO schemes up to eighth order are presented in the Appendix. We also develop new weight functions and derive constraints on their parameters, which provide consistency, much faster convergence of the high-order ESWENO schemes to their underlying linear schemes for smooth solutions with arbitrary number of vanishing derivatives, and better resolution near strong discontinuities than the conventional counterparts.     

Modulation instability,rogue waves and conservation laws in higher-order nonlinear Schrödinger equation

In this paper, the modulation instability(MI), rogue waves(RWs) and conservation laws of the coupled higher-order nonlinear Schr?dinger equation are investigated. According to MI and the2?×?2 Lax pair, Darboux-dressing transformation with an asymptotic expansion method, the existence and properties of the one-, second-, and third-order RWs for the higher-order nonlinear Schr?dinger equation are constructed. In addition, the main characteristics of these solutions are discussed through some graphics, which are draw widespread attention in a variety of complex systems such as optics, Bose–Einstein condensates, capillary flow, superfluidity, fluid dynamics,and finance. In addition, infinitely-many conservation laws are established.     

On symmetries,reductions, conservation laws and conserved quantities of optical solitons with inter-modal dispersion

This paper studies the compressional dispersive Alfvén (CDA) waves where Noether symmetries will be calculated from which the corresponding conservation laws will be obtained via Noether's theorem. Furthermore, one case of double reduction is performed via the association of a conserved vector with a Noether symmetry (with zero gauge). The conserved quantities of optical solitons in the presence of intermodal dispersion that is governed by the perturbed nonlinear Schrödinger's equation with Kerr law nonlinearity. The invariance-multiplier method is adopted to carry out the analysis, from which the conserved densities are then retrieved. Finally, the conserved quantities are obtained using the 1-soliton solution of the governing equation.