Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “root dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Transducer Resolution Effect on Pressure Fluctuations Beneath Hypersonic Turbulent Boundary Layers

The size of a pressure transducer is known to affect the accuracy of measurements of wall-pressure fluctuations beneath a turbulent boundary layer because of spatial averaging over the sensing area of the transducer. In this paper, the effect of finite transducer size is investigated by applying spatial averaging or wavenumber filters to a database of hypersonic wall pressure generated from a direct numerical simulation (DNS) that simulates the turbulent portion of the boundary layer over a sharp 7° half-angle cone at nominally Mach 8. Here, a good comparison between the DNS and the experiment in the Sandia Hypersonic Wind Tunnel at Mach 8 is achieved after spatial averaging is applied to the DNS data over an area similar to the sensing area of the transducer. The study shows that a finite sensor size similar to that of the PCB132 transducer can cause significant attenuation in the root-mean-square and power spectral density (PSD) of wall-pressure fluctuations, and the attenuation effect is identical between cone and flat plate configurations at the same friction Reynolds number. The Corcos theory is found to successfully compensate for the attenuated high-frequency components of the wall-pressure PSD.

42 ENGINEERING↗
DOE OSTI2025

Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings

This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Validation of Streaklines as Recorders of Synoptic Flow Direction in a Deltaic Setting

Knowledge of the flow patterns within distributary systems is key for understanding deltaic hydro- and morpho-dynamics, yet synoptic measurements of flow fields remain virtually nonexistent. As a means of overcoming this problem, a small number of studies have used biogenic surface films as synoptic flow tracers, under the assumption that biofilm streaklines are tangent to the local flow direction. Here we rigorously test this assumption and show that, despite flow patterns that change severely in space and time (over a range >270°), streaklines are relatively accurate synoptic flow tracers for the Wax Lake Delta, in Louisiana. When the incoming discharge was greater than 2400 m3/s with stable or falling tides, the streakline-derived flow direction departed from near bed flow direction measurements of 22.8° (root mean square). When the discharge was greater than 2400 m3/s and the tides were rising greater than 0.03 m/hr, they were accurate within 28.0°. Under conditions of discharge less than 2400 m3/s and tidal change less than a positive 0.03 m/hr, they were accurate within 33.3°, while during low discharge and rising tides they were accurate within 58.9°. Accuracy varied with distance from the delta, with proximal sites having greater precision. Our results demonstrate that a streakline-derived flow direction can characterize the spatiotemporal variability in the flow directions, but that the accuracy is significantly influenced by the hydrodynamic conditions and location within the network.

58 GEOSCIENCES↗
DOE OSTI2018

Gearbox high-speed-stage bearing slip induced by electric excitation in a test facility

The goal of this paper is to investigate the behaviour of the high-speed stage of a wind turbine gearbox during transient grid-loss events triggered at different power levels in a full-scale nacelle dynamometer test. The gear mesh behaviour and the bearing roller slip during these events are characterized by strain measurements placed in the gear teeth roots and embedded in the outer rings of the tapered roller bearing locating pair. From our experiments, we conclude that the high-speed stage experiences a multidimensional loading consisting of bending and torsional deformation during the transient event. This loading propagates throughout the gearbox and results in unfavourable loading conditions at the bearing that have the potential to go beyond the preload during the torque reversals that occur during the grid-loss event. During the torque reversal, roller slip can occur in the bearings. It is observed that the slip behaviour is directly linked to the initial power level of the generator. Electrically induced excitation that typically leads to dynamic torque events can therefore cause mechanical gearbox degradation.

17 WIND ENERGY↗
DOE OSTI2020

Ligand-Binding-Site Refinement to Generate Reliable Holo Protein Structure Conformations from Apo Structures

The first important step in a structure-based virtual screening is the judicious selection of a receptor protein. In cases where the holo protein receptor structure is unavailable, significant reduction in virtual screening performance has been reported. In this work, we present a robust method to generate reliable holo protein structure conformations from apo structures using molecular dynamics (MD) simulation with restraints derived from holo structure binding-site templates. We perform benchmark tests on two different datasets: 40 structures from a directory of useful decoy-enhanced (DUD-E) and 84 structures from the Gunasekaran dataset. Our results show successful refinement of apo binding-site structures toward holo conformations in 82% of the test cases. In addition, virtual screening performance of 40 DUD-E structures is significantly improved using our MD-refined structures as receptors with an average enrichment factor (EF), an EF1% value of 6.2 compared to apo structures with 3.5. Docking of native ligands to the refined structures shows an average ligand root mean square deviation (RMSD) of 1.97 Å (DUD-E dataset and Gunasekaran dataset) relative to ligands in the holo crystal structures, which is comparable to the self-docking (i.e., docking of the native ligand back to its crystal structure receptor) average, 1.34 Å (DUD-E dataset) and 1.36 Å (Gunasekaran dataset). On the other hand, docking to the apo structures yields an average ligand RMSD of 3.65 Å (DUD-E) and 2.90 Å (Gunasekaran). These results indicate that our method is robust and can be useful to improve virtual screening performance of apo structures.

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

Fourier-like Thermal Relaxation of Nanoscale Explosive Hot Spots

Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Accuracy assessment of the large-scale dynamic ocean topography from TOPEX/POSEIDON altimetry

The quality of TOPEX/POSEIDON determinations of the global scale dynamic ocean topography have been assessed by determining mean topography solutions for successive 10-day repeat cycles and by examining the temporal changes in the sea surface topography to identify known features. The assessment is based on the analysis of TOPEX altimeter data cycles 1 through 36. Important errors in the tide model used to correct the altimeter data have been identified. The errors were reduced significantly by use of a new tide model derived with the TOPEX/POSEIDON measurements. Maps of the global 1-year mean topography, produced using four of the most accurate of the marine geoid, show that the largest error in the dynamic ocean topography show expected features, such as the known annual hemispherical sea surface rise and fall and the seasonal variability due to monsoon influence in the Indian Ocean. Changes in the sequence of 10-day topography maps show the development and propagation of an equatorial Kelvin wave in the Pacific beginning in December 1992 with a propagation velocity of approximately 3 m/s. The observations are consistent with observed changes in the equatorial trade winds, and with tide gauge and other in situ observations of the strengthening of the El Nino. Comparison of TOPEX-determine sea surface height at points near oceanic tide gauges shows agreement at the 4 cm root-mean-square (RMS) level over the tropical Pacific. The results show that the TOPEX altimeter data set can be used to map the ocean surface with a temporal resolution of 10 days and an accuracy which is insonsistent with traditional in situ methods for the determination of sea level variations.

Tapley, B. D.↗
NASA NTRS2011

CFD Study of the J2-X LOX Two-Bladed Inducer in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of water flow experiments of the J-2X Oxidizer Turbopump (OTP) 2-bladed inducer with a circumferential groove that were conducted at Marshall Space Flight Center (MSFC). Sub-scale water flow testing results indicate that the circumferential groove greatly reduces synchronous cavitation and subsequent bearing loads at a minimal performance cost, but the energy reappears as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI-Stream CFD program were conducted in order to explore the root physical cause of the HOC. It was found that the axial recirculation pattern in the circumferential groove causes high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that the root cause of HOC is hydrodynamic. That is, the pearls are always present, even when cavitation is not. Comparisons to ongoing water flow experiments will be discussed, as well as predictions for the full-scale OTP inducer.

Coker, Robert↗
NASA NTRS2018

A Study of the Carbon Cycle Using NASA Observations and the GEOS Model

The Goddard Earth Observing System (GEOS) model has been developed in the Global Modeling and Assimilation Office (GMAO) at NASA's Goddard Space Flight Center. From its roots in chemical transport and as a General Circulation Model, the GEOS model has been extended to an Earth System Model based on a modular construction using the Earth System Modeling Framework (ESMF), combining elements developed in house in the GMAO with others that are imported through collaborative research. It is used extensively for research and for product generation, both as a free-running model and as the core of the GMAO's data assimilation system. In recent years, the GMAO's modeling and assimilation efforts have been strongly supported by Piers Sellers, building on both his earlier legacy as an observationally oriented model developer and his post-astronaut career as a dynamic leader into new territory. Piers' long-standing interest in the carbon cycle and the combination of models with observations motivates this presentation, which will focus on the representation of the carbon cycle in the GEOS Earth System Model. Examples will include: (i) the progression from specified land-atmosphere surface fluxes to computations with an interactive model component (Catchment-CN), along with constraints on vegetation distributions using satellite observations; (ii) the use of high-resolution satellite observations to constrain human-generated inputs to the atmosphere; (iii) studies of the consistency of the observed atmospheric carbon dioxide concentrations with those in the model simulations. The presentation will focus on year-to-year variations in elements of the carbon cycle, specifically on how the observations can inform the representation of mechanisms in the model and lead to integrity in global carbon dioxide simulations. Further, applications of the GEOS model to the planning of new carbon-climate observations will be addressed, as an example of the work that was strongly supported by Piers in the last months of his leadership of Earth Science at NASA Goddard.

Pawson, Steven↗
NASA NTRSDate not supplied

On the Solar Open Magnetic Flux Problem

The solar magnetic fields emerging from the photosphere into the chromosphere and corona are comprised of a combination of “closed” (field lines with both ends rooted at the Sun) and “open” (field lines with only one end at the Sun) fields. Since the early 2000’s, the magnitude of total unsigned open magnetic flux estimated by coronal models has been in significant disagreement with in situ spacecraft observations, especially during solar maximum. Estimates of total open unsigned magnetic flux using coronal hole observations (e.g., using extreme ultraviolet (EUV) or Helium (He) are in general agreement with the coronal model results and thus show similar disagreements with in situ observations. This paper provides a brief overview of the problem, summarizes the proposed explanations for the discrepancies, and presents results that strongly support the explanation that the discrepancy is due to dynamics at the open-closed boundary. These results are derived from determination of the total unsigned open magnetic flux, utilizing the Wang-Sheeley-Arge (WSA) model at a particular spatial resolution and different field-line tracing methods. One of these methods produces excellent agreement with in situ observations. Our results imply that active regions residing near the boundaries of mid-latitude coronal holes are the primary source of the missing open flux. Furthermore, they resolve many of the seemingly contradictory facts that have made resolving the open-flux problem so difficult.

Charles Arge↗
DOE OSTI2021

Declining carbohydrate content of Sitka-spruce treesdying from seawater exposure

Abstract Increasing sea levels associated with climate change threaten the survival of coastal forests, yet the mechanisms by which seawater exposure causes tree death remain poorly understood. Despite the potentially crucial role of nonstructural carbohydrate (NSC) reserves in tree survival, their dynamics in the process of death under seawater exposure are unknown. Here we monitored progressive tree mortality and associated NSC storage in Sitka-spruce (Picea sitchensis) trees dying under ecosystem-scale increases in seawater exposure in western Washington, USA. All trees exposed to seawater, because of monthly tidal intrusion, experienced declining crown foliage during the sampling period, and individuals with a lower percentage of live foliated crown (PLFC) died faster. Tree PLFC was strongly correlated with subsurface salinity and needle ion contents. Total NSC concentrations in trees declined remarkably with crown decline, and reached extremely low levels at tree death (2.4% and 1.6% in leaves and branches, respectively, and 0.4% in stems and roots). Starch in all tissues was almost completely consumed, while sugars remained at a homeostatic level in foliage. The decreasing NSC with closer proximity to death and near zero starch at death are evidences that carbon starvation occurred during Sitka-spruce mortality during seawater exposure. Our results highlight the importance of carbon storage as an indicator of tree mortality risks under seawater exposure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Non-Invasive Electrode Impedance Estimation for Optimized Charge Profile Parameterization of Lithium-Ion Batteries

This work presents a non-invasive method for parameterizing a physically motivated equivalent circuit model of lithium-ion batteries using operando electrochemical impedance spectroscopy and time-domain data. The proposed model consists exclusively of linear circuit elements, enabling computationally efficient simulation and real-time implementation on battery management system chips. By integrating frequency- and time-domain measurements, the model accurately estimates internal states such as the negative electrode potential, achieving a root mean square error of 12.3 mV during fast charging validation. Parameterization requires only rate tests with sinusoidal perturbations at three different ambient temperatures, making the approach experimentally accessible. The model reveals key insights into battery behavior, including rate-dependent overpotentials primarily governed by charge transfer kinetics at the positive electrode, and temperature-dependent impedance contributions from both charge transfer and solid-state diffusion processes. Validation using reference electrodes confirms the model’s ability to detect lithium plating onset and reproduce impedance behavior across a wide range of operating conditions. The approach enables in situ optimization of fast charging profiles and lays the foundation for future extensions incorporating aging effects and plating dynamics.

25 ENERGY STORAGE↗
DOE OSTI2023

A computational modeling framework for pre-clinical evaluation of cardiac mapping systems

There are a variety of difficulties in evaluating clinical cardiac mapping systems, most notably the inability to record the transmembrane potential throughout the entire heart during patient procedures which prevents the comparison to a relevant “gold standard”. Cardiac mapping systems are comprised of hardware and software elements including sophisticated mathematical algorithms, both of which continue to undergo rapid innovation. The purpose of this study is to develop a computational modeling framework to evaluate the performance of cardiac mapping systems. The framework enables rigorous evaluation of a mapping system’s ability to localize and characterize (i.e., focal or reentrant) arrhythmogenic sources in the heart. The main component of our tool is a library of computer simulations of various dynamic patterns throughout the entire heart in which the type and location of the arrhythmogenic sources are known. Our framework allows for performance evaluation for various electrode configurations, heart geometries, arrhythmias, and electrogram noise levels and involves blind comparison of mapping systems against a “silver standard” comprised of computer simulations in which the precise transmembrane potential patterns throughout the heart are known. A feasibility study was performed using simulations of patterns in the human left atria and three hypothetical virtual catheter electrode arrays. Activation times (AcT) and patterns (AcP) were computed for three virtual electrode arrays: two basket arrays with good and poor contact and one high-resolution grid with uniform spacing. The average root mean squared difference of AcTs of electrograms and those of the nearest endocardial action potential was less than 1 ms and therefore appears to be a poor performance metric. In an effort to standardize performance evaluation of mapping systems a novel performance metric is introduced based on the number of AcPs identified correctly and those considered spurious as well as misclassifications of arrhythmia type; spatial and temporal localization accuracy of correctly identified patterns was also quantified. This approach provides a rigorous quantitative analysis of cardiac mapping system performance. Proof of concept of this computational evaluation framework suggests that it could help safeguard that mapping systems perform as expected as well as provide estimates of system accuracy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2018

Particle production in ultrastrong-coupling waveguide QED

Understanding large-scale interacting quantum matter requires dealing with the huge number of quanta that are produced by scattering even a few particles against a complex quantum object. Prominent examples are found from high energy cosmic ray showers, to the optical or electrical driving of degenerate Fermi gases. We tackle this challenge in the context of many-body quantum optics, as motivated by the recent developments of circuit quantum electrodynamics at ultrastrong coupling. The issue of particle production is addressed quantitatively with a simple yet powerful concept rooted in the quantum superposition principle of multimode coherent states. This key idea is illustrated by the study of multi-photon emission from a single two-level artificial atom coupled to a high impedance waveguide, driven by a nearly-monochromatic coherent tone. Furthermore, we find surprisingly that the off-resonant inelastic emission lineshape is dominated by broadband particle production, due to the large phase space associated with contributions that do not conserve the number of excitations. Such frequency conversion processes produce striking signatures in time correlation measurements, which can be tested experimentally in quantum waveguides. Furthermore, these ideas open new directions for the simulation of a variety of physical systems, from polaron dynamics in solids to complex superconducting quantum architectures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

A Numerical Model of Unsteady, Subsonic Aeroelastic Behavior

A method for predicting unsteady, subsonic aeroelastic responses was developed. The technique accounts for aerodynamic nonlinearities associated with angles of attack, vortex-dominated flow, static deformations, and unsteady behavior. The fluid and the wing together are treated as a single dynamical system, and the equations of motion for the structure and flow field are integrated simultaneously and interactively in the time domain. The method employs an iterative scheme based on a predictor-corrector technique. The aerodynamic loads are computed by the general unsteady vortex-lattice method and are determined simultaneously with the motion of the wing. Because the unsteady vortex-lattice method predicts the wake as part of the solution, the history of the motion is taken into account; hysteresis is predicted. Two models are used to demonstrate the technique: a rigid wing on an elastic support experiencing plunge and pitch about the elastic axis, and an elastic wing rigidly supported at the root chord experiencing spanwise bending and twisting. The method can be readily extended to account for structural nonlinearities and/or substitute aerodynamic load models. The time domain solution coupled with the unsteady vortex-lattice method provides the capability of graphically depicting wing and wake motion.

Strganac, Thomas W.↗
DOE OSTI2024

Dynamic Modeling and Simulation of a Subcritical Coal-Fired Power Plant under Load-Following Conditions

Dynamic models for power plants that capture realistic general process trends and effects of manipulated variables are needed to improve load-following, while minimizing carbon footprint. In this work, a dynamic modeling approach and simulation results for subcritical coal-fired power plant components are presented. These encompass simulation of the dynamics in the fireside, including the effects of fuel, air combustion, and the dynamics of the entire waterside and power generation sections. This model development enables the simulation and analysis of the important short and long-time scale dynamics of components such as heaters, evaporative loop, and power generation units. Furthermore, additional variables in the power generation section are introduced to improve model accuracy, extending the prediction capability of subcritical power plant models and opening new opportunities for research in operator training, optimization, and advanced model-based controller design that are based on these models. The change in process gain for different ramp rates associated with disturbance signals that affect process variables is also explored and a correlation developed. This provides opportunities to study disturbance rejection control implementation and adaptation for scenarios with such variations in ramp rates. The prediction capabilities of selected components are compared to data available in literature, with the obtained root mean squared error ranges that reflect the model performance and quality of predictions.

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

The Principal Hugoniot of Forsterite to 950 GPa

Forsterite (Mg 2 SiO 4 ) single crystals were shock compressed to pressures between 200 and 950 GPa using independent plate–impact steady shocks and laser–driven decaying shock compression experiments. Additionally, we performed density functional theory–based molecular dynamics to aid interpretation of the experimental data and to investigate possible phase transformations and phase separations along the Hugoniot. We show that the experimentally obtained Hugoniot cannot distinguish between a pure liquid Mg 2 SiO 4 and an assemblage of solid MgO plus liquid magnesium silicate. In conclusion, the measured reflectivity is nonzero and increases with pressure, which implies that the liquid is a poor electrical conductor at low pressures and that the conductivity increases with pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Asymptotic Behavior of Memristive Circuits

The interest in memristors has risen due to their possible application both as memory units and as computational devices in combination with CMOS. This is in part due to their nonlinear dynamics, and a strong dependence on the circuit topology. We provide evidence that also purely memristive circuits can be employed for computational purposes. In the present paper we show that a polynomial Lyapunov function in the memory parameters exists for the case of DC controlled memristors. Such a Lyapunov function can be asymptotically approximated with binary variables, and mapped to quadratic combinatorial optimization problems. This also shows a direct parallel between memristive circuits and the Hopfield-Little model. In the case of Erdos-Renyi random circuits, we show numerically that the distribution of the matrix elements of the projectors can be roughly approximated with a Gaussian distribution, and that it scales with the inverse square root of the number of elements. This provides an approximated but direct connection with the physics of disordered system and, in particular, of mean field spin glasses. Using this and the fact that the interaction is controlled by a projector operator on the loop space of the circuit. We estimate the number of stationary points of the approximate Lyapunov function and provide a scaling formula as an upper bound in terms of the circuit topology only.

97 MATHEMATICS AND COMPUTING↗