The paper is considering of the basic differential equations of hydrodynamics: the equation of continuity and motion. On the simplest example, it is shown that the equation of continuity in a system with the equation of motion leads to contradictions and erroneous results of modeling. A more correct form of the continuity equation is described. It is shown that the equations of motion can be written in the form of complete differentials. Three possible integral forms of the equations of motion are presented. As a conclusion, the existence and smoothness of the solution of the Navier-Stokes equations are considered.

This paper attempts to clarify an long-standing issue about the number of unknowns in the Reynolds-Averaged Navier-Stokes equations (RANS). This study shows that all perspectives regarding the numbers of unknowns in the RANS stem from the misinterpretation of the Reynolds stress tensor. The current literature consider that the Reynolds stress tensor has six unknown components; however, this study shows that the Reynolds stress tensor actually has only three unknown components, namely the three components of fluctuation velocity. This understanding might shed a light to understand the well-known closure problem of turbulence.

This paper showed that turbulence closure problem is not an issue at all. All mistakes in the literature regarding the numbers of unknown quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor, i.e., all literature has stated that the symmetric Reynolds stress tensor has six unknowns; however, it actually has only three unknowns, i.e., the three components of fluctuation velocity. We showed the integral-differential equations of the Reynolds mean and fluctuation equations have exactly eight equations, which equal to the numbers of quantities in total, namely, three components of mean velocity, three components of fluctuation velocity, one mean pressure and one fluctuation pressure. With this understanding, the closed Reynolds Navier-Stokes turbulence equations of incompressible flows were formulated. This study may help to solve the puzzle that has eluded scientists and mathematicians for centuries.

The purpose of this paper is to bring to light a method through which the global in time existence for arbitrary large in H¹ initial data of a strong solution to 3D periodic Navier-Stokes equations follows. The method consists of subdividing the time interval of existence into smaller sub-intervals carefully chosen. These sub-intervals are chosen based on the hypothesis that for any wavenumber m, one can find an interval of time on which the energy quantized in low-frequency components (up to m) of the solution u is lesser than the energy quantized in high-frequency components (down to m) or otherwise the opposite. We associate then a suitable number m to each one of the intervals and we prove that the norm ||u(t)||_H1 is bounded in both mentioned cases. The process can be continued until reaching the maximal time of existence T_max which yields the global in time existence of strong solution.

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, for the Reynolds-averaged Navier-Stokes equations, a self-closed turbulence model without any adjustable parameter is formulated and a simple model is proposed. The validation rule for self-closed turbulence model is rigorously derived from the Reynolds-averaged Navier-Stokes equation. The rule is not effected by turbulence modelling on the Reynolds stresses. As an application, two dimensional turbulent boundary layers, and plate turbulent boundary layers have been formulated.

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.

A method of solution to solve the compressible unsteady 3D Navier-Stokes Equations in cylindrical co-ordinates coupled to the continuity equation in cylindrical coordinates is presented in terms of an additive solution of the three principle directions in the radial, azimuthal and z directions of flow. A dimensionless parameter is introduced whereby in the large limit case a method of solution is sought for in the boundary layer of the tube. A reduction to a single partial differential equation is possible and integral calculus methods are applied for the case of a body force directed to the centre of the tube to obtain an integral form of the Hunter-Saxton equation. Also an extension for a more general body force is shown where in addition there is a rotational force applied.

Based on author's previous work [Sun, B. The Reynolds Navier-Stokes Turbulence Equations of Incompressible Flow Are Closed Rather Than Unclosed. Preprints 2018, 2018060461 (doi: 10.20944/preprints201806.0461.v1)], this paper proposed an explicit representation of velocity fluctuation and formulated the Reynolds stress tensor in terms of the mean velocity field. The proposed closed Reynolds Navier-Stokes turbulence formulations reveal that the mean vorticity is the key source of producing turbulence.

This paper shown that turbulence closure problem is not an issue at all. All mistakes in the literature regarding the numbers of unknown quantities in the Reynolds turbulence equations stem from the misunderstandings of physics of the Reynolds stress tensor, i.e., all literatures have stated that the symmetric Reynolds stress tensor has six unknowns; however, it actually has only three unknowns, i.e., the three components of fluctuation velocity. We shown the integral-differential equations of the Reynolds mean and fluctuation equations have exactly eight equations, which equal to the numbers of quantities in total, namely, three components of mean velocity, three components of fluctuation velocity, one mean pressure and one fluctuation pressure. That is why we claim in this paper, that the Reynolds Navier-Stokes turbulence equations of incompressible flow are closed rather than unclosed. This study may help to solve the puzzle that has eluded scientists and mathematicians for centuries.

We studies the initial value problem for the fractional Navier-Stokes-Coriolis equations, which obtained by replacing the Laplacian operator in the Navier-Stokes-Coriolis equation by the more general operator $(-\Delta)^\alpha$ with $\alpha>0$. We introduce function spaces of the Besove type characterized by the time evolution semigroup associated with the general linear Stokes-Coriolis operator. Next, we establish the unique existence of global in time mild solutions for small initial data belonging to our function spaces characterized by semigroups in both the scaling subcritical and critical settings.

Into the frame of the French TANDEM project (Tsunamis in the Atlantic and the English ChaNnel: Definition of the Effects through numerical Modelling) Principia has been working on the development and qualification of two in-house CFD software’s: the 2D EOLE-SV (Saint Venant) model for simulation of large scale tsunami propagation from the source up to coastal scale and the 3D EOLE-NS (Navier-Stokes) model dedicated to tsunami coastal impact modelling. This paper presents a large range of test cases carried out into the frame of the project and dedicated to the validation of numerical codes in various tsunami wave conditions. The main aspects of phenomena such as wave generation, propagation and coastal impact are investigated on academic situations. A real case simulation is concerned as well, the devastating 2011 Tohoku event which is compared with in-situ data.

We address the dynamics of two-dimensional Navier-Stokes models with infinite delay and hereditary memory, whose kernels are a much larger class of functions than the one considered in the literature, on a bounded domain. We prove the existence and uniqueness of weak solutions by means of Faedo-Galerkin method. Moreover, we establish the existence of global attractor for the system with the existence of a bounded absorbing set and asymptotic compact property.

\abstract{The existence of vorticies is well known in fluid mechanics and particularly in turbulence theory. In the context of the millennium problem of the Clay Institute on the existence (or not) of Global solutions in time of the Periodic Navier Stokes equations (PNS) results are now proposed to confirm that Global solutions do not exist in the general case of PNS defined on $\mathbb{T}^3$. This corresponds to part ($D$) proposition in the Navier Stokes problem as posted by C.Fefferman. Recent work has shown that for the Euler equations defined on the 3-Torus there exists a H\"older continuous solution for the third component of velocity $v_3$ at the center of each cell or cube of the 3-Torus Lattice. It has also been determined that as the kinematic viscosity changes from the corresponding value at $\nu=0$ to the fully viscous case for the Periodic Navier Stokes equations (PNS) at the value $\nu=1$ the solution reaches a peak ($\frac{\partial v_3}{\partial y_3} = 0$ ) in $y_3$. (third component of $\vec{r}=(y_1,y_2,y_3)$ in Cartesian co-ordinates) Here in the viscous case there exists a dipole which is not centered at the center of the cell of the Lattice. This immediately implies that since the dipole by definition has an equal in magnitude positive and negative peak in the third component of velocity, then the dipole Riemann cut-off surface is covered by a closed surface which is the sphere $\|\vec{r}\|=1$ and where a given cell of dimensions $[-1,1]^3$ is circumscribed on a sphere of radius $1$. For such a closed surface containing a dipole it necessarily follows that the flux at the surface of the sphere of $v_3$ wrt to surface normal $\vec{n}$ is zero including at $y_i, i=1\dots 3$ is zero at the points where the surface of sphere touches the cube walls. At the finite time singularity on the sphere a rotation boundary condition is deduced. Here $\frac{\partial v_3}{\partial \vec{n}}= \vec{n}\cdot \nabla v_3$. Also a boundary condition on the sphere shows rotation of sphere. Finally as a result, on the sphere a solution for $v_3$ is obtained which is proven to be H\"older continuous and it is proven that it is possible to extend H\"older continuity on the sphere uniquely to all of the interior of the ball.}

This paper proposed an explicit and simple representation of velocity fluctuation and the Reynolds stress tensor in terms of the mean velocity field. The proposed turbulence equations are closed. The proposed formulations reveal that the mean vorticity is the key source of producing turbulence. It is found that there are no velocity fluctuation and turbulence if there were no vorticity. As a natural consequence, the laminar- turbulence transition condition was obtained in a rational way.

The problem of fluid dynamics can be greatly simplified if, for every point in space, the strain-rate tensor is diagonalized. This tensor is introduced into the Navier-Stokes equations via material law and divergence of the stress tensor. This article shows that local SO(3)xU(1) gauge fields can be used to locally diagonalize the diffusion components of the strain-rate tensor. The gauge fields resulting from the connection can be interpreted as convection components of the flow, they show properties of quasiparticles and can be interpreted as elementary vortices. Thus, the proposed approach not only offers new insights for the solution and situative simplification of the Navier-Stokes equations, it also uncovers hidden symmetries within the flow convection, allowing - depending on boundary conditions - further interpretation.

In this paper, we establish the existence and uniqueness results of the fractional Navier-Stokes evolution equation using the Banach fixed point theorem, where the fractional order β is in the form of the Atangana-Baleanu Caputo fractional order. The iterative method combined with the Laplace transform and Sumudu transform are employed to find the exact and approximate solutions of the fractional Navier-Stokes equation of a one-dimensional problem of unsteady flow of a viscous fluid in a tube. Mathematica software is used to present the solutions graphically.

For the solution existence condition of the Navier-Stokes equation, we propose a conjecture as follows: "\emph{The Navier-Stokes equation has a solution if and only if the determinant of flow velocity gradient is not zero, namely $\det (\bm \nabla \bm v)\neq 0$.}"

A careful examination of the fundamentals of electromagnetic theory shows that due to the underlying mathematical assumptions required for Stokes' Theorem, global charge conservation cannot be guaranteed in topologically non-trivial spacetimes. However, in order to break the charge conservation mechanism we must also allow the electromagnetic excitation fields D, H to possess a gauge freedom, just as the electromagnetic scalar and vector potentials $\varphi$ and A do. This has implications for the treatment of electromagnetism in spacetimes where black holes both form and then evaporate, as well as extending the possibilities for treating vacuum polarisation. Using this gauge freedom of D, H we also propose an alternative to the accepted notion that a charge passing through a wormhole necessarily leads to an additional (effective) charge on the wormhole's mouth.

The main object of this work is to establish the existence and uniqueness of a solution to the 3D Navier-Stokes (NS) system for an incompressible fluid with viscosity. The nonlinearity of the NS system, as well as the need to estimate velocity and pressure for every value of the viscosity parameter [1] make them challenging to solve. In this regard, in the present work, we study the Navier-Stokes system, which describe the flow of a viscous incompressible fluid and the solution was obtained for velocity and pressure in an analytical form. In addition, the found pressure distribution law, which is described by a Poisson type equation and plays a fundamental role in the theory of Navier-Stokes systems in constructing analytic smooth (conditionally smooth) solutions.

This article develops applications of the generalized method of lines to numerical solutions of the time-independent, incompressible Navier-Stokes system in fluid mechanics. We recall that for such a method, the domain of the partial differential equation in question is discretized in lines (or more generally in curves), and the concerning solutions are written on these lines as functions of the boundary conditions and the domain boundary shape.

An open source computational fluid dynamics (CFD) solver has been incorporated into the WindNinja modeling framework widely used by wildland fire managers as well as researchers and practitioners in other fields, such as wind energy, wind erosion, and search and rescue. Here we describe incorporation of the CFD solver and evaluate its performance compared to the conservation of mass (COM) solver in WindNinja and previously published large-eddy simulations (LES) for three field campaigns conducted over isolated terrain obstacles of varying terrain complexity: Askervein Hill, Bolund Hill, and Big Southern Butte. We also compare the effects of two important model settings in the CFD solver and provide guidance on model sensitivity to these settings. Additionally, we investigate the computational mesh and difficulties regarding terrain representation. Two important findings from this work are: (1) the choice of discretization scheme for advection has a significantly larger effect on the simulated winds than the choice of turbulence model and (2) CFD solver predictions are significantly better than the COM solver predictions at windward and lee side observation locations, but no difference was found in predicted speed-up at ridgetop locations between the two solvers.

Collisional mountain belts commonly develop intramontane basins from mechanical and isostatic subsidence during orogenic development. These frequently display a relict top surface, evidencing a change interval from basin infilling to erosion often via capture or overspill. Such surfaces provide markers that inform on orogenic growth patterns via climate and base level interplay. Here, we describe the top surface from the Sorbas Basin, a key intramontane basin within the Betic Cordillera (SE Spain). The surface is fragmentary comprising high elevation hilltops and discontinuous ridges developed onto the variably deformed final basin infill outcrop (Gochar Formation). We reconstruct surface configuration using DEM interpolation and apply ¹⁰Be/²⁶Al cosmonuclides to assess surface formation timing. The surface is a degraded Early Pleistocene erosional pediment developed via autogenic switching of alluvial fan streams under stable dryland climate and base level conditions. Base level lowering since the Middle Pleistocene focused headwards incision up interfan drainages, culminating in fan head capture and fan morphological preservation within the abandoned surface. Post abandonment erosion has lowered the basin surface by 31 m (average) and removed ~5.95 km³ of fill. Regional basin comparisons reveal a phase of Early Pleistocene surface formation, marking landscape stability following the most recent Pliocene-Early Pleistocene mountain building. Post-surface erosion rate quantification is low and in accordance with ¹⁰Be denudation rates typical of the low uplift Betic Cordillera.

Two-dimensional periodic Stokes flow of viscous incompressible fluid in a rectangular cavity with constant velocities applied to the top and bottom walls is considered. The study of chaotic advection regimes is reduced to the sequential solving of two problems. To solve the first problem, the analytical method of superposition, that allows to obtain any desired accuracy of the velocity field, is used. To solve the second problem, which is associated with obtaining the trajectories of individual fluid particles, numerical calculations of the Cauchy problem were performed. An analysis of the boundary conditions accuracy was performed based on the control of local integration. The advection of the selected volume of fluid in the rectangular cavity under the periodic motion of walls for a finite period time based on piecewise spline cubic interpolation is modelled. The obtained numerical results agree with well-known experimental data.

Unwanted proteins and metabolic waste in cerebral spinal fluid are cleared from the brain by meningeal and nasal lymphatics, and perineural sheath of cranial nerves; however, the distribution and clearance of CSF along the subarachnoid space of the entire spinal cord is not fully understood. It was hypothesized that the anatomical resolution of Cryo-Fluorescence Tomography (CFT) could provide the visual evidence of clearance across the spinal cord at each level of the vertebral column. To that end, isoflurane anesthetized mice were infused into the lateral cerebroventricle with 5.0 µl of quantum dots [QdotR 605 ITKTM amino (PEG)] over two mins. Mice were allowed to recover (ca 2-3 min) and remained awake and ambulatory for 5, 15, 30, 60 and 120 minutes after which they were euthanized, and the entire intact body frozen at -800. The entire mouse was sectioned and imaged to produce an isotropic voxel resolution of 35 µm. Both white light and fluorescent images were captured after each slice to produce high resolution three-dimensional volumes for each mouse. CFT highlighted the circulation of tracer throughout the ventricular system and central canal of the spinal cord and the entire subarachnoid space of the CNS. Signal could be visualized in the nasal cavity, deep cervical lymph nodes, thoracic lymph nodes and more superficial submandibular lymph nodes as early as 15 min post infusion. Fluorescent signal could be visualized along the dorsal root ganglia and down the proximal extension of the spinal nerves of the thoracic and lumbar segments at 30 min. There was significant accumulation of tracer in the lumbar and sacral lymph nodes between 15- 60 min. The dense fluorescent signal in the thoracic vertebrae noted at 5- and 15-mins post infusion was significantly reduced by 30 min. Indeed, all signal in the spinal cord was ostensibly absent by 120 min except for trace amounts in the coccyx. The brain still had some residual signal at 120 min. These data show that Qdots with a hydrodynamic diameter of 16-20 nm rapidly clear from the brain of awake mice. These data also clearly demonstrate the rapid distribution and efflux of trace along a major length of the vertebral column and the potential contribution of the spinal cord in the clearance of brain waste.

In this paper, we extend the KdVSKR equations in (1+1)-dimension and (2+1)-dimension to fractional KdVSKR equations with modified Riemann-Liouville derivative. The (2+1)-dimensional KdVSKR equation, which is a recent extension of (1+1)-dimensional KdVSKR equation, can model the resonances of solitons in shallow water. By means of the Hirota bilinear method, finite symmetry group method and consistent Riccati expansion method, many new interaction solutions have been derived. Soliton-cnoidal interaction solution of the (1+1)-dimensional fractional KdVSKR equation has been derived for the first time. For the (2+1)-dimensional frac-tional KdVSKR equation, two-wave interaction solutions and three-wave interaction solutions, including dark-soliton-sine interaction solution, bright-soliton-elliptic interaction solution, and lump-hyperbolic-sine interaction solution, have been derived. Impact of the fractional order on the shapes of these solutions has been illustrated by figures. Three-wave interaction solutions of fractional systems have not been reported in the existing references. The research idea in this paper can be applied to other fractional differential equations.

The COVID-19 pandemic underscores the significance of vaccine hesitancy in shaping vaccination outcomes. Understanding the factors underpinning COVID-19 vaccination hesitancy is crucial for tailoring effective vaccination strategies. This cross-sectional study, conducted in three communities across the United States and Lebanon, employed surveys to assess respondents' knowledge, attitudes, and perceptions regarding COVID-19 infection and vaccination. Among the 7,196 participants, comprising 6,775 from the US and 422 from Lebanon, vaccine hesitancy rates were comparable at 12.2% and 12.8%, respectively. Notably, a substantial proportion of respondents harbored misconceptions, such as attributing the potential to alter DNA (86.4%) or track individuals (92.8%) to COVID-19 vaccines, and believing in the virus's artificial origins (81%). Primary determinants of hesitancy included perceptions that the vaccine poses greater risk than the infection itself (aOR = 8.7 and 9.4, respectively) and negative recommendations from healthcare providers (aOR = 6.5 and 5.4, respectively). Conversely, positive endorsements from healthcare providers were associated with reduced hesitancy (aOR = 0.02 and 0.4, respectively). Targeting healthcare providers to dispel misinformation and elucidate COVID-19 vaccine risks holds promise for enhancing vaccination uptake.

The lab-on-a-chip concept, enabled by microfluidic technology, promises the integration of multi-ple discrete laboratory techniques into a miniaturised system. Research into microfluidics has generally focused on the development of individual elements of the total system (often with rela-tively limited functionality), without full consideration for integration into a complete fully opti-mised and miniaturised system. Typically, the operation of many of the reported lab-on-a-chip devices is dependent on the support of a laboratory framework. In this paper, a demonstrator platform for routine laboratory analysis is designed and built, which fully integrates a number of technologies into a single device with multiple domains such as fluidics, electronics, pneumatics, hydraulics, and photonics. This facilitates the delivery of breakthroughs in research, by incorpo-rating all physical requirements into a single device. To highlight this proposed approach, this demonstrator microsystem acts as a fully integrated biochemical assay reaction system. The re-sulting design determines enzyme kinetics in an automated process and combines reservoirs, three-dimensional fluidic channels, optical sensing, and electronics in a low-cost, low-power and portable package.

The oil and gas industry faces increasing pressure to remove people from dangerous offshore environments. Robots present a cost-effective and safe method for inspection, repair and maintenance of topside and marine offshore infrastructure. In this work, we introduce a new immobile multi-sensing robot, the Limpet, which is designed to be low-cost and highly manufacturable, and thus can be deployed in huge collectives for monitoring offshore platforms. The Limpet can be considered an instrument, where in abstract terms, an instrument is a device that transforms a physical variable of interest (measurand) into a form that is suitable for recording (measurement). The Limpet is designed to be part of the ORCA (Offshore Robotics for Certification of Assets) Hub System, which consists of the offshore assets and all the robots (UAVs, drones, mobile legged robots etc.) interacting with them. The Limpet comprises the sensing aspect of the ORCA Hub System. We integrated the Limpet with Robot Operating System (ROS), which allows it to interact with other robots in the ORCA Hub System. In this work, we demonstrate how the Limpet can be used to achieve real-time condition monitoring for offshore structures, by combining remote sensing with signal processing techniques. We show an example of this approach for monitoring offshore wind turbines. We demonstrate the use of four different communication systems (WiFi, serial, LoRa and optical communication) for the condition monitoring process. By processing the sensor data on-board, we reduce the information density of our transmissions, which allows us to substitute short-range high-bandwidth communication systems with low-bandwidth long-range communication systems. We train our classifier offline and transfer its parameters to the Limpet for online classification, where it makes an autonomous decision based on the condition of the monitored structure.

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

This article aims to explore the existence and stability of solutions to differential equations involving a regularized ψ−Hilfer fractional derivative, which, compared to standard ψ−Hilfer fractional derivatives, provide a clear physical interpretation when dealing with initial conditions. We discovered that the regularized ψ−Hilfer fractional derivative can be represented as the inverse operation of the ψ−Riemann–Liouville fractional integral, and used this property to prove the existence of solutions for linear differential equations with a regularized ψ−Hilfer fractional derivative. Additionally, we applied the Mönch’s fixed-point theorem and knowledge of non-compactness measures to demonstrate the existence of solutions for nonlinear differential equations with a regularized ψ−Hilfer fractional derivative, and further discussed Ulam-Hyers-Rassias stability and semi-Ulam-Hyers-Rassias stability of these solutions. Finally, we illustrated our results through case studies.

The polarised emission from AGN jets carries information about the physical conditions at the emitting plasma elements while its temporal evolution probes the physical processes that introduce variability and dynamically modify the local conditions. Here we present the analysis of multi-frequency radio linear and circular polarisation datasets with the aim to exactly quantify the conditions in blazar jets. Our analysis includes both the careful treatment of observational datasets and numerical modelling for the reproduction of synthetic polarisation curves that can be compared to the observed ones. In our approach the variability is attributed to traveling shocks. The emission from the cells of our jet model is computed with radiative transfer of all Stokes parameters. The model also accounts for Faraday effects which map the low-energy particle populations. We present two extreme cases in terms of the significance of Faraday conversion in the production of circular polarisation. As we show, in both regimes the model gives realistic reproduction of the observed emission.

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.

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

The study aims to explain the different effects between the multiplication, addition, division and subtraction equations. In this sense, multiplication and addition equations develop a new element or fact, meanwhile, division equations express information within the element or fact, subtraction is the selection of element or fact. Furthermore, for a multiplication and addition equation results in the new element or fact, it is necessary to depend on specific physical concepts, without extreme alteration, as it would result in the alteration of the element or fact. It is worth highlighting the importance of fixed values ​​in the equation to describe the specific nature of the element or fact. The study makes it possible to understand open equations (without solutions) in mathematics and physics, in addition to understanding existing equations. The study makes reference to the theory of obligatory necessity and the theory of differences between elements.

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

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

In this work, we study a boundary control problem for the evolutionary Stokes equations, under mixed boundary conditions, in two and three dimensions. The Stokes equation is a valuable tool in cardiovascular biomechanics, especially for modeling blood flow in low-velocity and high-viscosity conditions. Its application in the microvasculature, the design of medical devices, and the study of cardiovascular diseases helps improve the understanding and treatment of various pathological conditions. We provide a comprehensive theoretical framework to address the analysis and the derivation of a system of first-order optimality conditions that characterizes the solution of the control problem. Finally, solution-finding algorithm is proposed and illustrated with some simulations.

This paperfocusesontheoscillatorypropertiesofthesemi-canonicalnonlineardelaydifferential equation ofthird-order (δ2(t)(δ1(t)ζ′(t))′)′ + η(t)ζα(τ(t))= 0, t ≥ t0, under theconditions ∞t01δ2(t)dt = ∞ andZ ∞t01δ1(t)dt < ∞. Using the new canonical transform method, we transformed the studied equation into a canonical type equation and then we find the conditions for the oscillation of the studied equation from the canonical equations. The obtained results are new and complement to the existing results mentioned in the literature. Examples are provided to illustrate the importance and novelty of the main results.

The fundamental period is one of the most important parameters for the design of new structures , as well as the estimation of the capacity of existing ones. However, its calculation requires the solution of an eigenvalue problem. Especially for new buildings, this can involve material and geometric parameters that are not known beforehand. Thus, to estimate it, various design codes and researchers have adopted a several approximate analytical equations based on a number of key structural parameters. To this end, the present study introduces a novel methodology for the derivation of analytical equations for the fundamental period of Reinforced Concrete (RC) structures. The methodology is based on Machine Learning (ML) explainability techniques and specifically the so-called SHapley Additive exPlanations (SHAP) values. The novel methodology allows these equations to be constructed sequentially and in an informed manner, controlling the balance between accuracy and complexity. An extended dataset comprised of 4026 data points is employed. The results showed excellent accuracy, with a coefficient of determination R2≈0.95. This demonstrates the potential applicability of the proposed methodology in a wide array of similar engineering challenges.

The well-known equations for the process of powder compaction (PCP) in a rigid die published over the period from the beginning of the last century to the present are considered. Most of the considered equations are converted into the dependence of the densification pressure on the relative density. The equations were analyzed and their ability to describe PCP was assessed by determining the coefficient of determination when approximating experimental data on the compaction of various powders. It is shown that most of the equations contain two constants, the values of which are determined by fitting the mathematical dependence to the experimental curve. Such equations are able to describe PCP with high accuracy in the case of compaction of powders up to a relative density of 0.9-0.95. It is also shown that different equations can describe PCP in the range of density change from the initial to 0.9 with the same and sufficiently high accuracy, but when the process of compaction is extrapolated to higher values of density, the curves diverge. This indicates the importance of equations that can unambiguously describe PCP to a relative density equal to or close to 1.0. For an adequate description of PCP to a relative density greater than 0.95, equations containing three or four constants have proven themselves well.

Maxwell’s equations are valid only for a stationary observation point, therefore, to adequately describe real processes so far we have had to move to a moving reference frame. This paper presents the equations of electrodynamics for the moving observation point, it is shown that plane and spherical electromagnetic waves are their solutions, while the spherical wave propagates only outward, which cannot be said about Maxwell’s equations. The fields of uniformly moving charges are also solutions of the equations. Now there is no need to move to a moving reference frame, to use four-dimensional space and covariant form of equations. The question of finding a universal form of the equations that allows a solution in the form of the field of an arbitrarily moving charge remains open. This raises the question of the existence of a two-parameter group of transformations of electromagnetic fields along with the known one-parameter group has been posed. The phenomena derived from the equations, which make an additional contribution to the phase overrun in the Aharonov-Bohm effect are considered. The equation of motion of a charged particle in an electromagnetic field without simplifying approximations is considered, which allows us to take into account the radiation effects. It is shown that the fields in a moving observation point depend on its velocity and acceleration. In particular, although in a constant uniform electric field a force qE acts on a motionless charged particle, but on the same motionless but not fixed particle the force 4/3qE acts already, because it has a nonzero acceleration and the electric field at this point is larger. As the speed increases, the field decreases, and when it reaches the speed of light, when the particle stops accelerating, the force again becomes equal to qE The principle of operation of an unconventional alternator in a constant electric field and its corresponding engine, as well as new types of direct and impulse current generators, predicted by the equations, are described.

Boundary behavior of a given important function or its limit values are essential in the whole spectrum of mathematics and science. We consider some tractable cases of limit values in which either a difference of two ingredients or a difference equation is used coupled with the relevant functional equations to give rise to unexpected results. This involves the expression for the Laurent coefficients including the residue, the Kronecker limit formulas and higher order coefficients as well as the difference formed to cancel the inaccessible part, typically the Clausen functions. We also state Abelian results which yield asymptotic formulas for weighted summatory function from that for the original summatory function

In the paper, the authors view some ordinary differential equations and their solutions from the angle of (the generalized) derivative polynomials and simplify some known identities for the Bernoulli numbers and polynomials, the Frobenius-Euler polynomials, the Euler numbers and polynomials, in terms of the Stirling numbers of the first and second kinds.

In this paper, we compile the fractional power series method and the Laplace transform to design a new algorithm for solving the fractional Volterra integro-differential equation. For that, we assume the Laplace power series (LPS) solution in terms of power q=1m,m∈Z+, where the fractional derivative of order α=qγ for which γ∈Z+. This assumption will help us to write the integral, the kernel, and the nonhomogeneous terms as a LPS with the same power. The recurrence relations for finding the series coefficients can be constructed using this form. To demonstrate the algorithm's accuracy, the residual error is defined and calculated for several values of the fractional derivative. Two strongly nonlinear examples are discussed to provide the efficiency of the algorithm. The algorithm gains powerful results for this kind of problem. Under Caputo meaning the obtained results are illustrated numerically, and graphically. Geometrically, the behavior of the solution declares that the changing of the fractional derivative parameter values' in their domain alters the style of the attained solution in a symmetrical meaning and fully coinciding to the ordinary derivative value'. From these simulations, the results report that the recommended novel algorithm is a straightforward, accurate, and superb tool to generate analytic-approximate solutions for Integral, and integro-differental equations of fractional order.

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

Trail detection in mixed canopy ecosystems has important implications for forest management, monitoring, and conservation, although active sensor technology for sub-canopy trail detection is still developing. In order to assess the effectiveness of UAV(Unmanned Aerial Vehicle)-borne lidar (light detection and ranging) data for small trails (&lt; 2.5m width) in mixed forest canopy cover, we collected lidar data and trail characteristics (canopy cover and trail width) and created a high definition surface model map from the resulting lidar data, and also a high-resolution satellite imagery map using Google Earth. Through participatory mapping methods, seven respondents with limited prior geospatial experience completed a rapid identification of trails on both maps. Respondents’ trails were georeferenced in order to compare the rate of detectability between maps. We found greater detection on the lidar-derived map compared to the Google Earth map. Detectability in Google Earth maps was positively correlated with wider trails and trials with lower canopy. In lidar maps, trail detectability increased with wider trails, but canopy cover had no effect on detection rates. Our data indicate that a mixed-method approach that combines UAV-mounted lidar with high-resolution satellite imagery and participatory mapping increases rapid detection rates of small trails under varying canopy cover and trail widths.

The bark of Rhus verniciflua Stokes (RVS) is used as a food additive and herbal medicine for various inflammatory disorders and cancer in Eastern Asia. RVS exerted anti-inflammatory effects in lipopolysaccharide (LPS)-stimulated RAW264.7 cells, but whether this effect occurs in macrophages after oral administration has not been determined. We used a thioglycollate-induced peritonitis model to obtain macrophages from mice given RVS. We examined the systemic inflammatory response to intraperitoneal LPS. RVS-treated mice had an increased population of peritoneal exudate cells expressing CD11b and SRA. Increased uptake of Alexa Fluor 488-labeled acetylated lipoprotein was observed in monocyte-derived macrophages from RVS-treated mice. When these cells from the RVS group were stimulated with LPS, the levels of tumor necrosis factor (TNF)- and interleukin (IL)-6 in the supernatant decreased, but the level of IL-12 increased. The surface expression of CD86 was reduced, but surface expression of class II MHC molecules was increased. RVS suppressed the serum levels of LPS-induced TNF- and IL-6. RVS enhanced monocyte differentiation in thioglycollate-induced peritonitis by increasing scavenger receptor expression and activity. Macrophages isolated from mice given RVS responded differently to LPS. These findings demonstrate that RVS confers selective anti-inflammatory activity without causing the overall inhibitory effects on immune cells.

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.

We investigated the dynamics of a simple master equation describing the evolution of a probability distribution between two states, \( p_1 \) and \( p_2 \) (with \( p_1 + p_2 = 1 \)). Our focus is on the behavior of the entropy \( S \), the distance \( D \) from the uniform distribution (\( p_1 = p_2 = 1/2 \)), and the free energy \( F \). We define two ratios, \( \eta_S = \frac{dS/dt}{dF/dt} \) and \( \eta_D = \frac{dD/dt}{dF/dt} \), and find that both diverge to plus-minus infinity as the system evolves towards the uniform distribution. This divergence marks a critical transition point where the master equation temporarily loses physical meaning. However, after this divergence, the system settles into a new, well-behaved regime where the solutions stabilize at finite values, representing a new equilibrium or steady state. This two-regime behavior underscores the intricate dynamics of simple probabilistic systems as they approach and surpass equilibrium, providing insights into the complex interplay of entropy, distance, and free energy in statistical mechanics.

The research paper provides a mathematical framework for understanding phase shift in wave phenomena, bridging theoretical foundations with real-world applications. It emphasizes the importance of phase shift in physics and engineering, particularly in fields like telecommunications and acoustics. Key equations are introduced to explain phase angle, time delay, frequency, and wavelength relationships. The study also introduces the concept of time distortion due to a 1° phase shift, crucial for precise time measurements in precision instruments. The research also addresses infinitesimal wave energy loss related to phase shift, enriching our understanding of wave behavior, and impacting scientific and engineering disciplines.

Electrically charged particles such as Electrons and Protons carry electric, E, and magnetic, B, fields. In addition to these fields, Quantum Mechanics (QM) endows these particles with an `arcane and spooky' field --- the wavefunction. This wavefunction of QM is not only assumed to be separate but distinct from the electromagnetic field. We herein upend this view by demonstrating otherwise. That is, we demonstrate that the four components of the Dirac wavefunction, can be shown to not only be an intimate, but, a direct function of the electromagnetic field carried by the particle in question. Insofar as unity, depth in our understanding and insight into both Dirac and Maxwell's equations as major pillars of Modern Physics, we believe that this work may very well inch us one-step-closer to the truth.

Fluid mechanics of flow in hydrophobic, rectangular microchannels with finite aspect ratios is of paramount importance. In such microchannels, not only the effect of the side walls should be taken into account, but also the classical assumption of no-slip boundary condition (BC) is no longer valid at the solid-liquid interface. Accordingly, slip flow can occur in microchannels fabricated from surfaces with low wetting conditions, hydrophobic surfaces. Determining the interactions of liquid molecules adjacent to solid surface is still a challenging issue, and it is especially important in small scale domains. Herein, the fluid mechanics of flow through rectangular hydrophobic microchannels has been reconsidered by taking into account the general Navier-slip BCs at the solid-liquid interface. For fully developed incompressible flow in microchannels at low Reynolds number, partial differential equation (PDE) of the momentum equation simplifies to the classical Poisson equation. Accordingly, by analytically solving the Poisson equations with general Navier-slip BCs, the most general forms of velocity distributions, flow rate, friction factor and Poiseuille number have been obtained.

The objective of this study is to evaluate the effect of constitutive equations on the prediction accuracy for springback in cold stamping with various deformation modes. In this study, two types of yield functions—Hill’48 and Yld2000-2d—were considered to describe yield behavior. Isotropic and kinematic hardening models based on the Yoshida–Uemori model were also adopted to describe hardening behavior. Various material tests (such as uniaxial tension, tension- compression, loading-unloading, and hydraulic bulging tests) were carried out to determine the material parameters of the models. The obtained parameters were implemented in the finite element (FE) simulation to predict springback, and the results were compared with experimental data. U-bending and T-shape drawing were employed to evaluate the springback prediction accuracy. Obviously, the springback prediction accuracy was greatly influenced by constitutive equations. Therefore, it is important to choose appropriate constitutive equations for accurate description of material behaviors in FE simulation.

In this paper, we study the problem of controllability of impulsive neutral evolution integrodifferential equations with state dependent delay in Banach spaces. The main results are completely new and are obtained by using Sadovskii's fixed point theorem, theory of resolvent operators, and an abstract phase space. An example is given to illustrate the theory.

The second article in the series presents the application of the Smoothed Particle Hydrodynamics (SPH) method to modelling of batch crystallisation in stirred tanks. A methodology to integrate the population balance equations (PBE) in parallel and independently from the Navier-Stokes equations is demonstrated. The benefits of the proposed methodology in terms of computational requirements, accuracy and availability of the crystal size distribution are discussed. The specific formulation of the SPH equations where the resulting system of ordinary differential equations is solved using the weighted contributions rather than numerically by solving a linear system of equations allows for massive parallelisation and a very loose coupling of the population balance and the fluid dynamics. It has been demonstrated, that the population balance equations can be solved on a Shared Memory Architecture (SMA) system using the OpenMP interface while the fluid dynamics equations being computed independently on a General Purpose Graphics Processing Unit (GPGPU) using the NVidia CUDA technology. This way, a significant portion of the computational overhead due to the large number of additional transport equations resulting from the discretisation of the population balance was removed: the SPH simulation coupled with 200 population balance equations was only 40% slower compared to SPH-only simulation. Two methods for the solution of population balance equations that preserve full crystal size distribution were implemented: discretised population balance (DPB) and method of characteristics (MOCH). The DPB equations are solved using the high-resolution finite-volume method with flux limiter and the effect of a large number of different flux limiters have been investigated. Both methods were validated using the case studies from the literature where an analytical solution can be derived. The developed models were applied to a numerical solution of coupled computational fluid dynamics and population balance equations to model a batch crystallization process. The effect of the hydrodynamics on the local temperature/supersaturation and the resulting crystal size distribution was captured and compared to the ideal mixing case. The simulation results from the DPB and MOCH methods were compared in terms of computational requirements and accuracy and MOCH selected as computationally more efficient and accurate.

We explored an extension of the Cauchy-Riemann equations to the algebra of treons, recently described by Alejandro Bermejo, whose elements are ordered 3-tuples. We leveraged the isomorphism between the algebra of treons and algebra B, and deduced the Cauchy-Riemann equations for the algebra of treons, establishing the necessary conditions for analyticity in this algebraic structure. This work significantly broadened our horizons in complex analysis and introduced new possibilities for applications across various fields of advanced mathematics.

The paper considers the physical consequences of bimetric geodesically (pro jectively) invariant equations of gravitation. The correct expression for the energy- impulse tensor of the gravitational field is obtained and with its help the key role of the gravitational defect of mass-energy for the energy of the newborn universe is shown..

Most fully developed galaxies have vivid spiral structure, but the formation and evolution of spiral structure is still a mystery that is not fully understood in astrophysics. We find that the currently used equations of galactic dynamics contain some unreasonable components. In this paper, the following three working assumptions are introduced to simplify the galactic structural equations. 1. In the research of large-scale structure, the retarded potential of the gravitational field should be taken into account. The propagating time of the gravitational field from center to border is longer than the revolution periods of the stars near the center of galaxy. Newton's gravitational potential is unreasonable for such case, and the weak field and low velocity approximation of Einstein's field equation should be adopted. 2. The stars in a fully developed galaxy should be zero-pressure and inviscid fluid, and the equation of motion is different from that of ordinary continuum mechanics. Stars move along geodesics. 3. The structure of the galaxy is only related to the total mass density distribution. The equation of state of dark halo is different from that of ordinary luminous interstellar matter, so their trajectories are also very different. Dark halo and ordinary matter in galaxy are automatically separated. The total mass density distribution can be presupposed according to the observation data, and then it can be determined by comparing the solution of the equations with the observed data. These assumptions and treatments are supported by theory and observation. The variables of the equations of simplified galactic dynamics are separated from each other, and the equations are well-posed and can be solved according to a definite procedure. Therefore, this simplified dynamic equation system provides a more reasonable and practical framework for the further study of galactic structure, and can solve many practical problems. Besides, it is closely related to the study of dark energy and dark matter.

Traditional population-based metaheuristic algorithms are effective in solving complex real-world problems but require careful strategy selection and parameter tuning. Parameter-less population-based optimization algorithms have gained importance due to their simplicity and efficiency. However, research on their applicability for solving large systems of nonlinear equations is still incipient. In this paper, a systematic review and detailed description of metaphor-less optimization algorithms is made, including the Jaya algorithm, the main Jaya variants for continuous optimization problems, with a focus on the Enhanced Jaya (EJAYA) algorithm, the three Rao algorithms, the best-worst-play (BWP) algorithm, and the very recent max-min-greedy-interaction (MaGI) algorithm. This paper also discusses a few previous parallelizations of the Jaya optimization algorithm, as well as recent GPU-based massively parallel implementations of the Jaya, Enhanced Jaya (EJAYA), Rao, and BWP algorithms developed by the authors. In addition, a novel GPU-accelerated version of the MaGI algorithm is proposed. The GPU-accelerated versions of the metaphor-less algorithms developed by the authors were implemented using the Julia programming language and tested primarily on GeForce RTX 3090 and NVIDIA A100 GPUs, as well as on Tesla V100S and Tesla T4 GPUs, using a set of difficult, large-scale nonlinear equation system problems. The computational experiments carried out produced quite significant speedups, which highlights the efficiency of the GPU-based versions of the metaphor-less algorithms developed for solving large-scale systems of nonlinear equations.

In this paper we introduce the Adomian decomposition method for a class of first order fuzzy dynamic equations on arbitrary time scales for existence of solutions. The results are provided with suitable numerical examples.

We consider a Cauchy problem for differential equations in a Hilbert space $X$. The problem is stated in a time interval $I$, which can be finite or infinite. Weuse a fixed point argument for history-dependent operators to prove the unique solvability of the problem. Then, we state and prove convergence criteria for both a general fixed point problem and thecorresponding Cauchy problem. These criteria provide necessary and sufficient conditions on a sequence $\{u_n\}$ which guarantee its convergence to the solution of the corresponding problem, in the space of both continuous and continuously differentiable functions. We then specify our results in the study of a particular differential equation governed by two nonlinear operators. Finally, we provide an application in viscoelasticity and give a mechanical interpretation of the corresponding convergence result.

In the paper, the authors present an explicit form for a family of inhomogeneous linear ordinary differential equations, find a more significant expression for all derivatives of a function related to the solution to the family of inhomogeneous linear ordinary differential equations in terms of the Lerch transcendent, establish an explicit formula for computing all derivatives of the solution to the family of inhomogeneous linear ordinary differential equations, acquire the absolute monotonicity, complete monotonicity, the Bernstein function property of several functions related to the solution to the family of inhomogeneous linear ordinary differential equations, discover a diagonal recurrence relation of the Stirling numbers of the first kind, and derive an inequality for bounding the logarithmic function.

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

We explore the field equations in a 4-d complex space-time, in the same way, that general relativity does for our usual 4-d real space-time, forming this way, a new "general relativity" in C⁴ space-time, free of sources. Afterwards, by embedding our usual 4-d real space-time in C⁴ space-time, we describe geometrically the energy-momentum tensor T_μν as the lost geometric information of this embedding. We further give possible explanation of dark eld and dark energy.

When discussing the Diophantine quintic equation $p(a^5 + b^5) = q(c^5 + d^5)$ where p is a prime and q is an integer, there is a clear gap between the mathematical literature and online forums. Because of the intrinsic complexity of parameterizing fifth-degree equations, this equation re- mains largely unexplored. In this work, we approach this quintic problem through algebraic methods with the goal of providing numerical solutions that illuminate its properties and behaviors.Our research uncovers in- triguing patterns and connections that shed light on the enigmatic nature of Diophantine quintic equations in this particular form.

The purpose of this paper is to compare daily sunshine duration (SD) measured by Campbell-Stokes sunshine recorder (CS) with three automatic instruments (Kipp & Zonen pyranometers (CM6B), pyrheliometer (CHP1) and CSD3 sunshine recorders). The comparisons are essential, since in recent decades automatic sunshine duration sensors have been systematically replacing the manual measurements affecting therefore the continuity of the CS time series. The comparisons were applied to the daily SD data of three synoptic stations (Larnaca A/P, Pafos A/P and Athalassa) where parallel measurements of manual and automatic instruments are carried out during the periods 2009-2012 and 2012-2015. The CSD3 measurements were obtained from two automatic weather stations (AWS) (Mennogia and Pentakomo) which are closed to the synoptic station of Larnaca A/P. The comparisons cover the daily and monthly sums, their statistical distributions and their relationships. The daily differences of the pairs of the sunshine recorders installed at the same location are mostly within the range of -1.0 to +1.0 h. At Larnaca, in the period 2009-2012, the pyranometric method (Gl) is higher than the CS in the summer months, while during the winter months it is lower than CS. Pafos and Athalassa showed higher values of CS comparing to Gl and the pyrheliometric method (Pr), respectively, almost throughout the year. Generally, the annual average daily differences between the different pairs range between 0.3 and 0.5 h. Regression equations were determined allowing the estimation of daily SD under different sky conditions. The closest relations between different pairs were obtained in ‘clear sky’ and ‘all sky’ conditions, as it is indicated by the high coefficient of determinations. R2 reached the values of 0.92 and 0.96 under ‘all sky’ conditions, and between 0.80 and 0.91 under ‘clear’ sky conditions, respectively. An increase in cloud cover causes an increase in the SD differences recorded by the relevant devices. Comparing statistically the different pairs using the Mann-Whitney nonparametric test which is the alternative to the Two Sample T parametric test, it was found that the medians are not statistically different. However, significant differences were only detected in the case of the pair of CS and Pr at Athalassa (2012-2015) using both tests. In the rest of the cases the differences are not statistically significant at the 5% level.

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.

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.

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.

We have shown that the use of evolutionary gravitational constants G, the speed of light c, and cosmological constant Λ in the Einstein’s field equations leads to a very simple model that fits the supernovae Ia data with a single parameter as well as the standard ΛCDM model with two parameters, and has the predictive capability superior to the latter. We have shown that at the current time G and c both increase as dG/dt = 5.4GH₀ and dc/dt = 1.8cH₀ with H₀ as the Hubble constant, and Λ decreases as dΛ/dt = −1.2ΛH₀. Negative finding on the variation of G with time derived from the analysis of lunar laser ranging data in several studies has been attributed to the fact that these studies did not consider the variation of c simultaneously with G.

The autosomal recessive disorder Ataxia-Telangiectasia is caused by dysfunction of the stress response protein ATM. In the nucleus of proliferating cells, ATM senses DNA double-strand breaks and coordinates their repair. This role explains T-cell dysfunction and tumor risk. However, it remains unclear whether this function is relevant for postmitotic neurons and underlies the cerebellar atrophy, since ATM is cytoplasmic in postmitotic neurons. Here, we used ATM-null mice that survived early immune deficits by bone-marrow transplantation, and reached initial neurodegeneration stages at 12 months of age. Global cerebellar transcriptomics demonstrated ATM depletion to trigger upregulations in most neurotransmission and neuropeptide systems. Downregulated transcripts were found for the ATM interactome component Usp2, many non-coding RNAs, ataxia genes Itpr1, Grid2, immediate early genes and immunity factors. Allelic splice changes affected prominently neuropeptide machinery, e.g. Oprm1. Validation experiments with stressors were performed in human neuroblastoma cells, where ATM localized only to cytoplasm, similar to brain. Effect confirmation in SH-SY5Y cells occurred after ATM depletion and osmotic stress better than nutrient / oxidative stress, not after ATM kinase inhibition or DNA stressor bleomycin. Overall, we provide pioneer observations from a faithful A-T mouse model, which suggest general changes in synaptic and dense-core vesicle stress adaptation.

We prove general stability theorems for $n$-dimensional quartic-cubic-quadratic-additive type functional equations of the form \begin{eqnarray*} \sum_{i=1}^\ell c_i f \big( a_{i1}x_1 + a_{i2}x_2 + \cdots + a_{in}x_n \big) = 0 \end{eqnarray*} by applying the direct method. These stability theorems can save us much trouble of proving the stability of relevant solutions repeatedly appearing in the stability problems for various functional equations.

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

In this study, we focus on an abstract form of summation equations, which can be regarded as a particular class of nonlinear difference equations. We investigate sufficient conditions for the existence of discrete almost automorphic solutions by appealing to the fixed point theory and establish a linkage between the existence of discrete almost automorphic solutions and the existence of bounded solutions for the nonlinear difference equations in the spotlight. Consequentially, we propose a discrete counterpart of the Bohr-Neugebauer theorem for the qualitative theory of difference equations.

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.

We establish certain modular equations of degree three based on some results of Ramanujan's theta functions.

Considerable efforts have been devoted to modifying gravity, which aim to elucidate the possible existence or nature of dark matter and energy, describe observational data more effectively, and formulate quantum gravity. In addition, despite the immense success of the quantum field theory, the framework requires renormalization techniques and breaks down at high energies. Recently, the Planck legacy 2018 release has confirmed the existence of an enhanced lensing amplitude in the cosmic microwave background power spectra, which prefers a positively curved early Universe with a confidence level higher than 99%. This study considers the implied curvature of the early Universe as the curvature of “the background or 4D conformal bulk as a manifestation of vacuum energy” and distinguishes it from the localized curvature that is induced in the bulk by the presence of celestial objects that are regarded as ”4D relativistic cloud-worlds.” Analogously, because gravity appears to emerge owing to spacetime curvature and does not exhibit critical characteristics shared by other fields, it has been incorporated as the local curvature of the bulk affecting the embedded quantum fields that are regarded as propagating “4D relativistic quantum clouds.” To consider the effects of the bulk on embedded clouds, this paper presents interaction field equations in terms of brane-world modified gravity and the perspective of geometrization of the quantum mechanics wherein gravity is manifested by the curvature of 4D conformal bulk as an indicator of the field strength of vacuum energy on the embedded 4D relativistic clouds in addition to the boundary interactions, which could remove the singularities and satisfy a conformal invariance theory. A visualization of the evolution of the 4D relativistic cloud-worlds over the conformal spacetime of the 4D bulk is presented, whereas the standard said theories can be recovered from the interaction field equations.

Autonomous differential equations of fractional order and non-singular kernel are solved. While solutions can be obtained through numerical, graphical, or analytical solutions, we seek an implicit analytical solution.

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

In this paper we 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.

For the KdV-Burgers equations on cylindrical and spherical waves the development of a regular profile starting from an equilibrium under a periodic perturbation at the boundary is studied. The equations describe a medium which is both dissipative and dispersive. Symmetries, invariant solutions and conservation laws are investigated. For an appropriate combination of dispersion and dissipation the asymptotic profile looks like a periodical chain of shock fronts with a decreasing amplitude (sawtooth waves). The development of such a profile is preceded by a head shock of a constant height and equal velocity which depends on spatial dimension as well as on integral characteristics of boundary condition; an explicit asymptotic for this head shock and a median of the oscillating part is found.

The storm of November 12th-13th, 2019 provoked the displacements of boulders in a central Mediterranean rocky coast; with reference to a selected area, prone to the boulder production and geomorphologically monitored for years, a field-oriented study approach was applied for the phenomenon, by collating data concerning pre-storm locations and kinematics of these boulders. The number of displaced boulders is 11, that is, in terms of morphological imprint of a specific storm, one of the major study cases for the Mediterranean. In addition, based on widely used hydrodynamic equations, the minimum wave height required to displace the boulders is assessed. The values conform with the expected values for the wave climate dominating during the causative meteorological event and give a measure of the energy of the storm slamming the coast. Boulder dislodgements usually have plays a key role in determining the rate of the coastal recession, likely also in the investigated area. In view of an adverse climate evolution with a possible increase of energy and frequency of severe storms, the results deriving from the study of this morphodynamics should be considered for the hazard assessment and coastal management.

The prevalence of differential equations as a mathematical technique has refined the fields of control theory and constrained optimization due to the newfound ability to accurately model chaotic, unbalanced systems. However, in recent research, systems are increasingly more nonlinear and difficult to model using Differential Equations only. Thus, a newer technique is to use policy iteration and Reinforcement Learning, techniques that center around an action and reward sequence for a controller. Reinforcement Learning (RL) can be applied to control theory problems since a system can robustly apply RL in a dynamic environment such as the cartpole system (an inverted pendulum). This solution successfully avoids use of PID or other dynamics optimization systems, in favor of a more robust, reward-based control mechanism. This paper applies RL and Q-Learning to the classic cartpole problem, while also discussing the mathematical background and differential equations which are used to model the aforementioned system.

This paper presents machine learning methods for approximate solutions of reaction-diffusion equations with multivalued interaction functions. This approach addresses the challenge of finding all possible solutions for such equations, which often lack uniqueness. The proposed method utilizes physics-informed neural networks (PINNs) to approximate generalized solutions.

Recently several papers have been published on the Liénard equation x'' +L x+m x^3+n x^5=0, where the authors studied their explicit solutions and their applications. Here we describe the complete dynamics of these differential equations in the Poincar\'e disc. The Poincar\'e disc is the closed disc centered at the origin of coordinates of radius one, the whole plane R^2 is identified with the interior of this disc, and its boundary, the circle $\sss^1$ is identified with the infinity of R^2.

We explore the relationship between the exact renormalization group (RG) equation and quantum mechanics. The central idea is to rewrite the exact Wilsonian RG equation as a Schrödinger-type equation, opening avenues for analysis. By introducing a functional transformation, the RG equation takes on a form analogous to a quantum mechanical Schrödinger equation. The study investigates solutions to this equation and draws parallels with quantum mechanics, paving the way for understanding the infrared behavior of quantum field theories. Furthermore, we extends its insights to a one-dimensional toy model, providing illustrative examples. The discussion encompasses the implications of these findings in multiple dimensions and quantum field theory, with a brief exploration of their relevance in the context of inflationary cosmology.

In this paper, we consider with a class of system of differential equations whose argument transforms are involution. In this an initial value problem for a differential equation with involution is reduced to an initial value problem for a higher order ordinary differential equation. Than either two initial conditions are necessary for a solution, the equation is then reduced to a boundary value problem for a higher order ODE.

We present a class of efficient parametric closure models for 1D stochastic Burgers equations. Casting it as statistical learning of the flow map, we derive the parametric form by representing the unresolved high wavenumber Fourier modes as functionals of the resolved variables’ trajectory. The reduced models are nonlinear autoregression (NAR) time series models, with coefficients estimated from data by least squares. The NAR models can accurately reproduce the energy spectrum, the invariant densities, and the autocorrelations. Taking advantage of the simplicity of the NAR models, we investigate maximal and optimal space-time reduction. Reduction in space dimension is unlimited, and NAR models with two Fourier modes can perform well. The NAR model’s stability limits time reduction, with a maximal time step smaller than that of the K-mode Galerkin system. We report a potential criterion for optimal space-time reduction: the NAR models achieve minimal relative error in the energy spectrum at the time step where the K-mode Galerkin system’s mean CFL number agrees with the full model’s.

Based on previous research results, we propose a new preprocessing HSS iteration method (PHSS) for the generalized Lyapunov equation. At the same time, the corresponding inexact PHSS algorithm (IPHSS) is given from the angle of application. All the new methods presented in this paper have given the corresponding convergence proof. The numerical experiments are carried out to compare the new method with the existing methods, and the improvement effect is obvious. The feasibility and effectiveness of the proposed method are proved from two aspects of theory and calculation.

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

In the paper, by virtue of techniques in combinatorial analysis, the authors simplify two families of nonlinear ordinary differential equations in terms of the Stirling numbers of the first kind.

Maxwell's equations accurately describe the propagation of electromagnetic (EM) waves. However, the conventional formulation implies invariance to observers, posing challenges in explaining how the wave is perceived by different observers. For example, Doppler effect shows that observers perceive the same EM wave differently.By employing mathematical transformations, we derive a general form of Maxwell’s wave equations by incorporating a "time scaling" factor to account for observer perception. The original form is shown to be a special case when static. The Doppler effect and Sagnac effect are directly derived from Maxwell’s equations. All the results are consistent with established experiments.Our findings offer a fresh perspective, promising a deeper understanding and unification of EM phenomena. We extend Maxwell’s equations to describe not only the propagation of an EM wave, but also how it will be measured by different observers, say, antennas.

This manuscript presents a novel approach to estimating the potential within Schrödinger equations, with a particular focus on applications to the Navier-Stokes problem in fluid dynamics. By establishing new theoretical estimates, we delve into the intricate dynamics of fluid flow, aiming to unveil previously obscured aspects of the Navier-Stokes equations. Through a rigorous mathematical framework, we explore how these new potential estimates can provide a fresh perspective on fluid mechanics, contributing to the ongoing quest to solve some of its most persistent challenges.

The fractional derivatives are a generalization the derivatives of integer order and find applications in studying memory processes in various scientific fields. Numerical methods are used to solve and analyze fractional models of real world problems. In this paper, we consider an approximation of the Caputo fractional derivative and its asymptotic formula, whose generating function is the polylogarithmic function. In the paper, we prove the convergence of the approximation and derive an estimate for the error and order. We consider an application of the approximation for the construction of finite difference schemes for numerical solution of the two-term ordinary fractional differential equation and the time-fractional Black-Scholes equation for option pricing. The properties of the approximation are used to prove the convergence of the methods used for numerical solution of the fractional differential equations. The theoretical results for the error and order of the methods are illustrated by the results of numerical experiments.

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.

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

In this paper, we 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.

Good Boussinesq equations will be considered in this work. First we apply three combined compact schemes to approximate spatial derivatives of good Boussinesq equations. Then three fully discrete schemes are developed based on symplectic scheme in time direction, which are sympletic-structure preserving.Meanwhile, the convergence and conservation of the fully discrete schemes are analyzed. Finally, we present numerical experiments to confirm our theoretical analysis. Both our analysis and numerical test indicate that the fully discrete schemes are efficient in solving the spatial derivative mixed equation.

We propose a new hierarchy of the vector derivative nonlinear Schr{\"o}dinger equations and consider the simplest multiphase solutions of this hierarchy. The study of the simplest solutions of these equations led to the following results. First, the three-leaf spectral curves $\Gamma=\{(\mu,\lambda)\}$ of the simplest multiphase solutions have a quite simple symmetry. They are invariant with respect to holomorphic involution $\tau$. The type of this involution depends on the genus of the spectral curve. Or the involution has the form $\tau:(\mu,\lambda)\to(\mu,-\lambda)$, or $\tau:(\mu,\lambda)\to(-\mu,-\lambda)$. The presence of symmetry leads to the fact that the dynamics of the solution is determined not by the entire spectral curve $\Gamma$, but by its factor $\Gamma/\tau$, which has a smaller genus. Secondly, it turned out that the dynamics of the two-component vector $\bp=(p_1,p_2)^t$ is determined, first of all, by the dynamics of its length $\abs{\bp}$. Independent equations determine the dependence of the direction of the vector $\bp$ from its length. In cases where the direction of the vector $\bp$ is fixed, the corresponding spectral curve splits into separate components. In conclusion, we note that, as in the case of the Manakov system, the equation of the spectral curve is invariant with respect to the orthogonal transformation of the vector solutions. I.e., the solution can be found from the spectral curve up to the orthogonal transformation. This fact indicates that the spectral curve does not depend on the individual components of the solution, but on their symmetric functions. Thus, the spectral data of multiphase solutions have two symmetries. These symmetries make it difficult to reconstruct signals from their spectral data. The work contains examples illustrating these statements.