Time-dependent Solutions Research Articles

The Efficacy of Haar Wavelets in Addressing Discontinuities of McKean Equations with Heaviside Functions

This study explores the application of Haar wavelets to solve the McKean equation, a reaction-diffusion equation with discontinuous Heaviside step function. Haar wavelets, with their compact support and orthogonality, offer straightforward but yet powerful tool for addressing the equation’s nonlinear dynamics. We focus on the time-independent solution of the McKean equation, which is crucial for understanding the threshold phenomenon that determines the system’s behavior. Despite the existence of analytical time-independent solution to the McKean equation, achieving such solutions in closed form is uncommon for more complicated systems, highlighting the utility of the Haar wavelet approach. The proposed method integrates the Haar series expansion of the highest order derivative, enabling systematic solution derivation. Through a detailed comparison with analytical solution, we validate the Haar wavelet approach as a robust and computationally feasible tool for solving complex reaction-diffusion system. The results also demonstrate the method’s accuracy and efficiency, offering insights into its broader applicability to more complex reaction-diffusion system, especially those with discontinuity and sharp transitions.

Optimal local truncation error method on unfitted Cartesian meshes for solution of 3-D wave and heat equations for heterogeneous materials

In the paper we develop the optimal local truncation error method (OLTEM) with the non-diagonal and diagonal mass matrices on unfitted Cartesian meshes for the 3-D time-dependent wave and heat equations for heterogeneous materials with irregular interfaces. 27-point stencils that are similar to those for linear finite elements are used with OLTEM. There are no unknowns for OLTEM on interfaces between different materials; the structure of the global discrete equations is the same for homogeneous and heterogeneous materials. The calculation of the unknown stencil coefficients is based on the minimization of the local truncation error of the stencil equations, includes the use of the interface conditions and the entire PDE in derivations, and yields the optimal fourth order of accuracy of OLTEM with the non-diagonal mass matrix. This is two order higher than for linear finite elements with similar 27-point stencils. OLTEM with the diagonal mass matrix includes the rigorous calculation of the diagonal mass matrix and provides the optimal second order of accuracy. The 3-D numerical results for heterogeneous materials with irregular interfaces show that at the same number of degrees of freedom, OLTEM is even much more accurate than high-order (up to the 7-th order - the highest order in the commercial code COMSOL) finite elements with much wider stencils. Compared to linear finite elements with similar 27-point stencils, at an accuracy of 0.1% OLTEM decreases the number of degrees of freedom by 550–7300 times. This leads to a huge reduction in computation time.A new 3-D post-processing procedure has been developed with OLTEM for the calculation of the spatial derivatives of time-dependent numerical solutions. The spatial derivatives for each grid point are calculated with the help of one compact 27-point stencil (the same as that in basic computations). In contrast to known post-processing procedures, the new approach includes also the use of the time derivatives of the function, the interface conditions and the original PDE. The spatial derivatives of the OLTEM solutions calculated with the new post-processing procedure are much more accurate compared to those obtained by high-order (up to the 7-th order) finite elements with much wider stencils. At an accuracy of 0.1% for the spatial derivatives, OLTEM decreases the number of degrees of freedom by 1.9⋅106–6⋅109 times compared to linear finite elements. The new post-processing procedure can be equally applied to the calculation of the partial derivatives obtained by other numerical methods as well as to the numerical results for other PDEs.Due to the huge reduction in the computation time compared to existing methods and the use of trivial unfitted Cartesian meshes that are independent of irregular geometry, the proposed technique does not require remeshing for the shape change of irregular geometry and it will be effective for many design and optimization problems as well as for multiscale problems without the scale separation.

Modelling the evolution of an ice sheet’s weathering crust

Abstract The weathering crust is a layer of porous ice that can form at the surface of an ice sheet. It grows and decays in response changing weather and climate conditions, affecting the albedo, the melt rate and the transport of meltwater across the surface. To understand this behaviour, we seek time-dependent solutions to a continuum, thermodynamic model for the porosity, temperature and thickness of the weathering crust, and the internal and surface melt rates. We find solutions using a numerical enthalpy method, presented in this study. We use idealized ‘switching’ and sinusoidal forcings to explore the different dynamics exhibited during growth and decay, the timescales involved, and the impact of diurnal vs. annual variations. The results demonstrate qualitative agreement with observations, and provide insight into the relative importance of different surface heat fluxes during the growth and decay of the crust. The model therefore provides a useful tool for exploring the response of the weathering crust to climate change.

Exact Time-Dependent Thermodynamic Relationships for a Brownian Particle Navigating Complex Networks

The thermodynamic characteristics of systems driven out of equilibrium are examined for M Brownian ratchets organized in a complex network. The precise time-dependent solution reveals that the entropy S, entropy production ep(t), and entropy extraction hd(t) of the system in complex networks increase with system size, which is plausible as these thermodynamic quantities display extensive properties. In other words, as the number of lattice sites increases, the entropy S, entropy production ep(t), and entropy extraction hd(t) increase, demonstrating that these complex networks cannot be reduced to the corresponding one-dimensional lattice. Conversely, the rates for thermodynamic quantities such as velocity V, entropy production rate e˙p(t), and entropy extraction rate h˙ d(t) become independent of the network size in the long-term limit. The exact analytic results also indicate that the free energy decreases with system size. The model system is further analyzed by incorporating heat transfer via kinetic energy. Since heat exchange via kinetic energy does affect the energy extraction rate, the heat dumped to the cold reservoirs also contributes to internal entropy production. Consequently, such systems exhibit a higher degree of irreversibility. The thermodynamic features of a system operating between hot and cold baths are also compared and contrasted with a system functioning in a heat bath where temperature varies linearly along the reaction coordinate. Regardless of the network arrangements, the entropy, entropy production, and extraction rates are significantly larger for the linearly varying temperature case than for a system operating between hot and cold baths.

Cellular diffusion processes in singularly perturbed domains

There are many processes in cell biology that can be modeled in terms of particles diffusing in a two-dimensional (2D) or three-dimensional (3D) bounded domain Ω⊂Rd\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varOmega \\subset {\\mathbb {R}}^d$$\\end{document} containing a set of small subdomains or interior compartments Uj\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${\\mathcal {U}}_j$$\\end{document}, j=1,…,N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$j=1,\\ldots ,N$$\\end{document} (singularly-perturbed diffusion problems). The domain Ω\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varOmega $$\\end{document} could represent the cell membrane, the cell cytoplasm, the cell nucleus or the extracellular volume, while an individual compartment could represent a synapse, a membrane protein cluster, a biological condensate, or a quorum sensing bacterial cell. In this review we use a combination of matched asymptotic analysis and Green’s function methods to solve a general type of singular boundary value problems (BVP) in 2D and 3D, in which an inhomogeneous Robin condition is imposed on each interior boundary ∂Uj\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\partial {\\mathcal {U}}_j$$\\end{document}. This allows us to incorporate a variety of previous studies of singularly perturbed diffusion problems into a single mathematical modeling framework. We mainly focus on steady-state solutions and the approach to steady-state, but also highlight some of the current challenges in dealing with time-dependent solutions and randomly switching processes.

Axisymmetric consolidation behavior of multilayered unsaturated soils with transversely isotropic permeability

Time-dependent consolidation behavior of unsaturated soils is a vital problem in the geotechnical engineering. With the aid of the Fredlund consolidation theory, this work further assumes the total stress of soils skeleton freely change, and extends the Fredlund consolidation theory to a Biot-type theory, establishing the fully-coupled equation model of multilayered unsaturated poroelastic media with transversely isotropic permeability. To convert the partial differential governing equation into ordinary differential equations, the integration transform technology is applied. Subsequently, the precise integration method is used to acquire the time-dependent consolidation solution of multilayered unsaturated media with transversely isotropic permeability in the transformed domain, which is further solved in the actual domain by the inverse Hankel transform. A verification examples is provided to compare the present results with the existing work in the literature, showing a great coincidence and proving the feasibility of the present solution. Finally, numerous numerical examples are presented to investigate the evolution of excess pore pressure and settlement under quasi-static loads, revealing the consolidation behavior of unsaturated soils. The results demonstrates that the ramping time, stratification, permeability, depth and m1w have a significant effect on the consolidation behavior.

Well-posedness and stability of a stochastic neural field in the form of a partial differential equation

A system of partial differential equations representing stochastic neural fields was recently proposed with the aim of modelling the activity of noisy grid cells when a mammal travels through physical space. The system was rigorously derived from a stochastic particle system and its noise-driven pattern-forming bifurcations have been characterised. However, due to its nonlinear and non-local nature, standard well-posedness theory for smooth time-dependent solutions of parabolic equations does not apply. In this article, we transform the problem through a suitable change of variables into a nonlinear Stefan-like free boundary problem and use its representation formulae to construct local-in-time smooth solutions under mild hypotheses. Then, we prove that fast-decaying initial conditions and globally Lipschitz modulation functions imply that solutions are global-in-time. Last, we improve previous linear stability results by showing nonlinear asymptotic stability of stationary solutions under suitable assumptions.

Timelike Kasner singularities and Floquet states in 2+1d AdS/CFT

We consider a model of a holographic 2+1d CFT interacting with an oscillating background gauge field. It is solved by an AdS-Vaidya metric describing Ohmic heating of the boundary field theory. However, we also show that if timelike singularities of Kasner type are permitted then a time independent solution that may be interpreted as a Floquet state of the system can be constructed. In this state the system exhibits either Hall conductivity or kinetic induction, and we numerically evaluate the Kasner exponents for a range of boundary conditions. This model may contribute to the ongoing discussion on the validity and meaning of the Kasner metric in the AdS/CFT correspondence and its application in cosmology.

The modelling error in multi-dimensional time-dependent solute transport models

Starting from full-dimensional models of solute transport, we derive and analyze multi-dimensional models of time-dependent convection, diffusion, and exchange in and around pulsating vascular and perivascular networks. These models are widely applicable for modelling transport in vascularized tissue, brain perivascular spaces, vascular plants and similar environments. We show the existence and uniqueness of solutions to both the full- and the multi-dimensional equations under suitable assumptions on the domain velocity. Moreover, we quantify the associated modelling errors by establishing a-priori estimates in evolving Bochner spaces. In particular, we show that the modelling error decreases with the characteristic vessel diameter and thus vanishes for infinitely slender vessels. Numerical tests in idealized geometries corroborate and extend upon our theoretical findings.

On an Oberbeck-Boussinesq model relating to the motion of a viscous fluid subject to heating

Abstract This article surveys some results in the study of Iannelli [Su un modello di Oberbeck-Boussinesq relativo al moto di un fluido viscoso soggetto a riscaldamento, Fisica Matematica, Istituto Lombardo (rend. Sc.) A 121 (1987), 145–191], in which the motion of a viscous, compressible fluid in a two-dimensional domain, subject to heating at the walls, is studied. A global existence and uniqueness theorem for the time-dependent problem is given, and also, under more stringent assumptions, an existence and uniqueness theorem in the stationary case is given. A theorem on the asymptotic behavior for t → ∞ t\to \infty of the time-dependent solutions is proved.

Two-Layer Equilibrium Model of Miscible Inhomogeneous Fluid Flow

Two-layer flow of a density-stratified fluid with mass transfer between the layers is considered. In the Boussinesq approximation, the equations of motion are reduced to a homogeneous quasilinear system of partial differential equations of mixed type. The flow parameters in the intermediate mixed layer are determined from the equilibrium conditions in a more general model of three-layer flow of a miscible fluid. In particular, the equilibrium conditions imply the constancy of the interlayer Richardson number in velocity-shift flows. A self-similar solution to the problem of breakdown of an arbitrary discontinuity (the lock-exchange problem) in the domain of hyperbolicity of the system under consideration is constructed. The transcritical flow regimes over a local obstacle are studied and the conditions under which the obstacle determines the upstream flow are determined. A comparison of steady-state and time-dependent solutions with the solutions obtained for the original three-layer models of miscible fluid flow is carried out.

The combined effects of quantum and strong coupling on the nonlinear collective excitation in quantum dusty plasmas

The quantum hydrodynamic model for electrons and ions and the generalized hydrodynamic model for the strongly coupled dust particles are proposed in the strongly coupled quantum dusty plasma, where the combined quantum effects of quantum diffraction, quantum statistic pressure, as well as electron exchange and correlation effects are all considered in the quantum hydrodynamic model. The shear and bulk viscosity effects are included in the viscoelastic relaxation, which leads to the decay of the dust-ion-acoustic waves. The approximate time-dependent solitary solution is obtained by the momentum conservation law in the presence of viscosity.

Meridional circulation streamlined

Time-dependent meridional circulation and differential rotation in radiative zones are central open issues in stellar evolution theory. We streamline this challenging problem using the “downward control principle” of atmospheric science, under a geostrophic f-plane approximation. We recover the known stellar physics result that the steady-state meridional circulation decays on the length scale proportional to N / f × Pr , assuming molecular viscosity is the dominant drag mechanism. Prior to steady-state, the meridional circulation and the zonal wind (= differential rotation) spread together via radiative diffusion, under thermal wind balance. The corresponding (fourth-order) hyperdiffusion process is reasonably well approximated by regular (second-order) diffusion on scales of order a pressure scale-height. We derive an inhomogeneous diffusion (equiv. advection-diffusion) equation for the zonal flow which admits closed-form time-dependent solutions in a finite depth domain, allowing for rapid prototyping of differential rotation profiles. In the weak drag limit, we find that the time to rotational steady-state can be longer than the Eddington-Sweet time and be instead determined by the longer drag time. Unless strong enough drag operates, the internal rotation of main-sequence stars may thus never reach a steady-state. Streamlined meridional circulation solutions provide leading-order internal rotation profiles for studying the role of fluid/MHD instabilities (or waves) in redistributing angular momentum in the radiative zones of stars. Despite clear geometrical limitations and simplifying assumptions, one might expect our thin-layer geostrophic approach to offer qualitatively useful results to understand deep meridional circulation in stars.

Spherical symmetric solutions of conformal Killing gravity: black holes, wormholes, and sourceless cosmologies

Abstract The most general set of static and spherically symmetric solutions for conformal Killing gravity coupled to Maxwell fields is presented in closed form. These solutions, depending on six parameters, include non-asymptotically flat black holes or naked singularities, non-asymptotically flat traversable wormholes, and (possibly singularity-free) closed universes. We also consider the inverse problem, showing that the most general energy–momentum tensor generating a given static spherically symmetric metric depends on three parameters. Sourceless time-dependent isotropic solutions are also given. These solutions depending on the curvature, the cosmological constant, and a new integration constant α, present a rich variety, including singularity-free eternal cosmologies and universes evolving symmetrically from a big bang to a big crunch within a finite lapse of time.

The Effect of Nonlinear Sloshing Response of Water on Seismic Behavior of Reinforced Concrete Elevated Water Tanks

In this study, the fluid-structure interactive behavior by using the Westergaard approach and the Smoothed Particle Hydrodynamics (SPH) method of a 1000 m3 reinforced concrete (RC) elevated water tank under different seismic activities i.e., Kocaeli, Van, Kahramanmaraş and Kobe earthquakes was investigated. In the SPH method, the resulting hydrodynamic pressures were applied to the interior of the wall of the structure as an external force which is made as mass additions to the tank walls in Westergaard approach. With the help of software developed, the load carrying system of the elevated water tank was modeled using the MATLAB Partial Differential Equation Toolbox (PDE) software, based on the finite element method, and its verification was made by comparing with ANSYS software model. Dynamic condensation method was used to perform time history analyses for the cited earthquakes. The time-dependent solutions of partial differential equations are solved with the help of ODE45 functions based on Runge-Kutta method, by arranging the equation of motion in state-space format. The cited tank was analyzed linearly under different seismic activities for empty, half-filled (50%) and full-filled (100%) cases and the obtained results were compared with each other. It has been concluded that considering nonlinear behavior of the fluid during the earthquake with the SPH method provided consistency with the real behavior of RC elevated water tanks and gave more realistic results than the traditional Westergaard approach.

Fisher–KPP-type models of biological invasion: open source computational tools, key concepts and analysis

This review provides open-access computational tools that support a range of mathematical approaches to analyse three related scalar reaction–diffusion models used to study biological invasion. Starting with the classic Fisher–Kolomogorov (Fisher–KPP) model, we illustrate how computational methods can be used to explore time-dependent partial differential equation (PDE) solutions in parallel with phase plane and regular perturbation techniques to explore invading travelling wave solutions moving with dimensionless speed c ≥ 2 . To overcome the lack of a well-defined sharp front in solutions of the Fisher–KPP model, we also review two alternative modelling approaches. The first is the Porous–Fisher model where the linear diffusion term is replaced with a degenerate nonlinear diffusion term. Using phase plane and regular perturbation methods, we explore the distinction between sharp- and smooth-fronted invading travelling waves that move with dimensionless speed c ≥ 1 / 2 . The second alternative approach is to reformulate the Fisher–KPP model as a moving boundary problem on 0 < x < L ( t ) , leading to the Fisher–Stefan model with sharp-fronted travelling wave solutions arising from a PDE model with a linear diffusion term. Time-dependent PDE solutions and phase plane methods show that travelling wave solutions of the Fisher–Stefan model can describe both biological invasion ( c > 0 ) and biological recession ( c < 0 ) . Open source Julia code to replicate all computational results in this review is available on GitHub; we encourage researchers to use this code directly or to adapt the code as required for more complicated models.

Time-dependent dynamical energy analysis via convolution quadrature

Dynamical Energy Analysis was introduced in 2009 as a novel method for predicting high-frequency acoustic and vibrational energy distributions in complex engineering structures. In this paper we introduce the first time-dependent Dynamical Energy Analysis method. Time-domain models are important in numerous applications including sound simulation in room acoustics, predicting shock-responses in structural mechanics and modelling electromagnetic scattering from conductors. The first step is to reformulate Dynamical Energy Analysis in the time-domain by means of a convolution integral operator. We are then able to employ the Convolution Quadrature method to provide a link between the previous frequency-domain implementations of Dynamical Energy Analysis and fully time-dependent solutions by means of the Z-transform. By combining a modified multistep Convolution Quadrature approach for the time discretisation, together with Galerkin and Petrov-Galerkin methods for the space and momentum discretisations, respectively, we are able to accurately track the propagation of high-frequency transient signals through phase-space. The implementation here is detailed for finite two-dimensional spatial domains and we demonstrate the versatility of our approach by performing a range of numerical experiments for regular, non-convex and irregular geometries as well as different types of wave source.

C0–Semigroups Approach to the Reliability Model Based on Robot-Safety System

This paper considers a system with one robot and n safety units (one of which works while the others remain on standby), which is described by an integro-deferential equation. The system can fail in the following three ways: fails with an incident, fails safely and fails due to the malfunction of the robot. Using the C0–semigroups theory of linear operators, we first show that the system has a unique non-negative, time-dependent solution. Then, we obtain the exponential convergence of the time-dependent solution to its steady-state solution. In addition, we study the asymptotic behavior of some time-dependent reliability indices and present a numerical example demonstrating the effects of different parameters on the system.

Dynamical branes on expanding orbifold and complex projective space

We construct some new dynamical p-brane solutions to gravity theories on curved backgrounds. We discuss the relations between dynamical branes, a new time-dependent solution on complex projective space CPn, and the static p-branes on the orbifold Cn/Zn. Published by the American Physical Society 2024

Two Models on the Unsteady Heat and Fluid Flow Induced by Stretching or Shrinking Sheets and Novel Time-Dependent Solutions

Abstract This study sheds light on unsteady heat and fluid flow problems over stretching and shrinking surfaces, enriching our understanding of these complex phenomena. We derive two mathematical models using a rigorous approach. The first model aligns with the model commonly employed by researchers in this field, but its steady-state solution remains trivial. The second model, introduced in this work, demonstrably captures the physically relevant steady-state solutions well-studied in the literature. Notably, we introduce new similarity solutions for the temperature field specifically within the first model. We further demonstrate that a uniform wall temperature condition leads to the optimal heat transfer rate. While similarity solutions can be derived for specific cases with the second model, nonsimilar solutions may be necessary for more general scenarios. We discuss the implications of our analysis for stagnation-point flow and non-Newtonian viscoelastic fluid flow problems, illuminating future research directions in the open literature.