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

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

We consider the second order differential equation y'' + ω^(2)ρ(x)y = 0 where ω is a positive parameter. The principal concern here is to find conditions on the function ρ^(−1/2)(x), which ensure that the consecutive differences of sequences constructed from the zeros of a nontrivial solution of the equation are regular in sign for ω sufficiently large. In particular, if c_(νk)(α) denotes the k-th positive zero of the general Bessel (cylinder) function C_(ν)(x; α) = J_(ν)(x) cos α−Y_(ν)(x) sin α of order ν, and if |ν| < 1/2, we prove that (−1)^(m)∆^(m+2)c_(νk)(α) > 0 (m = 0, 1, 2, ...; k = 1, 2, ...), where ∆a_(k) = a_(k+1) − a_(k). This type of inequalities was conjectured by Lorch and Szego in 1963. We also show that the differences of the zeros of various orthogonal polynomials with higher degrees possess the sign-regularity.

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

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

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

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

The article considers an internal boundary value problem of the distribution of an electrostatic field in a lens formed by two identical semi-infinite circular cylinders coaxially located inside an infinite external cylinder. The problem is reduced to solving a system of singular Wiener-Hopf integral equations, which is further solved by the Wiener-Hopf method using factorized Bessel functions. Solutions to the problem for each region inside the infinite outer cylinder are presented as exponentially converging series in terms of eigenfunctions and eigenvalues.

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.

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

Hessenberg differential algebraic equations (Hessenberg-DAEs) with high index play a critical role in the modeling of mechanical systems and multibody dynamics. Motivated by the widely used Lie Group Differential Algebraic Equation (LGDAE) method which only handles index 2 systems, we propose a Modified Extended Lie Group Differential Algebraic Equation (MELGDAE) method for solving index 3 Hessenberg-DAEs, and provide theoretical analysis to deepen the foundation of the MELGDAE method.The performance of the MELGDAE method is compared with the standard methods RADAU and MEBDF on index 2 and 3 DAE systems, and it is demonstrated that the MELGDAE integrator exhibits competitive performance in terms of high accuracy and the preservation of algebraic constraints. In particular, all differential variables in index 3 Hessenberg DAEs achieve second-order convergence using the MELGDAE method, which suggests potential for extension to Hessenberg-DAEs with an index of 4 or higher.

This article proposes the multilayer and multilevel cosmological model, on average, balanced with respect to the Einsteinian vacuum. This hierarchical model is based on 10 possible Kottler solutions of the extended Einstein's field equations. The extension of the Einstein field equations is associated with the addition of an infinite number of lambda_j-terms to these equations with the condition that the series (i.e., the total sum) of these terms converges to zero. A simplified ten-layer and ten-level case is considered, in which ten spherical formations are sequentially nested into each other like "nesting dolls", while the radii of these spherical formations correspond to the characteristic sizes of a discrete sequence of observed spherical objects: metagalaxy core, galactic cores, stellar (planet) core, biological cells and nuclei of atoms and elementary particles, etc. As an example, from the general ten-level (hierarchical) solution of extended Einstein's field equations, two antipodal spherical formations with a core radius commensurate with the classical electron radius are singled out. Therefore, the metric-dynamic models of these formations obtained in this way are called «electron» and «positron». The article is aimed at the development of differential geometry and the program for the complete geometrization of physics by Clifford-Einstein-Wheeler.

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.

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

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

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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.

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.

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.

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.

Solving Linear equations with large number of variable contains many computations to be performed either iteratively or recursively. Thus it consumes more time when implemented in a sequential manner. There are many ways to solve the linear equations such as Gaussian elimination, Cholesky factorization, LU factorization, QR factorization. But even these methods when implemented on a sequential platform yield slower results as compared to a parallel platform where the time consumption is reduced considerably due to concurrent execution of instructions. The above mentioned linear equation solving methods can be implemented on the parallel platform using the direct approaches such as pipelining or 1D and 2D Partitioning approach. Vedic mathematics is a very ancient approach for solving mathematical problems. These Vedic mathematical approaches are well known for quicker and faster computation of mathematical problems. Vedic Mathematics provides a very different outlook towards the approach of solving linear Equation on parallel platform. It could be considered as a better approach for reducing space consumption and minimizing the number of algebraic operations involved in solving linear equation. In future Vedic Mathematics might serve as a viable solution for solving linear equation on parallel platform.

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

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

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

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

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

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

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

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

Improving the operational reliability of nuclear power plants, combined heat and power plants (CHP), as well as oil and gas pipelines is a priority task in the development of a variable-speed drive for lock valves used at these facilities. The paper analyzes technical requirements for such devices; the motor has been selected, its electrical equilibrium and moment equations have been obtained; recommendations for the selection of the kinematic drive scheme have been formulated. Based on the theoretical data obtained, a prototype has been developed, manufactured and tested.

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.

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

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

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

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

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.

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

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

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

In this 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

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

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

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

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

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

In this work, we investigate the behaviour of a relativistic scalar particle in the background of the Lorentz symmetry violation determined by a tensor (KF)µναβ out of the Standard Model Extension. A linear electric field and a uniform magnetic can be induced by the violation of the Lorentz symmetry breaking effects, and analyze the behaviour of the scalar particle. We see that the analytical solution to the KG-equation can be achieved, and a quantum effect characterized by the dependence of the magnetic field on the quantum numbers is observed

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

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

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

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

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.

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.

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

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.

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.

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

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

New class of conserved quantities is constructed. These quantities find direct application in mechanics of dissipative (generally non-conservative) dynamical systems. Approach demands formulation in the language of geometric mechanics, providing theoretical framework for situations with energy flow in and out of the system. As a by product, we suggest possibility of existence of Hamiltonian form for every autonomous ODE system, evolution of which is governed by non-potential generator of motion. Various examples are provided, ranging from physics and mathematics, to chemical kinetics and population dynamics in biology. Applications of these ideas in geometric integration techniques (GNI) of numerical analysis are discussed, and as an example of such, new discrete gradient-based numerical method is introduced.

The Maxwell equations of electrodynamics describe electrical and magnetic forces with great accuracy in the vacuum of space. But the equations of electrodynamics applied to material systems are usually written in a way that obscures their accuracy. The usual formulation of the Maxwell equations includes a dielectric constant as a single real number that grossly over-approximates the actual properties of matter, particularly liquid matter so important in applications of electrodynamics to biology and chemistry. We rewrite two Maxwell equations here to make clear the precision of the Maxwell equations in the presence or absence of matter. We discuss and derive two well known corollaries that are as universal and precise as the Maxwell equations themselves: (1) a continuity equation that relates charge and flux and (2) a conservation equation in which total current never accumulates at all. The total current is the right hand side of the Ampere-Maxwell law. Total current is perfectly incompressible. It is conserved exactly, everywhere and at every time, independent of the microphysics of charge conduction. The total current combines the flux of charges and the ethereal current. The ethereal current \varepsilon_{o\ }\sfrac{\partial\mathbf{E}}{\partial t}\ exists everywhere, including the interior of atoms, because it is a property of space, not matter. The ethereal current exists, and thus flows, in a vacuum devoid of mass. The ethereal current is a consequence of the relativistic (Lorentz) invariance of charge. Charge does not change even if it moves at speeds approaching the velocity of light. Total current has properties quite distinct from flux of mass because of the ethereal current. Most strikingly, in the one dimensional unbranched systems of electronic circuits and biological ion channels, the total current is independent of location, even if the flux of charges varies a great deal. Indeed, the Maxwell equations—and thus conservation of total current—act as a perfect low pass spatial filter, converting the infinite variation of (the Brownian model of) thermal motion of charges to the zero variation of the total current.

In this paper, we found a new model for compact star with charged anisotropic matter distribution considering an extended version of the Chaplygin equation of state. We specify a particular form of the metric potential Z(x) that allows us to solve the Einstein-Maxwell field equations. The obtained model satisfies all physical properties expected in a realistic star such that the expressions for the radial pressure, energy density, metric coefficients, measure of anisotropy and the mass are fully well defined and are regular in the interior of star. The solution obtained in this work can have multiple applications in astrophysics and cosmology.

Several problems of radiative transfer are yet unsolved because of the difficulties of the calculations involved in them, especially if the intervening shapes are geometrically complex. The main goal of our investigation in this domain is to convert the formulas that were previously derived, into a graphical interface based on the projected solid-angle principle. Such procedure is now feasible by virtue of several widely diffused programs for Algorithms Aided Design (AAD). Accuracy and reliability of the process is controlled by means of the analytical software DianaX developed at an earlier stage by the authors. With this new approach the often cumbersome procedure of lighting and thermal exchange calculations can be simplified and made available for the neophyte, with the undeniable advantage of reduced computer time.

This paper presents a new analytical method for obtaining new deformations compatibility equations or, new Saint-Venant’s identities, into the relative coordinate system with projected deformations by applying the hypothesis of the lineal shell theory in general flexion state. The method proposed generalizes the compatibility conditions established by A.L. Goldenveizer for the shell theory. On the other hand, the new equations include the deformations compatibility equations by other authors: Flügge, Saint-Venant, Love-Kirchhoff, Timoshenko, Goldenveizer, and Reissner-Mindlin. The results showed an increase of knowledge in general shell theory, and provide inverse and semi-inverse solutions, whose systems solution correspond to the hyper-statics degrees of the physical model, and not to their degrees of freedom.

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.

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.

In this paper, we investigate the dynamics of following higher order difference equation x_{n+1}=A+B((x_{n})/(x_{n-m}²)) with A,B and initial conditions are positive numbers, and m∈{2,3,⋯}. Especially we study the boundedness, periodicity, semi-cycles, global asymptotically stability and rate of convergence of solutions of related higher order difference equations.

In this article, author performs computational analysis for the generalized zeta functions by using computational software Mathematica. To achieve the purpose recently obtained difference equations are used. These difference equations have a computational power to compute these functions accurately while they can not be computed by using their known integral represenations. Several authors investigated such functions and their analytic properties, but no work has been reported to study the graphical representations and zeors of these functions. Author performs numerical computations to evaluate these functions for different values of the involved parameters. Taylor series expansions are also presented in this research.

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

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

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

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

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

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

Thermodynamics has been the foundation of many models of biological and technological systems. But thermodynamics is static and is misnamed. A more suitable name is thermostatics. Thermostatics does not include time as a variable and so has no velocity, flow or friction. Indeed, as usually formulated, thermostatics does not include boundary conditions. Devices require boundary conditions to define their input and output. They usually involve flow and friction. Thermostatics is an unsuitable foundation for understanding technological and biological devices. A time dependent generalization of thermostatics that might be called thermal dynamics is being developed by Chun Liu and collaborators to avoid these limitations. Electrodynamics is not restricted like thermostatics, but in its classical formulation involves drastic assumptions about polarization and an over-approximated dielectric constant. Once the Maxwell equations are rewritten without a dielectric constant, they are universal and exact. Conservation of total current, including displacement current, is a restatement of the Maxwell equations that leads to dramatic simplifications in the understanding of one dimensional systems, particularly those without branches, like the ion channel proteins of biological membranes and the two terminal devices of electronic systems. The Brownian fluctuations of concentrations and fluxes of ions become the spatially independent total current, because the displacement current acts as an unavoidable low pass filter, a consequence of the Maxwell equations for any material polarization. Electrodynamics and thermal dynamics together form a suitable foundation for models of technological and biological systems.

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.