Coupling the Maxwell tensor to the Riemann-Christoffel curvature tensor is shown to lead to a geometricized theory of electrodynamics. While this geometricized theory leads directly to the classical Maxwell equations, it also extends their physical interpretation by giving charge density, mass density and the four-velocity that describes their motions geometric definitions. Assuming the existence of a conserved energy-momentum tensor, all solutions to the geometricized theory of electrodynamics developed here are shown to be consistent with the emergence of gravity obeying the General Relativity field equation augmented by a term that mimics the properties of dark matter and/or dark energy. Finally, due to the symmetries of the theory, the properties and phenomenology of antimatter emerge in solutions.

Although Eulerian approaches are standard in computational acoustics, they are less effective for certain classes of problems like bubble acoustics and combustion noise. A different approach for solving acoustic problems is to compute with individual particles following particle motion. In this paper, a Lagrangian approach to model sound propagation in moving fluid is presented and implemented numerically, using three meshfree methods to solve the Lagrangian acoustic perturbation equations (LAPE) in the time domain. The LAPE split the fluid dynamic equations into a set of hydrodynamic equations for the motion of fluid particles and perturbation equations for the acoustic quantities corresponding to each fluid particle. Then, three meshfree methods, the smoothed particle hydrodynamics (SPH) method, the corrective smoothed particle (CSP) method, and the generalized finite difference (GFD) method, are introduced to solve the LAPE and the linearized LAPE (LLAPE). The SPH and CSP methods are widely used meshfree methods, while the GFD method based on the Taylor series expansion can be easily extended to higher orders. Applications to modeling sound propagation in steady or unsteady fluids in motion are outlined, treating a number of different cases in one and two space dimensions. A comparison of the LAPE and the LLAPE using the three meshfree methods is also presented. The Lagrangian approach shows good agreement with exact solutions. The comparison indicates that the CSP and GFD method exhibit convergence in cases with different background flow. The GFD method is more accurate, while the CSP method can handle higher Courant numbers.

Local polynomial smoother (LPS) is a weighted local least-squares nonparametric method. It provides a local Taylor series fit of the data at any location and can be directly used in a differential equation to provide a numerical scheme. In this article, we introduce this new nonparametric idea based on local polynomial smoother, for acquiring the numerical solution of the Bagley-Torvik fractional-order differential equations. Furthermore, this paper will present a numerical comparison with some methods, such as legendre operational matrix and pseudo-spectral method. The efficiency and accuracy of the LPS method are demonstrated by two numerical examples.

The study of different forms of preconditioners for solving a system of nonlinear equations, by using Newton’s method, is presented. The preconditioners provide numerical stability and rapid convergence with reasonable computation cost, whenever chosen accurately. Different families of iterative methods can be constructed by using a different kind of preconditioners. The multi-step iterative method consists of a base method and multi-step part. The convergence order of base method is quadratic and each multi-step add an additive factor of one in the previously achieved convergence order. Hence the convergence of order of an m-step iterative method is m + 1. Numerical simulations confirm the claimed convergence order by calculating the computational order of convergence. Finally, the numerical results clearly show the benefit of preconditioning for solving system of nonlinear equations.

This article focuses on obtaining the analytical solutions for parabolic Volterra integro- differential equations in d-dimensional with different types frictional memory kernel. Based on theories of Laplace transform, Fourier transform, the properties of Fox-H function and convolution theorem, analytical solutions of the equations in the infinite domain are derived under three frictional memory kernel functions respectively. The analytical solutions are expressed by infinite series, the generalized multi-parameter Mittag-Leffler function, Fox-H function and convolution form of Fourier transform. In addition, the graphical representations of the analytical solution under different parameters are given for one-dimensional parabolic Volterra integro-differential equation with power-law memory kernel. It can be seen that the solution curves subject to Gaussian decay at any given moment.

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.

In this paper, we study the semilinear pseudo-parabolic equations $\displaystyle u_{t} - \triangle_{\mathbb{B}}u - \triangle_{\mathbb{B}}u_{t} = \left|u\right|^{p-1}u$ on a manifold with conical singularity, where $\triangle _{\mathbb{B}}$ is Fuchsian type Laplace operator investigated with totally characteristic degeneracy on the boundary $x_{1} = 0$. Firstly, we discuss the invariant sets and the vacuum isolating behavior of solutions with the help of a family of potential wells. Then, we derive a threshold result of existence and nonexistence of global weak solution: for the low initial energy $J(u_{0})<d$, the solution is global in time with $I(u_{0}) >0$ or $\displaystyle\Vert \nabla_{\mathbb{B}}u_{0}\Vert_{L_{2}^{\frac{n}{2}}(\mathbb{B})} = 0$ and blows up in finite time with $I(u_{0}) < 0$; for the critical initial energy $J(u_{0}) = d$, the solution is global in time with $I(u_{0}) \geq0$ and blows up in finite time with $I(u_{0}) < 0$. The decay estimate of the energy functional for the global solution and the estimates of the lifespan of local solution are also given.

In this paper, by means of the Banach fixed point theorem and the Krasnoselskii's fixed point theorem, we investigate the existence of solutions for some fractional neutral functional integro-differential equations involving infinite delay. This paper deals with the fractional equations in the sense of Caputo fractional derivative and in the Banach spaces. Our results generalize the previous works on this issue. Also, an analytical example is presented to illustrate our results.

Splines and quasi-interpolation operators are important both in approximation theory and applications. In this paper, we construct a family of quasi-interpolation operators for the bivariate quintic spline spaces S₅³ (∆_mn⁽²⁾). Moreover, the properties of the proposed quasi-interpolation operators are studied, as well as its applications for solving two-dimensional Burgers’ equation and image reconstruction. Some numerical examples show that these methods, which are easy to implement, provide accurate results.

In this paper we use the conformable fractional derivative to discuss some fractional linear differential equations with constant coefficients. By applying some similar arguments to the theory of ordinary differential equations, we establish a sufficient condition to guarantee the reliability of solving constant coefficient fractional differential equations by the conformable Laplace transform method. Finally, we analyze the analytical solution for a class of fractional models associated with Logistic model, Von Foerster model and Bertalanffy model is presented graphically for various fractional orders and solution of corresponding classical model is recovered as a particular case.

Using four field equations, a recently proposed theory that covers the phenomenology of classical physics at the level of the Maxwell and Einstein Field Equations (M&EFEs) but then goes further by unifying electromagnetic and gravitational phenomena in a fundamentally new way is reviewed. Predictions of the field equations are shown to be consistent with those of the M&EFEs through specific solutions; a particle-like solution representing a point charge, and two radiative solutions representing electromagnetic and gravitational waves. A unique feature of the full set of field equations is that charge and mass are treated as dynamic fields instead of being introduced as external parameters as is done with the classical M&EFEs, a feature that enables a procedure for quantizing the mass, charge and angular momentum of particle-like solutions. Finally, antimatter is naturally accommodated by the theory and definite predictions regarding the interactions of matter and antimatter with gravitational fields are made.

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.

Conventionally, partial differential equations (PDE) problems are solved numerically through discretization process by using finite difference approximations. The algebraic systems generated by this process are then finalized by using an iterative method. Recently, scientists invented a short cut approach, without discretization process, to solve the PDE problems, namely by using machine learning (ML). This is potential to make scientific machine learning as a new sub-field of research. Thus, given the interest in developing ML for solving PDEs, it makes an abundance of an easy-to-use methods that allows researchers to quickly set up and solve problems. In this review paper, we discussed at least three methods for solving high dimensional of PDEs, namely PyDEns, NeuroDiffEq, and Nangs, which are all based on artificial neural networks (ANNs). ANN is one of the methods under ML which proven to be a universal estimator function. Comparison of numerical results presented in solving the classical PDEs such as heat, wave, and Poisson equations, to look at the accuracy and efficiency of the methods. The results showed that the NeuroDiffEq and Nangs algorithms performed better to solve higher dimensional of PDEs than the PyDEns.

The main bottleneck when performing computational fluid dynamics (CFD) simulations of combustion systems is the computation and integration of the highly non-linear and stiff chemical source terms. In recent times, machine learning has emerged as a promising tool to accelerate combustion chemistry, involving the use of regression models to predict the chemical source terms as functions of the thermochemical state of the system. However, combustion is a highly nonlinear phenomenon, and this often leads to divergence from the true solution when the neural network representation of chemical kinetics is integrated in time. This is because these approaches minimize the error during training without guaranteeing successful integration with ordinary differential equation (ODE) solvers. In this work, a novel neural ODE approach to combustion modeling, ChemNODE, is developed to address this issue. The source terms predicted by the neural network are integrated during training, and by backpropagating errors through the ODE solver, the neural network weights are adjusted accordingly to minimize the difference between the predicted and actual ODE solutions. It is shown that even when the dimensionality of the thermochemical manifold is trimmed to remove redundant species, the proposed approach accurately captures the correct physical behavior and reproduces the results obtained using the full chemical kinetic mechanism.

This paper is motivated by spatio-temporal pattern in the occurrence of Leishmaniasis in Afghanistan and the relatively high number of zero counts. We hold the view that correlations that arise from spatial and temporal sources are inherently distinct. Our method decouples these two sources of correlations, there are at least two advantages in taking this approach. First, it circumvents the need to inverting a large correlation matrix, which is a commonly encountered problem in spatio-temporal analyses. Second, it simplifies the modelling of complex relationships such as anisotropy, which would have been extremely difficult or impossible if spatio-temporal correlations were simultaneously considered. We identify three challenges in the modelling of a spatio-temporal process: (1) accommodation of covariances that arise from spatial and temporal sources; (2) choosing the correct covariance structure and (3) extending to situations where a covariance is not the natural measure of association. Moreover, because the data covers a period that overlaps with the US invasion of Afghanistan, the high number of zero counts may be the result of no disease incidence or lapse of data collection. To resolve this issue, a model truncated at zero built on a foundation of the generalized estimating equations was proposed.

In this paper, we study the fractional pseudo-parabolic equations u_t + (−△)^s u + (−△)^s u_t = u log |u|. Firstly, we recall the relationship between the fractional Laplace operator (−△)^s and the fractional Sobolev space H^s and discuss the invariant sets and the vacuum isolating behavior of solutions with the help of a family of potential wells. Then, we derive a threshold result of existence of global weak solution: for the low initial energy J(u₀) < d, the solution is global in time with I(u₀) > 0 or ∥u₀∥X_0(Ω) = 0 and blows up +∞ with I(u₀) < 0; for the critical initial energy J(u₀) = d, the solution is global in time with I(u₀) ≥ 0 and blows up at +∞ with I(u₀) < 0. The decay estimate of the energy functional for the global solution is also given.

The thymus hosts the development of a specific type of adaptive immune cells called T cells. T cells orchestrate the adaptive immune response through recognition of antigen by the highly variable T-cell receptor (TCR). T-cell development is a tightly coordinated process comprising lineage commitment, somatic recombination of Tcr gene loci and selection for functional, but non-self-reactive TCRs, all interspersed with massive proliferation and cell death. Thus, the thymus produces a pool of T cells throughout life capable of responding to virtually any exogenous attack while preserving the body through self-tolerance. The thymus has been of considerable interest to both immunologists and theoretical biologists due to its multiscale quantitative properties, bridging molecular binding, population dynamics and polyclonal repertoire specificity. Here, we review mathematical modelling strategies that were reported to help understand the flexible dynamics of the highly dividing and dying thymic cell populations. Furthermore, we summarize the current challenges to estimating in vivo cellular dynamics and to reaching a next-generation multiscale picture of T-cell development.

Dietary patterns established during adolescence might play a role in adulthood disease. We examined the stability of dietary patterns (DPs) from childhood through adolescence and into young adulthood (from age 8 to 34 years). Data from 130 participants (53 females) of Saskatchewan Pediatric Bone Mineral Accrual Study (aged 8-15 years, at baseline) were included. Multiple 24-hour recalls were collected annually from 1991 to 1997, 2002 to 2005 and 2010 and 2011. Using principal component analysis, “Vegetarian-style”, “Western-like”, “High-fat, high-protein”, “Mixed” and “Snack” DPs were derived at baseline. Applied DP scores for all annual measurements were calculated using factor loading of baseline DPs and energy-adjusted food group intakes. We analyzed data using generalized estimating equations. The tracking coefficient represents correlation between baseline dietary pattern scores and all other follow-up dietary pattern scores. We found a moderate tracking for the “Vegetarian-style” (ß=0.44, P<0.001) and “High-fat, high-protein” (ß=0.39, P<0.001) DPs in females and “Vegetarian-style” DP (ß=0.30, P<0.001) in males. Remaining DPs showed a poor-to-fair tracking in both sexes. No tracking for “Western-like” DP in females was observed. Assessing overall change in DP scores from childhood to young adulthood showed an increasing trend in adherence to “Vegetarian-style” DP and decreasing trend in adherence to “High-fat, high-protein” DP by age in both sexes (P<0.001). While “Western-like” and “Mixed” DP scores increased only in males (P<0.001). These findings suggest that healthy dietary habits established during childhood and adolescence moderately continue into adulthood.

The shallow water equations are widely applied for the simulation of flow routing in rivers and floodplains, as well as for flood inundation mapping. From a mathematical point of view, they are a hyperbolic system of nonlinear partial differential equations, whose numerical integration is sometimes computationally burdensome. For this reason, the interest of many researchers has been focused on the study of simplified forms of the original set of equations, which requires less computational effort. One of the most commonly applied simplifications consists in neglecting the inertial terms, which changes the hyperbolic model to a parabolic one. The effects of such a choice on the outputs of the simulations of flooding events are controversial and an important topic of debate. In the present paper, two numerical models, recently proposed for the solution of the complete and zero-inertia forms of the shallow waters equations, are applied to several unsteady flow routing scenarios. We simulate synthetic and laboratory studies, starting from very simple geometries and moving towards complex topographies. Analyzing the role of the terms in the momentum equations, we try to understand the effect, on the computed results, of neglecting the inertial terms in the zero-inertia formulation. We analyze the computational costs.

In recent years, applying deep learning to solve physics problems has attracted much attention. Data-driven deep learning methods produce operators that can learn solutions to the whole system of partial differential equations. However, the existing methods are only evaluated on simple flow equations (e.g., Burger's equation), and only consider the generalization ability on different initial conditions. In this paper, we construct CFDBench, a benchmark with four classic problems in computational fluid dynamics (CFD): lid-driven cavity flow, laminar boundary layer flow in circular tubes, dam flows through the steps, and periodic Karman vortex street. Each flow problem includes data with various boundary conditions, fluid physical properties, and domain geometries. Compared to existing datasets, the advantages of CFDBench are (1) comprehensive. It contains common physical parameters such as velocity, pressure, and cavity fraction. (2) realistic. It is very suitable for deep learning solutions of fluid mechanics equations. (3) challenging. It has a certain learning difficulty, prompting to find models with strong learning ability. (4) standardized. CFDBench facilitates a comprehensive and fair comparison of different deep learning methods for CFD. We make appropriate modifications to popular deep neural networks to apply them to CFDBench and enable the accommodation of more changing inputs. The evaluation on CFDBench reveals some new shortcomings of existing works and we propose possible directions for solving such problems.

For any positive integer n ≥ 2, the Erdős-Straus Conjecture claims that the Diophantine equation 4/n = 1/x + 1/y + 1/z has a solution where x, y, z are also positive integers. In this paper, a directed network based on this equation is generated, with properties such as its average clustering coefficient, average path length, degree distributions, and largest strongly connected component analysed to reveal some underlying trends about the nature of the conjecture. Potential connections between different numbers, that result from satisfying a source-solution relationship for this equation, are described using the appropriate number-theoretic interpretations wherever possible, while conjectures backed by these trends are made in other instances. Additionally, a directed configuration model is used to show that the origin of several results is the degree sequence of the network. Metrics relating to the prime number nodes, specifically their in and out degrees, are also explored to yield some intriguing observations. On the whole, the aim is to highlight the viability of complex networks as a computational tool to study this general class of problems pertaining to fixed-length unit fraction splits.

In this paper, a local meshless method (LMM) based on radial basis functions (RBFs) is utilized for the numerical solution of various types of PDEs, due to the flexibility with respect to geometry and high order of convergence rate. In case of global meshless methods, the two major deficiencies are the computational cost and the optimum value of shape parameter. Therefore, research is currently focus towards localized radial basis function approximations, as the local meshless method is proposed here. The local meshless procedures is used for spatial discretization whereas for temporal discretization different time integrators are employed. The proposed local meshless method is testified in terms of efficiency, accuracy and ease of implementation using regular and irregular domains.

Recently, the conformable derivative and its properties have been introduced. In this paper, we propose and prove some new results on conformable Boundary Value Problems. First, we introduce a conformable version of classical Sturm´s separation, and comparison theorems. For a conformable Sturm-Liouville problem, Green's function is constructed, and its properties are also studied. In addition, we propose the applicability of the Green´s Function in solving conformable inhomogeneous linear differential equations with homogeneous boundary conditions, whose associated homogeneous boundary value problem has only trivial solution. Finally, we prove the generalized Hyers-Ulam stability of the conformable inhomogeneous boundary value problem.

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.

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.

Recently, representation formulae and monotonicity properties of generalized k-Bessel functions, Wk v,c., were established and studied by SR Mondal [24]. In this paper, we pursue and investigate some of their image formulae. We then extract solutions for fractional kinetic equations, involving Wk v,c, by means of their Sumudu transforms. In the process, Important special cases are then revealed, and analyzed.

In this paper, we pursue and investigate the solutions for fractional kinetic equations, involving Bessel-Struve function by means of their Sumudu transforms. In the process, one Important special case is then revealed, and analyzed. The results obtained in terms of Bessel-Struve function are rather general in nature and can easily construct various known and new fractional kinetic equations.

In this paper, we employ fuzzy Sumudu transform for solving system of linear fuzzy differential equations with fuzzy constant coefficients. The system with fuzzy constant coefficients is interpreted under strongly generalized differentiability. For this purpose, new procedures for solving the system are proposed. A numerical example is carried out for solving system adapted from fuzzy radioactive decay model. Conclusion is drawn in the last section and some potential research directions are given.

The purpose of this paper is to demonstrate the existence of an artificial magnetic monopole and to introduce new electromagnetic equations by altering an electric field and a magnetic field vectors.As a principle device, a cylindrical condenser is prepared, and a superconducting loop is inserted into it. By this conduction, radial electric fields take a role as the centripetal force and both counterclockwise and clockwise motions are induced. As a result, a stationery wave is formed in which the nodes take a part in creating a monopole as follows.First, employing the Lorentz conservations and because node of the stationary wave has no phases, the momentum k and the vector potential A vanish and instead a magnetic potential appears in order to maintain the Lorentz conservation. This magnetic potential has relationship with an electric potential, and thus consequently, a dependent relationship is obtained between an electric field and a magnetic field vectors. Using this conclusive dependent relationship, we can derive new Maxwell equation assembly which are created by altering the electric field and the magnetic field vectors. In this process, we derive a divergent equation of magnetic fields which is not zero, i.e., the existence of a magnetic monopole. Employing these newly derived Maxwell equation, an electromagnetic wave is derived whose speed is the same as one the existing Maxwell equations provide. As a monopole configuration, this paper discusses the energy gap of the vacuum, which is a result of the Dirac equation and describes a monopole as pairs between two Cooper pairs (i.e. four electrons) whose interaction is a photon. As mentioned, because the total momenta and phases are zero, this paper defines the wave function as the Dirac function and demonstrate the condensation, employing the Bloch’s theorem. Moreover, using the macroscopic basic equations, we retrace the creation of the divergent magnetic field in view of macroscopic phenomenon., which provides results in this paper.In Result section in this paper, we succeeded in demonstrating the distribution of the divergent magnetic field of monopole in terms of both microscopic and macroscopic scales. Furthermore, Discussion section describes properties a magnetic monopole should follow.

In crystal structure prediction, Newton's second law can always be applied to particles (atoms or ions) in a cell to determine their positions. However the values of the period vectors (crystal cell edge vectors) should be determined as well. Here we applied the dynamical equations of period vectors derived recently based on Newton's laws (doi:/10.1139/cjp-2014-0518) for that purpose, where the period vectors are driven by the imbalance between the internal and external pressures/stresses. For equilibrium states, they became equations of state, which essentially turn out the equilibrium conditions of crystals from mechanical point of view. For external pressure, by ignoring some effects, they reduced to the Mie-Gruneisen equation. Since the internal stress has both a full kinetic energy term and a full interaction term, the influences of both external temperature and stress/pressure on crystal structures can be calculated, then thermodynamical properties and processes were presented. Contrary to usual ideas, the equations show that pure harmonic vibrating phonons can result in thermal expansion in crystals when the external temperature is changed. Finally, crystal system collapse due to temperature and/or stress/pressure change was discussed.

Complex and real valued exact solutions to some reaction-diffusion equations are suggested by using homogeneous balance and Sine-Gordon equation expansion method. The predicted solution of finite series of some hyperbolic functions are determined by using some relations between the hyperbolic functions and trigonometric functions based on Sine-Gordon equation and traveling wave transform. The Newel-Whitehead-Segel and Zeldovich equations are solved explicitly. Some real valued solutions are depicted for some particular choice of parameters.

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.

In this work, a mixture of Elzaki transform and projected differential transform method is applied to solve a nonlinear wave-like equations with variable coefficients. Nonlinear terms can be easily manipulated by using the projected differential transformation method. The method gives the results show that the proposed method is very efficient, simple and can be applied to other applications.

Empirical equations representing, interpolating and smoothing groups of measurement results by regression methods are in widespread use in metrology and other fields of science and engineering. Standard equations for the propagation of measurement uncertainties to values computed from such equations are available but suffer from a lack of general acceptance and are only infrequently applied in practice. One reason for the slow uptake is the lack of clear methods to account for systematic error. In this paper, uncertainty propagation equations in terms of covariance matrices are generalized to allow for systematic errors in least-squares regression, and effects of the resulting uncertainties are investigated analytically. A stochastic ensemble model of systematic error is proposed for the computation of the non-diagonal elements of the weight matrix of the Generalized Least-Squares method (GLS) from the measurement uncertainty. A GLS projector formalism is described which separates the effects of measurement scatter on values calculated by the equation from those on the related “residuals”, i.e., the residual errors in the fitted data. The same projectors act similarly on the associated uncertainties and covariance matrices and permit the effects of systematic errors on the simulation covariance to be quantified. It is demonstrated that, in order to include systematic errors in the uncertainty estimates of GLS, covariance matrices may be substituted by novel, specifically defined dispersion matrices that are not specified yet by the GUM[1]. Systematic measurement errors may be estimated from various sources; a particularly easy way suggested here involves analyzing the structure in the regression residuals. It is demonstrated that GLS projectors exhibit inherent features of “error blindness” also with respect to systematic errors. [1] GUM: Guide to the Expression of Uncertainty in Measurement, http://www.bipm.org/en/publications/guides

A cloud of relativistic material particles is considered, for which the distortions of length and time caused by the presence of the Lorentz factor an order of magnitude greater than curvature of space-time by gravity. The gravitational mass of the cloud is found in the region where its size is insignificant. It is established that the coefficient of dependence of the cloud mass on the total rest mass of its particles includes a complete elliptic integral of the 2nd kind.

The propeller jet from a ship has a significant component directed upwards towards the free surface of the water, which can be used for ice management. This paper describes a comprehensive laboratory experiment where the influences of operational factors affecting a propeller wake velocity field were investigated. The experiment was done on a steady wake field to investigate the characteristics of the axial velocity of the fluid in the wake and the corresponding variability downstream of the propeller. The axial velocities and the variability recorded were time-averaged. Propeller rotational speed was found to be the most significant factor, followed by propeller inclination. The experimental results also provide some idea about the change of the patterns of the mean axial velocity distribution against the factors considered for the test throughout the effective wake field, as well as the relationships to predict the axial velocity for known factors.

Nowadays, all the stakeholders, policy makers, regulators, aircraft designers, producers, operators, etc.) are intensively working on development of the aircraft with full electric and hybrid propulsion systems. However, the technical, technological constrains (like limit on accumulator energy density) require introducing a new approach to conceptual design of such aircraft. The new methods is based on energy and mass balance evaluation. This paper analyses the identified constrains; integrates the energy and mass balance equations into the preliminary definition and calculations of the aircraft performance. By this way, the technological constrains might be transferred into the limitation on the aircraft energy and mass breakdown, that initiates a new approach to aircraft conceptual design uses the knowledge based multidisciplinary optimization. The paper describes the developed methodology for conceptual design of aircraft. It show results of implementing this new development philosophy to conceptual design of a four-seat small electric/hybrid aircraft and a special hybrid cargo UAV. The discussion of the results including got by using the emerging and enabling new technologies and new methods and solutions (including for example distributed propulsion system, unconventional forms, morphing, biomimics, etc.), demonstrates the possible implementation of the new development philosophy, new approach to aircraft conceptual design.

Free vibration analysis of the porous functionally graded circular plates has been presented on the basis of classical plate theory. The three defined coupled equations of motion of the porous functionally graded circular/annular plate were decoupled to one differential equation of free transverse vibrations of plate. The one universal general solution was obtained as linear combination of the multiparametric special functions for the functionally graded circular and annular plates with even and uneven porosity distributions. The multiparametric frequency equations of functionally graded porous circular plate with diverse boundary conditions were obtained in the exact closed-form. The influences of the even and uneven distributions of porosity, power-law index, diverse boundary conditions and the negligibled effect of the coupling in-plane and transverse displacements on the dimensionless frequencies of the circular plate were comprehensively studied for the first time. The formulated boundary value problem, the exact method of solution and the numerical results for the perfect and imperfect functionally graded circular plates have not yet been reported.

In this paper we study the class of mixed-index time fractional differential equations in which different components of the problem have different time fractional derivatives on the left hand side. We prove a theorem on the solution of the linear system of equations, which collapses to the well-known Mittag-Leffler solution in the case the indices are the same, and also generalises the solution of the so-called linear sequential class of time fractional problems. We also investigate the asymptotic stability properties of this class of problems using Laplace transforms and show how Laplace transforms can be used to write solutions as linear combinations of generalisedMittag-Leffler functions in some cases. Finally we illustrate our results with some numerical simulations.

In this study, we investigate the form of solutions, stability character and asymptotic behavior of the following four rational difference equations x_{n+1} = (1/(x_{n}(x_{n-1}±1)±1)), x_{n+1} = ((-1)/(x_{n}(x_{n-1}±1)∓1)), such that their solutions are associated with Tribonacci numbers.

In this paper, optimum parameters of Tuned Mass Dampers (TMD) are considered to control the responses of 10-story shear building under harmonic loading and 22 set of seismic records of FEMA-P695. The criterion used to obtain the optimum parameters is to select mass ratio, the frequency (tuning) and damping ratio that would result in smallest lateral displacements. State-space equations of motion are presented to compute the structural responses by developing a MATLAB file. A 10-story shear building is presented as a case study to assess the effects of TMDs on the multi-story structures. The results indicate that using TMD can reduce structural responses up to the average 20% under earthquake excitation and up to 90% under harmonic loadings. TMDs are not always effective under any type of ground motion; therefore, being aware of the given location is significant to design TMDs properly.

Interaction of an electromagnetic field with matter in a laser cavity without the assumption of a fixed direction of the transverse electric field, described by the two-level Maxwell-Bloch equations, is studied. By using a perturbative nonlinear analysis, performed near the laser threshold, we report on the derivation of the laser (3+1)D vectorial complex cubic-quintic complex Ginzburg-Landau equation. Furthermore, we study the modulational instability of the plane waves both theoretically using the linear stability analysis, and numerically, using direct simulations via the split-step Fourier method. The linear theory predicts instability for any amplitude of the primary waves. Our numerical simulations confirm the theoretical predictions of the linear theory as well as the threshold of the amplitude of perturbations. The system understudy shows a deep dependence on the laser cavity parameters, for which there appear wave patterns in accordance with the predictions from the gain spectrum.

Due to failure of the continuum hypothesis for higher Knudsen numbers, rarefied gases and microflows of gases are particularly difficult to model. Macroscopic transport equations compete with particle methods, such as DSMC to find accurate solutions in the rarefied gas regime. Due to growing interest in micro flow applications, such as micro fuel cells, it is important to model and understand evaporation in this flow regime. Here, evaporation boundary conditions for the R13 equations, which are macroscopic transport equations with applicability in the rarefied gas regime, are derived. The new equations utilize Onsager relations, linear relations between thermodynamic fluxes and forces, with constant coefficients, that need to be determined. For this, the boundary conditions are fitted to DSMC data and compared to other R13 boundary conditions from kinetic theory and Navier-Stokes-Fourier (NSF) solutions for two one-dimensional steady-state problems. Overall, the suggested fittings of the new phenomenological boundary conditions show better agreement to DSMC than the alternative kinetic theory evaporation boundary conditions for R13. Furthermore, the new evaporation boundary conditions for R13 are implemented in a code for the numerical solution of complex, two-dimensional geometries and compared to NSF solutions. Different flow patterns between R13 and NSF for higher Knudsen numbers are observed.

This article is devoted to Feller's diffusion equation which arises naturally in probabilities and physics (e.g. wave turbulence theory). If discretized naively, this equation may represent serious numerical difficulties since the diffusion coefficient is practically unbounded and most of its solutions are weakly divergent at the origin. In order to overcome these difficulties we reformulate this equation using some ideas from the Lagrangian fluid mechanics. This allows us to obtain a numerical scheme with a rather generous stability condition. Finally, the algorithm admits an elegant implementation and the corresponding Matlab code is provided with this article under an open source license.

Nowadays, customers play a very important and vital role in the field of global economy. As a result, companies give special significance to the customers to survive and grow in the field of economic competition in the modern world and increase their relationship with the buyers of their products and services throughout their lifetime. Marketing growth through mobile phones has provided further motives for performing more researches in the field of customer behavior and attitude in mobile marketing. The goal of the current study is to analyze the effective factors on customer behavior in mobile marketing. The variables used in this study are perceived ease of use, individual knowledge, user’s mobile phone technology, customers’ negative attitudes, and customers’ positive attitudes. The current research is practical in terms of objective and is descriptive-analytic in terms of methodology. Data were gathered by distributing a questionnaire among 284 students of Tehran University. Data were analyzed by structural equation modeling using Lisrel and Expert Choice software. Test results showed that ease of use, individual knowledge, mobile phone technology, positive attitude, and negative attitude variables have a meaningful effect on customer behavior in mobile marketing. The confirmation of all the assumptions of the research supports the importance of the customer behavior analysis in mobile phone services.

This paper examines the indirect effect of democracy on economic growth using a dataset of 17 MENA countries from 1990 to 2015. Democracy is assumed to affect growth through a series of channels: education, health, physical capital accumulation per labor, government consumption, and trade openness. A system of six simultaneous equations, 3SLS, is used to estimate the effect of democracy on growth through these channels. For further analysis, the countries are classified into groups according to the democratic status on the one side, and the level of income on the other. The results indicate that democracy enhances growth through its positive effect on health in all classifications of countries within the MENA region. However, the effect of democracy on growth through education and physical capital/labor is non-monotonic. Democracy always hinders growth through government size and trade openness. Once all of these indirect effects are accounted for, the overall effect of democracy on growth is negative in less democratic countries and poor countries, but positive in more democratic countries and rich countries.

This paper deals with the form, the stability character, the periodicity and the global behavior of solutions of the following four rational difference equations x_{n+1} = ((±1)/(x_{n}(x_{n-1}±1)-1)) x_{n+1} = ((±1)/(x_{n}(x_{n-1}∓1)+1)).

Fractional differential equations have become a fundamental modelling approach for understanding and simulating the many aspects of nonlocality and spatial heterogeneity of complex materials and systems. Yet, while real-order fractional operators are nowadays widely adopted, little progress has been made in extending such operators to complex-order counterparts. In this work, we introduce new definitions for the complex-order fractional Laplacian, fully consistent with the theory of averaging of smooth functions over fractal sets, and present tailored spectral methods for their numerical treatment. The proposed complex-order operators exhibit a dual particle-wave behaviour, with solutions manifesting wave-like features depending on the magnitude of the imaginary part of the fractional order. Reaction-diffusion systems driven by the complex-order fractional Laplacian exhibit unique spatio-temporal dynamics, such as equilibrium of diffusion in random materials by interference of scattered waves, conduction block and highly fractionated propagation, or the generation of completely novel Turing patterns. Taken together, our results support that the proposed complex-order operators hold unparalleled capabilities to advance the description of multiscale transport phenomena in physical and biological processes highly influenced by the heterogeneity of complex media.

The ordinary continuity equation relating the current and density of a system is extended to incorporate systems with dual (longitudinal and transverse) currents. Such a system of equations is found to have the same mathematical structure as that of Maxwell equations. The horizontal and transverse currents and the densities associated with them are found to be coupled to each other. Each of these quantities are found to obey a wave equation traveling at the velocity of light in vacuum. London's equations of super-conductivity are shown to emerge from some sort of continuity equations. The new London's equations are symmetric and are shown to be dual to each other. It is shown that London's equations are Maxwell's equations with massive electromagnetic field (photon). These equations preserve the gauge invariance that is broken in other massive electrodynamics. The duality invariance may allow magnetic monopoles to be present inside superconductors. The new duality is called the comprehensive duality transformation.

In this paper, a further mathematical analysis of some previous results of the author’s ongoing investigation regarding Prandtl’s equations is accomplished. In particular, the objective of the present work is the performance of an approximate explicit solution to Prandtl’s system of equations for a laminar, isothermal, incompressible and steady boundary layer flow of a Newtonian fluid, over an infinite horizontal flat plate.

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.

We extended the usage of the expansion method based on Sine-Gordon equation to the two dimensional Fisher equation. The relation between the trigonometric and hyperbolic functions are derived from the Sine-Gordon equation dened in two space dimension. The complex-valued traveling wave solutions to the two dimensional Fisher and Nagumo equations are set in forms a nite series of multiplications of powers of sech(.) and tanh(.) functions.

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.

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.

Master equations define the dynamics that govern the time evolution of various physical processes on lattices. In the continuum limit, master equations lead to Fokker-Planck partial differential equations that represent the dynamics of physical systems in continuous spaces. Over the last few decades, nonlinear Fokker-Planck equations have become very popular in condensed matter physics and in statistical physics. Numerical solutions of these equations require the use of discretization schemes. However, the discrete evolution equation obtained by the discretization of a Fokker-Planck partial differential equation depends on the specific discretization scheme. In general, the discretized form is different from the master equation that has generated the respective Fokker-Planck equation in the continuum limit. Therefore, the knowledge of the master equation associated with a given Fokker-Planck equation is extremely important for the correct numerical integration of the latter, since it provides a unique, physically motivated discretization scheme. This paper shows that the Kinetic Interaction Principle (KIP) that governs the particle kinetics of many body systems, introduced in [G. Kaniadakis, Physica A 296, 405 (2001)], univocally defines a very simple master equation that in the continuum limit yields the nonlinear Fokker-Planck equation in its most general form.

Starting from the tenets of human imagination, i.e., the concepts of lines, points and infinity, we provide a biological demonstration that the skeptical claim “human beings cannot attain knowledge of the world” holds true. We show that the Euclidean account of the point as “that of which there is no part” is just a conceptual device, untenable in our physical/biological realm: terms like “lines, surfaces and volumes” label non-existent, arbitrary properties. We also elucidate the psychological and neuroscientific features hardwired in our brain that lead us humans to think to points and lines as truly occurring in our environment. Therefore, our current scientific descriptions of objects’ shapes, graphs and biological trajectories in phase spaces need to be revisited, leading to a proper portrayal of the real world’s events. In order to provide also a positive account, we view miniscule bounded physical surface regions as the basic objects in a biological context in a traversal of spacetime instead of the usual Euclidean points. Our account makes it possible to erase a painstaking problem that causes many theories to break down and/or being incapable of describing extreme events: the unwanted occurrence of infinite values in equations, such as singularity in the description of black holes. We propose a novel approach, based on point-free geometrical standpoints, that banishes infinitesimals and leads to a tenable physical/biological geometry. We conclude that points, lines, volumes and infinity do not describe the world, rather they are fictions introduced by ancient surveyors of land surfaces.

Digitalisation provides access to an integrated network of information that can benefit society, and business. Building digital network and society using digital means can create something unique opportunities to strategically address sustainable development challenges for the United Nations Targets (SDG) to ensure higher productivity, education and to equality oriented society. This point of view describes the potential of digitalisation for society and business of the future. The authors revise the links between digitalisation and sustainability in the European Union countries. The methodology for the research is suggested in the paper and linear regression method is applied. The results showed tiers with five SDG, focusing on society and business, and all these tiers are fixed in the constructed equations for each SDG. The suggested solution is statistically valid and proves the novelty of research. Among digitalisation indicators, only mobile-cellular subscriptions and fixed-broadband sub-basket prices in part have no effect on researched sustainable development indicators.

This article, we study the sufficient conditions for the controllability of second-order impulsive partial neutral evolution differential systems with infinite delay in Banach spaces by using the theory of cosine families of bounded linear operators and fixed point theorem.

This paper deals with the study of the existence and non existence of solutions of a three parameter's family of nonlinear fractional differential equation with mixed-integral boundary value conditions. We consider the $\alpha$-Riemann-Liouville fractional derivative, with $\alpha \in (1,2]$. In order to deduce the existence and non existence results, we first study the linear equation, by deducing the main properties of the related Green's functions. We obtain the optimal set of parameters where the Green's function has constant sign. After that, by means of the index theory, the nonlinear boundary value problem is studied. Some examples, at the end of the paper, are showed to illustrate the applicability of the obtained results.

This article deals with the conformable double Laplace transforms and their some properties with examples and also the existence Condition for the conformable double Laplace transform is studied. Finally, in order to obtain the solution of nonlinear fractional problems, we present a modified conformable double Laplace that we call conformable double Laplace decomposition methods (CDLDM). Then, we apply it to solve, Regular and singular conformable fractional coupled burgers equation illustrate the effectiveness of our method some examples are given.

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.

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.

In this paper, we established the existence of PC-mild solutions for non local fractional impulsive functional integro-differential equations with finite delay. The proofs are obtained by using the techniques of fixed point theorems, semi-group theory and generalized Bellman inequality. In this paper, we have used the distributed characteristic operators to define the mild solution of the system. Results obtained here improve and extend some known results.

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.

Recent astronomical observations with respect to measurements in distant supernovas, cosmic microwave background and weak gravitational lensing confirm that the Universe is undergoing a phase of accelerated expansion and it has been proposed that this cosmological behavior is caused by a hypothetical dark energy which has a strong negative pressure that allows explain the expanding universe. Several theoretical ideas and models related dark the energy includes the cosmological constant, quintessence, Chaplygin gas, braneworld and tachyonic scalar fields. In this paper, we have obtained new relativistic stellar configurations considering an anisotropic fluid distribution with a charge distribution which could represents a potential model of a dark energy star. In order to investigate the effect of a quadratic equation of state in this anisotropic model we specify particular forms for the gravitational potential that allow solving the Einstein-Maxwell field equations. For these new solutions we checked that the radial pressure, metric coefficients, energy density, anisotropy factor, charge density , mass function are well defined and are regular in the interior of the star. The solutions found can be used in the development of dark energy stars models satisfying all physical acceptability conditions but the causality condition and strong energy condition are violated. We expect that these models have multiple applications in astrophysics and cosmology.

In this study, we consider delta fractional Sturm--Liouville (DFSL) initial value problems in the sense of delta Caputo and delta Riemann-Liouville (R--L) operators. One of the properties of delta fractional difference operators which makes it different from nabla counterpart is to shift its domain. This feature makes it more complex than the nabla fractional operator. We obtain sum representation of solutions for DFSL initial value problems with the help of $\mathcal{Z}-$ transformation. Moreover, we get analytical solutions of homogeneous DFSL problem within Riemann-Liouville (R--L) and Caputo sense, discrete Sturm--Liouville (DSL) problem, continuous fractional Sturm--Liouville (FSL) problem in the sense of R--L and Caputo operators, and continuous Sturm--Liouville (SL) differential problem. From this point of view, we compare all the solutions with each other. Consequently, we show that all results for these four eigenvalue problems are compatible with each other and approach to each other while the orders tends to one, i.e. $\Delta^{\mu }\left( \Delta x\left( t-\mu \right) \right)\cong D_{0^{+}}^{\mu }\left( x^{\prime }\left( t\right) \right)\cong \Delta^2x(n-1) \cong x^{\prime \prime }\left( t\right) =\lambda x\left( t\right),\ \mu\rightarrow1 $ . We support our results comparatively by tables and simulations in detail.

From 1964 to the early 1970s, Orson Anderson led a research program at the Lamont Geological Observatory in the newly-emerging field of “mineral physics”. In collaboration with colleagues Edward Schreiber and Naohiro Soga, Orson exploited the techniques of physical acoustics to study the behavior of the sound velocities of minerals at elevated pressures and temperatures. This research program also included semi-empirical studies of relationships between the bulk modulus and the molar volume of solids and the use of lattice dynamics to calculate the elastic moduli of cubic structures as a function of pressure to predict instabilities, as well as theoretical investigations of the Lagrangian vs. Eulerian formulations of finite strain equations of state.

Numerical modeling of sedimentation and erosion in reservoirs is an active field of reservoir research. However, simulation of bed-load transport phenomena has rarely been applied to other water bodies, in particular, the fluctuating backwater area. This is because the complex morphological processes between hydrodynamics and sediment transport are generally challenging to accurately predict. In this study, the refinement and application of a two-dimensional shallow-water and bed-load transport model to the fluctuating backwater area is described. The model employs the finite volume method of the Godunov scheme and saturated sediment transport equations. The model was verified against experimental data of a scaled physical model. It was then applied to actual reservoir operation, including reservoir storage, reservoir drawdown and continuous flood process, to predict the morphology of reservoir sedimentation and sediment transport rates and bed level changes in the fluctuating backwater area. It was found that the location and morphology of sedimentation effected by the downstream water level results in random evolution of the river bed, and bed-load sedimentation is transported from upstream to downstream with the slope of the longitudinal section of the river bed generally reduced. Moreover, the sediment is mainly deposited in the main channel and the elevation difference between the riverbank and channel decreases gradually.

In this paper, we investigate the equilibrium points of following a system of difference equations x_n+1 = x_n 2y_n−1, y _n+1 = y_n 2x_n−1. We also study the asymptotic stability of related system of difference equations. Further we examine the periodic solutions of related system with period two. Additionally, we find out the invariant interval and periodic cycles of related system of difference equations.

The present paper discusses a study of a class of Charlier matrix polynomials and its generalized analogue. Certain generating matrix functions, recurrence matrix relations, matrix differential equation, summation formulas and many new results have been discussed for these matrix polynomials. Weisner's group theoretic method is used to obtain matrix generating relations for Charlier matrix polynomials and the details of this method were given in this paper. Finally, we will discuss only briefly the procedure followed.

In this paper, mean-field type games between two players with backward stochastic dynamics are defined and studied. They make up a class of non-zero-sum differential games where the players' state dynamics solve backward stochastic differential equations (BSDEs) that depend on the marginal distributions of player states. Players try to minimize their individual cost functionals, also depending on the marginal state distributions. Under some regularity conditions, we derive necessary and sufficient conditions for existence of Nash equilibria. Player behavior is illustrated by numerical examples, and is compared to a centrally planned solution where the social cost, the sum of player costs, is minimized. The inefficiency of a Nash equilibrium, compared to socially optimal behavior, is quantified by the so-called price of anarchy. Numerical simulations of the price of anarchy indicate how the improvement in social cost achievable by a central planner depends on problem parameters.

When solving the Einstein's equations for an isolated system of masses, V. Fock introduces harmonic reference frame and obtains an unambiguous solution. Further, he concludes that there exists a harmonic reference frame which is determined uniquely apart from a Lorentz transformation if suitable supplementary conditions are imposed. It is known that wave equations keep the same form under Lorentz transformations. Thus, we speculate that Fock's special harmonic reference frames may have provided us a clue to derive the Einstein's equations in some special class of non-inertial reference frames. Following this clue, generalized Einstein's equations in some special non-inertial reference frames are derived based on the theory of vacuum mechanics. If the field is weak and the reference frame is quasi-inertial, these generalized Einstein's equations reduce to Einstein's equations. Thus, this theory may also explain all the experiments which support the theory of general relativity. There exist some differences between this theory and the theory of general relativity.

In this paper, we investigate an analytical solution of a class of nonlinear fractional integro-differential equation base on a generalized fractional power series expansion. The fractional derivatives are described in the conformable's type. The new approach is a modified form of the well-known Taylor series expansion. The illustrative examples are presented to demonstrate the accuracy and effectiveness of the proposed method

We show that the quaternion quantum mechanics has well-founded mathematical roots and can be derived from the model of elastic continuum by French mathematician Augustin Cauchy, i.e., it can be regarded as representing physical reality of elastic continuum. Starting from the Cauchy theory (classical balance equations for isotropic Cauchy-elastic material) and using the Hamilton quaternion algebra we present a rigorous derivation of the quaternion form of the non- and relativistic wave equations. The family of the wave equations and the Poisson equation are a straightforward consequence of the quaternion representation of the Cauchy model of the elastic continuum. This is the most general kind of quantum mechanics possessing the same kind of calculus of assertions as conventional quantum mechanics. The problem of the Schrödinger equation, where imaginary ‘i’ should emerge, is solved. This interpretation is a serious attempt to describe the ontology of quantum mechanics, and demonstrates that, besides Bohmian mechanics, the complete ontological interpretations of quantum theory exists. The model can be generalized and falsified. To ensure this theory to be true, we specified problems allowing exposing its falsity.

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.

By defining a constant probabilistic orbit of and iterated functions, the stability dynamics of these functions in possible interactions through connectivity provides the formation of a dynamic fixed point as a metric space between both iterated functions. The presence of a dynamic fixed point identifies qualitatively phases of iteration time lengths and interaction orbits of the event. Qualitative results show that the greater the average distance from one of the functions to the fixed point of the other (all possible solutions), the higher the iteration expression on time (false asymptotic effect) of one of the functions and in the opposite hand, the lower the average distance, the higher orbit’s interactions proximity between iterated functions (stability). This feature reveals asymptotic (well-defined) behavior between functions f and g within a well-defined Lyapunov stability.

This paper is related to the fractional view analysis of Helmholtz equations, using innovative analytical techniques. The fractional analysis of the proposed problems has been done in terms of Caputo-operator sense. In the current methodology, first, we applied the r-Laplace transform to the targeted problem. The iterative method is then implemented to obtain the series form solution. After using the inverse transform of the r-Laplace, the desire analytical solution is achieved. The suggested procedure is verified through specific examples of the fractional Helmholtz equations. The present method is found to be an effective technique having a closed resemblance with the actual solutions. The proposed technique has less computational cost and a higher rate of convergence. The suggested methods are therefore very useful to solve other systems of fractional order problems.

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.

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, we study the third order semilocal convergence of the Newton-like method for finding the approximate solution of nonlinear operator equations in the setting of Banach spaces. First, we discuss the convergence analysis under ω-continuity condition, which is weaker than the Lipschitz and Hölder continuity conditions. Second, we apply our approach to solve Fredholm integral equations, where the first derivative of involved operator not necessarily satisfy the Hölder and Lipschitz continuity conditions. Finally, we also prove that the R-order of the method is 2p + 1 for any p $\in$ (0,1].

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.

By using the new recently proposed extended rational methods, we set some solitary and multi wave solutions of the coupled Higgs system in the present study. Some new families of hyperbolic and trigonometric solutions of the coupled Higgs equations are successfully obtained. Some solutions are real as some are complex valued functions. Particular forms of some solutions derived by choice of some parameters are demonstrated in three dimensional spaces. The complex valued solutions are also depicted by plotting both real and complex component.

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.

J. C. Maxwell, B. Riemann and H. Poincar$\acute{e}$ have proposed the idea that all microscopic particles are sink flows in a fluidic aether. Following this research program, a previous theory of gravitation based on a mechanical model of vacuum and a sink flow model of particles is generalized by methods of special relativistic continuum mechanics. In inertial reference frames, we construct a tensorial potential which satisfies the wave equation. Inspired by the equation of motion of a test particle, a definition of a metric tensor of a Riemannian spacetime is introduced. Applying Fock's theorem, generalized Einstein's equations in inertial systems are derived based on some assumptions. These equations reduce to Einstein's equations in case of weak field in harmonic reference frames. There exist some differences between this theory and Einstein's theory of general relativity.

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 paper, we proposed the numerical method called the variational iteration orthogonal collocation method (VIOCM), for the approximate solution of the deadly Corona virus model using Mamadu-Njoseh polynomials as basis functions. The proposed method is an elegant mixture of the variational iteration method (VIM) and the orthogonal collocation method (OCM). It was observed that the proposed method converges rapidly to the exact solution even as N increases. This suggests that the use of orthogonal polynomials as trial functions for the SEIR model is indeed an effective approximant as it produces the analytic solution at just few iterations. Resulting numerical evidences were compared with the standard variational iteration method as available in literature. All computational frameworks were executed with MAPLE software.

The technique of superposition of motions in the space of Lagrange variables is described, which allows us to obtain the equations of combined motion by replacing the Lagrange variables of superimposed (external) motion with Euler variables of nested (internal) motion. The components of velocity and acceleration in the combined motion obtained as a result of differentiating the equations of motion in time coincide with the results of vector addition of the velocities and accelerations of the particles involved in the superimposed motions at each moment of time. Examples of motion and superposition descriptions for absolutely solid and deformable bodies with equations for the main kinematic characteristics of motion, including for robot manipulators with three independent drives, pressing with torsion, bending with tension, and cross– helical rolling, are given. For example, given the fragment of calculation of forces in the kinematic pairs shown the advantages of the description of motion in Lagrangian form for the dynamic analysis of lever mechanisms, allows to determine the required external exposure when performing the energy conservation law at any time in any part of the mechanism.

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.

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.

In this paper, we study the dynamics of following system of nonlinear difference equations x_n+1 = x_n-1y_n - 1, y_n+1 = y_n-1z_n - 1, z_n+1 = z_n-1x_n - 1. Especially we investigate the periodicity, boundedness and stability of related system of difference equations.

The deep energy method (DEM), a type of physics-informed neural network, is evolving as an alternative to finite element analysis. This method employs the principle of minimum potential energy to predict deformations under static loading conditions. However, the model’s accuracy is contingent upon choosing the appropriate architecture for the model, which can be challenging due to the high interactions between hyperparameters, large search space, difficulty in identifying objective functions, and non-convex relationships with the objective functions. To improve DEM’s accuracy, we first introduce random Fourier feature (RFF) mapping. RFF mapping helps with the training of the model by reducing bias towards high frequencies. The effects of six hyperparameters are then studied under compression, tension, and bending loads in planar linear elasticity. Based on this study, a systematic automated hyperparameter optimization approach is proposed. Due to the high interaction between hyperparameters and the non-convex nature of the optimization problem, Bayesian optimization algorithms are used. The models trained using optimized hyperparameters and having Fourier feature mapping can accurately predict deflections compared to finite element analysis. Additionally, the deflections obtained for tension and compression load cases are more sensitive to variations in hyperparameters than bending.

Exact and soliton type solutions have great importance in propagation of surface waves, fluid dynamics, optics, and many other elds of nonlinear sciences. In this study, the explicit and exact soliton type solutions for two space-time fractional Equal- Width (FEW) equations with conformable derivative are procured via the hyperbolic function approach. The wave type solutions are represented in some hyperbolic and trigonometric functions.

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 this work, we investigate the behaviour of relativistic quantum oscillator under the effects of Lorentz symmetry violation determined by a tensor $(K_F)_{\mu\nu\alpha\beta}$ out of the Standard Model Extension. We analyze this relativistic system under an inverse radial electric field and a constant magnetic field induced by Lorentz symmetry violation. We see that the presence of Lorentz symmetry breaking terms modified the energy spectrum of the system, and a quantum effect arise due to the dependence of the linear charge density on the quantum numbers of the system

Physics-informed neural networks are a promising method to yield surrogate models of flow fields. We present a metamodeling technique for variable geometries based on physics-informed neural networks. The method was applied to the DU99W350 airfoil at a Reynolds number of 1×10⁵. The model predicted the Reynolds-averaged velocity and pressure field around the airfoil for arbitrary angles of attack between 10.0° and 17.5°. The model was trained with data from CFD simulations for a limited set of angles of attack. Additionally, satisfaction of the a priori known boundary conditions as well as the Reynolds-averaged Navier-Stokes equations were trained. A sensitivity analysis concerning the Reynolds number, the amount and distribution of training data, and the turbulence model was conducted showing the superiority of the pseudo-Reynolds stress method and the demand of labeled training data in the domain. The trained network was capable of predicting the developing flow separation on the suction surface and exhibited excellent agreement with CFD results even in the proximity to the wall for interpolations as well as extrapolations from the labeled data set.

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.

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$\%$.

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].

Conditional heteroskedastic financial time series are commonly modelled by ARCH and GARCH. ARCH(1) and GARCH processes were recently extended to the function spaces C[0,1] and L²[0,1], their probabilistic features were studied and their parameters were estimated. The projections of the operators on finite-dimensional subspace were estimated, as were the complete operators in GARCH(1,1). An explicit asymptotic upper bound of the estimation errors was stated in ARCH(1). This article provides sufficient conditions for the existence of strictly stationary solutions, weak dependence and finite moments of ARCH and GARCH processes in various L^p[0,1] spaces, C[0,1] and other spaces. In L²[0,1] we deduce explicit asymptotic upper bounds of the estimation errors for the shift term and the complete operators in ARCH and GARCH and for the projections of the operators on a finite-dimensional subspace in ARCH. The operator estimaton is based on Yule-Walker equations. The estimation of the GARCH operators also involves a result concerning the estimation of the operators in invertible, linear processes which is valid beyond the scope of ARCH and GARCH. Through minor modifications, all results in this article regarding functional ARCH and GARCH can be transferred to functional ARMA.

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.