This paper aims to investigate the global asymptotic stability of following system of higher order difference equations with quadratic terms: x_{n+1}=A+B((y_{n})/(y_{n-m}²)),y_{n+1}=A+B((x_{n})/(x_{n-m}²)) where A and B are positive numbers and the initial values are positive numbers. We also study the rate of convergence and oscillation behaviour of the solutions of related system.

In this paper we investigate the global asymptotic stability of following system ofhigher order difference equations with quadratic terms:xn+1=A+Byn/yn−m^2, yn+1=A+Bxn/xn−m^2, where A and B are positive numbers and the initial values are positive numbers.We also study the boundedness, rate of convergence and oscillation behaviour of thesolutions of related system.

A new type of ordinary differential equation is introduced and discussed, namely, the time-dependent order ordinary differential equations. These equations can be solved via fractional calculus and are mapped into Volterra integral equations of second kind with singular integrable kernel. The solutions of the time-dependent order differential equations smoothly deforms solutions of the classical integer order ordinary differential equations into one-another, and can generate or remove singularities. An interesting symmetry of the solution in relation to the Riemann zeta function and Harmonic numbers was also proved.

We employ another approach to quantize electromagnetic fields in the coordinate space, instead of the mode (or Fourier) space, such that local features of photons can be efficiently, physically, and more intuitively described. To do this, coordinate-ladder operators are defined from mode-ladder operators via the unitary transformation of systems involved in arbitrary inhomogeneous dielectric media. Then, one can expand electromagnetic field operators through the coordinate-ladder operators weighted by non-orthogonal and spatially-localized bases, which are propagators of initial quantum electromagnetic (complex-valued) field operators. Here, we call them QEM-CV-propagators. However, there are no general closed form solutions available for them. This inspires us to develop a quantum finite-difference time-domain (Q-FDTD) scheme to numerically time evolve QEM-CV-propagators. In order to check the validity of the proposed Q-FDTD scheme, we perform computer simulations to observe the Hong-Ou-Mandel effect resulting from the destructive interference of two photons in a 50/50 quantum beam splitter.

This paper is based on a talk given at the 14-th International Workshop on Differential Geometry and Its Applications, hosted by the Petroleum Gas University from Ploiesti, between July 9-th and July 11-th, 2019. After presenting some historical facts, we will consider some geometry problems related to unification approaches. Jordan algebras and Lie algebras are the main non-associative structures. Attempts to unify non-associative algebras and associative algebras led to UJLA structures. Another algebraic structure which unifies non-associative algebras and associative algebras is the Yang-Baxter equation. We will review topics relared to the Yang-Baxter equation and Yang-Baxter systems, with the goal to unify constructions from Differential Geometry.

The current health crisis derived from COVID 19 has become a serious threat to humanity. The collateral effects related to economic aspects are affecting companies, workers and, of course, families in general. The problems of financial insolvency are frequently present in debtors who have contracted a debt with their creditor. Therefore, the objective of the study is focused on modeling a debt restructuring proposal, through three hypothetical scenarios in which the debt is reconvened when some documents have expired and others are yet to expire. In each of the scenarios, the mathematical model is designed and subsequently corroborated in an Excel spreadsheet. The predominant variables in the restructuring are the adjustment rates used to update the past due promissory notes and the discount rates of those promissory notes that have not yet expired, capitalizations, transaction costs and, if applicable, whether the debt is in foreign currency. In this case, it is converted to the local currency, then the spot exchange rate is used.

It is universally accepted that Maxwell equations do not remain invariant under the Galilean transformation. This conflicts the principle of relativity which states that the physical law must remain invariant in the mathematical form in all inertial frames of reference. For this reason, the Lorentz transformation is invented, and the Galilean transformation is nowadays superseded by the Lorentz transformation. However, this paper challenges this widely held belief that the Maxwell equations are not invariant under the Galilean transformation. By applying the Galilean transformation to Lienard-Wiechert electromagnetic fields it is mathematically proven that the Maxwell equations indeed remain invariant under the Galilean transformation. In addition, the critical error in Lorentz's proof of Galilean non-invariance of Maxwell equations is pointed out, and it turns out that Lorentz's conclusion that Maxwell equations are not invariant under the Galilean transformation is the result of a mathematical error.

The Hubble parameter H(z), as a function of redshift, is modified by the presence of a new term originated from the extrinsic curvature in an embedded space-time. Assuming an asymptotic expansion factor $a\sim 0$ or, equivalently $z\rightarrow \infty$, it is obtained a nearly resemblance of the present model ($\beta$-model) with quintessence $w$CDM model with $H(z)\sim H_{wCDM}(z)$ at background level. In conjunction with $\Lambda$CDM, we test the models using a pack of recent datasets like that of the ``Gold 2018'' growth data, the best-fit Planck2018/$\Lambda$CDM parameters on the Cosmic Microwave Background (CMB), the Baryon Acoustic Oscillations (BAO) measurements, the Pantheon Supernovae type Ia and the Hubble parameter data with redshift ranging from $0.01 < z < 2.3$. Performing the Akaike Information Criterion (AIC) to ascertain the statistical viability of the model from Jeffreys' scale, we apply a joint likelihood analysis to the data with the Markov Chain Monte Carlo (MCMC) method. We find that the present model is in very good agreement with observations with a close statistical equivalence with the $\Lambda$CDM and $w$CDM cosmologies at 1-$\sigma$ level. We also show that a mild alleviation of the $\sigma$ tension between the growth amplitude factor and the matter content $(\sigma_8$-$\Omega_m)$ of the observations from CMB and Large Scale Structure (LSS) probes. A comparison of the aforementioned full pack of data in MCMC chains is made with the resulting MCMC from Pantheon SNIa+$H(z)$ in order to analyse the sensitivity of the models and how they respond to a cosmography analysis on the evolution of $H(z)$ and the deceleration parameter $q(z)$. In this sense, we find that only the $\beta$-model can be unaffected to the variations of the previous datasets due to its several minima in the likelihood (degeneracies) of the related parameters with an overall percentage relative difference only up to 4$\%$.

Deep learning machines are computational models composed of multiple processing layers of adaptive weights to learn representations of data with multiple levels of abstraction. Their structure is mainly reflecting the intuitive plausibility of decomposing a problem into multiple levels of computation and representation since it is believed that higher layers of representation allow a system to learn complex functions. Surprisingly, after decades of research, from learning and design perspectives these models are still deployed in a heuristic manner. In this paper, deep learning feed-forward machines are modeled from a statistical mechanics point of view as disordered physical systems where its macroscopic behavior is determined in terms of the interactions defined between the basic constituent of these models, namely, the artificial neuron. They are viewed as the equilibrium states of a theoretical body that is subject to the law of increase of the entropy. The study of the changes in energy of the body when passing from one equilibrium state to another is used to understand the structure and role of the phase space of the system, the stability of the equilibrium states, and the resulting degree of disorder. It is shown that the topology of these models is strongly linked to their stability and resulting level of disorder. Furthermore, the proposed theoretical characterization permit to assess the thermodynamic efficiency with which information can be processed by these models, and to provide a practical methodology to quantitatively estimate and compare their expected learning and generalization capabilities. These theoretical results provides new insights to the theory of deep learning and their implications are shown to be consistent through a set of benchmarks designed to experimentally assess their validity.

Using a joint statistical analysis, we test a four-dimensional FLWR model embedded in a five-dimensional bulk based on the Nash-Greene embedding theorem. Performing a Markov Chain Monte Carlo (MCMC) modelling, we combine observational data sets as those of the recent ``Gold 2018'' growth data, the best-fit Planck2018/$\Lambda$CDM parameters on the Cosmic Microwave Background (CMB), the Baryon Acoustic Oscillations (BAO) measurements, the Pantheon Supernovae type Ia and the Hubble parameter data with redshift ranging from $0.01 < z < 2.3$ to impose restrictions on the model. From linear Nash-Green fluctuations of the metric, we obtain the related perturbed equations in longitudinal Newtonian gauge to obtain the evolution of growth matter. We show that a mild alleviation is obtained from the $\sigma$ tension between the growth amplitude factor and the matter content $(\sigma_8$-$\Omega_m)$ of the observations from CMB and Large Scale Structure (LSS) probes with degeneracies on the parameters. We also apply the Om(z) diagnosis to distinguish the model from $\Lambda$CDM. The Akaike Information Criterion (AIC) is used and we find a relative statistical consistence of the present model with both $\Lambda$CDM and $w$CDM models up to 1$\%$ of percentage difference at early times.

We enhance Frobenius' method for solving linear ordinary differential equations about regular singular points. Key to Frobenius' approach is the exploration of the derivative with respect to a single parameter; this parameter is introduced through the powers of generalized power series. Extending this approach, we discover that tandem recurrence relations can be derived. These relations render coefficients for series occurring in logarithmic solutions. The method applies to the, practically important, exceptional cases in which the roots of the indicial equation are equal, or differ by a non-zero integer. We demonstrate the method on Bessel's equation and derive previously unknown tandem recurrence relations for coefficients of solutions of the second kind, for Bessel equations of all integer and half-integer order.

A transform approach based on a variable initial time (VIT) formulation is developed for discrete-time signals and linear time-varying discrete-time systems or digital filters. The VIT transform is a formal power series in z^(-1) which converts functions given by linear time-varying difference equations into left polynomial fractions with variable coefficients, and with initial conditions incorporated into the framework. It is shown that the transform satisfies a number of properties that are analogous to those of the ordinary z-transform, and that it is possible to do scaling of z^(- i) by time functions which results in left-fraction forms for the transform of a large class of functions including sinusoids with general time-varying amplitudes and frequencies. Using the extended right Euclidean algorithm in a skew polynomial ring with time-varying coefficients, it is shown that a sum of left polynomial fractions can be written as a single fraction, which results in linear time-varying recursions for the inverse transform of the combined fraction. The extraction of a first-order term from a given polynomial fraction is carried out in terms of the evaluation of z^(i) at time functions. In the application to linear time-varying systems, it is proved that the VIT transform of the system output is equal to the product of the VIT transform of the input and the VIT transform of the unit-pulse response function. For systems given by a time-varying moving average or an autoregressive model, the transform framework is used to determine the steady-state output response resulting from various signal inputs such as the step and cosine functions.

Formulating a relativistic equation for particles with arbitrary spin remains an open challenge in theoretical physics. In this study, the main algebraic approaches used to generalize the Dirac and Kemmer–Duffin equations for particles of arbitrary spin are investigated. It is proved that an irreducible relativistic equation formulated using spin matrices satisfying the commutation relations of the de Sitter group leads to inconsistent results, mainly as a consequence of violation of unitarity and the appearance of a mass spectrum that does not reflect the physical reality of elementary particles. However, the introduction of subsidiary conditions resolves the problem of unitarity and restores the physical meaning of the mass spectrum. The equations obtained by these approaches are solved and the physical nature of the solutions is discussed.

We review a modern differential geometric description of the fluid isotropic motion and featuring it the diffeomorphism group structure, modelling the related dynamics, as well as its compatibility with the quasi-stationary thermodynamical constraints. There is analyzed the adiabatic liquid dynamics, within which, following the general approach, there is explained in detail, the nature of the related Poissonian structure on the fluid motion phase space, as a semidirect Banach groups product, and a natural reduction of the canonical symplectic structure on its cotangent space to the classical Lie-Poisson bracket on the adjoint space to the corresponding semidirect Lie algebras product. We also present a modification of the Hamiltonian analysis in case of the isotermal liquid dynamics. We study the differential-geometric structure of the adiabatic magneto-hydrodynamic superfluid phase space and its related motion within the Hamiltonian analysis and invariant theory. In particular, we construct an infinite hierarchies of different kinds of integral magneto-hydrodynamic invariants, generalizing those before constructed in the literature, and analyze their differential-geometric origins. A charged liquid dynamics on the phase space invariant with respect to an abelian gauge group transformation is also investigated, some generalization of the canonical Lie-Poisson type bracket is obtained.

In this study, we investigate the formation of dynamical energy cascades in higher-order KdV-type equations. In the beginning, we recall what is known about the dynamic cascades for the classical KdV (quadratic) and mKdV (cubic) equations. Then, we investigate further the mKdV case by considering a richer set of initial perturbations in order to check the validity and persistence of various facts previously established for the narrow-banded perturbations. Afterwards, we focus on higher-order nonlinearities (quartic and quintic) which are found to be quite different in many respects from the mKdV equation. Throughout this study, we consider both the direct and double energy cascades. It was found that the dynamic cascade is always formed, but its formation is not necessarily accompanied by the nonlinear stage of the modulational instability. Direct cascade structure remains invariant regardless of the size of the spectral domain. In contrast, the double cascade shape can depend on the size of the spectral domain, even if the total number of cascading modes remains invariant. Results obtained in this study can be potentially applied to plasmas, free surface and internal wave hydrodynamics.

Queretaro is currently one of the states with the highest growth rate in the automotive industry. In the last five years, the number of companies has increased, reaching more than 300 small and medium-sized companies in the sector. However, they show a high degree of ignorance and lack of strategic foresight, which translates into poor planning in the short, medium and long term and low competitiveness, leading many of them to failure. This paper presents the results of a study within a strategic foresight evaluation with the development of a Structural Equation Model that allows the analysis of foresight and strategic planning within the automotive SMEs in the state of Queretaro. The study analyzes the necessary indicators to be evaluated and establishes the relationship of dependence between the variables, which is necessary to create the constructs of the model. It is confirmed that the adjustment of the model used is adequate for the evaluation of SMEs. The contributions of the research were: A theoretical contribution related to strategic foresight within SMEs and the construction of the structural equation model to evaluate strategic foresight in automotive SMEs in the state of Queretaro.

Integrating large-scale wind energy in modern power systems is demanding more efficient mathematical models to properly address classical assumptions in power system problems. In particular, there are two main assumptions in power system problems with wind integration that have not been adequately studied yet; First, non-linear AC power flow equations have been linearized in most of the literature. Such simplifications can lead to inaccurate power flow calculations that may result in other technical issues. Second, wind power uncertainties are inevitable and have been mostly modelled using the traditional uncertainty modelling approaches, that may not be suitable for large-scale wind power integration. In this paper, we address both challenges: we present a tight second-order conic relaxation (SOCR) for optimal power flow (OPF) problem, and simultaneously, implement the new effective budget of uncertainty approach for uncertainty modelling that determines the maximum wind power admissibility first and then addresses the uncertainty in the model. To the best of our knowledge, this is the first study that proposes an effective robust second-order conic programming (ERSOCP) model that simultaneously addresses the issues of power flow linearization and wind power uncertainty with the new paradigm on the budget of uncertainty approach. Our numerical results show the merit of the proposed model against traditional linearized power flow equations as well as traditional uncertainty modelling approaches.

One of the solutions of the Einstein equations, called McVittie solution, signifying a black hole embedded by the dynamic spacetime is studied. In the stationary spacetime the Mcvittie metric becomes the Schwarzschild-de Sitter metric (SdS). The geodesic of a freely falling test particle towards the black hole is examined in the SdS spacetime. It is found that unlike Schwarzschild case the potential of such particle becomes maximum at a point where it eventually stops to follow an unstable circular motion and then resumes its motion towards black hole center. It is shown that an observer or system of particles is spaghettified near the black hole singularity in the SdS spacetime. The dynamic of the universe in the framework of McVittie metric, being a generalized time dependent SdS solution, is represented in terms of that point, called stationary or turning point. The motion of the stationary point is studied in various regimes of the expanding universe and the possible outcomes are discussed in brief.

The constraint equations in Maxwell theory are investigated. In analogy with some recent results on the constraints of general relativity it is shown, regardless of the signature and dimension of the ambient space, that the "divergence of a vector field" type constraints can always be put into linear first order hyperbolic form for which global existence and uniqueness of solutions to an initial-boundary value problem is guaranteed.

Partial differential equations (PDEs) are used to describe a wide range of phenomena, such as heat transfer, fluid dynamics, and quantum mechanics. By solving PDEs, we can ob- tain insights into the behavior of the system and make predictions about its future evolution. Conventional numerical methods for obtaining the approximate solutions of PDEs may re- quire extensive computational resources and time, especially for complex PDEs and large- scale problems. The recently developed physics-informed neural network (PINN) has shown promise in a variety of scientific and engineering fields by incorporating physical laws into the loss functions of the neural network (NN). In addition, the training of PINN does not require ground truth data, but it has poor generalization ability to unseen domains. On the other hand, a convolutional neural network (CNN) has fast inference and better generalization ability, but it requires a large amount of training data. Taking the advantages of PINN and CNN by using Legendre multiwavelets (LMWs) as basis functions, we introduce a new method to approach the PDEs in this paper, namely Physics-Informed Legendre Multiwavelets CNN (PiLMWs-CNN), in order to continuously approximate a grid-based state representation that can be handled by a CNN. PiLMWs-CNN enable us to train our models using only physics-informed loss functions without any pre- computed training data, simultaneously providing fast and continuous solutions that gener- alize to previously unknown domains. In particular, the LMWs can simultaneously possess compact support, orthogonality, symmetry, high smoothness, and high approximation order. Compared to orthonormal polynomial (OP) bases, the approximation accuracy can be greatly increased and computation costs can be significantly reduced by using LMWs. We applied PiLMWs-CNN to approximate the damped wave equation, incompressible Navier-Stokes (N- S) equation, and two-dimensional heat conduction equation. The experimental results show that this method provides more accurate, efficient, and fast convergence with better stability when approximating the solution of PDEs.

We consider the Dirac electron as a nonperturbative particle-like solution consistent with its own Kerr-Newman (KN) gravitational and electromagnetic field. We develop the earlier models of the KN electron regularized by Israel and López, and consider the non-perturbative electron model as a bag model formed by Higgs mechanism of symmetry breaking. The López regularization determines the unique shape of the electron in the form of a thin disk with a Compton radius reduced by 4π. In our model this disk is coupled with a closed circular string which is placed on the border of the disk and creates the caused by gravitation frame-dragging string tension produced by the vector potential of the Wilson loop. Using remarkable features of the Kerr-Schild coordinate system, which linearizes the Dirac equation, we obtain solutions of the Dirac equation consistent with the KN gravitational and electromagnetic field, and show that this solution takes the form of a massless relativistic string. Parallelism of this model with quantum representations in Heisenberg and Schrodinger pictures explains remarkable properties of the stringy electron model in the relativistic scattering processes.

In this work, linear confinement of a relativistic scalar particle under the effects of Lorentz symmetry violation is investigated. We introduce a scalar potential by modifying the mass via transformation M → M + S(r) in the wave equation, and analyze the effects on the eigenvalues and the wave function. We see that the solution of the bound state to the wave equation can be achieved, and the energy eigenvalues and the wave function modified by the Lorentz symmetry breaking parameters as well as potential

A theoretical model to translate the evolution over time, in early stages, of growth and accumulation of biofilm bacterial mass is introduced. The model implies the solution of a system of differential-difference master equations. The application of an algorithm like Miller´s tree term recurrence, already known for Bessel functions of first kind, allows an exact calculation of the solutions of such equations, for a wide range of parameters values and time. For biofilm model a five term recurrence is deduced and applied in a backwards computation. A suitable normalisation condition completes the reach of the solution.

Considerable efforts have been devoted to modifying gravity, which aim to elucidate the possible existence or nature of dark matter and dark energy, achieve a better description of observation data, verify theoretical restrictions in the strong curvature regime and to formulate quantum gravity. In addition, despite the enormous success of the quantum field theory, the framework requires the so-called renormalization techniques and breaks down at high energies. Recently, the Planck Legacy 2018 release has confirmed the presence of an enhanced lensing amplitude in the cosmic microwave background power spectra, which prefers a positively curved early Universe with a confidence level greater than 99%. This study considers the implied positive curvature of the early Universe as the curvature of ‘the background or the 4D conformal bulk’ and distinguishes it from the localized curvature that is induced into the bulk by the presence of comparably smaller celestial objects that are regarded as ‘4D relativistic cloud-worlds’. To consider the interaction between the bulk and cloud-worlds, this paper presents interaction field equations in terms of the brane-world modified gravity counting for the bulk conformal curvature and the boundary terms, which could remove the singularities and satisfy a conformal invariance theory. Similarly, the quantum fields and clouds are regarded as ‘4D relativistic branes’ embedded in vacuum energy of a field strength reliant on the background curvature due to gravity. A visualization of the evolution of the 4D relativistic cloud-worlds over the conformal space-time of the 4D bulk is presented.

The mechanisms of natural oscillations and resonance are described, considering the peculiarities of the transformation of elastic and kinetic energy in the implementation of the law of conservation of energy in local and integral volumes of the body, using the concept of mechanics based on the concepts of space, time and energy. When describing the motion in the Lagrange form, the elastic deformation energy of the particles is determined by the quadratic invariant of the tensor, whose components are the partial derivatives of Euler variables with respect to Lagrange variables. The increment of the invariant due to elastic deformation is represented as the sum of two scalars, one of which depends on the average value of the relative lengths of the edges of the particles in the form of an infinitesimal parallelepiped, the second is equal to the standard deviation of these lengths from the average value. It is shown that each of the scalars can be represented in the form of two dimensionless kinematic parameters of elastic energy, which participate in different ways in the implementation of the law of conservation of energy. One part of the elastic energy passes into kinetic energy and participates in the implementation of the law of conservation of energy for the body as a whole, considering external forces. The second part is not converted into kinetic energy but changes the deformed state of the particles in accordance with the equations of motion while maintaining the same level of the part of the elastic energy of the particles used for this. The kinematic parameters differ from the volume density of the corresponding types of energy by a factor equal to the elastic modulus, which is directly proportional to the density and heat capacity of the material and inversely proportional to the volume compression coefficient. Transverse, torsional, and longitudinal vibrations are considered free and under resonance conditions. The mechanisms of transformation of forced vibrations into their own after the termination of external influences and resonance at the superposition of free and forced vibrations with the same or similar frequency are considered. The formation of a new free wave at each cycle with an increase in the amplitude, which occurs mainly due to internal energy sources, and not external forces, is justified.

In this paper, the sequential prime numbers are used as variables for the galactic spiral equations which were developed from the ROTASE model. Special spiral patterns are produced when prime numbers are treated with the unit of radian. The special spiral patterns produced with the first 1000 prime numbers have 20 spirals arranged in two groups. The two groups have perfect central symmetry with each other and are separated with two spiral gaps. The special spiral pattern produced with natural numbers from 1 to 7919 shows 6 spirals in the central area and 44 spirals in the outer area. The whole 7919 spiral points can be plotted with either 6-spiral pattern or 44-spiral pattern. For the spirals only produced by the prime numbers in the 6-spiral pattern plotting, the spiral 2 and spiral 3 each has only one spiral point produced by prime number 2 and 3, respectively, all other spiral points produced by other prime numbers are located on the spiral 1 and spiral 5. The special spiral pattern is well explained with careful analysis, it is concluded that all prime numbers greater than 3 must meet one of the equations: P₁ = 1 + 6 * n (n > 0) P₅ = 5 + 6 * n (n ≥ 0) In other words, every prime number greater than 3 is either a P₁ prime number or a P₅ prime number, no exception. Matching one of the equations is a necessary condition for a number to be a prime number, not a sufficient condition. Hope such sufficient condition can be found in the future. The number of P₁ prime numbers roughly equal the number of P₅ prime number in the first 2 billion prime numbers. The galactic spiral equations with golden angle can duplicate Vogel’s result for the simulation of sunflower seed head pattern, and a pinwheel pattern can be produced also with galactic spiral equations and 1 degree more than golden angle.

A new generation of universal tools and languages for modeling and simulation multi-physical domain applications emerged and became widely accepted, which generate large-scale systems of differential algebraic equations (DAEs) automatically. Motivated by the characteristics of DAEs systems with large dimension, high index or block structures, we first propose a modified Pantelides’ algorithm (MPA) for any high order DAEs based on its Σ matrix, which is similar to Pryce’s Σ method. By introducing a vital parameter vector, a modified Pantelides’ algorithm with parameter has been presented.It leads to a block Pantelides’ algorithm (BPA) naturally which can immediately compute the crucial canonical offsets for whole (coupled) systems with block-triangular form. We illustrate these algorithms by some examples. And numerical experiments show that the time complexity of BPA can be reduced by at least O(ℓ) compared to the MPA, which is mainly consistent with the results of our analysis.

In the present article we describe a few simple and efficient finite volume type schemes on moving grids in one spatial dimension for the sake of simplicity. The underlying finite volume scheme is conservative and it is accurate up to the second order in space. The main novelty consists in the motion of the grid. This new dynamic aspect can be used to resolve better the areas with high solution gradients or any other special features. No interpolation procedure is employed, thus an unnecessary solution smearing is avoided. Thus, our method enjoys excellent conservation properties. The resulting grid is completely redistributed according the choice of the so-called monitor function. Several more or less universal choices of the monitor function are provided. Finally, the performance of the proposed algorithm is illustrated on several examples stemming from the simple linear advection to the simulation of complex shallow water waves. We believe that the techniques described in our paper can be beneficially used to model tsunami wave propagation and run-up.

In this paper, the notion of θ∗-weak contraction is introduced, which is utilized to prove some fixed point results. These results are helpful to give a positive response to certain open question raised by Kannan [Amer. Math. Monthly 76:1969] and Rhoades [Contemp. Math. 72:1988] on the existence of contractive definition which does not force the mapping to be continuous at the fixed point. Some illustrative examples are also given to support our results. As applications of our result, we investigate the existence and uniqueness of a solution of non-linear matrix equations and integral equations of Volterra type as well.

In this work, we studied the Green’s functions of the second order differential operators with involution. Uniform equiconvergence of spectral expansions related to the second-order differential operators with involution is obtained. Basicity of eigenfunctions of the second-order differential operator operator with complex-valued coefficient is established.

In this study, we first present a time-fractional L$\hat{e}$vy diffusion equation of the exponential option pricing models of European option pricing and the risk-neutral parameter. Then, we modify a particular L$\hat{e}$vy-time fractional diffusion equation of European-style options. Introduce a more general model from the models based on the L$\hat{e}$vy-time fractional diffusion equation and review some recent findings regarding of the Europe option pricing of risk-neutral free.

In this paper, we prove some common fixed point theorems on complex partial metric space. The presented results gener- alize and expand some of the literature well-known results. We also explore some of the application of our key results.

In this paper defines the consecutive sum of the digits of a natural number, so far as it becomes less than ten, as an arithmetic function called and then introduces some important properties of this function by proving a few theorems in a way that they can be used as a powerful tool in many cases. As an instance, by introducing a test called test, it has been shown that we are able to examine many algebraic equalities in the form of in which and are arithmetic functions and to easily study many of the algebraic and diophantine equations in the domain of whole numbers. The importance of test for algebraic equalities can be considered equivalent to dimensional equation in physics relations and formulas. Additionally, this arithmetic function can also be useful in factorizing the composite odd numbers.

A multiscale mathematical model is proposed seeking to study the propagation dynamics of the Human Immunodeficiency Virus (HIV) in a group of young people between 15 and 24 years of age, through sexual contact without protection, considering the use of antiretroviral therapy (ART) and therapeutic failure. The model consists in a scale-free complex network that follows a power law, coupled with the immunological dynamics of each individual, that is, it considers the infection by the virus in the immune system of each HIV carrier, through a system of non-linear differential equations that govern the infection’s behavior in the immune system. Propagation of the virus in the network is modelled by taking into account information from the immunological status of each person. The study found that for a population to have high HIV prevalence, it is not necessary at the beginning of the simulation time for the virus to propagate rapidly. In addition, the study proves that with a higher number of sexual partners, there will be greater prevalence of HIV in the population and that the use of ART helps to control the propagation of the infection in the population. As an interesting result, it was also found that there is a higher number of HIV carriers who abandon ART than those who have access to it.

In the present paper, the effect of the non-linear thermal radiation on the neutral gas mixture in the unsteady state is investigated for the first time. The unsteady BGK technique of the Boltzmann kinetic equations for a neutral non-homogenous gas is solved. The solution of the unsteady case makes the problem more general significance than the stationary one. For this purpose, the moments' method, together with the traveling wave method, is applied. The temperature and concentration are calculated for each gas component and mixture for the first time.Furthermore, the study is held for aboard range of temperatures ratio parameter and a wide range of the molar fraction. The distribution functions are calculated for each gas component and the gas mixture. The significant non-equilibrium irreversible thermodynamic characteristics the entire system is acquired analytically. That technic allows us to investigate the consistency of Boltzmann's H-theorem, Le Chatelier principle, and thermodynamics laws. Moreover, the ratios among the different participation of the internal energy alteration are evaluated via the Gibbs formula of total energy. The final results are utilized to the argon-helium non-homogenous gas at different magnitudes of radiation force strength and molar fraction parameters. 3D-graphics are presented to predict the behavior of the calculated variables, and the obtained results are theoretically discussed.

In this work, we study a Klein-Gordon oscillator subject to Cornelltype potential in the background of the Lorentz symmetry violation determined by a tensor out of the Standard Model Extension. We introduce a Cornell-type potential S(r) = (η_L\,r + \frac{η_c}{r} ) by modifying the mass term via transformation $M → M + S(r)$ and then coupled oscillator with scalar particle by replacing the momentum operator $\vec{p}→ (\vec{p}+ i\,M\,ω\,\vec{r})$ in the relativistic wave equation. We see that the analytical solution to the Klein-Gordon oscillator equation can be achieved, and a quantum effect characterized by the dependence of the angular frequency of the oscillator on the quantum numbers of the relativistic system is observed

Recent studies have indicated that terahertz time-domain spectroscopy (THz-TDS) can stably and efficiently acquire output spectra using an affordable and compact multimode laser diode (MMLD) with delayed optical feedback as the light source. This research focused on a numerical analysis of the optimal conditions for employing an MMLD with delayed optical feedback (chaotic oscillating laser diode) in THz-TDS, utilizing multimode rate equations. The findings revealed that the intermittent chaotic output generated by the MMLD, characterized by concurrent picosecond pulse oscillations lasting several tens of picoseconds, proves to be highly effective for THz-TDS. By appropriately setting the amounts of injection current and optical feedback, and the delay time of optical feedback, intermittent chaotic oscillation can be attained within a considerably broad parameter range. Moreover, both the MMLD output spectrum and the THz-TDS output spectrum exhibit a consistently stable shape at the microsecond scale, demonstrating the attractor properties inherent in an MMLD with delayed optical feedback.

It is shown that the 5D geodetic equations and 5D Ricci identities give us a way to create a new viewpoint on some problems of Modern Physics, Astrophysics, and Cosmology. Specifically, the application of the 5D geodetic equations in (4+1) and (3+1+1) splintered forms obtained with the help of the monad and dyad methods made it possible to introduce a new, effective generalized concept of the rest mass of the elementary particle. The latter leads one to novel connections between the General Relativity and quantum field theories, as well as all of that, including the (4+1) splitting of the 5D Ricci identities, brings about a better understanding of the magnetic monopole problem and the vital difference in the origins of the Maxwell equations and gives rise to surprising connections between them. The obtained results also provide new insight into the mechanism of the 4D Universe’s expansion and its following acceleration.

In this paper, we introduced a generalized ∆-implicit locally contractive condition and give some examples to support it and to show its significance in fixed point theory. We prove that the mappings satisfying generalized ∆-implicit locally contractive condition admits a common fixed point, where, the ordered multiplicative GM−metric space is chosen as underlying space. The obtained fixed point theorems generalize many earlier fixed point theorems on implicit locally contractive mappings. In addition, some nontrivial and interesting examples are provided to support our findings. To demonstrate the originality of our new main result, we apply it to show existence of solutions to a system of nonlinear -Volterra type- integral equations.

In this paper it will be considered three deterministic models for the study of the interaction between the human immune system and a virus: the logistic model, the Gompertz model, and the generalized logistic model (or Richards model). A qualitative analysis of this models based on dynamical systems theory will be performed, by studying the local behavior of the equilibria and obtaining local dynamics properties from the linear stability point of view. Also, we will compare this models in order to understanding which is more appropriate to models the interaction between the human immune system and a virus. It will be obtained some natural medical interpretations which are available for all three models and can be useful to the medical community.

This paper results from our investigation into novel means of electromagnetic propulsion. It requires the basis of our claims to be put on a sound theoretical footing regarding the purported momentum exchange with the electromagnetic field. One of these concerns is the huge discrepancy between the energy density of the Zeropoint and its purported manifestation as the Cosmological Constant. Here we state that it is manifestly wrong to introduce the zeropoint at zero order into the stress-energy tensor, because it is something which describes zero particle count. As a fluctuation, it belongs in a higher order Taylor expansion in frequency of the stress-energy tensor. Furthermore in the 3^rd order in the Einstein constant our procedure is some 9 orders of magnitude too small. We make up this difference by suggesting that vacuum energy is much higher still and that more degrees of freedom exist in physics beyond the Standard Model or that there is interaction energy between the modes.

The relativistic quantum dynamics of a spin-0 scalar particle under the effects of the violation of Lorentz symmetry in the presence of a non-electromagnetic potential is analyzed. The central potential induced by the Lorentz symmetry violation is a linear electric and constant magnetic field and, analyze the effects on the eigenvalues and the wave function. We see there is a dependence of the linear charge density on the quantum numbers of the system

In this paper, we investigate the behaviour of a relativistic quantum oscillator under the effects of Lorentz symmetry violation determined by a tensor (KF)µναβ out of the Standard Model Extension. We analyze the quantum system under a Coulomb-type radial electric field and a uniform magnetic induced by Lorentz symmetry breaking effects under a Cornell-type potential, and obtain the bound states solution by solving the Klein-Gordon oscillator. We see a quantum effect due to the dependence of the angular frequency of the oscillator on the quantum numbers of the system, and the energy eigenvalues and the wave-function of the oscillator field get modified by the Lorentz symmetry breaking parameters as well as due to the presence of Cornell-type potential.

In this work, quantum dynamics of a spin-0 particle under the effects of Lorentz symmetry violation in the presence of Coulombtype non-electromagnetic potential $(S(r) ∝ \frac{1}{r})$ is investigated. The non-electromagnetic (or scalar) potential is introduced by modifying the mass term via transformation $M → M + \frac{η_c}{r}$ in the relativistic wave equation. The linear central potential induced by the Lorentz symmetry violation is a linear radial electric and constant magnetic field and, analyze the effects on the spectrum of energy and the wave function

In this work, the problem is considered for eliminating a non-uniform rectilinear smearing of an image by mathematical and computer-based way, for example, a picture of several cars taken with a fixed camera and moving with different speeds. The problem is described by a set of 1-dimensional Fredholm integral equations (IEs) of the first kind of convolution type with a 1-dimensional point spread function (PSF) at uniform smearing and by a set of new 1-dimensional IEs of a general type (not convolution type) with a 2-dimensional PSF at non-uniform smearing. The problem is also described by one 2-dimensional IE of convolution type with a 2-dimensional PSF at uniform smearing and by a new 2-dimensional IE of a general type with a 4-dimensional PSF at non-uniform smearing. The problem for solving Fredholm IE of the first kind is ill-posed (unstable). Therefore, IEs of convolution type are solved by the Fourier transform (FT) method and Tikhonov's regularization (TR), and IEs of general type are solved by the quadrature/cubature and TR methods. Moreover, the magnitude of the image smear Δ is determined by the original “spectral method”, which increases the accuracy of image restoration. It is shown that the use of a set of 1-dimensional IEs is preferable to one 2-dimensional IE in the case of non-uniform smearing. In the inverse problem (image restoration), the Gibbs effect (the effect of false waves) in the image may occur. It can be edge or inner. The edge effect is well suppressed by the proposed technique “diffusing the edges". In the case of an inner effect, it is eliminated with difficulty, but the image smearing itself plays the role of diffusing and suppresses the inner Gibbs effect to a large extent. It is shown (in the presence of impulse noise in an image) that the well-known Tukey median filter can distort the image itself, and the Gonzalez adaptive filter also distorts the image (but to a lesser extent). We propose a modified adaptive filter. A software package was developed in MatLab and illustrative calculations are performed.

This article deals primarily with the problem of determining the cutting force when machining hardened steels. Secondary issues are focused on the evaluation of surface quality on machined samples and the recommendation of cutting conditions. A wide variety of components are used in engineering, the final heat treatment of which is hardening. These components are usually critical in a particular product. The quality of these components determines the correct functioning of the entire technical equipment and ultimately its service life. In our case, these are the core parts of thrust bearings, specifically the rolling elements. The subject of the experiment is machining these components in the hardened state with cubic boron nitride tools and continuous measurement of the cutting force using a dynamometer. The following evaluation assesses the surface quality by both touch and non-touch methods. A structural equation with appropriate constant and exponents was then constructed from the data obtained using the dynamometer.

We investigate a scalar particle under Lorentz symmetry breaking effects determined by a tensor out of the Standard Model Extension (SME) in the presence of a Cornell-type potential by modifying the mass term M → M +S in the KG-equation. The field configuration is such that a Coulomb-type radial electric field and a constant magnetic field can be induced by Lorentz symmetry violation, and analyze the behaviour of a scalar particle. One can see that the bound states solution to the KG-equation under the consider effects can be achieved, and a quantum effect characterized by the dependence of charge density distribution parameter on the quantum numbers of the system is observed.

The relativistic quantum motion of a scalar particle under the effects of violation of the Lorentz symmetry in the presence of a linear confining potential is investigated. We see that the solution of the bound state to the modified Klein-Gordon equation can be obtained and a quantum effect characterized by the dependence of the magnetic field on the quantum numbers of the system is observed

This study delves into the intricate dynamics of sustainable fisheries through the lens of mathematical modeling, specifically employing the Lotka-Volterra equations. The Lotka-Volterra equations, originally conceived to understand predator-prey interactions, offer a potent framework to model and predict fish population dynamics within aquatic ecosystems. The aim of this research is to demonstrate the applicability and efficacy of the Lotka-Volterra equations in understanding the delicate balance of fish populations and aiding sustainable fisheries management. Through this mathematical approach, we seek to provide a structured methodology and present results that showcase the potential of this model in predicting and managing fish populations. In our analysis, we adjust the parameters of the Lotka-Volterra equations to simulate various scenarios, enabling a comprehensive understanding of sustainable fishing practices. Our findings underscore the significance of mathematical models, particularly the Lotka-Volterra equations, in informing strategies to maintain fish populations in a sustainable and ecologically balanced manner, thereby ensuring the longevity of marine ecosystems and the fisheries thatrely upon them.

In this paper, a relativistic scalar particle under Lorentz symmetry breaking effects in the presence of a scalar potential is investigated. We introduce the scalar potential by modifying the mass via transformation M → M+S(r) in the wave equation and analyze the behaviour of a scalar particle. We see that the analytical solution to the KleinGordon equation can be achieved, and the energy eigenvalues and the wave function depends on the Lorentz symmetry breaking parameters as well as potential

The negative energy solutions of the Klein-Gordon and Dirac equations have been used to define the antiparticle with the positive energy from the Feynman - Stueckelberg interpretation in the standard model based on the 4-D Minkowski space with the relative time. The negative energy solutions of the Klein-Gordon and Dirac equations are used to define the antiparticles with the negative energy in the present 3-D quantized space model. Note that all the present results are based on the 4-D Euclidean space with the absolute time (t) and relative time (t_l). Then the particle with the positive energy and antiparticle with the negative energy have the 4-D P₄ = TCP = T_cP symmetry. The continuity equations are derived from the Klein-Gordon, Dirac and Schrodinger equations. The Klein-Gordon and Schrodinger equations with the potential barriers are solved. The Klein paradox is explained by using the pair production of the electron with the positive energy and positron with the negative energy under the P₄ = _TcP symmetry, The effective velocity ( c_eff=√2 " c " ) and effective energy ("Etl =" √2 E ) along the relative time axis (ctl) give the relative time momentum of "p_tl"=E_tl/c_eff =(√2 E)/(√2 c)=E/c=p_t=p_x. This indicates that the photons with the energy of E have the constant photon velocity of c and the constant photon momentum of "ptl"=E_tl/c_eff =(√2 E)/(√2 c)=E/c=p_t=p_x. For the negative energy particles, the space momentum directions are opposite to the particle velocity directions because the energy is negative. And the matter collapse to the photons at the very high particle velocity is discussed from the length expansion. It is, for the first time, proposed that the gamma ray burst is originated from the matter collapse to the photons near the black hole. Also, the wave function collapse is explained from the length expansion. The wave function collapse takes place when the measurement makes the particle velocity to be zero. It is concluded that the negative energy solutions of the Klein-Gordon and Dirac equations support the existence of the partner antimatter universe with the negative energy and negative time direction. Our matter universe has the positive energy and the positive time direction.

A thermodynamic description of nano-porous media must handle the size- and shape-dependence of the media properties. Such dependencies are typically due to the presence of immiscible phases, contact areas and contact lines. We propose a way to obtain average densities suitable for integration on the course grained scale, applying Hill's thermodynamics for small systems to the subsystems. we argue that the average densities of the porous medium, when defined in a proper way, obey the Gibbs equation. All contributions are additive or weakly coupled. From the Gibbs equation and the balance equations, we derive the entropy production in the standard way, for transport of multi-phase fluids in a non-deformable, porous medium exposed to di¤erences in boundary pressures, temperatures, and chemical potentials. Linear relations between thermodynamic fluxes and forces follow for the control volume. Fluctuation- dissipation theorems are formulated for the first time, for the fluctuating contributions to fluxes in the porous medium. These give an added possibility for determination of porous media permeabilities. Practical possibilities are further discussed.

This short paper demonstrates that the galactic spiral equations developed from ROTASE model can be used to morphologically simulate the phyllotactic spiral-ring patterns of plants. It is also noticed during this project that although the numbers in the Fibonacci sequence dominate the phyllotactic numbers, there are still significant phyllotactic numbers are not in the Fibonacci sequence, however, most of those numbers are well fit into a Fibonacci sibling sequence starting with numbers 0 and 2. Many plants show the phyllotactic number pairs of 10-16, 16-26 and 26-42 from the Fibonacci sibling sequence, which is the second dominated sequence in the phyllotaxis.

The (numerical) solution of Ordinary Differential Equations (ODEs) problems is of paramount relevance, ODEs system being an ubiquitous mathematical formulation of many physical phenomena (such as those involved in fluid dynamics, chemistry, biology, evolutionary-anthropology, ...): almost every dynamici phenomenon can be modeled by means of an ODEs system. The present paper is the first manifesto of FOODIE, a library aimed to numerically solve ODEs problems by means of a clear, concise and efficient abstract interface. FOODIE, meaning Fortran Object oriented Ordinary Differential Equations integration library, has manifolds aims: to provide a set to built-in numerical schemes that are accurate, robust, validated and efficient and to allow easy application of these schemes to (almost) all ODEs problems by means of an effective Abstract Calculus Pattern. The key idea is to allow the same solver-implementation to be applied to all ODEs problems thus avoiding the re-implementation of the ODEs solver for each different ODEs problem: code re-usability is consequently maximized, FOODIE being a general robust framework. Besides, the same framework also allows rapid development of new ODEs solvers due to the high abstraction level of the library itself. The present paper is the first announcement of FOODIE project: the current implementation is extensively discussed and its capabilities are proved by means of tests and examples.

In this paper, we consider the effects of a radial electric field and a constant magnetic field induced by Lorentz symmetry violation on a generalized relativistic quantum oscillator by choosing a function f(r) = b1 r + b2/r in the equation subject to a Cornell-type potential S(r) = ηL r + ηc/ r introduce by modifying the mass term in the equation. We show that the analytical solutions to the Klein-Gordon oscillator can be achieved, and a quantum effect is observed due to the dependence of the angular frequency of the oscillator on the quantum numbers of the system

In the last years, a recent growing line of research has been developed with fruitful results by using fractional and conformable derivatives in the iterative procedures of classical methods for solving nonlinear equations. In that sense, the use of conformable derivatives has shown better behavior than fractional ones, not only in the theory, but also in the practice. In this work, we adapt the approximation of conformable derivatives in order to design the first conformable derivative-free iterative schemes to solve nonlinear equations: a Steffensen’s type method and a Secant type method; the latter with memory. Convergence analysis is made, preserving the order of classical cases, and the numerical performance is studied in order to confirm the results in the theory. It is shown that these methods can present some numerical advantages versus their classical partners, with wide sets of converging initial estimations.

In this work, we offer a novel and accurate method in order to find the solution of the linear differential equations over the intervals [0, 1) based on the generalization of Legendre wavelets. The mechanism is still upon workable implementation of the operational matrix of integration and its derivatives. This method reduces the problems into algebraic equations via the properties of generalized Legendre wavelet (GLW) together with the operational matrix of integration. The function approximation has been picked out in such a way so as to enumerate the connection coefficients in an facile manner. The proposed numerical technique, based on the GLW, has been examined on three linear problems as a consequence of this investigation. The outcomes have shown that this method, as opposed to some other existing numerical and analytical methods, is a very useful and advantageous for tackling such problems.

We suggest progress regarding the following six physics opportunities. List all elementary particles. Describe dark matter. Explain ratios of dark matter to ordinary matter. Explain eras in the history of the universe. Link properties of objects. Interrelate physics models. We use models based on Diophantine equations.

Reference evapotranspiration (ET0) is the first step in calculating crop irrigation demand, and numerous methods have been proposed to estimate this parameter. FAO-56 Penman-Monteith (PM) is the only standard method for defining and calculating ET0. However, it requires radiation, air temperature, atmospheric humidity, and wind speed data, limiting its application in regions where these data are unavailable; therefore, new alternatives are required. This study compared the accuracy of ET0 calculated with the Blaney-Criddle (BC) and Hargreaves-Samani (HS) methods versus PM using information from an automated weather station (AWS) and the NASA-POWER platform (NP) for different periods of time. The information collected corresponding at Module XII of the Lagunera Region Irrigation District 017, a semi-arid region in the North of Mexico. The HS method underestimated by 5.5 % the reference evapotranspiration (ET0) compared to the PM method during the period from March to August, and yielded the best fit in the different evaluation periods: daily, average, and 5-day cumulative; the latter showed the best values of inferential parameters. The information about maximum and minimum temperatures from the NP platform was suitable for estimating ET0 using the HS equation. This data source is a timely alternative, particularly in semi-arid regions where no data from weather stations are available.

The correspondence between classical and quaternion quantum equations, permits considering the universe (vacuum) as an ideal elastic solid. Elementary particles would have to be standing or soliton-like waves. Tension induced by the compression and twisting of the elastic medium would increase the energy density, consequently generate a gravity forcing and affect the wave speed. Consequently the gravity could be described by an index of refraction. Theory was created by combining the Cauchy model of the elastic continuum with the Planck-Kleinert crystal hypothesis. The quaternion-imaginary Lagrange’an, the quaternion motion equation and the quaternionic oscillator allowed deriving: - The Schrödinger equation from the functional integral, which identifies the quaternion-imaginary quantum Hamiltonian. - The 2nd order wave equation system describing both the bosons and the gravity in terms of quaternionic Poisson equation. - The first order quaternionic wave equation system. - The family of the second order wave equation systems describing both the particles and the generated quaternionic force-fields (four-potentials). - The fundamental Planck and gravity constants. - The quaternionic continuity equation in an ideal elastic solid.

This paper presented a two-dimensional, well-balanced hydrodynamic and sediment transport model based on the solutions of variable density shallow water equations (VDSWEs) for sediment-laden flows, and the Exner equation for bed changes. Those equations are solved in a coupled way by the first-order Godunov-type finite volume method. The Harten-Lax-van Leer-Contact (HLLC) Riemann solver is extended to find the local Riemann fluxes in order to maintain the exact balance between the momentum term and the bed slope term. The advantage of a well-balanced scheme over an unbalanced scheme is demonstrated by the synthetic standing contact-discontinuity test case. Then, the model is employed to simulate two laboratory experiments. At last, a field case, the 1996 Lake Ha! Ha! flood event (Canada), is simulated. Results of cross sectional geometries and profiles of longitudinal thalweg agree well with measurements. The accuracy and simplicity of the numerical model, together with the robust implementation, make the model a good candidate for practical engineering applications.

We review a number of preconditioners for the advection-diffusion operator and for the Schur complement matrix which in turn constitute the building blocks for Constraint and Triangular Preconditioners to accelerate the iterative solution of the discretized and linearized Navier-Stokes equations. An intensive numerical testing is performed onto the driven cavity problem with low values of the viscosity coefficient. We devise an efficient multigrid preconditioner for the advection-diffusion matrix which, combined with the commuted BFBt Schur complement approximation, and inserted in a 2 x 2 block preconditioner, provides convergence of the GMRES method in a number of iteration independent of the meshsize for the lowest values of the viscosity parameter. The low-rank acceleration of such preconditioner is also investigated showing its great potential.

An overview is presented of several diverse branches of work in the area of effectively 2D fluid equilibria which have in common that they are constrained by an infinite number of conservation laws. Broad concepts, and the enormous variety of physical phenomena that can be explored, are highlighted. These span, roughly in order of increasing complexity, Euler flow, nonlinear Rossby waves, 3D axisymmetric flow, shallow water dynamics, and 2D magnetohydrodynamics. The classical field theories describing these systems bear some resemblance to perhaps more familiar fluctuating membrane and continuous spin models, but the fluid physics drives these models into unconventional regimes exhibiting large scale jet and eddy structures. From a dynamical point of view these structures are the end result of various conserved variable forward and inverse cascades. The resulting balance between large scale structure and small scale fluctuations is controlled by the competition between energy and entropy in the system free energy, in turn highly tunable through setting the values of the conserved integrals. Although the statistical mechanical description of such systems is fully self-consistent, with remarkable mathematical structure and diversity of solutions, great care must be taken because the underlying assumptions, especially ergodicity, can be violated or at minimum lead to exceedingly long equilibration times. Generalization of the theory to include weak driving and dissipation (e.g., non-equilibrium statistical mechanics and associated linear response formalism) could provide additional insights, but has yet to be properly explored.

The temporal evolution of pandemics described by the susceptible-infectious-recovered (SIR)-compartment model is sensitively determined by the time dependence of the infection (a(t)) and recovery (u(t)) rates regulating the transitions from the susceptible to the infected and from the infected to the recovered compartment, respectively. Here approximated SIR-solutions for different time dependencies of the infection and recovery rates are derived which are based on the adiabatic approximation assuming time-dependent ratios k(t)=u(t)/a(t) varying slowly in comparison to the typical time characteristics of the pandemic wave. For such slow variations the available analytical approximations from the KSSIR-model, valid for a stationary value of the ratio k, are used to insert a-posteriori the adopted time-dependent ratio of the two rates. Instead of investigating endless different combinations of the time dependencies of the two rates a(t) and u(t) a suitably parameterized reduced time dependence of the ratio k(tau) is adopted. Together with the definition of the reduced time this parameterized ratio k(tau) allows us to cover a great variety of different time dependencies of the infection and recovery rates. The agreement between the solutions from the adiabatic approximation in its four different studied variants and the exact numerical solutions of the SIR-equations is remarkably good providing strong confidence in the accuracy of the proposed adiabatic approximation.

The present paper introduces a method of basis transformation of vector fields that is specifically applicable to polynomials space and differential equations with certain polynomials solutions such as Hermite, Laguerre and Legendre polynomials. The method based on separated transformation of vector space basis by a set of operators that are equivalent to the formal basis transformation and connected to it by linear combination with projection operators. Applying the Forbenius covariants yields a general method that incorporates the Rodrigues formula as a special case in polynomial space. Using the Lie algebra modules, specifically , on polynomial space results in isomorphic algebras whose Cartan sub-algebras are Hermite, Laguerre and Legendre differential operators. Commutation relations of these algebras and Baker-Campbell-Hausdorff formula gives new formulas for special polynomials

Many authors proposed analytical methods for solving fully fuzzy Sylvester matrix equation (FFSME) based on Vec-operator and Kronecker product. However, these methods are restricted to nonnegative fuzzy numbers and cannot be extended to FFSME with near-zero fuzzy numbers. The main intention of this paper is to develop a new numerical method for solving FFSME with near-zero trapezoidal fuzzy numbers that provides a wider scope of trapezoidal fully fuzzy Sylvester matrix equation (TrFFSME) in scientific applications. This numerical method can solve the trapezoidal fully fuzzy Sylvester matrix equation with arbitrary coefficients and find all possible finite arbitrary solutions for the system. In order to obtain all possible fuzzy solutions, the TrFFSME is transferred to a system of non-linear equations based on newly developed arithmetic fuzzy multiplication between trapezoidal fuzzy numbers. The fuzzy solutions to the TrFFSME are obtained by developing a new two-stage algorithm. To illustrate the proposed method numerical example is solved.

L\'evy walks represent important modeling tools for variety of real life processes. Their natural scaling limits are known to be described by the so-called material fractional derivatives. So far these scaling limits were derived for spatially homogeneous walks, where Fourier and Laplace transforms represent natural tools of analysis. Here we derive the corresponding limiting equations in the case of position depending times and velocities of walks, where Fourier transforms cannot be effectively applied. In fact, we derive three different limits (specified by the way the process is stopped at an attempt to cross the boundary), leading to three different multi-dimensional versions of Caputo-Dzherbashian derivatives, which correspond to different boundary conditions for the generators of the related Feller semigroups and processes. Some other extensions and generalisations are analysed.

Abstract:Shown is the derivation of Lorentz-Einstein k-factor in SRT as an amplitude-term of oscillation-differential equations of second order.This case is shown for classical Lorentz-factor as solution of an equation for undamped oscillation as well as the developed theorem as a second solution for advanced SRT of fourth order with an equation for damped oscillation-states.This advanced term allows a calculation for any velocities by real rest mass

The fractional diffusion equations with integrals play a significant role in describing anomalous diffusion phenomena. We first construct an appropriate resolvent family, and use the convolution theorem of Laplace transform, the probability density function, Cauchy integral formula and Fubini theorem to study its related equicontinuity, strongly continuity, compactness, etc. Then, a reasonable mild solution is defined. Finally, by combining some fixed point theorems, we obtain some conclusions on the existence and uniqueness of mild solutions.

Cell migration is essential for physiological, pathological and biomedical processes such as, in embryogenesis, wound healing, immune response, cancer metastasis, tumour invasion and inflammation. In light of this, quantifying mechanical properties during the process of cell migration is of great interest in experimental sciences, yet few theoretical approaches in this direction have been studied. In this work, we propose a theoretical and computational approach based on the optimal control of geometric partial differential equations to estimate cell membrane forces associated with cell polarisation during migration. Specifically, cell membrane forces are inferred or estimated by fitting a mathematical model to a sequence of images, allowing us to capture dynamics of the cell migration. Our approach offers a robust and potentially accurate framework to compute geometric mechanical membrane forces associated with cell polarisation during migration and also yields geometric information of independent interest, we illustrate one such example that involves quantifying cell proliferation levels which are associated with cell division, cell fusion or cell death.

In this study, the hot deformation of Cu-0.55Sn-0.08La (wt%) alloy was studied by gleeble-3180 testing machine at deformation temperatures of 400-700oC and various strain rates. The stress-strain curve showed that the hot deformation behavior of Cu-0.55Sn-0.08La (wt%) alloy was greatly affected by work hardening, dynamic recovery, and dynamic recrystallization. The activation energy Q was calculated as 287.141 kJ/mol , and hot compression constitutive equation was determined as ε ̇=[sinh(0.006388σ)]^12.123∙exp⁡(37.7505-287.141/RT). The microstructural evolution of the alloy during deformation at 400oC revealed the presence of both slip and shear bands in the grains. At 700oC, dynamic recrystallization grains were observed but recrystallization was incomplete. In sum, these results provided the theoretical basis for continuous extrusion process of alloys with promising future use.

The splitting method developed in a previous study for linear, multi – component, convection-diffusion mass transfer equations is extended to non-linear, multi-component, convection-diffusion mass transfer equations. The diffusion term is non-linear (the diffusion coefficients depend on chemical species concentrations). Theoretical results about the splitting error are presented, as well as numerical simulations for mixtures with n ≤ 5. The numerical experiments show second-order accuracy and stability for the present splitting algorithm.

This work emphasizes some geometrical properties of the Maxwell--Bloch equations. Based on these properties the closed-form solutions of their equations are established. Thus, the Maxwell-Bloch equations are reduced to a nonlinear differential equation depending on an auxiliary unknown function. The approximate analytical solutions were built using the Optimal Homotopy Asymptotic Method (OHAM). A good agreement between the analytical and corresponding numerical results was found. The accuracy of the obtained results is validated through the representative figures. This procedure could be successfully applied for more dynamical systems with geometrical properties.

Unsupervised physics-informed deep learning can be used to solve computational physics problems by training neural networks to satisfy the underlying equations and boundary conditions without labeled data. Parameters such as network architecture and training method determine the training success. However, the best choice is unknown a priori as it is case specific. Here, we investigated network shapes, sizes, and types for unsupervised physics-informed deep learning of the two-dimensional Reynolds-averaged flow around cylinders. We trained mixed-variable networks and compared them to traditional models. Several network architectures with different shape factors and sizes were evaluated. The models were trained to solve the Reynolds-averaged Navier-Stokes equations incorporating Prandtl's mixing length turbulence model. No training data were deployed to train the models. The superiority of the mixed-variable approach was confirmed for the investigated high Reynolds number flow. The mixed-variable models were sensitive to the network shape. For the two cylinders, differently deep networks showed superior performance. The best fitting models were able to capture important flow phenomena such as stagnation regions, boundary layers, flow separation, and recirculation. We also encountered difficulties when predicting high Reynolds number flows without training data

We provide methodology and numerical results for the Hattendorff differential equa- tion for the continuous time evolution of the variance of L(j)t , the loss at time t random variable for a multi-state process, given that the state at time t is j.

This article presents a concise survey of basic discrete and semi-discrete nonlinear models which produce two- and three-dimensional (2D and 3D) solitons, and a summary of main theoretical and experimental results obtained for such solitons. The models are based on the discrete nonlinear Schrödinger (DNLS) equations and their generalizations, such as a system of discrete Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang corrections, the 2D Salerno model (SM), DNLS equations with long-range dipole-dipole and quadrupole-quadrupole interactions, a system of coupled discrete equations for the second-harmonic generation with the quadratic (χ⁽²⁾) nonlinearity, a 2D DNLS equation with a superlattice modulation opening mini-gaps, a discretized NLS equation with rotation, a DNLS coupler and its PT-symmetric version, a system of DNLS equations for the spin-orbit-coupled (SOC) binary Bose-Einstein condensate, and others. The article presents a review of basic species of multidimensional discrete modes, including fundamental and vortex solitons, their bound states, gap solitons populating mini-gaps, symmetric and asymmetric solitons in the conservative and PT-symmetric couplers, cuspons in the 2D SM, discrete SOC solitons of the semi-vortex and mixed-mode types, 3D discrete skyrmions, and some others.

In this paper we study solutions and drift homogenization for a class of viscous lake equations by using the method of semigroups of bounded operators. Suppose that the initial value i.e.,for some Hölder continuous function onwith smooth function value satisfying and Then the initial value problem (2) for viscous lake equations has a unique smooth local strong solution. Using this result we study the drift homogenization for three-dimensional stationary Stokes equation in the usual sense

In this paper, we report orthogonal fuzzy versions of some celebrated iterative mappings. We provide various concrete conditions on the real valued functions J,S:(0,1]→(−∞,∞) for the existence of fixed-points of (J,S)-fuzzy iterative mappings. We obtain many fixed point theorems in orthogonal fuzzy metric spaces. We apply (J,S)-fuzzy version of Banach fixed point theorem to show the existence and uniqueness of the solution. We support these results with several non-trivial examples and applications to Volterra-type integral equations and fractional differential equations.

Using valid experimental parameters we quantify the magnitude of resonantly phonon-driven precession of exchange magnons in freestanding ferromagnetic nickel thin films on their thickness L. Analytical solutions of acoustically driven equations for magnon oscillators display a nonmonotonous dependence of the peak magnetization precession on the film thickness. It is explained by different L-dependence of multiple prefactors entering in the expression for the total magnetization dynamics. Depending on the ratio of acoustic and magnetic (Gilbert) damping constants, the magnetization precession is shown to be amplified by a Q-factor of either the phonon or the magnon resonance. The increase of the phonon mode amplitude for thinner membranes is also found to be significant. Focusing on the magnetization dynamics excited by the two first acoustic eigenmodes with p=1 and p=2 we predict the optimum thicknesses of nickel membranes to achieve large amplitude magnetization precession at multi 100 GHz frequencies at reasonably low values of an external magnetic field.

In this study we investigate the formation of dynamical energy cascades in higher order KdV-type equations. In the beginning we recall what is known about the dynamic cascades for the classical KdV (quadratic) and mKdV (cubic) equations. Then, we investigate further the mKdV case by considering a richer set of initial perturbations in order to check the validity and persistence of various facts previously established for the narrow-banded perturbations. Afterwards we focus on higher order nonlinearities (quartic and quintic) which are found to be quite different in many respects from the mKdV equation. Throughout this study we consider both the direct and double energy cascades. It was found that the dynamic cascade is always formed, but its formation is not necessarily accompanied by the nonlinear stage of the modulational instability. Direct cascade structure remains invariant regardless the size of the spectral domain. In contrast, the double cascade shape can depend on the size of the spectral domain, even if the total number of cascading modes remains invariant. Results obtained in this study can be potentially applied to plasmas, free surface and internal wave hydrodynamics.

The standard approach to integrable nonlinear evolution equations (NLEE) usually uses the following steps: 1) Lax representation $[L,M]=0$; 2)~construction of fundamental analytic solutions (FAS); 3) reducing the inverse scattering problem (ISP) to a Riemann-Hilbert problem (RHP) $\xi^+(x,t,\lambda) =\xi^-(x,t,\lambda) G(x,t\lambda)$ on a contour $\Gamma$ with sewing function $G(x,t,\lambda)$; 4) soliton solutions and possible applications.Step 1) involves several assumptions: the choice of the Lie algebra $\mathfrak{g}$ underlying $L$, as well as its dependence on the spectral parameter, typically linear or quadratic in $\lambda$. In the present paper we propose another approach which substantially extends the classes of integrable NLEE.Its first advantage is that one can effectively use any polynomial dependence in both $L$ and $M$. We use the following steps: A) Start with canonically normalized RHP with predefined contour $\Gamma$; B) Specify the $x$ and $t$ dependence of the sewing function defined on $\Gamma$; C) Introduce convenient parametrization for the solutions $ \xi^\pm(x,t,\lambda)$ of the RHP and formulate the Lax pair and the nonlinear evolution equations (NLEE); D) use Zakharov-Shabat dressing method to derive their soliton solutions. This needs correctly taking into account the symmetries of the RHP. E) Define the resolvent of the Lax operator and use it analyze its spectral properties.

The main goal of the paper is to use a generalized proportional Riemann-Liouville fractional derivative (GPRLFD) to model BAM neural networks and to study some stability properties of the equilibrium. Initially, several properties of the GPRLFD are proved such as the fractional derivative of a squared function. Also some comparison results for GPRLFD are provided. Two types of equilibrium of the BAM model with GPRLFD are defined. In connection with the applied fractional derivative and its singularity at the initial time the Mittag-Leffler exponential stability in time of the equilibrium is introduced and studied. An example is given illustrating the meaning of the equilibrium as well as its stability properties.

The gauge field equation for fluid mechanics established in Part 1 is developed into a first-order scattering theory in the simplified case of a two-dimensional incompressible flow over a flat plate. This is used to present a model for the origin of Tollmien-Schlichting (TS) waves based on scattering between fluid particles. As a result, analytical formulae for the maximum amplification factor and the transition point from laminar to turbulent flow in the boundary layer are obtained. The mathematical transformations from the stationary field equations in Part 1 to a scattering theory with time evolution along the flow axis using Wick rotation are elaborated in detail.

This paper develops a methodology to control the navigation of an autonomous ship robot from an initial state to a final. To solve the problem the following approach is used: The ship robot system is modelled as a control system of six ordinary differential equations involving six state variables and three control variables. After having computed the system Hamiltonian, the feasible controls for optimality are derived from the normal equations of optimality as functions of the state and costate variables. Such controls are substituted into the state and the costate equations and give a combined control-free state-costate system of ordinary differential equations which is solved numerically by using a developed Matlab computer program. Associated computational simulations are designed and provided to show the reliability of the approach.

Large values and gradients of stresses and deformations, triggering concentrations of stresses and deformations, arise in the corner areas of a structure. The action of forced deformations, leading to the finite rupture of the contact between the elements of a structure, also triggers the concentration of stresses, while the rupture reaches an irregular point, a line on the area boundary. The theoretical analysis of the stress-strain state (SSS) of areas with angular cutouts in the boundary under the action of discontinuous forced deformations is reduced to the study of singular solutions to the homogeneous problem of the elasticity theory that has power-related features. The calculation of stress concentration coefficients in the domain of a singular solution to the elastic problem makes no sense. It is experimentally proven that the zone, that is close to the vertex of the angular cutout in the area boundary, has substantial deformations, rotations, and it corresponds to rising values of the first and second derivatives of displacements along the radius in cases of sufficiently small radii in the neighbourhood of the irregular point of the boundary. For such areas, it is necessary to consider the plane problem of the elasticity theory, taking into account the geometric nonlinearity under the action of forced deformations. This will allow analyzing the effect of relations between orders of values of deformations, rotations, and forced deformations on the form of the equation of equilibrium. The purpose of this work is to analyze the effect of relations of deformation orders, rotations, forced deformations on the form of the equilibrium equation in the polar coordinate system for a V-shaped area under the action of forced temperature-induced deformations with regard for the geometrical non-linearity and physical linearity.

We have found an infinite dimensional manifold of exact solutions of the Navier-Stokes loop equation for the Wilson loop in decaying Turbulence in arbitrary dimension $d >2$. This solution family is equivalent to a fractal curve in complex space $\mathbb C^d$ with random steps parametrized by $N$ Ising variables $\sigma_i=\pm 1$, in addition to a rational number $\frac{p}{q}$ and an integer winding number $r$, related by $\sum \sigma_i = q r$. This equivalence provides a \textbf{dual} theory describing a strong turbulent phase of the Navier-Stokes flow in $\mathbb R_d$ space as a random geometry in a different space, like ADS/CFT correspondence in gauge theory. From a mathematical point of view, this theory implements a stochastic solution of the unforced Navier-Stokes equations. For a theoretical physicist, this is a \textbf{quantum} statistical system with integer-valued parameters, satisfying some number theory constraints. Its long-range interaction leads to critical phenomena when its size $N \rightarrow \infty$ or its chemical potential $\mu \rightarrow 0$. The system with fixed $N$ has different asymptotics at odd and even $N\rightarrow \infty$, but the limit $\mu \rightarrow 0$ is well defined. The energy dissipation rate is analytically calculated as a function of $\mu$ using methods of number theory. It grows as $\nu/\mu^2$ in the continuum limit $\mu \rightarrow 0$, leading to anomalous dissipation at $\mu \propto \sqrt{\nu} \to 0$. The same method is used to compute all the local vorticity distribution, which has no continuum limit but is renormalizable in the sense that infinities can be absorbed into the redefinition of the parameters. The small perturbation of the fixed manifold satisfies the linear equation we solved in a general form. This perturbation decays as $t^{-\lambda}$, where the anomalous dimensions $\lambda$ satisfy the spectral equation \eqref{spectralEq}. The spectrum becomes a continuum in the statistical limit $N \to \infty$, presumably leading to multifractal phenomena.

In the early 1980s it was well established that Leray solutions of the unforced Navier-Stokes equations in Rn decay in energy norm for large time. With the works of T. Miyakawa, M. Schonbek and others it is now known that the energy decay rate cannot in general be any faster than t^-(n+2)/4 and is typically much slower. In contrast, we show in this note that, given an arbitrary Leray solution u(.,t), the difference of any two Stokes approximations to the Navier-Stokes flow u(.,t) will always decay at least as fast as t^-(n+2)/4, no matter how slow the decay of || u(.,t) ||_L2 might happen to be.

We study the existence of positive solutions for a system of Hadamard fractional differential equations on an infinite interval with nonnegative nonlinearities, subject to coupled nonlocal boundary conditions which contain Hadamard fractional derivatives and Riemann-Stieltjes integrals. In the proof of the main results, we use the Guo-Krasnosel'skii fixed point theorem and the Leggett-Williams fixed point theorem.

This article develops dual variational formulations for a large class of models in variational optimization. The results are established through basic tools of functional analysis, convex analysis and duality theory. The main duality principle is developed as an application to a Ginzburg-Landau type system in superconductivity in the absence of a magnetic field. In the first sections, we develop new general dual convex variational formulations, more specifically, dual formulations with a large region of convexity around the critical points which are suitable for the non-convex optimization for a large class of models in physics and engineering. Finally, in the last section we present some numerical results concerning the generalized method of lines applied to a Ginzburg-Landau type equation.

We point out an inconsistency in Newton's equations of celestial mechanics. A set of differential equations implied by Newton's equations are shown to be free of this inconsistency. We then investigate the integrals of motion associated with this relative difference system.

Non-standard Lagrangians do not display any discernible energy-like terms, yet they give the same equations of motion as standard Lagrangians, which have easily identifiable energy-like terms. A new method to derive non-standard Lagrangians for second-order nonlinear differential equations with damping is developed and the limitations of this method are explored. It is shown that the limitations do not exist only for those nonlinear dynamical systems that can be converted into linear ones. The obtained results are applied to selected population dynamics models for which non-standard Lagrangians and their corresponding null Lagrangians and gauge functions are derived, and their roles in the population dynamics are discussed.

The famous Toms effect (1948) consists of a substantial increase of the critical Reynolds number when a small amount of soluble polymer is introduced into water. The most noticeable influence of polymer additives is manifested in the boundary layer near solid surfaces. The goal of the present paper is a group analysis of the boundary layer equations in two mathematical models of the flow of aqueous polymer solutions: the second grade fluid (Rivlin and Ericksen, 1955) and the model derived by Pavlovskii (1971). The equations of the unsteady two-dimensional boundary layer in the Pavlovskii and Rivlin-Ericksen models are analyzed for the first time here. These equations have no definite type so that finding their exact solutions is very important in order to understand the mathematical nature of the above mentioned models. The problem of group classification with respect to the arbitrary function of the longitudinal coordinate and time present in the equations, which sets the pressure gradient of the external flow, arises. All functions for which an extension of the admitted Lie group occurs are found. The task includes the ratio of two characteristic length scales. One of them is the Prandtl scale, and another is defined as the square root of the normalized coefficient of relaxation viscosity (Frolovskaya and Pukhnachev, 2018) and does not depend on the characteristics of the motion. The paper contains a number of exact solutions in the Pavlovskii model including a solution describing the flow near a critical point. Among the solutions of the new model of the boundary layer, a special place is taken by the solution of the stationary problem of flow around a rectilinear plate. Within the framework of the Prandtl theory of the boundary layer, such a solution was constructed by Blasius (1908). As is well-known, this solution has a non-removable defect: the transverse velocity near the edge of the plate increases without bound. The introduction of a relaxation term into the model makes it possible to eliminate this singularity.

Given any sequence a = (an)n1 of positive real numbers and any set E of complex sequences, we write Ea for the set of all sequences y = (yn)n1 such that y=a = (yn=an)n1 2 E. In this paper, we use the spaces w1, w0 and w of strongly bounded, summable to zero and summable sequences, that are the sets of all sequences y such that ? n?1Pn k=1 jykj
n is bounded, tends to zero and such that y ? le 2 w0, for some scalar l, respectively, (cf. [24, 22]). These sets where used in the statistical convergence, (cf. [17, Chapter 4]). Then we deal with the solvability of each of the (SSIE) F
 E + F0x where E is a linear space of sequences, F = c0, c, `1, w0, w, or w1 and F0 = c0, c, or `1. For instance, the solvability of the (SSIE) w
 w0+s(c) x consists in determining the set of all sequences x = (xn)n1 2 U+ that satisfy the following statement. For every sequence y that satisfy the condition limn!1 n?1Pn k=1 jyk ? yk?1 ? lj = 0, there are two sequences u and v, with y = u+v such that limn!1 n?1Pn k=1 jukj = 0 and limn!1 (vn=xn) = L for some scalars l and L. These results extend those stated in [11, 12, 10].

Calculus, as a starting point to everyone's study of mathematics, was developed along with classical mechanics. Though not rigorously formulated, the intuition for differential form and infinitesimal motion makes calculus not hard to understand. And hence it does not require great effort to grasp the principles of classical mechanics. However, the situation becomes quite different when it comes to the quantum world. Given a particle system, what of interest is always the physical quantity of this system, such as energy, position, momentum, entropy, etc. But it is no longer straightforward.

When 2−1/Q<p≤2, we establish the C^{0,1} and C^{1,α}-regularities of weak solutions to quasi-linear p-Laplacian type non-homogeneous equations in the Heisenberg group.

In this paper, a short survey of the existed results concerning the Jacobian Conjecture is first given. Then the 3-fold linear type polynomial map will be analyzed in detail. The expansion of the Jacobian condition is deduced to obtain its equivalent algebraic equations, and the Jacobian condition will be analyzed to derive two coordinate transformations that can maintain the invariance of the Jacobian condition. Finally, it is proved by mathematical induction method that one general chain expression presented in this paper is just the inverse polynomial map of 3-fold linear type polynomial map, i.e. LJC(n,[3]) holds such that the Jacobian Conjecture holds.

Maxwell equations governing electromagnetic effects are being shown to be equivalent to the compressible inviscid Navier-Stokes equations applicable in fluid dynamics and representing conservation of mass and linear momentum. The latter applies subject to a generalized Beltrami condition to be satisfied by the magnetic field. This equivalence indicates that the compressible inviscid Navier-Stokes equations are Lorentz invariant as they derive directly from the Lorentz invariant Maxwell equations subject to the same Beltrami condition. In addition, the derivation results provide support for the claim that electromagnetic potentials have physical significance as demonstrated by Aharonov-Bohm effect, and are not only a convenient mathematical formulation.

We develop a quantitative microscopic theory of decaying Turbulence by studying the dimensional reduction of the Navier-Stokes loop equation for the velocity circulation. We have found an infinite dimensional manifold of solutions of the Navier-Stokes loop equation\cite{M93, M23PR} for the Wilson loop in decaying Turbulence in arbitrary dimension $d >2$. This family of solutions corresponds to a fractal curve in complex space $\mathbb C^d$, described by an algebraic equation between consecutive positions plus a nonlinear periodicity condition. We derive the constrained SDE for the evolution of the fractal curve at a fixed moment of physical time as a function of an auxiliary stochastic time. We expect this stochastic process to cover our fixed manifold of the solutions of the decaying Turbulence. The energy density of the fluid decays as $\mathcal E_0/t$, where $\mathcal E_0$ is an initial dissipation rate. Presumably, we have found a new phase of extreme Turbulence yet to be observed in real or numerical experiments.

In this article, an iterative method of two step is proposed for finding all roots simultaneously of polynomial equations. The order of convergence of the proposed algorithm is 2m, by using any iterative scheme of order m. Numerical tests are performed to confirm the theoretical results and to compare the proposed scheme with existing methods for finding all roots simultaneously.

The so-called Born-Huang ansatz is a fundamental tool in the context of ab-initio molecular dynamics, viz., it allows to effectively separate fast and slow degrees of freedom and thus treating electrons and nuclei at different mathematical footings. Here we consider the use of a Born-Huang-like expansion of the three-dimensional time-dependent Schr\"odinger equation to separate transport and confinement degrees of freedom in electron transport problems that involve geometrical constrictions. The resulting scheme consists of an eigenstate problem for the confinement degrees of freedom (in the transverse direction) whose solution constitutes the input for the propagation of a set of coupled one-dimensional equations of motion for the transport degree of freedom (in the longitudinal direction). This technique achieves quantitative accuracy using an order less computational resources than the full dimensional simulation for a prototypical two-dimensional constriction.