Sparse-Grid Implementation of Fixed-Point Fast Sweeping WENO Schemes for Eikonal Equations

Fixed-point fast sweeping methods are a class of explicit iterative methods developed in the literature to efficiently solve steady-state solutions of hyperbolic partial differential equations(PDEs).As other types of fast sweeping schemes,fixed-point fast sweeping methods use the Gauss-Seidel iterations and alternating sweeping strategy to cover characteristics of hyperbolic PDEs in a certain direction simultaneously in each sweeping order.The resulting iterative schemes have a fast convergence rate to steady-state solutions.Moreover,an advantage of fixed-point fast sweeping methods over other types of fast sweeping methods is that they are explicit and do not involve the inverse operation of any nonlinear local system.Hence,they are robust and flexible,and have been combined with high-order accurate weighted essentially non-oscillatory(WENO)schemes to solve various hyperbolic PDEs in the literature.For multidimensional nonlinear problems,high-order fixed-point fast sweeping WENO methods still require quite a large amount of computational costs.In this technical note,we apply sparse-grid techniques,an effective approximation tool for multidimensional problems,to fixed-point fast sweeping WENO methods for reducing their computational costs.Here,we focus on fixed-point fast sweeping WENO schemes with third-order accuracy(Zhang et al.2006[41]),for solving Eikonal equations,an important class of static Hamilton-Jacobi(H-J)equations.Numerical experiments on solving multidimensional Eikonal equations and a more general static H-J equation are performed to show that the sparse-grid computations of the fixed-point fast sweeping WENO schemes achieve large savings of CPU times on refined meshes,and at the same time maintain comparable accuracy and resolution with those on corresponding regular single grids.

Pseudo S-Asymptotically(ω,c)-Periodic Solutions to Fractional Differential Equations of Sobolev Type

In this paper,we firstly recall some basic results on pseudo S-asymptotically(ω,c)-periodic functions and Sobolev type fractional differential equation.We secondly investigate some existence of pseudo S-asymptotically(ω,c)-periodic solutions for a semilinear fractional differential equations of Sobolev type.We finally present a simple example.

The Regularity of Solutions to Mixed Boundary Value Problems of Second-Order Elliptic Equations with Small Angles

This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of sufficient regularity of boundary conditions and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order.

THE STABILITY OF BOUSSINESQ EQUATIONS WITH PARTIAL DISSIPATION AROUND THE HYDROSTATIC BALANCE

This paper is devoted to understanding the stability of perturbations around the hydrostatic equilibrium of the Boussinesq system in order to gain insight into certain atmospheric and oceanographic phenomena.The Boussinesq system focused on here is anisotropic,and involves only horizontal dissipation and thermal damping.In the 2D case R^(2),due to the lack of vertical dissipation,the stability and large-time behavior problems have remained open in a Sobolev setting.For the spatial domain T×R,this paper solves the stability problem and gives the precise large-time behavior of the perturbation.By decomposing the velocity u and temperatureθinto the horizontal average(ū,θ)and the corresponding oscillation(ū,θ),we can derive the global stability in H~2 and the exponential decay of(ū,θ)to zero in H^(1).Moreover,we also obtain that(ū_(2),θ)decays exponentially to zero in H^(1),and thatū_(1)decays exponentially toū_(1)(∞)in H^(1)as well;this reflects a strongly stratified phenomenon of buoyancy-driven fluids.In addition,we establish the global stability in H^(3)for the 3D case R^(3).

THE GLOBAL EXISTENCE AND ANALYTICITY OF A MILD SOLUTION TO THE 3D REGULARIZED MHD EQUATIONS

In this paper,we study the three-dimensional regularized MHD equations with fractional Laplacians in the dissipative and diffusive terms.We establish the global existence of mild solutions to this system with small initial data.In addition,we also obtain the Gevrey class regularity and the temporal decay rate of the solution.

THE RIEMANN PROBLEM FOR ISENTROPIC COMPRESSIBLE EULER EQUATIONS WITH DISCONTINUOUS FLUX

We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.

Some Modified Equations of the Sine-Hilbert Type

Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.

On entire solutions of some Fermat type differential-difference equations

On one hand,we study the existence of transcendental entire solutions with finite order of the Fermat type difference equations.On the other hand,we also investigate the existence and growth of solutions of nonlinear differential-difference equations.These results extend and improve some previous in[5,14].

THE ASYMPTOTIC BEHAVIOR AND OSCILLATION FOR A CLASS OF THIRD-ORDER NONLINEAR DELAY DYNAMIC EQUATIONS

In this paper,we consider a class of third-order nonlinear delay dynamic equations.First,we establish a Kiguradze-type lemma and some useful estimates.Second,we give a sufficient and necessary condition for the existence of eventually positive solutions having upper bounds and tending to zero.Third,we obtain new oscillation criteria by employing the Potzsche chain rule.Then,using the generalized Riccati transformation technique and averaging method,we establish the Philos-type oscillation criteria.Surprisingly,the integral value of the Philos-type oscillation criteria,which guarantees that all unbounded solutions oscillate,is greater than θ_(4)(t_(1),T).The results of Theorem 3.5 and Remark 3.6 are novel.Finally,we offer four examples to illustrate our results.

THE GLOBAL EXISTENCE OF STRONG SOLUTIONS TO THERMOMECHANICAL CUCKER-SMALE-STOKES EQUATIONS IN THE WHOLE DOMAIN

We study the global existence and uniqueness of a strong solution to the kinetic thermomechanical Cucker-Smale(for short,TCS) model coupled with Stokes equations in the whole space.The coupled system consists of the kinetic TCS equation for a particle ensemble and the Stokes equations for a fluid via a drag force.In this paper,we present a complete analysis of the existence of global-in-time strong solutions to the coupled model without any smallness restrictions on the initial data.

THE EXACT MEROMORPHIC SOLUTIONS OF SOME NONLINEAR DIFFERENTIAL EQUATIONS

We find the exact forms of meromorphic solutions of the nonlinear differential equations■,n≥3,k≥1,where q,Q are nonzero polynomials,Q■Const.,and p_(1),p_(2),α_(1),α_(2)are nonzero constants withα_(1)≠α_(2).Compared with previous results on the equation p(z)f^(3)+q(z)f"=-sinα(z)with polynomial coefficients,our results show that the coefficient of the term f^((k))perturbed by multiplying an exponential function will affect the structure of its solutions.

THE SPARSE REPRESENTATION RELATED WITH FRACTIONAL HEAT EQUATIONS

This study introduces a pre-orthogonal adaptive Fourier decomposition(POAFD)to obtain approximations and numerical solutions to the fractional Laplacian initial value problem and the extension problem of Caffarelli and Silvestre(generalized Poisson equation).As a first step,the method expands the initial data function into a sparse series of the fundamental solutions with fast convergence,and,as a second step,makes use of the semigroup or the reproducing kernel property of each of the expanding entries.Experiments show the effectiveness and efficiency of the proposed series solutions.

Approximate solution of Volterra-Fredholm integral equations using generalized barycentric rational interpolant

It is well-known that interpolation by rational functions results in a more accurate approximation than the polynomials interpolation.However,classical rational interpolation has some deficiencies such as uncontrollable poles and low convergence order.In contrast with the classical rational interpolants,the generalized barycentric rational interpolants which depend linearly on the interpolated values,yield infinite smooth approximation with no poles in real numbers.In this paper,a numerical collocation approach,based on the generalized barycentric rational interpolation and Gaussian quadrature formula,was introduced to approximate the solution of Volterra-Fredholm integral equations.Three types of points in the solution domain are used as interpolation nodes.The obtained numerical results confirm that the barycentric rational interpolants are efficient tools for solving Volterra-Fredholm integral equations.Moreover,integral equations with Runge’s function as an exact solution,no oscillation occurrs in the obtained approximate solutions so that the Runge’s phenomenon is avoided.

Results Involving Partial Differential Equations and Their Solution by Certain Integral Transform

In this study,we aimto investigate certain triple integral transformand its application to a class of partial differentialequations.We discuss various properties of the new transformincluding inversion, linearity, existence, scaling andshifting, etc. Then,we derive several results enfolding partial derivatives and establish amulti-convolution theorem.Further, we apply the aforementioned transform to some classical functions and many types of partial differentialequations involving heat equations,wave equations, Laplace equations, and Poisson equations aswell.Moreover,wedraw some figures to illustrate 3-D contour plots for exact solutions of some selected examples involving differentvalues in their variables.

THE OPTIMAL LARGE TIME BEHAVIOR OF3D QUASILINEAR HYPERBOLIC EQUATIONS WITH NONLINEAR DAMPING

We are concerned with the large-time behavior of 3D quasilinear hyperbolic equations with nonlinear damping.The main novelty of this paper is two-fold.First,we prove the optimal decay rates of the second and third order spatial derivatives of the solution,which are the same as those of the heat equation,and in particular,are faster than ones of previous related works.Second,for well-chosen initial data,we also show that the lower optimal L^(2) convergence rate of the k(∈[0,3])-order spatial derivatives of the solution is(1+t)^(-(2+2k)/4).Therefore,our decay rates are optimal in this sense.The proofs are based on the Fourier splitting method,low-frequency and high-frequency decomposition,and delicate energy estimates.

Meta-Auto-Decoder:a Meta-Learning-Based Reduced Order Model for Solving Parametric Partial Differential Equations

Many important problems in science and engineering require solving the so-called parametric partial differential equations(PDEs),i.e.,PDEs with different physical parameters,boundary conditions,shapes of computational domains,etc.Typical reduced order modeling techniques accelerate the solution of the parametric PDEs by projecting them onto a linear trial manifold constructed in the ofline stage.These methods often need a predefined mesh as well as a series of precomputed solution snapshots,and may struggle to balance between the efficiency and accuracy due to the limitation of the linear ansatz.Utilizing the nonlinear representation of neural networks(NNs),we propose the Meta-Auto-Decoder(MAD)to construct a nonlinear trial manifold,whose best possible performance is measured theoretically by the decoder width.Based on the meta-learning concept,the trial manifold can be learned in a mesh-free and unsupervised way during the pre-training stage.Fast adaptation to new(possibly heterogeneous)PDE parameters is enabled by searching on this trial manifold,and optionally fine-tuning the trial manifold at the same time.Extensive numerical experiments show that the MAD method exhibits a faster convergence speed without losing the accuracy than other deep learning-based methods.

Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.

A Provable Positivity-Preserving Local Discontinuous Galerkin Method for the Viscous and Resistive MHD Equations

In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.

Factorized Smith Method for A Class of High-Ranked Large-ScaleТ-Stein Equations

We introduce a factorized Smith method(FSM)for solving large-scale highranked J-Stein equations within the banded-plus-low-rank structure framework.To effectively reduce both computational complexity and storage requirements,we develop techniques including deflation and shift,partial truncation and compression,as well as redesign the residual computation and termination condition.Numerical examples demonstrate that the FSM outperforms the Smith method implemented with a hierarchical HODLR structured toolkit in terms of CPU time.

Local existence and uniqueness of the incompressible Navier-Stokes-Landau-Lifshitz equations with the Dzyaloshinskii-Moriya interaction and V-flow term

In this paper,we prove that there exists a unique local solution for the Cauchy problem of a system of the incompressible Navier-Stokes-Landau-Lifshitz equations with the Dzyaloshinskii-Moriya interaction and V-flow term inR^(2) and R^(3).Our methods rely upon approximating the system with a perturbed parabolic system and parallel transport.