Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Consider parameters”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

CLPNets: Coupled Lie–Poisson neural networks for multi-part Hamiltonian systems with symmetries

To accurately compute data-based prediction of Hamiltonian systems, it is essential to utilize methods that preserve the structure of the equations over time. We consider a particularly challenging case of systems with interacting parts that do not reduce to pure momentum evolution. Such systems are essential in scientific computations, such as discretization of a continuum elastic rod, which can be viewed as the group of rotations and translations $SE(3)$. The evolution involves not only the momenta but also the relative positions and orientations of the particles. The presence of Lie group-valued elements, such as relative positions and orientations, poses a problem for applying previously derived methods for data-based computing. We develop a novel method of data-based computation and complete phase space learning of such systems. We follow the original framework of SympNets (Jin et al., 2020) and LPNets (Eldred et al., 2024), building the neural network from phase space mappings that preserve the Lie–Poisson structure. We derive a novel system of mappings that are built into neural networks describing the evolution of such systems. We call such networks Coupled Lie–Poisson Neural Networks, or CLPNets. We consider increasingly complex examples for the applications of CLPNets, starting with the rotation of two rigid bodies about a common axis, progressing to the free rotation of two rigid bodies, and finally to the evolution of two connected and interacting $SE(3)$ components, describing the discretization of an elastic rod into two elements. Our method preserves all Casimir invariants to machine precision, preserves energy to high accuracy, and shows good resistance to the curse of dimensionality, requiring only a few thousand data points for all cases studied (three to eighteen dimensions). Additionally, the method is highly economical in memory requirements, requiring only about 200 parameters for the most complex case considered.

Data-based modeling↗
DOE OSTI2025

Pair momentum dependence of a tilted source in heavy-ion collisions

In noncentral heavy-ion collisions, the particle-emitting source can be tilted away from the beam direction, an effect that becomes particularly significant at collision energies of a few GeV and lower. This phenomenon, manifest itself in many observables such as directed flow, polarization, and vorticity, is therefore important to investigate. Here, in this paper, we study the consistency between the tilt extracted directly from the freeze-out distribution of pions and the tilt parameter obtained using the azimuthally sensitive femtoscopy (asHBT) method. Using the Ultrarelativistic Quantum Molecular Dynamics model, we demonstrate a strong dependence of the tilt parameter extracted with asHBT on the momentum of the particle pair. Considering the experimental challenges in accessing low particle momenta—where the tilt parameter extracted with asHBT closely matches the tilt of the freeze-out distribution of pions—we propose an exponential extrapolation method to obtain the tilt of the entire freeze-out distribution. This approach aims to enhance the accuracy of experimental measurements of tilt in noncentral heavy-ion collisions.

particle correlations & fluctuations↗
DOE OSTI2024

Model-independent search for pair production of new bosons decaying into muons in proton-proton collisions at $\sqrt{s}$ = 13 TeV

The results of a model-independent search for the pair production of new bosons within a mass range of 0.21 < m < 60 GeV, are presented. This study utilizes events with a four-muon final state. We use two data sets, comprising 41.5 fb −1 and 59.7 fb −1 of proton-proton collisions at $\sqrt{s}$ = 13 TeV, recorded in 2017 and 2018 by the CMS experiment at the CERN LHC. The study of the 2018 data set includes a search for displaced signatures of a new boson within the proper decay length range of 0 < cτ < 100 mm. Our results are combined with a previous CMS result, based on 35.9 fb −1 of proton-proton collisions at $\sqrt{s}$ = 13 TeV collected in 2016. No significant deviation from the expected background is observed. Results are presented in terms of a model-independent upper limit on the product of cross section, branching fraction, and acceptance. The findings are interpreted across various benchmark models, such as an axion-like particle model, a vector portal model, the next-to-minimal supersymmetric standard model, and a dark supersymmetric scenario, including those predicting a non-negligible proper decay length of the new boson. In all considered scenarios, substantial portions of the parameter space are excluded, expanding upon prior results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Bayesian optimization of laser wakefield acceleration via spectral pulse shaping

In this paper, we investigate the effect of spectral pulse shaping of the laser driver on the performance of channel-guided, laser–plasma accelerators. The study was carried out with the assistance of Bayesian optimization using particle-in-cell simulations. We used a realistic plasma profile based on a novel optical-field-ionized channel technique with ionization injection and low, on-axis plasma densities to maximize the energy gain of the electron bunch trailing the laser. Spectral shaping allows us to modify the temporal profile of the laser driver while keeping the laser energy constant, affecting the acceleration and injection processes. In addition, we consider how modifying the plasma channel parameters may affect the target outputs. Given the complexity and breadth of the parameter space in question, we used numerical optimization to identify high-performing configurations. In particular, we found laser profiles with additional spectral content that, when used with optimal plasma channel parameters, result in charge content an order of magnitude higher than the baseline Gaussian case while also increasing the mean energy of the electron bunch.

Physics - Plasma physics↗
DOE OSTI2025

Performance of the spin-component-scaled methods for energy bands

The performance of various spin-component-scaled parameterisations is examined for the second-order many-body Green's-function [MBGF(2)] calculations of valence energy bands, taking three of the experimentally well-characterised polymers as examples: polyethylene, polytetrafluoroethylene, and polyacetylene. The parameterisations considered are Grimme's original SCS parameter set, Jung et al.'s original SOS set (retaining the opposite-spin component only), Śmiga et al.'s SCS(IP) set (calibrated specifically for ionization energies), and Śmiga et al.'s SOS(IP) set (calibrated for ionization energies with the opposite-spin component only; implicit in the os-D2 model of Opoku et al.). The SCS(IP) and SOS(IP) parameterisations are found to shift both outer and inner valence bands by up to a few electronvolts away from the experimental data. The original SCS and SOS parameter sets do not improve upon, but largely maintain the accuracy of the unscaled MBGF(2) methods. Given that the SOS-MBGF(2) method can be implemented in a quartic-scaling algorithm (for all roots), it is most promising for solid-state applications. Furthermore this observation is consistent with the success of the quartic-scaling GW methods without the vertex correction based on a density-functional theory reference.

Green’s-function theory↗
DOE OSTI2024

Fault Contribution of Grid-Following and Grid-Forming Inverters Considering Generic Inverter Controls and Ride-through Requirements

Gradual transformation of the synchronous generator (SG)-dominated power systems to inverter-based resource (IBR)-dominated power systems is bringing new challenges regarding power system protection. Increasing IBR’s penetration will reduce the reliability of protection coordination, originally developed for SG-dominated power systems, potentially leading to less accurate calculation or even mis-operation. Thus, new protection schemes compatible with IBR-dominated power systems should be established. Prior to establishing the protection schemes, research on the fault contribution of generic IBR models should be conducted considering interconnection requirements. This paper conducts fault current contribution analysis of generic grid-following (GFL) and grid-forming (GFM) inverters considering interconnection requirements using an electromagnetic transient (EMT) simulation tool. To do this, EMT-based generic GFL models have been developed with the ride-through requirements from IEEE-1547 standards. In addition, the major functions of the generic models are explained. Then, the fault current analysis was conducted using GFL and GFM inverters in a single-inverter-infinite-bus system and a modified IEEE 34-node system considering different grid, fault, inverter control parameters.

Kim, Jinho [BATTELLE (PACIFIC NW LAB)]↗
DOE OSTI2024

Aerosol-Jet Printed Transferable Millimeter-Wave Circuits

Printed millimeter wave (mmW) electronics have been of high interest in the field of communications for some time now due to the potential ability to fabricate mmW systems that utilize 3-Dimension Heterogeneous Integration (3DHI) to improve performance beyond traditional systems. Here, in this paper, a method for fabricating transferable mmW structures via aerosol-jet printing (AJP) is presented. A PDMS stamp assisted liftoff procedure is developed to separate the printed part from a rigid printing support surface and apply it to an adhesive target surface. Several microstrip (MS) line structures are demonstrated to characterize the effect of stamp transfer on the RF performance of the printed circuits. The effects of printed vias on RF probe pads and shifts in the resonant frequencies of a Beatty standard line after transfer are considered and characterized using S-parameter measurements.

42 ENGINEERING↗
DOE OSTI2025

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An Advanced Machine Learning and Artificial Intelligence System for Demonstrating Radiation Regulatory Compliance in DOE Accelerator Facilities

In this Phase II proposal, Applied Research LLC (ARLLC), Thomas Jefferson National Accelerator Facility (Jefferson Lab), and Old Dominion University (ODU) propose the combination of domain knowledge (beam characteristics, fixed structural shielding, earthen burden (the soil and foliage added to the dome of the experimental halls as additional shielding), etc.), machine learning (ML) and/or artificial intelligence (AI) to correlate a variety of multi-modal onsite signals and the radiation fields seen in accessible areas of the accelerator site and the site boundary. The ML/AI will consider the complex influence of environmental parameters affecting the radon contribution of the measurements, focusing on actual data obtained from Jefferson Lab. In Phase I, the coded beam and location data were fed into a deep learning model to predict doses at several designated locations in Jefferson Lab’s facility. Moreover, a dense radiation map was generated using only a sparse collection of the samples in a facility. In Phase II, we will develop a software prototype containing a radiation prediction algorithm, dense radiation map algorithms, and background noise prediction algorithms, with actual data used to evaluate the prototype. This work will provide a framework for evaluation of radiation measurement results around the site based on learned responses. In addition, the proposed approach allows more granular mapping of radiation levels. Better understanding and communication of these levels is related to the overall approach in keeping doses to personnel ALARA.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Evaluating the impact of anatomical and physiological variability on human equivalent doses using PBPK models

Abstract Addressing human anatomical and physiological variability is a crucial component of human health risk assessment of chemicals. Experts have recommended probabilistic chemical risk assessment paradigms in which distributional adjustment factors are used to account for various sources of uncertainty and variability, including variability in the pharmacokinetic behavior of a given substance in different humans. In practice, convenient assumptions about the distribution forms of adjustment factors and human equivalent doses (HEDs) are often used. Parameters such as tissue volumes and blood flows are likewise often assumed to be lognormally or normally distributed without evaluating empirical data for consistency with these forms. In this work, we performed dosimetric extrapolations using physiologically based pharmacokinetic (PBPK) models for dichloromethane (DCM) and chloroform that incorporate uncertainty and variability to determine if the HEDs associated with such extrapolations are approximately lognormal and how they depend on the underlying distribution shapes chosen to represent model parameters. We accounted for uncertainty and variability in PBPK model parameters by randomly drawing their values from a variety of distribution types. We then performed reverse dosimetry to calculate HEDs based on animal points of departure for each set of sampled parameters. Corresponding samples of HEDs were tested to determine the impact of input parameter distributions on their central tendencies, extreme percentiles, and degree of conformance to lognormality. This work demonstrates that the measurable attributes of human variability should be considered more carefully and that generalized assumptions about parameter distribution shapes may lead to inaccurate estimates of extreme percentiles of HEDs.

Toxicology↗
DOE OSTI2024

Uncertainty quantification of material parameters in modeling coupled metal and high explosive experiments

Experiments involving the coupling of metal and high explosives (HE) are of notable defense-related interest, and we seek to refine the uncertainty quantification associated with models of such experiments. In particular, our focus is on how uncertainty related to the metal constitutive model challenges our ability to infer high explosive model parameters when analyzing focused science experiments. We consider three focused experiments involving an HE accelerating metal: small plate tests with tantalum/LX-14 and tantalum/LX-17 pairings as well as a tantalum/LX-17 cylinder test. For all three models, we perform sensitivity analysis to ascertain the influence of metal strength on the coupled experimental response. Moreover, we calibrate each model in a Bayesian setting and study the quantification of metal strength on the inference of the HE parameters. Based on our results, we offer guidance for future metal/HE experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An international study on THM modelling of the full-scale heater experiment at Mont Terri laboratory

We present results from an international model comparison study of the Full-Scale Emplacement (FE) experiment in Opalinus Clay at the Mont Terri Laboratory, Switzerland. Based on a provided parameter set the teams decided which parameters they adopted for their models, whether they considered the excavation and the ventilation phase in addition to the heating phase and if they included technical features like the shotcrete or the EDZ. The teams were able to reproduce the measured parameters temperature, relative humidity and pore pressure. The modelled results for temperature agree very closely between the teams especially in the sensors in Opalinus Clay. All teams were able to reproduce the redistribution of water in the bentonite backfill due to heating. The evolution of the relative humidity showed similar trends with differences in the intensity of the dry out effect. To model the pore pressure evolution is more complex because it comprises the full interaction of the coupled THM processes. The spread between the pore pressure modelled by the teams was larger, with some teams overestimating the pressure increase due to heating and some teams overestimating the extent of drainage. The agreement of modelled results with measurements improves with larger distance to the heater. We conclude that the EDZ and the shotcrete potentially influence the behaviour of the rock causing higher differences closer to the heater. Further research is needed to better implement those influences into the models. Based on the calibrated models, the future evolution of temperature, relative humidity and pore pressure was predicted over the next 10 years following a change of the heat power applied in 2023 and 2024. Again, the predicted temperatures agree very closely between the teams. Most teams do not expect an increase in relative humidity during the next 10 years after the initial dry-out.

58 GEOSCIENCES↗
DOE OSTI2025

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U↗
DOE OSTI2025

Identifying high-impact and high-uncertainty parameters in MiniFuel model predictions

The MiniFuel irradiation platform at Oak Ridge National Laboratory's High Flux Isotope Reactor (HFIR) is a flexible, high-throughput separate effects test capability. Finite element thermal models are relied upon to design MiniFuel experiments and to achieve experimental objectives. Recent reports show good agreement in the model prediction of target fuel temperatures, but as the capability of the experiments is extended to higher temperatures, the uncertainty in the model predictions must be quantified. To that end, high-impact, high-uncertainty parameters that contribute the most uncertainty to the model are identified. The uncertainty quantification was accomplished through a series of screening and sensitivity analyses. The first analysis utilizes the method of Morris to perform a computationally efficient preliminary screening that considers uncertainty in a large number of the model inputs. The most important parameters identified in the Morris screening study were then considered in a Sobol sensitivity analysis that more robustly ranks and quantifies the uncertainty associated with each parameter. From these analyses, it was determined that thermal contact conductance between components is the parameter that contributes the highest uncertainty. The estimated uncertainty of the MiniFuel model fuel temperature predictions is ±80 °C in the removable beryllium and ±40 °C in the vertical experiment facilities. In conclusion, the framework established by the series of sensitivity analyses presented herein could easily be adapted to fit the needs of accelerated fuel qualification processes.

Fuel, Irradiation↗
DOE OSTI2024

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Photosynthetic responses to light levels in drought-tolerant novel peanut (Arachis hypogaea L) genotypes

Abstract Drought is a significant abiotic stressor that reduces peanut production because it alters photosynthetic activity and impacts crop growth. Therefore, developing drought-tolerant peanut genotypes capable of maintaining higher photosynthetic rates (A) under stress is crucial. This study assessed changes in photosynthetic and chlorophyll fluorescence responses to light (photosynthetic photon flux density, PPFD) in newly bred drought-tolerant peanut genotypes. Ten genotypes [NM-3, NM-5, NM-6, NM-23, NM-69, NM-70, NM-74, NM-77, V-C, and C-76–16] were evaluated under full irrigation (FC 100 ) and deficit irrigation (FC 50 ) in a split-plot design with four replications in a greenhouse. Under high PPFD levels, genotype NM-5 with deficit irrigation exhibited significantly higherA, stomatal conductance (gs), quantum efficiency of photosystem II (ΦPSII), and electron transport rate (ETR) by 40–59%, 135–525%, 31–212%, and 31–102%, respectively, than check varieties (V-C and C-76–16) and other genotypes. The NM-74 and NM-77 genotypes also performed well under deficit irrigations but with slightly lowerA,gs, ΦPSII, and ETR. Genotypes NM-5, NM-23, NM-74, and NM-77 exhibited significantly higher quantum efficiency of photosystem II (Fv’/Fm’) and photochemical quenching (qP) with higher light intensities in the daily cycle under deficit irrigation. The decline in ETR at the same PPFD levels in NM-3, NM 69, NM-70, and C-76–16 indicated photoinhibition or saturation of the photosynthetic apparatus compared to other genotypes. Concurrently, FC 100 irrigation level minimizes photoinhibition, enhancingA,gs, ΦPSII, and ETR in the genotypes than FC 50 . Therefore, we conclude that NM-5, NM-74, and NM-77 genotypes can perform better under water deficit environments. As such, chlorophyll fluorescence parameters Fv’/Fm’ and qP can be considered for selective breeding to enhance photosynthetic efficiencies.

Science & Technology - Other Topics↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2024

Uncertainties in greenhouse gas emission factors: A comprehensive analysis of switchgrass‐based biofuel production

Abstract This study investigates uncertainties in greenhouse gas (GHG) emission factors related to switchgrass‐based biofuel production in Michigan. Using three life cycle assessment (LCA) databases—US lifecycle inventory (USLCI) database, GREET, and Ecoinvent—each with multiple versions, we recalculated the global warming intensity (GWI) and GHG mitigation potential in a static calculation. Employing Monte Carlo simulations along with local and global sensitivity analyses, we assess uncertainties and pinpoint key parameters influencing GWI. The convergence of results across our previous study, static calculations, and Monte Carlo simulations enhances the credibility of estimated GWI values. Static calculations, validated by Monte Carlo simulations, offer reasonable central tendencies, providing a robust foundation for policy considerations. However, the wider range observed in Monte Carlo simulations underscores the importance of potential variations and uncertainties in real‐world applications. Sensitivity analyses identify biofuel yield, GHG emissions of electricity, and soil organic carbon (SOC) change as pivotal parameters influencing GWI. Decreasing uncertainties in GWI may be achieved by making greater efforts to acquire more precise data on these parameters. Our study emphasizes the significance of considering diverse GHG factors and databases in GWI assessments and stresses the need for accurate electricity fuel mixes, crucial information for refining GWI assessments and informing strategies for sustainable biofuel production.

Kim, Seungdo↗