Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic system”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Effect of Temperature on Thrombogenicity Testing of Biomaterials in an In Vitro Dynamic Flow Loop System

To develop and standardize a reliable in vitro dynamic thrombogenicity test protocol, the key test parameters that could impact thrombus formation need to be investigated and understood. In this study, we evaluated the effect of temperature on the thrombogenic responses (thrombus surface coverage, thrombus weight, and platelet count reduction) of various materials using an in vitro blood flow loop test system. Whole blood from live sheep and cow donors was used to assess four materials with varying thrombogenic potentials: negative-control polytetrafluoroethylene (PTFE), positive-control latex, silicone, and high-density polyethylene (HDPE). Blood, heparinized to a donor-specific concentration, was recirculated through a polyvinyl chloride tubing loop containing the test material at room temperature (22–24°C) for 1 hour, or at 37°C for 1 or 2 hours. The flow loop system could effectively differentiate a thrombogenic material (latex) from the other materials for both test temperatures and blood species ( p < 0.05). However, compared with 37°C, testing at room temperature appeared to have slightly better sensitivity in differentiating silicone (intermediate thrombogenic potential) from the relatively thromboresistant materials (PTFE and HDPE, p < 0.05). These data suggest that testing at room temperature may be a viable option for dynamic thrombogenicity assessment of biomaterials and medical devices.

Engineering↗
DOE OSTI2023

Density Functional Theory Guided Investigation of Ligand‐Induced Neptunyl‐Neptunyl Interactions

Abstract Actinyl‐actinyl interactions are particularly prevalent for the pentavalent neptunyl cation (Np(V)O 2 ) + where these interactions appear either as a T ‐ or D ‐shape (diamond‐shape). T ‐shaped interactions have been previously identified in high concentration Np(V) solutions containing simple anions (NO 3 − , ClO 4 − , Cl − ) whereas D ‐shaped have only been isolated in the solid‐state in the presence of carboxylate ligands. In this study, Density Functional Theory (DFT) calculations were paired with Raman spectroscopy to evaluate the formation of D ‐shaped interactions in the presence of aliphatic (R=H (formate), CH 3 (acetate), CH 2 CH 3 (propionate)) and aromatic (R=C 6 H 5 (benzoate), C 6 H 4 OH (4‐hydroxybenzoate), C 5 H 4 N (isonicotinate)) carboxylate ligands. DFT studies indicate that the Δ G to form hydrated T ‐ and D ‐shaped forms are not spontaneous but become so with the addition of the carboxylate ligands. Raman spectra of the Np(V) carboxylate solutions contained vibrational modes associated with the D ‐shaped interactions, but spectral changes observed over time indicate a dynamic system. Crystallization experiments from the Np(V) carboxylate systems confirmed the presence of D‐ shaped dimers for the aromatic carboxylates, suggesting that the choice of the anion in solution favors actinyl‐actinyl interactions even at low concentrations (≤20 mM) of Np(V).

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

Development of data‐driven exponential integrators with application to modeling of delay photocurrents

Abstract Radiation‐induced photocurrent effects represent a threat to microelectronic components operating in space, manmade and terrestrial radiation environments. Analysis of these threats by circuit simulations requires accurate and computationally efficient compact models. Most existing compact models are based on closed form analytic solutions of the governing equations and require empirical assumptions and idealizations that can limit their validity. In this paper we formulate an alternative numerical, data‐driven approach that learns a compact model from data representative of the type of measurements one can obtain in an experimental facility. To develop the model we start from a generic discrete‐time dynamical system and then use physics knowledge to refine its structure. Numerical studies demonstrate the potential of the model and establish some empirical guidelines for its training.

Bochev, Pavel↗
DOE OSTI2022

Dispersion-enhanced sequential batch sampling for adaptive contour estimation

In computer simulation and optimal design, sequential batch sampling offers an appealing way to iteratively stipulate optimal sampling points based upon existing selections and efficiently construct surrogate modeling. Nonetheless, the issue of near duplicates poses tremendous quandary for sequential learning. It refers to the situation that selected critical points cluster together in each sampling batch, which are individually but not collectively informative towards the optimal design. Near duplicates severely diminish the computational efficiency as they barely contribute extra information towards update of the surrogate. To address this issue, we impose a dispersion criterion on concurrent selection of sampling points, which essentially forces a sparse distribution of critical points in each batch, and demonstrate the effectiveness of this approach in adaptive contour estimation. Specifically, we adopt Gaussian process surrogate to emulate the simulator, acquire variance reduction of the critical region from new sampling points as a dispersion criterion, and combine it with the modified expected improvement (EI) function for critical batch selection. The critical region here is the proximity of the contour of interest. This proposed approach is vindicated in numerical examples of a two-dimensional four-branch function, a four-dimensional function with a disjoint contour of interest and a time-delay dynamic system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Trusted Simulation: Considering Model Quality in the Context of User Trust

A high‐quality simulation model should help its users to easily and appropriately calibrate their trust in the model. Traditional evaluation metrics such as validation and robustness are necessary but insufficient for this task. Trust calibration depends on factors like the model's transparency, applicability to intended use, usability, reputation, and consideration of potential bias. This article proposes a framework for designing and evaluating system dynamics models by considering factors that contribute to the proper calibration of user trust. This framework takes inspiration from trusted artificial intelligence, broadening our traditional concept of model quality and explicitly focusing on what users need to consider a model trustworthy and to understand the model's relevance to its intended purpose. The trusted simulation framework can improve our integration of model quality activities throughout the modeling process, leading to more impactful and better‐targeted model design, development, and evaluation.

Naugle, Asmeret Bier [Sandia National Laboratories↗
DOE OSTI2023

Robust multi‐loop control of a floating wind turbine

Abstract A principal challenge facing the control of floating offshore wind turbines (FOWTs) is the problem of instability, or “negative damping,” when using blade pitch feedback to control generator speed. This closed‐loop instability can be attributed to non‐minimum phase zeros in the transfer function from blade pitch to generator speed. Standard approaches to improving stability and performance include robust tuning of control gains and introducing multiple feedback loops to respond to platform motion. Combining these approaches is nontrivial because multiple control loops complicate the impact of coupling in the system dynamics. The single‐loop approach to analyzing stability robustness neglects inter‐loop coupling, while a simplistic multi‐loop approach is highly sensitive to dimensional scaling and overestimates the robustness of the single‐loop controller. This work proposes a sensitivity representation that separates some of the natural FOWT dynamic coupling into a parallel feedback loop in the sensitivity function loop to address both of these concerns. The modified robustness measure is used with a simplified linear FOWT model to optimize scheduled multi‐loop control parameters in an automated tuning procedure. This controller is implemented for the 10‐MW Ultraflexible Smart FLoating Offshore Wind Turbine (USFLOWT) and compared against conventional single‐ and multi‐loop controllers tuned using frequency‐domain analysis and high‐fidelity OpenFAST simulations. The multi‐loop robust controller shows the highest overall performance in generator speed regulation and tower load reduction, though consideration of power quality, actuator usage, and other structural loading leads to additional trade‐offs.

Stockhouse, David↗
DOE OSTI2025

Quantum computation of SU(2) lattice gauge theory with continuous variables

We present a quantum computational framework for pure SU(2) lattice gauge theory, using continuous variables instead of discrete qubits to represent the infinite-dimensional Hilbert space of the gauge fields. We consider a ladder as well as a two-dimensional grid of plaquettes, detailing the use of gauge fixing to reduce the degrees of freedom and simplify the Hamiltonian. We demonstrate how system dynamics, ground states, and energy gaps can be computed using the continuous-variable approach to quantum computing. Our results indicate that it is feasible to study non-Abelian gauge theories with continuous variables, providing new avenues for understanding the real-time dynamics of quantum field theories.

Gauge Symmetry↗
DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗
DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗
DOE OSTI2021

The MBARI-WEC: a power source for ocean sensing

Abstract Interest in wave energy converters to provide autonomous power to various ocean-bound systems, such as autonomous underwater vehicles, sensor systems, and even aquaculture farms, has grown in recent years. The Monterey Bay Aquarium Research Institute has developed and deployed a small two-body point absorber wave energy device suitable to such needs. This paper provides a description of the system to support future open-source access to the device and further the general development of similar wave energy systems. Additionally, to support future control design and system modification efforts, a set of hydrodynamic models are presented and cross-compared. To test the viability of using a linear frequency-domain admittance model for controller tuning, the linear model is compared against four WEC-Sim models of increasing complexity. The linear frequency-domain model is found to be generally adequate for capturing system dynamics, as the model agreement is good and the degree of nonlinearity introduced in the WEC-Sim models is generally less than 2.5%.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

From PINNs to PIKANs: recent advances in physics-informed machine learning

Physics-Informed Neural Networks (PINNs) have emerged as a key tool in Scientific Machine Learning since their introduction in 2017, enabling the efficient solution of ordinary and partial differential equations using sparse measurements. Over the past few years, significant advancements have been made in the training and optimization of PINNs, covering aspects such as network architectures, adaptive refinement, domain decomposition, and the use of adaptive weights and activation functions. A notable recent development is the Physics-Informed Kolmogorov-Arnold Networks (PIKANS), which leverage a representation model originally proposed by Kolmogorov in 1957, offering a promising alternative to traditional PINNs. In this review, we provide a comprehensive overview of the latest advancements in PINNs, focusing on improvements in network design, feature expansion, optimization techniques, uncertainty quantification, and theoretical insights. We also survey key applications across a range of fields, including biomedicine, fluid and solid mechanics, geophysics, dynamical systems, heat transfer, chemical engineering, and beyond. Lastly, we review computational frameworks and software tools developed by both academia and industry to support PINN research and applications.

Kolmogorov-Arnold networks↗
DOE OSTI2022

Conceptual core design of the molten metal fueled microreactor with self-regulating capability

The microreactor concept has received significant attention in the United States for its lower capital investment, siting flexibility and high mobility. A microreactor aims to provide reliable electricity in remote communities, mining areas, or isolated islands, which will save the high fuel transportation cost. Thus, minimizing the staff level for operating and maintaining a microreactor is desirable. The objective of the study is to design a self-regulated Molten Metal fueled microReactor core (MMR) and to demonstrate the self-regulating capability of the MMR core. The control of the reactor is solely relying on reactivity feedbacks from fuel temperature perturbation. Molten UMn (liquid form at the temperature over 720°C) was utilized as the fuel material in MMR due to the larger thermal expansion coefficient compared to the conventional solid fuel. The MMR core has a fast neutron spectrum, which corresponds to a small burnup reactivity swing. The core design parameters were selected based on DAKOTA, REBUS, and ANLHTP calculations. The MMR is designed to run for 10 years at a power level of 15 MWth. Furthermore, no refueling is needed during the reactor lifetime. The system dynamic analyses were performed, and it was found that the reactor can be self-regulating within a temperature range of 800 ± 80°C for a loss of heat sink transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Impact of climate change on water availability and its propagation through the Western U.S. power grid

Climate impact assessment for power systems has tended to focus on the effects of water availability on generation without explicitly considering how these impacts propagate through the interconnected grid to affect generation in other locations. We investigate the impact of projected water availability on the operations of the contemporary Western U.S. power grid, focusing specifically on regional interdependencies. Hydrologic simulations derived from three GCMs (CCSM4, INMCM4 and GFDL-CM3) and two radiative scenarios (RCP4.5 and RCP8.5) are used to force a large-scale, distributed water management model (MOSART-WM), which translates water availability into power generation constraints at hydropower plants and water-dependent thermo-electric plants. Power system dynamics are evaluated using the production cost model PLEXOS. We find that the inter-regional connections across the Western U.S. electricity infrastructure are key in managing drastic variations in regional water availability. Net generation in the Northwest and Northern California is substantially sensitive to local water availability. Generation from the Desert Southwest emerges as a critical resource to compensate for variations in water availability in those regions. Climate change impacts on Northwestern water availability therefore drive future changes in other regions’ generation and regional power flows. The large inter-annual hydrological variability of Northern California exacerbates grid stress arising from long-term changes in water availability in the Northwest. While the regional power flow directions seem insensitive to long-term variations in water availability, the analysis highlights the need to consider other compounding regional factors such as changes in Southern California net load or changes in regional fuel prices.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Plutonium mobilization from contaminated estuarine sediments, Esk Estuary (UK)

Since 1952, liquid radioactive effluent containing 238-242 Pu, 241 Am, 237 Np, 137 Cs, and 99 Tc has been released with authorization from the Sellafield nuclear complex (UK) into the Irish Sea. This represents the largest source of plutonium (Pu) discharged in all western Europe, with 276 kg having been released. In the Eastern Irish Sea, the majority of the transuranic activity has settled into an area of sediments (Mudpatch) located off the Cumbrian coast. Radionuclides from the Mudpatch have been re-dispersed via particulate transport in fine-grained estuarine and intertidal sediments to the North-East Irish Sea, including the intertidal saltmarsh located at the mouth of the Esk Estuary. Saltmarshes are highly dynamic systems which are vulnerable to external agents (sea level change, erosion, sediment supply, and freshwater inputs), and their stability remains uncertain under current sea level rise projections and possible increases in storm activity. In this work, we examined factors affecting Pu mobility in contaminated sediments collected from the Esk Estuary by conducting leaching experiments under both anoxic and oxic conditions. Here, leaching experiments were conducted over a 9-month period and were periodically sampled to determine solution phase Pu via multicollector-inductively coupled plasma-mass spectrometry (MC-ICP-MS), and to measure redox indicators (Eh, pH and extractable Fe(II)). Microbial community composition was also characterized in the sediments, and at the beginning and end of the anoxic/oxic experiments. Results show that: 1) Pu leaching is about three times greater in solutions leached under anoxic conditions compared to oxic conditions, 2) the sediment slurry microbial communities shift as conditions change from anoxic to oxic, 3) Pu leaching is enhanced in the shallow sediments (0–10 cm depth), and 4) the magnitude of Pu leached from sediments is not correlated with total Pu, indicating that the biogeochemistry of sediment-associated Pu is spatially heterogeneous. These findings provide constraints on the stability of redox sensitive Pu in biogeochemically dynamic/transient environments on a timescale of months and suggests that anoxic conditions can enhance Pu mobility in estuarine systems.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Exponential time differencing for the tracer equations appearing in primitive equation ocean models

The tracer equations are part of the primitive equations used in ocean modeling and describe the transport of tracers, such as temperature, salinity or chemicals, in the ocean. Depending on the number of tracers considered, several equations may be added to and coupled to the dynamics system. In many relevant situations, the time-step requirements of explicit methods imposed by the transport and mixing in the vertical direction are more restrictive than those for the horizontal, and this may cause the need to use very small time steps if a fully explicit method is employed. To overcome this issue, we propose an exponential time differencing (ETD) solver where the vertical terms (transport and diffusion) are treated with a matrix exponential, whereas the horizontal terms are dealt with in an explicit way. In this work, we investigate numerically the computational speed-ups that can be obtained over other semi-implicit methods, and we analyze the advantages of the method in the case of multiple tracers.

42 ENGINEERING↗
DOE OSTI2022

A comparison of neural network architectures for data-driven reduced-order modeling

The popularity of deep convolutional autoencoders (CAEs) has engendered new and effective reduced-order models (ROMs) for the simulation of large-scale dynamical systems. Despite this, it is still unknown whether deep CAEs provide superior performance over established linear techniques or other network-based methods in all modeling scenarios. To elucidate this, the effect of autoencoder architecture on its associated ROM is studied through the comparison of deep CAEs against two alternatives: a simple fully connected autoencoder, and a novel graph convolutional autoencoder. Through benchmark experiments, it is shown that the superior autoencoder architecture for a given ROM application is highly dependent on the size of the latent space and the structure of the snapshot data, with the proposed architecture demonstrating benefits on data with irregular connectivity when the latent space is sufficiently large.

42 ENGINEERING↗
DOE OSTI2022

Operator inference for non-intrusive model reduction with quadratic manifolds

This paper proposes a novel approach for learning a data-driven quadratic manifold from high-dimensional data, then employing this quadratic manifold to derive efficient physics-based reduced-order models. The key ingredient of the approach is a polynomial mapping between high-dimensional states and a low-dimensional embedding. This mapping consists of two parts: a representation in a linear subspace (computed in this work using the proper orthogonal decomposition) and a quadratic component. The approach can be viewed as a form of data-driven closure modeling, since the quadratic component introduces directions into the approximation that lie in the orthogonal complement of the linear subspace, but without introducing any additional degrees of freedom to the low-dimensional representation. Combining the quadratic manifold approximation with the operator inference method for projection-based model reduction leads to a scalable non-intrusive approach for learning reduced-order models of dynamical systems. Further, applying the new approach to transport-dominated systems of partial differential equations illustrates the gains in efficiency that can be achieved over approximation in a linear subspace.

42 ENGINEERING↗