Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turb...Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.展开更多
The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study o...The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.展开更多
Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishe...Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.展开更多
Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer...Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.展开更多
Spillover of trypanosomiasis parasites from wildlife to domestic livestock and humans remains a major challenge world over.With the disease targeted for elimination by 2030,assessing the impact of control strategies i...Spillover of trypanosomiasis parasites from wildlife to domestic livestock and humans remains a major challenge world over.With the disease targeted for elimination by 2030,assessing the impact of control strategies in communities where there are human-cattle-wildlife interactions is therefore essential.A compartmental framework incorporating tsetse flies,humans,cattle,wildlife and various disease control strategies is developed and analyzed.The reproduction is derived and its sensitivity to different model parameters is investigated.Meanwhile,the optimal control theory is used to identify a combination of control strategies capable of minimizing the infected human and cattle population over time at minimal costs of implementation.The results indicates that tsetse fly mortality rate is strongly and negatively correlated to the reproduction number.It is also established that tsetse fly feeding rate in strongly and positively correlated to the reproduction number.Simulation results indicates that time dependent control strategies can significantly reduce the infections.Overall,the study shows that screening and treatment of humans may not lead to disease elimination.Combining this strategy with other strategies such as screening and treatment of cattle and vector control strategies will result in maximum reduction of tsetse fly population and disease elimination.展开更多
Malaria is a significant global health challenge.This devastating disease continues to affect millions,especially in tropical regions.It is caused by Plasmodium parasites transmitted by female Anopheles mosquitoes.Thi...Malaria is a significant global health challenge.This devastating disease continues to affect millions,especially in tropical regions.It is caused by Plasmodium parasites transmitted by female Anopheles mosquitoes.This study introduces a nonlinear mathematical model for examining the transmission dynamics of malaria,incorporating both human and mosquito populations.We aim to identify the key factors driving the endemic spread of malaria,determine feasible solutions,and provide insights that lead to the development of effective prevention and management strategies.We derive the basic reproductive number employing the next-generation matrix approach and identify the disease-free and endemic equilibrium points.Stability analyses indicate that the disease-free equilibrium is locally and globally stable when the reproductive number is below one,whereas an endemic equilibrium persists when this threshold is exceeded.Sensitivity analysis identifies the most influential mosquito-related parameters,particularly the bite rate and mosquito mortality,in controlling the spread of malaria.Furthermore,we extend our model to include a treatment compartment and three disease-preventive control variables such as antimalaria drug treatments,use of larvicides,and the use of insecticide-treated mosquito nets for optimal control analysis.The results show that optimal use of mosquito nets,use of larvicides for mosquito population control,and treatment can lower the basic reproduction number and control malaria transmission with minimal intervention costs.The analysis of disease control strategies and findings offers valuable information for policymakers in designing cost-effective strategies to combat malaria.展开更多
In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase dra...In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.展开更多
The increasing demand due to development and advancement in every field of life has caused the depletion of fossil fuels.This depleting fossil fuel reserve throughout the world has enforced to get energy from alternat...The increasing demand due to development and advancement in every field of life has caused the depletion of fossil fuels.This depleting fossil fuel reserve throughout the world has enforced to get energy from alternative/renewable sources.One of the economicalways to get energy is through the utilization of solar ponds.In this study,a mathematical model of a salt gradient solar pond under the Islamabad climatic conditions has been analyzed for the first time.The model uses a one-dimensional finite difference explicit method for optimization of different zone thicknesses.The model depicts that NCZ(Non-Convective Zone)thickness has a significant effect on LCZ(Lower Convective Zone)temperature and should be kept less than 1.7mfor the optimal temperature.It is also observed that for long-termoperation of a solar pond,heat should be extracted by keeping the mass flowrate of 17.3 kg/m^(2)/day.Themodel also suggests that when the bottom reflectivity is about 0.3,then only 24%of the radiation is absorbed in the pond.展开更多
Improving the specific,technical,economic,and environmental characteristics of piston engines(ICE)operating on alternative gaseous fuels is a pressing task for the energy and mechanical engineering industries.The aim ...Improving the specific,technical,economic,and environmental characteristics of piston engines(ICE)operating on alternative gaseous fuels is a pressing task for the energy and mechanical engineering industries.The aim of the study was to optimize the parameters of the ICE working cycle after replacing the base fuel(propane-butane blend)with syngas from wood sawdust to improve its technical and economic performance based on mathematical modeling.The modeling results were verified through experimental studies(differences for key parameters did not exceed 4.0%).The object of the study was an electric generator based on a single-cylinder spark ignition engine with a power of 1 kW.The article describes the main approaches to creating a mathematical model of the engine working cycle,a test bench for modeling verification,physicochemical properties of the base fuel(propane-butane blend),and laboratory syngas.It was shown that replacing the fuel from a propane-butane blend to laboratory syngas caused a decrease in engine efficiency to 33%(the efficiency of the base ICE was 0.179 vs.the efficiency of 0.119 for the converted ICE for the 0.59 kW power mode).Engine efficiency was chosen as the key criterion for optimizing the working cycle.As a result of optimization,the efficiency of the converted syngas engine was 6.1%higher than that of the base engine running on the propane-butane blend,and the power drop did not exceed 8.0%.Thus,careful fine-tuning of the working cycle parameters allows increasing the technical and economic characteristics of the syngas engine to the level of ICEs running on traditional types of fuel.展开更多
The work considers the problem of gas hydrate dissociation in a porous medium using the two-term Forchheimer law,corresponding to high flow rates of reservoir fluids.Such rates can arise during the decomposition of ga...The work considers the problem of gas hydrate dissociation in a porous medium using the two-term Forchheimer law,corresponding to high flow rates of reservoir fluids.Such rates can arise during the decomposition of gas hydrates,since a large amount of gas is released.Intensive emissions of gases from the earth’s interior are observed on the ocean floor.They are also associated with a large number of subvertical geological structures under the ocean floor,coming to the surface in the formof local ring funnels(pockmarks).Many similar objects have also been found on land.Particular interest in this problemis caused by climate threats associated with the release of greenhouse gases.The movement of gas released fromthe hydrate to the breakthrough channel is similar to the gas inflow to the well(without hydrate),which is usually described by a two-term filtration law.In this work,a mathematical model of gas hydrate dissociation with a nonlinear Forchheimer-type law ofmotion is developed.The systemis split in two blocks by physical processes,taking into account the quadratic correction to the velocity in the filtration law.The first block is responsible for the convective transfer of saturation parameters in the model,water,gas and hydrate saturations are taken into account.The second block corresponds to the equation of dissipative piezoconductivity with a different number of thermodynamic degrees of freedom,taking into account heat and mass transfer in a porous medium.The performed splitting allows using explicit-implicit difference schemes when solving problems and avoiding strong refinement of the step in time and space.For numerical modeling,the support operator method is used,which makes it possible to discretize partial differential equations on irregular grids,which allows taking into account the complex geometry and lithology of the reservoir.A difference scheme based on the support operator method is developed,which,due to the mutually consistent approximation of vector analysis operations(divergence and gradient),allows to take into account the various flux laws between adjacent grid cells,including quadratic corrections to the velocity.Based on the developed numerical algorithms and their program implementations,calculations of gas hydrate dissociation are performed both in a reservoir of simple geometric structure and in a heterogeneous reservoir of complex configuration.The results obtained correspond to the physics of the processes under consideration.展开更多
The mathematical model of 4He quantum interferometer gyroscope is presented. The model includes the driven equation, the current equation and the position equation. Therefore, it can sufficiently describe the gyro- sc...The mathematical model of 4He quantum interferometer gyroscope is presented. The model includes the driven equation, the current equation and the position equation. Therefore, it can sufficiently describe the gyro- scope system. The driven equation shows the thermally driven gyroscope can work for a long time but the pres- sure driven one cannot. From the current equation, the superfluid currents passing through the weak link contain the AC currents which show the rotation flux, and other currents caused by drive. As shown in the position equa- tion, the displacement of diaphragm is the only detectable parameter in the gyroscope system. The model is tested by the simulations based on experimental parameters, and can be used to research performance of the gyroscope and analyse the gyroscope error.展开更多
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wa...Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.展开更多
Based on the combined hydraulic calculation for the eastern network region at the Pearl River estuary and several outlets to the Lingdingyang Bay, the sediment calculation modelling was introduced in the establishment...Based on the combined hydraulic calculation for the eastern network region at the Pearl River estuary and several outlets to the Lingdingyang Bay, the sediment calculation modelling was introduced in the establishment of the sediment mathematical model for Lingdingyang Bay and the eastern region with one and two dimensional flow calculation. Model adjustment and verification were performed in conjunction with field data. The simulated results coincide well with measured data.In addition the model is applied to predict the shore-line planning scheme of Lingdingyang Bay.The theoretical criterion is provided for the shore line plan in the model.And a new mathematical simulated method is put out to research the planning engineering concerned with one-dimensional net rivers and two-dimensional estuary.展开更多
The most prospective method for certain structural failures and damages that cannot employ redundancy is self-repairing techniques, to ensure especially the maximum flight safety. Based on the characters of self-repai...The most prospective method for certain structural failures and damages that cannot employ redundancy is self-repairing techniques, to ensure especially the maximum flight safety. Based on the characters of self-repairing aircraft, this paper states some basic assumptions of the self-repairing aircraft, and puts forward some special new conceptions concerning the self-repairing aircraft: control input, operating input, command input, repair input and operating and control factor as well as their relationships. Thus it provides a simple and reliable mathematical model structure for the research on the self-repairing control of the aircraft.展开更多
Mathematical models of tire-longitudinal road adhesion for use in the study of road vehicle dynamics are set up so as to express the relations of longitudinal adhesion coefficients with the slip ratio. They perfect th...Mathematical models of tire-longitudinal road adhesion for use in the study of road vehicle dynamics are set up so as to express the relations of longitudinal adhesion coefficients with the slip ratio. They perfect the Pacejka's models in practical use by taking into account the influences of all essential parameters such as road surface condition. vehicle velocity. slip angle. vertical load and slip ratio on the longitudinal adhesion coefficients. The new models are more comprehensive more concise. simpler and more convenient in application in all kinds of simulations of car dynamics in various sorts of braking modes.展开更多
To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of ...To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.展开更多
Mathematical models of tire-lateral mad adhesion for use in mad vehicle dynamics studies are set up to express the relations of adhesion coefficients with slip ratio in lateral direction.The models of tire-lateral mad...Mathematical models of tire-lateral mad adhesion for use in mad vehicle dynamics studies are set up to express the relations of adhesion coefficients with slip ratio in lateral direction.The models of tire-lateral mad adhesion revolutionize the Pacejka's model in concept and therefore make it possible for applications in vehicle dynamics studies by the expression of lateral adhesion coefficient as a function of wheel slip ratio,instead of the wheel slip angle,taking into account in the mean time the influences of mad surface condition, vehicle velocity,vertical load,tire slip angle,and wheel camber angle.展开更多
AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the ra...AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.展开更多
In the chemical vapor deposition（CVD） process of C/C composites,the dynamics and mechanism of precursor gas flowing behavior were analyzed mathematically,in which the precursor gas was infiltrated by the pressure di...In the chemical vapor deposition（CVD） process of C/C composites,the dynamics and mechanism of precursor gas flowing behavior were analyzed mathematically,in which the precursor gas was infiltrated by the pressure difference of the gas flowing through felt.Differential equations were educed which characterized the relations among the pressure inside the felt,the pressure outside the felt of the precursor gas and the porosity of the felt as a function of CVD duration.The gas residence time during the infiltration process through the felt was obtained from the differential equations.The numerical verification is in good agreement with the practical process,indicating the good reliability of the current mathematical model.展开更多
According to the research theory of improved black oil simulator, a practical mathematical model for CO2 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas...According to the research theory of improved black oil simulator, a practical mathematical model for CO2 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance（fast,） streamline method is employed to solve this model as an alternative to traditional finite（difference） （methods.） Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method （BEM）, an explicit total variation diminishing （TVD） method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less（information） need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an（effective） tool for early stage screening and reservoir dynamic management of the CO2（miscible） oil field.展开更多
文摘Based on the tensor analysis of water-sediment two-phase how, the basic model equations for clear water flow and sediment-laden flow are deduced in the general curve coordinates for natural water variable-density turbulent how. Furthermore, corresponding boundary conditions are also presented in connection with the composition and movement of non-uniform bed material. The theoretical results are applied to the calculation of the float open caisson in the construction period and good results are obtained.
基金supported by the National Natural Science Foundation of China(11722104,11671150)supported by the National Natural Science Foundation of China(11571280,11331005)+3 种基金supported by the National Natural Science Foundation of China(11331005,11771150)by GDUPS(2016)the Fundamental Research Funds for the Central Universities of China(D2172260)FANEDD No.201315
文摘The two-phase flow models are commonly used in industrial applications, such as nuclear, power, chemical-process, oil-and-gas, cryogenics, bio-medical, micro-technology and so on. This is a survey paper on the study of compressible nonconservative two-fluid model, drift-flux model and viscous liquid-gas two-phase flow model. We give the research developments of these three two-phase flow models, respectively. In the last part, we give some open problems about the above models.
基金funded by the National Key R&D Program of China,China(Grant No.2023YFB4005500)National Natural Science Foundation of China,China(Grant Nos.52379113 and 52379114).
文摘Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications,yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions.This study establishes a pore-scale numerical framework integratingMonte Carlo-generated fractal porousmedia with Volume of Fluid(VOF)simulations to unravel the coupling among pore distribution characterized by fractal dimension(Df),flow dynamics,and displacement efficiency.A pore-scale model based on the computed tomography(CT)microstructure of Berea sandstone is established,and the simulation results are compared with experimental data.Good agreement is found in phase distribution,breakthrough behavior,and flow path morphology,confirming the reliability of the numerical simulation method.Ten fractal porous media models with Df ranging from 1.25~1.7 were constructed using a Monte-Carlo approach.The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05-5m/s,inwhich the time-resolved two-phase distribution patternswere systematically recorded.The results reveal that smaller fractal dimensions(Df=1.25~1.45)accelerate fingering breakthrough(peak velocity is 1.73 m/s at Df=1.45)due to a bimodal pore size distribution dominated by narrow channels.Increasing Df amplifies vorticity generation by about 3 times(eddy viscosity is 0.033 Pa⋅s at Df=1.7)through reduced interfacial curvature,while tortuosity-driven pressure differentials transition from sharp increases(0.4~6.3 Pa at Df=1.25~1.3)to inertial plateaus(4.8 Pa at Df=1.7).A nonlinear increase in equilibrium gas volume fraction(fav=0.692 at Df=1.7)emerges from residual gas saturation and turbulence-enhanced dispersion.This behavior is further modulated by flow velocity,with fav peaking at 0.72 under capillary-dominated conditions(0.05 m/s),but decreasing to 0.65 in the inertial regime(0.5 m/s).The work quantitatively links fractal topology to multiphase flow regimes,demonstrating the critical role of Df in governing preferential pathways,energy dissipation,and phase distribution.
基金supported by the National Natural Science Foundation of China(52104049)the Young Elite Scientist Sponsorship Program by Beijing Association for Science and Technology(BYESS2023262)。
文摘Polymer flooding is an important means of improving oil recovery and is widely used in Daqing,Xinjiang,and Shengli oilfields,China.Different from conventional injection media such as water and gas,viscoelastic polymer solutions exhibit non-Newtonian and nonlinear flow behavior including shear thinning and shear thickening,polymer convection,diffusion,adsorption,retention,inaccessible pore volume,and reduced effective permeability.However,available well test model of polymer flooding wells generally simplifies these characteristics on pressure transient response,which may lead to inaccurate results.This work proposes a novel two-phase numerical well test model to better describe the polymer viscoelasticity and nonlinear flow behavior.Different influence factors that related to near-well blockage during polymer flooding process,including the degree of blockage(inner zone permeability),the extent of blockage(composite radius),and polymer flooding front radius are explored to investigate these impacts on bottom hole pressure responses.Results show that polymer viscoelasticity has a significant impact on the transitional flow segment of type curves,and the effects of near-well formation blockage and polymer concentration distribution on well test curves are very similar.Thus,to accurately interpret the degree of near-well blockage in injection wells,it is essential to first eliminate the influence of polymer viscoelasticity.Finally,a field case is comprehensively analyzed and discussed to illustrate the applicability of the proposed model.
文摘Spillover of trypanosomiasis parasites from wildlife to domestic livestock and humans remains a major challenge world over.With the disease targeted for elimination by 2030,assessing the impact of control strategies in communities where there are human-cattle-wildlife interactions is therefore essential.A compartmental framework incorporating tsetse flies,humans,cattle,wildlife and various disease control strategies is developed and analyzed.The reproduction is derived and its sensitivity to different model parameters is investigated.Meanwhile,the optimal control theory is used to identify a combination of control strategies capable of minimizing the infected human and cattle population over time at minimal costs of implementation.The results indicates that tsetse fly mortality rate is strongly and negatively correlated to the reproduction number.It is also established that tsetse fly feeding rate in strongly and positively correlated to the reproduction number.Simulation results indicates that time dependent control strategies can significantly reduce the infections.Overall,the study shows that screening and treatment of humans may not lead to disease elimination.Combining this strategy with other strategies such as screening and treatment of cattle and vector control strategies will result in maximum reduction of tsetse fly population and disease elimination.
基金supported by the Deanship of Scientific Research,Vice Presidency for Graduate Studies and Scientific Research,King Faisal University,Saudi Arabia[Grant No.KFU252959].
文摘Malaria is a significant global health challenge.This devastating disease continues to affect millions,especially in tropical regions.It is caused by Plasmodium parasites transmitted by female Anopheles mosquitoes.This study introduces a nonlinear mathematical model for examining the transmission dynamics of malaria,incorporating both human and mosquito populations.We aim to identify the key factors driving the endemic spread of malaria,determine feasible solutions,and provide insights that lead to the development of effective prevention and management strategies.We derive the basic reproductive number employing the next-generation matrix approach and identify the disease-free and endemic equilibrium points.Stability analyses indicate that the disease-free equilibrium is locally and globally stable when the reproductive number is below one,whereas an endemic equilibrium persists when this threshold is exceeded.Sensitivity analysis identifies the most influential mosquito-related parameters,particularly the bite rate and mosquito mortality,in controlling the spread of malaria.Furthermore,we extend our model to include a treatment compartment and three disease-preventive control variables such as antimalaria drug treatments,use of larvicides,and the use of insecticide-treated mosquito nets for optimal control analysis.The results show that optimal use of mosquito nets,use of larvicides for mosquito population control,and treatment can lower the basic reproduction number and control malaria transmission with minimal intervention costs.The analysis of disease control strategies and findings offers valuable information for policymakers in designing cost-effective strategies to combat malaria.
基金funded by Sichuan Natural Science Foundation Outstanding Youth Science Foundation(No.2024NSFJQ0012)Key project of Regional Innovation and Development Joint Fund of National Natural Science Foundation(No.U23A20669)Sichuan Science and Technology Program(2022ZDZX0041).
文摘In multiphase pumps transporting gas-liquid two-phase flows,the high-speed rotation of the impeller induces complex deformations in bubble shapes within the flow domain,making the prediction of gasliquid two-phase drag forces highly challenging in numerical simulations.To achieve precise prediction of the drag forces on irregular bubbles within multiphase pumps,this study modifies the existing bubble drag force model and applies the revised model to the prediction of gas-liquid two-phase flow within multiphase pumps.The research findings indicate that the modified drag force model significantly enhances the accuracy of predicting flow characteristics within the pump,particularly under high gas volume fraction conditions.The simulation results for gas phase distribution and vorticity exhibit strong agreement with experimental data.The modified drag model better captures the accumulation of the gas phase at the suction side of the impeller outlet.It also accurately predicts the vortex characteristics induced by bubble backflow from the trailing edges of the diffuser.Additionally,the adjustment of the drag coefficient enhances the model’s ability to represent local flow field characteristics,thereby optimizing the performance simulation methods of multiphase pumps.Compared to traditional drag force models,the modified model reduces prediction errors in head and efficiency by 36.4%and 27.5%,respectively.These results provide important theoretical foundations and model support for improving the accuracy of gas-liquid two-phase flow simulations and optimizing the design of multiphase pumps under high gas volume fraction conditions.
文摘The increasing demand due to development and advancement in every field of life has caused the depletion of fossil fuels.This depleting fossil fuel reserve throughout the world has enforced to get energy from alternative/renewable sources.One of the economicalways to get energy is through the utilization of solar ponds.In this study,a mathematical model of a salt gradient solar pond under the Islamabad climatic conditions has been analyzed for the first time.The model uses a one-dimensional finite difference explicit method for optimization of different zone thicknesses.The model depicts that NCZ(Non-Convective Zone)thickness has a significant effect on LCZ(Lower Convective Zone)temperature and should be kept less than 1.7mfor the optimal temperature.It is also observed that for long-termoperation of a solar pond,heat should be extracted by keeping the mass flowrate of 17.3 kg/m^(2)/day.Themodel also suggests that when the bottom reflectivity is about 0.3,then only 24%of the radiation is absorbed in the pond.
基金the Ministry of Science and Higher Education of the Russian Federation(Ural Federal University Program of Development within the Priority-2030 Program)is gratefully acknowledged.
文摘Improving the specific,technical,economic,and environmental characteristics of piston engines(ICE)operating on alternative gaseous fuels is a pressing task for the energy and mechanical engineering industries.The aim of the study was to optimize the parameters of the ICE working cycle after replacing the base fuel(propane-butane blend)with syngas from wood sawdust to improve its technical and economic performance based on mathematical modeling.The modeling results were verified through experimental studies(differences for key parameters did not exceed 4.0%).The object of the study was an electric generator based on a single-cylinder spark ignition engine with a power of 1 kW.The article describes the main approaches to creating a mathematical model of the engine working cycle,a test bench for modeling verification,physicochemical properties of the base fuel(propane-butane blend),and laboratory syngas.It was shown that replacing the fuel from a propane-butane blend to laboratory syngas caused a decrease in engine efficiency to 33%(the efficiency of the base ICE was 0.179 vs.the efficiency of 0.119 for the converted ICE for the 0.59 kW power mode).Engine efficiency was chosen as the key criterion for optimizing the working cycle.As a result of optimization,the efficiency of the converted syngas engine was 6.1%higher than that of the base engine running on the propane-butane blend,and the power drop did not exceed 8.0%.Thus,careful fine-tuning of the working cycle parameters allows increasing the technical and economic characteristics of the syngas engine to the level of ICEs running on traditional types of fuel.
基金the framework of the state assignment of Keldysh Institute of Applied Mathematics of RAS(Project No.125020701776-0)the Ministry of Education and Science of Russia for IO RAS(Project No.FMWE-2024-0018).
文摘The work considers the problem of gas hydrate dissociation in a porous medium using the two-term Forchheimer law,corresponding to high flow rates of reservoir fluids.Such rates can arise during the decomposition of gas hydrates,since a large amount of gas is released.Intensive emissions of gases from the earth’s interior are observed on the ocean floor.They are also associated with a large number of subvertical geological structures under the ocean floor,coming to the surface in the formof local ring funnels(pockmarks).Many similar objects have also been found on land.Particular interest in this problemis caused by climate threats associated with the release of greenhouse gases.The movement of gas released fromthe hydrate to the breakthrough channel is similar to the gas inflow to the well(without hydrate),which is usually described by a two-term filtration law.In this work,a mathematical model of gas hydrate dissociation with a nonlinear Forchheimer-type law ofmotion is developed.The systemis split in two blocks by physical processes,taking into account the quadratic correction to the velocity in the filtration law.The first block is responsible for the convective transfer of saturation parameters in the model,water,gas and hydrate saturations are taken into account.The second block corresponds to the equation of dissipative piezoconductivity with a different number of thermodynamic degrees of freedom,taking into account heat and mass transfer in a porous medium.The performed splitting allows using explicit-implicit difference schemes when solving problems and avoiding strong refinement of the step in time and space.For numerical modeling,the support operator method is used,which makes it possible to discretize partial differential equations on irregular grids,which allows taking into account the complex geometry and lithology of the reservoir.A difference scheme based on the support operator method is developed,which,due to the mutually consistent approximation of vector analysis operations(divergence and gradient),allows to take into account the various flux laws between adjacent grid cells,including quadratic corrections to the velocity.Based on the developed numerical algorithms and their program implementations,calculations of gas hydrate dissociation are performed both in a reservoir of simple geometric structure and in a heterogeneous reservoir of complex configuration.The results obtained correspond to the physics of the processes under consideration.
基金Supported by the National Natural Science Foundation of China(61074162)the Ph.D.Program Foundation of Ministry of Education of China(200802870011)~~
文摘The mathematical model of 4He quantum interferometer gyroscope is presented. The model includes the driven equation, the current equation and the position equation. Therefore, it can sufficiently describe the gyro- scope system. The driven equation shows the thermally driven gyroscope can work for a long time but the pres- sure driven one cannot. From the current equation, the superfluid currents passing through the weak link contain the AC currents which show the rotation flux, and other currents caused by drive. As shown in the position equa- tion, the displacement of diaphragm is the only detectable parameter in the gyroscope system. The model is tested by the simulations based on experimental parameters, and can be used to research performance of the gyroscope and analyse the gyroscope error.
基金National Natural Science Foundation (Project number 40604013).
文摘Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement.
文摘Based on the combined hydraulic calculation for the eastern network region at the Pearl River estuary and several outlets to the Lingdingyang Bay, the sediment calculation modelling was introduced in the establishment of the sediment mathematical model for Lingdingyang Bay and the eastern region with one and two dimensional flow calculation. Model adjustment and verification were performed in conjunction with field data. The simulated results coincide well with measured data.In addition the model is applied to predict the shore-line planning scheme of Lingdingyang Bay.The theoretical criterion is provided for the shore line plan in the model.And a new mathematical simulated method is put out to research the planning engineering concerned with one-dimensional net rivers and two-dimensional estuary.
文摘The most prospective method for certain structural failures and damages that cannot employ redundancy is self-repairing techniques, to ensure especially the maximum flight safety. Based on the characters of self-repairing aircraft, this paper states some basic assumptions of the self-repairing aircraft, and puts forward some special new conceptions concerning the self-repairing aircraft: control input, operating input, command input, repair input and operating and control factor as well as their relationships. Thus it provides a simple and reliable mathematical model structure for the research on the self-repairing control of the aircraft.
文摘Mathematical models of tire-longitudinal road adhesion for use in the study of road vehicle dynamics are set up so as to express the relations of longitudinal adhesion coefficients with the slip ratio. They perfect the Pacejka's models in practical use by taking into account the influences of all essential parameters such as road surface condition. vehicle velocity. slip angle. vertical load and slip ratio on the longitudinal adhesion coefficients. The new models are more comprehensive more concise. simpler and more convenient in application in all kinds of simulations of car dynamics in various sorts of braking modes.
文摘To improve the accuracy of the conventional finite-difference method, finitedifference numerical modeling methods of any even-order accuracy are recommended. We introduce any even-order accuracy difference schemes of any-order derivatives derived from Taylor series expansion. Then, a finite-difference numerical modeling method with any evenorder accuracy is utilized to simulate seismic wave propagation in two-phase anisotropic media. Results indicate that modeling accuracy improves with the increase of difference accuracy order number. It is essential to find the optimal order number, grid size, and time step to balance modeling precision and computational complexity. Four kinds of waves, static mode in the source point, SV wave cusps, reflection and transmission waves are observed in two-phase anisotropic media through modeling.
文摘Mathematical models of tire-lateral mad adhesion for use in mad vehicle dynamics studies are set up to express the relations of adhesion coefficients with slip ratio in lateral direction.The models of tire-lateral mad adhesion revolutionize the Pacejka's model in concept and therefore make it possible for applications in vehicle dynamics studies by the expression of lateral adhesion coefficient as a function of wheel slip ratio,instead of the wheel slip angle,taking into account in the mean time the influences of mad surface condition, vehicle velocity,vertical load,tire slip angle,and wheel camber angle.
基金supported by the National Natural Science Foundation of China(Grant Nos.41404101,41174114,41274130,and 41404102)
文摘AVO forward modeling is based on two-phase medium theory and is considered an effective method for describing reservoir rocks and fluids. However, the method depends on the input matrix mineral bulk modulus and the rationality of the two-phase medium model. We used the matrix mineral bulk modulus inversion method and multiple constraints to obtain a two-phase medium model with physical meaning. The proposed method guarantees the reliability of the obtained AVO characteristicsin two-phase media. By the comparative analysis of different lithology of the core sample, the advantages and accuracy of the inversion method can be illustrated. Also, the inversion method can be applied in LH area, and the AVO characteristics can be obtained when the porosity, fluid saturation, and other important lithology parameters are changed. In particular, the reflection coefficient amplitude difference between the fast P wave and S wave as a function of porosity at the same incidence angle, and the difference in the incidence angle threshold can be used to decipher porosity.
基金Projects (50702078,50874123) supported by the National Natural Science Foundation of ChinaProject (2009AA03Z536) supported by the National High-tech Research and Development Program of China+1 种基金Project (2011CB606306) supported by the National Research Program of ChinaProject supported by the Program for New Century Excellent Talents in University of China
文摘In the chemical vapor deposition（CVD） process of C/C composites,the dynamics and mechanism of precursor gas flowing behavior were analyzed mathematically,in which the precursor gas was infiltrated by the pressure difference of the gas flowing through felt.Differential equations were educed which characterized the relations among the pressure inside the felt,the pressure outside the felt of the precursor gas and the porosity of the felt as a function of CVD duration.The gas residence time during the infiltration process through the felt was obtained from the differential equations.The numerical verification is in good agreement with the practical process,indicating the good reliability of the current mathematical model.
文摘According to the research theory of improved black oil simulator, a practical mathematical model for CO2 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance（fast,） streamline method is employed to solve this model as an alternative to traditional finite（difference） （methods.） Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method （BEM）, an explicit total variation diminishing （TVD） method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less（information） need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an（effective） tool for early stage screening and reservoir dynamic management of the CO2（miscible） oil field.
