scispace - formally typeset
Search or ask a question

Showing papers on "Hyperbolic partial differential equation published in 2020"



Journal ArticleDOI
TL;DR: A new family of very high order accurate direct Arbitrary-Lagrangian-Eulerian Finite Volume and Discontinuous Galerkin schemes for the solution of nonlinear hyperbolic PDE systems on moving two-dimensional Voronoi meshes that are regenerated at each time step and which explicitly allow topology changes in time.

64 citations


Book ChapterDOI
10 Mar 2020
TL;DR: In this paper, the authors describe some results concerning oscillatory integrals and their application to Radon-like transforms and show that Radon transforms have an impact on the estimates for averaging operators involving integration over lower dimensional manifolds and connections with hyperbolic equations.
Abstract: This chapter describes some results concerning oscillatory integrals and, in particular, their application to Radon-like transforms. It briefly discusses three classes of oscillatory integrals. The first class consists of maximal oscillatory integrals. The second is made up of the oscillatory integrals arising in restrictions theorems and Bochner–Riesz summability. The third class contains the oscillatory integrals related to the Radon-like transforms, which is most closely related to Fourier integral operators. The chapter shows that the study of Radon-like transforms has an impact on the estimates for averaging operators involving integration over lower dimensional manifolds and connections with hyperbolic equations, relations of Radon transforms with Fourier integral operators, the initial study of singular Radon transforms as “Hilbert transforms along curves,” and the relevance to several complex variables.

38 citations


Journal ArticleDOI
TL;DR: A novel thermoelasticity model based on the Moore–Gibson–Thompson equation is considered in cases where some of these models fail to be positive, and is applied to investigate the wave propagation in an isotropic and infinite body subjected to a continuous thermal line source.
Abstract: Many attempts have been made to investigate the classical heat transfer of Fourier, and a number of improvements have been implemented. In this work, we consider a novel thermoelasticity model based on the Moore-Gibson-Thompson equation in cases where some of these models fail to be positive. This thermomechanical model has been constructed in combination with a hyperbolic partial differential equation for the variation of the displacement field and a parabolic differential equation for the temperature increment. The presented model is applied to investigate the wave propagation in an isotropic and infinite body subjected to a continuous thermal line source. To solve this problem, together with Laplace and Hankel transform methods, the potential function approach has been used. Laplace and Hankel inverse transformations are used to find solutions to different physical fields in the space-time domain. The problem is validated by calculating the numerical calculations of the physical fields for a given material. The numerical and theoretical results of other thermoelastic models have been compared with those described previously.

37 citations


Journal ArticleDOI
TL;DR: This paper presents a control design for the one-phase Stefan problem under actuator delay via a backstepping method, and robustness analysis with respect to the delay mismatch between the plant and the controller is studied, which provides analogous results to the exact compensation by restricting the control gain.
Abstract: This paper presents a control design for the one-phase Stefan problem under actuator delay via a backstepping method. The Stefan problem represents a liquid-solid phase change phenomenon which describes the time evolution of a material's temperature profile and the interface position. The actuator delay is modeled by a first-order hyperbolic partial differential equation (PDE), resulting in a cascaded transport-diffusion PDE system defined on a time-varying spatial domain described by an ordinary differential equation (ODE). Two nonlinear backstepping transformations are utilized for the control design. The setpoint restriction is given to guarantee a physical constraint on the proposed controller for the melting process. This constraint ensures the exponential convergence of the moving interface to a setpoint and the exponential stability of the temperature equilibrium profile and the delayed controller in the  1 norm. Furthermore, robustness analysis with respect to the delay mismatch between the plant and the controller is studied, which provides analogous results to the exact compensation by restricting the control gain.

35 citations


Journal ArticleDOI
TL;DR: Several numerical examples, including the long-time evolution of a stable neutron star in anti-Cowling approximation, are presented in order to show the obtained improvements with respect to the standard formulation without special treatment of the curl involution constraints.

33 citations


Journal ArticleDOI
TL;DR: The (2 + 1)-dimensional Schrodinger complex equations are essential physical models that describe the short-range pulse spread in nonlinear media fiber optics as discussed by the authors, and two different techniques are used to construct novel complex solutions to Schroffinger's complex hyperbolic model.
Abstract: The (2 + 1)-dimensional Schrodinger complex equations are essential physical models that describe the short-range pulse spread in nonlinear media fiber optics. We construct novel complex solutions to Schrodinger's complex hyperbolic model by using two different techniques. One method is characterized by the efficient algebraic equations that eventually form. Meanwhile, it uses the dependency variable expressions and its derivatives in the differential equation of the polynomial of a solitary wave. New acquired solutions are rational and exponential solutions expressed by periodic solutions. The solutions are illustrated through 3D- and 2D- plots to clarify the physical features for this model.

30 citations


Journal ArticleDOI
TL;DR: In this article, an analogue of the Gellerstedt problem for a loaded parabolic-hyperbolic equation of the third order in an infinite three-dimensional domain is studied, and sufficient conditions are found such that all differentiation operations are legal.
Abstract: In this paper, we study an analogue of the Gellerstedt problem for a loaded parabolic-hyperbolic equation of the third order in an infinite three-dimensional domain. The main method to study this Gellerstedt problem is the Fourier transform. Based on the Fourier transform, we reduce the considering problem to a planar analogue of the Gellerstedt spectral problem with a spectral parameter. The uniqueness of the solution of this problem is proved by the new extreme principle for loaded third-order equations of the mixed type. The existence of a regular solution of the Gellerstedt spectral problem is proved by the method of integral equations. In addition, the asymptotic behavior of the solution of the Gellerstedt spectral problem is studied for large values of the spectral parameter. Sufficient conditions are found such that all differentiation operations are legal in this work.

28 citations


Journal ArticleDOI
TL;DR: In this paper, the authors considered the energy-critical defocusing nonlinear wave equation and established almost sure global existence and scattering for randomized radially symmetric initial data in the random data setting.
Abstract: We consider the energy-critical defocusing nonlinear wave equation on $\mathbb{R}^4$ and establish almost sure global existence and scattering for randomized radially symmetric initial data in $H^s_x(\mathbb{R}^4) \times H^{s-1}_x(\mathbb{R}^4)$ for $\frac{1}{2} < s < 1$. This is the first almost sure scattering result for an energy-critical dispersive or hyperbolic equation with scaling super-critical initial data. The proof is based on the introduction of an approximate Morawetz estimate to the random data setting and new large deviation estimates for the free wave evolution of randomized radially symmetric data.

26 citations


Journal ArticleDOI
TL;DR: GRAMSES as mentioned in this paper is a pipeline for nonlinear cosmological N-body simulations in General Relativity (GR), with constant mean curvature and minimum distortion gauge fixings, which provides a fully nonlinear and background independent framework for relativistic cosmology.
Abstract: We present GRAMSES, a new pipeline for nonlinear cosmological N-body simulations in General Relativity (GR). This code adopts the Arnowitt-Deser-Misner (ADM) formalism of GR, with constant mean curvature and minimum distortion gauge fixings, which provides a fully nonlinear and background independent framework for relativistic cosmology. Employing a fully constrained formulation, the Einstein equations are reduced to a set of ten elliptical equations which are solved using multigrid relaxation with adaptive mesh refinements (AMR), and three hyperbolic equations for the evolution of tensor degrees of freedom. The current version of GRAMSES neglects the latter by using the conformal flatness approximation, which allows it to compute the two scalar and two vector degrees of freedom of the metric. In this paper we describe the methodology, implementation, code tests and first results for cosmological simulations in a ΛCDM universe, while the generation of initial conditions and physical results will be discussed elsewhere. Inheriting the efficient AMR and massive parallelisation infrastructure from the publicly-available N-body and hydrodynamic simulation code RAMSES, GRAMSES is ideal for studying the detailed behaviour of spacetime inside virialised cosmic structures and hence accurately quantifying the impact of backreaction effects on the cosmic expansion, as well as for investigating GR effects on cosmological observables using cosmic-volume simulations.

25 citations


Journal ArticleDOI
TL;DR: A variational framework for the distance-minimizing data-driven computing method is proposed that renders the high-order polynomial interpolation to straightforward implementation and spectral element methods are employed to reduce the computational cost while assuring the high accuracy of the data- driven solution.

Journal ArticleDOI
TL;DR: In this article, a new limiting method for bicompact schemes is proposed that preserves them conservative, based upon a finite-element treatment of the bic-compact approximation.

Journal ArticleDOI
TL;DR: In this article, a modified argument for proving conditional stability for inverse problems of determining spatially varying functions in evolution equations by Carleman estimates is provided, which can simplify the existing proofs.
Abstract: In this article, we provide a modified argument for proving conditional stability for inverse problems of determining spatially varying functions in evolution equations by Carleman estimates. Our method needs not any cut-off procedures and can simplify the existing proofs. We establish the conditional stability for inverse source problems for a hyperbolic equation and a parabolic equation, and our method is widely applicable to various evolution equations.

Journal ArticleDOI
TL;DR: A novel first-order reformulation of the most well-known Boussinesq-type systems that are used in ocean engineering, based on a relaxed augmented system in which the divergence constraints of the velocity flow variables are coupled with the other conservation laws via an evolution equation for the depth-averaged non-hydrostatic pressures.
Abstract: In this paper, we propose a novel first-order reformulation of the most well-known Boussinesq-type systems that are used in ocean engineering. This has the advantage of collecting in a general framework many of the well-known systems used for dispersive flows. Moreover, it avoids the use of high-order derivatives which are not easy to treat numerically, due to the large stencil usually needed. These first-order PDE dispersive systems are then approximated by a novel set of first-order hyperbolic equations. Our new hyperbolic approximation is based on a relaxed augmented system in which the divergence constraints of the velocity flow variables are coupled with the other conservation laws via an evolution equation for the depth-averaged non-hydrostatic pressures. The most important advantage of this new hyperbolic formulation is that it can be easily discretized with explicit and high-order accurate numerical schemes for hyperbolic conservation laws. There is no longer need of solving implicitly some linear system as it is usually done in many classical approaches of Boussinesq-type models. Here a third-order finite volume scheme based on a CWENO reconstruction has been used. The scheme is well-balanced and can treat correctly wet–dry areas and emerging topographies. Several numerical tests, which include idealized academic benchmarks and laboratory experiments are proposed, showing the advantage, efficiency and accuracy of the technique proposed here.

Journal ArticleDOI
TL;DR: B bound estimates for the two fastest wave speeds emerging from the solution of the Riemann problem for three well-known hyperbolic systems, namely the Euler equations of gas dynamics, the shallow water equations and the blood flow equations for arteries are provided.

Journal ArticleDOI
TL;DR: In this article, a two-phase model and its application to wavefields numerical simulation are discussed in the context of modeling of compressible fluid flows in elastic porous media, and the results of solving some numerical test problems are presented and discussed.

Journal ArticleDOI
TL;DR: In this paper, a linear one-dimensional Bresse system consisting of three hyperbolic equations coupled in a certain manner under mixed homogeneous Dirichlet-Neumann boundary conditions is considered.
Abstract: In this paper, we consider a linear one-dimensional Bresse system consisting of three hyperbolic equations coupled in a certain manner under mixed homogeneous Dirichlet-Neumann boundary conditions. Here, we consider that only the longitudinal displacement is damped, and the vertical displacement and shear angle displacement are free. We prove the well-posedness of the system and some exponential, lack of exponential and polynomial stability results depending on the coefficients of the equations and the smoothness of initial data. At the end, we use some numerical approximations based on finite difference techniques to validate the theoretical results. The proof is based on the semigroup theory and a combination of the energy method and the frequency domain approach.

Journal ArticleDOI
TL;DR: In this article, the authors propose a solution to the problem of nonequilibrium thermodynamics in the framework of SHTC equations, which does not rely on postulates of equilibrium irreversible thermodynamics but treats irreversible processes from the non-equilibrium point of view.
Abstract: The lack of formulation of macroscopic equations for irreversible dynamics of viscous heat-conducting media compatible with the causality principle of Einstein's special relativity and the Euler-Lagrange structure of general relativity is a long-lasting problem. In this paper, we propose a possible solution to this problem in the framework of SHTC equations. The approach does not rely on postulates of equilibrium irreversible thermodynamics but treats irreversible processes from the non-equilibrium point of view. Thus, each transfer process is characterized by a characteristic velocity of perturbation propagation in the non-equilibrium state, as well as by an intrinsic time/length scale of the dissipative dynamics. The resulting system of governing equations is formulated as a first-order system of hyperbolic equations with relaxation-type irreversible terms. Via a formal asymptotic analysis, we demonstrate that classical transport coefficients such as viscosity, heat conductivity, etc., are recovered in leading terms of our theory as effective transport coefficients. Some numerical examples are presented in order to demonstrate the viability of the approach. This article is part of the theme issue 'Fundamental aspects of nonequilibrium thermodynamics'.

Book ChapterDOI
01 Jan 2020
TL;DR: It is argued that the reduced system, constructed with the new method, can be identified by a reduced energy that mimics the energy of the high-fidelity system, and the loss in energy, associated with the model reduction, remains constant in time.
Abstract: In the past decade, model order reduction (MOR) has been successful in reducing the computational complexity of elliptic and parabolic systems of partial differential equations (PDEs). However, MOR of hyperbolic equations remains a challenge. Symmetries and conservation laws, which are a distinctive feature of such systems, are often destroyed by conventional MOR techniques which result in a perturbed, and often unstable reduced system. The importance of conservation of energy is well-known for a correct numerical integration of fluid flow. In this paper, we discuss model reduction, that exploits skew-symmetry of conservative and centered discretization schemes, to recover conservation of energy at the level of the reduced system. Moreover, we argue that the reduced system, constructed with the new method, can be identified by a reduced energy that mimics the energy of the high-fidelity system. Therefore, the loss in energy, associated with the model reduction, remains constant in time. This results in an, overall, correct evolution of the fluid that ensures robustness of the reduced system. We evaluate the performance of the proposed method through numerical simulation of various fluid flows, and through a numerical simulation of a continuous variable resonance combustor model.

Journal ArticleDOI
TL;DR: It is proved that u ( x, t ) = 0 is asymptotic stable in terms of natural energy associated with the solution of the above equation.
Abstract: In this paper, we consider the following quasilinear hyperbolic equation involving variable sources: u t t − d i v ( | ∇ u | p ( x ) − 2 ∇ u ) + | u t | m ( x ) − 2 u t = | u | q ( x ) − 2 u . Some energy estimates and K o m o r n i k inequality are used to prove some uniform estimates of decay rates of the solution. And then, we prove that u ( x , t ) = 0 is asymptotic stable in terms of natural energy associated with the solution of the above equation. As we know, such results are seldom seen for the variable exponent case. At last, we give some numerical examples to illustrate our results.

Journal ArticleDOI
TL;DR: In this article, the authors considered the semilinear wave equation with the multiplication of logarithmic and polynomial nonlinearities and established the global existence and finite time blow up of solutions at three different energy levels using potential well method.
Abstract: In this paper we consider the semilinear wave equation with the multiplication of logarithmic and polynomial nonlinearities. We establish the global existence and finite time blow up of solutions at three different energy levels (E(0) < d, E(0) = d and E(0) > 0) using potential well method. The results in this article shed some light on using potential wells to classify the solutions of the semilinear wave equation with the product of polynomial and logarithmic nonlinearity.

Book ChapterDOI
TL;DR: A new numerical method for solving time dependent Maxwell equations, which is also suitable for general linear hyperbolic equations, based on an unstructured partitioning of the spacetime domain into tent-shaped regions that respect causality is presented.
Abstract: We present a new numerical method for solving time dependent Maxwell equations, which is also suitable for general linear hyperbolic equations. It is based on an unstructured partitioning of the spacetime domain into tent-shaped regions that respect causality. Provided that an approximate solution is available at the tent bottom, the equation can be locally evolved up to the top of the tent. By mapping tents to a domain which is a tensor product of a spatial domain with a time interval, it is possible to construct a fully explicit scheme that advances the solution through unstructured meshes. This work highlights a difficulty that arises when standard explicit Runge Kutta schemes are used in this context and proposes an alternative structure-aware Taylor time-stepping technique. Thus explicit methods are constructed that allow variable time steps and local refinements without compromising high order accuracy in space and time. These Mapped Tent Pitching (MTP) schemes lead to highly parallel algorithms, which utilize modern computer architectures extremely well.

Journal ArticleDOI
TL;DR: A version of the convexification numerical method for a Coefficient Inverse Problem for a 1D hyperbolic PDE is presented and the global convergence of the gradient projection method is established.
Abstract: A version of the convexification numerical method for a Coefficient Inverse Problem for a 1D hyperbolic PDE is presented. The data for this problem are generated by a single measurement event. This method converges globally. The most important element of the construction is the presence of the Carleman Weight Function in a weighted Tikhonov-like functional. This functional is strictly convex on a certain bounded set in a Hilbert space, and the diameter of this set is an arbitrary positive number. The global convergence of the gradient projection method is established. Computational results demonstrate a good performance of the numerical method for noisy data.

Posted Content
TL;DR: This work demonstrates how to make DNNs satisfy boundary and initial conditions automatically in a systematic manner and ensures that the loss function in MIM is free of the penalty term and does not have any modeling error.
Abstract: In theory, boundary and initial conditions are important for the wellposedness of partial differential equations (PDEs). Numerically, these conditions can be enforced exactly in classical numerical methods, such as finite difference method and finite element method. Recent years have witnessed growing interests in solving PDEs by deep neural networks (DNNs), especially in the high-dimensional case. However, in the generic situation, a careful literature review shows that boundary conditions cannot be enforced exactly for DNNs, which inevitably leads to a modeling error. In this work, based on the recently developed deep mixed residual method (MIM), we demonstrate how to make DNNs satisfy boundary and initial conditions automatically in a systematic manner. As a consequence, the loss function in MIM is free of the penalty term and does not have any modeling error. Using numerous examples, including Dirichlet, Neumann, mixed, Robin, and periodic boundary conditions for elliptic equations, and initial conditions for parabolic and hyperbolic equations, we show that enforcing exact boundary and initial conditions not only provides a better approximate solution but also facilitates the training process.

Posted Content
TL;DR: In this article, a two time-scale scheme was introduced which allows to extend the method of minimizing movements to hyperbolic problems, which is used to show the existence of weak solutions to a fluid-structure interaction problem between a nonlinear, visco-elastic, $n$-dimensional bulk solid and an incompressible fluid governed by the Navier-Stokes equations.
Abstract: We introduce a two time-scale scheme which allows to extend the method of minimizing movements to hyperbolic problems. This method is used to show the existence of weak solutions to a fluid-structure interaction problem between a nonlinear, visco-elastic, $n$-dimensional bulk solid governed by a hyperbolic evolution and an incompressible fluid governed by the ($n$-dimensional) Navier-Stokes equations for $n\geq 2$.

Posted Content
TL;DR: A numerical method is proposed to solve an inverse source problem of computing the initial condition of hyperbolic equations from the measurements of Cauchy data by the quasi-reversibility method, and it is rigorously proved that the convergence of this method as the noise level tends to 0.
Abstract: We propose a numerical method to solve an inverse source problem of computing the initial condition of hyperbolic equations from the measurements of Cauchy data. This problem arises in thermo- and photo- acoustic tomography in a bounded cavity, in which the reflection of the wave makes the widely-used approaches, such as the time reversal method, not applicable. In order to solve this inverse source problem, we approximate the solution to the hyperbolic equation by its Fourier series with respect to a special orthogonal basis of $L^2$. Then, we derive a coupled system of elliptic equations for the corresponding Fourier coefficients. We solve it by the quasi-reversibility method. The desired initial condition follows. We rigorously prove the convergence of the quasi-reversibility method as the noise level tends to 0. Some numerical examples are provided. In addition, we numerically prove that the use of the special basic above is significant.

Journal ArticleDOI
TL;DR: This paper considers the Runge--Kutta discontinuous Galerkin (RKDG) method to solve linear constant-coefficient hyperbolic equations, where the fourth-order explicit Runge -- Kutta time-marching method is considered.
Abstract: In this paper we consider the Runge--Kutta discontinuous Galerkin (RKDG) method to solve linear constant-coefficient hyperbolic equations, where the fourth-order explicit Runge--Kutta time-marching...

Journal ArticleDOI
TL;DR: In this paper, the integrability of nonlinear lattices is studied in the sense of reductions of a chain to a system of hyperbolic equations of an arbitrarily high order that are integrable in the Darboux sense.
Abstract: We study the problem of the integrable classification of nonlinear lattices depending on one discrete and two continuous variables. By integrability, we mean the presence of reductions of a chain to a system of hyperbolic equations of an arbitrarily high order that are integrable in the Darboux sense. Darboux integrability admits a remarkable algebraic interpretation: the Lie—Rinehart algebras related to both characteristic directions corresponding to the reduced system of hyperbolic equations must have a finite dimension. We discuss a classification algorithm based on the properties of the characteristic algebra and present some classification results. We find new examples of integrable equations.

Journal ArticleDOI
TL;DR: The result demonstrates that the superconvergence property of the semi-discrete DG method is preserved, and the optimal order in time is provided under the smoothness assumption that is independent of the number of stages.
Abstract: In this paper, we shall establish the superconvergence property of the Runge–Kutta discontinuous Galerkin (RKDG) method for solving a linear constant-coefficient hyperbolic equation. The RKDG method is made of the discontinuous Galerkin (DG) scheme with upwind-biased numerical fluxes coupled with the explicit Runge–Kutta algorithm of arbitrary orders and stages. Superconvergence results for the numerical flux, cell averages as well as the solution and derivative at some special points are shown, which are based on a systematical study of the $$\hbox {L}^2$$ -norm stability for the RKDG method and the incomplete correction techniques for the well-defined reference functions at each time stage. The result demonstrates that the superconvergence property of the semi-discrete DG method is preserved, and the optimal order in time is provided under the smoothness assumption that is independent of the number of stages. As a byproduct of the above superconvergence study, the expected order of the post-processed solution is obtained when a special initial solution is used. Some numerical experiments are also given.

Journal ArticleDOI
TL;DR: A numerical scheme, based on power series approximations of the kernels is formulated, taking into account the fact that the kernels are only piecewise differentiable, and is applicable to the dual problem of boundary stabilization of coupled linear reaction–diffusion PDEs.