Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamics”

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 433 records · Page 24

Dynamics of plasma atomic layer etching: Molecular dynamics simulations and optical emission spectroscopy

Atomic layer etching is intrinsically dynamic as it involves sequential and repeated exposures of a surface to be etched with different species at different energies. The composition and structure of the near surface region change in both time and depth. Full understanding of this process requires resolving both temporal and spatial variations. In this work, we consider silicon (Si) atomic layer etching (ALE) by alternating exposure to chlorine gas (Cl 2 ) and argon ions (Ar + ). Molecular dynamics (MD) simulations are compared to experimental measurements with the aim of better understanding the dynamics of ALE and to test the simulation procedure. The simulations help to more fully interpret the experimental measurements. Further, optical emission measured just above the surface being etched can be related to etch products and can, therefore, be directly compared to simulation predictions. The simulations capture the measured initial product distribution leaving the surface and match the measured etch per cycle reasonably well. While simulations demonstrate the importance of ion-induced surface damage and mixing into a layer below the surface, the depth of which depends mainly on ion energy, the experiments suggest there is more Cl mixed into the layer than the MD procedure predicts.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dynamical Aquaplanet Experiments with Uniform Thermal Forcing: System Dynamics and Implications for Tropical Cyclone Genesis and Size

Existing hypotheses for the dynamical dependence of tropical cyclone genesis and size on latitude depend principally on the Coriolis parameter $f$. These hypotheses are tested via dynamical aquaplanet experiments with uniform thermal forcing in which planetary rotation rate and planetary radius are varied relative to Earth values; the control simulation is also compared to a present-day Earth simulation. Storm genesis rate collapses to a quasi-universal dependence on f that attains its maximum at the critical latitude, where the inverse-$f$ scale and Rhines scale are equal. Minimum genesis distance from the equator is set by the equatorial Rhines (or deformation) scale and not by a minimum value of $f$. Outer storm size qualitatively follows the smaller of the two length scales, including a slow increase with latitude equatorward of 45° in the control simulation, similar to the Earth simulation. The latitude of peak size scales with the critical latitude for varying planetary radius but not rotation rate, possibly because of the dependence of genesis specifically on $f$. The latitudes of peak size and peak packing density scale closely together. Results suggest that temporal effects and interstorm interaction may be significant for size dynamics. More generally, the critical latitude separates two regimes: 1) a mixed wave–cyclone equatorial belt, where wave effects are strong and the Rhines scale likely limits storm size, and 2) a cyclone-filled polar cap, where wave effects are weak and cyclones persist. Finally, the large-planet limit predicts a world nearly covered with long-lived storms, equivalent to the f plane. Overall, spherical geometry is likely important for understanding tropical cyclone genesis and size on Earthlike planets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data Projection of the High Temperature Electrolysis System in the Dynamic Energy Transport and Integration Laboratory using Dynamic System Scaling

For nuclear power to be flexible in a functioning Integrated Energy System (IES), excess produced heat must be stored or utilized during times of low power demand to ensure a load factor of 1 while load balancing. The Dynamic Energy Transport and Integration Laboratory (DETAIL) is one facility that is under development to emulate IES conditions on the engineering-scale, planned to conduct virtual real time operations with industry-scale facilities, and is currently testing thermal storage and high temperature electrolysis. As part of the study to develop a method to preprocess input signals or postprocess output signals between systems of different scales via Dynamical System Scaling (DSS), the current research is one of the continued efforts branching from the data projection activity conducted for the Thermal Energy Distribution System and currently engages the High Temperature Electrolysis (HTE) System in DETAIL. The HTE SOEC electrical, fluid, and thermal dynamics Figure of Merits (FOM) were identified, governing equations and closure relations were successfully scaled, and relations between FOM scaling ratios were determined. Setting the scaling objectives to reform existing data to project a data set that doubly accelerated the electrolysis process while preserving the produced amount of hydrogen was generated for the full transient. The calculated boundary conditions were inlet temperature, stack current, and inlet steam mass flow rate at 1470 K, 121.1 A, and 1.886 g/s, respectively. The research outcomes demonstrated an output signal postprocessing case accelerating the hydrogen production without changing geometry, number of cells, and partial pressures.

08 HYDROGEN↗

Computational Fluid Dynamics Simulations of the Transonic Dynamics Tunnel Airstream Oscillator System

This paper presents a Computational Fluid Dynamics (CFD) model of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT) with the Airstream Oscillator System (AOS) in operation. The TDT is a continuous-flow, closed circuit wind tunnel with a 16- by 16-foot slotted test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas (R-134a) as the test medium. The TDT was also fitted with an AOS which can be used to create a gust field in the test section. To date, no computational analyses of the TDT involving the AOS have been performed. In this study, experimental data acquired of an oscillating airstream in the tunnel will be used to calibrate the computational analyses. A significant motivation of this work is to attempt to compare computational gust velocities to those recorded in the TDT literature. In addition to a validation of this model with experimental data, it may also be possible to supplement the experimental data with computational data. The experimental data is rather sparse and at only a few Mach numbers. Computational results may be able to expand the AOS data set. Another motivation for this work is that the Integrated Adaptive Wing Technology Maturation (IAWTM) semi-span model of the high-aspect-ratio CRM with 10 trailing edge control surfaces is currently being fabricated and will be delivered to the NASA Langley Transonic Dynamics Tunnel (TDT) for testing in late 2020. Among tests to be conducted, gust load alleviation (GLA) will be demonstrated at transonic conditions using the AOS.

Computational Fluid Dynamics↗

Testing thermal gradient driving force for grain boundary migration using molecular dynamics simulations

Strong thermal gradients in low-thermal-conductivity ceramics may drive extended defects, such as grain boundaries and voids, to migrate in preferential directions. In this work, molecular dynamics simulations are conducted to study thermal gradient driven grain boundary migration and to verify a previously proposed thermal gradient driving force equation, using uranium dioxide as a model system. It is found that a thermal gradient drives grain boundaries to migrate up the gradient and the migration velocity increases under a constant gradient owing to the increase in mobility with temperature. Different grain boundaries migrate at very different rates due to their different intrinsic mobilities. The extracted mobilities from the thermal gradient driven simulations are compared with those calculated from two other well-established methods and good agreement between the three different methods is found, demonstrating that the theoretical equation of the thermal gradient driving force is valid, although a correction of one input parameter should be made. The discrepancy in the grain boundary mobilities between modeling and experiments is also discussed.

36 MATERIALS SCIENCE↗

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↗

Differentiation of hepatic abscess from metastasis on contrast-enhanced dynamic computed tomography in patients with a history of extrahepatic malignancy: emphasis on dynamic change of arterial rim enhancement

Objectives: The objective of the study is to identify computed tomography (CT) findings that differentiate hepatic abscess from hepatic metastasis in a patient with a history of extrahepatic malignancy. Materials and methods: This retrospective study included 30 patients with 93 hepatic abscesses and 40 patients with 125 hepatic metastases who had a history of extrahepatic malignancy and underwent contrast-enhanced dynamic CT with arterial phase (AP) and portal venous phase (PVP). The diagnosis of hepatic abscess and hepatic metastasis was made using pathological confirmation or clinical diagnosis. Margin, patchy parenchymal enhancement, arterial rim enhancement, dynamic change of arterial rim enhancement, size discrepancy of lesions between arterial and portal phases, bile duct dilatation, perilesional hyperemia, and perilesional low density were evaluated by two radiologists independently. Significant findings for differentiating two groups were identified at univariate and multivariate analysis with nomogram for predicting hepatic abscess. Interobserver agreement was also analyzed for each variable. Results: Multivariate analysis revealed that patchy parenchymal enhancement (P < 0.001), arterial rim enhancement persistent through PVP (P < 0.001), and perilesional hyperemia (P = 0.013) were independent significant findings to predict hepatic abscess than metastasis. Among them, arterial rim enhancement persistent through PVP showed a highest odds ratio (OR 33.73) on multivariate analysis and a highest predictor point on a nomogram for predicting hepatic abscess. When two of these three criteria were combined, 80.7% (75/93) of hepatic abscess were correctly identified, with a specificity of 85.6% (107/125). When all three criteria were satisfied, specificity was up to 100% (125/125). Conclusions: At contrast-enhanced dynamic CT, patchy parenchymal enhancement, arterial rim enhancement persistent through PVP, perilesional hyperemia, and their combinations may be reliable CT features for differentiating hepatic abscess from metastasis in patients with a history of primary extrahepatic malignancy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Observational constraints on the tilted flat-XCDM and the untilted nonflat XCDM dynamical dark energy inflation parameterizations

In this work, we constrain tilted spatially-flat and untilted nonflat XCDM dynamical dark energy inflation parameterizations using Planck 2015 cosmic microwave background (CMB) anisotropy data and recent baryonic acoustic oscillations distance measurements, Type Ia supernovae data, Hubble parameter observations, and growth rate measurements. Inclusion of the four non-CMB data sets results in a significant strengthening of the evidence for nonflatness in the nonflat XCDM model from 1.1σ for the CMB data alone to 3.4σ for the full data combination. In this untilted nonflat XCDM case the data favor a spatially-closed model in which spatial curvature contributes a little less than a percent of the current cosmological energy budget; they also mildly favor dynamical dark energy over a cosmological constant at 1.2σ. These data are also better fit by the flat-XCDM parameterization than by the standard ΛCDM model, but only at 0.3σ significance. Current data is unable to rule out dark energy dynamics. The nonflat XCDM parameterization is compatible with the Dark Energy Survey limits on the present value of the rms mass fluctuations amplitude (σ 8 ) as a function of the present value of the nonrelativistic matter density parameter (Ω m ). Furthermore, it does not provide as good a fit to the higher multipole CMB temperature anisotropy data as does the standard tilted flat-ΛCDM model. A number of measured cosmological parameter values differ significantly when determined using the tilted flat-XCDM and the nonflat XCDM parameterizations, including the baryonic matter density parameter and the reionization optical depth.

79 ASTRONOMY AND ASTROPHYSICS↗

Structure and dynamics of a compact state of a multidomain protein, the mercuric ion reductase

Here, the functional efficacy of colocalized, linked protein domains is dependent on linker flexibility and system compaction. However, the detailed characterization of these properties in aqueous solution presents an enduring challenge. Here, we employ a novel, to our knowledge, combination of complementary techniques, including small-angle neutron scattering, neutron spin-echo spectroscopy, and all-atom molecular dynamics and coarse-grained simulation, to identify and characterize in detail the structure and dynamics of a compact form of mercuric ion reductase (MerA), an enzyme central to bacterial mercury resistance. MerA possesses metallochaperone-like N-terminal domains (NmerA) tethered to its catalytic core domain by linkers. The NmerA domains are found to interact principally through electrostatic interactions with the core, leashed by the linkers so as to subdiffuse on the surface over an area close to the core C-terminal Hg(II)-binding cysteines. How this compact, dynamical arrangement may facilitate delivery of Hg(II) from NmerA to the core domain is discussed.

59 BASIC BIOLOGICAL SCIENCES↗

Elucidating the Solvation Structure and Dynamics of Lithium Polysulfides Resulting from Competitive Salt and Solvent Interactions

Designing optimal electrolytes is key to enhancing the performance of energy storage devices,especially relating to cycle life, efficiency, and stability.1 Specifically, for future beyond-Li ion energy storage, there is still ample room for electrolyte improvements. Among the candidates for higher gravimetric energy storage, the Li-S battery is considered quite promising, owing to its theoretical specific energy density (2600 Wh/kg) and specific capacity (1675 mAh/g) and significantly lower cost as compared to state-of-art lithium-ion batteries.

36 MATERIALS SCIENCE↗

Dynamic Process Intensification via Data-Driven Dynamic Optimization: Concept and Application to Ternary Distillation

Process intensification is a design philosophy aimed at making chemical processes safer and more efficient. Its implementation often results in significant modifications to the design and structure of the process, with several conventional unit operations occurring in the same physical device. Traditionally, process intensification has focused on steady-state operation. In our previous works, we introduced dynamic process intensification (DPI) as a new intensification paradigm based on operational changes for conventional or intensified units. DPI is predicated on switching operation between two auxiliary steady states selected via a steady-state optimization calculation that ensures that the system generates, on average and over time, the same products as in nominal steady-state operation, but with favorable economics. This paper extends the DPI concept and introduces a novel dynamic optimization-based DPI strategy (Do-DPI) that involves imposing a true cyclic operation rather than switching between two discrete states. We discuss its implementation using surrogate dynamic models learned via system identification. Here, an extensive case study concerning a ternary distillation column separating a canonical hydrocarbon mixture shows that Do-DPI can reduce energy use by more than 4% relative to steady-state operation, with no significant deviations in product quality and production rate.

42 ENGINEERING↗

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17 O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and 17 O NMR quadrupolar relaxation rates are investigated in the series of complexes UO 2 2+ , UO 2 (CO 3 ) 3 4– , and (UO 2 ) 3 (CO 3 ) 6 6– . Car–Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO 2 2+ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl 17 O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO 2 (CO 3 ) 3 4– and (UO 2 ) 3 (CO 3 ) 6 6– complexes show a solvent dependence of the terminal carbonate 17 O and 13 C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin–orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique 13 C chemical shifts is correctly reproduced for (UO 2 ) 3 (CO 3 ) 6 6– . NMR relaxation rate data support large 17 O peak widths, but remain below those noted in the experimental literature. Furthermore, comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonadiabatic Dynamics in a Continuous Circularly Polarized Laser Field with Floquet Phase-Space Surface Hopping

Non-adiabatic chemical reactions involving continuous circularly polarized light (cw CPL) have not attracted as much attention as dynamics in unpolarized/linearly polarized light. However, including circularly (in contrast to linearly) polarized light allows one to effectively introduce a complex-valued time-dependent Hamiltonian, which offers a new path for control or exploration through the introduction of Berry forces. Here, we investigate several inexpensive semiclassical approaches for modeling such nonadiabatic dynamics in the presence of a time-dependent complex-valued Hamiltonian, beginning with a straightforward instantaneous adiabatic fewest-switches surface hopping (IA-FSSH) approach (where the electronic states depend on position and time), continuing to a standard Floquet fewest switches surface hopping (F-FSSH) approach (where the electronic states depend on position and frequency), and ending with an exotic Floquet phase-space surface hopping (F-PSSH) approach (wherethe electronic states depend on position, frequency, and momentum). Using a set of model systems with time-dependent complex-valued Hamiltonians, we show that the Floquet phase- space adiabats are the optimal choice of basis as far as accounting for Berry phase effects and delivering accuracy. Furthermore, the F-PSSH algorithm sets the stage for modeling nonadiabatic dynamics under strong externally pumped circular polarization in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗