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.

AENO:a Novel Reconstruction Method in Conjunction with ADER Schemes for Hyperbolic Equations

In this paper,we present a novel spatial reconstruction scheme,called AENO,that results from a special averaging of the ENO polynomial and its closest neighbour,while retaining the stencil direction decided by the ENO choice.A variant of the scheme,called m-AENO,results from averaging the modified ENO(m-ENO)polynomial and its closest neighbour.The concept is thoroughly assessed for the one-dimensional linear advection equation and for a one-dimensional non-linear hyperbolic system,in conjunction with the fully discrete,high-order ADER approach implemented up to fifth order of accuracy in both space and time.The results,as compared to the conventional ENO,m-ENO and WENO schemes,are very encouraging.Surprisingly,our results show that the L_(1)-errors of the novel AENO approach are the smallest for most cases considered.Crucially,for a chosen error size,AENO turns out to be the most efficient method of all five methods tested.

Efficient Sparse-Grid Implementation of a Fifth-Order Multi-resolution WENO Scheme for Hyperbolic Equations

High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.

A Combination of Residual Distribution and the Active Flux Formulations or a New Class of Schemes That Can Combine Several Writings of the Same Hyperbolic Problem:Application to the 1D Euler Equations

We show how to combine in a natural way(i.e.,without any test nor switch)the conservative and non-conservative formulations of an hyperbolic system that has a conservative form.This is inspired from two different classes of schemes:the residual distribution one(Abgrall in Commun Appl Math Comput 2(3):341–368,2020),and the active flux formulations(Eyman and Roe in 49th AIAA Aerospace Science Meeting,2011;Eyman in active flux.PhD thesis,University of Michigan,2013;Helzel et al.in J Sci Comput 80(3):35–61,2019;Barsukow in J Sci Comput 86(1):paper No.3,34,2021;Roe in J Sci Comput 73:1094–1114,2017).The solution is globally continuous,and as in the active flux method,described by a combination of point values and average values.Unlike the“classical”active flux methods,the meaning of the point-wise and cell average degrees of freedom is different,and hence follow different forms of PDEs;it is a conservative version of the cell average,and a possibly non-conservative one for the points.This new class of scheme is proved to satisfy a Lax-Wendroff-like theorem.We also develop a method to perform nonlinear stability.We illustrate the behaviour on several benchmarks,some quite challenging.

A Posteriori Stabilized Sixth-Order Finite Volume Scheme with Adaptive Stencil Construction:Basics for the 1D Steady-State Hyperbolic Equations

We propose an adaptive stencil construction for high-order accurate finite volume schemes a posteriori stabilized devoted to solve one-dimensional steady-state hyperbolic equations.High accuracy(up to the sixth-order presently)is achieved,thanks to polynomial recon-structions while stability is provided with an a posteriori MOOD method which controls the cell polynomial degree for eliminating non-physical oscillations in the vicinity of dis-continuities.We supplemented this scheme with a stencil construction allowing to reduce even further the numerical dissipation.The stencil is shifted away from troubles(shocks,discontinuities,etc.)leading to less oscillating polynomial reconstructions.Experimented on linear,Burgers',and Euler equations,we demonstrate that the adaptive stencil technique manages to retrieve smooth solutions with optimal order of accuracy but also irregular ones without spurious oscillations.Moreover,we numerically show that the approach allows to reduce the dissipation still maintaining the essentially non-oscillatory behavior.

Actively tuning anisotropic light-matter interaction in biaxial hyperbolic materialα-MoO_(3) using phase change material VO_(2) and graphene

Anisotropic hyperbolic phonon polaritons(PhPs)in natural biaxial hyperbolic materialα-MoO_(3) has opened up new avenues for mid-infrared nanophotonics,while active tunability ofα-MoO_(3) PhPs is still an urgent problem necessarily to be solved.In this study,we present a theoretical demonstration of actively tuningα-MoO_(3) PhPs using phase change material VO_(2) and graphene.It is observed thatα-MoO_(3) PhPs are greatly dependent on the propagation plane angle of PhPs.The insulator-to-metal phase transition of VO_(2) has a significant effect on the hybridization PhPs of theα-MoO_(3)/VO_(2) structure and allows to obtain actively tunableα-MoO_(3) PhPs,which is especially obvious when the propagation plane angle of PhPs is 900.Moreover,when graphene surface plasmon sources are placed at the top or bottom ofα-MoO_(3) inα-MoO_(3)/VO_(2)structure,tunable coupled hyperbolic plasmon-phonon polaritons inside its Reststrahlen bands(RB s)and surface plasmonphonon polaritons outside its RBs can be achieved.In addition,the above-mentionedα-MoO_(3)-based structures also lead to actively tunable anisotropic spontaneous emission(SE)enhancement.This study may be beneficial for realization of active tunability of both PhPs and SE ofα-MoO_(3),and facilitate a deeper understanding of the mechanisms of anisotropic light-matter interaction inα-MoO_(3) using functional materials.

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.

Influence of substrate effect on near-field radiative modulator based on biaxial hyperbolic materials

Relative rotation between the emitter and receiver could effectively modulate the near-field radiative heat transfer(NFRHT)in anisotropic media.Due to the strong in-plane anisotropy,natural hyperbolic materials can be used to construct near-field radiative modulators with excellent modulation effects.However,in practical applications,natural hyperbolic materials need to be deposited on the substrate,and the influence of substrate on modulation effect has not been studied yet.In this work,we investigate the influence of substrate effect on near-field radiative modulator based onα-MoO_(3).The results show that compared to the situation without a substrate,the presence of both lossless and lossy substrate will reduce the modulation contrast(MC)for different film thicknesses.When the real or imaginary component of the substrate permittivity increases,the mismatch of hyperbolic phonon polaritons(HPPs)weakens,resulting in a reduction in MC.By reducing the real and imaginary components of substrate permittivity,the MC can be significantly improved,reaching 4.64 forε_(s)=3 at t=10 nm.This work indicates that choosing a substrate with a smaller permittivity helps to achieve a better modulation effect,and provides guidance for the application of natural hyperbolic materials in the near-field radiative modulator.

Dimension by Dimension Finite Volume HWENO Method for Hyperbolic Conservation Laws

In this paper,we propose a finite volume Hermite weighted essentially non-oscillatory(HWENO)method based on the dimension by dimension framework to solve hyperbolic conservation laws.It can maintain the high accuracy in the smooth region and obtain the high resolution solution when the discontinuity appears,and it is compact which will be good for giving the numerical boundary conditions.Furthermore,it avoids complicated least square procedure when we implement the genuine two dimensional(2D)finite volume HWENO reconstruction,and it can be regarded as a generalization of the one dimensional(1D)HWENO method.Extensive numerical tests are performed to verify the high resolution and high accuracy of the scheme.

An Extended Numerical Method by Stancu Polynomials for Solution of Integro-Differential Equations Arising in Oscillating Magnetic Fields

In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.

STARLIKENESS ASSOCIATED WITH THE SINE HYPERBOLIC FUNCTION

Let q_(λ)(z)=1+λsinh(ζ),0<λ<1/sinh(1)be a non-vanishing analytic function in the open unit disk.We introduce a subclass S^(*)(q_(λ))of starlike functions which contains the functions f such that zf'/f is subordinated by q_(λ).We establish inclusion and radii results for the class S^(*)(q_(λ))for several known classes of starlike functions.Furthermore,we obtain sharp coefficient bounds and sharp Hankel determinants of order two for the class S^(*)(q_(λ)).We also find a sharp bound for the third Hankel determinant for the caseλ=1/2.

Theoretical study of particle and energy balance equations in locally bounded plasmas

In this study,new particle and energy balance equations have been developed to predict the electron temperature and density in locally bounded plasmas.Classical particle and energy balance equations assume that all plasma within a reactor is completely confined only by the reactor walls.However,in industrial plasma reactors for semiconductor manufacturing,the plasma is partially confined by internal reactor structures.We predict the effect of the open boundary area(A′_(L,eff))and ion escape velocity(u_(i))on electron temperature and density by developing new particle and energy balance equations.Theoretically,we found a low ion escape velocity(u_(i)/u_(B)≈0.2)and high open boundary area(A′_(L,eff)/A_(T,eff)≈0.6)to result in an approximately 38%increase in electron density and an 8%decrease in electron temperature compared to values in a fully bounded reactor.Additionally,we suggest that the velocity of ions passing through the open boundary should exceedω_(pi)λ_(De)under the condition E^(2)_(0)?(Φ/λ_(De))^(2).

Efficient Finite Difference WENO Scheme for Hyperbolic Systems withNon-conservativeProducts

Higher order finite difference weighted essentially non-oscillatory(WENO)schemes have been constructed for conservation laws.For multidimensional problems,they offer a high order accuracy at a fraction of the cost of a finite volume WENO or DG scheme of the comparable accuracy.This makes them quite attractive for several science and engineering applications.But,to the best of our knowledge,such schemes have not been extended to non-linear hyperbolic systems with non-conservative products.In this paper,we perform such an extension which improves the domain of the applicability of such schemes.The extension is carried out by writing the scheme in fluctuation form.We use the HLLI Riemann solver of Dumbser and Balsara(J.Comput.Phys.304:275-319,2016)as a building block for carrying out this extension.Because of the use of an HLL building block,the resulting scheme has a proper supersonic limit.The use of anti-diffusive fluxes ensures that stationary discontinuities can be preserved by the scheme,thus expanding its domain of the applicability.Our new finite difference WENO formulation uses the same WENO reconstruction that was used in classical versions,making it very easy for users to transition over to the present formulation.For conservation laws,the new finite difference WENO is shown to perform as well as the classical version of finite difference WENO,with two major advantages:(i)It can capture jumps in stationary linearly degenerate wave families exactly.(i)It only requires the reconstruction to be applied once.Several examples from hyperbolic PDE systems with non-conservative products are shown which indicate that the scheme works and achieves its design order of the accuracy for smooth multidimensional flows.Stringent Riemann problems and several novel multidimensional problems that are drawn from compressible Baer-Nunziato multiphase flow,multiphase debris flow and twolayer shallow water equations are also shown to document the robustness of the method.For some test problems that require well-balancing we have even been able to apply the scheme without any modification and obtain good results.Many useful PDEs may have stiff relaxation source terms for which the finite difference formulation of WENO is shown to provide some genuine advantages.

Existence of Entropy Solution for Degenerate Parabolic-Hyperbolic Problem Involving p(x)-Laplacian with Neumann Boundary Condition

We consider a strongly non-linear degenerate parabolic-hyperbolic problem with p(x)-Laplacian diffusion flux function. We propose an entropy formulation and prove the existence of an entropy solution.

A New Efcient Explicit Deferred Correction Framework:Analysis and Applications to Hyperbolic PDEs and Adaptivity

The deferred correction(DeC)is an iterative procedure,characterized by increasing the accuracy at each iteration,which can be used to design numerical methods for systems of ODEs.The main advantage of such framework is the automatic way of getting arbitrarily high order methods,which can be put in the Runge-Kutta(RK)form.The drawback is the larger computational cost with respect to the most used RK methods.To reduce such cost,in an explicit setting,we propose an efcient modifcation:we introduce interpolation processes between the DeC iterations,decreasing the computational cost associated to the low order ones.We provide the Butcher tableaux of the new modifed methods and we study their stability,showing that in some cases the computational advantage does not afect the stability.The fexibility of the novel modifcation allows nontrivial applications to PDEs and construction of adaptive methods.The good performances of the introduced methods are broadly tested on several benchmarks both in ODE and PDE contexts.

Wave propagation of a functionally graded plate via integral variables with a hyperbolic arcsine function

Several studies on functionally graded materials(FGMs)have been done by researchers,but few studies have dealt with the impact of the modification of the properties of materials with regard to the functional propagation of the waves in plates.This work aims to explore the effects of changing compositional characteristics and the volume fraction of the constituent of plate materials regarding the wave propagation response of thick plates of FGM.This model is based on a higher-order theory and a new displacement field with four unknowns that introduce indeterminate integral variables with a hyperbolic arcsine function.The FGM plate is assumed to consist of a mixture of metal and ceramic,and its properties change depending on the power functions of the thickness of the plate,such as linear,quadratic,cubic,and inverse quadratic.By utilizing Hamilton’s principle,general formulae of the wave propagation were obtained to establish wave modes and phase velocity curves of the wave propagation in a functionally graded plate,including the effects of changing compositional characteristics of materials.

Matrix Riccati Equations in Optimal Control

In this paper, the matrix Riccati equation is considered. There is no general way for solving the matrix Riccati equation despite the many fields to which it applies. While scalar Riccati equation has been studied thoroughly, matrix Riccati equation of which scalar Riccati equations is a particular case, is much less investigated. This article proposes a change of variable that allows to find explicit solution of the Matrix Riccati equation. We then apply this solution to Optimal Control.

Convergence of Hyperbolic Neural Networks Under Riemannian Stochastic Gradient Descent

We prove,under mild conditions,the convergence of a Riemannian gradient descent method for a hyperbolic neural network regression model,both in batch gradient descent and stochastic gradient descent.We also discuss a Riemannian version of the Adam algorithm.We show numerical simulations of these algorithms on various benchmarks.

Analytical solutions fractional order partial differential equations arising in fluid dynamics

This article describes the solution procedure of the fractional Pade-Ⅱ equation and generalized Zakharov equation(GSEs)using the sine-cosine method.Pade-Ⅱ is an important nonlinear wave equation modeling unidirectional propagation of long-wave in dispersive media and GSEs are used to model the interaction between one-dimensional high,and low-frequency waves.Classes of trigonometric and hyperbolic function solutions in fractional calculus are discussed.Graphical simulations of the numerical solutions are flaunted by MATLAB.

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).