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

o estudo visa explicar os efeitos diferentes entre as equações de multiplicação, soma, divisão e subtração. Nesse sentido, as equações de multiplicação e soma desenvolvem um novo elemento ou fato, enquanto isso, as equações de divisão expressar informações dentro do elemento ou fato, a subtração é a seleção de elemento ou fato Alem disso, para uma equação de multiplicação e soma resulta no novo elemento ou fato é necessario dependência de conceitos físicos específicos, sem alteração extrema, pois resultaria na alteração do elemento ou fato. vale ressaltar, a importância de valores fixos na equação para descreve a natureza específica do elemento ou fato. O estudo possibilita o entendimento das equações em aberta(sem solução) na matematica e fisica, além de compreender equações existentes. O estudo faz referência à teoria da necessidade obrigatória e teoria da diferenças entre elementos.

We present three reduced integrable hierarchies of nonlocal integrable NLS type equations from a vector integrable hierarchy associated with a matrix Lie algebra, not being A type. Three similarity transformations are taken to keep the invariance of the transformed zero curvature equations. The key step is to formulate a solution to a reduced stationary zero curvature equation so that the zero curvature formulation works for a reduced case.

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

For a linear ordinary differential equation (ODE in short) of the third order, results are presented that supplement the theory of conjugate points and extremal solutions by W. Leighton, Z. Nehari, M. Hanan. It is especially noted the sensitivity of solutions to the initial data, which makes their numerical study difficult. Similar results were obtained for the third-order nonlinear equations of the Emden-Fowler type.

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

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

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

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

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

The creep response of commercially pure titanium was investigated in air at 550, 600 and 650°C, to assess the effect of oxidation on the mechanical response. Experiments demonstrated that exposure at high temperature produced a marked reduction of the minimum creep rate under a given applied stress. Microhardness measurements showed that a hardened zone forms in proximity of the surface, due to oxygen penetration into metal. A simplified composite-model was then used to describe the creep response. In this model, the sample consisted of two zones, the hard case, enriched in oxygen, and the soft pure-titanium core, both creeping with similar strain rates. Calculation led to an estimation of the dependence of the minimum creep rate on stress and temperature for the hard high-oxygen zone. The simplified composite-model here presented gave an excellent description of experimental creep data for pure titanium tested in air and also provided a reliable picture of the effect of oxidation on complex Ti-alloys.

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

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

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

Multi-parameter families of Lax pairs for the modified Korteweg-de Vries (mKdV) equation are defined by applying a direct method developed in the present study. The gauge transformations converting the defined Lax pairs to some simpler forms and reducing the variety of the Lax pairs this way are found. The direct method and its possible applications to other types of evolution equations are discussed.

Antiangiogenic therapy (AAT) is an indirect oncological modality which is aimed at disruption of cancer cells nutrient supply. Invasive tumors have been shown to possess inherent resistant to this treatment, while compactly growing benign tumors react to it by shrinkage. It is generally accepted that AAT by itself is not curative. This study presents a mathematical model of non-invasive tumor growth with physiologically justified account of alteration of microvasculature and biomechanical aspects during tumor growth and AAT. In untreated setting the model reproduces tumor growth with saturation, where the maximum tumor volume depends on the level of angiogenesis. The outcomes of simulations of AAT depend on the tumor size at the moment of treatment initiation. If it is close to the stable size of avascular tumor grown in absence of angiogenesis, then the tumor is rapidly stabilized by AAT. Treatment of large tumors is accompanied by displacement of normal tissue due to the tumor shrinkage. During it, microvasculature undergoes distortion which degree depends on the displacement distance. As it affects tumor nutrient supply, the stable size of a tumor that undergoes AAT negatively correlates with its size at the beginning of treatment. For sufficiently large initial tumors, the long-term survival of tumor cells is compromised by competition with normal cells for severely limited inflow of nutrients, which makes AAT effectively curative.

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

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

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

Fokker-Planck PDEs (incl. diffusions) for stable Lévy processes (incl. Wiener processes) on the joint space of positions and orientations play a major role in mechanics, robotics, image analysis, directional statistics and probability theory. Exact analytic designs and solutions are known in the 2D case, where they have been obtained using Fourier transform on $SE(2)$. Here we extend these approaches to 3D using Fourier transform on the Lie group $SE(3)$ of rigid body motions. More precisely, we define the homogeneous space of 3D positions and orientations $\mathbb{R}^{3}\rtimes S^{2}:=SE(3)/(\{\mathbf{0}\} \times SO(2))$ as the quotient in $SE(3)$. In our construction, two group elements are equivalent if they are equal up to a rotation around the reference axis. On this quotient we design a specific Fourier transform. We apply this Fourier transform to derive new exact solutions to Fokker-Planck PDEs of $\alpha$-stable Lévy processes on $\mathbb{R}^{3}\rtimes S^{2}$. This reduces classical analysis computations and provides an explicit algebraic spectral decomposition of the solutions. We compare the exact probability kernel for $\alpha = 1$ (the diffusion kernel) to the kernel for $\alpha=\frac12$ (the Poisson kernel). We set up SDEs for the Lévy processes on the quotient and derive corresponding Monte-Carlo methods. We verify that the exact probability kernels arise as the limit of the Monte-Carlo approximations.

This experimental research presents an easier model for fitting experimental results that can be simulate the dynamical behaviour of a single degree freedom system, SDFS, subjected in a forced vibration in a shaking table. The fractional damping using a Kelvin-Voigt fractional derivative rheological model, KVFRdm, to represent a fractional oscillator and compare with the experimental model.

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

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

Electrodynamics is usually written using polarization fields to describe changes in distribution of charge as electric fields change. This approach does not specify polarization fields uniquely from electrical measurements. Many polarization fields will produce the same electrodynamic forces and flows because only divergence of polarization enters Maxwell’s first equation, relating charge and electric field. The curl of any function can be added to a polarization field without changing the electric field at all. The divergence of the curl is always zero. To be unique, models must describe the charge distribution and how it varies. I propose a different paradigm to describe field dependent charge, i.e., the phenomenon of polarization. This operational definition of polarization has worked well in biophysics for fifty years, where a field dependent, time dependent polarization provides gating current that makes neurons respond sensitively to voltage. Theoretical estimates of polarization computed with this definition fit experimental data. I propose that operational definition be used to define polarization charge in general. Charge movement needs to be computed from a combination of electrodynamics and mechanics because ‘everything interacts with everything else’. The classical polarization field need not enter into that treatment at all. When nothing is known about polarization, it is necessary to use an approximate representation with a dielectric constant that is a single real positive number. This approximation allows important results in some cases, e.g., design of integrated circuits in silicon semiconductors, but can be seriously misleading in other cases, e.g., ionic solutions.

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

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

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

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

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

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

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

This short paper briefly explores the interplay between quantum mechanics and gravity, drawing from the principle of gravitational redshift and bringing it into quantum mechanical frameworks. The exploration looks through quantum mechanical relations and equations such as Planck-Einstein relation, De-Broglie’s relation and the relativistic wave equations, and propagators, accounting for the presence of gravitational fields through gravitational redshift of its momentum.

The susceptible-infected-recovered-vaccinated-deceased (SIRVD) epidemic compartment model extends the SIR model to include the effects of vaccination campaigns and time-dependent fatality rates on epidemic outbreaks. It encompasses the SIR, SIRV, SIRD, and SI models as special cases, with individual time-dependent rates governing transitions between different fractions. We investigate a special class of exact solutions and accurate analytical approximations for the SIRVD and SIRD compartment models. While the SIRVD and SIRD equations pose complex integro-differential equations for the rate of new infections and the fractions as a function of time, a simpler approach considers determining equations for the sum of ratios for given variations. This approach enables us to derive fully exact analytical solutions for the SIRVD and SIRD models. For nonlinear models with a high-dimensional parameter space, such as the SIRVD and SIRD models, analytical solutions, exact or accurately approximative, are of high importance and interest: not only as suitable benchmarks for numerical codes, but especially as they allow us to understand the critical behavior of epidemic outbursts as well as the decisive role of certain parameters. In the second part of our study, we apply a recently developed analytical approximation for the SIR and SIRV models to the more general SIRVD model. This approximation offers accurate analytical expressions for epidemic quantities, such as the rate of new infections and the fraction of infected persons, particularly when the cumulative fraction of infections is small. The distinction between recovered and deceased individuals in the SIRVD model affects the calculation of the death rate, which is proportional to the infected fraction in the SIRVD/SIRD cases but often proportional to the rate of new infections in many SIR models using a posteriori approach. We demonstrate that the temporal dependence of the infected fraction and the rate of new infections differs when considering the effects of vaccinations and when the real-time dependence of fatality and recovery rates diverge. These differences are highlighted for stationary ratios and gradually decreasing fatality rates. \Revised{The case of stationary ratios allows one to construct a new powerful diagnostics method to extract analytically all SIRVD model parameters from measured COVID-19 data of a completed pandemic wave.

We consider two abstract systems of parabolic-hyperbolic type that model thermoelastic problems. We study the influence of the physical constants and the initial data on the nonexistence of global solutions that in our framework are produced by the blow-up in finite time of the norm of the solution in the phase space. We employ a differential inequality to find sufficient conditions that produce the blow-up. To that end, we construct a set that is positive invariant for any positive value of the initial energy. As a result we found that the coupling with the parabolic equation stabilizes the system, as well as the damping term in the hyperbolic equation. Moreover, for any pair of positive values $(\xi,\epsilon)$, there exist initial data such that the corresponding solution with initial energy $\xi$ blows-up at a finite time less than $\epsilon$. Our purpose is to improve results previously published in the literature.

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

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

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

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

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

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

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

A mathematical model, composed of two non-linear differential equations that describe the population dynamics of CD4 T cells in the human immune system, as well as viral HIV particles, is proposed. The invariance region is determined, classical equilibria stability analysis is performed using the basic reproduction number, and numerical simulations are carried out, in order to illustrate stability results. Later, the model is modified with a delay term, which describes the time that cells require for immunological activation. This generates a two-dimensional integro-differential system, which is transformed into a system with three ordinary differential equations, via auxiliary variable use. For the new model, equilibrium points are determined, their local stability is examined, and results are studied by way of numerical simulation.

Periodic travelling waves are observed in various brain activities including visual, motor, language, sleep, and so on. There are several neural field models describing periodic waves assuming nonlocal interaction and, possibly, inhibition, time delay, or some other properties. In this work we study the influence of asymmetric connectivity functions and of time delay on the emergence of periodic waves and on their properties. Nonlinear wave dynamics is studied, including modulated and aperiodic waves. Multiplicity of waves for the same values of parameters is observed. External stimulation in order to restore wave propagation in a damaged tissue is discussed.

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

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

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

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

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

Use of MS Excel as a platform for teaching simple iterative solution of partial differential equations runs into issues of reliability of the computation because of the large number of calculations required. This article documents the problem and an empirical solution. The cause of this unreliability was found by detailed examination of the numerical results from three simultaneous iterations. Computers of several types from different manufacturers were found to fail occasionally in updating variables in the numerical iteration. A simple work-around is suggested and was tested on several platforms.

This paper considers the Hermitian solutions of a new system of commutative quaternion matrix equations, we establish both necessary and sufficient conditions for the existence of solutions. Furthermore, we derive an explicit general expression when it is solvable. To illustrate our main findings, we present two numerical algorithms and examples in this paper.

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.

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

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

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.

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

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.

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

This article deals with special curves similar to those known in mathematics as Lamé curves. Their equations have a certain similarity with the equations of Lamé curves. But these curves are more like ellipses in appearance. Like ellipses, they have major and minor axes. These are such curves that their major and minor semi-axis coincide with an ellipse but are different than an ellipse. The main feature of these curves is that for any triangle whose base coincides with the axis of a given curve and the vertex is on a given curve, the condition aⁿ+bⁿ=cⁿ is obtained. The number “n” is called the degree or rank of this curve. To distinguish this group of curves, we called them Elbe curves. The equations of these lines are also provided in the paper.

This article describes an approach for optimising sensor-controlled systems through minimal intervention, utilising Fuzzy Linear Systems of Equations (FLSE). Staring with a generalised model of the system behaviour, incorporating an array of control units, environmental sensors, and an expert knowledge base. The described problems of detecting the level of intervention needed to change the system state to another is handled with the help of developed methods for solving the inverse problem for FLSE. By achieving minimal intervention, we ensure that system adjustments effective, economically optimal and non-intrusive. A MATLAB-based implementation is presented.

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 dynamics of the main photosynthetic structures is studied by mathematical modeling methods in this article. Chlorophyll portion variability in phytoplankton and the formation of energy-intensive substances in the process of photosynthesis underlie the models. These cellular components are considered in terms of their participation in the growth of specific biomass. Computational experiments are presented that simulate various degrees of environmental friendliness. The corresponding functions are built in accordance with seasonal fluctuations throughout the year in the Far East region of Russia. The stability of model solutions in long-term dynamics is also investigated. The models were tested for biological adequacy and their effectiveness was compared. For the model selected as a result of the comparison, the optimal control problem has been formulated and solved. It is shown that in this way it is possible to reduce the space of the initial components of the model system.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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

RNA transcripts play a crucial role as witnesses of gene expression health. Identifying disruptive short sequences in RNA transcription and regulation is essential for potentially treating diseases. Let’s delve into the mathematical intricacies of these sequences. We’ve previously devised a mathematical approach for defining a “healthy” sequence (10.14218/GE.2023.00079). This sequence is characterized by having at most four distinct nucleotides (denoted as nt≤4). It serves as the generator of a group denoted as fp. The desired properties of this sequence are as follows: fp should be close to a free group of rank nt−1, it must be aperiodic, fp should not have isolated singularities within its SL2(C) character variety (specifically within the corresponding Groebner basis). Now, let’s explore the concept of singularities. There are cubic surfaces associated with the character variety of a four-punctured sphere denoted as S24. When we encounter these singularities, we find ourselves dealing with some algebraic solutions of a dynamical second-order differential (and transcendental) equation known as the Painlevé VI equation (10.3390/dynamics4010001). In certain cases, S24 degenerates, in the sense that two punctures collapse, resulting in a “wild” dynamics governed by the Painlevé equations of an index lower than VI. In our paper, we provide examples of these fascinating mathematical structures within the context of miRNAs. Specifically, we find a clear relationship between decorated character varieties of Painlevé equations and the character variety calculated from the seed of oncomirs. These findings should find applications in cancer research.

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.

The work is devoted to the development of methods for constructing asymptotic formulas as x→∞ of a fundamental system of solutions of linear differential equations generated by a symmetric two-term differential expression of odd order. The coefficients of the differential expression belong to classes of functions that allow oscillation (for example, those that do not satisfy the classical Titchmarsh-Levitan regularity conditions). As a model equation, the 5th order equation i2p(x)y′′′′′+p(x)y′′′′′+q(x)y=λy, for which various cases of behavior of the coefficients p(x),q(x), is investigated. New asymptotic formulas are obtained for the case when the function h(x)=−1+p−1/2(x)∉L1[1,∞) significantly influences the asymptotics of solutions to the equation. The case when the equation contains a nontrivial bifurcation parameter is studied.

Applications of gradient method for nonlinear optimization in development of Gradient Neural Network (GNN) and Zhang Neural Network (ZNN) are investigated. Particularly, the solution of the matrix equation AXB=D which changes over time is studied using the novel GNN model, termed as GGNN(A,B,D). The GGNN model is developed applying GNN dynamics on the gradient of the error matrix used in the development of the GNN model. The convergence analysis shows that the neural state matrix of the GGNN(A,B,D) design converges asymptotically to the solution of the matrix equation AXB=D, for any initial state matrix. It is also shown that the convergence result is the least square solution which is defined depending on the selected initial matrix. A hybridization of GGNN with analogous modification GZNN of the ZNN dynamics is considered. The Simulink implementation of presented GGNN models is carried out on the set of real matrices.

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

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