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

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

The main purpose of this article is to consider the smoothness control of a weak solution after some moment if there is known solution regularity until this moment. The necessary tools can be varied. It is possible to control kinetic energy dissipation to fix moment or changing of velocity square or summability of acceleration square to fix at point time.

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

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.

We introduce the linear subgroup of volume-preserving diffeomorphism as the underlying symmetry to construct an action within the effective field theory framework and the Schwinger-Keldish formalism. The formulated action is posited to effectively incorporate dissipative effects into the Navier-Stokes equations.

Monge-Ampere equations have important research significance in many fields such as geometry, convex geometry and mathematical physics. In this paper, under some superlinear and sublinear conditions, the existence of nontrivial solutions for a system arising from Monge-Ampere equations with two parameters is investigated based on Guo-Krasnosel'skii fixed point theorem. In the end, two examples are given to illustrate our main results.

Recently, Li et al [1] investigated the nonlinear Langevin equations existence including two fractional orders with infinite-point boundary conditions. It has shown that the outcomes given it does count on the solution form for which it has boundary values. It means that their solution needs more to be in the final step. In current work, we will have discussion about the existence and uniqueness for the same boundary conditions by investigating the closely explicit solution form that has no boundary values. Numerical example is supported to have more vision to compare the between new and previous results. It turns out that our results are superior.

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

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

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.

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.

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

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

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

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.

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.

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

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

The asymptotic dynamics of the classes of rational difference equations (RDEs) of third order defined over the positive real-line as $$\displaystyle{x_{n+1}=\frac{x_{n}}{ax_n+bx_{n-1}+cx_{n-2}}}, \displaystyle{x_{n+1}=\frac{x_{n-1}}{ax_n+bx_{n-1}+cx_{n-2}}}, \displaystyle{x_{n+1}=\frac{x_{n-2}}{ax_n+bx_{n-1}+cx_{n-2}}}$$ and $$\displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n}}}, \displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n-1}}}, \displaystyle{x_{n+1}=\frac{ax_n+bx_{n-1}+cx_{n-2}}{x_{n-2}}}$$ is investigated computationally with theoretical discussions and examples. It is noted that all the parameters $a, b, c$ and the initial values $x_{-2}, x_{-1}$ and $x_0$ are all positive real numbers such that the denominator is always positive. Several periodic solutions with high periods of the RDEs as well as their inter-intra dynamical behaviours are studied.

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.

In this paper, we compare the performance between systems of ordinary and (Caputo) fractional differential equations depicting the susceptible-exposed-infectious-recovered (SEIR) models of diseases. In order to understand the origins of both approaches as mean-field approximations of integer and fractional stochastic processes, we introduce the fractional differential equations as approximations of some type of fractional nonlinear birth--death processes. Then, we examine validity of the two approaches against empirical courses of epidemics; we fit both of them to case counts of three measles epidemics that occurred during the pre-vaccination era in three different locations. While FDEs appear more flexible in fitting empirical data, our ODEs offered better fits to two out of three data sets. Important differences in transient dynamics between these modeling approaches are discussed.

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

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

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

This paper provides the numerical investigation of two-dimensional stochastic Navier-Stokes equation driven by additive white noises. The vorticity-stream function method is presented to calculate the incompressible fluid flow. In order to solve the formulation of steam function equation, we propose a Crank-Nicolson Fourier pseudo-spectral method. Numerical results obviously show that the presence of stochastic forcing can have an impact upon the shapes of vortex in fluid flow. To a certain extent, the phenomenon of eddy shedding and the formation of some new vortex in a random medium can be observed, which make it possible for us to conduct optimal control experiments in the development of vortex structures.

New possibilities of Gramian computation by using canonical transformations into diagonal, controllable and observable canonical forms are shown. With the help of such a technique the Gramian matrices can be represented in the form of products of Hadamard matrices of multipliers and matrices of the transformed right-hand side of Lyapunov equations. It is shown that the multiplier matrices are invariant under various canonical transformations of linear continuous systems. The modal Lyapunov equations for continuous SISO LTI systems in diagonal form are obtained and their new solutions based on Hadamard decomposition are proposed. New algorithms for element-by-element computation of Gramian matrices for stable continuous MIMO LTI systems are developed. For continuous SISO LTI systems given by equations of state in controllable and observable canonical forms, new algorithms for the computation of controllability Gramians and their traces in the form of Hadamard products in the form of Xiao matrices are developed. The application of transformations to the canonical forms of controllability and observability allowed to simplify the formulas of spectral decompositions in the form of Xiao matrices. In the paper new spectral decompositions in the form of Hadamard’s products for solutions of Sylvester algebraic and differential equations of MIMO LTI systems, including spectral decompositions of finite and infinite cross-Gramians of continuous MIMO LTI systems. Recommendations on the use of the obtained results are given.

This study focuses on the optimal deployment problem, and determines the optimal size of a military force to send to the battle field. The decision is optimized, based on an objective function, that considers the cost of deployment, the cost of the time it takes to win the battle, and the costs of killed and wounded soldiers with equipment. The cost of deployment is modeled as an explicit function of the number of deployed troops and the value of a victory with access to a free territory, is modeled as a function of the length of the time it takes to win the battle. The cost of lost troops and equipment, is a function of the size of the reduction of these lives and resources. An objective function, based on these values and costs, is optimized, under different parameter assumptions. The battle dynamics is modeled via the Lanchester differential equation system based on the principles of directed fire. First, the deterministic problem is solved analytically, via derivations and comparative statics analysis. General mathematical results are reported, including the directions of changes of the optimal deployment decisions, under the influence of alternative types of parameter changes. Then, the first order optimum condition from the analytical model, in combination with numerically specified parameter values, is used to determine optimal values of the levels of deployment in different situations. A concrete numerical case, based on documented facts from the Battle of Iwo Jima, during WW II, is analyzed, and the optimal US deployment decisions are determined under different assumptions. The known attrition coefficients of both armies, from USA and Japan, and the initial size of the Japanese force, are parameters. The analysis is also based on some parameters without empirical documentation, that are necessary to include to make optimization possible. These parameter values are motivated in the text. The optimal solutions are found via Newton- Raphson iteration. Finally, a stochastic version of the optimal deployment problem is defined. The attrition parameters are considered as stochastic, before the deployment decisions have been made. The attrition parameters of the two armies have the same expected values as in the deterministic analysis, are independent of each other, have correlation zero, and have relative standard deviations of 20%. All possible deployment decisions, with 5000 units intervals, from 0 to 150000 troops, are investigated, and the optimal decisions are selected. The analytical, and the two numerical, methods, all show that the optimal deployment level is a decreasing function of the marginal deployment cost, an increasing function of the marginal cost of the time to win the battle, an increasing function of the marginal cost of killed and wounded soldiers and lost equipment, an increasing function of the initial size of the opposing army, an increasing function of the efficiency of the soldiers in the opposing army and a decreasing function of the efficiency of the soldiers in the deployed army. With stochastic attrition parameters, the stochastic model also shows that the probability to win the battle is an increasing function of the size of the deployed army. When the optimal deployment level is selected, the probability of a victory is usually less than 100%, since it would be too expensive to guarantee a victory with 100% probability. Some of many results of relevance to the Battle of Iwo Jima, are the following: In the deterministic Case 0 analysis, the optimal US deployment level is 66200, the time to win the battle is 30 days and 14000 US soldiers are killed or wounded. If the marginal cost of the time it takes to win a victory doubles, the optimal deployment increases to 75400, the time to win a victory is reduced to 26 days, and less than 12000 soldiers are killed or wounded. In the stochastic Case 0 analysis, the optimal US deployment level is 65000, the expected time to win the battle is 46 days and almost 25000 US soldiers are expected to be killed or wounded. If the cost per killed or wounded soldier increases from 0 to 5 M $US, the optimal deployment level increases to 75000. Then, the victory is expected to appear after 35 days and 19900 US soldiers are expected to be killed or wounded.

Dual quaternions have wide applications in automatic differentiation, computer graphics, mechanics and others. Due to its application in control theory, the matrix equation AXB=C has been extensively studied. However, there is currently limited information on matrix equation AXB=C over the dual quaternion algebra. In this paper, we provide the necessary and sufficient conditions for the solvability of the dual quaternion matrix equation AXB=C, and present the expression for the general solution when it is solvable. As an application, we derive the ϕ-Hermitian solutions for the dual quaternion matrix equation AXAϕ=C, where the ϕ-Hermitian extends the concepts of Hermiticity and η-Hermiticity. Lastly, we present a numerical example to verify the main research results of this paper.

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

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

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

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.

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.

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

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

This is the second paper in a cycle investigating the exact solution of loop equations in decaying turbulence. We perform numerical simulations of the Euler ensemble, suggested in the previous work, as a solution to the loop equations. We designed novel algorithms for simulation, which take a small amount of computer RAM so that only the CPU time grows linearly with the size of the system and its statistics. This algorithm allows us to simulate systems with billions of discontinuity points on our fractal curve, dual to the decaying turbulence. For the vorticity correlation function, we obtain \textbf{quantum fractal laws} with regime changes between universal discrete values of indexes. The traditional description of turbulence with fractal or multifractal scaling laws does not apply here. The measured conditional probabilities of fluctuating variables are smooth functions of the logarithm of scale, with statistical errors being negligible at our large number of random samples $T = 167,772,160$ with $N = 200,000,000$ points on a fractal curve. The quantum jumps arise from the analytical distribution of scaling variable $X = \frac{\cot^2(\pi p/q)}{q^2}$ where $\frac{p}{q}$ is a random fraction. This distribution is related to the totient summatory function and is a discontinuous function of $X$. In particular, the energy spectrum is a devil's staircase with tilted uneven steps with slopes between $-1$ and $-2.03$. The logarithmic derivative of energy decay $n(t)$ as a function of time is jumping down the stairs of universal levels between $1.50$ and $1.53$. The quantitative verification of this quantization would require more precise experimental data.

A cosmological model described by the perturbation metric is proposed, in which space is entrained by matter expanding as a result of the Big Bang. In the absence of an input of external energy and momentum, the solution of the Einstein equations assumes a pressure equal to zero and an accelerated expansion. The cosmological perturbation model explains the difference in measurements of the Hubble constant using gravitational lensing and a distance ladder.

To accurately predict the matching relationships between the various components and the engine performance in the whole aero-engine environment, this study introduces a two-dimensional throughflow simulation method for the whole aero-engine. This method is based on individual throughflow solvers for the turbo-machinery and the combustor. It establishes a throughflow simulation model for the whole engine by integrating with the compressor-turbine co-operating equations and boundary conditions. The turbo-machinery throughflow solver employs a circumferentially averaged form of the time-dependent Navier-Stokes equations (N-S) as the governing equation. The combustor solver uses the Reynolds Average Navier-Stokes (RANS) method to solve flow and chemical reaction processes by constructing turbulence, combustion, and radiation models. The accuracy of the component solver is validated using Pratt &amp; Whitney’s three-stage axial compressor (P&amp;W3S1) and General Electric’s high-pressure turbine (GE-EEE HPT), and the predicted results are consistent with the experimental data. Finally, the developed throughflow method is applied to simulate the throttling characteristics of the WZ-X turboshaft engine. The results predicted by the throughflow program are consistent with the GasTurb calculations, including the trends of shaft power delivered, specific fuel consumption (SFC), inlet airflow, and total pressure ratio of the compressor. The developed method to perform throughflow simulation of the whole aero-engine eliminates the dependence on a general component map. It can quickly obtain the meridian flow field parameters and overall engine characteristics, which is expected to guide the design and modification of the engine in the future.

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

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

Exercising surds to represent figures is a common practice in scientific and Engineering fields, especially in scripts where calculators are banned or unapproachable. Peter Chew Theorem make result becomes simple when dealing with converting Quadratic Surds. The substance of the Peter Chew Theorem lies in enabling the forthcoming generation to simple break problems related to Quadratic Surds more effectively, easing a direct comparison with contemporary results. By employing the Peter Chew Theorem, one can streamline the tutoring and literacy of math, particularly concerning second- order direct equations with constant portions. This theorem's objective aligns with Albert Einstein's famed quotation Everything should be made as simple as possible, but not simpler.

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

This study aims to reveal the revolutionary Gamma Distribution Manifold (GDM) info-geometric analysis combined with potential Gamma distribution (GD) applications to advance Satellite Earth Observations (SEOs). An exposition of the GDM Fisher Information Matrix (FIM) was undertaken, together with the analytic investigation of the existence of FIM’s inverse. The derivation of GDM’s Geodesic Equations (GEs) is undertaken. The solution of the obtained Geodesic Equations is optimally complicated, and it provides more extensions to numerous open problems. Notably, the current paper provides a first-time unification between three multi-interdisciplinary fields, namely, information geometry theory, probability theory, and Satellite earth observations. Looking at the other side of the coin, the provided revolutionary approach in this current work reports the significant impact of info-geometric analysis on the probabilistic dynamics as well as the influential role of information geometry on satellite and space industry. This poses numerous open issues that initiates a new form of contemporary unification between these inter-disciplinaries to discover another visionary thinking out of traditional statistical frames. In other words, refraining ourselves from using classical statistical frames of thinking and delving into a new info-geometric vision of statistical distributions as well as examining real life applications of information geometry to advance satellite earth observations. The curtain is drawn on this paper by providing several challenging open problems in combination with a concluding remark and the next phase of research.

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

This paper brings to fruition recently found \cite{migdal2023exact} exact reduction of decaying turbulence in the Navier-Stokes equation in $3 + 1$ dimensions to a Number Theory problem of finding the statistical limit of the Euler ensemble. We reformulate the Euler ensemble as a Markov chain and show the equivalence of this formulation to the \textbf{quantum} statistical theory of $N$ fermions on a ring, with an external field related to the random fractions of $\pi$. We find the solution of this system in the turbulent limit $N\to \infty, \nu \to 0$ in terms of a complex trajectory (instanton) providing a saddle point to the path integral over the density of these fermions. This results in an analytic formula for the observable correlation function of vorticity in wavevector space. This is a full solution of decaying turbulence from the first principle without approximations or fitted parameters. The energy spectrum decays as $k^{-\frac{7}{2}}$. We compute resulting integrals in \Mathematica{} and present effective index $n(t) = -t\partial_t \log E $ for the energy decay as a function of time Fig.\ref{fig::NPlot}. The asymptotic value of the effective index in energy decay $n(\infty) = \frac{5}{4}$, but the universal function $n(t)$ is neither constant nor linear.

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.

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

In this paper, we introduce as a pedagogical strategy an internal division-free, straightforward, and symmetrically progressing algorithm in manually computing matrix inverse and solving systems of linear equations by revisiting the application of elementary row operations in the Gauss-Jordan reduction method and connecting it to the determinant method. The proposed cross-multiplication-flip (CMF) algorithm employs cross-multiplication similar to the butterfly movement in computing determinants as a strategic application of elementary row operations to efficiently reduce the rows and then applies flipping of rows and entries to put an upper triangular matrix into lower triangular form to continue the reduction process.

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.

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

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.

Partial Differential Equations (PDEs) constitute a fundamental framework for modeling various physical phenomena across diverse fields such as physics, engineering, and finance. Solving PDEs accurately and efficiently remains a significant challenge, particularly when dealing with complex systems or high-dimensional data. In this paper, we present a novel approach based on deep learning techniques to investigate the effective factors in solving PDEs. Our methodology leverages the expressive power of deep neural networks to learn and represent the underlying relationships between input parameters and the solutions of PDEs. By incorporating domain knowledge into the network architecture and training process, our approach aims to provide insights into the key factors influencing the solution of PDEs. We demonstrate the effectiveness of our method through numerical experiments on various PDEs, showcasing its capability to identify important features and improve predictive accuracy. Our findings suggest that deep learning techniques offer a promising avenue for understanding and analyzing the complex dynamics inherent in PDEs, paving the way for enhanced computational methods in scientific and engineering applications.

The article introduces the reader to the new thermodynamic potential of classical thermodynamics, called the J-potential. For the case of a homogeneous system, it vanishes, but it perfectly serves multiphase (heterogeneous) systems. The J-potential is especially important for the thermodynamics of capillary systems in colloidal science. The J-potential simplifies the derivation of the Neumann and Young equations, and has found application in the thermodynamics of thin films and the theory of strength of materials. In this work, for example, it was used to study the dependence of the chemical potential of a substance of a small bubble on the bubble size.

This paper presents a numerical method for solving a two-dimensional subdiffusion equation with a Caputo fractional derivative. The problem considered assumes symmetry in both the equation’s solution domain and the boundary conditions, allowing for a reduction of the two-dimensional equation to a one-dimensional one. The proposed method is an extension of the fractional Crank-Nicolson method, based on the discretization of the equivalent integral-differential equation. To validate the method, the obtained results were compared with a solution obtained through Laplace transform. The analytical solution in the image of the Laplace transform was inverted using the Gaver-Wynn-Rho algorithm implemented in the specialized mathematical computing environment, Wolfram Mathematica. The results clearly show the mutual convergence of the solutions obtained by the two methods.

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.

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

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

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

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

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

A formulation of nonequilibrium thermo field dynamics has been performed using the nonequilibrium statistical operator method by D.N. Zubarev. Generalized transfer equations for a consistent description of kinetics and hydrodynamics of the dense quantum-field system with strongly coupled states are derived.

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.

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

Recently, a paper was written to establish semi-structured complex numbers ℍ (a new number set invented to enable division by zero) as a useful mathematical tool providing novel results from its application in engineering and science. Whilst the paper was a milestone, the application of this number set to quantum and classical physics had not been fully explored. Consequently, the aim of this research was to use semi-structured complex numbers to develop a new mathematical model for the mass of an object and in course unite classical and quantum physics in a simple yet effective manner. As its major contributions this paper: (1) develops a new mathematical model for the mass of an object called the wave-mass function using De Broglie’s mass and frequency relation and semi-structured complex numbers; (2) uses the new wave-mass function to determine what happens to an object traveling at the speed of light; (3) derives a relation between Schrödinger’s wave function and the new wave-mass function; (4) derives a relation between Einstein’s Field Equation and the new wave-mass function; (5) uses the wave-mass function to combine Einstein’s Field Equation with Schrödinger’s wave function to form a new quantum gravity equation. These results and their applications provide a firm foundation to advance the number set ℍ as a useful mathematical tool.

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.

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.

This paper describes developing and analyzing an analytical model of Microelectromechanical Systems (MEMS) internal vibrating ring gyroscope. The internal ring structure consists of eight semicircular beams attached to the externally placed anchors. The paper analyses the vibrating ring gyroscope inplane displacement behavior and the resulting elliptical vibrational modes. The elliptical vibrational modes appear as pairs with the same resonance frequency due to the symmetric nature of the design. The analysis commences by conceptualizing the ring as a geometric structure with a circular shape, possessing specific dimensions such as thickness, height, and radius. We construct a linear model that characterizes the vibrational dynamics of the internal vibrating ring. The analysis encompasses comprehensive mathematical formulations for the radial and tangential displacements in local polar coordinates, considering the inextensional displacement of the ring structure. By utilizing the derived equations, we highlight the underlying relationships driving the vibrational characteristics of the MEMS vibrating ring gyroscope. These dynamic vibrational relationships are essential in enabling the vibrating ring gyroscope's future utilization in accurate navigation and motion sensing technologies.

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

The 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

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.

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.

In this paper, initially a mathematical model is formulated for transient frequency of power system considering time delays which occur while transmitting the control signals in open communication infrastructure. Time delay negligence in a power system leads to improper measurement of frequency variation in power system. The study of impact of time delays on the stability of power system is performed by estimating the decay rate of frequency wave form using Kalman Filter (KF). In power system, there is a possibility of multiple time delays. This paper also focusses on developing Interacting Multiple Model (IMM) Algorithm with multiple model space using Kalman Filter (KF) as state estimator tool. The multiple time delays in power system is considered as multiple model space. The result shows that KF provides better estimate of correct model for a particular input-set. The qualitative properties of Riccati difference equation (RDE) in terms of state error covariance of IMMKF are also analyzed and presented.

We primarily investigate the existence of solutions for fractional neutral integro-differential equations subjected to Neumann-type boundary conditions, which is crucial for understanding natural phenomena. Taking into account factors such as neutral type, fractional-order integrals, and fractional-order derivatives, we employ probability density functions, Laplace transforms, and resolvent operators to formulate a well-defined concept of a mild solution for the specified equation. Following this, by integrating fixed point theorems, we establish the existence of mild solutions under more relaxed conditions.

Our starting point is the covariant coordinate transformation equation of a relativistic positioning system in Minkowski space-time, which maps the receiver’s emission coordinates (proper times broadcast by the emitters) to its coordinates in some inertial reference frame. Bancroft’s analytical (closed-form) solution to the basic pseudorange navigation equations with four emitters is recovered and the subjacent elements are geometrically interpreted. The case of four static beacons is analized as a clarifying situation.

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.

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

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

The article proposes asymptotic methods for calculating Sommerfeld integrals, which enable us to calculate the integral using the expansion of a function into an infinite power series at the saddle point, where the role of a rapidly oscillating function under the integral can be fulfilled either by an exponential or by its product by the Hankel function. The proposed types of Sommerfeld integrals are generalized on the basis of integral representations of the Hertz radiator fields in the form of the inverse Hankel transform with the subsequent replacement of the Bessel function by the Hankel function. It is shown that the numerical values of the saddle point are complex. During integration, auxiliary or so-called standard integrals, which contain the main features of the integrand function, were used. As a demonstration of the accuracy of the technique, a previously known asymptotic formula for the Hankel functions was obtained in the form of an infinite series. The proposed method for calculating Sommerfeld integrals can be useful in solving the half-space Sommerfeld problem. An example is the expression in the form of an infinite series for the magnetic field of reflected waves, obtained directly through the Sommerfeld integral (SI).

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

Euler-Poisson equations of a charged symmetrical body in external constant and homogeneous electric and magnetic fields are deduced starting from the variational problem, where the body is considered as a system of charged point particles subject to holonomic constraints. The final equations are written for the center-of mass-coordinate, rotation matrix and angular velocity. General solution to the equations of motion is obtained for the case of a charged ball. For the case of a symmetrical charged body (solenoid), the task of obtaining the general solution is reduced to the problem of a one-dimensional cubic pseudo-oscillator. Besides, we present a one-parametric family of solutions to the problem in elementary functions.

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

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

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

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

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

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

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

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.

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.

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.

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.

This study suggests a novel cooperative multi-simplex algorithm that generalizes a local search optimizer to design a novel global search heuristic algorithm. The proposed algorithm exploits the vertex sharing strategy to enhance the search abilities of the working simplexes. The vertex sharing among the simplexes is carried out through cooperative step that is based on fitness of the underlying simplex. The proposed algorithm is applied to solve some systems of nonlinear equations by transforming them to optimization problems. Comparative analysis of results shows that the proposed method is practical and effective.

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

In this paper, existence and stability results of non-local integro-differential equation with Hilfer derivative of order ω∈(1,2) have been investigated. In this derivative is applied in RLC circuit equation of various fractional order. To proposed problem the existence solutions are derived using Schaefer’s fixed point theorem and the Banach contraction principle is using for uniqueness results. To demonstrate the effectiveness and applicability of our theoretical conclusions and two-step Lagrange polynomial interpolation were used to solve four numerical results.

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.

This paper discusses how mathematical models are used to predict environmental conditions. While these models are successful for larger scale phenomena, it becomes more complex for local weather or wave processes. Several approaches have been proposed to address this problem, including neural networks and statistical and mathematical models. The novelty in this approach is the use of statistical manifolds, which provide a more accurate way to estimate distances and cost functions. This work focuses on advancements in differential and information geometry for creating different types of models that can be used for environmental simulation and forecasting. These fields are vital for many institutes and centres worldwide, as they help solve challenges in areas such as meteorology, renewable energy, and aviation safety. Some potential IG applications to advance energy were given. Closing remarks with new open problems combined with the next research phase are provided.

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

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