Search

This paper aims to reconstruct the initial condition of a hyperbolic equation with an unknown damping coefficient. Our approach involves approximating the hyperbolic equation's solution by its truncated Fourier expansion in the time domain and using the recently developed polynomial-exponential basis. This truncation process facilitates the elimination of the time variable, consequently, yielding a system of quasi-linear elliptic equations. To globally solve the system without needing an accurate initial guess, we employ the Carleman contraction principle. We provide several numerical examples to illustrate the efficacy of our method. The method not only delivers precise solutions but also showcases remarkable computational efficiency. [ABSTRACT FROM AUTHOR]

This document is an editorial for a special issue of the International Journal of Robust & Nonlinear Control on the topic of time-delay systems. It provides a historical overview of functional differential equations and highlights the contributions of Nikolay Nikolayevich Krasovskii and Aristide Halanay to the field. The special issue includes 21 research papers covering stability, robustness, observer theory, predictors, and delay-based approaches to partial differential equation models. Each paper is briefly described, covering topics such as stabilization of uncertain systems, sampled-data control for nonlinear systems, observer theory for differential delay systems, sensor fault estimation, and control approaches for vehicular platoons. The summaries aim to provide library patrons with a quick understanding of the papers' contents and their relevance to their research interests. [Extracted from the article]

The article starts from the challenge of the critical cases in difference operator stability for neutral functional differential equations (NFDE). Such cases occur in the NFDE associated to 1D$$ 1D $$ hyperbolic partial differential equations (PDE) dynamics in mechanical and hydraulic engineering. For some of such applications it resulted that the aforementioned critical (nonasymptotic) stability is connected to the character and level of the energy losses. It is shown that suitable choice of the losses to be taken into account can remove the critical stability properties and give the difference operator the asymptotic stability thus ensuring asymptotic stability for the system's dynamics and also other asymptotic properties. [ABSTRACT FROM AUTHOR]

For population systems modeled by age‐structured hyperbolic partial differential equations (PDEs), we redesign the existing feedback laws, designed under the assumption that the dilution input is directly actuated, to the more realistic case where dilution is governed by actuation dynamics (modeled simply by an integrator). In addition to the standard constraint that the population density must remain positive, the dilution dynamics introduce constraints of not only positivity of dilution, but possibly of given positive lower and upper bounds on dilution. We present several designs, of varying complexity, and with various measurement requirements, which not only ensure global asymptotic (and local exponential) stabilization of a desired positive population density profile from all positive initial conditions, but do so without violating the constraints on the dilution state. To develop the results, we exploit the relation between first‐order hyperbolic PDEs and an equivalent representation in which a scalar input‐driven mode is decoupled from input‐free infinite‐dimensional internal dynamics represented by an integral delay system. [ABSTRACT FROM AUTHOR]

We consider for a small parameter > 0 a parabolic convection–diffusion problem with Péclet number of order ( − 1) in a three-dimensional graph-like junction consisting of thin curvilinear cylinders with radii of order () connected through a domain (node) of diameter (). Inhomogeneous Robin type boundary conditions, that involve convective and diffusive fluxes, are prescribed both on the lateral surfaces of the thin cylinders and the boundary of the node. The asymptotic behavior of the solution is studied as → 0 , i.e. when the diffusion coefficients are eliminated and the thin junction is shrunk into a three-part graph connected in a single vertex. A rigorous procedure for the construction of the complete asymptotic expansion of the solution is developed and the corresponding energy and uniform pointwise estimates are proved. Depending on the directions of the limit convective fluxes, the corresponding limit problems (= 0) are derived in the form of first-order hyperbolic differential equations on the one-dimensional branches with novel gluing conditions at the vertex. These generalize the classical Kirchhoff transmission conditions and might require the solution of a three-dimensional cell-like problem associated with the vertex to account for the local geometric inhomogeneity of the node and the physical processes in the node. The asymptotic ansatz consists of three parts, namely, the regular part, node-layer part, and boundary-layer one. Their coefficients are classical solutions to mixed-dimensional limit problems. The existence and other properties of those solutions are analyzed. [ABSTRACT FROM AUTHOR]

We study nonlinear inverse problems for first order hyperbolic equations. We find a solution and an unknown coefficient depending on the time variable. For additional conditions we consider integral and boundary overdetermination conditions. We establish the existence of a regular solution, i.e., the solution possesses all generalized derivatives in the sense of Sobolev entering the equation. [ABSTRACT FROM AUTHOR]

Summary: This paper studies the adaptive synchronization of complex spatio‐temporal networks modeled by semi‐linear hyperbolic partial differential equations (CSTNSLHPDEs) as well as considering time‐invariant and time‐varying delays in a one‐dimensional space. Firstly, a distributed adaptive controller is proposed, where different nodes are with different adaptive gains. Secondly, four cases, CSTNSLHPDEs with time‐invariant delays and one single weight, with time‐invariant delays and multi‐weights, with time‐varying delays and one single weight, and with time‐varying delays and multi‐weights, are successively analyzed, and synchronization conditions of these four cases are obtained by using the proposed distributed adaptive controller. In the end, examples illustrate the effectiveness of the proposed distributed adaptive controller. [ABSTRACT FROM AUTHOR]

A body may have a structural, thermal, electromagnetic or optical role. In wave propagation, many models are described for transmission problems, whose solutions are defined in two or more domains. In this paper, we consider an inverse source hyperbolic problem on disconnected intervals, using solution point constraints. Applying a transform method, we reduce the inverse problems to direct ones, which are studied for well-posedness in special weighted Sobolev spaces. This means that the inverse problem is said to be well posed in the sense of Tikhonov (or conditionally well posed). The main aim of this study is to develop a finite difference method for solution of the transformed hyperbolic problems with a non-local differential operator and initial conditions. Numerical test examples are also analyzed. [ABSTRACT FROM AUTHOR]

This article studies the problem of prescribed‐time constrained tracking control for a class of 2×2$$ 2\times 2 $$ hyperbolic partial differential equation (PDE) systems with actuator dynamics, which are described by a set of nonlinear ordinary differential equations (ODEs). Since the control input only appears in the ODE subsystem rather than directly on the boundary of PDE subsystem, the control task becomes quite challenging. The most important is that for the control of such ODE controlled PDE systems we mainly make the following two contributions: (1) the controlled output of the PDE system tracks the reference signal within the prescribed time; (2) the controlled output and all the actuator states are constrained. It is the first time that such a prescribed‐time constrained tracking control problem is addressed for the PDE‐ODE coupled system considered in this article. Through rigorous theoretical proof, it is demonstrated that all the system states and control signals are bounded and sufficiently continuous by configuring appropriate design parameters. Finally, the performance is investigated via numerical simulation. [ABSTRACT FROM AUTHOR]

By using the power-exponential function method and the extended hyperbolic auxiliary equation method, we present the explicit solutions of the subsidiary elliptic-like equation. With the aid of the subsidiary elliptic-like equation and a simple transformation, we obtain the exact solutions of Hirota equation which include bright solitary wave, dark solitary wave, bell profile solitary wave solutions and Jacobian elliptic function solutions. Some of these solutions are found for the first time, which may be useful for depicting nonlinear physical phenomena. This approach can also be applied to solve the other nonlinear partial differential equations. [ABSTRACT FROM AUTHOR]

Local wellposedness for a nonlinear parabolic-hyperbolic coupled system modeling Micro-Electro-Mechanical System (MEMS) is studied. The particular device considered is a simple capacitor with two closely separated plates, one of which has motion modeled by a semilinear hyperbolic equation. The gap between the plates contains a gas and the gas pressure is taken to obey a quasilinear parabolic Reynolds' equation. Local-in-time existence of strict solutions of the system is shown, using well-known local-in-time existence results for the hyperbolic equation, then Hölder continuous dependence of its solution on that of the parabolic equation, and finally getting local-in-time existence for a combined abstract parabolic problem. Semigroup approaches are vital for the local-in-time existence results. [ABSTRACT FROM AUTHOR]

We consider the Cauchy problem for a one-dimensional hyperbolic equation whose lower-order coefficient and right-hand side oscillate in time with a high frequency and the amplitude of the lower-order coefficient is small. Under study is the reconstruction of the cofactors of these rapidly oscillating functions independent of the space variable from a partial asymptotics of a solution at some point of the space. The classical theory of inverse problems examines the numerous problems of determining unknown sources, and coefficients without rapid oscillations for various evolutionary equations, where the exact solution to the direct problem appears in the additional overdetermination condition. Equations with rapidly oscillating data are often encountered in modeling the physical, chemical, and other processes that occur in media subjected to high-frequency electromagnetic, acoustic, vibrational, and others fields, which demonstrates the topicality of perturbation theory problems on the reconstruction of unknown functions in high-frequency equations. We give some nonclassical algorithm for solving such problems that lies at the junction of asymptotic methods and inverse problems. In this case the overdetermination condition involves a partial asymptotics of solution of a certain length rather than the exact solution. [ABSTRACT FROM AUTHOR]

We consider a hyperbolic equation with variable leading part and nonlinearity in the lower-order term. The coefficients of the equation are smooth functions constant beyond some compact domain in the three-dimensional space. A plane wave with direction falls to the heterogeneity from the exterior of this domain. A solution to the corresponding Cauchy problem for the original equation is measured at boundary points of the domain for a time interval including the moment of arrival of the wave at these points. The unit vector is assumed to be a parameter of the problem and can run through all possible values sequentially. We study the inverse problem of determining the coefficient of the nonlinearity on using this information about solutions. We describe the structure of a solution to the direct problem and demonstrate that the inverse problem reduces to an integral geometry problem. The latter problem consists of constructing the desired function on using given integrals of the product of this function and a weight function. The integrals are taken along the geodesic lines of the Riemannian metric associated with the leading part of the differential equation. We analyze this new problem and find some stability estimate for its solution, which yields a stability estimate for solutions to the inverse problem. [ABSTRACT FROM AUTHOR]

Integral identities that are fulfilled for measure of difference between the exact solution of a parabolic-hyperbolic equation and any functions from a corresponding energy class are proved. These identities make it possible to derive two-sided a posteriori estimates for approximate solution to the corresponding Cauchy problem. The left-hand side of such an estimate is a natural measure of deviation from the solution, and the right-hand side depends on the problem data and the approximate solution itself and, therefore, it can be explicitly calculated. These estimates can be used to control the accuracy of approximate solutions and to compare solutions to Cauchy problems with different initial conditions. These estimates also allow one to quantitatively assess the effects occurring due to inaccuracies in the initial data and in the coefficients of the equation. [ABSTRACT FROM AUTHOR]

The influence of cleat permeability and porosity on the cumulative gas produced is studied. The cleat parameters are weighed for producing reserve in the least time and the minimum breakdown pressure requirement to create fractures. The reservoir is stimulated. The breakdown pressure required is more to create fractures at closer spacing than needed to increase the width of new fractures. The non-linear hyperbolic partial differential equations (PDE) are solved by the implicit finite difference method. Cleat permeability has a more significant influence on the gas produced than the porosity. Optimized values of 300 mD and 0.03 porosity sweep more area in a short time. Due to the faster propagation of drainage, the adsorbed gas is not recovered to nearly complete before the entire reservoir/pay zone length is drained. Fracturing a formation can be successful if the same adsorbed gas is spread over a larger area, and length is insufficient to accommodate the gas in place. The pay zone length needs to be larger than the present value. The novel approach of checking both the economic feasibility of fracturing a formation and its relation with the length of the reservoir, is discovered and solved here. [ABSTRACT FROM AUTHOR]

This paper is dedicated to investigating the chaos of a initial-boundary value (IBV) problem of a multi-dimensional weakly hyperbolic equation subject to two general nonlinear boundary conditions (NBCs). The existence and uniqueness of solution for the IBV problem are established. By employing the snap-back repeller and heteroclinic cycle theories, it has been proven that the IBV problem with a linear and a general NBCs exhibits chaos in the sense of both Devaney and Li–Yorke. Furthermore, these chaotic results are extended to the IBV problem with two general NBCs. Two stability criteria of the IBV problem are established, respectively, for the corresponding two cases of boundary conditions. Finally, numerical simulations are presented to illustrate the theoretical results. [ABSTRACT FROM AUTHOR]

The iterated Crank-Nicolson (ICN) method is a successful numerical algorithm for solving partial differential equations. When unequal weights are used in the ICN predictor-corrector process, the convergence rate is reduced to first order. In this paper, we propose two modified ICN algorithms that achieve second order convergence rate, based on two different ways of choosing the weights. The first approach employs geometrically averaged weights in two consecutive iterations, and the second one uses arithmetically averaged weights for two consecutive time steps. The stability and second order accuracy of our methods are verified using stability and truncation error analysis and numerically demonstrated by solving linear and semi-linear hyperbolic partial differential equations, Burgers' equation, and the Gross-Pitaevskii equation. [ABSTRACT FROM AUTHOR]

In this paper, we address the approximation of the coupling problem for the wave equation and Maxwell's equations of electromagnetism in the time domain in terms of electric field by means of a nodal linear finite element discretization in space, combined with a classical explicit finite difference scheme for time discretization. Our study applies to a particular case where the dielectric permittivity has a constant value outside a subdomain, whose closure does not intersect the boundary of the domain where the problem is defined. Inside this subdomain, Maxwell's equations hold. Outside this subdomain, the wave equation holds, which may correspond to Maxwell's equations with a constant permittivity under certain conditions. We consider as a model the case of first-order absorbing boundary conditions. First-order error estimates are proven in the sense of two norms involving first-order time and space derivatives under reasonable assumptions, among which lies a CFL condition for hyperbolic equations. The theoretical estimates are validated by numerical computations, which also show that the scheme is globally of the second order in the maximum norm in time and in the least-squares norm in space. [ABSTRACT FROM AUTHOR]

The potential of a finite-element technique based on an egalitarian meshing of the space–time domain D of physical problems described by parabolic or hyperbolic differential equations is explored. A least-squares minimization technique is applied in the meshed domain D to obtain stiffness-like matrices associated with various linear differential operators. Applications discussed include problems of boundary growth, and diffusive coalescence, in which D cannot be regarded as the Cartesian product of two independent domains in space and time. [ABSTRACT FROM AUTHOR]

Different from the construction process of the traditional finite difference method, we derive a large class of high-accuracy methods for approximating derivatives. Since Taylor expansion is avoided, the requirement for function smoothness in the new methods is greatly reduced. We analyze the approximation errors of the proposed methods and compare their approximation effects. As a typical application, we use the proposed methods to solve Elliptic, Hyperbolic and Parabolic problems numerically. For time-dependent problems, a class of new fully implicit difference schemes and the fourth-order Runge-Kutta method are both used for discretization. When using fully implicit difference schemes, unconditional stability and convergence are proved for hyperbolic and parabolic problems with the periodic boundary conditions. A large number of numerical experiments are provided to demonstrate the theoretical results. [ABSTRACT FROM AUTHOR]

This paper deals with the initial value problem for a semilinear wave equation on a bounded domain and solutions are required to vanish on the boundary of this domain. The essential feature of the situation considered here is that the ellipticity of the spatial part of the differential operator degenerates like the square of the distance from given point in the domain and so hyperbolicity is lost at this point. The assumptions ensure that u ≡ 0 is a stationary solution of the problem and the object is to study the stability of this solution with respect to perturbations of the initial data. Stability, instability and asymptotic stability are all considered. The assumptions about the nonlinear terms ensure that the problem has a well-defined linearization at u ≡ 0. There are simple cases where this linearization is asymptotically stable but u ≡ 0 is an unstable solution of the nonlinear problem. We also establish conditions under which the stability of a stationary solution u ≢ 0 can be determined using our results. The quadratic degeneracy at a point treated here is typical of what is required in models for acoustic (or sonic) black holes. It also occurs in a simplified Wheeler–DeWitt model which we discuss in some detail. [ABSTRACT FROM AUTHOR]

In this paper, we deal with the problem of stabilisation for a class of distributed semilinear systems. The considered systems present a state time delay and evolve in Banach spaces. A new stabilising feedback control is proposed. Two kinds of stabilizability, namely exponential and weak stabilizability, are investigated. An explicit estimate of the energy decay rate is given. Some illustrating applications to parabolic and hyperbolic equations with numerical simulations are considered. [ABSTRACT FROM AUTHOR]

The formulations and methods of solution of the problems of diagnostics of thermoelasticity on a stationary boundary of a solid body are considered. A weak formulation of the boundary-value problem is presented for a nonlinear hyperbolic equation (a nonstationary wave equation with density, elasticity modulus, and thermal stress depending on temperature) with mixed boundary conditions. An approximate solution of this problem in the space W 2 1 has been obtained, and an example of solution is given. Based on the weak formulation of the boundary-value problem, problems of diagnostics have been formulated in the sense of the theory of function traces. According to the theorem about the uniqueness of the relationship between the function and its trace in the space W 2 1 , a conclusion has been drawn about the need to control the convergence of the calculated wave function to the experimental one only on boundaries with unstable conditions in solving the problem of diagnostics. [ABSTRACT FROM AUTHOR]

We study hyperbolic equations with positive coefficients at potentials undergoing translations with respect to the spatial independent variable. The qualitative novelty of the investigation is that the real part of the symbol of the differential-difference operator contained in the equation is allowed to change its sign. Earlier, only the case where the said sign is constant was investigated. We find a condition relating the coefficient at the nonlocal term of the investigated equation and the length of the translation, guaranteeing the global solvability of the investigated equation. Under this condition, we explicitly construct a three-parametric family of smooth global solutions of the investigated equation. [ABSTRACT FROM AUTHOR]

For a linear hyperbolic equation of fourth order, a Dirichlet type boundary problem in an orthogonally convex domain is investigated. Sharp sufficient conditions guaranteeing solvability and well-posedness of the problem under consideration are established. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the global asymptotic stability of an age-structured population dynamics model with a Ricker's type of birth function. This model is a hyperbolic partial differential equation with a nonlinear and nonlocal boundary condition. We prove a uniform persistence result for the semiflow generated by this model. We obtain the existence of global attractors and we prove the global asymptotic stability of the positive equilibrium by using a suitable Lyapunov functional. Furthermore, we prove that our global asymptotic stability result is sharp, in the sense that Hopf bifurcation may occur as close as we want from the region global stability in the space of parameter. [ABSTRACT FROM AUTHOR]

A PDE‐constrained model predictive control (MPC) algorithm for open‐channel systems based on the Saint‐Vevant(S‐V) equations is investigated in this paper. The S‐V equations, which precisely model the dynamics of open‐channel systems, are quasi‐linear hyperbolic partial differential equations (PDEs) without analytical solutions. Directly applying the S‐V equations to an MPC controller design becomes sophisticated. In this work, the calculus of variation is used to obtain the adjoint equations and the adjoint analysis method is utilized to deduce the gradients of the MPC optimization problem. Particularly, the physical constraints involving both the state and control variables are also considered. A gradient‐based optimization algorithm in combination with the numerical computation of Preissmann implicit scheme is proposed to solve the constrained MPC optimization problem. The control performances of the developed PDE‐constrained MPC algorithm with respect to the controlled water levels and gate openings are compared with those of the MPC controller designed for the linearized model. All the simulation tests are carried out on an aqueduct reach in Yehe Irrigation District in Hebei Province, China. The results show that the proposed PDE‐constrained MPC algorithm is a promising method in dealing with the constraints in terms of hyperbolic PDEs, control variables and state variables simultaneously. [ABSTRACT FROM AUTHOR]

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach, and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. We demonstrate the method's capabilities and performance characteristics using a standard benchmark test network. • Formulation of a stochastic hyperbolic PDE as a high-dimensional parametric PDE. • Novel SFV method for quantification of intertemporal uncertainty in systems of conservation laws posed on graphs. • Formulation of the stochastic junction Riemann problem to compute numerical fluxes at the nodes of the graph. • Application of the SFV method to a real-world problem of stochastic gas flows on networks. • Convergence analysis for the statistical moments of the solution. [ABSTRACT FROM AUTHOR]

In this paper, we focus on the Riemann problem for two‐layered blood flow model, which is represented by a system of quasi‐linear hyperbolic partial differential equations (PDEs) derived from the Euler equations by vertical averaging across each layer. We consider the Riemann problem with varying velocities and equal constant density through arteries. For instance, the flow layer close to the wall of vessel has a slower average speed than the layer far from the vessel because of the viscous effect of the blood vessel. We first establish the existence and uniqueness of the corresponding Riemann solution by a thorough investigation of the properties of elementary waves, namely, shock wave, rarefaction wave, and contact discontinuity wave. Further, we extensively analyze the elementary wave interaction between rarefaction wave and shock wave with contact discontinuity and rarefaction wave and shock wave. The global structure of the Riemann solutions after each wave interaction is explicitly constructed and graphically illustrated towards the end. [ABSTRACT FROM AUTHOR]

Since self-gravity is crucial in the structure formation of the Universe, many hydrodynamics simulations with the effect of self-gravity have been conducted. The multigrid method is widely used as a solver for the Poisson equation of the self-gravity; however, the parallelization efficiency of the multigrid method becomes worse when we use a massively parallel computer, and it becomes inefficient with more than 104 cores, even for highly tuned codes. To perform large-scale parallel simulations (>104 cores), developing a new gravity solver with good parallelization efficiency is beneficial. In this article, we develop a new self-gravity solver using the telegraph equation with a damping coefficient, κ. Parallelization is much easier than the case of the elliptic Poisson equation since the telegraph equation is a hyperbolic partial differential equation. We analyse convergence tests of our telegraph equations solver and determine that the best non-dimensional damping coefficient of the telegraph equations is |$\tilde{\kappa } \simeq 2.5$|⁠. We also show that our method can maintain high parallelization efficiency even for massively parallel computations due to the hyperbolic nature of the telegraphic equation by weak-scaling tests. If the time-step of the calculation is determined by heating/cooling or chemical reactions, rather than the Courant–Friedrichs–Lewy (CFL) condition, our method may provide the method for calculating self-gravity faster than other previously known methods such as the fast Fourier transform and multigrid iteration solvers because gravitational phase velocity determined by the CFL condition using these time-scales is much larger than the fluid velocity plus sound speed. [ABSTRACT FROM AUTHOR]

The model of Timoshenko's shell theory of shells was used to analyze the dynamic characteristics of conical shells of variable thickness on a Pasternak elastic bed under nonstationary loading. Based on the Hamilton–Ostrogradsky variational principle, the equations of motion of a conical shell of variable thickness on a Pasternak elastic bed were derived. This system of hyperbolic differential equations is solved by the finite difference method. The numerical algorithm for solving the obtained equations is based on applying the integral-interpolation method for constructing difference schemes in the spatial coordinate and an explicit finite difference scheme for integration in the time coordinate. The influence of geometric dimensions, taper angle, and elastic media on the natural frequencies and other dynamic characteristics of a conical shell of variable thickness under the action of a pulsed load is analyzed using specific examples. New mechanical effects are revealed. [ABSTRACT FROM AUTHOR]

Dynamic perturbations reveal unconventional nonlinear behavior in rocks, as evidenced by field and laboratory studies. During the passage of seismic waves, rocks exhibit a decrease in elastic moduli, slowly recovering after. Yet, comprehensive physical models describing these moduli alterations remain sparse and insufficiently validated against observations. Here, we demonstrate the applicability of two physical damage models—the internal variable model (IVM) and the continuum damage model (CDM)—to provide quantitative descriptions of nonlinear co‐seismic elastic wave propagation observations. While the IVM uses one internal variable to describe the evolution of elastic material moduli, the CDM damage variable is a mathematical representation of microscopic defects. We recast the IVM and CDM models as nonlinear hyperbolic partial differential equations and implement 1D and 2D numerical simulations using an arbitrary high‐order discontinuous Galerkin method. We verify the modeling results with co‐propagating acousto‐elastic experimental measurements. Subsequently, we infer the parameters for these nonlinear models from laboratory experiments using probabilistic Bayesian inversion and 2D simulations. By adopting the Adaptive Metropolis Markov chain Monte Carlo method, we quantify the uncertainties of inferred parameters for both physical models, investigating their interplay in 70,000 simulations. We find that the damage variables can trade off with the stress‐strain nonlinearity in discernible ways. We discuss physical interpretations of both damage models and that our CDM quantitatively captures an observed damage increase with perturbation frequency. Our results contribute to a more holistic understanding of co‐seismic damage and post‐seismic recovery after earthquakes bridging the worlds of theoretical analysis and laboratory findings. Plain Language Summary: Rocks react to earthquakes by softening when seismic waves—the energy released by earthquakes—pass through them. Observations of such rock softening during the passage of seismic waves are common both in the laboratory and in the field. Interestingly, rocks gradually harden again once the shaking stops. Different physical mechanisms have been proposed to explain the observations. In this study, we put two existing theories to the test. One assumes that a term in the internal energy of the material increases with damage accumulation, while the other incorporates the opening and closing of micro‐cracks. We implement them into a powerful simulation program called ExaHyPE. This allows us to model how nonlinear waves move through rocks. When we compare the computer simulation outcomes with real laboratory tests, we find that both models match what we see in reality. Studying thousands of simulations with different model parameters, we find some intriguing insights. For instance, the initial state of strain and the tiny cracks that open and close within the rock may be key to understanding the hardening and softening process. We hope to use these physical models in future earthquake simulations, offering more accurate predictions of how our Earth's crust reacts to earthquakes. Key Points: We analyze two physical models suitable for simulations of nonlinear elastic wave propagation observed in the laboratoryThe experimentally observed co‐seismic acoustic modulus drop can be explained with nonlinear damage modelsWe use the Markov chain Monte Carlo method to explore connections and uncertainties of nonlinear parameters [ABSTRACT FROM AUTHOR]

This document is an editorial for a special issue of the journal Mathematics titled "Differential Equations of Mathematical Physics and Related Problems of Mechanics." The special issue covers a range of topics related to differential equations in mathematical physics and mechanics, including wave equations, spectral theory, scattering, and inverse problems. The editorial provides a summary of the published papers in the special issue, highlighting their contributions to the field. The document emphasizes the importance of the special issue in covering both applied and fundamental aspects of mathematics, physics, and their applications in various fields. The author expresses gratitude to the authors, reviewers, assistants, associate editors, and editors for their contributions to the special issue. The report does not provide specific details about the content of the papers or the nature of the special issue. [Extracted from the article]

A model of a two-phase flow of compressible immiscible fluids is presented. Its derivation is based on the use of the theory of symmetric hyperbolic thermodynamically compatible systems. The model is an extension of the previously proposed thermodynamically compatible model of compressible two-phase flows due to the inclusion of new state variables of a medium associated with surface-tension forces. The governing equations of the model form a hyperbolic system of differential equations of the first order and satisfy the laws of thermodynamics (energy conservation and entropy increase). The properties of the model equations are studied, and it is shown that the Young–Laplace law of capillary pressure is fulfilled in the asymptotic approximation at the continuum level. [ABSTRACT FROM AUTHOR]

In this article, we derive error estimates for two second-order numerical schemes for solving the viscous Cahn-Hilliard equation with hyperbolic relaxation, one based on the second-order Crank-Nicolson time marching method and the other on the backward differentiation formula. In both schemes, the nonlinear potential is discretized by the Invariant Energy Quadratization (IEQ) method, which employs an auxiliary variable to produce a linear and unconditional energy-stable structure. After assuming some reasonable boundedness and continuity conditions for the nonlinear function, the optimal error estimate is rigorously derived using mathematical induction. Finally, several numerical experiments are carried out to verify the theoretical predictions of the convergence rate and energy stability of the algorithms. [ABSTRACT FROM AUTHOR]

The paper considers nonlocal problems with integral conditions for hyperbolic equations with the Bessel differential operator whose statement substantially depends on the intervals where the parameter occurring in this operator varies. The well-posedness of these problems is studied according to a unified scheme based on the classical method of separation of variables, which is also used to study nonclassical problems with integral conditions for equations of elliptic–hyperbolic type containing the Bessel operator in one or two variables as well. [ABSTRACT FROM AUTHOR]

In this paper, we consider a coupled system of hyperbolic and biharmonic-wave equations with variable exponents in the damping and coupling terms. In each equation, the damping term is modulated by a time-dependent coefficient a(t) (or b(t)). First, we state and prove a well-posedness theorem of global weak solutions, by exploiting Galerkin's method and some compactness arguments. Then, using the multiplier method, we establish the decay rates of the solution energy, under suitable assumptions on the time-dependent coefficients and the range of the variable exponents. We end our work with some illustrative examples. [ABSTRACT FROM AUTHOR]

A model of a two-phase flow of compressible immiscible fluids is presented. Its derivation is based on the use of the theory of symmetric hyperbolic thermodynamically compatible systems. The model is an extension of the previously proposed thermodynamically compatible model of compressible two-phase flows due to the inclusion of new state variables of a medium associated with surface-tension forces. The governing equations of the model form a hyperbolic system of differential equations of the first order and satisfy the laws of thermodynamics (energy conservation and entropy increase). The properties of the model equations are studied, and it is shown that the Young–Laplace law of capillary pressure is fulfilled in the asymptotic approximation at the continuum level. [ABSTRACT FROM AUTHOR]

This research article addresses a nonclassical initial boundary value problem characterized by a non-local constraint within the framework of a pseudo-hyperbolic equation. Employing rigorous analytical techniques, the paper establishes the existence, uniqueness, and continuous dependence of a strong solution to the problem at hand. With respect to the associated linear problem, the uniqueness of its solution is ascertained through an energy inequality, which provides an a priori bound for the solution. Moreover, the solvability of this linear problem is verified by proving that the operator range engendered by the problem is indeed dense. Extending the analysis to the nonlinear problem, an iterative methodology is utilized. This approach is predicated on the insights gained from the linear problem and facilitates the demonstration of both the existence and uniqueness of a solution for the nonlinear problem under study. Consequently, the paper contributes a robust mathematical framework for solving both linear and nonlinear variants of complex initial boundary value problems with non-local constraints. [ABSTRACT FROM AUTHOR]

In this article, we discuss an inverse problem of determining unknown coefficients of the first-order derivatives in an ultrahyperbolic equation. By a finite set of measurements, we prove the uniqueness of solution of the problem in semi-geodesic coordinates under some conditions on the principal coefficients of the equation. Our main tool is a Carleman estimate. [ABSTRACT FROM AUTHOR]

We deal with an optimal design problem governed by an initial‐boundary value problem for a hyperbolic variational inequality describing the perpendicular vibrations of a simply supported anisotropic elastic plate against a rigid obstacle. A variable thickness of a plate plays the role of a control variable. The admissible set of states for the design problem contains solutions of a state problem gained as limits of the sequences of functions solving penalized problems. We verify the existence of the optimal thickness function using the optimal controls achieved as solutions of optimal design problems connected with penalized state problems. [ABSTRACT FROM AUTHOR]

The particle size distribution (PSD) in emulsion polymerization (EP) has been modeled in the past using either the pseudo bulk (PB) or the 0-1/0-1-2 approaches. There is some controversy on the proper type of model to be used to simulate the experimental PSDs, which are apparently broader than the theoretical ones. Additionally, the numerical technique employed to solve the model equations, involving hyperbolic partial differential equations (PDEs) with moving and possibly steep fronts, has to be precise and robust, which is not a trivial matter. A deterministic kinetic model for the PSD evolution of ab initio EP of vinyl monomers was developed to investigate these issues. The model considers three phases, micellar nucleation, and particles that can contain n ≥ 0 radicals. Finite volume (FV) and weighted-residual methods are used to solve the system of PDEs and compared; their limitations are also identified. The model was validated by comparing predictions with data of monomer conversion and PSD for the batch emulsion homopolymerization of styrene (Sty) and methyl methacrylate (MMA) using sodium dodecyl sulfate (SDS)/potassium persulfate (KPS) at 60 °C, as well as the copolymerization of Sty-MMA (50/50; mol/mol) at 50 and 60 °C. It is concluded that the PB model has a structural problem when attempting to adequately represent PSDs with steep fronts, so its use is discouraged. On the other hand, there is no generalized evidence of the need to add a stochastic term to enhance the PSD prediction of EP deterministic models. [ABSTRACT FROM AUTHOR]

We prove global Lipschitz stability for an inverse source problem of a system of wave equations on a Lorentzian manifold. The method used in this paper is widely known as the Bukhgeim–Klibanov method. However, the conventional method is not sufficient for the application to the hyperbolic partial differential equation with time-dependent coefficients to obtain the Lipschitz stability, and further innovations are needed. In this paper, we present an improved global Carleman estimate and an energy estimate to obtain the Lipschitz stability. [ABSTRACT FROM AUTHOR]

We propose a novel framework for model-order reduction of hyperbolic differential equations. The approach combines a relaxation formulation of the hyperbolic equations with a discretization using shifted base functions. Model-order reduction techniques are then applied to the resulting system of coupled ordinary differential equations. On computational examples including in particular the case of shock waves we show the validity of the approach and the performance of the reduced system. [ABSTRACT FROM AUTHOR]

The existence and uniqueness theorems of the solution to the boundary-value problem for a parabolic-hyperbolic fractional-order equation with the gluing condition are proved. [ABSTRACT FROM AUTHOR]