Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrostatic”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Role of electrostatic perturbation on kinetic resistive wall mode with application to spherical tokamak

Abstract A more complete non-perturbative magnetohydrodynamic (MHD)-kinetic hybrid formulation is developed by including the perturbed electrostatic potential δφ in the particle Lagrangian. The fluid-like counter-parts of the hybrid equations, in the Chew-Goldberger-Low high-frequency limit, are also derived and utilized to test the new toroidal implementation in the MARS-K code. Application of the updated non-perturbative hybrid model for a high- β spherical tokamak plasma in MAST finds that the perturbed electrostatic potential generally plays a minor role in the n = 1 ( n is the toroidal mode number) resistive wall mode instability. The effect of δφ is largely destabilizing, with the growth rate of the instability increased by several (up to 20) percent as compared to the case without including δφ . A similar relative change is also obtained for the kinetic-induced resonant field amplification effect at high- β in the MAST plasma considered. The updated capability of the MARS-K code allows quantitative exploration of drift kinetic effects on various MHD instabilities and the antenna-driven plasma response where the electrostatic perturbation, coupled to magnetic perturbations, may play important roles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A guide on modeling electrostatics of semiconductor detectors in COMSOL Multiphysics ® for DRiFT

Adding a new semiconductor detector into the detector response function toolkit (DRiFT) requires a model of the electric potential and the electric field. To get an accurate electrostatics model, the geometry of the detector should be modeled as closely to the real geometry as feasible, including the semiconductor materials, doping layers and concentrations, and contacts. Since the electric field and potential are used to calculate the induced signal, the results of the detector response functions are largely influenced by the electrostatics models. Initial electrostatics models were completed using Silvaco, however, this document provides guidance on using COMSOL Mulitphysics ® (COMSOL) to model the detectors to give users additional flexibility. Guidance on using COMSOL and its user interface are mostly left out of this document, but best practices for geometry, modeling methods, and data format/exporting are included to help users generate models that are compatible with DRiFT, prevents bugs, issues, and inaccuracies during charge collection calculations. The following summary provides an overview of the document.

42 ENGINEERING↗
DOE OSTI2023

An implicit, conservative and asymptotic-preserving electrostatic particle-in-cell algorithm for arbitrarily magnetized plasmas in uniform magnetic fields

Here, we introduce a new electrostatic particle-in-cell algorithm capable of using large timesteps compared to particle gyro-period under a uniform external magnetic field. The algorithm extends earlier electrostatic fully implicit PIC implementations with a new asymptotic-preserving particle-push scheme that allows timesteps much larger than particle gyroperiods. In the large-timestep limit, the integrator preserves all particle drifts, while recovering the full orbit for small timesteps. The scheme allows for a seamless, efficient treatment of particles with coexisting magnetized and unmagnetized species, and conserves energy and charge exactly without spoiling implicit solver performance. We demonstrate by numerical experiment with several problems of variable species magnetization (diocotron instability, modified two-stream instability, and drift instability) that orders of magnitude wall-clock-time speedups vs. the standard fully implicit electrostatic PIC algorithm are possible without sacrificing solution accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electrostatic Interactions Control the Nanostructure of Conjugated Polyelectrolyte–Polymeric Ionic Liquid Blends

Polyelectrolyte complexation offers unique opportunities to compatibilize polymers with very different backbone chemistries and to control the morphology of the resulting blend via electrostatic manipulation. In this study, we demonstrate the ability to formulate homogeneous complexes of a conjugated polyelectrolyte with a polymeric ionic liquid, utilizing the electrostatic attraction among their oppositely charged side chains. Variation of electrostatic parameters, such as counterion concentration or polymer charge fraction, tunes the morphology of these polymer complexes from homogeneously disordered blend to weakly structured microemulsion where the local ordering arises from backbone-immiscibility-induced microphase segregation. Our experimental observations are in qualitative agreement with both field-theoretic simulation and random-phase approximation calculations. Simulated morphology snapshots suggest and experimental evidence also indicates that the microphase-segregated complex likely takes on a cocontinuous microemulsion structure. Our findings show that ionic interactions are an effective pathway to compatibilize polymers at macroscopic length scales while achieving controlled nanostructures in these ionic blends. Such systems have great potential for engineering the nanostructure of polymers to tailor applications such as nanofiltration, catalysis, and energy storage, where local ordering can enhance the physical properties of an otherwise macroscopically homogeneous structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Pivotal Roles of Triple Screening-Topological, Electrostatic, and Hydrodynamic-On Dynamics in Semidilute Polyelectrolyte Solutions

For semidilute polyelectrolyte solutions, it is generally assumed that topological, electrostatic, and hydrodynamic interactions are screened (called triple screening). Despite a large body of research focused on polyelectrolyte solutions, the concept of triple screening has never been rigorously verified. In this work, we test the concept by probing concentration fluctuations in aqueous solutions containing a well-studied polyelectrolyte, sodium poly(styrenesulfonate) (NaPSS) with neutron scattering, theory, and molecular dynamics simulations. Neutron spin–echo (NSE) and small-angle neutron scattering (SANS) data from semidilute solutions of NaPSS are presented at different polymer and salt (NaCl) concentrations. A combined theory for structure (J. Chem. Phys. 105, 5183 (1996)) and dynamics (J. Chem. Phys. 107, 2619 (1997)), which captures effects of hydrodynamic, topological, and electrostatic screening, is used to interpret the experimental results. Here, the theory quantitatively predicts the decay rate obtained from the NSE measurements while capturing the shape and concentration dependencies of the polyelectrolyte peak observed in the SANS spectra. Detailed comparisons of the theory and the experiments reveal that the wavevector-dependent decay rate of concentration fluctuations in semidilute solutions of polyelectrolytes is dictated by the screening of hydrodynamic, topological, and electrostatic interactions. This conclusion is corroborated by coarse-grained molecular dynamics simulations, executed without any hydrodynamic interactions, which fail to capture the correct wavevector dependence of the decay rate. These results highlight that the theories based on the concept of triple screening provide a quantitative framework for predicting a relation between the structure and dynamics of polyelectrolyte solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Two-Dimensional Binary Superlattice of BNNT-Surfactant Vesicle Complex Induced by Electrostatic Interaction

For a wide range of practical applications of boron nitride nanotubes (BNNTs), it is essential to achieve their highly ordered self-assembled structures. This study reports on a two-dimensional (2D) binary superlattice of individually exfoliated BNNTs with a negative surface charge (p-BNNT25) and cationic surfactant vesicles (CTAT/SDBS vesicles, prepared by mixing cetyltrimethylammonium tosylate (CTAT) and sodium dodecylbenzenesulfonate (SDBS)) complexes through electrostatic interactions. Depending on the surface charge density of the CTAT/SDBS vesicles and the mass ratio between the CTAT/SDBS vesicle and p-BNNT25, the CTAT/SDBS-BNNT complexes formed highly ordered superstructures. These structures include an intercalated lamellar phase with a centered rectangular structure (ICLP), in which a 2D array of p-BNNT25 is inserted into the multilamellar structure, and an AB 3 structure, in which the BNNTs are surrounded by surfactant micelles in a triangular arrangement. To the best of our knowledge, this is the first demonstration of the fabrication of highly ordered superstructures of individually exfoliated and negatively charged BNNTs with positively charged surfactant vesicles through electrostatic interactions. This approach for the 2D binary superlattices of CTAT/SDBS-BNNT complexes induced by electrostatic interactions is expected to be beneficial for a wide range of one-dimensional (1D) nanoparticle applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Electrostatically controlled spin polarization in Graphene-CrSBr magnetic proximity heterostructures

Abstract The magnetic proximity effect can induce a spin dependent exchange shift in the band structure of graphene. This produces a magnetization and a spin polarization of the electron/hole carriers in this material, paving the way for its use as an active component in spintronics devices. The electrostatic control of this spin polarization in graphene has however never been demonstrated so far. We show that interfacing graphene with the van der Waals antiferromagnet CrSBr results in an unconventional manifestation of the quantum Hall effect, which can be attributed to the presence of counterflowing spin-polarized edge channels originating from the spin-dependent exchange shift in graphene. We extract an exchange shift ranging from 27 – 32 meV, and show that it also produces an electrostatically tunable spin polarization of the electron/hole carriers in graphene ranging from − 50% to + 69% in the absence of a magnetic field. This proof of principle provides a starting point for the use of graphene as an electrostatically tunable source of spin current and could allow this system to generate a large magnetoresistance in gate tunable spin valve devices.

Science & Technology - Other Topics↗
DOE OSTI2024

A contactless in situ EFISH method for measuring electrostatic potential profile of semiconductor/electrolyte junctions

In photoelectrochemical cells, promising devices for directly converting solar energy into storable chemical fuels, the spatial variation of the electrostatic potential across the semiconductor–electrolyte junction is the key parameter that determines the cell performance. In principle, electric field induced second harmonic generation (EFISH) provides a contactless in situ spectroscopic tool to measure the spatial variation of electrostatic potential. However, the total second harmonic generation (SHG) signal contains the contributions of the EFISH signals of semiconductor space charge layer and the electric double layer, in addition to the SHG signal of the electrode surface. The interference of these complex quantities hinders their analysis. In this work, to understand and deconvolute their contributions to the total SHG signals, bias-dependent SHG measurements are performed on the rutile TiO2(100)–electrolyte junction as a function of light polarization and crystal azimuthal angle (angle of the incident plane relative to the crystal [001] axis). A quadratic response between SHG intensity and the applied potential is observed in both the accumulation and depletion regions of TiO2. The relative phase difference and amplitude ratio are extracted at selected azimuthal angles and light polarizations. At 0° azimuthal angle and s-in–p-out polarization, the SHG intensity minimum has the best match with the TiO2 flatband potential due to the orthogonal relative phase difference between bias-dependent and bias-independent SHG terms. We further measure the pH-dependent flatband potential and probe the photovoltage under open circuit conditions using the EFISH technique, demonstrating the capability of this contactless method for measuring electrostatic potential at semiconductor–electrolyte junctions.

Chemistry↗
DOE OSTI2022

Optimizing the Internal Cell Structure of the Radon Monitor Based on Electrostatic Collection Method

Abstract As Rn-222 decays, an alpha particle is emitted and the residual polonium nucleus recoils in the opposite direction. At the end of the recoil path, 88% of the polonium atoms have a positive charge and 12% are neutral. The electric potential distribution in the 60 ml hemispherical internal cell of the radon monitor based on electrostatic collection is studied for reducing the combined probability of the positively charged Po-218 and the OH− produced by the ionization of water vapour in the air. The COMSOL software is used to simulate the electric potential distribution in the internal cell of the radon monitor based on the electrostatic collection method at 27°C, a pressure of 0.1 Mpa. For improving the collection efficiency of Po-218 ions, the average collection time along vertical and oblique lines is calculated when the upper surface of the internal cell is plastic, uncharged metal and charged metal, respectively. Assuming that the gas in the internal cell is uniformly distributed, the results show that if the upper surface of the hemispherical internal cell is plastic or uncharged metal, the electric potential formed in the internal cell is more uniform, and it is beneficial to reduce the total collection time of the positively charged Po-218 ions, thereby improving collection efficiency. The simulation results can be used as an effective reference for optimizing the design of the internal cell structure of the radon monitor based on electrostatic collection method.

Fan, Zhongkai (ORCID:0000000181247459)↗
DOE OSTI2023

GUI Control System for the Mu2e Electrostatic Septum High Voltage at Fermilab

The Mu2e Experiment has stringent beam structure requirements; namely, its proton bunches with a time structure of 1.7 $\mu$s in the Fermilab Delivery Ring. This beam structure will be delivered using the Fermilab 8-GeV Booster, the 8-GeV Recycler Ring, and the Delivery Ring. The 1.7-$\mu$s period of the Delivery Ring will generate the required beam structure by means of a third order resonant extraction system operating on a single circulating bunch. The electrostatic septum (ESS) for this system is particularly challenging, requiring mechanical precision in a ultra high vacuum of 1 x 10$^-8$ Torr to generate 100 kV across 15 mm. This paper describes a graphical user interface that has been developed to automate the conditioning and commissioning process for the electrostatic septa. It is based on an interface to the Fermilab ACNET system using the ACSys Python Data Pool Manager (DPM) Client produced and maintained by Fermilab Accelerator Controls. Network interfacing between data pool managers made by the application and ACNET devices introduce an inherent (approximately 1 s) latency in throughput of the readouts. This delay is utilized to process and graph incoming data events of devices crucial to conditioning of a electrostatic septum (ESS). 'Ramping' and 'Monitoring' modes adjust settings of the power supply based on internal logic to efficaciously increase and maintain the high voltage (HV) in the ESS, easing the voltage setting on incidence of sparking or other possibly damaging events. A timestamped log file is produced as the application runs.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Thermodynamic Model and Initial Experimental Investigation of Air Dehumidification through Electrically Charged Vapor Capturing Electrostatic Droplets

The removal of water vapor from the air to reduce relative humidity is a well known indoor environmental comfort requirement. Common dehumidification approaches require a substantial amount of energy and usually involve the cooling of atmospheric humid air below its dew point or the use of absorbent/adsorbent materials to extract water vapor out of the air. More recently, researchers investigated the effect of electrostatic forces for enhancing water vapor condensation with the goal to reduce the energy consumption associated with dehumidification. However, the studies are limited, and there is a lack of correlations that can predict the dehumidification rate. In this thesis, a broad theoretical investigation to study the electrostatically enhanced condensation processes was carried out. These processes consist of the use of highly charged particles, preferably highly charged water droplets to attract polar water vapor molecules to their surfaces and promote condensation, a phenomenon known as dielectrophoresis. The electrical charge promotes the reduction of the vapor pressure on the droplets' surface with respect to the saturated pressure predicted by the Kelvin equation on curved surfaces and, consequently, the equilibrium between evaporation and condensation is shifted towards condensation. This investigation resulted in the development of a thermodynamics model, which was able to predict an effective size range of the charged droplets for optimal dehumidification, under ideal conditions. The range resulted in about 2 to 4 m in diameter, while the electrical charge was kept to the maximum limit predicted by the Rayleigh model. A sensitivity analysis on the variation of the size and the charge of the charged droplets was also performed. In terms of dehumidification rates, when six electrospray heads, i.e. the maximum number of heads tested in the experimental campaign, were considered, the model predicted very limited rates, way lower than the target of this research, i.e. a dehumidification rate in terms of relative humidity of 5 % with an air flow rate of 5 cfm. While ways to increase this rate existed, their implementation proved difficult for this initial investigation. According to the requirements predicted by the model, the use of electrosprays appeared the most suitable solution for the production of small but highly charged droplets. The electrospray features and operational modes were studied in detail and an electrospray assembly was designed to be tested in this initial experimental investigation. For this initial investigation, an experimental system was designed and built to validate the thermodynamics model and with the goal to achieve a dehumidification rate of 5 % in terms of relative humidity with an air flow rate of 5 cfm. The experimental system consisted of a wind tunnel test apparatus, dew point sensors, thermocouples for temperature measurements, pressure sensors and other devices. The experimental system was equipped with a computerized data acquisition and storage system. The system was suitable to evaluate the air water content differential before and after the test section where the electrospray heads were installed. The air water content differential was evaluated in terms of the ∆ω/ω_1 ratio, to eliminate the dependence on the dry bulb temperature along the test apparatus. This work presented initial experimental data for electrostatically enhanced dehumidification processes. Working parameters, such as air and water flow rates, high voltage potentials and polarity, deionized water types with different electrical conductivity, were varied to find the best combination for an improved dehumidification. With the set up used and the conditions considered, it was never possible to achieve a 5 % dehumidification rate with 5 cfm. In general, the dehumidification rate was always limited and lower than 1 % for all the air flow rates considered and mainly within the uncertainty of the dew point sensors. These results confirmed the already limited predictions of the thermodynamics model, for ideal conditions and the interesting and promising results obtained in some experimental investigations available in the open literature were not achievable in this thesis. The more probable reasons were the higher air flow rates considered and the complete absence of the use of cooling power, which was mainly implemented to facilitate the vapor condensation.

Morcelli, Stefano↗
DOE OSTI2023

Saddle-Shaped Third Component with Out-of-Plane Electrostatic Dipole for Realizing High-Performance Photovoltaic Donor Terpolymers

Terpolymer fabrication is an effective methodology for molecular engineering and generating high-performance organic photovoltaic materials to construct highly efficient polymer solar cells. Modification of the polymer PM6 by incorporating a third component resulting in the formation of a ternary copolymer is reported to outperform PM6 in achieving enhanced device performances. However, one of the major challenges in constructing high-performance terpolymers is to counter the molecular disorder caused by the backbone entropy induced by the third moiety. Here, in this work, double B←N bridged bipyridine (BNBP) is used as the third component, which possesses a strong out-of-plane electrostatic dipole owing to the saddle-shaped B←N fused ring structure. The out-of-plane dipole moment introduced in the modified PM6 terpolymer can be used as a means for tuning and optimizing the nanostructures of the blended films. The prepared PM6-BNBP-4 blend polymer with 4% of the benzodithiophene dione monomers replaced by BNBP results in excellent power conversion efficiency of 19.13%. This work demonstrates that the out-of-plane electrostatic dipole moment in saddle-shaped molecules is valuable for achieving high-performance organic photovoltaic donor materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electrostatic powder spreading for metal powder bed fusion applications

Typical powder bed fusion processes spread each layer of powder by pushing it across the build area with a spreader bar or roller. While simple, these methods require good flowing powder and do not have the ability to place more than one powder type per layer. This naturally limits the ability to incorporate more than one material onto each powder layer in powder bed additive manufacturing processes. Here we demonstrate electrostatic powder spreading, a method to deposit a powder layer without mechanical contact with the previous layers. Electrostatic powder spreading (ESPS) uses an electric field formed between a powder container and a counter electrode to move powder particles from the container to the powder bed. We show that the rate at which powder is deposited can be controlled by adjusting the strength of the electric field. Correspondingly, the deposited layer thickness is controlled by the strength of the electric field and the recoating speed. To demonstrate the feasibility of ESPS for laser powder bed fusion powder spreading, nearly fully dense (> 99.8%) coupons were built in a commercial L-PBF machine using the ESPS setup. Finally, we achieve 2D patterning of a single layer by using an array of electrodes and switching the applied electric field on and off. ESPS represents a method to produce an arbitrary patterned multi-material powder bed without complicated powder removal steps and easy implementation into existing L-PBF machines. Furthermore, ESPS presents an opportunity to build parts with 3D gradients in material properties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Interplay of Pauli Repulsion, Electrostatics, and Field Inhomogeneity for Blueshifting and Redshifting Vibrational Probe Molecules

Many molecules' vibrational frequencies are sensitive to intermolecular electric fields, enabling them to probe the field in complex molecular environments. However, it is often unclear whether the probe is responding to the local electric field or other types of intermolecular interactions, inhibiting interpretation of the frequency and effectiveness as probes. This is especially true for molecules whose vibrational frequencies blueshift instead of the more typical redshift in hydrogen bonding configurations. In this work, we computationally investigate the causes of redshifting versus blueshifting over a range of vibrational reporters. First, we apply adiabatic energy decomposition analysis to a paradigmatic set of probes, finding that redshifting only occurs when electrostatic interactions are strong enough to overcome the dominant and large blueshifting contribution of Pauli repulsion. Furthermore, we demonstrate that field inhomogeneity can further shift the frequency of many probes substantially to either reinforce or counteract the shift expected from a homogeneous field. We find that redshifting is reinforced by electric field inhomogeneity, otherwise field inhomogeneity further weakens the electrostatic contribution relative to Pauli repulsion, leading to blueshifting. Further calculations indicate that the probe's response to field inhomogeneity can be understood by considering the mass of the atoms involved in the stretching mode and sign of the electric field. In explaining the interplay of different intermolecular interactions and field inhomogeneity for many probes, our results should enable the use and interpretation of spectroscopic probes and their connection to electric fields in more complex systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Revealing Local and Directional Aspects of Catalytic Active Sites by the Nuclear and Surface Electrostatic Potential

This work examines the prospects of using the electrostatic potential, V(r), as a descriptor in heterogeneous catalysis. In particular, the subatomic spatial resolution of the property allows for analysis of both directionality and confinement effects in surface adsorption. This feature of V(r) is used to identify adsorption sites, orientations, and energetics for metal surfaces, particles, and nanoclusters upon interactions with catalytically relevant intermediates. The use of V(r) in assessing the 3D nature of catalytic sites in low-temperature and electrocatalysis is highlighted, and future directions in catalysis design are discussed. Ultimately, we provide a critical analysis of the use of V(r) in the predictions of local adsorption susceptibilities, and we address its limitations. The link between V(r) and other established descriptors in catalysis are motivated via physical relations and theoretical derivations; close ties are established between V(r) and the d-band center (ε d ), as well as the surface site stability (BE M ). In conclusion, by comparing the performance of V(r) evaluated on isodensity contours, i.e., the surface electrostatic potential, to that of V(r) evaluated at the nucleus of an atom, we investigate the application space for a directional and an atom-localized version of the V(r) descriptor for catalyst design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Direct Imaging of Asymmetric Interfaces and Electrostatic Potentials inside a Hafnia–Zirconia Ferroelectric Nanocapacitor

In hafnia-based thin-film ferroelectric devices, chemical phenomena during growth and processing, such as oxygen vacancy formation and interfacial reactions, appear to strongly affect device performance. However, the correlation between the structure, chemistry, and electrical potentials at the nanoscale in these devices is not fully known, making it difficult to understand their influence on device properties. Here, we directly image the composition and electrostatic potential with nanometer resolution in the cross section of a nanocrystalline W/Hf 0.5 Zr 0.5 O 2−δ (HZO)/W ferroelectric capacitor using multimodal electron microscopy. This reveals a 1.4 nm wide tungsten suboxide interfacial layer formed at the bottom interface during fabrication, which introduces a potential dip and leads to asymmetric switching fields. Additionally, we compare the measured potentials to DFT calculations and find it is nearly 3 V lower than expected in the HZO, which appears to be caused by oxygen vacancies and a resulting negative built-in potential. In conclusion, these chemical and electrostatic details are important to characterize and tune to achieve high-performance ferroelectric devices.

Defects in solids↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗