Search NASA⌕ Search

SEARCH · Search NASA

Results for “MD simulation”

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 991 records · Page 55

Cyber-Physical Reconfiguration for Disaster Resilience of Power Distribution Systems

Cyber-physical distribution systems (CPDS) have emerged from the integration of information technology into distribution systems. While offering substantial benefits, this integration also introduces vulnerabilities. The interaction between cyber networks and distribution systems renders CPDS susceptible to disasters. To ensure critical load supply and system resilience, rapid post-disaster load restoration is required. The paper proposes a critical load restoration (CLR) framework in CPDS using a network reconfiguration approach that exploits the existing post-disaster resources to restore critical loads within the shortest possible time. Using graph theory, the cyber network and distribution system are integrated into a single digraph, minimizing the CLR complexity in CPDS. A cost metric is also defined to satisfy network-specific objectives and constraints. A heuristic is proposed to guide the load restoration process using the cost metric within the integrated digraph. Simulation results confirm the framework's superiority over existing literature, which either overlooks cyber components or prolongs restoration with additional resource deployment.

cyber-physical system↗

Transition from Vehicular to Structural Ionic Transport in Electrified Alkali Aqueous Solutions

A molecular understanding of the solvation and dynamics of ions under static electric fields is crucial for modelling a wide range of natural and technological processes. Yet, traditional simulation methods suffer from a trade-off that has to be made between accuracy and statistical convergence. To bridge this gap, herein we extend our recently introduced Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) approach to investigate the solvation structures and ionic transport mechanisms of electrified alkali cationic solutions. We obtain ionic conductivities for Li+, Na+ and Cs+ from the field dependence of the ionic current density in good agreement with experiment. Surprisingly, the migration mechanism is found to be strikingly different for the three ions despite their similar ionic conductivities. While Li+ conducts predominantly through vehicular migration of a stable 4-fold coordinated ion at all field strengths, Cs+ conducts strictly through a structural diffusion mechanism, where 9-12 transient first shell water coordination bonds are broken and renewed. Notably, aqueous Na+ emerges as a “Goldilocks” ion: its ion-water interactions are strong enough to maintain distinct 5-6 fold coordination shells at zero field (unlike Cs+) yet labile enough to be strongly perturbed by electric fields (unlike Li+). As a consequence, we observe an electric field-induced transition from vehicular to structural ionic transport for Na+ that is accompanied by a marked increase in the ionic current density. Our results imply that the conductance mechanism of ions with moderate ion-solvent interactions can be effectively tuned by external electric fields.

Joll, Kit↗

Design, analysis, and control of large transport aircraft utilizing engine thrust as a backup system for the primary flight controls

A review of accidents that involved the loss of hydraulic flight control systems serves as an introduction to this project. In each of the accidents--involving transport aircraft such as the DC-10, the C-5A, the L-1011, and the Boeing 747--the flight crew attempted to control the aircraft by means of thrust control. Although these incidents had tragic endings, in the absence of control power due to primary control system failure, control power generated by selective application of engine thrust has proven to be a viable alternative. NASA Dryden has demonstrated the feasibility of controlling an aircraft during level flight, approach, and landing conditions using an augmented throttles-only control system. This system has been successfully flown in the flight test simulator for the B-720 passenger transport and the F-15 air superiority fighter and in actual flight tests for the F-15 aircraft. The Douglas Aircraft Company is developing a similar system for the MD-11 aircraft. The project's ultimate goal is to provide data for the development of thrust control systems for mega-transports (600+ passengers).

Gerren, Donna S.↗

Oil in Alaska North Slope Gas Hydrate Reservoir: Micro-CT and Flow Simulation Insights into Permeability

Gas hydrate-bearing sands on the Alaska North Slope (ANS) host minor volumes of crude oil whose impact on formation permeability has never been quantified relative to a gas hydrate reservoir system. Here, we combine in situ pressure-core microcomputed-tomography (μ-CT), thermogravimetric analysis (TGA), gas-chromatography–mass-spectrometry (GC–MS), scanning-electron microscopy coupled with energy-dispersive spectroscopy (SEM-EDS), and Stokes-flow simulation to (i) quantify oil saturation, (ii) infer its source, and (iii) evaluate its influence on permeability after gas-hydrate dissociation. μ-CT and SEM-EDS imaging identify the host sediment as silty, with some grain-coating clay present. Eight representative 3003-voxel subvolumes extracted from a preserved pressure core (Core 15P-3c, from the B1 sand (unit B), HYDRATE 02 Geo Data Well (GDW), 886.96–887.29 m measured depth, MD) exhibit porosities of 0.34–0.37 and an average oil saturation of 0.09 ± 0.03, in comparison to an independent TGA-based estimate of 0.16. GC-MS fingerprints obtained from an adjacent depressurized core (Core 17P-1, 891.37–891.44 m MD), together with oil pore habit revealed by μ-CT, suggest that the oil is partly native but also includes contributions from mineral oil-based drilling fluid contamination. Further analysis of oil–porewater interactions and flushing tests confirms that native oil saturation in the analyzed samples ranges from 0.04 to 0.08 with an average of approximately 0.06, and the oil is highly immobile. Flow simulations demonstrate that the native oil saturation of ≈0.06 reduces permeability by approximately 50% due to both pore blockage and increased flow-path tortuosity. These findings reveal the potential presence of native oil within the ANS gas hydrate reservoir and suggest that future simulation models may need to account for its impact on permeability to improve long-term performance predictions of gas and water production.

02 PETROLEUM↗

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗

Protecting and Defending against Autonomous Control Systems and Digital Twin Cyber Attacks: Response Strategy for Hyperparameter attacks of Digital Twin Machine Learning Models in Nuclear Power Plants (Final)

Navigating through the complex tapestry of technological advancements, "Response Strategy for Hyperparameter attacks of Digital Twin Machine Learning Model in Nuclear Power Plants" stands at the intersection of cybersecurity and nuclear power plant operations, embarking on a journey through the intricacies of securing digital twins against malicious cyber activities. As nuclear power plants progressively integrate digital twin technology and machine learning models to optimize operations and ensure system reliability, they inadvertently expose themselves to a new spectrum of vulnerabilities, notably in the realm of hyperparameter attacks. Hyperparameters, integral in machine learning model tuning and optimal performance of digital twins, have emerged as a target for adversaries aiming to destabilize the predictive capabilities and therefore, the operational accuracy of these digital entities within critical infrastructures like nuclear plants. This paper, therefore, meticulously threads the needle through the development of a robust response strategy, poised to shield these digital reflections against calculated hyperparameter manipulations, ensuring that the digital twin can effectively and securely function as a reliable proxy for its physical counterpart. The ensuing sections delve into the orchestrated maelstrom of multi-rate time-changing intelligent coordinated hyperparameter attacks and the implementation of event-triggered predictive control, laying down a structured, predictive, and responsive framework that safeguards the nexus where the digital and physical realms of nuclear power plants coalesce. The operational integrity of digital twins in nuclear power plants depends critically on the security of machine learning hyperparameters. This study makes two different contributions. First, a decision-based idea known as a multi-rate time changing intelligent coordinated hyperparameter attack is put forth. In this attack, many hyperparameters are repeatedly changed using both random and intelligent optimal techniques by the attacker. These assaults introduce varied rates at different attack steps, compromise various amounts of hyperparameters, and improve stealth and flexibility. Second, a technique is developed for event triggered predictive control to rapidly respond to potential hyperparameter attacks. This control integrates a sliding window framework, retaining a history of previous data points and employing linear regression to predict the next data point from the current dataset. The control gain K is determined using the Lyapunov-Krasovskii method, and subsequently, an action is developed. Finally, the outcome of the simulation demonstrates the viability of the proposed method for defending nuclear power plant digital twins from hyperparameter attacks.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Tough Errors are no Match (TEAM): Optimizing the Quantum Compiler for Noise Resilience

This project builds toward a comprehensive error-mitigating toolkit that makes quantum programming more robust and adaptive to the noisy, resource-limited nature of today’s quantum hardware. To that end, it integrates established error-mitigation methods — such as zero-noise extrapolation and dynamical decoupling — directly into compiler infrastructures. These techniques will be packaged as modules that can automatically adjust and combine based on performance analysis, enabling compilers to explore large design spaces and produce optimized, low-noise quantum programs with minimal manual intervention. In parallel, this project also explores new approaches to analog quantum programming or quantum simulation, and has developed the programming language SimuQ which treats quantum Hamiltonian evolution as the central object.

97 MATHEMATICS AND COMPUTING↗

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗

Pathfinding quantum simulations of neutrinoless double- β decay

We present results from co-designed quantum simulations of the neutrinoless double- β decay of a simple nucleus in 1+1D quantum chromodynamics using IonQ’s Forte-generation trapped-ion quantum computers. Electrons, neutrinos, and up and down quarks are distributed across two lattice sites and mapped to 32 qubits, with an additional 4 qubits used for flag-based error mitigation. A four-fermion interaction is used to implement weak interactions, and lepton-number violation is induced by a neutrino Majorana mass. Quantum circuits that prepare the initial nucleus and time evolve with the Hamiltonian containing the strong and weak interactions are executed on IonQ Forte Enterprise. Enabled by tuned model parameters, lepton-number violation is observed in real time, providing a clear signal of neutrinoless double- β decay. This was made possible by co-designing the simulation to maximally utilize the all-to-all connectivity and native gate-set available on IonQ’s quantum computers. Quantum circuit compilation techniques and co-designed error-mitigation methods, informed from executing benchmarking circuits with up to 2,356 two-qubit gates, enabled observables to be extracted with high precision. We discuss the potential of future quantum simulations to provide yocto-second resolution of the reaction pathways in these, and other, nuclear processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An International Round-Robin Study on Thermoelectric Module Testing and Development of Standard Power Generation Modules

An international round-robin study on thermoelectric power generation modules was conducted with nine participating laboratories. Two types of commercially available bismuth telluride modules, 30 mm × 30 mm and 40 mm × 40 mm, were used. A test protocol was followed with five temperature set points from 50°C to 150°C. Graphite sheets were used as thermal interface materials with test pressure at 100 psi (0.69 MPa). The results showed large lab-to-lab variations and the key source of uncertainty for module efficiency was identified as the heat flux measurement. In the meantime, significant uncertainty was also found in maximum electrical power (P max ) measurements. As a result of the round-robin, a “standard module” with 4 × 4 legs on a 20 mm × 20 mm platform was suggested. A skutterudite module and a half-Heusler module were produced with identical geometry and 4 mm × 4 mm × 8 mm legs. All transport properties to calculate the figure-of-merit, zT, were measured from ambient temperature to 500°C. Module performance was measured by two laboratories. Two finite-element-analysis (FEA)-based models were developed independently to simulate and predict the module performance. In conclusion, the standard modules eliminated significant test uncertainties and are aimed at assisting device design and achieving more accurate performance predictions.

Round-robin↗

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗

Interpreting the Strain Tensor During Injection into a Reservoir

A suite of well tests was conducted to evaluate the feasibility of measuring and interpreting the strain tensor during injection. Water was injected at 1 liter/s for a few days to a few weeks into a well completed at a depth of 530 m in the Bartlesville sandstone in the North Avant Field in Oklahoma. The strain tensor was measured using strainmeters at a depth of 30m and a distance of 220m from the well. We used a Gladwin strainmeter, a state-of-the-art geodetic strainmeter, as well as novel borehole strainmeters developed for the project. After processing, the horizontal strains were tensile and the vertical strain was compressive during injection. The radial strain was less than the circumferential strain. Strain magnitudes increased during injection and recovered afterward. The strain rate was as fast as 100 nanostrain/d at the start of injection, and decreased to 10 nanostrain/d after a few days. The normal strain signals were approximately proportional to the log of time during the first few days of injection, and the semi-log slope was 50 to 75 nanostrain. Data from the Gladwin strainmeter were essentially identical to data from a strainmeter we developed using optical fiber sensors. Tilt was measured using two tiltmeters, but the tilts were erratic compared to the normal strain and were of limited value. Strain data were analyzed using a strain type-curve method, an analytical solution to a poroelastic inhomogeneity, and a stochastic inversion method. Poroelastic simulations explain the data well. Results from the interpretation methods are consistent with each other, and with independent information characterizing the area. The hydraulic diffusivity is approximately 0.5 m2/s and the permeability is between 100 and 500 mD, according to analysis of the strain signals. This is consistent with measurements of permeability from cores, and with analyses of pressure transients. The analyses indicate that Young’s modulus of the reservoir is between 2 and 6 GPa. Strain interpretations indicate the well intersects a permeable lens with a maximum dimension of 500 to 1,000 m. An isopach map of coarse-grained sand confirms a feature roughly 1,000 m in lateral extent. The results demonstrate that the strain signal measured at shallow (30m) depth can be used to infer properties and heterogeneities in an underlying reservoir or aquifer.

Murdoch, Larry↗

Analysis of interior noise ground and flight test data for advanced turboprop aircraft applications

Interior noise ground tests conducted on a DC-9 aircraft test section are described. The objectives were to study ground test and analysis techniques for evaluating the effectiveness of interior noise control treatments for advanced turboprop aircraft, and to study the sensitivity of the ground test results to changes in various test conditions. Noise and vibration measurements were conducted under simulated advanced turboprop excitation, for two interior noise control treatment configurations. These ground measurement results were compared with results of earlier UHB (Ultra High Bypass) Demonstrator flight tests with comparable interior treatment configurations. The Demonstrator is an MD-80 test aircraft with the left JT8D engine replaced with a prototype UHB advanced turboprop engine.

Simpson, M. A.↗

Comparison ofdvanced turboprop interior noise control ground and flight test data

Interior noise ground tests conducted on a DC-9 aircraft test section are described. The objectives were to study ground test and analysis techniques for evaluating the effectiveness of interior noise control treatments for advanced turboprop aircraft, and to study the sensitivity of the ground test results to changes in various test conditions. Noise and vibration measurements were conducted under simulated advanced turboprop excitation, for two interior noise control treatment configurations. These ground measurement results were compared with results of earlier UHB (Ultra High Bypass) Demonstrator flight sts with comparable interior treatment configurations. The Demonstrator is an MD-80 test aircraft with the left JT8D engine replaced with a prototype UHB advanced turboprop engine.

Simpson, Myles A.↗

Kinetics and Mechanisms of ZnO to ZIF-8 Transformations in Supercritical CO2 Revealed by in situ X-ray Diffraction

ZIF-8 was synthesized in supercritical carbon dioxide (scCO2). In situ powder X-ray diffraction, ex situ microscopy and simulations, provides a encompassing view of the formation of ZIF-8 and intermediary ZnO@ZIF-8 composites in this nontraditional solvent. Time-resolved imaging exposed divergent physicochemical reaction pathways from previous studies of the growth of anisotropic ZIF-8 core@shell structures in traditional solvents. Synthetically relevant physio-chemical properties of scCO2 were integrated into classical nucleation theory, relating interfacial forces with 3D nu-cleation outcomes. The kinetics of crystallization were examined and displayed a characteristic signature of time- and temperature-dependent mechanisms over the extent of the reaction. Lastly, we show that subtle factors, such as the extent of reaction and the size/shape of sacrificial templates can tailor ZIF-8 composition and size, eliciting control over hierarchical porosity in a nonconventional green solvent.

Sinnwell, Michael A.↗

Electronic Structure and Bonding of US, SUO, and US 2

Anion photoelectron spectra of US – and US 2 – were recorded using the third (355 nm) and fourth (266 nm) harmonics of an Nd:YAG laser, which yielded vertical detachment energies (VDEs) of 1.71 and 2.02 eV, respectively. The experimental results are supported by extensive relativistic ab initio calculations, primarily at the coupled cluster level of theory, with systematic sequences of correlation consistent basis sets. Calculations include the closely related SUO and SUO – molecules, as well as the oxide congeners UO/UO – and UO 2 /UO 2 – which are well-known experimentally and provide benchmark systems for the sulfide calculations. Adiabatic electron detachment energies (ADEs) are computed for UO – , UO 2 – , US – , SUO – , and US 2 – using the Feller–Peterson–Dixon (FPD) composite approach. Additionally, ADEs are determined for UO – and US – using a spinor-based coupled cluster approach where spin–orbit coupling is included at the orbital level. VDEs are derived from the ab initio results from Franck–Condon simulations of the photoelectron spectra.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Influence of Surface Aerosol Injection on Stratocumulus‐to‐Cumulus Transition: Cloud‐Surface Coupling and Background Aerosol Concentrations

The influence of surface aerosol injection on the stratocumulus‐to‐cumulus transition (SCT) is explored using large‐eddy simulations. We examine how cloud‐surface coupling (or the strength of the marine boundary layer (MBL) stratification that limits vertical turbulent mixing and convection) impacts the vertical transport of aerosols, and how injected aerosols influence cloud properties and associated cloud radiative effects during the SCT. By injecting aerosols at different stages of the SCT, noting that cloud is more decoupled from the surface over time due to entrainment warming, we find that cloud‐surface coupling significantly affects aerosol vertical transport. However, injection timing (before drizzle if any) does not notably affect the SCT and the efficiency of marine cloud brightening, because aerosol number concentrations due to injections at different times rapidly converge before the transition onset. By varying the background aerosol concentration, we find that injected aerosols can significantly extend the persistence of stratocumulus decks by suppressing precipitation in clean environments but have little impact on stratocumulus breakup with higher background aerosol concentrations due to saturated aerosol effects. In clean MBLs, the SCT‐delay‐induced increase in cloud fraction dominates the overall cooling effects in response to aerosols, followed by Twomey effects. These cooling effects are slightly offset by decreased liquid water path (LWP) due to entrainment drying. In polluted MBLs, the Twomey effect is more dominant, followed by cloud fraction adjustments, and these coolings are also partly offset by LWP adjustments. All the simulations are made in relatively small domains in which injected aerosols are homogenized over a short time scale.

Zhang, Haipeng [Univ. of Maryland, College Park, M↗

Omnigenous stellarators with improved ideal and kinetic ballooning stability

Omnigenity is a property of a magnetic field which ensures confinement of trapped particles. It is a necessary requirement for any high-performance stellarator. After creating an omnigenous equilibrium, one must also ensure reduced transport resulting from kinetic and magnetohydrodynamic (MHD) instabilities. To this end, we leverage the GPU-accelerated DESC optimization suite, which is used to design stable, finite-β omnigenous equilibria with poloidal, toroidal, and helical symmetry, achieving Mercier, ideal ballooning, and as a consequence, improved kinetic ballooning stability. We discover stellarators with second stability, a regime of large pressure gradient where an equilibrium becomes ideal ballooning stable, and demonstrate and explore both using theory and gyrokinetic simulations the connection between ideal and kinetic ballooning stability.

optimization↗