Numerical Simulation of Deflagration to Detonation Transition in a Straight Duct: Effects of Energy Release and Detonation Stability

Numerical simulation based on the Euler equation and one-step reaction model is carried out to investigate the process of deflagration to detonation transition (DDT) occurring in a straight duct. The numerical method used includes a high resolution fifth-order weighted essentially non-oscillatory (WENO) scheme for spatial discretization, coupled with a third order total variation diminishing Runge-Kutta time stepping method. In particular, effect of energy release on the DDT process is studied. The model parameters used are the heat release at $q=50, 30, 25, 20, 15, 10$ and $5$, the specific heat ratio at $1.2$, and the activation temperature at $Ti=15$, respectively. For all the cases, the initial energy in the spark is about the same compared to the detonation energy at the Chapman-Jouguet (CJ) state. It is found from the simulation that the DDT occurrence strongly depends on the magnitude of the energy release. The run-up distance of DDT occurrence decreases with the increase of the energy release for $q$=50~20, and increases with the increase of the energy release for $q$=20~5. This phenomenon is found to be in agreement with the analysis of mathematical stability theory. It is suggested that the factors to strengthen the DDT would make the detonation more stable, and vice versa. Finally, it is concluded from the simulations that the interaction of the shock wave and the flame front is the main reason for leading to DDT.     

Nuclear hyperfine coupling interactions in the rotational spectra of the linear and bent isomers of HF-N2O

The rotational spectra of six isotopomers of the linear and bent isomers of HF-N₂O have been collected in the 7-18 GHz region with a Fourier transform microwave spectrometer. The nuclear hyperfine structure in the spectra produced by HF spin-spin coupling interaction and nuclear quadrupole coupling interactions due to the D nucleus of DF and the nuclei of N₂O have been resolved and analyzed. In the linear isomer, H in HF is bonded to the terminal N in N₂O. The NF bond lengths are 2.9808(2) Å for the HF-containing isotopomers and 2.9732(2) Å for the DF-containing isotopomers. The zero point angles are 23.1° for HF and 31-34° for N₂O. The hyperfine constants suggest that the HF bond is lengthened by 0.0105 Å upon complexation and that the electric field gradients of the two nitrogen nuclei in N₂O are perturbed differently in the complex. In the bent isomer, the hydrogen bond is formed between HF and O in N₂O. The intermolecular distances are 3.4942(2) Å for the HF-containing isotopomers and 3.4436(2) Å for the DF-containing isotopomers, with HF and N₂O forming angles of 34° and 46°, respectively, with the intermolecular axis. The nuclear quadrupole coupling constants of the two nitrogen nuclei do not indicate electric field gradient perturbation in this isomer.     

Many particle magnetic dipole-dipole and hydrodynamic interactions in magnetizable stent assisted magnetic drug targeting

The implant assisted magnetic targeted drug delivery system of Avilés, Ebner and Ritter is considered both experimentally (in vitro) and theoretically. The results of a 2D mathematical model are compared with 3D experimental results for a magnetizable wire stent. In this experiment a ferromagnetic, coiled wire stent is implanted to aid collection of particles which consist of single domain magnetic nanoparticles (radius ). In order to model the agglomeration of particles known to occur in this system, the magnetic dipole-dipole and hydrodynamic interactions for multiple particles are included. Simulations based on this mathematical model were performed using open source C++ code. Different initial positions are considered and the system performance is assessed in terms of collection efficiency. The results of this model show closer agreement with the measured in vitro experimental results and with the literature. The implications in nanotechnology and nanomedicine are based on the prediction of the particle efficiency, in conjunction with the magnetizable stent, for targeted drug delivery.     

One-dimensional optical bright screening photovoltaic photorefractive solitons, soliton pairs and incoherent interaction between them in BaTiO3 crystal

By the use of theory of optical solitons, nonlinear equations governed by screening photorefractive photovoltaic solitons are introduced. By means of numerical methods, intensity profile, refractive index change of the media and also stability under propagation are investigated. Self-bending of solitons under first-order diffusion is also studied. Besides, incoherent coupled bright-bright soliton pairs and incoherent interaction are studied.     

Numerical studies on autoignition and detonation development from a hot spot in hydrogen/air mixtures

Detonation development inside spark ignition engines can result in the so called super-knock with extremely high pressure oscillation above 200?atm. In this study, numerical simulations of autoignitive reaction front propagation in hydrogen/air mixtures are conducted and the detonation development regime is investigated. A hot spot with linear temperature distribution is used to induce autoignitive reaction front propagation. With the change of temperature gradient or hot spot size, three typical autoignition reaction front modes are identified: supersonic reaction front; detonation development and subsonic reaction front. The effects of initial pressure, initial temperature, fuel type and equivalence ratio on detonation development regime are examined. It is found that the detonation development regime strongly depends on mixture composition (fuel and equivalence ratio) and thermal conditions (initial pressure and temperature). Therefore, to achieve the quantitative prediction of super-knock in engines, we need use the detonation development regime for specific fuel at specific initial temperature, initial pressure, and equivalence ratio.     

Water behavior in the neighborhood of hydrophilic and hydrophobic membranes: Lessons from molecular dynamics simulations

The study of properties of water in the vicinity of surfaces poses a fascinating challenge. In this article we studied the behavior of water molecules in the neighborhood of membranes. We addressed the question of how these water molecules are influenced by the membranes’ hydrophilicity. Three systems were studied through molecular dynamics simulations: water in the presence of a hydrophilic membrane (PL), water in the presence of a hydrophobic (PB) one and water in the absence of membranes (BULK). Additionally, in order to study the dependence of the effect of the membrane on the behavior of neighboring water molecules with temperature, each system was simulated at three different temperatures (K): 250, 300 and 350. For each condition, kinetic and structural features were studied. The first feature involved the calculation of diffusion coefficients and activation energy. The second feature was evaluated through the study of water density and hydrogen bond distribution. From the present study we concluded that: (1) density studies underestimate the influence of both hydrophilic and hydrophobic membranes on the neighboring water molecules; (2) the hydrophilic and hydrophobic membranes disturb the hydrogen bond network within distances ranging from 1 to 8 nm, depending on the nature of the membrane and the temperature conditions; (3) the presence of a hydrophobic surface results in an enhancement of the natural hydrogen bond network present in liquid water, to a greater extent than what even an ordered Ih ice structure is able to achieve (i.e. PL membrane); (4) the structural enhancement due to the presence of a hydrophobic surface involves roughly 18 to 24 water hydration layers, for ambient and above temperature conditions.     

Molecular dynamics study on water self-diffusion in aqueous mixtures of methanol,ethylene glycol and glycerol: investigations from the point of view of hydrogen bonding

Molecular dynamics simulations have been performed to investigate the aqueous binary mixtures of alcohols, including methanol, ethylene glycol (EG) and glycerol of molalities ranging from 1 to 5 m at the temperatures of 273, 288 and 298 K, respectively. The primary purpose of this paper is to investigate the mechanism of water self-diffusion in water-alcohol mixtures from the point of view of hydrogen bonding. The effects of temperature and concentration on water self-diffusion coefficient are evaluated quantitatively in this work. Temperature and concentration to some extent affect the hydrogen bonding statistics and dynamics of the binary mixtures. It is shown that the self-diffusion coefficient of water molecules decreases as the concentration increases or the temperature decreases. Moreover, calculations of mean square displacements of water molecules initially with different number n of H-bonds indicate that the water self-diffusion coefficient decreases as n increases. We also studied the aggregation of alcohol molecules by the hydrophobic alkyl groups. The largest cluster size of the alkyl groups clearly increases as the concentration increases, implying the emergence of a closely connected network of water and alcohols. The clusters of water and alcohol that interacted could block the movement of water molecules in binary mixtures. These findings provide insight into the mechanisms of water self-diffusion in aqueous binary mixtures of methanol, EG and glycerol.