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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Combined Day-ahead and Real-time approach for Scheduling Real and Reactive power dispatch of Battery Energy Storage System

A combined day-ahead (look-ahead) and real-time optimization approach is proposed to schedule BESS real and reactive power dispatch. The look ahead scheduling agent maximizes profit and essentially performs the function of energy arbitrage by looking at forecasted prices, network constraints and BESS operating conditions to generate real power dispatch schedule. While this schedule is followed, the real time optimization agent computes an optimum reactive power set-point (for a particular objective function) using optimal power flow framework. The concepts are implemented by co-simulating legacy tools and demonstrated on a large real world distribution system. In addition to the profits from energy arbitrage, optimum reactive power set point computation results in achieving desirable network attributes like flat voltage profiles and minimizing voltage deviations from desired value.

Bhatti, Bilal A.↗
NASA NTRS2022

P- and L-Band Retrieval of Subsurface Soil Moisture and Temperature Profiles as First-Order Polynomial Function

This paper demonstrates the potential use of P and L band passive measurements to determine root zone soil moisture (SM) and soil temperature(ST). SM and ST data have been taken as a function of depth during the NASA GSFC PLEX19 experiment in the summer of 2019 at Beltsville, MD, USA. Using these data, a coherent model has been used to compute H and V brightness temperatures at frequencies of 0.8 and 1.4 GHz with an observation angle of 35 degrees. These synthetic brightness data are then used to estimate the SM and ST profiles which are represented by linear polynomials. The inversion problem is formulated as a least square problem that is solved by a global optimization method known as the Adaptive Simulated Annealing(ASA) method. Four inversion examples having different SM and ST profiles are presented. Selected results show that the standard deviation between the retrieved and measured data is less than 0.077 cm3/cm3 for SM, and 2.245 °C for ST.

Ming Li↗
DOE OSTI2023

A Membrane Contactor Enabling Energy-Efficient CO 2 Capture from Point Sources with Deep Eutectic Solvents

We demonstrate a scalable and energy-efficient hollow fiber membrane contactor (HFMC)-based process using a green solvent for CO 2 capture. This process uses a deep eutectic solvent (DES) in an HFMC to provide close interfacial interactions and contact between the DES and CO 2 . This approach overcomes disadvantages associated with direct absorption in DES and could potentially be applied to a variety of solvent-based CO 2 capture methods. Commercial low-cost polymer hollow fiber membranes (e.g., microporous polypropylene) were evaluated for CO 2 capture with reline, a prototypical DES. Single-gas measurements showed that the DES-based polypropylene HFMC can capture and separate CO 2 while rejecting N2. From a mixed gas containing 50 mol % N 2 and 50 mol % CO ν , the DES-based HFMC separated CO 2 with a purity of 96.9 mol %. The effect of several process parameters including solvent flow rate, pressure, and temperature on the CO 2 separation performance was studied. The flux of the recovered CO 2 was 67.43 mmole/m 2 /h at a feed pressure of 4 bar. In situ Fourier transform infrared (FTIR) measurements combined with density functional theory (DFT)-based molecular dynamics simulations revealed that reline absorbs CO 2 by physical absorption without forming a new chemical compound, and CO 2 separation by reline occurs via the pressure swing mechanism. This research provides fundamental insights about physical solvent-based separation processes and a pathway toward practical deployment.

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

KRAS4a and KRAS4b show distinct lipid-dependent regulation of RAS-RAF membrane dynamics

KRAS4a and KRAS4b are important regulators of signaling, and their interactions with the plasma membrane are dynamic and influenced by lipid composition. KRAS 4a and 4b have nearly identical globular domains but differ in their membrane-associated hyper variable region (HVR). The functional distinctions between these isoforms remain unclear, particularly with regards to their dependence on specific lipids and the membrane environment. Previous work showed that the membrane orientation of KRAS4b affects its ability to bind to RAF kinase RBDCRD and that the KRAS–RBDCRD complex adopts different poses on the membrane as well as influences the size and composition of the lipid environment. To model differences between KRAS 4a and 4b protein–lipid interactions, we extended the Multiscale Machine-Learned Modeling Infrastructure (MuMMI) to incorporate continuum simulations in the grand canonical ensemble, enabling sampling across macroscopic, coarse-grained, and all-atom resolutions. Using this framework, we systematically altered PIP2 concentrations, KRAS 4a versus 4b, and RAF RBDCRD complexation to assess impacts on membrane–protein interactions and dynamics. Our results reveal that reducing PIP2 shifts and broadens the membrane orientational preference of both KRAS 4b and 4a, with stronger effects on 4b HVR localization versus 4a. We demonstrate that with depletion of the strong negatively charged PIP2 lipid, the less charged phosphatidylserine replaces PIP2. Our findings highlight similarities and distinctions in the dynamics and lipid dependency of KRAS isoforms and suggest that ordering of the local lipid composition by HVRs is a shared property and key modulator of RAS-mediated signaling at the plasma membrane.

Biological and medical sciences↗
DOE OSTI2025

Interface‐Induced Stability of Nontrivial Topological Spin Textures: Unveiling Room‐Temperature Hopfions and Skyrmions

Topological spin configurations, such as soliton-like spin texture and Dirac electron assemblies, have recently emerged in fundamental science and technology. Achieving stable topological spin textures at room temperature is crucial for their use as long-range information carriers. However, their creation and manipulation are hindered by multi-step field training and competing interactions. Thus, a spontaneous ground state for multidimensional topological spin textures is desirable, with skyrmions forming swirling, hedgehog-like spin structures in two dimensions and hopfions as their twisted 3D counterparts. Here, the first observation of robust and reproducible topological spin textures of hopfions and skyrmions observed at room temperature and in zero magnetic field is reported, which are stabilized by geometric confinement and protected by interfacial magnetism in a ferromagnet/topological insulator/ferromagnet trilayer heterostructure. These skyrmion-hopfion configurations are directly observed at room temperature with Lorenz transmission electron microscopy. Using micromagnetic modeling, the experimental observations of hopfion-skyrmion assemblies are reproduced. This model reveals a complete picture of how spontaneously organized skyrmion lattices encircled by hopfion rings are controlled by surface electrons, uniaxial anisotropy, and Dzyaloshinskii-Moriya interaction. This study provides evidence that topological chiral spin textures can facilitate the development of magnetic topological carriers, paving the way for ultralow-power and high-density information processing.

hopfions↗
DOE OSTI2024

Real-time hybrid controls of energy storage and load shedding for integrated power and energy systems of ships

This paper presents an original energy management methodology to enhance the resilience of ship power systems. The integration of various energy storage systems (ESS), including battery energy storage systems (BESS) and super-capacitor energy storage systems (SCESS), in modern ship power systems poses challenges in designing an efficient energy management system (EMS). The EMS proposed in this paper aims to achieve multiple objectives. The primary objective is to minimize shed loads, while the secondary objective is to effectively manage different types of ESS. Considering the diverse ramp-rate characteristics of generators, SCESS, and BESS, the proposed EMS exploits these differences to determine an optimal long-term schedule for minimizing shed loads. Furthermore, the proposed EMS balances the state-of-charge (SoC) of ESS and prioritizes the SCESS’s SoC levels to ensure the efficient operation of BESS and SCESS. For better computational efficiency, we introduce the receding horizon optimization method, enabling real-time EMS implementation. Further, a comparison with the fixed horizon optimization (FHO) validates its effectiveness. Simulation studies and results demonstrate that the proposed EMS efficiently manages generators, BESS, and SCESS, ensuring system resilience under generation shortages. Additionally, the proposed methodology significantly reduces the computational burden compared to the FHO technique while maintaining acceptable resilience performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Impacts of swinging door lossy compression of synchrophasor data

As the number of operational Phasor Measurement Units (PMUs) is increasing, compression is becoming crucial for efficiently transferring and storing enormous amount of synchrophasor data continuously measured from a power system. Some utilities in the U.S. have already adopted a nonlinear compression technique, called Swinging Door Algorithm (SDA), assuming compressed data meets technical and regulatory requirements without resulting in any degradation of the measurements. This paper studies impacts of this compression algorithm on field measured and simulated PMU data, and presents results from spectral analyses of raw and decompressed measurements. The results show blind application of the SDA compression alters frequency content in the measurements and introduces false harmonics and aliasing. These can lead to false conclusions being drawn from the measurements about the power system condition.

42 ENGINEERING↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2023

Change of Electrical and Transport Properties of Nickel Oxide by Carrier Concentration and Temperature through First-Principle Calculations

Abstract Nickel is typically used as one of the main components in electrical contact devices or connectors. Nickel oxide (NiO) is usually formed on the surfaces of electrodes and can negatively impact system performance by introducing electrical contact resistance. The thermal, electrical, and transport properties of NiO, as a Mott insulator or a p-type semiconductor, can be altered by operating and environmental conditions such as temperature and stress/strain by contact. In this study, we investigate the fundamental material properties of NiO through the first-principle calculations. First, we obtain and compare the lattice parameter, magnetic moment, and electronic structure for NiO via the WIEN2K simulations with four different potentials (i.e., GGA, GGA + U, LSDA, and LSDA + U). Then, using the WIEN2K simulation results with LSDA + U potential that produces a highly accurate bandgap for NiO, we calculate the electrical conductivity and electrical part of the thermal conductivity of nickel and NiO as a function of temperature and carrier concentration through the BoltzTraP simulations. Systematic simulation results revealed that the electrical conductivity relative to the relaxation time for NiO increases with the carrier concentration, while it shows a slightly decreasing trend with temperature under a fixed carrier concentration. By contrast, the electrical part of the thermal conductivity shows an increasing trend considering carrier concentration and temperature.

Shah, Nayem Md. Reza↗
DOE OSTI2020

Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of V o complex

Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a Vo proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo–electron microscopy (cryo-EM) structure of yeast V o proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of V o that shows the c-ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c-ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free V o that is generated as a result of V-ATPase regulation by reversible disassembly in vivo.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Low viscosity of solid MgO at high pressures and strain rates measured using the laser-driven Richtmyer-Meshkov instability

Solids are often assumed to behave as viscous fluids under high-strain rates. This behavior has been studied experimentally in metals but largely unexplored in brittle ceramic materials. Here, in this study, we present a technique for measuring the viscosity of MgO using time-resolved velocimetry to track the growth rate of a perturbation caused by the Richtmyer-Meshkov instability at the OMEGA EP laser facility. To interpret the results, we use an in-house Eulerian hydrocode to simulate our experiments and model the plastic deformation of solid MgO as a viscous fluid. Results indicate that MgO has a surprisingly low upper bound to its effective viscosity of ∼10 2 ⁢ Pa s at 175 ± 15⁢ GPa, ∼3500⁢K, and 10 6 −10 7 s −1 strain rate.

Perez, Tyler M. [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2025

New particle pusher with hadronic interactions for modeling multimessenger emission from compact objects

We propose novel numerical schemes based on the Boris method in curved spacetime, incorporating both hadronic and radiative interactions for the first time. Once the proton has lost significant energy due to radiative and hadronic losses, and its gyroradius has decreased below typical scales on which the electromagnetic field varies, we apply a guiding center approximation (GCA). We fundamentally simulate collision processes either with a Monte-Carlo method or, where applicable, as a continuous energy loss, contingent on the local optical depth. To test our algorithm for the first time combining the effects of electromagnetic, gravitational, and radiation fields including hadronic interactions, we simulate highly relativistic protons traveling through various electromagnetic fields and proton backgrounds. We provide unit tests in various spatially dependent electromagnetic and gravitational fields and background photon and proton distributions, comparing the trajectory against analytic results. We propose that our method can be used to analyze hadronic interactions in black hole accretion disks, jets, and coronae to study the neutrino abundance from active galactic nuclei.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Maximizing dynamic range and performance of anatase TiO 2 ECRAM through structure and programming

Here, in this study, we investigate the structure-dependent modulation characteristics of all-solid-state three-terminal electrochemical random-access memory (ECRAM) based on an anatase Li x TiO 2 channel. By directly comparing “asymmetric” and “symmetric” ECRAM device architectures, we reveal significant insight into the impact of a non-zero gate-drain open-circuit voltage and its influence on voltage vs. current-controlled gating. We also explore the impact of potentiation/depression write parameters on the symmetry, linearity, and dynamic range of the device response. Together, initial results from optimizing structure and programming approaches yielded unprecedented G max /G min ratios of >1,000 for ECRAM and hundreds of tunable memory states with excellent linearity and symmetry. Simulations based on these ECRAM devices further illustrate the promise of this analog memory technology, achieving near 2% classification error in the MNIST digit recognition benchmark for a range of training parameters compared to a theoretical best of 1.66% and outperforming other device models extracted from the literature.

AIHWKit↗
DOE OSTI2025

Dispersion Leverage Coronagraph: a nulling coronagraph for use on primary objective grating telescopes

Here, we present the Dispersion Leverage Coronagraph (DLC), a variation of the Achromatic Interfero Coronagraph (AIC) that is designed for optical systems featuring large, dispersive primary objective gratings. DLC was originally designed for the Diffractive Interfero Coronagraph Exoplanet Resolver (DICER), a notional 20 m class infrared space telescope utilizing the enhanced one-dimensional angular resolution of large diffraction gratings to discover and characterize near-Earth exoplanets. Here we develop the theoretical foundation for DLC, and apply it to DICER as an example use case. We derive important properties of the DLC system, including focal plane transmission maps, stellar leakage, residual optical path difference tolerance, and pointing error/jitter considerations. Ultimately, we found that DLC effectively nulls an on-axis target across the entire spectrum in the focal plane, allowing for 2D/λ diffraction-limited imaging. It requires asymmetrical fine-guidance tolerances on pointing error/jitter. We work through a benchmark DICER design, explaining the need for a second disperser to reduce background from Zodiacal light, and showing that it could plausibly find and characterize ∼4 nearby, habitable exoplanets around Sun-like stars in a 7-year mission; about 30% of the habitable exoplanets within 8 pc were found in our simulation. The DLC may be useful for any application requiring extremely high-resolution, close companion spectroscopy.

AIC↗
DOE OSTI2023

Design of a Multi-Chemistry Battery Pack System for Behind-the-Meter Storage Applications

Battery management systems (BMS) are essential for a battery pack's safe operation and longevity. This paper presents an active balancing method-based BMS for different cell chemistry structures to be used in behind-the-meter storage (BTMS) applications. The proposed system utilizes modular isolated dual active bridge (DAB) DC/DC converters to actively balance the battery pack through a low voltage (LV) bus. A supervisory controller monitors all the cell voltage, current, and state of charge (SOC) values. Based on the estimation of the SOCs, reference currents for the DAB converters are generated by the supervisory controller. Detailed modeling and the control approach of the modular DAB converters are presented in the paper. Moreover, the control strategy of the supervisory control is also analyzed. The proposed method and structure can be extended to any combination of the number of cells to design the battery pack. Simulation results are provided for a system consisting of three cells in parallel to form a cell block and three cell blocks in series to form the battery module. Experimental results are provided for three modular DAB converters operating with a LiFeMnPO4 prismatic cell with 3.2V, 20Ah rated values.

battery management systems↗
DOE OSTI2022

Magnetic dipole emission in resonant metal-dielectric-metal structures

We explore a possibility to control magnetic dipole emission with plasmonic cavities, placing Eu 3+ emitters inside profile-modulated metal-dielectric-metal structures. Significant variations in the branching ratio of the magnetic and electric dipole transitions are observed as the function of the thickness of the intermediate layer. The experimental results are confirmed with numerical simulations which account for cavity and gap plasmon resonances and predict modifications in the spontaneous emission spectrum as the function of the gap size and a strong directionality of the emission for small thicknesses of the intermediate layer. The implications of having a competition between electric and magnetic dipole relaxation channels in Eu 3+ are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Development of a DC Distribution Testbed for High-Power EV Charging

This paper explores the design and implementation of a power hardware-in-the-loop (P-HIL) setup for DC distribution infrastructure integrated with high-power charging (HPC) of electric vehicles (EVs), DC loads, and sources. The utilization of DC distribution holds significant potential for enhancing the operation of a HPC station architecture. However, there are challenges establishing a DC charging hub including interoperability, commoditization, distributed energy resource integration, stability, DC protection, and lack of common system level controllers. To address these challenges, a testing setup is required that accommodates commercial off-the-shelf (COTS) products to evaluate different use cases at rated power and voltage levels. The developed P-HIL setup features a dedicated DC charging hub, DC-coupled chargers, DC loads/sources, DC protection, and a communication architecture. The integrated P-HIL system provides a versatile testing environment to address technology and interoperability gaps and implements a smart energy management system (SEMS). This platform enables comprehensive and robust testing of COTS devices, charger prototypes, SEMS controllers and protection schemes, which together will accelerate transition to EVs at scale. The setup is tested for various use-cases at full-scale, integrating 950 V DC bus voltage, 660 kW grid-tied inverter, 150 kW COTS charger, and 100 kW energy storage system within an open-source SEMS platform.

ADVANCED PROPULSION SYSTEMS,POWER TRANSMISSION AND↗