Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Modeling”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Winter precipitation predictability in Central Southwest Asia and its representation in seasonal forecast systems

In Central Southwest Asia (CSWA; 22°N to 40°N and 30°E to 70°E), winter (November to February) precipitation contributes up to 70% of the annual mean, but substantial interannual variations exist. Dynamical models exhibit subpar predictability in this region, but the limits of their skills are not well established. Here, we identify the tropical and extratropical forcings that explain ~75% of area-averaged seasonal variability in CSWA winter precipitation. Tropical forcing comes from the indirect El Niño-Southern Oscillation (ENSO) pathway, the leading mode of tropical Indian Ocean precipitation variability. This mode is coupled with ENSO-related Pacific Ocean sea surface temperature variability. A direct ENSO influence on CSWA does not extend beyond its Indian Ocean connection. Extratropical forcing comes from a large-scale mode of internal atmospheric variability. The spatial structure, variability of tropical forcing, and teleconnection with CSWA winter precipitation are skillfully depicted in two seasonal forecasting systems: the fifth-generation seasonal forecasting system (SEAS5) and Seamless System for Prediction and Earth System Research (SPEAR). Extratropical forcing’s spatial structure is also produced skillfully in the two modelling systems; however, the representation of its interannual variability and teleconnection requires improvement. While SEAS5 displays skill in representing extratropical forcing influence on CSWA winter precipitation and marginal skill in reproducing interannual variability, SPEAR has negligible ability in both areas. Consequently, these models have limited predictive skills over CSWA in winter. While improvements in representing extratropical forcing may be inherently limited as it arises from internal atmospheric variability, further research is needed to establish its predictability limits fully.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Deep transfer operator learning for partial differential equations under conditional shift

Transfer learning enables the transfer of knowledge gained while learning to perform one task (source) to a related but different task (target), hence addressing the expense of data acquisition and labelling, potential computational power limitations and dataset distribution mismatches. Here, we propose a new transfer learning framework for task-specific learning (functional regression in partial differential equations) under conditional shift based on the deep operator network (DeepONet). Task-specific operator learning is accomplished by fine-tuning task-specific layers of the target DeepONet using a hybrid loss function that allows for the matching of individual target samples while also preserving the global properties of the conditional distribution of the target data. Inspired by conditional embedding operator theory, we minimize the statistical distance between labelled target data and the surrogate prediction on unlabelled target data by embedding conditional distributions onto a reproducing kernel Hilbert space. We demonstrate the advantages of our approach for various transfer learning scenarios involving nonlinear partial differential equations under diverse conditions due to shifts in the geometric domain and model dynamics. Our transfer learning framework enables fast and efficient learning of heterogeneous tasks despite considerable differences between the source and target domains.

42 ENGINEERING↗
DOE OSTI2020

Programmable topotaxis of magnetic rollers in time-varying fields

We discuss how spatially uniform, time-periodic magnetic fields can be designed to power and direct the migration of ferromagnetic spheres up (or down) local gradients in the topography of a solid substrate. Our results are based on a dynamical model that considers the time-varying magnetic torques on the particle and its motion through the fluid at low Reynolds number. We use both analytical theory and numerical simulation to design magnetic fields that maximize the migration velocity up (or down) an inclined plane. We show how “topotaxis” of spherical particles relies on differences in the hydrodynamic resistance to rotation about axes parallel and perpendicular to the plane. Importantly, the designed fields can drive multiple independent particles to move simultaneously in different directions as determined by gradients in their respective environments. Experiments on ferromagnetic spheres provide evidence for topotactic motions up inclined substrates. The ability to program the autonomous navigation of driven particles within anisotropic environments is relevant to the design of colloidal robots.

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

Positional ordering induced by dynamic steric interactions in superparamagnetic rods

The dynamic motion produced by precessing magnetic fields can drive matter into far-from-equilibrium states. We predict 1D periodic ordering in systems of precessing rods when magnetic interactions between rods remain insignificant. The precession angle of the rods is completely determined by the field's precession angle and the ratio of the field's precession frequency and the characteristic response frequency of the rods. Here, we develop a molecular dynamics model that explicitly calculates magnetic interactions between particles, and we also simulate rods in the limit of a strong and fast precessing magnetic field where inter-rod magnetic interactions are negligible, using a purely steric model. Our simulations show how steric interactions drive the rods from a positionally disordered phase (nematic) to a layered (smectic) phase. As the rod precession angle increases, the nematic–smectic transition density significantly decreases. The minimization of unfavorable steric interactions also induces phase separation in binary mixtures of rods of different lengths. This effect is general to any force that produces precession in elongated particles. This work will advance the understanding and control of out-of-equilibrium soft matter systems.

42 ENGINEERING↗
DOE OSTI2022

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Time propagation of electronic wavefunctions using nonorthogonal determinant expansions

The use of truncated configuration interaction in real-time time-dependent simulations of electron dynamics provides a balance of computational cost and accuracy, while avoiding some of the failures associated with real-time time-dependent density functional theory. However, low-order truncated configuration interaction also has limitations, such as overestimation of polarizability in configuration interaction singles, even when perturbative doubles are included. Increasing the size of the determinant expansion may not be computationally feasible, and so, in this work, we investigate the use of nonorthogonality in the determinant expansion to establish the extent to which higher-order substitutions can be recovered, providing an improved description of electron dynamics. Model systems are investigated to quantify the extent to which different methods accurately reproduce the (hyper)polarizability, including the high-harmonic generation spectrum of H 2 , water, and butadiene.

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

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Hydrodynamic and frictional modulation of deformations in switchable colloidal crystallites

Displacive transformations in colloidal crystals may offer a pathway for increasing the diversity of accessible configurations without the need to engineer particle shape or interaction complexity. To date, binary crystals composed of spherically symmetric particles at specific size ratios have been formed that exhibit floppiness and facile routes for transformation into more rigid structures that are otherwise not accessible by direct nucleation and growth. There is evidence that such transformations, at least at the micrometer scale, are kinetically influenced by concomitant solvent motion that effectively induces hydrodynamic correlations between particles. Here, we study quantitatively the impact of such interactions on the transformation of binary bcc-CsCl analog crystals into close-packed configurations. We first employ principal-component analysis to stratify the explorations of a bcc-CsCl crystallite into orthogonal directions according to displacement. We then compute diffusion coefficients along the different directions using several dynamical models and find that hydrodynamic correlations, depending on their range, can either enhance or dampen collective particle motions. These two distinct effects work synergistically to bias crystallite deformations toward a subset of the available outcomes.

42 ENGINEERING↗
DOE OSTI2024

Command of three-dimensional solitary waves via photopatterning

Multidimensional solitons are prevalent in numerous research fields. In orientationally ordered soft matter system, three-dimensional director solitons exemplify the localized distortion of molecular orientation. However, their precise manipulation remains challenging due to unpredictable and uncontrolled generation. Here, we utilize preimposed programmable photopatterning in nematics to control the kinetics of director solitons. This enables both unidirectional and bidirectional generation at specific locations and times, confinement within micron-scaled patterns of diverse shapes, and directed propagation along predefined trajectories. A focused dynamical model provides insight into the origins of these solitons and aligns closely with experimental observations, underscoring the pivotal role of anchoring conditions in soliton manipulation. Our findings pave the way for diverse fundamental research avenues and promising applications, including microcargo transportation and optical information processing.

Science & Technology - Other Topics↗
DOE OSTI2025

Integrative mapping reveals molecular features underlying the mechanism of nucleocytoplasmic transport

Nuclear pore complexes (NPCs) enable rapid, selective, and robust nucleocytoplasmic transport. To explain how transport emerges from the system components and their interactions, we used experimental data and theoretical information to construct an integrative Brownian dynamics model of transport through an NPC, coupled to a kinetic model of transport in the cell. The model recapitulates key aspects of transport for a wide range of molecular cargoes, including preribosomes and viral capsids. Our model quantifies how flexible phenylalanine-glycine (FG) repeat proteins create an entropic barrier to passive diffusion and how this barrier is selectively lowered in facilitated diffusion by the many transient interactions of nuclear transport receptors with the FG repeats. Selective transport is enhanced by “fuzzy” multivalent interactions, redundant FG repeat mass, coupling to the energy-dependent RanGTP concentration gradient, and exponential dependence of transport kinetics on the transport barrier. Our model will facilitate rational modulation of the NPC and its artificial mimics.

Raveh, Barak [School of Computer Science and Engin↗
DOE OSTI2020

Machine learning for collective variable discovery and enhanced sampling in biomolecular simulation

Classical molecular dynamics simulates the time evolution of molecular systems through the phase space spanned by the positions and velocities of the constituent atoms. Molecular-level thermodynamic, kinetic, and structural data extracted from the resulting trajectories provide valuable information for the understanding, engineering, and design of biological and molecular materials. The cost of simulating many-body atomic systems makes simulations of large molecules prohibitively expensive, and the high-dimensionality of the resulting trajectories presents a challenge for analysis. Driven by advances in algorithms, hardware, and data availability, there has been a flare of interest in recent years in the applications of machine learning – especially deep learning – to molecular simulation. Furthermore, these techniques have demonstrated great power and flexibility in both extracting mechanistic understanding of the important nonlinear collective variables governing the dynamics of a molecular system, and in furnishing good low-dimensional system representations with which to perform enhanced sampling or develop long-timescale dynamical models. It is the purpose of this article to introduce the key machine learning approaches, describe how they are married with statistical mechanical theory into domain-specific tools, and detail applications of these approaches in understanding and accelerating biomolecular simulation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Expanded Accident Scenarios and Frequency Characteristics for a Fast Spectrum Molten Salt Reactor in an Integrated Energy System

Integration of advanced nuclear reactors to industrial processes in an integrated energy system (IES) has many potential advantages. The coupling of nuclear power to varied industrial processes introduces additional accident scenarios unique to the specific systems included in the IES. Similarly, noise in one system of the IES could potentially create negative effects in the other systems. Because of the limited amount of comparable nuclear IES operating experience, there is added value in scoping the inherent behavior of these systems in response to these transients. Here, this paper investigates the safety-related behavior of an open-source IES dynamic model of a fast spectrum molten salt reactor (MSR) powering a regenerative Rankine cycle for electricity production and a hybrid sulfur (HyS) cycle for hydrogen production. The simulations include a loss of cooling water to the Rankine cycle, a loss of sulfur in the HyS cycle, and the frequency characteristics of the entire IES across six orders of magnitude of frequency. The results show safe behavior in response to the two accident scenarios, wherein a disruption in heat transfer leads to increased salt temperatures and a decrease in reactor power. The frequency analysis shows that at higher frequencies, temperature changes propagating through an IES are increasingly damped to the point of negligibility.

Dunkle, Nicholas [Univ. of Tennessee, Knoxville, T↗
DOE OSTI2024

Co-Active Subspace Methods for the Joint Analysis of Adjacent Computer Models

Active subspace (AS) methods are a valuable tool for understanding the relationship between the inputs and outputs of a Physics simulation. In this article, an elegant generalization of the traditional ASM is developed to assess the co-activity of two computer models. This generalization, which we refer to as a Co-Active Subspace (Co-AS) Method, allows for the joint analysis of two or more computer models allowing for thorough exploration of the alignment (or non-alignment) of the respective gradient spaces. We define co-active directions, co-sensitivity indices, and a scalar “concordance” metric (and complementary “discordance” pseudo-metric) and we demonstrate that these are powerful tools for understanding the behavior of a class of computer models, especially when used to supplement traditional AS analysis. Details for efficient estimation of the Co-AS and an accompanying R package (concordance) are provided. Practical application is demonstrated through analyzing a set of simulated rate stick experiments for PBX 9501, a high explosive, offering insights into complex model dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A hardware-in-the-loop (HIL) testbed for cyber-physical energy systems in smart commercial buildings

In recent years, there has been a growing trend toward the development of smart buildings that rely on cyber-physical systems (CPS) to optimize occupant comfort, safety, and energy efficiency. To ensure the reliable and efficient operation of CPS with designed control strategies, it is important to evaluate their performance under various scenarios before deploying them in the real world. This is where a Hardware-in-the-loop (HIL) testbed designed for studying sensor and control-related studies in smart buildings can be highly valuable. With the growing threat of cyber-attacks and physical faults targeting smart buildings, it is essential to ensure the security of building operations. A HIL testbed can emulate cyber-attack and physical fault scenarios, allowing researchers to develop and test threat detection and mitigation algorithms. This enables researchers to identify potential issues and optimize the algorithms in a safe and controlled environment before they are deployed in real-world settings, reducing the risk of failures that can negatively impact occupant comfort, safety, and energy efficiency. Therefore, this paper developed a HIL testbed designed for cyber-physical energy systems (e.g. buildings automation system (BAS)) in smart commercial buildings. The HIL testbed is comprised of a real-time building and Heating, Ventilation, and Air-Conditioning (HVAC) emulator using Modelica-based dynamic models, a set of BAS controllers, and a BAS computer server. The data generation capability of the HIL testbed is demonstrated by tracking normal and faulty operating data in the BAS, as well as monitoring detailed network traffic in the local BAS network. Here, this study further demonstrates the HIL testbed’s capability by conducting case studies on real-time physical fault and cyber-attack experiments using a Department of Energy (DOE) prototype commercial building. It is anticipated that the fully functional HIL testbed will be utilized for a variety of sensor and control-related studies, including but not limited to testing, developing, validating of different HVAC control strategies, fault detection & diagnosis, energy monitoring and analysis, cyber security study, etc.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Review: Turbulence dynamics during the pedestal evolution between edge localized modes in magnetic fusion devices

Fusion performance has been shown to be related to the H-mode pedestal structure. The pedestal is associated with steep gradients that are the source of free energy for microinstabilities. A variety of instabilities have been shown to co-exist in the pedestal. This paper reviews the experimentally observed signatures of instabilities during the pedestal parameters' evolutions, with a focus on turbulence observations, made in between edge localized modes (ELMs), in multiple tokamaks. ELMs are cyclic events associated with bursty relaxations of the pedestal. The multiple machine results point to very similar pedestal localized modes suggesting the existence of a unifying mechanism governing the modes' onset and dynamics. Modeling efforts to identify the generation mechanism of these instabilities are introduced. Several remaining challenges include the study of the instability saturation mechanisms and impact on the pedestal structure formation. Finally, potential future research avenues will require multiscale-non-linear-gyrokinetic analyses to study the non-linear interaction between these instabilities, and which pedestal conditions facilitate such non-linear interactions, as well as the relationship with the formation of the pedestal through transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Precision gamma-ray constraints for sub-GeV dark matter models

The indirect detection of dark matter particles with mass below the GeV scale has recently received significant attention. Future space-borne gamma-ray telescopes, including All-Sky-ASTROGAM, AMEGO, and GECCO, will probe the MeV gamma-ray sky with unprecedented precision, offering an exciting test of particle dark matter in the MeV-GeV mass range. While it is typically assumed that dark matter annihilates into only one Standard Model final state, this is not the case for realistic dark matter models. Here we analyze existing indirect detection constraints and the discovery reach of future detectors for the well-motivated Higgs and vector-portal models using our publicly-available code Hazma. In particular, we show how to leverage chiral perturbation theory to compute the dark matter self-annihilation cross sections into final states containing mesons, the strongly-interacting Standard Model dynamical degrees of freedom below the GeV scale. We find that future telescopes could probe dark matter self-annihilation cross sections orders of magnitude smaller than those presently constrained by cosmic microwave background, gamma-ray and terrestrial observations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗