Search

Molecular dynamics simulations have been performed to investigate the behavior of the EMIM- BF4 ionic liquid for conditions similar to those of electrospray thrusters. To study the physics of ion extrusion a large system composed of approximately 2160 ion pairs was placed inside a platinum capillary and equilibrated. Diffusion coefficients and electrical conductivity were calculated and compared with literature values, showing good agreement and validating the use of a constrained dynamics model with a coarse-grained potential. An electric field was applied to the system in the longitudinal direction and the fraction of solvated and non-solvated emitted ions was analyzed. The threshold electric field for particle emission was found to be 1.2 V/nm, consistent with other reported work. [ABSTRACT FROM AUTHOR]

Using molecular dynamics (MD), we have studied the mechanism of heat accommodation between carbon dioxide clusters and monomers for temperatures and cluster size conditions that exist in homogeneous condensing supersonic expansion plumes. The work was motivated by our meso-scale direct simulation Monte Carlo and Bhatnagar-Gross-Krook based condensation simulations where we found that the heat accommodation model plays a key role in the near-field of the nozzle expansion particularly as the degree of condensation increases [R. Kumar, Z. Li, and D. Levin, Phys. Fluids 23, 052001 (2011)]. The heat released by nucleation and condensation and the heat removed by cluster evaporation can be transferred or removed from either the kinetic or translational modes of the carbon dioxide monomers. The molecular dynamics results show that the time required for gas-cluster interactions to establish an equilibrium from an initial state of non-equilibrium is less than the time step used in meso-scale analyses [R. Kumar, Z. Li, and D. Levin, Phys. Fluids 23, 052001 (2011)]. Therefore, the good agreement obtained between the measured cluster and gas number density and gas temperature profiles with the meso-scale modeling using the second energy exchange mechanism is not fortuitous but is physically based. Our MD simulations also showed that a dynamic equilibrium is established by the gas-cluster interactions in which condensation and evaporation processes take place constantly to and from a cluster. [ABSTRACT FROM AUTHOR]

A coalescence model for homogeneous condensation of ammonia in supersonic expansions to vacuum has been developed using molecular dynamics trajectory calculations. The MD calculations show that the sticking probability increases as the ammonia cluster size increases or the cluster temperature decreases. In addition, the sensitivity of the sticking probability to cluster size decreases as the temperature decreases. Comparison of the Ashgriz-Poo semiempirical coalescence model with MD simulations show that for cluster sizes larger than 100 the former model may be used. To model homogeneous nucleation in an ammonia jet, direct simulation Monte Carlo (DSMC) simulations were performed for different stagnation pressure conditions using the MD simulation outcomes for smaller cluster-cluster collisions and the Ashgriz-Poo model for cluster sizes larger than 100. We found that, by including the combined coalescence model, the average cluster sizes and size distributions predicted by DSMC agree reasonably well with experiment. [ABSTRACT FROM AUTHOR]

In the present work, we have simulated the homogeneous condensation of carbon dioxide and ethanol using the Bhatnagar-Gross-Krook based approach. In an earlier work of Gallagher-Rogers et al. [J. Thermophys. Heat Transfer 22, 695 (2008)], it was found that it was not possible to simulate condensation experiments of Wegener et al. [Phys. Fluids 15, 1869 (1972)] using the direct simulation Monte Carlo method. Therefore, in this work, we have used the statistical Bhatnagar-Gross-Krook approach, which was found to be numerically more efficient than direct simulation Monte Carlo method in our previous studies [Kumar et al., AIAA J. 48, 1531 (2010)], to model homogeneous condensation of two small polyatomic systems, carbon dioxide and ethanol. A new weighting scheme is developed in the Bhatnagar-Gross-Krook framework to reduce the computational load associated with the study of homogeneous condensation flows. The solutions obtained by the use of the new scheme are compared with those obtained by the baseline Bhatnagar-Gross-Krook condensation model (without the species weighting scheme) for the condensing flow of carbon dioxide in the stagnation pressure range of 1-5 bars. Use of the new weighting scheme in the present work makes the simulation of homogeneous condensation of ethanol possible. We obtain good agreement between our simulated predictions for homogeneous condensation of ethanol and experiments in terms of the point of condensation onset and the distribution of mass fraction of ethanol condensed along the nozzle centerline. [ABSTRACT FROM AUTHOR]

Recent direct simulation Monte Carlo (DSMC) simulations of homogeneous condensation in free expansion water plumes [Z. Li, J. Zhong, D. A. Levin, and B. Garrison, AIAA J. 47, 1241 (2009)] show that the nucleation rate is a key factor for accurately modeling condensation phenomenon. In this work, we use molecular dynamics (MD) simulations of a free expansion to explore the microscopic mechanisms of water dimer formation and develop collision models required by DSMC. Bimolecular and termolecular dimer cluster formation mechanisms are considered and the former is found to be the main mechanism in expanding flows to vacuum. MD simulations between two water molecules using the simple point charge intermolecular potential were performed to predict the bimolecular dimer formation probability and the probability was found to decrease with collision energy. The formation probabilities and postcollisional velocity and energy distributions were then integrated into DSMC simulations of a free expansion of an orifice condensation plume with different chamber stagnation temperatures and pressures. The dimer mole fraction was found to increase with distance from the orifice and become constant after a distance of about two orifice diameters. Similar to experiment, the terminal dimer mole fraction was found to decrease with chamber stagnation temperatures and increase linearly with chamber stagnation pressures which is consistent with a bimolecular nucleation mechanism. [ABSTRACT FROM AUTHOR]

In this work, the use of molecular dynamics as a predictive tool for modeling the atomistic behavior of electrospray propulsion is discussed. 1-ethyl-3-methylimidazolium tetrafluoroborate (EMIM-BF4) and ethylammonium nitrate (EAN) were considered as two limits of ionic liquid (IL) propellants that tend to operate in an ion versus a droplet mode. The emission modes were found to depend on the electro-chemical properties of the IL propellant. The aprotic EMIM-BF4-based electrospray emitted primarily monomers and trimers as the dominant species and only small quantities of droplets. In contrast, trimers were the dominant emitted species in the protic EAN emissions with a significantly large contribution from droplets towards the total emission current, suggesting that EMIM-BF4-based colloid thrusters operate in ion mode and EAN-based devices operate in the droplet mode. Furthermore, the formation of the Taylor cone was found to depend on the mass flow rate and the external electric field strength. This paper provides a framework that can be extended for use to simulate any other ILs or their combinations. [ABSTRACT FROM AUTHOR]