Diesel engines have proven over the years important in terms of efficiency and fuel consumption to power generation ratio. Many research works show the potential of biodiesel as a substitute for conventional gasoil. M...Diesel engines have proven over the years important in terms of efficiency and fuel consumption to power generation ratio. Many research works show the potential of biodiesel as a substitute for conventional gasoil. Mainly, previous and recent researches have focused on experimental investigation of diesel engine performance fuelled by biodiesel. Researches on the mathematical description of diesel engine process running on biodiesel are scarce, and mostly about chemical and thermodynamic description of the combustion process of biodiesel rather than performance studies. This work describes a numerical investigation on the performance analysis of a diesel engine fuelled by palm oil biodiesel. The numerical investigation was made using a semi empirical 0D model based on Wiebe’s and Watson’s model which was implemented via the open access numerical calculation software Scilab. The model was validated first by comparing with experimental pressure and performance data of a one cylinder engine at rated speed and secondly by comparing with a six cylinders engine performance data at various crankshaft rotational speeds. Simulations were then made to analyze the engine performance when running on biodiesel. The calculations were made at constant combustion duration and constant coefficient of excess air. Results showed that the model matches the overall experimental data, such as the power output and peak cylinder pressure. The ignition delay was somehow underestimated by the model for the first experiment, which caused a slight gap on in cylinder pressure curve, whereas it predicted the average ignition delay fairly well for the second set of validation. The simulations of engine performance when running on biodiesel confirmed results obtained in previous experimental researches on biodiesel. The model will be further investigated for engine control when shifting to biodiesel fuel.展开更多
Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an applic...Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.展开更多
A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during an...A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during and after electrical discharges are verified in the experiments. A modified 0D internal ballistic model is established to simulate the ETC launch. According to the measured pressure and electrical parameters, a transient burning rate law including the influence of EGGR coefficient by electric power and pressure gradient(dp/dt) is added into the model. The EGGR coefficient of 4/7 high-nitrogen solid propellant is equal to 0.005 MW-1. Both simulated breech pressure and projectile muzzle velocity accord with the experimental results well. Compared with Woodley's modified burning rate law, the breech pressure curves acquired by the transient burning rate law are more consistent with test results. Based on the parameters calculated in the model, the relationship among propellant burning rate, pressure gradient(dp/dt) and electric power is analyzed. Depending on the transient burning rate law and experimental data, the burning of solid propellant under the condition of plasma is described more accurately.展开更多
According to the physics of tokamak start-up,this study constructs a zero-dimensional(0D)model applicable to electron cyclotron(EC)wave assisted start-up in NCST spherical torus(spherical tokamak)and CN-H1 stellarator...According to the physics of tokamak start-up,this study constructs a zero-dimensional(0D)model applicable to electron cyclotron(EC)wave assisted start-up in NCST spherical torus(spherical tokamak)and CN-H1 stellarators.Using the constructed 0D model,the results obtained in this study under the same conditions are compared and validated against reference results for pure hydrogen plasma start-up in tokamak.The results are in good agreement,especially regarding electron temperature,ion temperature and plasma current.In the presence of finite Ohmic electric field in the spherical tokamak,a study on the EC wave assisted start-up of the NCST plasma at frequency of 28 GHz is conducted.The impact of the vertical magnetic field B_(v)on EC wave assisted start-up,the relationship between EC wave injection power P_(inj),Ohmic electric field E,and initial hydrogen atom density n_(H0)are explored separately.It is found that under conditions of Ohmic electric field lower than ITER(~0.3 V m^(-1)),EC wave can expand the operational space to achieve better plasma parameters.Simulating the process of28 GHz EC wave start-up in the CN-H1 stellarator plasma,the plasma current in the zerodimensional model is replaced with the current in the poloidal coil of the stellarator.Plasma startup can be successfully achieved at injection powers in the hundreds of kilowatts range,resulting in electron densities on the order of 10^(17)-10^(18)m^(-3).展开更多
In this paper a new recursive method for ARMA model estimation is given. Same as in [1], theorder's estimator is strongly consistent, and the parameter's estimators defer to CLT and LILunder a natural conditio...In this paper a new recursive method for ARMA model estimation is given. Same as in [1], theorder's estimator is strongly consistent, and the parameter's estimators defer to CLT and LILunder a natural condition. Compared with the previous metheds suggested by Hannan & Kavalieris(1984), Wang Shouren & Chen Zhaoguo (1985) and Franke (1985), this methed has some advantages:the amount of calculat on work is smaller, the minimum-phase property of coeffcient estimators canbe guaranteed,the BAN estimators for MA or AR model can be obtained directly,and the simulationshows that this method is more accurate in estimating the order and parameters.展开更多
The type B aortic dissection(TBAD)is a perilous disease with high morbidity and mortality rates.The hemodynamics of TBAD in different scenarios has been widely studied by computational fluid dynamics(CFD)research.Howe...The type B aortic dissection(TBAD)is a perilous disease with high morbidity and mortality rates.The hemodynamics of TBAD in different scenarios has been widely studied by computational fluid dynamics(CFD)research.However,the flow pattern and wave propagation characteristics in the cardiovascular system with TBAD are not yet clear,and the effect of the distal tear is still unknown.In this work,a onedimensional(1D)cardiovascular system model coupling with a zero-dimensional(0D)lumped-parameter model is introduced to study the hemodynamics and wave propagation in the cardiovascular system.The results show that the proposed 0D-1D method well captures the oscillation and retrograde characteristics for the flow in the false lumen(FL),and the smaller distal tear damps the retrograde flow.Besides,the distal tear should also be paid attention to,and the wave intensity(WI)can be used as an access mark of the degree of the aortic dissection(AD).展开更多
We present the physics program of the open charm and heavy cc states above the DD production energy threshold, which will be studied with the BES-Ⅲ detector at the BEPC- Ⅱ collider in the coming years. Based on some...We present the physics program of the open charm and heavy cc states above the DD production energy threshold, which will be studied with the BES-Ⅲ detector at the BEPC- Ⅱ collider in the coming years. Based on some full Monte Carlo simulations with the BES-Ⅲ detector, we predict the accuracy levels on measuring some physical quantities related to D^0, D^＋ and Ds^＋ decays as well as some non-charmed decays of the heavy cc states.展开更多
文摘Diesel engines have proven over the years important in terms of efficiency and fuel consumption to power generation ratio. Many research works show the potential of biodiesel as a substitute for conventional gasoil. Mainly, previous and recent researches have focused on experimental investigation of diesel engine performance fuelled by biodiesel. Researches on the mathematical description of diesel engine process running on biodiesel are scarce, and mostly about chemical and thermodynamic description of the combustion process of biodiesel rather than performance studies. This work describes a numerical investigation on the performance analysis of a diesel engine fuelled by palm oil biodiesel. The numerical investigation was made using a semi empirical 0D model based on Wiebe’s and Watson’s model which was implemented via the open access numerical calculation software Scilab. The model was validated first by comparing with experimental pressure and performance data of a one cylinder engine at rated speed and secondly by comparing with a six cylinders engine performance data at various crankshaft rotational speeds. Simulations were then made to analyze the engine performance when running on biodiesel. The calculations were made at constant combustion duration and constant coefficient of excess air. Results showed that the model matches the overall experimental data, such as the power output and peak cylinder pressure. The ignition delay was somehow underestimated by the model for the first experiment, which caused a slight gap on in cylinder pressure curve, whereas it predicted the average ignition delay fairly well for the second set of validation. The simulations of engine performance when running on biodiesel confirmed results obtained in previous experimental researches on biodiesel. The model will be further investigated for engine control when shifting to biodiesel fuel.
基金funded by CEA,EDF and Framatomefinancial and scientific support of CEA Cadarache.
文摘Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.
文摘A 30 mm electrothermal-chemical(ETC) gun experimental system is employed to research the burning rate characteristics of 4/7 high-nitrogen solid propellant. Enhanced gas generation rates(EGGR) of propellants during and after electrical discharges are verified in the experiments. A modified 0D internal ballistic model is established to simulate the ETC launch. According to the measured pressure and electrical parameters, a transient burning rate law including the influence of EGGR coefficient by electric power and pressure gradient(dp/dt) is added into the model. The EGGR coefficient of 4/7 high-nitrogen solid propellant is equal to 0.005 MW-1. Both simulated breech pressure and projectile muzzle velocity accord with the experimental results well. Compared with Woodley's modified burning rate law, the breech pressure curves acquired by the transient burning rate law are more consistent with test results. Based on the parameters calculated in the model, the relationship among propellant burning rate, pressure gradient(dp/dt) and electric power is analyzed. Depending on the transient burning rate law and experimental data, the burning of solid propellant under the condition of plasma is described more accurately.
基金supported by the National Key Research and Development Program of China(Nos.2022YFE03070000 and 2022YFE03070003)National Natural Science Foundation of China(Nos.12375220 and 12075114)。
文摘According to the physics of tokamak start-up,this study constructs a zero-dimensional(0D)model applicable to electron cyclotron(EC)wave assisted start-up in NCST spherical torus(spherical tokamak)and CN-H1 stellarators.Using the constructed 0D model,the results obtained in this study under the same conditions are compared and validated against reference results for pure hydrogen plasma start-up in tokamak.The results are in good agreement,especially regarding electron temperature,ion temperature and plasma current.In the presence of finite Ohmic electric field in the spherical tokamak,a study on the EC wave assisted start-up of the NCST plasma at frequency of 28 GHz is conducted.The impact of the vertical magnetic field B_(v)on EC wave assisted start-up,the relationship between EC wave injection power P_(inj),Ohmic electric field E,and initial hydrogen atom density n_(H0)are explored separately.It is found that under conditions of Ohmic electric field lower than ITER(~0.3 V m^(-1)),EC wave can expand the operational space to achieve better plasma parameters.Simulating the process of28 GHz EC wave start-up in the CN-H1 stellarator plasma,the plasma current in the zerodimensional model is replaced with the current in the poloidal coil of the stellarator.Plasma startup can be successfully achieved at injection powers in the hundreds of kilowatts range,resulting in electron densities on the order of 10^(17)-10^(18)m^(-3).
文摘In this paper a new recursive method for ARMA model estimation is given. Same as in [1], theorder's estimator is strongly consistent, and the parameter's estimators defer to CLT and LILunder a natural condition. Compared with the previous metheds suggested by Hannan & Kavalieris(1984), Wang Shouren & Chen Zhaoguo (1985) and Franke (1985), this methed has some advantages:the amount of calculat on work is smaller, the minimum-phase property of coeffcient estimators canbe guaranteed,the BAN estimators for MA or AR model can be obtained directly,and the simulationshows that this method is more accurate in estimating the order and parameters.
基金supported by the National Natural Science Foundation of China(No.51976026)the Fundamental Research Funds of Central Universities of China(Nos.DUT22YG206 and DUT21JC25)。
文摘The type B aortic dissection(TBAD)is a perilous disease with high morbidity and mortality rates.The hemodynamics of TBAD in different scenarios has been widely studied by computational fluid dynamics(CFD)research.However,the flow pattern and wave propagation characteristics in the cardiovascular system with TBAD are not yet clear,and the effect of the distal tear is still unknown.In this work,a onedimensional(1D)cardiovascular system model coupling with a zero-dimensional(0D)lumped-parameter model is introduced to study the hemodynamics and wave propagation in the cardiovascular system.The results show that the proposed 0D-1D method well captures the oscillation and retrograde characteristics for the flow in the false lumen(FL),and the smaller distal tear damps the retrograde flow.Besides,the distal tear should also be paid attention to,and the wave intensity(WI)can be used as an access mark of the degree of the aortic dissection(AD).
基金Supported by National Natural Science Foundation of China (10935007)National Key Basic Research Program (973 by MOST) (2009CB825200XX)
文摘We present the physics program of the open charm and heavy cc states above the DD production energy threshold, which will be studied with the BES-Ⅲ detector at the BEPC- Ⅱ collider in the coming years. Based on some full Monte Carlo simulations with the BES-Ⅲ detector, we predict the accuracy levels on measuring some physical quantities related to D^0, D^＋ and Ds^＋ decays as well as some non-charmed decays of the heavy cc states.
