Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic parameter”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Water dynamics in C–S–H and M-S-H cement pastes: A revised jump-diffusion and rotation-diffusion model

The Quasi-Elastic Neutron Scattering (QENS) spectra from four cement pastes are re-analyzed by a new revised jump-diffusion and rotation-diffusion model (rJRM). From the QENS fit, it can be seen that the rJRM is an improved model to fit QENS spectra within the whole detected neutron energy transfer and scattering vector. By the rJRM fitting, the structure parameters extracted from QENS spectra show that both the additives aluminum-silicate nanotubes (ASN) and carboxyl group functionalized ASN (ASN-COOH) can improve magnesium-silicate-hydrate (M-S-H) toward calcium-silicate-hydrate (C–S–H) direction in mechanical properties, but the improvement is weakened with decreasing temperatures. The extracted dynamical parameters show that there is a dynamic anomaly near 230 K in not only translational but also rotational diffusion of water confined in all the investigated samples. In conclusion, the anomaly in rotational diffusion is new compared with those results obtained by the QENS fit using other models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Parallel Grid-Forming Inverter-Driven Black Start: Power-Hardware-in-the-Loop Validation: Preprint

With power systems encountering increasing deployment levels of inverter-based resources (IBRs), system restoration using grid-forming (GFM) IBRs has gained attention. Engineered to establish grid voltages in the absence of a stiff grid, black-start-capable GFM IBRs are expected to enhance power system resilience by playing a critical role in bottom-up system restoration. This paper experimentally studies the feasibility of the novel approach through power-hardware-in-the-loop (PHIL) testing and validation with a commercial GFM inverter. The PHIL test setup demonstrates an inverter-driven black start of a 5-MW unbalanced distribution feeder where two GFM inverters collectively black start the feeder: one is a commercial hardware GFM inverter interacting through the PHIL interface, and the other is a software GFM emulated by a real-time simulation with full electromagnetic transient (EMT) models. Both GFM inverters are equipped with negative-sequence voltage control to suppress the voltage imbalance resulting from the unbalanced loading, allowing us to study the dynamic interactions between the GFMs with their control parameters unknown. To evaluate the dynamic behavior of the GFM inverters under the entire black-start process, the EMT model of the distribution system details the transformer and motor dynamics to emulate their inrush and startup behaviors. It abstracts conventional grid-connected inverter dynamics, i.e., grid-following inverters. Oscilloscope measurements of the hardware GFM inverter are also presented for critical steps. Takeaways for further study and field deployment are provided.

grid-forming inverter↗
DOE OSTI2020

RAS Protein Molecular Dynamics Simulations

This dataset contains 6401 files that describe parameters and other inputs for running molecular dynamics computer simulations of Ras proteins associated with a lipid-based cell membrane model in water, as well as resulting time-trajectories of atomic motion in these simulations. Specifically, the content includes: (a) 21 files containing generalized force field parameters, and (b) 22 files of similar format and content for each of 290 separate molecular dynamics simulations (22 ´ 290 = 6,380). The total size of this data set is 1.0 TB.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Parameterized Neural Ordinary Differential Equations: Applications to Computational Physics Problems

This work proposes an extension of neural ordinary differential equations (NODEs) by introducing an additional set of ODE input parameters to NODEs. This extension allows NODEs to learn multiple dynamics specified by the input parameter instances. Our extension is inspired by the concept of parameterized ordinary differential equations, which are widely investigated in computational science and engineering contexts, where characteristics of the governing equations vary over the input parameters. We apply the proposed parameterized NODEs (PNODEs) for learning latent dynamics of complex dynamical processes that arise in computational physics, which is an essential component for enabling rapid numerical simulations for time-critical physics applications. For this, we propose an encoder-decoder-type framework, which models latent dynamics as PNODEs. We demonstrate the effectiveness of PNODEs with important benchmark problems from computational physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

On Projection of Safe Operation for Grid-Following Inverters - Grid Parameter Estimation

This work investigates the projection of safe operation for grid-following inverters using a reference model. This work applies recursive-least square and model reference adaptive estimation techniques to estimate the unknown grid parameters. The estimated grid parameters are used in a dynamic reference model to project the safe operation of the inverter. This work also demonstrates that the controller nonlinearity and parameter variations can cause unsafe operations, such as unexpected instability at low power levels, namely hidden mode of instability. This instability issue in the nominal operating range may occur for given control parameters in a weak-grid condition. The dynamic reference model for PQ-controlled inverters is applied to analyze this instability issue using the estimated grid parameters. Both adaptive identification techniques can estimate the unknown parameters accurately, and the dynamic reference model can safeguard inverters from an unsafe operation, i.e., hidden instability, by examining the incoming new power setpoints before engaging them to the local controller. The findings are experimentally verified using a small-scale two-level 208 V, 5 kVA inverter feeding a 12 kW NHR 9410 power grid emulator.

Hossen, Tareq↗
DOE OSTI2021

Dynamic Probabilistic Safety Assessment Studies for Advanced Reactor Using RAVEN

Probabilistic Safety Assessment (PSA) is used extensively to evaluate the risks associated with complex engineering systems like Nuclear Power Plants (NPPs). Current PSA models are based on the Event-Tree/Fault-Tree (ET/FT) methodology. ET and FT models are static and are based on Boolean logic approaches. In the past, concerns have been raised in the literature regarding the capability of the traditional static modelling approaches to adequately account for the impact of process, hardware, software, firmware and human interactions on the stochastic system behaviour. To overcome the limitations of the traditional approach to PSA, several dynamic PSA methodologies have been proposed. One of the dynamic PSA methodologies used for dynamic evaluations is Dynamic Event Tree (DET) framework which can be used to assess the impact of the parameter variability and scenario dynamics on the PSA model for the initiating event. The DET framework couples the stochastic model (number of component/trains that start on demand, operator action timing, etc.) with a Thermal-Hydraulic (TH) model of the plant. This paper explores the use of DET along with a case study on advanced reactor. The initiating event selected for the study was Class IV power supply failure event. The TH analysis considering uncertainty in various parameters was performed using RELAP5 and Reactor Analysis and Virtual control ENvironment (RAVEN) tool. Based on the uncertainty analysis, it is concluded that the peak clad temperatures (PCT) are within the limits in all the code runs implying a high-degree of safety margin. However, variation in time to reach the PCT was observed among the code runs and the mean time to reach the PCT was found to be around 8590sec (approximately 2.4 hours). Hence, sufficient time margin is available for human intervention and the operator might have a relatively stress-free state during such an accident scenario. Due to the static nature of the traditional PSA models, the safety margin available was lesser, whereas, with the help of dynamic PSA models, one can demonstrate that the actual available safety margin is more in the present case study and is valuable input from the design point of view.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

DART-PFLOTRAN: An ensemble-based data assimilation system for estimating subsurface flow and transport model parameters

Ensemble-based Data Assimilation (EDA), based on the Monte Carlo approach, has been effectively applied to estimate model parameters through inverse modeling in subsurface flow and transport problems. However, implementation of EDA approach involves a complicated workflow that include setting up and executing ensemble forward model simulations, processing observations and model simulation results for parameter updates, and repeat for sequential or iterative EDA. To facilitate the management of such workflow and lower the barriers for adopting EDA-based parameter estimation in subsurface science, we develop a generic software frame-work linking the Data Assimilation Research Testbed (DART) with a massively parallel subsurface FLOw and TRANsport code PFLOTRAN. The new DART-PFLOTRAN leverages both the core data assimilation engines in DART and the computational power afforded by PFLOTRAN. In addition to the standard smoother and filtering options, DART-PFLOTRAN enables an iterative EDA workflow based on the Ensemble Smoother for Multiple Data Assimilation method (ES-MDA) to improve estimation accuracy for nonlinear forward problems. Here, we verify the implementation of ES-MDA in DART-PFLOTRAN using two synthetic cases designed to estimate static permeability and dynamic exchange fluxes across the riverbed, respectively, from continuous temperature measurements made across a depth profile. One-dimensional hydro-thermal simulations are performed in both cases to relate temperature responses with the parameters of interest. In the case of estimating dynamic parameters, we demonstrate the flexibility of DART-PFLOTRAN in automating sequential ES-MDA workflow, which will significantly reduce the time researchers spend on managing complex workflows in similar applications. Both studies yield accurate estimations of the parameters compared to their synthetic truth, while ES-MDA leads to more accurate estimation when a high level of nonlinearity exist between observed responses and unknown parameters. With a code base in Python and Fortran, DART-PFLOTRAN paves the way for applications in large-scale subsurface inverse modeling by automating the complex workflow of sequential ES-MDA that can be executed on various computing platforms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Dynamic Model Reduction for Large-Scale Power Systems Using Wide-Area Measurements

To perform faster than the real-time dynamic simulation of large-scale power systems, it is necessary to reduce the simulated system size by using equivalents for surrounding areas of the study area, and existing dynamic model reduction approach could provide the needed structure of the reduced area. However, further parameter optimization is required to achieve the desired accuracy. In this paper, a particle swarm optimization (PSO) based approach is used to solve the above problem. Parameters for the individual dynamic elements in the reduced system are calibrated repeatedly until the wide-area measurements of the reduced model and the original model are very similar to each other with satisfactory accuracy. Results indicate that after optimization, the dynamic response of the reduced model matches better with that of the original one than using existing methods. Under both the generator-trip event and the bus-fault event, the reduced model has a higher frequency match and less power mismatch.

42 ENGINEERING↗
DOE OSTI2022

The influence of random packed column parameters on the liquid holdup and interfacial area

Abstract Carbon dioxide capture via solvent absorption in packed columns has emerged as a potential technology to mitigate coal‐fired power plant CO 2 emissions. Parameters, including packing types, solvent properties, and operating conditions, could potentially affect the packed column CO 2 capture efficiency. To understand the importance of those parameters and help packed column optimization, a design of experiments (DoEs) method was proposed to generate input parameter matrix. Combined with multiphase computational fluid dynamics (CFD), the random packed column parameter influence on the liquid holdup and interfacial area can be efficiently investigated. Surrogate‐based sensitivity analysis shows that the solvent flow rate and contact angle are key factors dictating liquid holdup and interfacial area. On the other hand, solvent viscosity has a marginal impact on the interfacial area. The sensitivity scores were calculated for each input parameter to guide the selection of dimensionless numbers for the liquid holdup and interfacial area correlation development.

42 ENGINEERING↗
DOE OSTI2022

Distinguishing Bulk and Surface Recombination By Measuring and Modeling Ultrafast Carrier Dynamics

Distinguishing Bulk and Surface Recombination By Measuring and Modeling Ultrafast Carrier Dynamics The combination of modeling with TRTS and TRPL under multiple excitation conditions improves accuracy and precision of parameter estimation. Ultrafast carrier dynamics can link materials, processing, and device performance to direct the design of more efficient solar cells.

14 SOLAR ENERGY↗
DOE OSTI2021

Toward nonthermal control of excited quantum materials: framework and investigations by ultrafast electron scattering and imaging

Quantum material systems upon applying ultrashort laser pulses provide a rich platform to access excited material phases and their transformations that are not entirely like their equilibrium counterparts. The addressability and potential controls of metastable or long-trapped out-of-equilibrium phases have motivated interests both for the purposes of understanding the nonequilibrium physics and advancing the quantum technologies. Thus far, the dynamical spectroscopic probes eminently focus on microscopic electronic and phonon responses. For characterizing the long-range dynamics, such as order parameter fields and fluctuation effects, the ultrafast scattering probes offer direct sensitivity. Bridging the connections between the microscopic dynamics and macroscopic responses is central toward establishing the nonequilibrium physics behind the light-induced phases. Here, we present a path toward such understanding by cross-examining the structure factors associated with different dynamical states obtained from ultrafast electron scattering, imaging, and modeling. We give the basic theoretical framework on describing the non-equilibrium scattering problems and briefly describe how such framework relates to the out-of-equilibrium phenomena. We give effective models outlining the emergences of nonthermal critical points, hidden phases, and non-equilibrium relaxational responses from vacuum-suspended rare-earth tritellurides, tantalum disulfides thin films, and vanadium dioxide nanocrystalline materials upon light excitations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

Here we analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

TDCOSMO - IV. Hierarchical time-delay cosmography – joint inference of the Hubble constant and galaxy density profiles

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H0 = 73.3−1.8+1.7 km s−1 Mpc−1, describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H0. We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H0. Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H0 = 74.5−6.1+5.6 km s−1 Mpc−1. Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H0 = 67.4−3.2+4.1 km s−1 Mpc−1. This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H0. While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data.Key words: gravitational lensing: strong / galaxies: general / galaxies: kinematics and dynamics / distance scale / cosmological parameters / cosmology: observations⋆ The full analysis is available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Parametrized Hamiltonian simulation using quantum optimal control

Analog quantum simulation offers a hardware-specific approach to studying quantum dynamics, but mapping a model Hamiltonian onto the available device parameters requires matching the hardware dynamics. In this work, we introduce a paradigm for quantum Hamiltonian simulation that leverages digital decomposition techniques and optimal control to perform analog simulation. We validate this approach by constructing the optimal analog controls for a superconducting transmon device to emulate the dynamics of an extended Bose-Hubbard model. We demonstrate the role of control time, digital error, and pulse complexity, and we explore the accuracy and robustness of these controls. We conclude by discussing the opportunity for implementing this paradigm in near-term quantum devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spherically expanding flame in silane–hydrogen–nitrous oxide–argon mixtures

The effect of silane addition on the laminar flame speed ($S^0_u$) of flames propagating in hydrogen–nitrous oxide–argon mixtures has been investigated experimentally for the first time using the spherically expanding flame technique in a constant volume combustion chamber. Replacing hydrogen by silane and maintaining the equivalence ratio constant, much higher flame speeds, explosion peak pressures, and pressure rise coefficients were measured. A previously developed detailed reaction model has been updated based on ab initio thermodynamic properties calculations and collision limit violation analysis. The improved reaction model demonstrates encouraging performance in predicting the flame speed, with a mean absolute error below 11%. To explain the effect of silane addition on the flame dynamics, a number of parameters have been calculated including OH and H rate of production, heat release rate per reaction, and sensitivity coefficient on . The dynamics of freely propagating flames in SiH 4 –H 2 –N 2 O–Ar mixtures is essentially controlled by reactions of the H–O–N chemical system: N 2 O+H=N 2 +OH, OH+H 2 =H 2 O+H, and N 2 O(+M)=N 2 +O(+M). Whereas silane addition does not influence much the rate of production of OH, it significantly modifies that of H with a number of pyrolytic chemical pathways of silicon hydrides, such as SiH+H 2 =SiH 2 +H and Si+H 2 =SiH+H, which act as sink of H atom as they proceed in the backward direction. The reactions forming SiO(s) and SiO 2 (s), such as SiO+OH=SiO 2 (s)+H and 2SiO=2SiO(s), are exothermic and significantly contribute to the temperature increase. The adiabatic, constant pressure flame temperature for mixture containing silane is significantly higher, up to several 100’s K. The increase of induced by silane addition seems to be mostly related to the large increase of the flame temperature which leads to higher energy release rate.

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

The effect of hydroxyl spacing in diols on the solvation structure, dynamics, and transport properties of choline chloride-based deep eutectic solvents

Deep eutectic solvents (DESs) are a class of liquids that offer great potential in alleviating some of the challenges present in today's long-term energy storage methods because they have physical properties that are favorable for storable electrolyte solutions. In this work, a series of glycols (ethylene glycol, 1,3-propanediol, 1,4-butanediol, and 1,5-pentanediol) were studied as potential hydrogen bond donors (HBD) with a common choline chloride (ChCl) as the hydrogen bond acceptor (HBA). The solvation dynamics of the prepared systems were studied by measuring the solvent reorganization response using femtosecond transient absorption spectroscopy (fs-TA). Conductivity, viscosity, density, ET(30) polarity, and dynamics of the prepared DESs were analyzed, with a particular interest in determining the effect of HBD chain length on these parameters. Here, classical molecular dynamics simulations were employed to investigate how the local liquid structure, solvent dynamics, and bulk solvent properties vary with changes in glycol chain length.

Conductivity↗
DOE OSTI2024

An experimentally validated electro-thermal EV battery pack model incorporating cycle-life aging and cell-to-cell variations

Lithium-ion batteries are used in a wide variety of applications. To meet the power and energy demands of these applications battery packs are composed of hundreds to thousands of cells. The electrical and thermal interactions between cells introduce additional complexity in the pack dynamics. To capture these effects, a battery pack model composed of 192 cells based on a first-generation (2012) Nissan Leaf battery pack is developed in MATLAB/Simulink/Simscape. Here, with this model, we simulate the electrical dynamics (using a first-order equivalent-circuit model), the thermal dynamics (using a first-order lumped-parameter thermal model), and the aging dynamics (using a semi-empirical severity factor-based model) of every cell in the pack and we also create a pack thermal model that explicitly captures the heat exchange between the modules, and the cells contained within, during operation. The models are calibrated and validated, both at the cell and pack level, with experimental data. Two different case studies of this pack model are investigated. In the first case study, an initial, normally-distributed, cell-to-cell capacity variation is introduced and its effect on the pack voltage and module temperatures is studied. In the second case study, we deliberately insert cells with lower than nominal capacity into the pack and we investigate how this type of initial cell-to-cell capacity variation affects the pack’s ability to deliver energy over time. Finally, we also study how parallel-connected cells can reduce the effects of cell-to-cell variations at the expense of increased aging of the pack overall.

25 ENERGY STORAGE↗
DOE OSTI2024

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗