Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Performance 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 1,495 records · Page 83

End-to-end protocol for high-quality quantum approximate optimization algorithm parameters with few shots

The quantum approximate optimization algorithm (QAOA) is a quantum heuristic for combinatorial optimization that has been demonstrated to scale better than state-of-the-art classical solvers for some problems. For a given problem instance, QAOA performance depends crucially on the choice of the parameters. While average-case optimal parameters are available in many cases, meaningful performance gains can be obtained by fine-tuning these parameters for a given instance. This task is especially challenging, however, when the number of circuit executions (shots) is limited. In this work, we develop an end-to-end protocol that combines multiple parameter settings and fine-tuning techniques. We use large-scale numerical experiments to optimize the protocol for the shot-limited setting and observe that optimizers with the simplest internal model (linear) perform best. We implement the optimized pipeline on a trapped-ion processor using up to 32 qubits and 5 QAOA layers, and we demonstrate that the pipeline is robust to small amounts of hardware noise. To the best of our knowledge, these are the largest demonstrations of QAOA parameter fine-tuning on a trapped-ion processor in terms of two-qubit gate count.

quantum algorithms & computation↗

Magnetic excitations of the hybrid multiferroic (ND 4 ) 2 FeCl 5 ·D 2 O

We report a comprehensive inelastic neutron scattering study of the hybrid molecule-based multiferroic compound ( ND 4 ) 2 FeCl 5 · D 2 O in the zero-field incommensurate cycloidal phase and the high-field quasicollinear phase. The spontaneous electric polarization changes its direction concurrently with the field-induced magnetic transition, from mostly aligned with the crystallographic a axis to the c axis. To account for such a change in polarization direction, the underlying multiferroic mechanism was proposed to switch from the spin-current model induced via the inverse Dzyaloshinskii-Moriya interaction to the p - d hybridization model. We perform a detailed analysis of the inelastic neutron data of ( ND 4 ) 2 FeCl 5 · D 2 O using linear spin-wave theory to quantify magnetic interaction strengths and investigate the possible impact of different multiferroic mechanisms on the magnetic couplings. Overall our result reveals that the spin dynamics of both multiferroic phases can be well described by a Heisenberg Hamiltonian with easy-plane anisotropy. We do not find notable differences between the optimal model parameters of the two phases. The hierarchy of exchange couplings and the balance among frustrated interactions remain the same between two phases, suggesting that magnetic interactions in ( ND 4 ) 2 FeCl 5 · D 2 O are much more robust than the electric polarization in response to delicate reorganizations of the electronic degrees of freedom in an applied magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Transactive Electric Water Heater Agent: Design and Performance Evaluation

Electric water heaters (EWHs) are usually equipped with inbuilt thermostats to measure water temperature near the installed positions. Most of the existing EWHs have only two thermostats and often do not have water flow sensor installed. Estimation of state of heat energy (SOHE) inside the hot water tank, which is crucial for efficient control and operation of the EWH, under such imperfect system conditions is very challenging. This paper, therefore, designs a transactive EWH agent (TEWHA) to better estimate SOHE in presence of measurement errors and imperfect system knowledge. The TEWHA, built-in Python, approximates the EWH physics using limited measurements obtained from the ground truth GridLABD model. The performance of the TEWHA is validated against ground truth system for three scenarios: a) ideal condition where TEWHA has complete knowledge of the system, b) imperfect condition where TEWHA has incomplete system information, and c) noised condition where TEWHA measurements contain random errors. Finally, an optimization problem is formulated and solved to demonstrate how TEWHA could use the estimated and measured parameters for a transactive control of EWH. Index Terms—Demand flexibility

demand flexibility, distribution system, electric ↗

1.1 eV GaInAs Cell Development for Dual-Use Solar and 1070 nm Laser Power Converters

Dual-use power converter cells can receive solar and laser power simultaneously to generate current, an application relevant to space and terrestrial industries. This study investigates two concepts of solar cells optimized for dual-use 1070 nm laser and solar power conversion, a single-junction and triple-junction cell. Because it is bandgap-tuned for a 1070 nm laser, the 1.1 eV junction is incorporated into each of the two designs, making its development key to the success of both concepts. Efficiency data for a one junction GaInAs cell demonstrates a laser conversion efficiency of 38% at 1070nm wavelength without an anti-reflection coating. However, the single-junction device requires optimization to reduce short wavelength absorption of the broad solar spectrum. In both devices, the graded buffer layers in the GaInAs cell affects the cell's performance by reducing threading dislocations in the active junction. However, the buffer in the three-junction device also acts as a lateral transport layer and so affects the fill factor depending on its sheet resistance. By varying the buffer thickness, we demonstrate a direct relationship between buffer thickness and sheet resistance reduction, while considering implications to open-circuit voltages. We also performed resistance modeling to determine the optimal grid spacing and thickness of the grid fingers to minimize losses due to sheet resistance and grid shading.

junctions↗

Resilience Analytics for Interdependent Power and Water Distribution Systems

The occurrence of large-scale outages in power distribution systems (PDSs) caused, for instance, by extreme weather events has raised concerns on the cascading impact of the power outages on the resilience of other critical infrastructures, including water distribution systems (WDSs), whose healthy operation highly relies on the availability of power from PDSs. Here, this paper proposes an analytics model for quantifying the interdependence between the resilience of power and water distribution systems. The proposed model first performs spatio-temporal outage analysis on the PDS and WDS to determine the PDS power serving capability and the resulting cascading impact on the WDS water serving capability during the occurrence of disruptions. The resilience interdependence is then determined by a proposed set of six metrics calculated using the power and water serving capability profiles, which trace and quantify the cascading impacts of power outages in the operation of WDSs during the degradation and recovery stages. The proposed analytics model is implemented on the IEEE 33-bus power distribution system supplying power for pumping stations of a 15-node test WDS, considering multiple cases of water storage capacity, single and multiple outage scenarios, and availability of distributed generation. The numerical results demonstrate that the proposed metrics track the temporal and spatial interdependence between the resilience of power and water distribution systems, providing a mechanism to identify weak spots, and assess options to enhance the resilience of the critical water infrastructure.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Pyrolysis of Cyclohexane and 1-Hexene at High Temperatures and Pressures—A Photoionization Mass Spectrometry Study

Cycloalkanes are important components of a wide range of fuels. However, there are few experimental data at simultaneously high temperatures and pressures similar to those found in practical systems. Such data are necessary for developing and testing chemical kinetic models. In this study, data relevant to cycloalkane pyrolysis were obtained from high repetition rate shock tube experiments coupled with synchrotron-based photoionization mass spectrometry diagnostics. The pyrolysis of cyclohexane was studied over 1270–1550 K and ~9 bar, while the more reactive primary decomposition product, 1-hexene, was studied at 1160–1470 K and ~5 bar. Insights into the decomposition of the parent molecules, the formation of primary products and the production of aromatic species were gained. Simulations were performed with models for cyclohexane and 1-hexene that were based on literature models. The results indicate that over several hundred microseconds reaction time at high pressures and temperatures the pyrolysis of cyclohexane is largely dominated by reactions initiated by cyclohexyl radicals. Furthermore, good agreement between the simulations and the experiments were observed for cyclohexane and 1-hexene with a modified version of the cyclohexane model. Conversely, the 1-hexene model did not reproduce the experimental observations.

1-alkenes↗

Experimental evaluation of thermal-fluids performance of Helium-Cooled Flat Plate (HCFP) divertor

Various helium (He)-cooled solid-tungsten (W) divertor concepts have been proposed for long-pulse magnetic fusion energy reactors. Among these concepts, the He-cooled flat plate divertor (HCFP) modules have the largest plasma-facing surface area of ~0.2 m 2 . Simulations have shown that the most recent version of the design can withstand heat fluxes as great as 8 MW/m 2 . Earlier experimental studies of a single shortened HCFP cooling unit with a slot length of 7.6 cm used air at ambient temperature and pressures below 0.6 MPa. Here, we present initial experimental studies of a copper alloy and steel test section modeling a shortened HCFP cooling unit using He at prototypical pressure of 10 MPa, inlet temperatures T i ≤ 200 °C and steady-state incident heat fluxes q" ≤ 1.2 MW/m 2 . Results for Nusselt number Nu as a function of Reynolds number Re were obtained for Re = 1.2×10 4 –3.4×10 4 , and used to develop a Nu(Re) correlation and validate numerical models of the test section using commercial computational fluid dynamics (CFD) software. Simulations with this model are performed to evaluate the effect of the shortened slot. Furthermore, these analyses are used to estimate the thermal-fluids performance of the HCFP under prototypical conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fan and Pump Efficiency in Modelica based on the Euler Number

Simulation programs often assume constant hydraulic efficiency for fan or pump models when performance curves are unavailable. This is inaccurate because the hydraulic efficiency varies with the operation condition. It therefore consistently underestimates the power draw at off-design conditions at which the hydraulic efficiency drops. Use of a modified Euler number allows computing the hydraulic efficiency and shaft power with limited data. Others showed the validity of the modified Euler number for fan efficiency calculations. We show that it is also applicable for pumps, and present its implementation in Modelica for a fan or pump model. Here, the only input required from the user is one data point at which the hydraulic efficiency is at its maximum. The reported method is applicable regardless of the type, size, or operational region of the fan or pump. Across a sample of eighteen sets of pump data and seven sets of fan data, the errors of the computed power from interpolated data were within 15% for the range of 20% - 70% of maximum flow rate and 40% - 90% of maximum pressure rise, excluding outliers.

component model↗

Reduced‐Order Probabilistic Emulation of Physics‐Based Ring Current Models: Application to RAM‐SCB Particle Flux

Abstract In this work, we address the computational challenge of large‐scale physics‐based simulation models for the ring current. Reduced computational cost allows for significantly faster than real‐time forecasting, enhancing our ability to predict and respond to dynamic changes in the ring current, valuable for space weather monitoring and mitigation efforts. Additionally, it can also be used for a comprehensive investigation of the system. Thus, we aim to create an emulator for the Ring current‐Atmosphere interactions Model with Self‐Consistent magnetic field (RAM‐SCB) particle flux that not only improves efficiency but also facilitates forecasting with reliable estimates of prediction uncertainties. The probabilistic emulator is built upon the methodology developed by Licata and Mehta (2023), https://doi.org/10.1029/2022sw003345 . A novel discrete sampling is used to identify 30 simulation periods over 20 years of solar and geomagnetic activity. Focusing on a subset of particle flux, we use Principal Component Analysis for dimensionality reduction and Long Short‐Term Memory (LSTM) neural networks to perform dynamic modeling. Hyperparameter space was explored extensively resulting in about 5% median symmetric accuracy across all data sets for one‐step dynamic prediction. Using a hierarchical ensemble of LSTMs, we have developed a reduced‐order probabilistic emulator (ROPE) tailored for time‐series forecasting of particle flux in the ring current. This ROPE offers accurate predictions of omnidirectional flux at a single energy with no pitch angle information, providing robust predictions on the test set with an error score below 11% and calibration scores under 8% with bias under 2% providing a significant speed up as compared to the full RAM‐SCB run.

79 ASTRONOMY AND ASTROPHYSICS↗

An ultra-compact x-ray free-electron laser

In the field of beam physics, two frontier topics have taken center stage due to their potential to enable new approaches to discovery in a wide swath of science. These areas are: advanced, high gradient acceleration techniques, and x-ray free electron lasers (XFELs). Further, there is intense interest in the marriage of these two fields, with the goal of producing a very compact XFEL. In this context, recent advances in high gradient radio-frequency cryogenic copper structure research have opened the door to the use of surface electric fields between 250 and 500 MV m -1 . Such an approach is foreseen to enable a new generation of photoinjectors with six-dimensional beam brightness beyond the current state-of-the-art by well over an order of magnitude. This advance is an essential ingredient enabling an ultra-compact XFEL (UC-XFEL). In addition, one may accelerate these bright beams to GeV scale in less than 10 m. Such an injector, when combined with inverse free electron laser-based bunching techniques can produce multi-kA beams with unprecedented beam quality, quantified by 50 nm-rad normalized emittances. The emittance, we note, is the effective area in transverse phase space (x, p x /m e c) or (y, p y /m e c) occupied by the beam distribution, and it is relevant to achievable beam sizes as well as setting a limit on FEL wavelength. These beams, when injected into innovative, short-period (1–10 mm) undulators uniquely enable UC-XFELs having footprints consistent with university-scale laboratories. We describe the architecture and predicted performance of this novel light source, which promises photon production per pulse of a few percent of existing XFEL sources. We review implementation issues including collective beam effects, compact x-ray optics systems, and other relevant technical challenges. To illustrate the potential of such a light source to fundamentally change the current paradigm of XFELs with their limited access, we examine possible applications in biology, chemistry, materials, atomic physics, industry, and medicine—including the imaging of virus particles—which may profit from this new model of performing XFEL science.

47 OTHER INSTRUMENTATION↗

Bead–Spring Simulation of Ionomer Melts—Studying the Effects of Chain-Length and Associating Group Fraction on Equilibrium Structure and Extensional Flow Behavior

Ionomers are associative polymers with diverse applications ranging from selective membranes and high-performance adhesives to abrasion- and chemical-resistant coatings, insulation layers, vacuum packaging, and foamed sheets. Within equilibrium melt, the ionic or associating groups are known to form thermally reversible, associative clusters whose presence can significantly affect the system’s mechanical, viscoelastic, and transport properties. It is, thus, of great interest to understand how to control such clusters’ size distribution, shape, and stability through the designed choice of polymer architecture and the ionic groups’ fraction, arrangement, and interaction strength. In this work, we represent linear associating polymers using a Kremer–Grest type bead–spring model and perform large-scale MD simulations to explore the effect of polymer chain-length (l) and fraction (f s ) of randomly placed associating groups on the size distribution and stability of formed clusters. We consider different chain-lengths (below and above entanglement), varying fractions of associating groups (represented by ‘sticky’ beads) between 5 and 20%, and a fixed sticky–sticky nonbond interaction strength of four times that between regular non-associating beads. For all melts containing associating groups the equilibrium structure factor S(q) displays a signature ionomer peak at low wave vector q whose intensity increases with increasing f s and l. The average cluster size Nc increases with f s . However, the effect of chain-length on Nc appears to be pronounced only at higher values of f s . Under extensional flows, the computed stress (and viscosity) is higher at higher f s and l regardless of strain rate. Beyond a critical strain rate, we observe fragmentation of the associative clusters, which has interesting effects on the stress/viscous response.

36 MATERIALS SCIENCE↗

Magnetic properties of the Ising-like rare earth pyrosilicate: D-Er 2 Si 2 O 7

Ising-like spin-1/2 magnetic materials are of interest for their ready connection to theory, particularly in the context of quantum critical behavior. In this work we report detailed studies of the magnetic properties of a member of the rare earth pyrosilicate family, D-Er 2 Si 2 O 7 , which is known to display a highly anisotropic Ising-like g-tensor and effective spin-1/2 magnetic moments. We used powder neutron diffraction, powder inelastic neutron spectroscopy (INS), and single crystal AC susceptibility to characterize its magnetic properties. Neutron diffraction enabled us to determine the magnetic structure below the known transition temperature (T N = 1.9 K) in zero field, confirming that the magnetic state is a four-sublattice antiferromagnetic structure with two non-collinear Ising axes, as was previously hypothesized. Our powder INS data revealed a gapped excitation at zero field, consistent with anisotropic (possibly Ising) exchange. An applied field of 1 T produces a mode softening, which is consistent with a field-induced second order phase transition. To assess the relevance of D-Er 2 Si 2 O 7 to the transverse field Ising model, we performed AC susceptibility measurements on a single crystal with the magnetic field oriented in the direction transverse to the Ising axes. This revealed a transition at 2.65 T at 0.1 K, a field significantly higher than the mode-softening field observed by powder INS, showing that the field-induced phase transitions are highly field-direction dependent as expected. Furthermore, these measurements suggest that D-Er 2 Si 2 O 7 may be a candidate for further exploration related to the transverse field Ising model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Community Modeling and Flux Balance Analysis of Synthetic Bacterial Community in MOPS Medium

This narrative was used for the 4-member community modeling using MOPS medium for the bacteria isolated from Populus deltoides (PD). 10 bacterial strains isolated from PD were co-cultured as a synthetic microbial community in MOPS minimal medium. After 15 passages there were 4 members survived in MOPS medium. The genome of each community member was sequenced and uploaded into KBase for genome annotation and metabolic modeling. To simulate the metabolic interactions among the finally surviving community members, the individual metabolic models of each strain were merged into a compartmentalized community model, and then the community model was gapfilled using MOPS medium. The FBA modeling was performed and the metabolites exchange among community members was proposed. The publication by Jia Wang, Dana L. Carper, Leah H. Burdick, Him Shrestha, Manasa Appidi, Paul E. Abraham, Collin M. Timm, Robert L. Hettich, Dale A. Pelletier, Mitchel J. Doktycz can be found here: https://doi.org/10.​1016/​j.​csbj.​2021.​03.​034

59 BASIC BIOLOGICAL SCIENCES↗

Community Modeling and Flux Balance Analysis of Synthetic Bacterial Community in R2A Medium

This narrative was used for the 3-member community modeling using R2A medium for the bacteria isolated from Populus deltoides (PD). 10 bacteiral strains isolated from PD were co-cultured as a synthetic microbial community in R2A complex medium. After 15 passages, there were 3 members survived in R2A medium. The genome of each community member was sequenced and uploaded into KBase for genome annotation and metabolic modeling. To simulate the metabolic interactions among the finally surviving community members, the individual metabolic models of each strain were merged into a compartmentalized community model, and then the community model was gapfilled using R2A medium. The FBA modeling was performed and the metabolites exchange among community members was proposed. The publication by Jia Wang, Dana L. Carper, Leah H. Burdick, Him Shrestha, Manasa Appidi, Paul E. Abraham, Collin M. Timm, Robert L. Hettich, Dale A. Pelletier, Mitchel J. Doktycz can be found here: https://doi.org/10.​1016/​j.​csbj.​2021.​03.​034

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation of precipitation across the contiguous United States, Alaska, and Puerto Rico in multi-decadal convection-permitting simulations

Abstract This study is an early effort to generate a multi-decadal convection-permitting regional climate dataset that covers nearly the entire North American continent. We assessed a 20 year dynamically downscaled regional climate simulation at a 4 km spatial resolution with explicit convection across the contiguous United States (CONUS), Alaska, and Puerto Rico. Specifically, we evaluated the model’s performance in representing mean, 95th percentile, and extreme precipitation across regions. Our findings indicate that when compared with ERA5 reanalysis, the forcing data, convection-permitting simulation improves representations of seasonal, 95th percentile, and extreme precipitation over a large portion of the CONUS, Alaska, and Puerto Rico, particularly in areas where precipitation is heaviest. The simulation adds value over its forcing data (ERA5) in up to 53% of all grid cells in the CONUS, 68.8% in Alaska, and 84.0% in Puerto Rico. It is important to note that, however, despite improvements, model errors in Puerto Rico remain large. Similar improvements are observed in extreme indices, including consecutive dry days, maximum 5 days precipitation, and extreme precipitation. Analysis of the diurnal cycle of mean hourly precipitation suggests that representations of convective processes—including onset, dissipation, suppression, downstream propagation, and local circulation—improved overall.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Two-Dimensional to One-Dimensional Site Response for Idaho National Laboratory

We perform two-dimensional (2D) site response analyses accounting for spatial variability of soil properties and subsurface geometry of the Eastern Snake River Plane (ESRP), and quantify their effects on ground surface motion relative to one-dimensional (1D) site response analyses at the Idaho National Laboratory (INL). We first present the development of random field idealizations of the repeated basalt lava flows, heterogeneously inter-layered with sediments, from seismic velocity data collected over four decades in the ESRP. Using realizations of the stochastic fields mapped on 2D deterministic finite element models, we perform 2D viscoelastic and equivalent-linear wave propagation simulations, and quantify the mean and variance of site response aggravation factors, defined as the response spectral ratio of 2D to 1D analyses on the ground surface. Results are shown to be insensitive to the constitutive material behavior considered here, for strains induced by rock outcrop peak ground acceleration (PGA) as high as 0.7g: viscoelastic and equivalent-linear analyses predict peak mean 2D/1D aggravation factor 1.05 at period T=0.075 sec (i.e. the 2D response spectrum is 5% higher than the corresponding 1D at that period, on average), which corresponds to the wavelength of the horizontal correlation length of the random field (50m). For periods longer than the fundamental period of the site (here, T 1 =0.3125 sec), the propagating wavelengths are too long to be affected by the 1D site response and the aggravation factor becomes equal to 1. The standard deviation of the natural logarithms of the 2D/1D aggravation factors is ~0.15 for periods shorter than the fundamental period of the site, and decays thereafter at a steady rate.

58 GEOSCIENCES↗

A multi-year short-range hindcast experiment with CESM1 for evaluating climate model moist processes from diurnal to interannual timescales

Abstract. We present a multi-year short-range hindcast experiment and its experimental design for better evaluation of both the mean state and variability of atmospheric moist processes in climate models from diurnal to interannual timescales and facilitate model development. We used the Community Earth System Model version 1 as the base model and performed a suite of 3 d hindcasts initialized every day starting at 00:00 Z from 1997 to 2012. Three processes – the diurnal cycle of clouds during different cloud regimes over the central US, precipitation and diabatic heating associated with the Madden–Julian Oscillation (MJO), and the response of precipitation, surface radiative and heat fluxes, as well as zonal wind stress to sea surface temperature anomalies associated with the El Niño–Southern Oscillation – are evaluated as examples to demonstrate how one can better utilize simulations from this experiment to gain insights into model errors and their connection to physical parameterizations or large-scale state. This is achieved by comparing the hindcasts with corresponding long-term observations for periods based on different phenomena. These analyses can only be done through this multi-year hindcast approach to establish robust statistics of the processes under well-controlled large-scale environment because these phenomena are either a result of interannual climate variability or only happen a few times in a given year (e.g., MJO, or cloud regime types). Furthermore, comparison of hindcasts to the typical simulations in climate mode with the same model allows one to infer what portion of a model's climate error directly comes from fast errors in the parameterizations of moist processes. As demonstrated here, model biases in the mean state and variability associated with parameterized moist processes usually develop within a few days and manifest within weeks to affect the simulations of large-scale circulation and ultimately the climate mean state and variability. Therefore, model developers can achieve additional useful understanding of the underlying problems in model physics by conducting a multi-year hindcast experiment.

54 ENVIRONMENTAL SCIENCES↗