Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Density”

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 1,585 records · Page 88

Carbon Nanotubes: Molecular Electronic Components

The carbon Nanotube junctions have recently emerged as excellent candidates for use as the building blocks in the formation of nanoscale molecular electronic networks. While the simple joint of two dissimilar tubes can be generated by the introduction of a pair of heptagon-pentagon defects in an otherwise perfect hexagonal graphene sheet, more complex joints require other mechanisms. In this work we explore structural characteristics of complex 3-point junctions of carbon nanotubes using a generalized tight-binding molecular-dynamics scheme. The study of pi-electron local densities of states (LDOS) of these junctions reveal many interesting features, most prominent among them being the defect-induced states in the gap.

Srivastava, Deepak↗

Scalability of a Low-Cost Multi-Teraflop Linux Cluster for High-End Classical Atomistic and Quantum Mechanical Simulations

Scalability of a low-cost, Intel Xeon-based, multi-Teraflop Linux cluster is tested for two high-end scientific applications: Classical atomistic simulation based on the molecular dynamics method and quantum mechanical calculation based on the density functional theory. These scalable parallel applications use space-time multiresolution algorithms and feature computational-space decomposition, wavelet-based adaptive load balancing, and spacefilling-curve-based data compression for scalable I/O. Comparative performance tests are performed on a 1,024-processor Linux cluster and a conventional higher-end parallel supercomputer, 1,184-processor IBM SP4. The results show that the performance of the Linux cluster is comparable to that of the SP4. We also study various effects, such as the sharing of memory and L2 cache among processors, on the performance.

Kikuchi, Hideaki↗

Modeling the Extended Neutral Atmosphere and Plasma Environment near Saturn

In the three years of this study we have published five papers in refereed journals. We have first examined satellite sources and their contribution to the observed neutral cloud. Based on the total calculated satellite sources and the spatial distribution of sputtered H20, we concluded that they cannot produce observed OH cloud. The E-ring contribution has been also studied in detail.In order to produce observed OH cloud we suggested that the E-ring might be the dominant source in inner Saturnian magnetosphere. We proposed a possible resupply mechanism which is needed to keep both E-ring and OH cloud in the present state: collisions between E-ring grains and remains of a disrupted satellite near Enceladus Lagrangian points. In this scenario a large amount of material, ranging from a few micrometers to hundred of meters, which is optically invisible at present, is likely to exist there. The fourth paper compares the magnetosheaths of the outer planets. A surprising result is that the hot proton component comprises about 40% of the total density, much larger than predicted by shock theory. Gas dynamic models of the boundaries show that the magnetospheres of Jupiter and Saturn are flattened at the poles. The last paper was published in GRL and is the first based of the model of neutrals developed as a main goal of this grant and which is now operational.This Monte Carlo collision code self- consistently determines the neutral distribution. from the rings and satellites until they are lost by ionization, by collisions with rings, moons, or Saturn, or by escape from Saturn. Our model is unique in that it includes the effects of plasma chemistry and both plasma-neutral and neutral neutral collisions to determine the dynamical evolution of the water group neutrals in Saturn's magnetosphere. The lifetimes of the neutrals against loss to photoionization, charge exchange, electron dissociation and electron impact dissociation are based on the model given by (Richardson et al. 1998) and vary with position in the magnetosphere. The dominant neutral dissociation channels H20->OH+H, H20-> O+H2, and OH->O+H2 are considered.

Richardson, John D.↗

Surface Time-Variable Gravity Signals and Possible Sources Including Core Mass Flow

Over two decades of geodetic satellite-laser-ranging (SLR) data show that the variation of the Earth's oblateness parameter J2 has a clear seasonal signal of amplitude of about 3e-10 and a secular decrease of about -2.8e-11/year, superimposed on some interesting interannual fluctuations. Physically, any change in mass distribution or/inside the Earth will be reflected in the time-variable gravity signal obtained outside the Earth, according to Newton s gravitational law. Therefore, such signal contains contributions from all geophysical sources that redistribute mass, on all temporal and spatial scales, including those from the core. Besides Earth rotation and geomagnetic field variations, the time-variable gravity also contains information linking Earth surface observations with internal core dynamical processes. The time scales of the gravity signal are critical in helping differentiate different contributions. The atmosphere and hydrosphere are responsible for the seasonal and much of the interannual and intraseasoanl fluctuations, while the secular trend is due mainly to the post-glacial rebound but possibly core mass flow. To estimate the latter effect, we use our MoSST (Modular, Scalable, Self-consistent, Three-dimensional) core dynamics model to forward simulate the core flow, and density variation due to the core convection. Our results suggest that, when upward continued to the surface, the J2 component of the core mass redistribution can reach an overall amplitude of e-11/year, approaching the SLR detectability and significant in geophysical terms. We also find a general westward drift of the mass flow, with a speed comparable to that of the geomagnetic westward drift.

Chao, Benjamin F.↗

LDSD POST2 Simulation and SFDT-1 Pre-Flight Launch Operations Analyses

The Low-Density Supersonic Decelerator (LDSD) Project's first Supersonic Flight Dynamics Test (SFDT-1) occurred June 28, 2014. Program to Optimize Simulated Trajectories II (POST2) was utilized to develop trajectory simulations characterizing all SFDT-1 flight phases from drop to splashdown. These POST2 simulations were used to validate the targeting parameters developed for SFDT- 1, predict performance and understand the sensitivity of the vehicle and nominal mission designs, and to support flight test operations with trajectory performance and splashdown location predictions for vehicle recovery. This paper provides an overview of the POST2 simulations developed for LDSD and presents the POST2 simulation flight dynamics support during the SFDT-1 launch, operations, and recovery.

Bowes, Angela L.↗

Magnus Effect on a Spinning Satellite in Low Earth Orbit

A spinning body in a flow field generates an aerodynamic lift or Magnus effect that displaces the body in a direction normal to the freestream flow. Earth orbiting satellites with substantial body rotation in appreciable atmospheric densities may generate a Magnus force to perturb orbital dynamics. We investigate the feasibility of using this effect for spacecraft at a perigee of 80km using the Systems Tool Kit (STK). Results show that for a satellite of reasonable properties, the Magnus effect doubles the amount of time in orbit. Orbital decay was greatly mitigated for satellites spinning at 10000 and 15000RPM. This study demonstrates that the Magnus effect has the potential to sustain a spacecraft's orbit at a low perigee altitude and could also serve as an orbital maneuver capability.

Orbit↗

Simultaneous Conventional and Plenoptic Background Oriented Schlieren Imaging

Plenoptic Background Oriented Schlieren (BOS) is an emerging schlieren technique that is capable of providing 3D qualitative and quantitative information about density gradients present in a wide range of fluid dynamics problems. In this work, the fundamental concepts of plenoptic BOS are reviewed before discussing an open-air experiment with a buoyant plume where both conventional BOS and plenoptic BOS measurements were acquired simultaneously. Both cameras had the same field-of-view for all experiments, and three different focal plane arrangements were explored: (1) the focal plane was set to the background positions and the plume varied between 11 different positions relative to this focal plane, (2) the focal plane was set to 635- millimeters in front of the background position, and (3) the nominal focal plane varied while the position of the plume remained fixed. Such discussion will provide insight on how the two techniques compare, and what additional work is required to better understand the results provided by these two imaging systems.

Klemkowsky, Jenna N.↗

Redshifting the Cosmological Constant in Unimodular Gravity via Nonlinear Quantum Mechanics

The cosmological constant problem represents a profound conflict between quantum field theory and general relativity. Unimodular gravity offers a compelling starting point by de-gravitating the vacuum energy of the Standard Model, but this framework traditionally trades the problem of vacuum energy for a fine-tuning of initial conditions, which manifest as a ``shadow" cosmological constant. In this paper, we resolve this initial conditions problem by proposing a novel modification to gravity based on nonlinear quantum mechanics. We introduce specific state-dependent terms to the Hamiltonian, constructed from expectation values of the metric such as the average Ricci scalar. These terms alter the dynamical equations of gravity such that the shadow energy density associated with unconstrained initial conditions redshifts away with cosmic expansion, rendering it negligible at late times. The resulting cosmology is naturally dominated by matter and radiation without fine-tuning. We demonstrate that this significant infrared modification of gravity is consistent with local and cosmological tests of gravity. We comment on the possibility of testing this solution in cosmological measurements of Newton's constant.

Kaplan, David E. [Johns Hopkins U.; Tokyo U., IPMU↗

POLAR Observations of Topside Field-Aligned O+ Flows and Auroral Forms

Measurements of thermal O (sup +) ion densities, field-aligned velocities, and fluxes from the Thermal Ion Dynamics Experiment (TIDE) on POLAR obtained near 5000 km altitude over the Southern hemisphere are compared with auroral images from the Ultra Violet Imager (UVI). We find upward O (sup +) flows in the cleft region, but subsonic O (sup +) downflows in the polar cap region. Also, the O (sup +) ion density follows a decreasing trend from the poleward side of the cusp region into the nightside aurora region. The magnitude of the downward O (sup +) parallel velocities increases from dayside to nightside across the polar cap boundary. The upflows tend to occur over or near auroral forms, while the downflows are seen in relatively dark regions, such as the polar cap. These results are consistent with a cleft ion fountain source for the polar cap O (sup +) ions. In the nightside polar cap, the results indicate a transition from downward to upflowing field-aligned O (sup +) ions near boundaries of bright auroral arcs.

Stevenson, B. A.↗

Plasmon-induced excitation energy transfer in silver nanoparticle dimers: A real-time TDDFTB investigation

Here, using real-time quantum dynamics calculations, we perform theoretical investigations of light-induced interactions and electronic excitation transfer in a silver nanoparticle dimer. Real-time time-dependent density functional tight-binding (RT-TDDFTB) calculations provide details of the quantum dynamical processes at an electronic/atomistic level with attosecond resolution. The computational efficiency of RT-TDDFTB allows us to examine electronic dynamics up to picosecond time scales. With time scales varying over six orders of magnitude, we provide insight into interactions between the nanoparticle and laser and between nanoparticles. Our results show that the coupling between nanoparticle monomers is dependent on the separation distance between the nanoparticles in the dimer. As the interparticle distance is varied, the dipole–dipole interactions and electronic excitation transfer mechanisms are markedly different. At large distances (from 50 to 20 Å), the energy transfer from NP1 to NP2 becomes more efficient as the interparticle distance decreases. The total dipole moment of the Ag 14 nanoparticle dimer increases linearly at an interparticle distance of 20 Å and reaches its maximum after 1.2 ps. The electronic excitation transfer is also the most efficient at 20 Å. At short distances, back-transfer effects reduce the ability of the dimer and NP1 to accept energy from the incident electric field. We attribute the distance-dependent features of the nanoparticle dimer to the beating between the laser acting on NP1 and the back transfer from NP2 to NP1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Zero-cost corrections to influence functional coefficients from bath response functions

Recent work has shown that it is possible to circumvent the calculation of the spectral density and directly generate the coefficients of the discretized influence functionals using data from classical trajectory simulations. However, the accuracy of this procedure depends on the validity of the high temperature approximation. In this work, an alternative derivation based on the Kubo formalism is provided. This enables the calculation of additional correction terms that increases the range of applicability of the procedure to lower temperatures. Because it is based on the Kubo-transformed correlation function, this approach allows the direct use of correlation functions obtained from methods such as ring-polymer molecular dynamics and centroid molecular dynamics in determining the influence functional coefficients for subsequent system-solvent simulations. The accuracy of the original procedure and the corrected procedure is investigated across a range of parameters. It is interesting that the correction term comes at zero additional cost. Furthermore, it is possible to improve upon the correction using zero-cost physical intuition and heuristics making the method even more accurate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Learning compact physics‐aware delayed photocurrent models using dynamic mode decomposition

Abstract Radiation‐induced photocurrent in semiconductor devices can be simulated using complex physics‐based models, which are accurate, but computationally expensive. This presents a challenge for implementing device characteristics in high‐level circuit simulations where it is computationally infeasible to evaluate detailed models for multiple individual circuit elements. In this work we demonstrate a procedure for learning compact delayed photocurrent models that are efficient enough to implement in large‐scale circuit simulations, but remain faithful to the underlying physics. Our approach utilizes dynamic mode decomposition (DMD), a system identification technique for learning reduced‐order discrete‐time dynamical systems from time series data based on singular value decomposition. To obtain physics‐aware device models, we simulate the excess carrier density induced by radiation pulses by solving numerically the ambipolar diffusion equation, then use the simulated internal state as training data for the DMD algorithm. Our results show that the significantly reduced‐order delayed photocurrent models obtained via this method accurately approximate the dynamics of the internal excess carrier density—which can be used to calculate the induced current at the device boundaries—while remaining compact enough to incorporate into larger circuit simulations.

Hanson, Joshua↗

Mechanical and Thermal Properties for Uranium and U–6Nb Alloy from First-Principles Theory

Elasticity, lattice dynamics, and thermal expansion for uranium and U–6Nb alloy (elastic moduli) are calculated from density functional theory that is extended to include orbital polarization (DFT+OP). Introducing 12.5 at.% of niobium, substitutionally, in uranium softens all the cii elastic moduli, resulting in a significantly softer shear modulus (G). Combined with a nearly invariant bulk modulus (B), the quotient B/G increases dramatically for U–6Nb, suggesting a more ductile material. Lattice dynamics from a harmonic model coupled with a DFT+OP electronic structure is applied for α uranium, and the obtained phonon density of states compares well with inelastic neutron-scattering measurements. The Debye temperature associated with the lattice dynamics falls within the range of experimentally observed Debye temperatures and it also validates our quasi-harmonic (QH) phonon model. The QH Debye–Grüneisen phonon method is combined with a DFT+OP electronic structure and used to explore the anisotropic thermal expansion in α uranium. The anomalous negative thermal expansion (contraction) of the b lattice parameter of the α-phase orthorhombic cell is relatively well reproduced from a free-energy model consisting of QH-phonon and DFT+OP electronic structure contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

Environmental controls on the kinetics of iron-sulfur cluster nucleation and nanoparticle formation

Anoxic, sulfidic conditions have been prevalent since the early Proterozoic and favor aqueous iron-sulfur (FeS aq ) clusters as a major fraction of the soluble, reduced iron and sulfur pool. FeS aq cluster formation and nucleation is driven by the high affinity between ferrous iron (Fe(II)) and sulfide (HS − ), ultimately yielding particles that precipitate as iron sulfide minerals. FeS aq clusters were recently shown to be bioavailable sources of iron and sulfur for a variety of anaerobes, yet little is known of the factors that influence the kinetics of their formation and nucleation. Here we apply computational and spectroscopic approaches to investigate the dynamics of FeS aq nucleation, cluster growth, precipitation, and redissolution as a function of Fe(II)/HS − concentration, temperature, and pH. Experiments were conducted under excess HS − to mimic euxinic conditions common to contemporary anaerobic aquatic ecosystems and those of the Proterozoic. Density functional theory calculations reveal the key role of water oxygen-iron interactions in stabilizing small FeS aq clusters and promoting solubility. Dynamic light scattering revealed a concentration-dependent increase in the kinetics of FeS aq nucleation and cluster aggregation. Increasing temperature promoted FeS aq cluster nucleation and aggregation while also enhancing dissolution. Alkaline pH also promoted FeS aq nucleation and cluster aggregation. At 25 °C, pH 7.0, and at reactant concentrations of 30 µM, FeS aq clusters < 10 nm in diameter remained in solution for > 2 h. These results underscore the importance of temperature, pH, and reactant concentration in the kinetics of FeS aq nucleation and cluster growth that, in turn, influence their bioavailability in anaerobic ecosystems.

Aquatic ecosystems↗

Gas Density Discontinuities in Merging Clusters

In a number of galaxy clusters, Chandra has discovered a new phenomenon, the sharp gas density edges. Depending on the sign of the temperature jump across the edge, these features may either be bow shocks or cold fronts. While the merger origin of the bow shocks is unambiguous, what causes cold fronts is not entirely clear, as they are observed both in mergers and in relaxed clusters. The purpose of the study of A3376, an interesting cluster with density edges, is to understand the origin of cold fronts and to look for possible shocks. This cluster also has a radio halo, and comparison of an X-ray temperature map and radio image may shed light on the nature of the cluster radio halos, which is our secondary goal. The XMM data for A3376 have been obtained, processed and cleaned. Most of the preliminary analysis, including deriving a temperature map, is complete. We also have Chandra data for this cluster, and intend to compare and combine the two datasets. We are well underway toward determining the dynamical state of A3376 and the nature of its gas density discontinuities.

Mushotzky, Richard↗

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗