Classical Equations Research Articles

On coupled semilinear evolution systems: Techniques on fractional powers of 4×4$4\times 4$ matrices and applications

AbstractIn this paper, we provide several techniques to explicitly calculate fractional powers of operator matrices focusing on creating a theory that can be applied to distinct situations. To illustrate the abstract results developed, we consider its application in systems of coupled reaction–diffusion equations and in (strongly damped) wave equations. We also discuss how these techniques can be applied to higher order matrices and we specifically calculate the fractional powers of a operator matrix associated to a weakly coupled system of wave equation. In addition, we deal with the applicability of this analysis with respect to solvability, stabilization, regularity, smooth dynamics, and connection with evolutionary classic equation and its fractional counterpart.

Well-Posedness and Regularity of Solutions to Neural Field Problems with Dendritic Processing

We study solutions to a recently proposed neural field model in which dendrites are modelled as a continuum of vertical fibres stemming from a somatic layer. Since voltage propagates along the dendritic direction via a cable equation with nonlocal sources, the model features an anisotropic diffusion operator, as well as an integral term for synaptic coupling. The corresponding Cauchy problem is thus markedly different from classical neural field equations. We prove that the weak formulation of the problem admits a unique solution, with embedding estimates similar to the ones of nonlinear local reaction–diffusion equations. Our analysis relies on perturbing weak solutions to the diffusion-less problem, that is, a standard neural field, for which weak problems have not been studied to date. We find rigorous asymptotic estimates for the problem with and without diffusion, and prove that the solutions of the two models stay close, in a suitable norm, on finite time intervals. We provide numerical evidence of our perturbative results.

A peridynamic differential operator modeling approach for ceramic matrix composites microstructure with uniform or non-uniform discretization

Based on the theory of peridynamic differential operator (PDDO), a novel three-dimensional (3D) peridynamic (PD) model for isotropic and orthotropic material is proposed to investigate the elastic mechanical behaviors and stress distribution of Ceramic Matrix Composites (CMC) microstructure. The non-local expressions of strain and stress tensors are derived by employing PDDO and the classical constitutive equations. The bond forces in interface region crossing two different constituent materials are obtained by converting the partial differential terms into the non-local forms. The weak form of PD equations of motion is derived by using the principle of virtual work and PDDO. The validity of this approach is verified by predicting the elastic response and stress distributions of isotropic and orthotropic materials under uniaxial tension and pure shear deformation. Finally, simulations are conducted on a CMC microstructure with fiber and matrix under periodic boundary conditions. It is demonstrated that the current approach can effectively build 3D CMC microstructure with uniform or non-uniform discretization, and accurately predict the stress fields and effective elastic properties.

DECAY RATES AND INITIAL VALUES FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATIONS

We consider a solution u(·,t) to an initial boundary value problem for time-fractional diffusion-wave equation with the order α ∈ (0,2)\ {1} where t is a time variable. We first prove that a suitable norm of u(·,t) is bounded by the rate t −α for 0 α 1 and t 1−α for 1 α 2 for all large t > 0. Second, we characterize initial values in the cases where the decay rates are faster than the above critical exponents. Differently from the classical diffusion equation α = 1, the decay rate can keep some local characterization of initial values. The proof is based on the eigenfunction expansions of solutions and the asymptotic expansions of the Mittag-Leffler functions for large time.

Nonplanar ion acoustic waves in e-p-i plasmas with non-Maxwellian distributed electrons and positrons

The process of propagation of cylindrical and spherical ion acoustic solitary waves in an electron, positron, and ion plasma system is studied. The electrons are considered to have the generalized (r, q) distribution, while the positrons are supposed to have the kappa distribution, also known as the Lorentzian distribution. The classic nonlinear equation, Korteweg de-Vries (KdV), has been numerically solved by employing appropriate mathematical transformations. The role of various plasma parameters and nonplanar geometry on the widths and amplitudes of nonlinear ion acoustic structures has been numerically investigated. The current study could be helpful to understand the significance of cylindrical and spherical solitary waves in laboratory, astrophysical as well as in space plasmas.

Galerkin-finite difference method for fractional parabolic partial differential equations

The fractional form of the classical diffusion equation embodies the super-diffusive and sub-diffusive characteristics of any flow, depending on the fractional order. This study aims to approximate the solution of parabolic partial differential equations of fractional order in time and space. For this, firstly, we briefly discussed on some existing methods to solve Partial Differential Equations (PDEs) and fractional Differential Equations (DEs), then introduce a combined technique such as the Galerkin weighted residual method for the space fractional term with modified Bernoulli polynomials as basis functions, and the finite difference approximation for the time fractional term, respectively. The mathematical formulation of the proposed method is explained elaborately. Then we describe the order of convergence for the time fractional term only, as the convergence of the Galerkin method is obvious. We impose this technique on the fractional Black-Scholes model subsequently. Finally, we experimented with our proposed technique on some numerical problems. All the results are depicted in both tabular and 3D visualizations as well. We compare our results with the available methods in the literature, and our accuracy is considerable. To summarize: •The paper introduces an approach that integrates the Galerkin weighted residual method with modified Bernoulli polynomials to handle space fractional terms, alongside employing a finite difference approximation for time fractional terms.•The convergence analysis is focused.•The technique is implemented on the fractional Black-Scholes model and other numerical problems, with outcomes depicted through tables and 3D visualizations.

Heterotic integrable deformation of the principal chiral model

A novel classically integrable model is proposed. It is a deformation of the two-dimensional principal chiral model, embedded into a heterotic σ model by a particular heterotic gauge field. This is inspired by the bosonic part of the heterotic σ model and its recent Hamiltonian formulation in terms of O(d,d+n)-generalized geometry in Hatsuda [Gauged double field theory, current algebras and heterotic sigma models, ]. Classical integrability is shown by the construction of a Lax pair and a classical R matrix. The latter is almost of the canonical form with a twist function and solves the classical Yang-Baxter equation. Published by the American Physical Society 2024

On the use of spherical harmonic approximations in nonclassical particle transport problems

Presented here is the application of the spherical harmonic expansion to the nonclassical transport equation to derive a system of equations for the nonclassical flux moments, namely nonclassical spherical harmonic approximations (NSHA). The nonclassical transport equation is a recently developed mathematical model that allows the modeling of transport problems wherein the particle flux does not undergo exponential attenuation. We employ a Spectral Approach technique to depict the nonclassical flux by expressing it as a truncated series of Laguerre polynomials with respect to the free-path variable s. This yields a system of equations for the nonclassical flux moments structured similarly to the classical PN equations. We show that the NSHA is reduced to the classical PN equation in a special case. Numerical results for slab-geometry test problems are given to verify the derivation.

Imaginary past of a quantum particle moving on imaginary time

The analytical continuation of classical equations of motion to complex times suggests that a tunneling particle spends in the barrier an imaginary duration i|T|. Does this mean that it takes a finite time to tunnel, or should tunneling be seen as an instantaneous process? It is well known that examination of the adiabatic limit in a small additional AC field points towards |T| being the time it takes to traverse the barrier. However, this is only half the story. We probe the transmitted particle's history, and find that it very little of the field's past behavior, as if the transit time were close to zero. The ensuing contradiction suggests that the question is ill posed, and we explain why. Published by the American Physical Society 2024

Quantum simulation of dissipation for Maxwell equations in dispersive media

The dissipative character of an electromagnetic medium breaks the unitary evolution structure that is present in lossless, dispersive optical media. In dispersive media, dissipation appears in the Schrödinger representation of classical Maxwell equations as a sparse diagonal operator occupying an r-dimensional subspace. A first order Suzuki-Trotter approximation for the evolution operator enables us to isolate the non-unitary operators (associated with dissipation) from the unitary operators (associated with lossless media). The unitary operators can be implemented through qubit lattice algorithm (QLA) on n qubits, based on the discretization and the dimensionality of the pertinent fields. However, the non-unitary-dissipative part poses a challenge both physically and computationally on how it should be implemented on a quantum computer. In this paper, two probabilistic dilation algorithms are considered for handling the dissipative operators. The first algorithm is based on treating the classical dissipation as a linear amplitude damping-type completely positive trace preserving (CPTP) quantum channel where an unspecified environment interacts with the system of interest and produces the non-unitary evolution. Therefore, the combined system-environment is now closed, and must undergo unitary evolution in the dilated space. The unspecified environment can be modeled by just one ancillary qubit, resulting in an implementation scaling of O(2n−1n2) elementary gates for the total system-environment unitary evolution operator. The second algorithm approximates the non-unitary operators by the Linear Combination of Unitaries (LCU). On exploiting the diagonal structure of the dissipation, we obtain an optimized representation of the non-unitary part, which requires O(2n) elementary gates. Applying the LCU method for a simple dielectric medium with homogeneous dissipation rate, the implementation scaling can be further reduced into O[poly(n)] basic gates. For the particular case of weak dissipation we show that our proposed post-selective dilation algorithms can efficiently delve into the transient evolution dynamics of dissipative systems by calculating the respective implementation circuit depth. A connection of our results with the non-linear-in-normalization-only (NINO) quantum channels is also presented.

[formula omitted]-player game formulation of the majority-vote model of opinion dynamics

From a self-centered perspective, it can be assumed that people only hold opinions that can benefit them. If opinions have no intrinsic value, and acquire their value when held by the majority of individuals in a discussion group, then we have a situation that can be modeled as an N-player game. Here we explore the dynamics of (binary) opinion formation using a game-theoretic framework to study an N-player game version of Galam’s local majority-vote model. The opinion dynamics is modeled by a stochastic imitation dynamics in which the individuals copy the opinion of more successful peers. In the infinite population limit, this dynamics is described by the classical replicator equation of evolutionary game theory. The equilibrium solution shows a threshold separating the initial frequencies that lead to the fixation of one opinion or the other. A comparison with Galam’s deterministic model reveals contrasting results, especially in the presence of inflexible individuals, who never change their opinions. In particular, the N-player game predicts a polarized equilibrium consisting only of extremists. Using finite-size scaling analysis, we evaluate the critical exponents that determine the population size dependence of the opinion’s fixation probability and mean fixation times near the threshold. The results underscore the usefulness of combining evolutionary game theory with opinion dynamics and the importance of statistical physics tools to summarize the results of Monte Carlo simulations.

Transmission of waves and particles through the interface: Reversibility and coherence

We examine the transmission of quantum particles (phonons, electrons, and photons) across interfaces, identifying universal patterns in diverse physical scenarios. Starting with classical wave equations, we quantize them and derive kinetic equations. Those are matching conditions for the distribution functions of particles at the interface. We note the time irreversibility of the derived kinetic equations — an essential feature for accurately describing irreversible processes like heat transport. We identify the juncture in our derivation where the time symmetry of wave equations is disrupted, it is the assumption of the non-coherence of incident waves. Consequently, we infer that non-coherent transmission through the interface exhibits time irreversibility. We propose an experiment to validate this hypothesis.

Two-wave interaction solutions of perturbation and CKdVE integrability for (2+1)-D CDGKS equation

By the Hirota bilinear method, some new interaction solutions with the complex perturbation for CDGKS equation are obtained. Meanwhile, with the help of the classical nonlinear KdV equation, many new exact solutions of CDGKS equation are derived through the CKdVE method, since it satisfies the CKdVE solvability. Two typical examples also show the local geometric characteristics of the parameter perturbation of these exact solutions.

Comment on “Hydraulic characteristics of labyrinth sluice gates” by T. Hashem, A.Y. Mohammed, T.J. Alfatlawi

In this paper, the stage-discharge relationship of a labyrinth sluice gate is theoretically deduced by the Π-Theorem of dimensional analysis and the self-similarity condition. This approach was then compared with the classical approach using the discharge coefficient Cd proposed by Hashem et al. (2024). The classical stage-discharge equation (Eq. 4), with the discharge coefficient estimated by the empirical relationship (Eq. (1)) proposed by Hashem et al. (2024), resulted in errors E that are less than or equal to ±5 % for 74.3 % of the investigated cases and a mean square error MSE of 0.0078. The new stage-discharge relationship, deduced by dimensional analysis and self-similarity theory, was calibrated using the dataset by Hashem et al. (2024). The obtained stage-discharge equation (Eq. 20) is characterized by errors in the discharge estimate E that are less than or equal to ±5 % for 82.2 % of the investigated cases and a mean square error MSE of 0.0051. In conclusion, the stage-discharge relationship deduced by dimensional analysis allows for the best performance in the discharge estimate, which can be obtained without knowledge of the empirical discharge coefficient.

Stochastic microbial dispersal drives local extinction and global diversity.

Airborne dispersal of microorganisms is a ubiquitous migration mechanism, allowing otherwise independent microbial habitats to interact via biomass exchange. Here, we study the ecological implications of such advective transport using a simple spatial model for bacteria-phage interactions: the population dynamics at each habitat are described by classical Lotka-Volterra equations; however, species populations are taken as integer, that is, a discrete, positive extinction threshold exists. Spatially, species can spread from habitat to habitat by stochastic airborne dispersal. In any given habitat, the spatial biomass exchange causes incessant population density oscillations, which, as a consequence, occasionally drive species to extinction. The balance between local extinction events and dispersal-induced migration allows species to persist globally, even though diversity would be depleted by competitive exclusion, locally. The disruptive effect of biomass dispersal thus acts to increase microbial diversity, allowing system-scale coexistence of multiple species that would not coexist locally.

A phase-splitting approach to describe macroscopically non-equilibrium transport in porous media

Classical macro-scale dispersion equation cannot correctly represent non-local, non-equilibrium effects, and alternatives had to be proposed to overcome these discrepancies. Two-equation models have been widely used and proved useful to incorporate, to some extent, local non-equilibrium effects. More accurate descriptions may be obtained through the introduction of additional complex terms in the one-equation dispersion equation, e.g. spatial and time convolutions. An alternative more practical route rely on the introduction of a N-equations description. For instance, instead of a two-temperature model, transfers in the solid and/or fluid phase maybe represented by several equations, each one defined by some splitting framework. Splitting of the solid phase has been proposed in the past (multi-rate models). Splitting of the flowing phase has a rustic illustration in the case of the Coats and Smith model, i.e., ”flowing” and ”stagnant” fluid equations. In this paper, we explore further the possibilities offered by N-equations splitting of the flowing phase. The N-equations macro-scale model is developed applying an averaging technique to N-pseudo phases obtained by splitting the flowing phase based on different criteria, e.g. the velocity histogram. Most important model properties are the matrices of dispersion tensors and exchange coefficients, which are provided by a system of specific ”closure” problems. To illustrate this approach, the methodology is applied to Taylor’s dispersion problem, in which non-local effects are produced by pore-scale spatially distributed inlet conditions. The proposed N-equations description was compared to estimates of the N-temperatures obtained from direct numerical simulations. The N-equations model was able to reproduce accurately the dynamic of pore-scale computations, for such reputedly non-homogenizable problems, with a relatively small number of equations.

Identification of the effective crane hook's cross-section by incorporating finite element method and programming language

The crane hook is a widely utilized component in several industries for the purpose of lifting things. The crane hook must possess the capacity to withstand the intended load without encountering any complications, hence ensuring the safety of both personnel and the objects being lifted. The process of analysis is crucial for the effective utilization of a crane hook. The primary aim of this study was to determine the most efficient cross-sectional crane hook among five distinct geometric profiles. This was achieved through the application of finite element analysis using Solidworks software. Subsequently, the identified cross-sectional profile was further examined using the Python programming language, taking into account the classical equation of a curved beam. The five cross-sectional shapes seen in the study were circular, rectangular, trapezoidal, I-shaped, and T-shaped. For the purposes of this investigation, the chosen material for each cross-sectional crane hook model was 34CrMo4 steel. Despite the identical boundary constraints imposed on all the chosen cross-sectional crane hook profiles, it was observed that the trapezoidal cross-sectional crane hook exhibited superior performance compared to the others. The trapezoidal cross-sectional crane hook model exhibited a Von Mises stress of 203 MPa, with a corresponding factor of safety of 3.20. Further experimentation was conducted using Python to examine the trapezoidal profile. The results indicated that an increased level of parallelism in the inner side of the trapezoidal shape corresponded to a higher factor of safety. Hence, it is advisable to maintain the trapezoidal cross-sectional profile of the crane hook, with due consideration given to maximizing the length of the inner parallel side. The enhancement of design leads to a decrease in the likelihood of failure and the occurrence of undesirable accidents.

Lie Symmetries and the Invariant Solutions of the Fractional Black–Scholes Equation under Time-Dependent Parameters

In this paper, we consider the time-fractional Black–Scholes model with deterministic, time-varying coefficients. These time parametric constituents produce a model with greater flexibility that may capture empirical results from financial markets and their time-series datasets. We make use of transformations to reduce the underlying model to the classical heat transfer equation. We show that this transformation procedure is possible for a specific risk-free interest rate and volatility of stock function. Furthermore, we reverse these transformations and apply one-dimensional optimal subalgebras of the infinitesimal symmetry generators to establish invariant solutions.

Higher Time-Derivative Theories from Space–Time Interchanged Integrable Field Theories

We compare a relativistic and a nonrelativistic version of Ostrogradsky’s method for higher-time derivative theories extended to scalar field theories and consider as an alternative a multi-field variant. We apply the schemes to space–time rotated modified Korteweg–de Vries systems and, exploiting their integrability, to Hamiltonian systems built from space–time rotated inverse Legendre transformed higher-order charges of these systems. We derive the equal-time Poisson bracket structures of these theories, establish the integrability of the latter theories by means of the Painlevé test and construct exact analytical period benign solutions in terms of Jacobi elliptic functions to the classical equations of motion. The classical energies of these partially complex solutions are real when they respect a certain modified CPT-symmetry and complex when this symmetry is broken. The higher-order Cauchy and initial-boundary value problem are addressed analytically and numerically. Finally, we provide the explicit quantization of the simplest mKdV system, exhibiting the usual conundrum of having the choice between having to deal with either a theory that includes non-normalizable states or spectra that are unbounded from below. In our non-Hermitian system, the choice is dictated by the correct sign in the decay width.

Noise-induced transport in the Aubry-André-Harper model

We study quantum transport in a quasiperiodic Aubry-André-Harper (AAH) model induced by the coupling of the system to a Markovian heat bath. We find that coupling the heat bath locally does not affect transport in the delocalized and critical phases, while it induces logarithmic transport in the localized phase. Increasing the number of coupled sites at the central region introduces a transient diffusive regime, which crosses over to logarithmic transport in the localized phase and in the delocalized regime to ballistic transport. On the other hand, when the heat bath is coupled to equally spaced sites of the system, we observe a crossover from ballistic and logarithmic transport to diffusion in the delocalized and localized regimes, respectively. We propose a classical master equation, which in the localized phase, captures our numerical observations for both coupling configurations on a qualitative level and for some parameters, even on a quantitative level. Using the classical picture, we show that the crossover to diffusion occurs at a time that increases exponentially with the spacing between the coupled sites, and the resulting diffusion constant decreases exponentially with the spacing.