Scalar Equation Research Articles

On stability of equations with an infinite distributed delay

For scalar equations of population dynamics with an infinite distributed delay x′(t)=r(t)[∫−∞tf(x(s))dsR(t,s)−x(t)], x(t)=φ(t), t⩽t0, where f is the delayed production function, we consider asymptotic stability of the zero and a positive equilibrium K. It is assumed that the initial distribution is an arbitrary continuous function. Introducing conditions on the memory decay, we characterize functions f such that any solution with nonnegative nontrivial initial conditions tends to a positive equilibrium. The differences between finite and infinite delays are outlined, in particular, we present an example when the weak Allee effect (meaning that f′(0)=1 together with f(x)>x , x∈(0,K) ) which has no effect in the finite delay case (all solutions are persistent) can lead to extinction in the case of an infinite delay.

Wild solutions to scalar Euler-Lagrange equations

We study very weak solutions to scalar Euler-Lagrange equations associated with quadratic convex functionals. We investigate whether W 1 , 1 W^{1,1} solutions are necessarily W loc 1 , 2 W^{1,2}_{\operatorname {loc}} , which would make the theories by De Giorgi-Nash and Schauder applicable. We answer this question positively for a suitable class of functionals. This is an extension of Weyl’s classical lemma for the Laplace equation to a wider class of equations under stronger regularity assumptions. Conversely, using convex integration, we show that outside this class of functionals, there exist W 1 , 1 W^{1,1} solutions of locally infinite energy to scalar Euler-Lagrange equations. In addition, we prove an intermediate result which permits the regularity of a W 1 , 1 W^{1,1} solution to be improved to W loc 1 , 2 W^{1,2}_{\operatorname {loc}} under suitable assumptions on the functional and solution.

Periodic Solutions in a Simple Delay Differential Equation

A simple-form scalar differential equation with delay and nonlinear negative periodic feedback is considered. The existence of several types of slowly oscillating periodic solutions is shown with the same and double periods of the feedback coefficient. The periodic solutions are built explicitly in the case with piecewise constant nonlinearities involved. The periodic dynamics are shown to persist under small perturbations of the equation, which make it smooth. The theoretical results are verified through extensive numerical simulations.

The classical double copy in curved spacetimes: perturbative Yang-Mills from the bi-adjoint scalar

We formulate a version of the double copy for classical fields in curved spacetimes. We provide a correspondence between perturbative solutions to the bi-adjoint scalar equations and those of the Yang-Mills equations in position space. At the linear level, we show that there exists a map between these solutions for maximally symmetric spacetime backgrounds, that provides every Yang-Mills solution by the action of an appropriate differential operator on a bi-adjoint scalar solution. Given the existence of a linearized map, we show that it is possible to cast the solutions of the Yang-Mills equations at arbitrary perturbation order in terms of the corresponding bi-adjoint scalar solutions. This all-order map is reminiscent of the flat space BCJ double copy, and works for any curved spacetime where the perturbative expansion holds. We show that these results have the right flat space limit, and that the correspondence is agnostic to the choice of gauge.

Finslerian extension of an anisotropic strange star in the domain of modified gravity

In this article, we apply the Finsler spacetime to develop the Einstein field equations in the extension of modified geometry. Following Finsler geometry, which is focused on the tangent bundle with a scalar function, a scalar equation should be the field equation that defines this structure. This spacetime maintains the required causality properties on the generalized Lorentzian metric manifold. The matter field is coupled with the Finsler geometry to produce the complete action. The developed Einstein field equations are employed on the strange stellar system to improve the study. The interior of the system is composed of a strange quark matter, maintained by the MIT bag equation of state. In addition, the modified Tolman–Oppenheimer–Volkov (TOV) equation is formulated. In particular, the anisotropic stress attains the maximum at the surface. The mass-central density variation confirms the stability of the system.

A finite element framework for fluid–structure interaction of turbulent cavitating flows with flexible structures

We present a finite element framework for the numerical prediction of cavitating turbulent flows interacting with flexible structures. The vapor-fluid phases are captured through a homogeneous mixture model, with a scalar transport equation governing the spatio-temporal evolution of cavitation dynamics. High-density gradients in the two-phase cavitating flow motivate the use of a positivity-preserving Petrov–Galerkin stabilization method in the variational framework. A mass transfer source term introduces local compressibility effects arising as a consequence of phase change. The turbulent fluid flow is modeled through a dynamic subgrid-scale method for large eddy simulations. The flexible structure is represented by a set of eigenmodes, obtained through a modal decomposition of the linear elasticity equations. While a partitioned iterative approach is adopted to couple the structural dynamics and cavitating fluid flow, the deforming flow domain is described by an arbitrary Lagrangian-Eulerian frame of reference. We establish the fidelity of the proposed framework by comparing it against experimental and numerical studies for both rigid and flexible hydrofoils in cavitating flows. Under unstable partial cavitating conditions, we identify specific vortical structures leading to cloud cavity collapse. We further explore features of cavitating flow past a rigid body such as re-entrant jet and turbulence-cavity interactions during cloud cavity collapse. Based on the validation study conducted over a flexible NACA66 rectangular hydrofoil, we elucidate the role of cavity and vortex shedding in governing the structural dynamics. Subsequently, we identify a broad spectrum frequency band whose central peak does not correlate to the frequency content of the cavitation dynamics or the natural frequencies of the structure, indicating the induction of unsteady flow patterns around the hydrofoil. Finally, we discuss the coupled fluid–structure dynamics during a cavitation cycle and the underlying mechanism associated with the promotion and mitigation of cavitation.

Port-Hamiltonian discontinuous Galerkin finite element methods

Abstract A port-Hamiltonian (pH) system formulation is a geometrical notion used to formulate conservation laws for various physical systems. The distributed parameter port-Hamiltonian formulation models infinite dimensional Hamiltonian dynamical systems that have a nonzero energy flow through the boundaries. In this paper, we propose a novel framework for discontinuous Galerkin (DG) discretizations of pH-systems. Linking DG methods with pH-systems gives rise to compatible structure preserving semidiscrete finite element discretizations along with flexibility in terms of geometry and function spaces of the variables involved. Moreover, the port-Hamiltonian formulation makes boundary ports explicit, which makes the choice of structure and power preserving numerical fluxes easier. We state the Discontinuous Finite Element Stokes–Dirac structure with a power preserving coupling between elements, which provides the mathematical framework for a large class of pH discontinuous Galerkin discretizations. We also provide an a priori error analysis for the port-Hamiltonian discontinuous Galerkin Finite Element Method (pH-DGFEM). The port-Hamiltonian discontinuous Galerkin finite element method is demonstrated for the scalar wave equation showing optimal rates of convergence.

Numerical Investigation of Some Reductions for the Gatenby–Gawlinski Model

Two (consecutive) reductions of the complete Gatenby–Gawlinski model for cancer invasion are proposed in order to investigate the mathematical framework, mainly from a computational perspective. After a brief overview of the full model, we proceed by examining the case of a two-equations-based and one-equation-based reduction, both obtained by means of a quasi-steady-state assumption. We focus on invasion fronts, exploiting a numerical strategy based on a finite volume approximation, and perform corresponding computational simulations to study the sharpness/smoothness of the traveling waves. Then, we employ a space-averaged wave speed estimate—referred to as the LeVeque–Yee formula—to quantitatively approach the propagation phenomenon. Concerning the one-equation-based model, we propose a scalar degenerate reaction-diffusion equation, which proves to be effective in order to qualitatively recover the typical trends arising from the Gatenby–Gawlinski model. Finally, we carry out some numerical tests in a specific case where the analytical solution is available.

Numerical Analysis for Sturm–Liouville Problems with Nonlocal Generalized Boundary Conditions

For the generalized Sturm–Liouville problem (GSLP), a new formulation is undertaken to reduce the number of unknowns from two to one in the target equation for the determination of eigenvalue. The eigenparameter-dependent shape functions are derived for using in a variable transformation, such that the GSLP becomes an initial value problem for a new variable. For the uniqueness of eigenfunction an extra condition is imposed, which renders the right-end value of the new variable available; a derived implicit nonlinear equation is solved by an iterative method without using the differential; we can achieve highly precise eigenvalues. For the nonlocal Sturm–Liouville problem (NSLP), we consider two types of integral boundary conditions on the right end. For the first type of NSLP we can prove sufficient conditions for the positiveness of the eigenvalue. Negative eigenvalues and multiple solutions may exist for the second type of NSLP. We propose a boundary shape function method, a two-dimensional fixed-quasi-Newton method and a combination of them to solve the NSLP with fast convergence and high accuracy. From the aspect of numerical analysis the initial value problem of ordinary differential equations and scalar nonlinear equations are more easily treated than the original GSLP and NSLP, which is the main novelty of the paper to provide the mathematically equivalent and simpler mediums to determine the eigenvalues and eigenfunctions.

Scalar BSDEs of iterated-logarithmically sub-linear generators with integrable terminal values

We establish a general existence and uniqueness result of integrable adapted solutions to scalar backward stochastic differential equations with integrable parameters, where the generator g has an iterated-logarithmically uniform continuity in the second unknown variable z. The result improves our previous one in Fan et al. (2023), and can be applied to the (conditional) g-expectation and the star-shaped risk measure.

On hierarchical competition through reduction of individual growth

We consider a population organised hierarchically with respect to size in such a way that the growth rate of each individual depends only on the presence of larger individuals. As a concrete example one might think of a forest, in which the incidence of light on a tree (and hence how fast it grows) is affected by shading by taller trees. The classic formulation of a model for such a size-structured population employs a first order quasi-linear partial differential equation equipped with a non-local boundary condition. However, the model can also be formulated as a delay equation, more specifically a scalar renewal equation, for the population birth rate. After discussing the well-posedness of the delay formulation, we analyse how many stationary birth rates the equation can have in terms of the functional parameters of the model. In particular we show that, under reasonable and rather general assumptions, only one stationary birth rate can exist besides the trivial one (associated to the state in which there are no individuals and the population birth rate is zero). We give conditions for this non-trivial stationary birth rate to exist and analyse its stability using the principle of linearised stability for delay equations. Finally, we relate the results to the alternative, partial differential equation formulation of the model.

Frozen Gaussian sampling for scalar wave equations

Frozen Gaussian sampling for scalar wave equations

Bubble dynamics under the influence of the Marangoni force induced by a stratified field of contamination

AbstractThe Marangoni effect assumes significance in bubbly flows when temperature or concentration gradients exist in the domain. This study investigated the hydrodynamics of single bubbles under the influence of the Marangoni force induced by stratified fields of dissolved sugar, providing a numerical framework for examining these phenomena. A laboratory-scale bubble column and high-speed imaging were utilized to analyze the bubble behavior. The OpenFOAM-based geometric volume of the fluid solver was extended by incorporating the solutocapillary Marangoni effect, and a passive scalar transport equation for the sugar concentration was solved. The results revealed that small bubbles entering regions with elevated sugar concentrations experienced deceleration, transitioning into linear paths, while those departing from regions with high sugar concentrations exhibited fluctuations and meandering. Furthermore, the concentration gradient leads larger bubbles to meander throughout the entire column, without a notable increase in their velocity. The intensity of these behaviors is governed by the magnitude of the Marangoni force. The findings provide a better understanding of single bubble hydrodynamics in complex environments.

Spreading speeds for time heterogeneous prey-predator systems with nonlocal diffusion on a lattice

We investigate the spreading behaviour for the solutions of a non-autonomous prey-predator system on a discrete lattice. These time variations are assumed to enjoy an averaging property. This includes periodicity, almost periodicity and unique ergodicity as special cases. The spatial motion of individuals from one site to another is modelled by a discrete convolution operator. In order to take into account external fluctuations such as seasonality, daily variations and so on, the convolution kernels and reaction terms may vary with time. Our analysis of the spreading speeds of invasion of the species is based on the careful and detailed study of the hair-trigger effect and spreading speed for a non-autonomous scalar Fisher-KPP equation on a lattice. Then, we are able to compare the solutions of the prey-predator system with those of a suitable scalar Fisher-KPP equation and derive the invasion speeds of the prey and of the predator.

Normal forms of a class of partial functional differential equations

In this paper, we investigate the normal form a class of partial functional differential equations, which takes account of delays in the diffusion terms as well as a wider scope of nonlinear terms. We first study the associated linear theory, mainly including the spectral properties of infinitesimal generator, formal adjoint and decomposition of phase space. Based on these results, the normal form theory for nonlinear equation is established, which can be used to study the local dynamics near the steady state for such equations. As an application, we consider the Hopf bifurcation problem of a scalar diffusive equation with delay not only involved in the diffusive term but also in reaction terms. The normal form, depending on the original coefficients, up to the third order term is calculated, which allows us to determine the direction of Hopf bifurcation and stability of bifurcated periodic solutions. The results are then applied to study the scalar population model with memory-based diffusion for modeling the movement of highly-developed animals.

Parallel and energy conservative/dissipative schemes for sine–Gordon and Allen–Cahn equations

The sine–Gordon and Allen–Cahn equations are two typical models in the fields of conservative Hamiltonian systems and dissipative gradient flows, respectively. As the demand for numerical methods that respect intrinsic energy conservation/dissipation laws turns into a fundamental principle, the subsequent computational efficiency is getting more and more desired. Linearly implicit methods, which only require solutions of linear algebraic systems at each time step, have become a popular way to design efficient schemes. However, the overall computational cost of these methods depends heavily on the speed at which the linear system can be solved. In this paper, we propose an alternative approach for developing highly efficient energy-conservative or dissipative schemes for the sine–Gordon and Allen–Cahn equations. Our approach is based on spatial finite difference approximations and is fully implicit at first glance, but actually it is completely decoupled point by point, allowing for the implementation of a scalar nonlinear equation successively with a lower complexity that is comparable only to the degrees of freedom. We further establish the connections between our approach and the classic Itoh–Abe discrete gradient and the partitioned averaged vector field methods, and then propose the generalized Itoh–Abe discrete gradient method, which offers great flexibility in the ordering of updating each unknown point. One immediate benefit is that we can specifically choose an ordering to achieve a naturally parallel computation of the resulting scheme, significantly improving the computational efficiency. Various numerical experiments are presented to illustrate the performance of the proposed schemes.

Nonparaxial accelerating waves as a superposition of nondiffracting Bessel-lattice optical fields.

In the first part of this work, we introduce a monochromatic solution to the scalar wave equation in free space, defined by a superposition of monochromatic nondiffracting half Bessel-lattice optical fields, which is determined by two scalar functions; one is defined on frequency space, and the other is a complete integral to the eikonal equation in free space. We obtain expressions for the geometrical wavefronts, the caustic region, and the Poynting vector. We highlight that this solution is stable under small perturbations because it is characterized by a caustic of the hyperbolic umbilical type. In the second part, we introduce the corresponding solution to the Maxwell equations in free space.

Implicit-explicit time integration for the immersed wave equation

Efficient solution strategies for wave propagation problems with complex geometries are essential in many engineering fields. Immersed boundary methods simplify mesh generation by embedding the domain of interest into an extended domain that is easy to mesh, introducing the challenge of dealing with cells that intersect the domain boundary. We consider the finite cell method that extends the weak form from the physical into the extended domain, multiplied by a small value, to stabilize badly cut cells. Combined with explicit time integration methods, the presence of finite cells with very little support in the physical domain results in tiny critical time step sizes. While the finite cell stabilization limits how small the critical time step size can become, this limit is still restrictive for many applications. Explicit transient analyses commonly use the spectral element method due to its natural way of obtaining diagonal mass matrices through nodal lumping. The resulting method is very efficient since nodal lumping renders the solution of the equation systems trivial while not reducing the accuracy or the critical time step size. The combination of the spectral element and finite cell methods is called the spectral cell method. Unfortunately, a direct application of nodal lumping in the spectral cell method is impossible due to the special quadrature necessary to treat the discontinuous integrand inside the cut cells. The existing approaches to lump the mass matrices of cut cells significantly reduce the accuracy of the approximation.We analyze an implicit-explicit (IMEX) time integration method to exploit the advantages of the nodal lumping scheme for uncut cells on one side and the unconditional stability of implicit time integration schemes for cut cells on the other. In this hybrid, immersed Newmark IMEX approach, we use explicit second-order central differences to integrate the uncut degrees of freedom that lead to a diagonal block in the mass matrix and an implicit trapezoidal Newmark method to integrate the remaining degrees of freedom (those supported by at least one cut cell). The immersed Newmark IMEX approach preserves the high-order convergence rates and the geometric flexibility of the finite cell method while retaining the efficiency of the nodal lumping scheme for uncut cells without compromising the critical time step size. We analyze a simple system of spring-coupled masses to highlight some of the essential characteristics of Newmark IMEX time integration. We then solve the scalar wave equation with immersed boundaries on two- and three-dimensional examples with significant geometric complexity to show that our approach is more efficient than state-of-the-art time integration schemes when comparing accuracy and runtime. While we focus on the finite and spectral cell methods, we expect other immersed methods to benefit equally from Newmark IMEX time integration.

Asymptotic Errors in the Superconvergence of Discontinuous Galerkin Methods Based on Upwind-Biased Fluxes for 1D Linear Hyperbolic Equations

AbstractIn this paper, we study the superconvergence of the semi-discrete discontinuous Galerkin (DG) method for linear hyperbolic equations in one spatial dimension. The asymptotic errors in cell averages, downwind point values, and the postprocessed solution are derived for the initial discretization by Gaussian projection (for periodic boundary condition) or Cao projection Cao et al. (SIAM J. Numer. Anal. 5, 2555–2573 (2014)) (for Dirichlet boundary condition). We proved that the error constant in the superconvergence of order $$2k+1$$ 2 k + 1 for DG methods based on upwind-biased fluxes depends on the parity of the order k. The asymptotic errors are demonstrated by various numerical experiments for scalar and vector hyperbolic equations.

The existence and multiplicity of L 2-normalized solutions to nonlinear Schrödinger equations with variable coefficients

Abstract The existence of L 2–normalized solutions is studied for the equation − Δ u + μ u = f ( x , u ) in R N , ∫ R N u 2 d x = m . $-{\Delta}u+\mu u=f\left(x,u\right)\quad \quad \text{in} {\mathbf{R}}^{N},\quad {\int }_{{\mathbf{R}}^{N}}{u}^{2} \mathrm{d}x=m.$ Here m > 0 and f(x, s) are given, f(x, s) has the L 2-subcritical growth and (μ, u) ∈ R × H 1(R N ) are unknown. In this paper, we employ the argument in Hirata and Tanaka (“Nonlinear scalar field equations with L 2 constraint: mountain pass and symmetric mountain pass approaches,” Adv. Nonlinear Stud., vol. 19, no. 2, pp. 263–290, 2019) and find critical points of the Lagrangian function. To obtain critical points of the Lagrangian function, we use the Palais–Smale–Cerami condition instead of Condition (PSP) in Hirata and Tanaka (“Nonlinear scalar field equations with L 2 constraint: mountain pass and symmetric mountain pass approaches,” Adv. Nonlinear Stud., vol. 19, no. 2, pp. 263–290, 2019). We also prove the multiplicity result under the radial symmetry.