A novel cooling system combining ultra fast cooling rigs with laminar cooling devices was investigated.Based on the different cooling mechanisms,a serial of mathematic models were established to describe the relations...A novel cooling system combining ultra fast cooling rigs with laminar cooling devices was investigated.Based on the different cooling mechanisms,a serial of mathematic models were established to describe the relationship between water flow and spraying pressure and the relationship between water spraying heat flux and layout of nozzles installed on the top and bottom cooling headers.Model parameters were validated by measured data.Heat transfer models including air convection model,heat radiation model and water cooling capacity model were detailedly introduced.In addition,effects on cooling capacity by water temperature and different valve patterns were also presented.Finally,the comparison results from UFC used or not have been provided with respect to temperature evolution and mechanical properties of Q235B steel grade with thickness of 7.8 mm.Since online application of the sophisticated CTC process control system based on these models,run-out table cooling control system has been running stably and reliably to produce resource-saving,low-cost steels with smaller grain size.展开更多
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.展开更多
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.展开更多
An experimental study was performed to determine the characteristics and drying process of mushroom （Lentinus edodes） by 6 different hot-air drying methods namely isothermal drying, uniform raise drying, non-uniform...An experimental study was performed to determine the characteristics and drying process of mushroom （Lentinus edodes） by 6 different hot-air drying methods namely isothermal drying, uniform raise drying, non-uniform raise drying, uniform intermittent drying, non-uniform intermittent drying and combined drying. The chemical composition （dry matter, ash, crude protein, crude fat, total sugars, dietary fiber, and energy）, color parameters （L, a＊, b＊, c＊, and h~） and rehydration capacities were determined. Among all the experiments, non-uniform intermittent drying reached a better comprehensive results due to the higher chemical composition, better color quality associated with high bright （26.381＋5.842）, high color tone （73.670＋2.975）, low chroma （13.349a：3.456） as well as the highest rehydration （453.76% weigh of dried body）. Nine kinds of classical mathematical model were used to obtained moisture data and the Midili-kucuk model can be described by the drying process with the coefficient （R2 ranged from 0.99790 to 0.99967）, chi-square （X2 ranged from 0.00003 to 0.00019） and root mean square error （RMSE ranged from 0.000486 to 0.0012367）.展开更多
The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displa...The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displacement changes of the piston are ignored, even experiment verification is not conducted. Therefore, in view of deficiencies above, a nonlinear mathematical model is established in this paper, including dynamic characteristics of servo valve, nonlinear characteristics of pressure-flow, initial displacement of servo cylinder piston and friction nonlinearity. The transfer function block diagram is built for the hydraulic drive unit closed loop position control, as well as the state equations. Through deriving the time-varying coefficient items matrix and time-varying free items matrix of sensitivity equations respectively, the expression of sensitivity equations based on the nonlinear mathematical model are obtained. According to structure parameters of hydraulic drive unit, working parameters, fluid transmission characteristics and measured friction-velocity curves, the simulation analysis of hydraulic drive unit is completed on the MATLAB/Simulink simulation platform with the displacement step 2 mm, 5 mm and 10 mm, respectively. The simulation results indicate that the developed nonlinear mathematical model is sufficient by comparing the characteristic curves of experimental step response and simulation step response under different constant load. Then, the sensitivity function time-history curves of seventeen parameters are obtained, basing on each state vector time-history curve of step response characteristic. The maximum value of displacement variation percentage and the sum of displacement variation absolute values in the sampling time are both taken as sensitivity indexes. The sensitivity indexes values above are calculated and shown visually in histograms under different working conditions, and change rules are analyzed. Then the sensitivity indexes values of four measurable parameters, such as supply pressure, proportional gain, initial position of servo cylinder piston and load force, are verified experimentally on test platform of hydraulic drive unit, and the experimental research shows that the sensitivity analysis results obtained through simulation are approximate to the test results. This research indicates each parameter sensitivity characteristics of hydraulic drive unit, the performance-affected main parameters and secondary parameters are got under different working conditions, which will provide the theoretical foundation for the control compensation and structure optimization of hydraulic drive unit.展开更多
The available studies in the literature on physical and mathematical modeling of the argon oxygen decarburization (AOD) process of stainless steel have briefly been reviewed. The latest advances made by the author wi...The available studies in the literature on physical and mathematical modeling of the argon oxygen decarburization (AOD) process of stainless steel have briefly been reviewed. The latest advances made by the author with his research group have been summarized. Water modeling was used to investigate the fluid flow and mixing characteristics in the bath of an 18 t AOD vessel, as well as the 'back attack' action of gas jets and its effects on the erosion and wear of the refractory lining, with sufficiently full kinematic similarity. The non rotating and rotating gas jets blown through two annular tuyeres, respectively of straight tube and spiral flat tube type, were employed in the experiments. The geometric similarity ratio between the model and its prototype (including the straight tube type tuyeres) was 1:3. The influences of the gas flow rate, the angle included between the two tuyeres and other operating parameters, and the suitability of the spiral tuyere as a practical application, were examined. These latest studies have clearly and successfully brought to light the fluid flow and mixing characteristics in the bath and the overall features of the back attack phenomena of gas jets during the blowing, and have offered a better understanding of the refining process. Besides, mathematical modeling for the refining process of stainless steel was carried out and a new mathematical model of the process was proposed and developed. The model performs the rate calculations of the refining and the mass and heat balances of the system. Also, the effects of the operating factors, including adding the slag materials, crop ends, and scrap, and alloy agents; the non isothermal conditions; the changes in the amounts of metal and slag during the refining; and other factors were all considered. The model was used to deal with and analyze the austenitic stainless steel making (including ultra low carbon steel) and was tested on data of 32 heats obtained in producing 304 grade steel in an 18 t AOD vessel. The changes in the bath composition and temperature during the refining process with time can be accurately predicted using this model. The model can provide some very useful information and a reliable basis for optimizing the process practice of the refining of stainless steel and control of the process in real time and online.展开更多
Recent advances in modeling of tornadoes and twisters consist of significant achievements in mathematical calculation of occurrence and evolution of a violent F5-class tornado on the Fujita scale, and four-dimensional...Recent advances in modeling of tornadoes and twisters consist of significant achievements in mathematical calculation of occurrence and evolution of a violent F5-class tornado on the Fujita scale, and four-dimensional mathematical modeling of a tornado with the fourth coordinate time multiplied by its characteristic velocity. Such a tornado can arise in a thunderstorm supercell filled with turbulent whirlwinds. A theory of the squall storms is proposed. The squall storm is modeled by running pertur- bation of the temperature inversion on the lower boundary of cloudiness. This perturbation is induced by the action of strong, hurricane winds in the upper and middle troposphere, and looks like a running solitary wave （soliton）; which is developed also in a field of pressure and velocity of a wind. If a soliton of a squall storm gets into the thunderstorm supercell then this soliton is captured by supercell. It leads to additional pressure fall of air inside a storm supercell and stimulate amplification of wind velocity here. As a result, a cyclostrophic balance inside a storm supercell generates a tornado. Comparison of the radial distribution of wind velocity inside a tornado calculated by using the new formulas and equations with radar observations of the wind velocity inside Texas Tornado Dummit in 1995 and inside the 3 May 1999 Oklahoma City Tornado shows good correspondence.展开更多
The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The ...The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The field of displacements is for- mulated using the classical broken line hypothesis and the proposed nonlinear hypothesis that generalizes the classical one. Using both hypotheses, the strains are determined as well as the stresses of each layer. The kinetic energy, the elastic strain energy, and the work of load are also determined. The system of equations of motion is derived using Hamilton＇s principle. Finally, the system of three equations is reduced to one equation of motion, in particular, the Mathieu equation. The Bubnov-Galerkin method is used to solve the system of equations of motion, and the Runge-Kutta method is used to solve the second-order differential equation. Numerical calculations are done for the chosen family of beams. The critical loads, unstable regions, angular frequencies of the beam, and the static and dynamic equilibrium paths are calculated analytically and verified numerically. The results of this study are presented in the forms of figures and tables.展开更多
Inverted pendulums are important objects of theoretical investigation and experiment in the area of control theory and engineering. The researches concentrate on the rigid finite dimensional models which are described...Inverted pendulums are important objects of theoretical investigation and experiment in the area of control theory and engineering. The researches concentrate on the rigid finite dimensional models which are described by ordinary differential equations (ODEs) .Complete rigidity is the approximation of practical models ; Elasticity should be introduced into mathematical models in the analysis of system dynamics and integration of highly precise controller. A new kind of inverted pendulum, elastic inverted pendulum was proposed, and elasticity was considered. Mathematical model was derived from Hamiltonian principle and variational methods, which were formulated by the coupling of partial differential equations (PDE) and ODE. Because of infinite dimensional, system analysis and control of elastic inverted pendulum is more sophisticated than the rigid one.展开更多
A mathematical model was developed to monitor the facilitated transport of germanium(IV) from oxalic acid solutions through a flat sheet supported liquid membrane(FSSLM) containing four trialkylphosphine oxides(Cyanex...A mathematical model was developed to monitor the facilitated transport of germanium(IV) from oxalic acid solutions through a flat sheet supported liquid membrane(FSSLM) containing four trialkylphosphine oxides(Cyanex 923). The FSSLM modeling was based on the extraction constant(Kext) calculated from the liquid-liquid extraction(LLX) modeling. The LLX model presented a reliable calculation of the extraction constant(Kex= 2.057×103 L/mol4). The FSSLM model was solved using Matlab■ software according to extraction constant, Fick’s law, and diffusional principles. The model predicts the overall mass transfer coefficient(Korg) to be 3.84 cm/s. Using this value, diffusion coefficients(Dm) for various Cyanex 923 concentrations of 0.126, 0.252, 0.378, 0.505, 0.631 and 0.757 mol/L are found to be 8.50×10^-4, 4.30×10^-4, 1.87×10^-4, 5.87×10^-5, 2.57×10^-5, 2.09×10^-5 cm2/s, respectively. The results show that the diffusion rate of the current study is approximately more than that of similar FSSLM systems containing Cyanex 923 used to transport various metals. The modeling values are in good agreement with the experimental data, showing the good reliability of the mathematical model.展开更多
The oxidation of polyunsaturated fatty acids （PUFA） in emulsion with stirring and limited oxygen compensation was studied. A mathematical model of diffusion-oxidation was developed considering the mass transfer resi...The oxidation of polyunsaturated fatty acids （PUFA） in emulsion with stirring and limited oxygen compensation was studied. A mathematical model of diffusion-oxidation was developed considering the mass transfer resistance of a gas-liquid boundary, the resistance of the boundary layer from the emulsifier membrane, and the autocatalytic-type autoxidation reaction of PUFA. The dynamic mass transfer coefficient of the emulsifier membrane, k0, was introduced. The model was verified by comparing the predictions of the model with the experi- mental data. The results indicated that the model was in good agreement with the oxygen diffusion and linoleic acid oxidation in the emulsion, and showed good applicability in the prediction of the effect of the emulsifier type on the oxidation of PUFA in the emulsion. It indicated that the oxidation of PUFA in emulsions, with stirring and limited oxygen compensation from the atmosphere, was controlled mostly by mass transfer resistance from the emulsifier membrane.展开更多
The study on the fluid flow, meniscus oscillation, slag entrapment in continuous casting mould was conducted mathematically and experimentally. The results show that the injection of argon into submerged nozzle enhan...The study on the fluid flow, meniscus oscillation, slag entrapment in continuous casting mould was conducted mathematically and experimentally. The results show that the injection of argon into submerged nozzle enhances the meniscus oscillation, thus increases the probability of slag entrapment, and the critical argon blowing flow rate, which will give rise to slag entrapment, is around 10l/min. The trajectory of bubble is affected by the bubble diameter and the molten steel flow, and the bubble diameter is dominant. The bubble with diameter 1.4mm floats fastest with 0.47m/s terminal velocity.展开更多
Objective： The progression of human cancer is characterized by the accumulation of genetic instability. An increasing number of experimental genetic molecular techniques have been used to detect chromosome aberration...Objective： The progression of human cancer is characterized by the accumulation of genetic instability. An increasing number of experimental genetic molecular techniques have been used to detect chromosome aberrations. Previous studies on chromosome abnormalities often focused on identifying the frequent loci of chromosome alterations, but rarely addressed the issue of interrelationship of chromosomal abnormalities. In the last few years, several mathematical models have been employed to construct models of carcinogenesis, in an attempt to identify the time order and cause-and-effect relationship of chromosome aberrations. The principles and applications of these models are reviewed and compared in this paper. Mathematical modeling of carcinogenesis can contribute to our understanding of the molecular genetics of tumor development, and identification of cancer related genes, thus leading to improved clinical practice of cancer.展开更多
Foam injection is a promising solution for control of mobility in oil and gas field exploration and development,including enhanced oil recovery,matrix-acidization treatments,contaminated-aquifer remediation and gas le...Foam injection is a promising solution for control of mobility in oil and gas field exploration and development,including enhanced oil recovery,matrix-acidization treatments,contaminated-aquifer remediation and gas leakage prevention.This study presents a numerical investigation of foam behavior in a porous medium.Fractional flow method is applied to describe steady-state foam displacement in the entrance region.In this model,foam flow for the cases of excluding and including capillary pressure and for two types of gas,nitrogen(N2)and carbon dioxide(CO2)are investigated.Effects of pertinent parameters are also verified.Results indicate that the foam texture strongly governs foam flow in porous media.Required entrance region may be quite different for foam texture to accede local equilibrium,depending on the case and physical properties that are used.According to the fact that the aim of foaming of injected gas is to reduce gas mobility,results show that CO2 is a more proper injecting gas than N2.There are also some ideas presented here on improvement in foam displacement process.This study will provide an insight into future laboratory research and development of full-field foam flow in a porous medium.展开更多
The dynamic performance of a nozzle-flapper servo valve can be affected by several factors such as the disturbance of the input signal,the motion of the armature assembly and the oscillation of the jet force.As the pa...The dynamic performance of a nozzle-flapper servo valve can be affected by several factors such as the disturbance of the input signal,the motion of the armature assembly and the oscillation of the jet force.As the part of vibrating at high frequency,the armature assembly plays a vital role during the operation of the servo valve.In order to accurately predict the transient response of the armature assembly during the vibration,a mathematical model of armature assembly is established based on the distributed parameters method(DPM)and Hamilton principle.The new mathematical model is composed of three main parts,the modal eigenfunction,modal mechanical response expressions of the spring tube and the motion equation of the other armature assembly.After programing,the purpose of using the DPM to predict the dynamic response of different positions located on the armature assembly is achieved.For verifying the validity of the mathematical model,the finite element method(FEM)and classic model(CM)of armature assembly are applicated by commercial software under the same condition.The comparison results prove that the DPM can effectively predict the axial and tangential deflection of the armature assembly different positions which the CM can’t duing to its over-simplification.A certain error is generated when predicting the axial deformation at different heights by DPM,which is caused by an approximate method to simulate the torsion of the spring tube.The comparison results of the spring tube deflection at different vibration frequencies shows that the adaptability of DPM is significantly higher than the classic model,which verify the model is more adaptable for predicting the dynamic response of the armature assembly.展开更多
Optimization of mathematical model of flow field in slab continuous casting mold was performed by means of industrial measurement and mathematical modeling.The rod deflection method was used to quantitatively measure ...Optimization of mathematical model of flow field in slab continuous casting mold was performed by means of industrial measurement and mathematical modeling.The rod deflection method was used to quantitatively measure the velocities near the mold surface at high temperature.The measurement results were compared with the simulation results of three mathematical models at different argon gas flow rates of 6,10 and 14 L min^(−1).The model 1 neglects the mold powder layer,thermal effect and solidified shell.The model 2 only considers the influence of mold powder layer.The model 3 considers the influence of mold powder layer,thermal effect and solidified shell on the flow field.In all three models,the diameter of argon bubbles obeys Rosin-Rammler distribution fitted according to the experimental data of others’previous work.With increasing the argon gas flow rate,the velocity of liquid steel near the mold surface decreases.The model 1 seriously underestimates the shear stress of liquid steel near the mold surface,and its calculation results show higher velocity near the mold surface,lower turbulent kinetic energy and wider distribution of argon gas bubbles in the mold.The simulation results of model 2 only considering the viscous resistance of the mold powder layer to liquid steel makes the velocity near the surface lower than the measurement results obviously.The calculated velocities near the mold surface with model 3 are in best agreement with the measured results,showing the reasonable spatial distribution range of argon bubbles in the mold and the moderate turbulent kinetic energy.In the present conditions,the best argon gas flow rate is 10 L min^(−1) due to the moderate velocity near the mold surface,the appropriate distribution of argon gas bubbles in the mold and the smallest fluctuation amplitude on the mold surface.展开更多
The salt-gradient operation mode used in ion-exchange simulated moving bed chromatography (SMBC) can improve the efficiency of protein separations. A detailed model that takes into account any kind of adsorption/ion-e...The salt-gradient operation mode used in ion-exchange simulated moving bed chromatography (SMBC) can improve the efficiency of protein separations. A detailed model that takes into account any kind of adsorption/ion-exchange equilibrium, salt gradient, size exclusion, mass transfer resistance, and port periodic switching mechanism, was developed to simulate the complex dynamics. The model predictions were verified by the experimental data on upward and downward gradients for protein separations reported in the literature. All design and operating parameters (number, configuration, length and diameter of columns, particle size, switching period, flow rates of feed, raffinate, desorbent and extract, protein concentrations in feed, different salt concentrations in desorbent and feed) can be chosen correctly by numerical simulation. This model can facilitate the design, operation, optimization, control and scale-up of salt-gradient ion-exchange SMBC for protein separations.展开更多
Structure of porous media and fluid distribution in rocks can significantly affect the transport characteristics during the process of microscale tracer flow.To clarify the effect of micro heterogeneity on aqueous tra...Structure of porous media and fluid distribution in rocks can significantly affect the transport characteristics during the process of microscale tracer flow.To clarify the effect of micro heterogeneity on aqueous tracer transport,this paper demonstrates microscopic experiments at pore level and proposes an improved mathematical model for tracer transport.The visualization results show a faster tracer movement into movable water than it into bound water,and quicker occupancy in flowing pores than in storage pores caused by the difference of tracer velocity.Moreover,the proposed mathematical model includes the effects of bound water and flowing porosity by applying interstitial flow velocity expression.The new model also distinguishes flowing and storage pores,accounting for different tracer transport mechanisms(dispersion,diffusion and adsorption)in different types of pores.The resulting analytical solution better matches with tracer production data than the standard model.The residual sum of squares(RSS)from the new model is 0.0005,which is 100 times smaller than the RSS from the standard model.The sensitivity analysis indicates that the dispersion coefficient and flowing porosity shows a negative correlation with the tracer breakthrough time and the increasing slope,whereas the superficial velocity and bound water saturation show a positive correlation.展开更多
The available studies in the literature on mathematical modeling of the vacuum circulation (RH) refining process of molten steel have briefly been reviewed. The latest advances obtained by the author with his research...The available studies in the literature on mathematical modeling of the vacuum circulation (RH) refining process of molten steel have briefly been reviewed. The latest advances obtained by the author with his research group have been summarized. On the basis of the mass and momentum balances in the system, a new mathematical model for decarburization and degassing during the RH and RH KTB refining processes of molten steel was proposed and developed. The refining roles of the three reaction sites, i.e. the up snorkel zone, the droplet group and steel bath in the vacuum vessel, were considered in the model. It was assumed that the mass transfer of reactive components in the molten steel is the rate control step of the refining reactions. And the friction losses and drags of flows in the snorkels and vacuum vessel were all counted. The model was applied to the refining of molten steel in a multifunction RH degasser of 90 t capacity. The decarburization and degassing processes in the degasser under the RH and RH KTB operating conditions were modeled and analyzed using this model. Besides, proceeded from the two resistance mass transfer theory and the mass balance of sulphur in the system, a kinetic model for the desulphurization by powder injection and blowing in the RH refining of molten steel was developed. Modeling and predictions of the process of injecting and blowing the lime based powder flux under assumed operating modes with the different initial contents of sulphur and amounts of powder injected and blown in a RH degasser of 300 t capacity were carried out using the model. It was demonstrated that for the RH and RH KTB refining processes, and the desulphurization by powder injection and blowing in the RH refining, the results predicted by the models were all in good agreement respectively with data from industrial experiments and practice. These models may be expected to offer some useful information and a reliable basis for determining and optimizing the technologies of the RH and RH KTB refining and desulphurization by powder injection and blowing in the RH refining and for controlling the processes.展开更多
基金Project(51034009)supported by the National Natural Science Foundation of China
文摘A novel cooling system combining ultra fast cooling rigs with laminar cooling devices was investigated.Based on the different cooling mechanisms,a serial of mathematic models were established to describe the relationship between water flow and spraying pressure and the relationship between water spraying heat flux and layout of nozzles installed on the top and bottom cooling headers.Model parameters were validated by measured data.Heat transfer models including air convection model,heat radiation model and water cooling capacity model were detailedly introduced.In addition,effects on cooling capacity by water temperature and different valve patterns were also presented.Finally,the comparison results from UFC used or not have been provided with respect to temperature evolution and mechanical properties of Q235B steel grade with thickness of 7.8 mm.Since online application of the sophisticated CTC process control system based on these models,run-out table cooling control system has been running stably and reliably to produce resource-saving,low-cost steels with smaller grain size.
基金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.
文摘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 National High-Tech R&D Program of China(863 Program,2011AA100805-2)the Project from Chongqing Science and Technology Committee(CSTC2011AC1010)supported by the National Natural Science Foundation of China(31271825)
文摘An experimental study was performed to determine the characteristics and drying process of mushroom （Lentinus edodes） by 6 different hot-air drying methods namely isothermal drying, uniform raise drying, non-uniform raise drying, uniform intermittent drying, non-uniform intermittent drying and combined drying. The chemical composition （dry matter, ash, crude protein, crude fat, total sugars, dietary fiber, and energy）, color parameters （L, a＊, b＊, c＊, and h~） and rehydration capacities were determined. Among all the experiments, non-uniform intermittent drying reached a better comprehensive results due to the higher chemical composition, better color quality associated with high bright （26.381＋5.842）, high color tone （73.670＋2.975）, low chroma （13.349a：3.456） as well as the highest rehydration （453.76% weigh of dried body）. Nine kinds of classical mathematical model were used to obtained moisture data and the Midili-kucuk model can be described by the drying process with the coefficient （R2 ranged from 0.99790 to 0.99967）, chi-square （X2 ranged from 0.00003 to 0.00019） and root mean square error （RMSE ranged from 0.000486 to 0.0012367）.
基金Supported by National Key Basic Research Program of China(973 Program,Grant No.2014CB046405)Hebei Provincial Applied Basic Research Program(Grant No.12962147D)National Natural Science Foundation of China(Grant No.51375423)
文摘The previous sensitivity analysis researches are not accurate enough and also have the limited reference value, because those mathematical models are relatively simple and the change of the load and the initial displacement changes of the piston are ignored, even experiment verification is not conducted. Therefore, in view of deficiencies above, a nonlinear mathematical model is established in this paper, including dynamic characteristics of servo valve, nonlinear characteristics of pressure-flow, initial displacement of servo cylinder piston and friction nonlinearity. The transfer function block diagram is built for the hydraulic drive unit closed loop position control, as well as the state equations. Through deriving the time-varying coefficient items matrix and time-varying free items matrix of sensitivity equations respectively, the expression of sensitivity equations based on the nonlinear mathematical model are obtained. According to structure parameters of hydraulic drive unit, working parameters, fluid transmission characteristics and measured friction-velocity curves, the simulation analysis of hydraulic drive unit is completed on the MATLAB/Simulink simulation platform with the displacement step 2 mm, 5 mm and 10 mm, respectively. The simulation results indicate that the developed nonlinear mathematical model is sufficient by comparing the characteristic curves of experimental step response and simulation step response under different constant load. Then, the sensitivity function time-history curves of seventeen parameters are obtained, basing on each state vector time-history curve of step response characteristic. The maximum value of displacement variation percentage and the sum of displacement variation absolute values in the sampling time are both taken as sensitivity indexes. The sensitivity indexes values above are calculated and shown visually in histograms under different working conditions, and change rules are analyzed. Then the sensitivity indexes values of four measurable parameters, such as supply pressure, proportional gain, initial position of servo cylinder piston and load force, are verified experimentally on test platform of hydraulic drive unit, and the experimental research shows that the sensitivity analysis results obtained through simulation are approximate to the test results. This research indicates each parameter sensitivity characteristics of hydraulic drive unit, the performance-affected main parameters and secondary parameters are got under different working conditions, which will provide the theoretical foundation for the control compensation and structure optimization of hydraulic drive unit.
文摘The available studies in the literature on physical and mathematical modeling of the argon oxygen decarburization (AOD) process of stainless steel have briefly been reviewed. The latest advances made by the author with his research group have been summarized. Water modeling was used to investigate the fluid flow and mixing characteristics in the bath of an 18 t AOD vessel, as well as the 'back attack' action of gas jets and its effects on the erosion and wear of the refractory lining, with sufficiently full kinematic similarity. The non rotating and rotating gas jets blown through two annular tuyeres, respectively of straight tube and spiral flat tube type, were employed in the experiments. The geometric similarity ratio between the model and its prototype (including the straight tube type tuyeres) was 1:3. The influences of the gas flow rate, the angle included between the two tuyeres and other operating parameters, and the suitability of the spiral tuyere as a practical application, were examined. These latest studies have clearly and successfully brought to light the fluid flow and mixing characteristics in the bath and the overall features of the back attack phenomena of gas jets during the blowing, and have offered a better understanding of the refining process. Besides, mathematical modeling for the refining process of stainless steel was carried out and a new mathematical model of the process was proposed and developed. The model performs the rate calculations of the refining and the mass and heat balances of the system. Also, the effects of the operating factors, including adding the slag materials, crop ends, and scrap, and alloy agents; the non isothermal conditions; the changes in the amounts of metal and slag during the refining; and other factors were all considered. The model was used to deal with and analyze the austenitic stainless steel making (including ultra low carbon steel) and was tested on data of 32 heats obtained in producing 304 grade steel in an 18 t AOD vessel. The changes in the bath composition and temperature during the refining process with time can be accurately predicted using this model. The model can provide some very useful information and a reliable basis for optimizing the process practice of the refining of stainless steel and control of the process in real time and online.
基金supported by the Presidium of the Russian Academy of Sciences in the framework of the theme 'Disasters'
文摘Recent advances in modeling of tornadoes and twisters consist of significant achievements in mathematical calculation of occurrence and evolution of a violent F5-class tornado on the Fujita scale, and four-dimensional mathematical modeling of a tornado with the fourth coordinate time multiplied by its characteristic velocity. Such a tornado can arise in a thunderstorm supercell filled with turbulent whirlwinds. A theory of the squall storms is proposed. The squall storm is modeled by running pertur- bation of the temperature inversion on the lower boundary of cloudiness. This perturbation is induced by the action of strong, hurricane winds in the upper and middle troposphere, and looks like a running solitary wave （soliton）; which is developed also in a field of pressure and velocity of a wind. If a soliton of a squall storm gets into the thunderstorm supercell then this soliton is captured by supercell. It leads to additional pressure fall of air inside a storm supercell and stimulate amplification of wind velocity here. As a result, a cyclostrophic balance inside a storm supercell generates a tornado. Comparison of the radial distribution of wind velocity inside a tornado calculated by using the new formulas and equations with radar observations of the wind velocity inside Texas Tornado Dummit in 1995 and inside the 3 May 1999 Oklahoma City Tornado shows good correspondence.
基金Project supported by the Ministry of Science and Higher Education of Poland(Nos.04/43/DSPB/0085and 02/21/DSPB/3464)
文摘The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The field of displacements is for- mulated using the classical broken line hypothesis and the proposed nonlinear hypothesis that generalizes the classical one. Using both hypotheses, the strains are determined as well as the stresses of each layer. The kinetic energy, the elastic strain energy, and the work of load are also determined. The system of equations of motion is derived using Hamilton＇s principle. Finally, the system of three equations is reduced to one equation of motion, in particular, the Mathieu equation. The Bubnov-Galerkin method is used to solve the system of equations of motion, and the Runge-Kutta method is used to solve the second-order differential equation. Numerical calculations are done for the chosen family of beams. The critical loads, unstable regions, angular frequencies of the beam, and the static and dynamic equilibrium paths are calculated analytically and verified numerically. The results of this study are presented in the forms of figures and tables.
基金This worie was supported by Ningbo Institute of Technology, Zhejiang University (No. 1051157G301).
文摘Inverted pendulums are important objects of theoretical investigation and experiment in the area of control theory and engineering. The researches concentrate on the rigid finite dimensional models which are described by ordinary differential equations (ODEs) .Complete rigidity is the approximation of practical models ; Elasticity should be introduced into mathematical models in the analysis of system dynamics and integration of highly precise controller. A new kind of inverted pendulum, elastic inverted pendulum was proposed, and elasticity was considered. Mathematical model was derived from Hamiltonian principle and variational methods, which were formulated by the coupling of partial differential equations (PDE) and ODE. Because of infinite dimensional, system analysis and control of elastic inverted pendulum is more sophisticated than the rigid one.
文摘A mathematical model was developed to monitor the facilitated transport of germanium(IV) from oxalic acid solutions through a flat sheet supported liquid membrane(FSSLM) containing four trialkylphosphine oxides(Cyanex 923). The FSSLM modeling was based on the extraction constant(Kext) calculated from the liquid-liquid extraction(LLX) modeling. The LLX model presented a reliable calculation of the extraction constant(Kex= 2.057×103 L/mol4). The FSSLM model was solved using Matlab■ software according to extraction constant, Fick’s law, and diffusional principles. The model predicts the overall mass transfer coefficient(Korg) to be 3.84 cm/s. Using this value, diffusion coefficients(Dm) for various Cyanex 923 concentrations of 0.126, 0.252, 0.378, 0.505, 0.631 and 0.757 mol/L are found to be 8.50×10^-4, 4.30×10^-4, 1.87×10^-4, 5.87×10^-5, 2.57×10^-5, 2.09×10^-5 cm2/s, respectively. The results show that the diffusion rate of the current study is approximately more than that of similar FSSLM systems containing Cyanex 923 used to transport various metals. The modeling values are in good agreement with the experimental data, showing the good reliability of the mathematical model.
基金Supported by the National Natural Science Foundation of China （20401007）.
文摘The oxidation of polyunsaturated fatty acids （PUFA） in emulsion with stirring and limited oxygen compensation was studied. A mathematical model of diffusion-oxidation was developed considering the mass transfer resistance of a gas-liquid boundary, the resistance of the boundary layer from the emulsifier membrane, and the autocatalytic-type autoxidation reaction of PUFA. The dynamic mass transfer coefficient of the emulsifier membrane, k0, was introduced. The model was verified by comparing the predictions of the model with the experi- mental data. The results indicated that the model was in good agreement with the oxygen diffusion and linoleic acid oxidation in the emulsion, and showed good applicability in the prediction of the effect of the emulsifier type on the oxidation of PUFA in the emulsion. It indicated that the oxidation of PUFA in emulsions, with stirring and limited oxygen compensation from the atmosphere, was controlled mostly by mass transfer resistance from the emulsifier membrane.
文摘The study on the fluid flow, meniscus oscillation, slag entrapment in continuous casting mould was conducted mathematically and experimentally. The results show that the injection of argon into submerged nozzle enhances the meniscus oscillation, thus increases the probability of slag entrapment, and the critical argon blowing flow rate, which will give rise to slag entrapment, is around 10l/min. The trajectory of bubble is affected by the bubble diameter and the molten steel flow, and the bubble diameter is dominant. The bubble with diameter 1.4mm floats fastest with 0.47m/s terminal velocity.
基金supported by a grant from the Education Department of Zhejiang Province (No.Y200803235)
文摘Objective： The progression of human cancer is characterized by the accumulation of genetic instability. An increasing number of experimental genetic molecular techniques have been used to detect chromosome aberrations. Previous studies on chromosome abnormalities often focused on identifying the frequent loci of chromosome alterations, but rarely addressed the issue of interrelationship of chromosomal abnormalities. In the last few years, several mathematical models have been employed to construct models of carcinogenesis, in an attempt to identify the time order and cause-and-effect relationship of chromosome aberrations. The principles and applications of these models are reviewed and compared in this paper. Mathematical modeling of carcinogenesis can contribute to our understanding of the molecular genetics of tumor development, and identification of cancer related genes, thus leading to improved clinical practice of cancer.
文摘Foam injection is a promising solution for control of mobility in oil and gas field exploration and development,including enhanced oil recovery,matrix-acidization treatments,contaminated-aquifer remediation and gas leakage prevention.This study presents a numerical investigation of foam behavior in a porous medium.Fractional flow method is applied to describe steady-state foam displacement in the entrance region.In this model,foam flow for the cases of excluding and including capillary pressure and for two types of gas,nitrogen(N2)and carbon dioxide(CO2)are investigated.Effects of pertinent parameters are also verified.Results indicate that the foam texture strongly governs foam flow in porous media.Required entrance region may be quite different for foam texture to accede local equilibrium,depending on the case and physical properties that are used.According to the fact that the aim of foaming of injected gas is to reduce gas mobility,results show that CO2 is a more proper injecting gas than N2.There are also some ideas presented here on improvement in foam displacement process.This study will provide an insight into future laboratory research and development of full-field foam flow in a porous medium.
基金supported by National Natural Science Foundation of China(No.51675119)。
文摘The dynamic performance of a nozzle-flapper servo valve can be affected by several factors such as the disturbance of the input signal,the motion of the armature assembly and the oscillation of the jet force.As the part of vibrating at high frequency,the armature assembly plays a vital role during the operation of the servo valve.In order to accurately predict the transient response of the armature assembly during the vibration,a mathematical model of armature assembly is established based on the distributed parameters method(DPM)and Hamilton principle.The new mathematical model is composed of three main parts,the modal eigenfunction,modal mechanical response expressions of the spring tube and the motion equation of the other armature assembly.After programing,the purpose of using the DPM to predict the dynamic response of different positions located on the armature assembly is achieved.For verifying the validity of the mathematical model,the finite element method(FEM)and classic model(CM)of armature assembly are applicated by commercial software under the same condition.The comparison results prove that the DPM can effectively predict the axial and tangential deflection of the armature assembly different positions which the CM can’t duing to its over-simplification.A certain error is generated when predicting the axial deformation at different heights by DPM,which is caused by an approximate method to simulate the torsion of the spring tube.The comparison results of the spring tube deflection at different vibration frequencies shows that the adaptability of DPM is significantly higher than the classic model,which verify the model is more adaptable for predicting the dynamic response of the armature assembly.
基金supported by the National Natural Science Foundation of China(U1960202)and HBIS Handan Iron and Steel Group Co.,Ltd.
文摘Optimization of mathematical model of flow field in slab continuous casting mold was performed by means of industrial measurement and mathematical modeling.The rod deflection method was used to quantitatively measure the velocities near the mold surface at high temperature.The measurement results were compared with the simulation results of three mathematical models at different argon gas flow rates of 6,10 and 14 L min^(−1).The model 1 neglects the mold powder layer,thermal effect and solidified shell.The model 2 only considers the influence of mold powder layer.The model 3 considers the influence of mold powder layer,thermal effect and solidified shell on the flow field.In all three models,the diameter of argon bubbles obeys Rosin-Rammler distribution fitted according to the experimental data of others’previous work.With increasing the argon gas flow rate,the velocity of liquid steel near the mold surface decreases.The model 1 seriously underestimates the shear stress of liquid steel near the mold surface,and its calculation results show higher velocity near the mold surface,lower turbulent kinetic energy and wider distribution of argon gas bubbles in the mold.The simulation results of model 2 only considering the viscous resistance of the mold powder layer to liquid steel makes the velocity near the surface lower than the measurement results obviously.The calculated velocities near the mold surface with model 3 are in best agreement with the measured results,showing the reasonable spatial distribution range of argon bubbles in the mold and the moderate turbulent kinetic energy.In the present conditions,the best argon gas flow rate is 10 L min^(−1) due to the moderate velocity near the mold surface,the appropriate distribution of argon gas bubbles in the mold and the smallest fluctuation amplitude on the mold surface.
文摘The salt-gradient operation mode used in ion-exchange simulated moving bed chromatography (SMBC) can improve the efficiency of protein separations. A detailed model that takes into account any kind of adsorption/ion-exchange equilibrium, salt gradient, size exclusion, mass transfer resistance, and port periodic switching mechanism, was developed to simulate the complex dynamics. The model predictions were verified by the experimental data on upward and downward gradients for protein separations reported in the literature. All design and operating parameters (number, configuration, length and diameter of columns, particle size, switching period, flow rates of feed, raffinate, desorbent and extract, protein concentrations in feed, different salt concentrations in desorbent and feed) can be chosen correctly by numerical simulation. This model can facilitate the design, operation, optimization, control and scale-up of salt-gradient ion-exchange SMBC for protein separations.
基金funded by National Science and Technology Major Projects(2017ZX05009004,2016ZX05058003)Beijing Natural Science Foundation(2173061)and State Energy Center for Shale Oil Research and Development(G5800-16-ZS-KFNY005).
文摘Structure of porous media and fluid distribution in rocks can significantly affect the transport characteristics during the process of microscale tracer flow.To clarify the effect of micro heterogeneity on aqueous tracer transport,this paper demonstrates microscopic experiments at pore level and proposes an improved mathematical model for tracer transport.The visualization results show a faster tracer movement into movable water than it into bound water,and quicker occupancy in flowing pores than in storage pores caused by the difference of tracer velocity.Moreover,the proposed mathematical model includes the effects of bound water and flowing porosity by applying interstitial flow velocity expression.The new model also distinguishes flowing and storage pores,accounting for different tracer transport mechanisms(dispersion,diffusion and adsorption)in different types of pores.The resulting analytical solution better matches with tracer production data than the standard model.The residual sum of squares(RSS)from the new model is 0.0005,which is 100 times smaller than the RSS from the standard model.The sensitivity analysis indicates that the dispersion coefficient and flowing porosity shows a negative correlation with the tracer breakthrough time and the increasing slope,whereas the superficial velocity and bound water saturation show a positive correlation.
文摘The available studies in the literature on mathematical modeling of the vacuum circulation (RH) refining process of molten steel have briefly been reviewed. The latest advances obtained by the author with his research group have been summarized. On the basis of the mass and momentum balances in the system, a new mathematical model for decarburization and degassing during the RH and RH KTB refining processes of molten steel was proposed and developed. The refining roles of the three reaction sites, i.e. the up snorkel zone, the droplet group and steel bath in the vacuum vessel, were considered in the model. It was assumed that the mass transfer of reactive components in the molten steel is the rate control step of the refining reactions. And the friction losses and drags of flows in the snorkels and vacuum vessel were all counted. The model was applied to the refining of molten steel in a multifunction RH degasser of 90 t capacity. The decarburization and degassing processes in the degasser under the RH and RH KTB operating conditions were modeled and analyzed using this model. Besides, proceeded from the two resistance mass transfer theory and the mass balance of sulphur in the system, a kinetic model for the desulphurization by powder injection and blowing in the RH refining of molten steel was developed. Modeling and predictions of the process of injecting and blowing the lime based powder flux under assumed operating modes with the different initial contents of sulphur and amounts of powder injected and blown in a RH degasser of 300 t capacity were carried out using the model. It was demonstrated that for the RH and RH KTB refining processes, and the desulphurization by powder injection and blowing in the RH refining, the results predicted by the models were all in good agreement respectively with data from industrial experiments and practice. These models may be expected to offer some useful information and a reliable basis for determining and optimizing the technologies of the RH and RH KTB refining and desulphurization by powder injection and blowing in the RH refining and for controlling the processes.
