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,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Current-Driven Symmetry Breaking and Spin–Orbit Polarization in Chiral Wires

The spin dynamics of electrons in chiral molecular systems remains a topic of intense interest, particularly regarding whether geometric chirality inherently induces spin polarization in current-carrying electrons. In this work, we employ ab initio real-time time-dependent density functional theory (rt-TDDFT) to directly simulate the interplay among charge current, spin, and orbital. This realtime tracking extends beyond perturbative treatments, and we analyze how nonequilibrium currents effectively lift the symmetry constraints of screw rotation and time-reversal symmetry. We find that the emergence of spin and orbital angular momenta is dynamically correlated with a concomitant loss of translational (linear) momentum, which we interpret as an intrinsic consequence of current-driven symmetry lowering. The implications of this mechanism for chirality-induced spin selectivity and spintronics device design are discussed.

chiral wire↗

NASA NTRS1992

Correlation and prediction of dynamic human isolated joint strength from lean body mass

A relationship between a person's lean body mass and the amount of maximum torque that can be produced with each isolated joint of the upper extremity was investigated. The maximum dynamic isolated joint torque (upper extremity) on 14 subjects was collected using a dynamometer multi-joint testing unit. These data were reduced to a table of coefficients of second degree polynomials, computed using a least squares regression method. All the coefficients were then organized into look-up tables, a compact and convenient storage/retrieval mechanism for the data set. Data from each joint, direction and velocity, were normalized with respect to that joint's average and merged into files (one for each curve for a particular joint). Regression was performed on each one of these files to derive a table of normalized population curve coefficients for each joint axis, direction, and velocity. In addition, a regression table which included all upper extremity joints was built which related average torque to lean body mass for an individual. These two tables are the basis of the regression model which allows the prediction of dynamic isolated joint torques from an individual's lean body mass.

Pandya, Abhilash K.↗

DOE OSTI2025

Nonlinear Topological Photonics: Capturing Nonlinear Dynamics and Optical Thermodynamics

Combining multiple optical resonators or engineering dispersion of complex media has provided an effective method for demonstrating topological physics controlling photons in unprecedented ways such as unidirectional light propagation and spatially localized modes between an interface or on a corner. Further, adding nonlinear responses to those topological photonic systems has enabled achieving diverse phases of photons in both space and time, allowing for more functionalities in photonic devices that provide a new playground for studying dynamic features of nonlinear topological systems. However, most methods for describing nonlinear topological photonic systems rely on linear topological theories, making it challenging to accurately characterize the topology of nonlinear systems. Thus, substantial efforts have focused on rigorously describing nonlinear topological phases and developing effective tools to analyze nonlinear topological effects. Meanwhile, coupled multimode optical waveguides with nonlinear dynamic responses provide an excellent platform for the statistical description of photons, opening a new paradigm called “optical thermodynamics”. This review will introduce the basic concepts of nonlinear topological photonics and the recent development of theoretical approaches focusing on data-driven approaches for creating phase diagrams as well as the spectral localizer framework and the pseudospectrum method for understanding optical nonlinearities in topological systems. In addition, the new concept of optical thermodynamics will be introduced with some recent theoretical works.

Topological photonics↗

DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Diversity in the Common Fold: Structural Insights into Class D β-Lactamases from Gram-Negative Pathogens

Class D β-lactamases (DBLs) represent a major threat to antibiotic efficacy by hydrolyzing β-lactam drugs, including last-resort carbapenems, thereby driving antimicrobial resistance in Gram-negative bacteria. The enzymes share a structurally conserved two-domain α/β architecture with seven active-site motifs and three flexible extended loops (the P-loop, Ω-loop, and newly designated B-loop) that surround the active site. While each of these loops is known to influence enzyme function, their coordinated roles have not been fully elucidated. To investigate the significance of their interplay, we compared the sequences and crystal structures of 40 DBLs from clinically relevant Gram-negative pathogens and performed molecular dynamics simulations on selected representatives. Combined structural and dynamical analyses revealed a strong correlation between B-loop architecture and carbapenemase activity in the pathogens Klebsiella and Acinetobacter, particularly regarding loop length and spatial organization. These findings emphasize the B-loop’s critical contribution, in concert with the P- and Ω-loops, in tuning active site versatility, substrate recognition, catalytic activity, and structural stability. A deeper understanding of how these motifs and loops govern DBL function may inform the development of novel antibiotics and inhibitors targeting this class of enzymes.

antibiotic resistance↗

DOE OSTI2025

The Two Arctic Wintertime Boundary Layer States: Disentangling the Role of Cloud and Wind Regimes in Reanalysis and Observations During MOSAiC

The wintertime central Arctic atmosphere comprises a radiatively clear and a radiatively opaque state, which are linked to synoptic forcing and mixed-phase clouds. Weather and climate models often lack process representations surrounding these states, but prior work mostly treated the problem as an aggregate of synoptic conditions, resulting in partially overlapping biases. Here, we disaggregate the Arctic states and confront ERA5 reanalysis with observations from the MOSAiC campaign over the central Arctic sea ice during winter 2019/2020. Low-level winds and liquid water path (LWP) are combined to derive different synoptic classes. Results show that the clear state is primarily formed by weak/moderate winds and the absence of liquid-bearing clouds, while strong winds and enhanced LWP primarily form the radiatively opaque state. ERA5 struggles to reproduce these basic statistics, shows too weak sensitivity of thermal radiation to synoptic forcing, and overestimates thermal radiation for similar LWP amounts. The latter is caused by a warm bias, which has a pronounced inversion structure and is largest in clear and calm conditions. Under strong synoptic forcing, the warm bias is constant with height and discrepancies in mixed-phase cloud altitude appear. Separating synoptic conditions is regarded as useful for process-oriented evaluation of the Arctic troposphere in models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

What are the Impacts on the HVAC System when it Provides Frequency Regulation? – A Comprehensive Case Study with a Multi-Zone Variable Air Volume (VAV) System

Heating, ventilation, and air-conditioning (HVAC) systems in commercial buildings have proven to be a great asset for providing secondary frequency regulation (FR) service to the power grid. However, how the provision of FR service affects the HVAC system Quality of Service (QOS) is still largely unknown. In particular, there exists a controversy regarding the energy efficiency impact. This paper investigates how a HVAC system performance and its associated controls are impacted during and after the FR service. For the assessment, we developed a dynamic Modelica-based medium office building model as a virtual testbed with the Air Handling Unit (AHU) fans providing the FR service. Different regulation capacities (demand levels) and standard test signals under various load profiles in the U.S. are considered. The results show that FR has little or no impact on the HVAC operation and occupant thermal comfort while the building is providing the service if an appropriate regulation capacity is pre-determined. However, if the regulation capacity is overestimated, the quality of the FR service provided by the HVAC system and the HVAC energy efficiency will be affected due to the equipment operating constraints. In addition, the AHU supply air temperature fluctuates, and the controller cannot maintain the desired setpoint as the FR magnitude is increased. After the FR service time, the system takes a certain period to return to normal operation. The recovery period extends with the increase in the FR magnitude. For this case study, there exists a control conflict between the FR control and the existing HVAC control (from ASHRAE Guideline 36-2018) during high load conditions, causing the fan speed to saturate. The fan power consumption could increase significantly; up to 77% increase compared to the baseline was observed. This case study demonstrates that the HVAC system controls need to be enhanced and adapted to minimize the interaction between the FR controller and the indoor environmental controls in order to preserve the HVAC system QOS while (during and after) providing satisfactory FR service.

Lu, Xing↗

DOE OSTI2024

Unraveling the multi-step crystallization mechanism of polytetrafluoroethylene, modified polytetrafluoroethylene, and their nanocomposites with boron nitride nanobarbs: Experimental insights and theoretical analysis

The non-isothermal crystallization behavior and kinetics of polytetrafluoroethylene (PTFE) composites with boron nitride nanobarb (BNNB), a new generation nanostructure with unique surface morphology and mechanical “barbs” have been analyzed, understanding these properties is essential for their high-end applications as thermal interface materials (TIM) for microwave, 5G and microelectronic devices. The analysis of the crystallization parameters includes crystallization onset, peak and end temperatures, crystallization half-life and overall crystallinity of PTFE, modified PTFE and their BNNB composites. The results were further analyzed using theoretical models such as the combined Avrami-Ozawa model. It was found that BNNB supports crystallization in the modified PTFE but shows minimal effect on the crystallization of PTFE. Due to the limitation of the classical theoretical models used above in fully characterizing the multi-step crystallization process of PTFE, an in-depth analysis using the model-free advanced isoconversional computation was used to characterize the PTFE crystallization based on the evolution of activation energy with fractional crystallinity and for the first time with temperature. Three kinetic regions were identified in the crystallization mechanism. Here, this study investigated the molecular organization and microstructural evolution of PTFE, modified PTFE and their composites during non-isothermal crystallization using advanced X-ray scattering measurements. An insight into the changes undergone by the material's microstructural units including crystallite size and morphology, lamellar thickness and lamellar interfacial layer thickness, and crystallographic phase dynamics during non-isothermal cooling from the melt, was provided here in this work. The effect of copolymer modification of PTFE and the inclusion of pristine and functionalized BNNB (a thermally conductive and electrically insulating ceramic) are both new investigations that provide valuable knowledge for the development of materials with strong matrix-nanofiller interaction and guidance for optimizing sintering and cooling cycles, two key steps in PTFE processing that largely affect the material microstructural features. Overall, the result of the three-part investigation demonstrates that BNNB supports crystallization in the modified PTFE up to 20 wt% concentration and at low and high cooling rates typically used in the industrial processing of PTFE.

36 MATERIALS SCIENCE↗

DOE OSTI2026

A fast and robust computational modeling approach for density and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near-net-shape parts with high material utilization and uniform microstructures. PM-HIP is frequently used for producing small-scale parts with complicated geometries and is potentially economical for producing large-scale parts. However, excessive post-HIP shape distortions can reduce its effectiveness and economic advantage, especially for larger parts. A PM-HIP computational model can predict and help mitigate these distortions. However, due to complex deformation mechanisms and thermo-mechanical coupling present in PM-HIP processes, these non-linear computational models sometimes become numerically unstable. The numerical instabilities in these models can lead to very slow convergence or no convergence at all, which often translates to slow and unreliable models. These limitations are more pronounced in large models with complicated geometries. Hence, in this work, an alternative modeling approach is presented that improves numerical stability and computational performance. The presented approach achieves these improvements through approximating the fully coupled thermo-mechanical PM-HIP model as a decoupled model and adding inertial damping to the model’s mechanical part. In conclusion, a comparison with the fully coupled model indicated a slight dip in prediction accuracy (<5% error) but significant improvements in numerical stability (>20 times larger time step size) and computational performance (5-10 times speed-up with less computational resource usage) when using the presented approach.

Hot isostatic pressing↗

DOE OSTI2021

Elevated levels of radiocarbon in methane dissolved in seawater reveal likely local contamination from nuclear powered vessels

Measurements of the natural radiocarbon content of methane ( 14 C-CH 4 ) dissolved in seawater and freshwater have been used to investigate sources and dynamics of methane. However, during investigations along the Atlantic, Pacific, and Arctic Ocean Margins of the United States, as well as in the North American Great Lakes, some samples revealed highly elevated 14 C-CH 4 values, as much as 4–5 times above contemporary atmospheric 14 C-CH 4 levels. Natural production of the 14 CH 4 isotopologue is too low to cause these observations nor can it explain the variations in location and depth. Numerous lab and field validation tests and blanks, as well as the relatively small number of samples that display these elevated values, all suggest that these signals are not derived from an unknown procedural issue. Here, in this work, we suggest that the byproducts of nuclear power generation include localized discharges of the 14 CH4 isotopologue into marine and aquatic environments, severely altering the measured 14 C-CH 4 isotopic signals. Since several of our sample sites are distant from on-land nuclear powerplants, we conduct further calculations concluding that the most elevated anomalies in 14 C-CH 4 likely originate with discharge from nuclear-powered vessels.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Permeability Hysteresis From Microchannels Opening During Dissolution/Reprecipitation Cycle

Abstract Permeability is a critical parameter for geological resources characterization. Its evolution with respect to porosity is particularly interesting and many research initiatives focus on deriving such relationships, to understand some hydraulic impacts of microstructure alteration. Permeability evolution from chemical reactions, for instance, can become complex as flow channels may open during rock dissolution. In this contribution, we show that this phenomenon can lead to irregular porosity‐permeability curves and permeability hysteresis after reprecipitation. Current approaches describing permeability as a simple function of porosity can therefore not capture this behavior, and we advocate instead the use of dynamic modeling for such scenarios. We demonstrate our approach by modeling a dissolution/precipitation cycle for a unit cell pore channel and quantify the process at larger scale on three different rock samples, whose microstructures are reconstructed from segmented micro‐Computerised Tomography scans.

Lesueur, Martin↗

DOE OSTI2025

Cloud Feedback Uncertainty in the Equatorial Pacific Across CMIP6 Models

Cloud feedback is the largest uncertainty in estimating Equilibrium Climate Sensitivity. In this study we focus on the equatorial Pacific, where CMIP6 model cloud feedback spread is notably large. Cloud radiative effects in this region are relevant for the global climate. Our findings show that models predict a consistent shift towards the ascent regime in response to El Nino-like sea surface warming. Models diverge in terms of the radiative impact due to differences in cloud characteristics in ascent and subsidence regimes. Using the observed relationship between circulation regime and cloud radiative effect, we find a reduction in the regional mean cloud feedback estimate from 0.77 to 0.22 W m -2 K -1 , though this does not substantially lessen the model spread in total feedback. Pathways to reduce this spread include: improving confidence in the regional ocean warming pattern, and using observations and models to understand cloud type and circulation interactions.

CMIP6↗

DOE OSTI2025

Atmospheric wind energization of ocean weather

Ocean weather comprises vortical and straining mesoscale motions, which play fundamentally different roles in the ocean circulation and climate system. Vorticity determines the movement of major ocean currents and gyres. Strain contributes to frontogenesis and the deformation of water masses, driving much of the mixing and vertical transport in the upper ocean. While recent studies have shown that interactions with the atmosphere damp the ocean’s mesoscale vortices O(100) km in size, the effect of winds on straining motions remains unexplored. Here, we derive a theory for wind work on the ocean’s vorticity and strain. Using satellite and model data, we discover that wind damps strain and vorticity at an equal rate globally, and unveil striking asymmetries based on their polarity. Subtropical winds damp oceanic cyclones and energize anticyclones outside strong current regions, while subpolar winds have the opposite effect. A similar pattern emerges for oceanic strain, where subtropical convergent flow is damped along the west-equatorward east-poleward direction and energized along the east-equatorward west-poleward direction. These findings reveal energy pathways through which the atmosphere shapes ocean weather.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Soil moisture gradients strengthen mesoscale convective systems by increasing wind shear

Mesoscale convective systems are a class of storm linked to extensive flooding and other destructive hazards in many regions globally. In West Africa, soil moisture impacts provide a valuable source of predictability for mature storm hazards, but little is known about mature storm sensitivity to soil moisture in other climatic regions. Here we use a storm track dataset, satellite observations and reanalysis fields to investigate the response of mature storms to soil moisture in seven global storm hotspots—West Africa, India, South America, South Africa, Australia and the United States Great Plains. We demonstrate that mesoscale soil moisture gradients (~500 km) can enhance storms by driving increased vertical wind shear conditions, a crucial ingredient for storm organization, through the strengthening of atmospheric temperature gradients. This is evidenced by a 10–30% increase in precipitation feature size and rainfall for the largest storms on days with favourable soil moisture gradients compared with unfavourable gradients. Global simulations confirm that soil moisture gradients influence wind shear. The results demonstrate the importance of soil moisture feedbacks for accurate forecasting of mesoscale convective systems and future projections of extreme events under climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

CAMELSH: A Large-Sample Hourly Hydrometeorological Dataset and Attributes at Watershed-Scale for CONUS

We present CAMELSH (Catchment Attributes and Hourly HydroMeteorology for Large-Sample Studies), the first large-sample hydrometeorological dataset at the hourly scale for the contiguous United States. CAMELSH intergrates hourly meteorological time series, catchment attributes and boundaries from GAGES-II and HydroATLAS for 9,008 catchments across diverse climatic, hydrological, and anthropogenic conditions. In addition, hourly streamflow time series is provided for 3,166 catchments. The dataset spans 45 years (1980–2024) with 11 meteorological variables from the NLDAS-2 forcing dataset, from which we compute nine climate indices related to precipitation, evapotranspiration, seasonality, and snow fraction. Additionally, CAMELSH includes two sets of catchment attributes: 439 from GAGES-II and 195 derived from HydroATLAS. These attributes include factors related to climate, geology, hydrology, river/stream morphology, landscape, nutrient, soil, topography, and anthropogenic influences. Developed in accordance with FAIR (Findability, Accessibility, Interoperability, and Reusability) principles, CAMELSH is the first large-sample dataset at an hourly timescale, supporting machine learning applications for short-term streamflow (flood) prediction and advancing data-driven hydrological research across multiple timescales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

How the choice of exchange–correlation functional affects DFT-based simulations of the hydrated electron

Hydrated electrons are anionic species that are formed when an excess electron is introduced into liquid water. Building an understanding of how hydrated electrons behave in solution has been a long-standing effort of simulation methods, of which density functional theory (DFT) has come to the fore in recent years. The ability of DFT to model the reactive chemistry of hydrated electrons is an attractive advantage over semi-classical methodologies; however, relatively few density functional approximations (DFAs) have been used for the hydrated electron simulations presented in the literature. Here, we simulate hydrated electron systems using a series of exchange–correlation (XC) functionals spanning Jacob’s ladder. We calculate a variety of experimental and other observables of the hydrated electron and compare the XC functional dependence for each quantity. We find that the formation of a stable localized hydrated electron is not necessarily limited to hybrid XC functionals and that some hybrid functionals produce delocalized hydrated electrons or electrons that react with the surrounding water at an unphysically fast rate. Here, we further characterize how different DFAs impact the solvent structure and predicted spectroscopy of the hydrated electron, considering several methods for calculating the hydrated electron’s absorption spectrum for the best comparison between structures generated using different density functionals. None of the dozen or so DFAs that we investigated are able to correctly predict the hydrated electron’s spectroscopy, vertical detachment energy, or molar solvation volume.

Ab-initio molecular dynamics↗

DOE OSTI2024

Verification and Validation of Spectral Element Code for Supercritical CO2 Flow in Vertical Heated Tubes

The investigation of heat transfer in supercritical CO2 (sCO2) has garnered considerable attention in recent decades, given sCO2's potential as a promising working fluid for advanced power conversion cycles. Despite previous research efforts, there are still gaps in our understanding of sCO2 heat transfer, particularly in conditions associated with heat transfer deterioration. To delve into sCO2 heat transfer more comprehensively, we propose employing the high-fidelity computational fluid dynamics code NekRS to simulate sCO2 flow using the large eddy simulation technique. Through graphics processing unit acceleration, NekRS achieves a higher computational speed than traditional CPU-based systems. However, before using NekRS in practical applications involving sCO2, it is imperative to perform verification and validation. Here, this paper presents our efforts to verify and validate the NekRS code's capability for simulating sCO2 using heated vertical tubes, where heat transfer deterioration usually happens. To accommodate the unique properties of sCO2, we have modified the NekRS code by integrating third-party property modules, such as REFPROP and PROPATH. Our simulations are compared with experimental and numerical data from the literature, instilling confidence in leveraging NekRS for future engineering applications. Our simulations also reveal that the accuracy of the property module significantly impacts the results, with REFPROP outperforming PROPATH for sCO2 properties. Additionally, we observed that, depending on the flow direction, buoyancy can either enhance or suppress turbulence in sCO2 flow. In upward flow, under certain conditions, the suppressed turbulence leads to heat transfer deterioration, resulting in elevated wall temperatures.

NekRS↗

DOE OSTI2025

Wave topology in Hall magnetohydrodynamics

Hall magnetohydrodynamics (HMHD) extends ideal MHD by incorporating the Hall effect via the induction equation, making it more accurate for describing plasma behavior at length scales below the ion skin depth. Despite its importance, a comprehensive description of the eigenmodes in HMHD has been lacking. In this work, we derive the complete spectrum and eigenvectors of HMHD waves and identify their underlying topological structure. We prove that the HMHD wave spectrum is homotopic to that of ideal MHD, consisting of three distinct branches: the slow magnetosonic-Hall waves, the shear Alfvén-Hall waves, and the fast magnetosonic-Hall waves, which continuously reduce to their ideal MHD counterparts in the limit of vanishing Hall parameter. Contrary to a recent claim [Mahajan, Sharma, and Lingam, Phys. Plasmas 31, 090701 (2024)], we find that HMHD does not admit any additional wave branches beyond those in ideal MHD. In conclusion, the key qualitative difference lies in the topological nature of the HMHD wave structure: it exhibits nontrivial topology characterized by a Weyl point—an isolated eigenmode degeneracy point—and associated nonzero Chern numbers of the eigenmode bundles over a 2-sphere in 𝐤-space surrounding the Weyl point.

Alfvén waves↗