A CombustionModel for Explosive Charge Affected by a Bottom Gap in the Launch Environment

Launch safety of explosive charges has become an urgent problem to be solved by all countries in the world aslaunch situation of ammunition becomes consistentlyworse.However, the existing numericalmodels have differentdefects. This paper formulates an efficient computational model of the combustion of an explosive charge affectedby a bottom gap in the launch environment in the context of the material point method. The current temperatureis computed accurately from the heat balance equation, and different physical states of the explosive charges areconsidered through various equations of state. Microcracks in the explosive charges are described with respectto the viscoelastic statistical crackmechanics (Visco–SCRAM) model. Themethod for calculating the temperatureat the bottomof the explosive charge with respect to the bottomgap is described. Based on this combustionmodel,the temperature history of a Composition B (COMB) explosive charge in the presence of a bottom gap is obtainedduring the launch process of a 155-mm artillery. The simulation results show that the bottom gap thickness shouldbe no greater than 0.039 cm to ensure the safety of the COM B explosive charge in the launch environment. Thisconclusion is consistent with previous results and verifies the correctness of the proposed model. Ultimately, thispaper derives amathematical expression for themaximumtemperature of the COMB explosive chargewith respectto the bottomgap thickness (over the range of 0.00–0.039 cm), and establishes a quantitative evaluationmethod forthe launch safety of explosive charges.The research results provide some guidance for the assessment and detectionof explosive charge safety in complex launch environments.

Combustion crack-network reaction evolution model for highly-confined explosives

The evolution behavior of combustion crack reaction of highly confined solid explosives after non-shock ignition is governed by multiple dynamic processes,including intrinsic combustion of explosives,crack propagation,and rapid growth of combustion surface area.Here,the pressure increase can accelerate the combustion rate of explosives,and the crack propagation can enlarge the combustion surface area.The coupling between these two effects leads to the self-enhanced combustion of explosive charge system,which is the key mechanism for the reaction development after ignition.In this study,combustion cracknetwork(CCN) model is established to describe the evolution of combustion crack reaction of highly confined solid explosives after non-shock ignition and quantify the reaction violence.The feasibility of the model is verified by comparing the computational and experimental results.The results reveal that an increase in charge structure size causes an increase in the time of crack pressurization and extension of cracks due to the high temperature-generated gas flow and surface combustion during the initial stage of explosive reaction,but when the casing is fractured,the larger the charge structure,the more violent the late reaction and the larger the charge reaction degree.The input pressure has no obvious influence on the final reaction violence.Further,a larger venting hole area leads to better pressure relief effect,which causes slower pressure growth inside casing.Larger reserved ullage volume causes longer lowpressure induction stage,which further restrains the internal pressure growth.Furthermore,the stronger the casing constraint,the more rapid the self-enhanced combustion of the high temperaturegenerated gas,which results in more violent charge reaction and larger charge reaction degree during casing break.Overall,the proposed model can clarify the effects of intrinsic combustion rate of explosives,charge structure size,input pressure,relief area,ullage volume,and constraint strength on the reaction evolution,which can provide theoretical basis for violence evaluation and safety design for ammunition under accident stimulus.

Kinetic models of natural gas combustion in an internal combustion engine

In this study, combustion of methane was simulated using four kinetic models of methane in CHEMKIN 4.1.1 for 0-D closed internal combustion （IC） engine reactor. Two detailed （GRIMECH3.0 ＆ UBC MECH2.0） and two reduced （One step ＆ Four steps） models were examined for various IC engine designs. The detailed models （GRIMECH3.0, ＆ UBC MECH2.0） and 4-step models successfully predicted the combustion while global model was unable to predict any combustion reaction. This study illustrated that the detailed model showed good concordances in the prediction of chamber pressure, temperature and major combustion species profiles. The detailed models also exhibited the capabilities to predict the pollutants formation in an IC engine while the reduced schemes showed failure in the prediction of pollutants emissions. Although, there are discrepancies among the profiles of four considered model, the detailed models （GRIMECH3.0 ＆ UBC MECH2.0） produced the acceptable agreement in the species prediction and formation of pollutants.

Estimating the Operation Status of Steam Cracking Furnace Using Numerical Simulation with Combustion Models

An accurate and complete geometric model was constructed to simulate the combustion, flow and temperature environment in the radiant section of the steam cracking furnace. Simulation of flow and radiation status has utilized the standard k-ε model and P1 model. The finite-rate/eddy-dissipation (finite-rate/ED) combustion model and non-premixed combustion model were both used to simulate accurately the combustion and the operation status of the steam cracking furnace. Three different surfaces of the steam cracking furnace were obtained from the simulation, namely:the flue gas temperature field of the entrance surface in long flame burners, the central surface location of tubes, and the crossover section surface. Detailed information on the flue gas temperature and the mass concentration fraction of these different surfaces in the steam cracking furnace can also be obtained by the simulation. This paper analyzed and compared the simulation results with the two combustion models, estimated the operation status of the steam cracking furnace, and reported that the finite-rate/ED model is appropriate to simulate the steam cracking furnace by comparing key simulation data with actual test data. This work has also provided a theoretical basis for simulating and operating the steam cracking furnace.

Three-Dimensional Models for Analyzing the Cyclic Variations in a Lean Burn CNG Engine

Three-dimensional models, consisting of the flame kernel formation model, flame kernel development model and natural gas single step reaction model, are used to analyze the contribution of cyclic equivalence ratio variations to cyclic variations in the compressed natural gas （CNG） lean burn spark ignition engine. Computational results including the contributions of equivalence ratio cyclic variations to each combustion stage and effects of engine speed to the extent of combustion variations are discussed. It is concluded that the equivalence ratio variations affect mostly the main stage of combustion and hardly influence initial kernel development stage.

DYNAMICS MODEL AND SIMULATION OF FLAT VALVE SYSTEM OF INTERNAL COMBUSTION WATER PUMP

The dynamics differential equations are constructed, and the initial conditions are also given. Simulation shows the following conclusions： The water pressure in cylinder has great instantaneous pulsation and phase step when outlet valve or inlet valve opens, but is more gently in other time; The volume efficiency is influenced by the output pressure slightly, and decreases as the working rotational speed increases; When the inherent frequency of the valves is integer multiple of the working frequency, the volume efficiency of system will decrease evidently.

Three-Dimensional Tidal Model and Its Application to Numerical Simulation of Water Quality in Coastal Waters

The turbulence mechanism plays an important part in the mixing process and momentum transfer of turbulence. A three-dimensional Prandtl mixing length tidal model has been developed to simulate tidal flows and water quality. The eddy viscosities and diffusivities are computed from the Prandtl mixing length model. In order to model the water quality of an estuary or coastal area many interdependent processes need to be simulated. These may be conveniently separated into three main groups: transport and mixing processes, biochemical interaction of water quality variables and the utilization and re-cycling of nutrients by living matter. The model simulates full oxygen and nutrient balance, primary productivity and the transport, reaction mechanism and fate of pollutants over tidal time-scales. The model is applied to numerical simulation of tidal flows and water quality in Dalian Bay. The model has been calibrated against a limited data set of historical water quality observations and in general demonstrates excellent agreement with all available data.

Mathematical Model of Combustion in Blunt Annular Ceramic Burner

The computer simulation of the combustion process in blast furnace(BF) stove has been studied by using the k-ε-g turbulent diffusion flame model.The combustion process in blunt annular ceramic burner was calculated by using the software.The profiles of gas and air velocity,temperature of the combustion products,concentration of the components,and the shape and length of the flame during combustion have been researched.Compared with the original annular ceramic burner,the new design of the blunt one improves the mixing of the gas and the air significantly,and shortened the length of the flame.

The three-dimension model for the rock-breaking mechanism of disc cutter andanalysis of rock-breaking forces

To study the rock deformation with three- dimensional model under rolling forces of disc cutter, by car- rying out the circular-grooving test with disc cutter rolling around on the rock, the rock mechanical behavior under rolling disc cutter is studied, the mechanical model of disc cutter rolling around the groove is established, and the the- ory of single-point and double-angle variables is proposed. Based on this theory, the physics equations and geometric equations of rock mechanical behavior under disc cutters of tunnel boring machine （TBM） are studied, and then the bal- ance equations of interactive forces between disc cutter and rock are established. Accordingly, formulas about normal force, rolling force and side force of a disc cutter are de- rived, and their validity is studied by tests. Therefore, a new method and theory is proposed to study rock- breaking mech- anism of disc cutters.

Large eddy simulation of turbulent premixed and stratified combustion using flame surface density model coupled with tabulation method

Large eddy simulations(LESs) are performed to investigate the Cambridge premixed and stratified flames, SwB1 and SwB5, respectively. The flame surface density(FSD) model incorporated with two different wrinkling factor models, i.e., the Muppala and Charlette2 wrinkling factor models, is used to describe combustion/turbulence interaction, and the flamelet generated manifolds(FGM) method is employed to determine major scalars. This coupled sub-grid scale(SGS) combustion model is named as the FSD-FGM model. The FGM method can provide the detailed species in the flame which cannot be obtained from the origin FSD model. The LES results show that the FSD-FGM model has the ability of describing flame propagation, especially for stratified flames. The Charlette2 wrinkling factor model performs better than the Muppala wrinkling factor model in predicting the flame surface area change by the turbulence.The combustion characteristics are analyzed in detail by the flame index and probability distributions of the equivalence ratio and the orientation angle, which confirms that for the investigated stratified flame, the dominant combustion modes in the upstream and downstream regions are the premixed mode and the back-supported mode, respectively.

Numerical Modeling of NO Formation during Packed-bed Combustion of Coke Granules

A comprehensive kinetic model of NO formation during coke combustion in packed-bed in presence of noncombustible particles was developed. The detailed homogeneous gas-phase chemistry (including 102 chemical reactions), heterogeneous gas-solid chemistry (including 11 reactions) of coke combustion and NO formation, and the heat and mass transfer were taken into account in the present model. The governing equations which are strongly coupled, non-linear and unsteady with 26 unknowns in total, were dispersed into differential equations with the finite differential method. Meanwhile, all the differential equations were numerically solved to give the time-histories and space-distributions oftemperatures of the bed and gas phase as well as the concentrations of all the gaseous species. By comparison, the experimental data were explained well by the calculated results. Based on the kinetic and mathematical model, the effects of O2 content of inlet gas, the initial chemical analysis of coke, bed-temperature and local reductive atmosphere (CO/O2) on NO formation during packed-bed coke combustion were numerically discussed. It was found that coke samples with a higher initial content of [N] and volatile matters, combusted under a suitable O2-containing atmosphere produced less NO emission. The reactions between CO and NO, catalyzed by high temperature surface of coke particles may be responsible for efficient reduction of NO.

THE GENERALIZED RIEMANN PROBLEM FOR A SCALAR NONCONVEX COMBUSTION MODEL-THE PERTURBATION ON INITIAL BINDING ENERGY

In this article, we study the generalized Riemann problem for a scalar non- convex Chapman-Jouguet combustion model in a neighborhood of the origin （t 〉 0） on the （x, t） plane. We focus our attention to the perturbation on initial binding energy. The solutions are obtained constructively under the entropy conditions. It can be found that the solutions are essentially different from the corresponding Riemann solutions for some cases. Especially, two important phenomena are observed： the transition from detonation to deflagration followed by a shock, which appears in the numerical simulations [7, 27]; the transition from deflagration to detonation （DDT）, which is one of the core problems in gas dynamic combustion.

A Simplified Model for SO_(2) Generation during Spontaneous Combustion of Coal Gangue

A simplified model for SO_(2) generation during spontaneous combustion of coal gangue was put forward and validated using the measured data.Using the proposed model,the effects of initial temperature inside the gangue and fresh air supply on SO_(2) generation were discussed.The results showed that,higher initial temperature inside the gangue could accelerate the oxidation rate of FeS_(2) and increase the maximum concentration of SO_(2).If initial temperature inside the gangue increased by about 37%,the total SO_(2) generation increased by 166%.Fresh air supply had less significant effect on the oxidation rate of FeS_(2).However,the higher the fresh air supply was,the more FeS_(2) could be oxidized,which ultimately produced more SO_(2).Although the computed results and the measured data concerning the inner locations inside the gangue had a certain degree of error,the proposed model can provide a relatively precise total release of SO_(2) within acceptable accuracy.Besides,this method provides a useful prototype to predict the generation of hazardous materials,such as CO,NO_(x),and chlorine during the spontaneous combustion of coal gangue.

Response Surface Modeling of Fuel Rich and Fuel Lean Catalytic Combustion of the Stabilized Confined Turbulent Gaseous Diffusion Flames

The Response Surface Methodology (RSM) has been applied to explore the thermal structure of the experimentally studied catalytic combustion of stabilized confined turbulent gaseous diffusion flames. The Pt/γAl2O3 and Pd/γAl2O3 disc burners were situated in the combustion domain and the experiments were performed under both fuel-rich and fuel-lean conditions at a modified equivalence (fuel/air) ratio (ø) of 0.75 and 0.25 respectively. The thermal structure of these catalytic flames developed over the Pt and Pd disc burners were inspected via measuring the mean temperature profiles in the radial direction at different discrete axial locations along the flames. The RSM considers the effect of the two operating parameters explicitly (r), the radial distance from the center line of the flame, and (x), axial distance along the flame over the disc, on the measured temperature of the flames and finds the predicted maximum temperature and the corresponding process variables. Also the RSM has been employed to elucidate such effects in the three and two dimensions and displays the location of the predicted maximum temperature.

Kinetic analysis and modeling of maize straw hydrochar combustion using a multi-Gaussian-distributed activation energy model

Combustion kinetics of the hydrochar was investigated using a multi-Gaussian-distributed activation energy model(DAEM)to ex-pand the knowledge on the combustion mechanisms.The results demonstrated that the kinetic parameters calculated by the multi-Gaussian-DAEM accurately represented the experimental conversion rate curves.Overall,the feedstock combustion could be divided into four stages:the decomposition of hemicellulose,cellulose,lignin,and char combustion.The hydrochar combustion could in turn be divided into three stages:the combustion of cellulose,lignin,and char.The mean activation energy ranges obtained for the cellulose,lignin,and char were 273.7-292.8,315.1-334.5,and 354.4-370 kJ/mol,respectively,with the standard deviations of 2.1-23.1,9.5-27.4,and 12.1-22.9 kJ/mol,re-spectively.The cellulose and lignin contents first increased and then decreased with increasing hydrothermal carbonization(HTC)temperature,while the mass fraction of char gradually increased.

Modeling and Simulation of Coke Combustion Regeneration for Coked Cr2O3/Al2O3 Propane Dehydrogenation Catalyst

A heterogeneous model is developed for the regeneration of the Cr2O3/Al2O3 catalyst for the propane dehydrogenation process by considering the internal mass transfer and external mass/heat transfer during the coke combustion.Simulation shows that under practical operating conditions,multi-steady states exist for the catalyst pellets and the catalyst temperature is sensitive to gas temperature.However,at increased mass flow rate or lowered oxygen concentration,multi-steady states will not appear.Under the strong influences of film diffusion,the coke in the packed bed reactor will first be exhausted at the inlet,while if the film diffusion resistance is decreased,the position of first coke exhaustion moves toward the outlet of the reactor.

THE WEAK DETONATION SOLUTIONS TO THE SIMPLEST COMBUSTION MODEL

In the present paper the authors prove that all the generalised entropy,solutions of the CJ-model, which for the given Riemann initial data are a oneparameter family u(eta) characterised by weak detonation discontinuity point eta, are just the limits of the admissible solutions u(ek) Of the selfsimilar combustion model. In fact, the authors prove that for any possible eta, there exists a constant B>0 s.t.

Combustion modeling in a model combustor

The flow-field of a propane-air diffusion flame combustor with interior and exterior conjugate heat transfers was numerically studied.Results obtained from four combustion models,combined with the re-normalization group(RNG) k-ε turbulence model,discrete ordinates radiation model and enhanced wall treatment are presented and discussed.The results are compared with a comprehensive database obtained from a series of experimental measurements.The flow patterns and the recirculation zone length in the combustion chamber are accurately predicted,and the mean axial velocities are in fairly good agreement with the experimental data,particularly at downstream sections for all four combustion models.The mean temperature profiles are captured fairly well by the eddy dissipation(EDS),probability density function(PDF),and laminar flamelet combustion models.However,the EDS-finite-rate combustion model fails to provide an acceptable temperature field.In general,the flamelet model illustrates little superiority over the PDF model,and to some extent the PDF model shows better performance than the EDS model.

Modelling of a tubular membrane contactor for pre-combustion CO_2 capture using ionic liquids:Influence of the membrane configuration, absorbent properties and operation parameters

A membrane contactor using ionic liquids(ILs) as solvent for pre-combustion capture CO_2 at elevated temperature(303-393 K) and pressure(20 bar) has been studied using mathematic model in the present work. A comprehensive two-dimensional(2 D) mass-transfer model was developed based on finite element method. The effects of liquid properties, membrane configurations, as well as operation parameters on the CO_2 removal efficiency were systematically studied. The simulation results show that CO_2 can be effectively removed in this process. In addition, it is found that the liquid phase mass transfer dominated the overall mass transfer. Membranes with high porosity and small thickness could apparently reduce the membrane resistance and thus increase the separation efficiency. On the other hand, the membrane diameter and membrane length have a relatively small influence on separation performance within the operation range.

THE ASYMPTOTIC LIMIT FOR A COMBUSTION MODEL IN REGARD TO INFINITE REACTION RATE

The Zeldovich-von Neumann-Doring model and the Chapman-Jouguet model for a simplified combustion model-Majda＇s model is studied.The author proves a uniform maximum norm estimate,then proves that as the rate of chemical reaction tends to infinity the solutions to the Zeldovich-von Neumann-Doring model tend to that of the Chapman-Jouguet model.The type of combustion waves is studied.This result is compared with the result of the projection and finite difference method for the same model.