System Of Differential Equations Research Articles

Analytical simulation of Darcy–Forchheimer nanofluid flow over a curved expanding permeable surface

Abstract This research paper presents an analytical simulation of Darcy–Forchheimer flow over a porous curve stretching surface. In fluid dynamics, the Darcy–Forchheimer model combines Forchheimer adjustment and high-velocity effects with Darcy’s formula for porous media flow: two nanofluid particles, molybdenum disulphide, and graphene oxide, form nanofluid with the base fluid blood. The governing partial differential equations for momentum and energy are converted into a nonlinear ordinary differential equations system by applying the appropriate similarity transformations. The homotopy analysis method is used to solve the transform equations analytically. The impact of essential factors includes the Forchheimer parameter, porosity parameter, slip parameter, Eckert number, nanoparticle volume friction, magnetic field parameter, and curvature parameter. The results have applications in the design of sophisticated cooling systems, where effective thermal control is essential.

SDIR Mathematical Approach with Time Delay Factor to Reveal Indonesian Students' Interest in Government Organizations and Internships: A Case in FMIPA UNM

This study is to build a time-delay SDIR model on the case of student interest in Organizations and Government Internship Programs, analyze the model and conduct a model simulation to predict the level of student interest cases in Organizations. This study is a theoretical and applied study. Analysis of the time-delay SDIR model of the case of student interest in Organizations, while the model simulation uses Maple Software. The population of the study was active students of FMIPA UNM, with a sample size of 1029 students obtained using the Sloving technique. The results of this study are a mathematical model of SDIR on the case of student interest in Organizations which is a system of differential equations. The model analysis provides a value of the free equilibrium point of the case of student interest in Organizations and a stable endemic equilibrium point. The results also found that the basic reproduction number value for cases without a solution would produce R0n = 9.912507841, which means that in social cases where one individual can influence 9-10 people in their environment to hesitate to participate in organizational activities or internship programs, but on the other hand if the case is given a solution, it will produce R0s = 0.2737372211, which means that there is no psychological spread, where each individual does not influence other individuals.

A mathematical model to study the role of dystrophin protein in tumor micro-environment.

In this research work, the authors have developed a mathematical model to examine the interaction between dystrophin protein and tumor. The authors formulated a system of ordinary differential equations to describe the dynamics of the dystrophin-tumor interaction system. Jacobian matrix and Routh-Hurwitz stability techniques were used to determine equilibrium points, perform stability and bifurcation analysis, and establish the conditions required for the stability of the proposed model. Numerical simulations are performed using Euler's method to investigate the temporal evolution of the proposed model under different parameter values, such as tumor growth rate and feedback strength of dystrophin protein. The numerical results are presented in tables, and corresponding to each table, a graphical analysis is done. The graphical analysis includes creating phase portraits to visually represent stability regions around the equilibrium points, bifurcation diagrams to identify critical points, and time series analysis to highlight the behavior of the proposed model. The authors explore how variations in dystrophin expression impact tumor progression, identifying potential therapeutic implications of maintaining higher dystrophin levels. This comprehensive analysis enhances our understanding of the dystrophin-tumor interaction, providing a basis for further experimental validation and potential therapeutic strategies.

A study for analysing predator–prey ecological system using Legendre polynomial

This paper presents a mathematical model to examine the effects of the coexistence of predators on single prey. Based on Caputo operators, we present a newly developed system of differential equations for the predator–prey system using wavelet method. It is well known that a system of nonlinear singular models cannot operate smoothly since they are singular and nonlinear. Therefore, with the help of this numerical approach, we have converted the system into a nonlinear system of algebraic equations by extending it through operational matrix of Legendre wavelets. Using the wavelet collocation scheme, we have calculated these unknown coefficients. It has been demonstrated in tables and graphs that the developed approach is consistent and proficient. Further bifurcation diagrams, as well as phase portraits, have been used to study the proposed system numerically and to analyse its behaviour. In addition, a nonlinear functional analysis have used to establish uniformly boundedness for the proposed model. Also we have discussed residual error analysis and Lyapunov exponent. The applicability and efficacy of this methodology have been demonstrated through this nonlinear system. Additionally, a comparison with existing results highlights the advantages of our numerical approach. All calculations have been done using MATLAB.

Lie symmetry method to discuss the optimal system of Lie subalgebras and similarity solutions for shock propagation in a perfectly conducting dusty gas with magnetic field in rotating medium

In this paper, the similarity solutions for the shock wave propagation in a perfectly conducting dusty gas under the effect of axial or azimuthal magnetic field in the case of rotating medium by using Lie symmetry method have been discussed. The optimal system of one-dimensional Lie subalgebras related to the hyperbolic system of nonlinear partial differential equations had been constructed by using general invariant function and general adjoint transformation matrix. The optimal system is the collection of inequivalent optimal classes, which provides crucial information regarding different types of the similarity solutions. On the basis of the optimal system, we have discussed all possible similarity solutions for inequivalent optimal classes individually. Also, we have found that the similarity solutions exist in two cases for power law shock path and in one case with exponential law shock path. The shock wave decay is due to an increment in the physical parameter such as shock Cowling number or non-idealness parameter or adiabatic exponent or rotational parameter. For smaller value of the ratio of density of solid particles to initial species density of the gas i.e. [Formula: see text], the shock strength decreases with an increment in the mass fraction of solid particles in the mixture [Formula: see text], but the shock strength decreases as [Formula: see text] increases for [Formula: see text]. Also, the shock strength increases with an increment in [Formula: see text] at constant [Formula: see text]. It is observed that the shock wave is stronger in the presence of axial magnetic field as compared to azimuthal magnetic field.

Vector-borne disease outbreak control via instant releases.

This paper is devoted to the study of optimal release strategies to control vector-borne diseases, such as dengue, Zika, chikungunya and malaria. Two techniques are considered: the sterile insect one (SIT), which consists in releasing sterilized males among wild vectors in order to perturb their reproduction, and the Wolbachia one (presently used mainly for mosquitoes), which consists in releasing vectors, that are infected with a bacterium limiting their vectorial capacity, in order to replace the wild population by one with reduced vectorial capacity. In each case, the time dynamics of the vector population is modeled by a system of ordinary differential equations in which the releases are represented by linear combinations of Dirac measures with positive coefficients determining their intensity. We introduce optimal control problems that we solve numerically using ad-hoc algorithms, based on writing first-order optimality conditions characterizing the best combination of Dirac measures. We then discuss the results obtained, focusing in particular on the complexity and efficiency of optimal controls and comparing the strategies obtained. Mathematical modeling can help testing a great number of scenarios that are potentially interesting in future interventions (even those that are orthogonal to the present strategies) but that would be hard, costly or even impossible to test in the field in present conditions.

Analytical solutions and instability analysis of truncated M-fractional coupled dispersionless equations

Abstract This paper comprehensively investigates the truncated M-fractional coupled dispersionless equations, a nonlinear system of partial differential equations characterized by its M-fractional derivative. The Jacobi elliptic function expansion method is employed to derive analytical solutions for the coupled system. In addition, the modulation instability of the solutions is thoroughly explored, providing a detailed exposition of the mathematical framework governing the system. The analytical solutions are graphically illustrated and analyzed to highlight their physical significance. These findings have significant applications in nonlinear optics, offering new insights into wave propagation and stability within such systems.

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

Non-axisymmetric coupled non-stationary problem of thermoelectroelasticity for a long piezoceramic cylinder

A new closed solution to the non-axisymmetric coupled non-stationary problem of thermoelectroelasticity was constructed for a long piezoceramic cylinder for the case of satisfaction of the first and the third kind boundary conditions. Cylindrical surfaces were made as electrodes and connected to a measurement device with large input resistance. Limitation of a temperature change “load” rate made it possible to include equations of statics, electrostatics and thermal conductivity in the initial formula. The finite biorthogonal transforms are applying to explore a non-selfadjoint system of differential equations and to develop a closed solution. The obtained relations made it possible to determine the temperature and electric fields, and the stress-strain state in the piezoceramic cylinder, as well as the potential difference between cylindrical surfaces (electrodes) under non-stationary non-axisymmetric temperature impact.

New aspects of the development of Kantorovich’s one-product dynamic model of replacement of production funds

In the context of the development of scientific and technological progress and the growth of the capital stock, an important task is modeling and optimization of the terms of operation of production funds. The development of technologies and automation of production lead to a reduction in the workforce employed in production. In the article, for the previously developed single-product dynamic model of the replacement of production assets, taking into account the inertial properties of the funds being introduced, the case of a reduction in labor resources under the condition of an increase in output and capital investment is investigated. A solution was obtained that allows to determine the optimal strategy for the withdrawal of obsolete funds and the introduction of the new ones in case of decrease in labor resources. As the optimization criterion the principle of differential optimization is used. The theorem of equivalence of a trajectory satisfying the principle of differential optimization, the obtained solution is given. A system of functional differential equations is presented, which should be satisfied by variable models both in terms of growth and reduction of labor resources. Numerical methods are used to solve such a system. The variants of the system development are considered under various assumptions regarding changes in the total amount of labor resources. The results of numerical implementation of the model are presented.

PARASIAN OVER PARISIAN, HOW MUCH EARLIER SHOULD ONE EXERCISE?

In this paper, an integral equation approach used for pricing American-style Parisian options is extended to pricing Parasian options, after overcoming an additional difficulty of losing the “reset” feature of the latter in addition to still dealing with the high nonlinearity associated with American-style exercises. More specifically, compared with an American Parisian option, the accumulative feature of the “tracking clock” time results in a pair of coupled three-dimensional (3D) Partial Differential Equation (PDE) systems instead of coupled two-dimensional (2D) and 3D PDE systems. Upon successfully establishing the couple integral equation systems, we have successfully computed the option price and the optimal exercise price and compared the results with those of its Parisian counterpart. Through the comparison, we can directly observe how the accumulativeness of the “tracking clock” time affects the option and the optimal exercise prices and provide some meaningful financial explanations.

Evaluating the risk of malaria transmission within the Central African Republic with the goal of stabilising and eliminating the infection

In this paper, we formulate a multi-compartmental mathematical model for humans and mosquitoes. We construct the system of differential equations for an SEITVR for the human compartment and an SEI for the mosquito compartment. We investigate the outbreak of malaria and its effect on the Central African Republic. The analysis of the compartmental model is carried out using stability analysis and Routh Hurwitz Criterion technique is used to indicate the major impact of the model and to improve the model through minor modifications in the transformation of disease in the population. Our model exhibits two equilibrium points, disease free equilibrium points and endemic equilibrium points. The next generation matrix is used to determine the basic reproduction number R 0. A new compartmental model was framed and estimated the malaria spread after 2023 in the Central African Republic, which is the novelty of this research. Our main motivation is to make the Central African Republic a malaria free country. A Numerical example are provided to validate our results for both the disease free state and endemic state of each model. We believe that this investigation will be more effective in reducing malaria infection and stop spreading.

Numerical Simulation Study of Thermal Performance in Hot Water Storage Tanks with External and Internal Heat Exchangers

This paper presents a numerical analysis of two hot water storage tank configurations—one equipped with an external heat exchanger (Tank-1) and the other with an internal heat exchanger (Tank-2). The objective is to evaluate and compare their thermal performance during charging and discharging processes. The numerical model is developed by solving a system of ordinary differential equations using the 4th-order Runge–Kutta method, implemented in the Python programming language. The results indicate that Tank-1 demonstrated a higher charging efficiency of 94.6%, achieving full charge in approximately 2 h and 20 min. In comparison, Tank-2 required 3 h and 47 min to reach full charge, with a charging efficiency of 85.9%. During discharge, both configurations exhibited similar behavior, with an efficiency of 13.63% over approximately 33 min. The analysis showed that the external heat exchanger configuration led to more effective thermal stratification, supported by the Richardson number analysis, which indicated a significant effect of buoyancy during charging. This design advantage makes Tank-1 particularly suitable for applications requiring rapid heating and minimal heat loss, such as in cold climates or intermittent demand systems. The numerical model demonstrated reliable predictive accuracy, achieving an RMSE of 6.1% for the charging process and 6.8% for the discharging process, thereby validating the model’s reliability. These findings highlight the superior performance of the external heat exchanger configuration for fast and efficient energy storage, particularly for applications in cold climates.

Numerical solutions and stability analysis of hybrid Casson nanofluid flow with MHD and heat transfer effects

AbstractThis research addresses the complex dynamics of hybrid Casson nanoliquids in flow and heat transfer applications, focusing on the interaction between fluid dynamics and thermal phenomena in the presence of magnetohydrodynamics (MHD), viscous dissipation, and Joule heating. The motivation behind this study stems from the need to enhance the efficiency of heat transfer processes in various engineering applications, such as cooling systems and electronic devices, where hybrid nanofluids can offer superior thermal performance compared to conventional fluids. The novelty of this work lies in its comprehensive numerical investigation of a hybrid Casson nanoliquid over a moving permeable surface, incorporating a detailed analysis of MHD effects, viscous dissipation, and Joule heating. Using the Tiwari and Das model to formulate the governing equations, the study employs similarity transformations to convert these equations into a system of ordinary differential equations (ODEs). MATLAB is then used to derive numerical solutions, with a stability analysis ensuring the physical dependability of these solutions. Key findings reveal that the temperature distribution of the hybrid nanoliquid shows a positive correlation with both Prandtl and Hartmann numbers. Additionally, a positive relationship between temperature and the Eckert number is observed. These insights offer valuable guidance for engineers aiming to optimize heat transfer processes using hybrid nanofluids, highlighting their potential for improved thermal management in practical applications.

The influence of lymphatic vessels on nanoparticle distribution and heat transfer within tissue

AbstractThis study analytically investigates the dynamics of nanoparticle transport within a three‐dimensional porous cylinder simulating a lymphatic vessel, without external heat sources. The governing equations and boundary conditions are transformed to yield a system of ordinary differential equations, which are solved numerically using MATLAB built‐in function, bvp4c. Key parameters are visually examined and physically interpreted in relation to temperature, velocity, concentration, and Nusselt number profiles. The study reveals that the distribution of temperature and Nusselt number are maximized by increasing the heat transfer coefficient, whereas NP concentration is increased by decreasing it. Furthermore, the Brownian motion parameter enhances both heat transmission and NP concentration. It is also observed that simpler extravasation into lymphatics decreases tissue nanoparticle levels and heat conduction. Ultimately, optimal intra‐lymphatic nanoparticle distribution pathways are achieved by specifically varying heat transfer and interstitial mass flux patterns. By simulating biological barriers and lymphatic drainage, this model enhances our understanding of the underlying transport mechanisms controlling nanoparticle mobilization.

Mathematical Approach for Directly Solving Air–Water Interfaces in Water Emptying Processes

Emptying processes are operations frequently required in hydraulic installations by water utilities. These processes can result in drops to sub-atmospheric pressure pulses, which may lead to pipeline collapse depending on soil characteristics and the stiffness of a pipe class. One-dimensional mathematical models and 3D computational fluid dynamics (CFD) simulations have been employed to analyse the behaviour of the air–water interface during these events. The numerical resolution of these models is challenging, as 1D models necessitate solving a system of algebraic differential equations. At the same time, 3D CFD simulations can take months to complete depending on the characteristics of the pipeline. This presents a mathematical approach for directly solving air–water interactions in emptying processes involving entrapped air, providing a predictive tool for water utilities. The proposed mathematical approach enables water utilities to predict emptying operations in water pipelines without needing 2D/3D CFD simulations or the resolution of a differential algebraic equations system (1D model). A practical application is demonstrated in a case study of a 350 m long pipe with an internal diameter of 350 mm, investigating the influence of air pocket size, friction factor, polytropic coefficient, pipe diameter, resistance coefficient, and pipe slope. The mathematical approach is validated using an experimental facility that is 7.36 m long, comparing it with 1D mathematical models and 3D CFD simulations. The results confirm that the derived mathematical expression effectively predicts emptying operations in single water installations.

From interacting agents to Boltzmann-Gibbs distribution of money

Abstract We investigate the unbiased model for money exchanges: agents give at random time a dollar to one another (if they have one). Surprisingly, this dynamics eventually leads to a geometric distribution of wealth (shown empirically by Dragulescu and Yakovenko, and rigorously by several follow-up papers). We prove a uniform-in-time propagation of chaos result as the number of agents goes to infinity, which links the stochastic dynamics to a deterministic infinite system of ordinary differential equations. This deterministic description is then analyzed by taking advantage of several entropy–entropy dissipation inequalities and we provide a quantitative almost-exponential rate of convergence toward the equilibrium (geometric distribution) in relative entropy.

Unraveling the role of exercise in cancer suppression: insights from a mathematical model

Recent experimental studies have shown that physical exercise has the potential to suppress tumor progression. Such suppression has been reported to be mediated by the exercise-induced activation of natural killer (NK) cells through the release of IL-6, a cytokine. Aimed at shedding light on how exercise-induced NK cell activation helps in the suppression of cancer, we developed a coarse-grained mathematical model based on a system of ordinary differential equations describing the interaction between IL-6, NK-cells, and tumor cells. The model is then used to study how exercise duration and exercise intensity affect tumor suppression. Our results show that increasing exercise intensity or increasing exercise duration leads to greater and sustained tumor suppression. Furthermore, multi-bout exercise patterns hold promise for improving cancer treatment strategies by adjusting exercise intensity and frequency. Thus, the proposed mathematical model provides insights into the role of exercise in tumor suppression and can be instrumental in guiding future experimental studies, potentially leading to more effective exercise interventions.

Computational Technology for Shell Models of Magnetohydrodynamic Turbulence Constructing

The paper discusses the computational technology for constructing one type of small-scale magnetohydrodynamic turbulence models – shell models. Any such model is a system of ordinary quadratic nonlinear differential equations with constant coefficients. Each phase variable is interpreted in absolute value as a measure of the intensity of one of the fields of the turbulent system in a certain range of spatial scales (scale shell). The equations of any shell model must have several quadratic invariants, which are analogues of conservation laws in ideal magnetohydrodynamics. The derivation of the model equations consists in obtaining such expressions for constant coefficients for which the predetermined quadratic expressions will indeed be invariants. Derivation of these expressions «manually» is quite cumbersome and the likelihood of errors in formula transformations is high. This is especially true for non-local models in which large-scale shells that are distant in size can interact. The novelty and originality of the work lie in the fact that the authors proposed a computational technology that allows one to automate the process of deriving equations for shell models. The technology was implemented using computer algebra methods, which made it possible to obtain parametric classes of models in which the invariance of given quadratic forms is carried out absolutely accurately – in formula form. The determination of the parameter values in the resulting parametric class of models is further carried out by agreement with the measures of the interaction of shells in the model with the probabilities of their interaction in a real physical system. The idea of the described technology and its implementation belong to the authors. Some of its elements were published by the authors earlier, but in this work, for the first time, its systematic description is given for models with complex phase variables and agreement of measures of interaction of shells with probabilities. There have been no similar works by other authors previously. The technology allows you to quickly and accurately generate equations for new non-local turbulence shell models and can be useful to researchers involved in modeling turbulent systems.

A generalized differential scheme for the effective conductivity of percolating microinhomogeneous materials with the Hall effect

In this paper, we propose a self-consistent scheme for the calculation of the components of the effective electrical conductivity tensor. The calculations were fulfilled for a microinhomogeneous material, the components of which have the Hall effect. The presence of the Hall effect leads to appearance of asymmetry of the components of the conductivity tensor and to dependence of these components on the magnitude of the magnetic field applied to the material. Our approach is based on the Generalized Differential Effective Medium (GDEM) method. This method generalizes the classical differential scheme (DEM) for the case of several inclusion types instead of one. In this case, the GDEM scheme leads to a system of matrix differential equations that were solved numerically. This solution was obtained for materials containing spherical or cylindrical inclusions (3D and 2D-problems). In the case of cylindrical inclusions, the results were obtained for inclusions with the symmetry axes orthogonal to the magnetic field. The application of the GDEM method allows us to consider the percolation effect for 2D and 3D-microheterogeneous materials. The results obtained are compared to the experimental data and the calculation results obtained by other self-consistent schemes.