SEARCH · Engineering Papers

Results for “unit 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,495 records · Page 83

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ ≈ 1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ ≈ 10 –2 –10 2 . Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. As a result, the generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω ≳ 2⁢ω p at strongly coupled conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Large-Angle Rocking Beam Electron Diffraction of Large Unit Cell Crystals Using Direct Electron Detector

We report a large-angle rocking beam electron diffraction (LARBED) technique for electron diffraction analysis. Diffraction patterns are recorded in a scanning transmission electron microscope (STEM) using a direct electron detector with large dynamical range and fast readout. We use a nanobeam for diffraction and perform the beam double rocking by synchronizing the detector with the STEM scan coils for the recording. Using this approach, large-angle convergent beam electron diffraction (LACBED) patterns of different reflections are obtained simultaneously. By using a nanobeam, instead of a focused beam, the LARBED technique can be applied to beam-sensitive crystals as well as crystals with large unit cells. Here, this paper describes the implementation of LARBED and evaluates the performance using silicon and gadolinium gallium garnet crystals as test samples. We demonstrate that our method provides an effective and robust way for recording LARBED patterns and paves the way for quantitative electron diffraction of large unit cell and beam-sensitive crystals.

4DSTEM↗

DOE OSTI2024

Dynamic Transformation of High-Architectural Nanocrystal Superlattices upon Solvent Molecule Exposure

The cluster-based body-centered-cubic superlattice (cBCC SL) represents one of the most complicated structures among reported nanocrystal assemblies, comprised of 72 truncated tetrahedral quantum dots per unit cell. Our previous report revealed that truncated tetrahedral quantum dots within cBCC SLs possessed highly controlled translational and orientational order owing to an unusual energetic landscape based on the balancing of entropic and enthalpic contributions during the assembly process. However, the cBCC SL’s structural transformability and mechanical properties, uniquely originating from such complicated nanostructures, have yet to be investigated. Herein, we report that cBCC SLs can undergo dynamic transformation to face-centered-cubic SLs in response to post-assembly molecular exposure. We monitored the dynamic transformation process using in situ synchrotron-based small-angle X-ray scattering, revealing a dynamic transformation involving multiple steps underpinned by interactions between incoming molecules and TTQDs’ surface ligands. Furthermore, our mechanistic study demonstrated that the precise configuration of TTQDs’ ligand molecules in cBCC SLs was key to their high structural transformability and unique jelly-like soft mechanical properties. While ligand molecular configurations in nanocrystal SLs are often considered minor features, our findings emphasize their significance in controlling weak van der Waals interactions between nanocrystals within assembled SLs, leading to previously unremarked superstructural transformability and unique mechanical properties. Finally, our findings promote a facile route toward further creation of soft materials, nanorobotics, and out-of-equilibrium assemblies based on nanocrystal building blocks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Coarse-grained molecular dynamics simulation of solvent-dependent cellulose nanofiber interactions

Associations between cellulose are important both in biofuel production and in the use of cellulose for biomaterials. Cellulose nanofibers (CNFs) are sustainable, strong, light-weight alternatives to traditional materials in manufacturing, but are challenging to obtain due to irreversible aggregation in solution during preparative fibrillation. Therefore, it is imperative to understand the underlying factors driving aggregation with a view to designing solvents that can effectively compete with interfiber interactions, hence reducing aggregation. Molecular dynamics (MD) simulation at atomic detail can provide useful information on local interactions. However, the length and timescales accessible are too short to fully capture association processes. Here, we provide a method for accessing the longer length and timescales required using coarse-grained (CG) MD simulations with a MARTINI force field to calculate the interaction behavior of CNFs in three selected solvents: NaOH-urea-water, acetone, and neat water. The CG results are consistent with our prior all-atom MD and with previous experimental results. While acetone is found not to be an effective solvent, urea and ionic moieties in NaOH-urea-water not only solvate the fibrils but also improve the confinement of water molecules around them as shown by the solvent residence times and mean-square displacements. Overall, the presence of urea and ions reduces the likelihood of aggregation in multi-CNF systems relative to neat water irrespective of whether the hydrophobic or hydrophilic CNF surfaces are interacting. In conclusion, the CG method shows clear promise for selecting potential high-performance solvents for experimental prioritization in bioenergy and biomaterials research in a relatively fast manner as well as for understanding the aggregation and rheological behavior of CNF-solvent systems.

aggregation↗

DOE OSTI2021

How Do Ionic Liquids “Fold” Ionenes? Computational and Experimental Analysis of Imidazolium Polymers Based on Ether and Alkyl Chain Variations Dissolved in an Ionic Liquid

Molecular dynamics (MD) simulations and complementary experiments are used here to understand the inter- and intramolecular structure and conformational properties of ionenes (i.e., cationic polymers formed from condensation reactions) dissolved in an ionic liquid (IL), 1-ethyl-3-methylimidazolium bistriflimide [C 2 mim + ][Tf 2 N - ]. The simulated structural properties are benchmarked against experimental analyses of these same polymers, mainly using dynamic light scattering (DLS) experiments. Four different imidazolium ionenes are considered, corresponding to variations in the chemistry and structure of the repeat units including: (1) poly (decylimidazolium) (PD 10 ); (2) poly (tetraethyleneglycolimidazolium) (PE 10 ); (3) alternating copolymer (P(ED) 5 ); and (4) a block copolymer (PE 5 D 5 ). Detailed computational analyses of the polymer structure and conformational properties were performed, including the radius of gyration, end-to-end distance, torsional distributions, and site-site and spatial distribution functions. Overall, there is a competition between intramolecular associations between the imidazolium groups and the ether sites in the polymers versus the intermolecular interactions of the polymer imidazolium groups with the surrounding anion molecules. The polymer with only ether linkages (PE 10 ) results in the most intramolecular interactions, leading to significant coiling behavior and chain contraction. These strong interactions reduce chain flexibility, but they also result in much more linear chain configurations and alignment of the imidazolium groups. The polymer structure is strongly affected as the concentration of alkyl groups is increased, and we find that the polymer architecture (alternating versus block copolymer) also has an important influence. Furthermore, these detailed observations are important for understanding the connection between the molecular design of imidazolium polymers and their emergent structural properties, which could lead to unique opportunities for creating polymer composites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Understanding the structural mechanics of ligated DNA crystals via molecular dynamics simulation

DNA self-assembly is a highly programmable method to construct arbitrary architectures based on sequence complementarity. Among various constructs, DNA crystals are macroscopic crystalline materials formed by assembling motifs via sticky end association. Due to their high structural integrity and size ranging from tens to hundreds of micrometers, DNA crystals offer unique opportunities to study the structural properties and deformation behaviors of DNA assemblies. For example, enzymatic ligation of sticky ends can selectively seal nicks resulting in more robust structures with enhanced mechanical properties. However, the research efforts have been mostly on experiments involving different motif designs, structural optimization, or new synthesis methods, while their mechanics are not yet fully understood. The complex properties of DNA crystals are difficult to study via experiments alone, and numerical simulation can complement and aid the experiments. The coarse-grained molecular dynamics (MD) simulation is a powerful tool that can probe the mechanics of DNA assemblies. Here, we investigate DNA crystals made of four different motif lengths with various ligation patterns (full ligation, major directions, connectors, and in-plane) using oxDNA, an open-source, coarse-grained MD platform. We found that several distinct deformation stages emerge in response to mechanical loading and that the number and the location of ligated nucleotides can significantly modulate structural behaviors. These findings should be useful for predicting crystal properties and thus improving the design.

DNA crystal↗

DOE OSTI2025

Dynamical phase transition in a strongly hybridized phonon-triplon chain

In this article, we study a dimerized spin-1/2 chain, such as CuGeO 3 , hosting triplon excitations coupled to optical phonons under weak terahertz laser driving. Both phonons and triplons weakly lose energy into the surrounding baths, forming a non-equilibrium steady state. In the strong phonon-triplon coupling regime, phonons near the two-triplon continuum hybridize strongly with triplons. Using mean-field Lindblad dynamics, we show that strong hybridization induces sharp first-order phase transitions—either single or simultaneous double—in the emission spectrum, mainly due to nonlinearities in the chain and the instantaneous breaking of spatial inversion symmetry, as captured by the phonon displacement. Using mean-field Floquet analysis of harmonic modes in both sectors, we analytically confirm the existence of these phase transitions. Furthermore, we map the complete steady-state phase diagram by varying key control parameters and provide experimentally relevant parameters for observing these transitions in laser-driven CuGeO 3 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Topological structure and dynamics of three-dimensional active nematics

Topological structures are effective descriptors of the nonequilibrium dynamics of diverse many-body systems. For example, motile, point-like topological defects capture the salient features of two-dimensional active liquid crystals composed of energy-consuming anisotropic units. We dispersed force-generating microtubule bundles in a passive colloidal liquid crystal to form a three-dimensional active nematic. Light-sheet microscopy revealed the temporal evolution of the millimeter-scale structure of these active nematics with single-bundle resolution. The primary topological excitations are extended, charge-neutral disclination loops that undergo complex dynamics and recombination events. Our work suggests a framework for analyzing the nonequilibrium dynamics of bulk anisotropic systems as diverse as driven complex fluids, active metamaterials, biological tissues, and collections of robots or organisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

A bilevel multistage stochastic self-scheduling model with indivisibilities for trading in the continuous intraday electricity market

In this paper, we study the profit maximization problem of a virtual power plant trading in the continuous intraday electricity market. Our virtual power plant model is compatible with renewable, and thermal assets, covering a range of virtual power plants currently participating in energy markets. We model the trading problem as a bilevel multistage stochastic program. The upper level of the problem accounts for the profit maximization of the virtual power plant with explicit modeling of the technical constraints of the operational status of the thermal power plant including minimum start-up and shut-down times, ramp-up and ramp-down rates, and minimum generation level. The upper level also decides which continuous and indivisible (fill-or-kill) orders are submitted to the market. The lower-level problem accounts for the clearing of the continuous intraday market, i.e., matching of buy and sell orders. Because of the presence of fill-or-kill orders, the lower-level problem is mixed-integer, which prevents its direct conversion to a single-level problem using duality. In order to solve this challenging problem, we develop a convex-hull extended formulation for the lower-level problem, apply duality theory to obtain a single-level stochastic equivalent formulation, and employ McCormick envelopes to turn the problem into a multistage stochastic mixed-integer linear problem, which we solve using the stochastic dual dynamic integer programming algorithm. We conduct numerical experiments and analyze the optimal trading behavior of a virtual power plant trading in an ideal continuous market without arbitrage.

Bilevel multistage stochastic programming problem↗

DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

DOE OSTI2025

Application of a temporal multiscale method for efficient simulation of degradation in PEM Water Electrolysis under dynamic operating conditions

Hydrogen is emerging as a vital energy carrier, driven by the need to reduce carbon emissions. Proton Electrolyte Membrane Water Electrolysis (PEMWE) enables hydrogen production under fluctuating renewable power conditions but requires improved understanding and stability of the anode catalyst layer under dynamic operating conditions, especially with low noble metal loadings. Long-term degradation experiments are both time-consuming and costly; therefore, a systematic, model-aided approach is essential. In the present work, a temporal multiscale method is applied to reduce the computational effort of simulating long-term degradation processes in PEMWE, with an exemplary focus on catalyst dissolution. A mechanistic model incorporating the oxygen evolution reaction, catalyst dissolution, and hydrogen permeation from the cathode to the anode was hypothesized and implemented. In this way, the local periodicity of transport and reaction processes in dynamic PEMWE operation, which influence the gradual degradation of the catalyst layer, is captured. The temporal multiscale method significantly reduces the computational effort of simulation, decreasing processing time from hours to mere minutes. This efficiency gain is attributed to the limited evolution of Slow-Scale variables during each period of time P of the Fast-Scale variables. Consequently, simulation is required only until local periodicity is achieved within each Slow-Scale time step. Hence, the fully resolved dynamic problem is decoupled into these two scales, employing a heterogeneous multiscale technique. The developed approach effectively accelerates parameter estimation and predictive simulations, supporting systematic modeling of PEMWE degradation under dynamic conditions.

08 HYDROGEN↗

NASA NTRS1993

Architecture for Adaptive Intelligent Systems

We identify a class of niches to be occupied by 'adaptive intelligent systems (AISs)'. In contrast with niches occupied by typical AI agents, AIS niches present situations that vary dynamically along several key dimensions: different combinations of required tasks, different configurations of available resources, contextual conditions ranging from benign to stressful, and different performance criteria. We present a small class hierarchy of AIS niches that exhibit these dimensions of variability and describe a particular AIS niche, ICU (intensive care unit) patient monitoring, which we use for illustration throughout the paper. We have designed and implemented an agent architecture that supports all of different kinds of adaptation by exploiting a single underlying theoretical concept: An agent dynamically constructs explicit control plans to guide its choices among situation-triggered behaviors. We illustrate the architecture and its support for adaptation with examples from Guardian, an experimental agent for ICU monitoring.

ACTIVE CONTROL↗

DOE OSTI2026

Electromagnetic Transient Modeling of Large Data Centers for Grid-Level Studies

The magnitude and complexity of electricity usage patterns from large data centers are having significant impacts on the operation and dynamics of the power grid; grid operators and planners require a range of specialized data center models to properly evaluate these impacts and specify technical solutions as needed. Towards addressing this need, Pacific Northwest National Laboratory (PNNL) has developed a library of electromagnetic transient (EMT) models for grid-level studies of data centers called the data center model library (DML). This report describes how the DML was created and how it may properly be used. The models present in the DML are generic models; subject matter expertise and additional technical data are needed to modify these models before they can represent any real data center. However, they will significantly reduce the level of effort required to develop site-specific models and can serve as a common starting point to guide industry towards a more refined consensus. Most of the models within DML are dedicated to representing the power electronics interfaces commonly used in modern data centers, such as double-conversion uninterruptible power supplies and single-phase power factor correction converters. These models are intended for use in grid-level studies and are a simplified aggregation of many small components. That said, background material on the physical and electrical design of large data centers is provided as companion material so that users can be aware of many of the details which have been omitted or streamlined as a matter of practical necessity. Additionally, guidance on the application of EMT analysis for data center interconnection studies is provided, which aids users in identifying when the DML is necessary and what sort of additional model development may be necessary for conducting real-world studies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2026

Validated Reactive Force Field Quantifies MXene Interfacial Properties, Mechanics, and Thermal Transport

MXenes combine rich surface chemistry, mechanical strength, and high conductivity for a multitude of emerging applications. Predictive modeling supports accelerated materials designs and has been limited by the absence of validated and transferable force fields. Here, we introduce an interpretable, reactive INTERFACE force field (IFF and IFF-R) for Ti 3 C 2 T x MXenes that is trained based on chemical knowledge and achieves quantitative agreement with experiments across lattice parameters (<0.5%), density (<0.2%), liquid contact angles, Raman spectra, and the in-plane elastic modulus (∼320 GPa). The models cover surface terminations from hydroxyl (−OH) to fluorine (−F) groups and are extensible to other chemistries. We introduce pH-resolved surface chemistry and identify dopamine adsorption mechanisms at MXene–aqueous interfaces supported by QCM-D and UV–Vis experiments. The data reveal coplanar and perpendicular binding modes and concentration-dependent multilayer assembly. We predict previously inaccessible properties, including termination-dependent cleavage energies, interlayer shear moduli and dynamic shear failure, nanoindentation and brittle fracture, anisotropic in-plane and out-of-plane thermal conductivities, including the role of defects. Agreement with available experimental data is consistently close and exceeds DFT accuracy across the benchmark properties examined. The IFF/IFF-R model is compatible with CHARMM, AMBER, OPLS, and CVFF force fields for simulations of MXenes with diverse surface terminations, electrolyte interfaces, biointerfaces, and polymer composites without additional parameters. Parameter sets, 3D models, and analysis scripts are provided for community use. The validated, reactive, and transferable IFF framework facilitates predictive design of MXene-based films, membranes, sensing interfaces, and composites.

MXene↗

DOE OSTI2024

Mo 3 S 13 Chalcogel: A High-Capacity Electrode for Conversion-Based Li-Ion Batteries

Despite large theoretical energy densities, metal-sulfide electrodes for energy storage systems face several limitations that impact the practical realization. Here, we present the solution-processable, room temperature (RT) synthesis, local structures, and application of a sulfur-rich Mo 3 S 13 chalcogel as a conversion-based electrode for lithium-sulfide batteries (LiSBs). The structure of the amorphous Mo 3 S 13 chalcogel is derived through operando Raman spectroscopy, synchrotron X-ray pair distribution function (PDF), X-ray absorption near edge structure (XANES), and extended X-ray absorption fine structure (EXAFS) analysis, along with ab initio molecular dynamics (AIMD) simulations. A key feature of the three-dimensional (3D) network is the connection of Mo 3 S 13 units through S–S bonds. Li/Mo 3 S 13 half-cells deliver initial capacity of 1013 mAh g –1 during the first discharge. After the activation cycles, the capacity stabilizes and maintains 312 mAh g –1 at a C/3 rate after 140 cycles, demonstrating sustained performance over subsequent cycling. Such high-capacity and stability are attributed to the high density of (poly)sulfide bonds and the stable Mo–S coordination in Mo 3 S 13 chalcogel. Importantly, these findings showcase the potential of Mo 3 S 13 chalcogels as metal-sulfide electrode materials for LiSBs.

25 ENERGY STORAGE↗

DOE OSTI2022

Tradeoffs among indoor air quality, financial costs, and CO 2 emissions for HVAC operation strategies to mitigate indoor virus in U.S. office buildings

Adapting building operation during the COVID-19 pandemic to improve indoor air quality (IAQ) while ensuring sustainable solutions in terms of costs and CO 2 emissions is challenging and limited in literature. Our previous study investigated different HVAC operation strategies, including increased filtration using MERV 10, MERV 13, or HEPA filters, as well as supplying 100% outdoor air into buildings for a system initially sized for MERV 10 filtration. This paper significantly extends that research by systematically analyzing the potential financial and environmental impact for different locations in the U.S. The previous medium office building system model is improved to account for operation in different climates. New evaluation metrics are created to consider the comprehensive impact of improving IAQ on costs and CO 2 emissions, using dynamic emission factors for electricity generation depending on the location. HVAC operation strategies are studied in five different locations across the United States, with distinct climates and electricity sources. In four of the five locations, MERV 13 filtration offers the best improvement in IAQ per increase in costs and emissions relative to MERV 10. The exception is the mildest climate of San Diego, where use of 100% outdoor air provides the best IAQ with a limited increase in costs and emissions. Finally, a system not sized for HEPA filtration can lead to increased costs and emissions without much improvement in IAQ.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Quantitative multiplexed analysis of gene and protein expression patterns in Yarrowia lipolytica

In this report, we present coordinated observations of protein and mRNA transcript counts at the single-cell level in the oleaginous yeast model Yarrowia lipolytica. The transcription factor Xbp1p regulates entry into a quiescent state, representing a shift of resources to sequestration of nutrients rather than cell division. We observed the responses of wild-type and Δxbp1 cells to protein (by fluorescence) and transcript quantification and localization at both single-cell and population-averaged levels. Data were collected via single-molecule fluorescence in situ hybridization (smFISH) and qPCR under nitrogen depletion, a condition that drives lipid accumulation. These techniques reveal a complex and heterogeneous population of Xbp1p dynamics and downstream regulation. Our findings highlight the need for single-cell resolution analyses to describe cellular dynamics and regulatory processes.

Yarrowia lipolytica↗

DOE OSTI2023

Diffusion of small anti-Schottky clusters in UO 2

The aggregation of irradiation-induced defect clusters in UO 2 leads to the formation of dislocation loops and cavities which contribute swelling and ultimately produce deleterious effects on the fuel. Despite their fundamental role in fuel evolution, the kinetics of these defect clusters are not currently well understood. In this work, we investigate the diffusive behavior of interstitial clusters via Molecular Dynamics (MD) simulations. In this work, our investigation considers a range of defect cluster sizes, each composed of N UO 2 units or anti-Schottky defects. We report a complex cluster size - mobility relation; increasing cluster size corresponds to an increase in mobility up until a critical size (N=4) where diffusivity reaches a maximum, after which the trend reverts and further increasing size corresponds to a decrease in mobility. The rapid migration observed correlates well with the very low barriers reported in a few historic experimental studies. Further analysis of the cluster shape and orientation reveals a strong structural preference for near-planar configurations oriented normal to the $\langle$100$\rangle$ direction. Rotational energy barriers are found to be similar in both magnitude and trend to the observed diffusion barriers. Lastly, connections are drawn to the behavior of large defect clusters including the formation of (1/3) $\langle$111$\rangle$ Frank dislocation loops.

36 MATERIALS SCIENCE↗