Investigation of Vaporized Kerosene Injection and Combustion in a Supersonic Model Combustor

Injection and combustion of vaporized kerosene was experimentally investigated in a Mach 2.5 model combustor at various fuel temperatures and injection pressures. A unique kerosene heating and delivery system, which can prepare heated kerosene up to 820 K at a pressure of 5.5 MPa with negligible fuel coking, was developed. A three-species surrogate was employed to simulate the thermophysical properties of kerosene. The calculated thermophysical properties of surrogate provided insight into the fuel flow control in experiments. Kerosene jet structures at various preheat temperatures injecting into both quiescent environment and a Mach 2.5 crossflow were characterized. It was shown that the use ofvaporized kerosene injection holds the potential of enhancing fuel-air mixing and promoting overall burning. Supersonic combustion tests further confirmed the preceding conjecture by comparing the combustor performances of supercritical kerosene with those of liquid kerosene and effervescent atomization with hydrogen barbotage. Under the similar flow conditions and overall kerosene equivalence ratios, experimental results illustrated that the combustion efficiency of supercritical kerosene increased approximately 10-15% over that of liquid kerosene, which was comparable to that of effervescent atomization.

QM/MM and Classical Molecular Dynamics Simulation of Histidine-Tagged Peptide Immobilization on Nickel Surface

The hybrid quantum mechanics (QM) and molecular mechanics (MM) method is employed to simulate the His-tagged peptide adsorption to ionized region of nickel surface. Based on the previous experiments, the peptide interaction with one Ni ion is considered. In the QM/MM calculation, the imidazoles on the side chain of the peptide and the metal ion with several neighboring water molecules are treated as QM part calculated by “GAMESS”, and the rest atoms are treated as MM part calculated by “TINKER”. The integrated molecular orbital/molecular mechanics (IMOMM) method is used to deal with theQMpart with the transitional metal. By using the QM/MM method, we optimize the structure of the synthetic peptide chelating with a Ni ion. Different chelate structures are considered. The geometry parameters of the QM subsystem we obtained by QM/MM calculation are consistent with the available experimental results. We also perform a classical molecular dynamics (MD) simulation with the experimental parameters for the synthetic peptide adsorption on a neutral Ni(1 0 0) surface. We find that half of the His-tags are almost parallel with the substrate, which enhance the binding strength. Peeling of the peptide from the Ni substrate is simulated in the aqueous solvent and in vacuum, respectively. The critical peeling forces in the two environments are obtained. The results show that the imidazole rings are attached to the substrate more tightly than other bases in this peptide.

Variation of drag coefficient with wind wave development

Turbulent air flows over developing wind waves in the air-sea boundary layer are numerically simulated without considering wave breaking. Influences of wind waves on air flows are considered using a model of significant wave and surface roughness, with a formula proposed for calculating the surface roughness, k - epsilon model is adopted to simulate turbulent flows. The results of the drag coefficient and turbulence characteristics agree well with the observations.

Numerical Simulation of Pulsed Liquid Fluidized Bed and Its Experimental Validation

The pulsed liquid fluidized bed was studied using numerical simulation and experimental methods, The area-averaged two-fluid model (TFM) was used to simulate the pulsed fluidization. The bed expansion and collapse processes were simulated first and the phenomena obtained from the calculation were consistent with our previous experiments and observations. In the pulsed fluidization, the variation of bed height, the variations of particle velocity and concentration distribution were obtained and analyzed. Experiments were carried out to validate the simulation results. The pressure variation with time at different locations was measured using pressure transducers and compared with the simulated results. The variations of bed height and particle concentration distribution were recorded using a digital video camera recorder. The results were consistent with the simulation results as a whole.

Molecular Dynamics Simulation of Barnacle Cement

Barnacle cement is an underwater adhesive that is used for permanent settlement. Its main components are insoluble protein complexes that have not been fully studied. In present article, we chose two proteins of barnacle cement for study, 36-KD protein and Mrcp-100K protein. In order to investigate the characteristic of above two proteins, we introduced the method of molecular modeling. And the simulation package GROMACS was used to simulate the behavior of these proteins. In this article, before the simulations, we introduce some theories to predict the time scale for polymer relaxation. During the simulation, we mainly focus on two properties of these two proteins: structural stability and adhesive force to substrate. First, we simulate the structural stability of two proteins in water, and then the stability of 36-KD protein in seawater environment is investigated.We find that the stability varies in the different environments. Next, to study adhesive ability of two proteins, we simulate the process of peeling the two proteins from the substrate (graphite). Then, we analyze the main reasons of these results. We find that hydrogen bonds in proteins play an important role in the protein stability. In the process of the peeling, we use Lennard–Jones 12-6 potential to calculate the van der Waals interactions between proteins and substrate.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Adaptive Mesh Refinement FEM for Damage Evolution of Heterogeneous Brittle Media

Damage evolution of heterogeneous brittle media involves a wide range of length scales. The coupling between these length scales underlies the mechanism of damage evolution and rupture. However, few of previous numerical algorithms consider the effects of the trans-scale coupling effectively. In this paper, an adaptive mesh refinement FEM algorithm is developed to simulate this trans-scale coupling. The adaptive serendipity element is implemented in this algorithm, and several special discontinuous base functions are created to avoid the incompatible displacement between the elements. Both the benchmark and a typical numerical example under quasi-static loading are given to justify the effectiveness of this model. The numerical results reproduce a series of characteristics of damage and rupture in heterogeneous brittle media.

Deformation and Failure of Polymer Bonded Explosives under Diametric Compression Test

The tensile deformation and failure of polymer bonded explosives (PBXs), a particulate composite, is studied in this paper. Two HMX-based PBXs with different binder were selected for study. A diametric compression test, in which a disc-shaped specimen is loaded diametrically, was chosen to generate tensile failure in the materials. The quasi-static tensile properties and the tensile creep properties were studied by using conventional displacement transducers to measure the lateral strain along the horizontal diameter. The whole-field in-plane creep deformation was measured by using the technique of high resolution moire´ interferometry. Real time microscopic examination was conducted to monitor the process of deformation and failure of PBXs by using a scanning electron microscope equipped with a loading stage. A manifold method (MM) was used to simulate the deformation and failure of PBX samples under the diametric compression test, including the crack initiation, crack propagation and final cleavage fracture. The mechanisms of deformation and failure of PBXs under diametric compression were analyzed. The diametric compression test and the techniques developed in this research have proven to be applicable to the study of tensile properties of PBXs.

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Behaviour analysis of numerical solutions and simulation of coherent structures in compressible mixing layers

For understanding the correctness of simulations the behaviour of numerical solutions is analysed, Tn order to improve the accuracy of solutions three methods are presented. The method with GVC (group velocity control) is used to simulate coherent structures in compressible mixing layers. The effect of initial conditions for the mixing layer with convective Mach number 0.8 on coherent structures is discussed. For the given initial conditions two types of coherent structures in the mixing layer are obtained.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Studies of the Effect of a Coal Concentrator on NO Formation in Swirling Coal Combustion

To develop low-pollution burners, the effect of a coal concentrator on NO formation in swirling coal combustion is studied using both numerical simulation and experiments. The isothermal gas-particle two-phase velocities and particle concentration in a cold model of swirl burners with and without coal concentrators were measured using the phase Doppler particle anemometer (PDPA). A full two-fluid model of reacting gas-particle flows and coal combustion with an algebraic unified second-order moment (AUSM) turbulence-chemistry model for the turbulent reaction rate of NO formation are used to simulate swirling coal combustion and NO formation with different coal concentrators. The results give the turbulent kinetic energy, particle concentration, temperature and NO concentration in cases of with and without coal concentrators. The predicted results for cold two-phase flows are in good agreement with the PDPA measurement results, showing that the coal concentrator increases the turbulence and particle concentration in the recirculation zone. The combustion modeling results indicate that although the coal concentrator increases the turbulence and combustion temperature, but still can remarkably reduce the NO formation due to creating high coal concentration in the recirculation zone.

The investigation on dynamic fracture behavour of materials under compressive-shear combined stress waves

In this paper, an improved plate impact experimental technique is presented for studying dynamic fracture mechanism of materials, under the conditions that the impacting loading is provided by a single pulse and the loading time is in the sub-microsecond range. The impacting tests are carried out on the pressure-shear gas gun. The loading rate achieved is dK/dt similar to 10(8) MPa m(1/2) s(-1). With the elimination of influence of the specimen boundary, the plane strain state of a semi-infinite crack in an infinite elastic plate is used to simulate the deformation fields of crack tip. The single pulses are obtained by using the "momentum trap" technique. Therefore, the one-time actions of the single pulse are achieved by eradicating the stress waves reflected from the specimen boundary or diffracted from the crack surfaces. In the current study, some important phenomena have been observed. The special loading of the single pulse can bring about material damage around crack tip, and affect the material behavior, such as kinking and branching of the crack propagation. Failure mode transitions from mode I to mode II crack are observed under asymmetrical impact conditions. The mechanisms of the dynamic crack propagation are consistent with the damage failure model.

A Semi-Empirical Equation of Penetration Depth on Concrete Target Impacted by Ogive-Nose Projectiles

In this paper, the penetration process of ogive-nose projectiles into the semi-infinite concrete target is investigated by the dimensional analysis method and FEM simulation. With the dimensional analysis, main non-dimensional parameters which control the penetration depth are obtained with some reasonable hypothesis. Then, a new semi-empirical equation is present based on the original work of Forrestal et al., has only two non-dimensional combined variables with definite physical meanings. To verify this equation, prediction results are compared with experiments in a wide variation region of velocity. Then, a commercial FEM code, LS-DYNA, is used to simulate the complex penetration process, that also show the novel semi-empirical equation is reasonable for determining the penetration depth in a concrete target.

Modelling of Laminar Plasma Jet Impinging on a Flatplate with Approximate Box Relaxation Method

Approximate Box Relaxation method was used t'o simulate a plasma jet flow impinging on a flatplate at atmospheric pressure, to achieve a better understanding of the characteristics of plasma jet in materials surface treating. The flow fields under different conditions were simulated and analyzed. The distributions of temperature, velocity and pressure were obtained by modelling. Computed results indicate that this numerical method is suitable for simulation of the flow characteristics of plasma jet: and is helpful for understanding of the mechanism of the plasma-material processing.