Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic simulations”

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 127 records · Page 7

Simulation of fundamental atomization mechanisms in fuel sprays

Growth of instabilities on the liquid/gas interface in the initial region of fuel sprays is studied by means of numerical simulations. The simulations are based on solutions of the variable-density incompressible Navier-Stokes equations, which are obtained with a new numerical algorithm. The simulations give good agreement with analytical results for the instabilities on a liquid cylinder induced by surface tension and wind-induced instabilities. The effects of boundary layers on the wind-induced instabilities are investigated. It is found that a boundary layer reduces the growth rate for a single interface, and a comparison with inviscid theory suggests that boundary layer effects may be significantly more important than surface tension effects. The results yield a better estimate than inviscid theory for the drop sizes as reported for diesel sprays. Results for the planar jet show that boundary layer effects hasten the growth of Squire's 'symmetric' mode, which is responsible for jet disintegration. This result helps explain the rapid atomization which occurs in swirl and air-blast atomizers.

Childs, Robert, E.↗

Kinetic Simulation and Energetic Neutral Atom Imaging of the Magnetosphere

Advanced simulation tools and measurement techniques have been developed to study the dynamic magnetosphere and its response to drivers in the solar wind. The Comprehensive Ring Current Model (CRCM) is a kinetic code that solves the 3D distribution in space, energy and pitch-angle information of energetic ions and electrons. Energetic Neutral Atom (ENA) imagers have been carried in past and current satellite missions. Global morphology of energetic ions were revealed by the observed ENA images. We have combined simulation and ENA analysis techniques to study the development of ring current ions during magnetic storms and substorms. We identify the timing and location of particle injection and loss. We examine the evolution of ion energy and pitch-angle distribution during different phases of a storm. In this talk we will discuss the findings from our ring current studies and how our simulation and ENA analysis tools can be applied to the upcoming TRIO-CINAMA mission.

Fok, Mei-Ching H.↗

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser Transport System Vacuum Simulations and LED Atom Tracker

MAGIS-100 is an experiment using atom interferometry on strontium atoms in a 100m vacuum chamber. The laser travels through a laser transport system (LTS) before going down the chamber. Two projects were completed over the course of the internship. One project was conducting vacuum simulations on the LTS to ensure it meets the experimental pressure requirement of 10-11 torr. This was done by characterizing the pressure profile with variations in pump size, pump spacing, and orifice size. The initial design of the LTS included one 30 L/s ion pump in the high vacuum (HV) region, a lens mount with holes of diameter 0.788”, and three 150 L/s ion pumps in the ultra-high vacuum (UHV) region. After adjusting these variables, it was concluded that the lens mount holes could be reduced to 0.200” in diameter and one 150 L/s ion pump in the UHV could be removed and still meet the experimental requirements. The other project was to design an LED atom tracker that shows the position of the strontium atoms in the 100m shaft. This will help with public outreach by allowing visitors to have a visual of the experiment. An Arduino Uno and Arduino Due were programmed using classical physics formulas to direct the RGB LED strips to flash in sync with the falling atoms. The requirements for powering and connecting the hardware were calculated and recommendations for scaling the system to 100m were made.

43 PARTICLE ACCELERATORS↗

Atomic Oxygen (ATOX) simulation of Teflon FEP and Kapton H surfaces using a high intensity, low energy, mass selected, ion beam facility

A high intensity (greater than 10(exp 15) ions/sq cm) low energy (down to 5 eV) mass selected ion beam (MSIB) facility was used to study the effects of ATOX on two polymers commonly used for space applications (Kapton H and Teflon FEP). The polymers were exposed to O(+) and Ne(+) fluences on 10(exp 15) - 10(exp 19) ions/sq cm, using 30eV ions. A variety of analytical methods were used to analyze the eroded surfaces including: (1) atomic force microscopy (AFM) for morphology measurements; (2) total mass loss measurements using a microbalance; (3) surface chemical composition using x-ray photoelectron spectroscopy (XPS), and (4) residual gas analysis (RGA) of the released gases during bombardment. The relative significance of the collisional and chemical degradation processes was evaluated by comparing the effects of Ne(+) and O(+) bombardment. For 30 eV ions it was found that the Kapton is eroded via chemical mechanisms while Teflon FEP is eroded via collisional mechanisms. AFM analysis was found very powerful in revealing the evolution of the damage from its initial atomic scale (roughness of approx. 1 nm) to its final microscopic scale (roughness greater than 1 micron). Both the surface morphology and the average roughness of the bombarded surfaces (averaged over 1 micron x 1 micron images by the system's computer) were determined for each sample. For 30 eV a non linear increase of the Kapton roughness with the O(+) fluence was discovered (a slow increase rate for fluences phi less than 5 x 10(exp 17) O(+)/sq cm, and a rapid increase rate for phi greater than 5 x 10(exp 17) O(+)/sq cm). Comparative studies on the same materials exposed to RF and DC oxygen plasmas indicate that the specific details of the erosion depend on the simulation facility emphasizing the advantages of the ion beam facility.

Vered, R.↗

All-Atom Molecular Dynamics Simulations of the Temperature Response of Densely Grafted Polyelectrolyte Brushes

Water acts a “good” solvent for most polyelectrolyte (PE) systems due to its high dielectric strength. Theoretical and experimental studies have highlighted some remarkably unconventional properties of strongly water-swollen PE brushes such as a reduction in brush height with increase in temperature. However, an understanding of the atomistic-scale response of the grafted PE chains as well as the local arrangement and organization of the brush-supported counterions and water molecules to changes in temperature is missing. Here, we conduct an all-atom molecular dynamics (MD) study to probe the influence of temperature on the “microstructure” of densely grafted, fully ionized polyacrylic acid (PAA) brushes. The atomistic insights obtained from our study shed light on prior experimental observations and elucidate the key role played by changes in hydrogen bonding network in dictating the hydrophilicity of the brushes. In addition to the temperature-mediated reduction in height of the PE brush layer, we study the influence of temperature on several properties of the brush-supported water molecules such as the number distribution and kinetics of water-water and water-PE hydrogen bonds (HBs), mass density, orientational tetrahedral order parameter, and the activation energy associated with the self-diffusion of water. The effect of temperature on the properties of the brush-trapped counterions is also quantified via changes in their solvation structure, counterion-water radial distribution function (RDF), and the translational mobility. Our findings unravel several hitherto unknown phenomena such as the enhanced propensity for attaining PE brush-induced “water-in-salt”-like scenarios at elevated temperatures, and an increase in the activation energy for the self-diffusion of water with increase in the degree of PE brush-induced nanoconfinement (from ‘bulk water’ to ‘sparsely grafted PE brushes’ to ‘densely grafted PE brushes’). Further, we anticipate that this work will be helpful in improving our predictive capabilities of the temperature response of water-swollen PE brushes, which would be instrumental for better design of several PE brush systems used in a multitude of applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Grain Boundary Misorientation on Spall Strength in Ta via Shock-Free Simulations with Relatively Few Atoms

A suite of 37 molecular dynamics simulations is conducted at two system sizes to systematically characterize the role of grain boundary (GB) misorientation on spall strength in pure BCC tantalum (Ta). The systems studied consist of bicrystals with a single [110] symmetric tilt grain boundary. Two loading conditions are compared: (i) homogeneous extension under uniaxial strain simulated in this study and (ii) piston/flyer impact of sample, which induces heterogeneous deformation via shockwave propagation along the length of the sample. The piston/flyer impact is taken from the literature and run on the same set of GB misorientation angles using LAMMPS. The major finding here is that both methods result in similar spall strength predictions, but the homogeneous extension method generally requires two to three orders of magnitude fewer atoms and similar reductions in computational costs. Spall strength results systematically overpredict using this method, by about 10% for the dataset three orders of magnitude smaller than piston/flyer simulations, and 5% for the dataset two orders of magnitude smaller. Lastly, the effect of system size and pre-compression magnitude on spall strength is systematically characterized.

36 MATERIALS SCIENCE↗

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗

All-Atom Modeling and Simulation of Biopolymer Interface: Dual Role of Antifouling Polymer Brushes

Antifouling polymer brushes are well-known for their exceptional resistance to unwanted protein adsorption. While experimental studies have extensively characterized protein–polymer brush interactions, computational investigations remain limited, largely due to the challenges in accurately modeling and integrating these complex, multicomponent systems without resorting to oversimplification. To address this challenge, this study presents one of the most comprehensive and realistic model systems to date, comprising the substrate, grafted polymers, and proteins. Our work interprets the interactions between polymer brush and protein based on realistic modeling without simplification. In particular, this study utilizes molecular modeling and simulation of polycationic and polyzwitterionic brushes─poly(dimethylaminoethyl methacrylate) (PDMAEMA), poly(2-(N-oxide-N,N-dimethylamino)ethyl methacrylate) (PNOMA), and poly(2-(N-3-sulfopropyl-N,N-dimethylammonium)ethyl methacrylate) (PSBMA)─grafted onto α-quartz substrates via polymerization initiator linkers. The brush models were developed to closely replicate experimentally synthesized samples and to provide detailed insights into structural and dynamical changes at the molecular level during protein adsorption. Using steered molecular dynamics simulations, we show that the PSBMA brush, due to its high local density, exhibits the greatest resistance to protein insertion. Cα root-mean-square deviation and interaction pattern analyses further reveal that the PSBMA brush also induces the most significant destabilization of lysozyme, while the PDMAEMA brush enhances protein stability through ion-mediated interactions. The PNOMA brush, while requiring the lowest force for protein adsorption, induces greater protein destabilization than the PDMAEMA brush, primarily due to electrostatic repulsion caused by a short carbon spacer length. Hydration analysis reveals that both the PSBMA brush and the lysozyme interacting with it exhibit the most rapid dehydration, attributed to the brush’s high local chain density, which results in the greatest lysozyme destabilization and the highest adsorption force. These findings highlight the dual role of antifouling polymer brushes: resisting protein adsorption and modulating protein structural dynamics. In conclusion, this study provides valuable insights for the rational design of next-generation antifouling materials and offers a framework for realistic model development in complex multicomponent systems.

Adsorption↗

Solvent-induced membrane stress in biofuel production: molecular insights from small-angle scattering and all-atom molecular dynamics simulations

The disruptive effect of organic solvents on microbial membranes represents a significant challenge to the economical production of green fuels and value-added chemicals from lignocellulosic feedstocks. One route to overcoming this challenge is to engineer microbes with membranes capable of resisting organic solvent stresses. In this regard, it is useful to understand the mechanisms by which organic solvents disrupt typical biomembranes. In this study, molecular dynamics (MD) simulation, complemented by small-angle X-ray and neutron scattering (SANS/SAXS), provide a molecular-scale view of the disruption of a microbial model membrane by 1-butanol and tetrahydrofuran (THF), two common water–organic cosolvent mixtures of importance in biofuel production. Solvent interactions at the interface between the head-group and fatty acid tail regions lead to more dramatic membrane changes than interactions solely at the head-groups or tails. Although both organic solvents are found to partition into the membrane, the depth of solvent penetration into the membrane is quite different. Specifically, 1-butanol localizes near the interface between the lipid heads and tails at low concentrations, but partitions into both the head and tail regions at high concentrations. In contrast, THF, overall, partitions less than 1-butanol and prefers the lipid tail regions. Importantly, the presence of 1-butanol near the head/tail interface introduces drastic membrane changes not seen with THF. The organic solvent interactions with the lipids lead to membrane thinning and fluidization, but more so for 1-butanol than for THF. These results suggest that an aim for the future engineering of robust membranes could be to design lipid head groups that reduce the accumulation of organic solvents at the head–tail interface and that rational designs need also be cognizant of the different solvent-specific mechanisms responsible for membrane disruption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Laser Transport Vacuum Simulations and LED Atom Tracker

MAGIS-100 uses atom interferometry to measure the distance between the Sr atoms traversing the shaft. The laser must travel through vacuum to measure accurately. Project ensures the vacuum system has low enough pressure to meet experimental requirements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ultra-high energy density relativistic plasmas from nanostructures: scaling to ultra-high intensities (Final Report)

This project investigated the ultra-high energy density (UHED) regime of matter found in the center of stars, using a compact PW-class laser. These extreme conditions are typically only obtained in the laboratory in the central hot-spot of spherically imploded capsules in inertial confinement fusion experiments driven by the world’s largest lasers. We have shown that near-solid density arrays of aligned nanostructures can be volumetrically heated to multi-keV temperatures by irradiation at relativistic intensity with ultrafast laser pulses of modest energy, opening a path for the generation of UHED plasmas with compact lasers. This new UHED plasma generation approach promises to create an environment with extreme energy densities and degrees of ionization, record conversion of optical laser light into ultrafast x-ray pulses, gigantic magnetic fields and pressures, and directed beams of high energy particles . In this project we achieved record degree of ionization in volumetric heating solid density and near-solid density plasmas. Gold plasmas which spectra is characterized by L shell transition emission from ions with charge up to the Ne-like state, Au+ 72 were generated using ultrafast laser pulses of less than 10 J of energy from a compact laser focused to an intensity of ~ 3x10 21 Wcm -2 . We also conducted measurements to determine the heat penetration depth in Ni nanowire arrays as compared to Ni foil targets by monitoring the line emission of a Co buried tracer underneath a variable amount of Ni. The measurements revealed that the nanowire plasmas are roughly six times larger in depth than the solid density target plasmas. A result of this increased plasma in nanowire arrays is a greatly increased conversion of optical laser light into > 1 KeV x-rays, a record conversion efficiency of 20 %. Critical to the realization of the proposed experiments was the generation of ultrafast laser pulses with ultra-high contrast that can deposit the energy deep into the nanowire arrays before the nanowires explode to create a continuous plasma. Supporting the proposed experiments was the recent demonstration at Colorado State University of a Petawatt-class laser that emits 30 fs pulses at high repetition rates. The experiments combined this unique laser tool with a variety of tailored nanowire arrays fabricated in house and with an extensive suite of diagnostics. The experiments were accompanied by 3-dimensional particle-in-cell simulations and detailed atomic physics simulations with transient kinetics and radiation transport. The proposed research allowed us to continue training Ph.D students and post-docs with broad experience in HEDP Science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗