Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model initialization”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Inverse design of photonic surfaces via multi fidelity ensemble framework and femtosecond laser processing

We demonstrate a multi-fidelity (MF) machine learning ensemble framework for the inverse design of photonic surfaces, trained on a dataset of 11,759 samples that we fabricate using high throughput femtosecond laser processing. The MF ensemble combines an initial low fidelity model for generating design solutions, with a high fidelity model that refines these solutions through local optimization. The combined MF ensemble can generate multiple disparate sets of laser-processing parameters that can each produce the same target input spectral emissivity with high accuracy (root mean squared errors < 2%). SHapley Additive exPlanations analysis shows transparent model interpretability of the complex relationship between laser parameters and spectral emissivity. Finally, the MF ensemble is experimentally validated by fabricating and evaluating photonic surface designs that it generates for improved efficiency energy harvesting devices. Our approach provides a powerful tool for advancing the inverse design of photonic surfaces in energy harvesting applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Probing remdesivir nucleotide analogue insertion to SARS-CoV-2 RNA dependent RNA polymerase in viral replication

Remdesivir (RDV) prodrug can be metabolized into a triphosphate form nucleotide analogue (RDV-TP) to bind and insert into the active site of viral RNA dependent RNA polymerase (RdRp) to further interfere with viral genome replication. Here in this work, we computationally studied how RDV-TP binds and inserts to the SARS-CoV-2 RdRp active site, in comparison with natural nucleotide substrate adenosine triphosphate (ATP). To do that, we first constructed atomic structural models of an initial binding complex (active site open) and a substrate insertion complex (active site closed), based on high-resolution cryo-EM structures determined recently for SARS-CoV-2 RdRp or non-structural protein (nsp) 12, in complex with accessory protein factors nsp7 and nsp8. By conducting all-atom molecular dynamics simulation with umbrella sampling strategies on the nucleotide insertion between the open and closed state RdRp complexes, our studies show that RDV-TP can initially bind in a comparatively stabilized state to the viral RdRp active site, as it primarily forms base stacking with the template uracil nucleotide (nt +1), which under freely fluctuations supports a low free energy barrier of the RDV-TP insertion (~1.5 kcal mol -1 ). In comparison, the corresponding natural substrate ATP binds initially to the RdRp active site in Watson–Crick base pairing with the template nt, and inserts into the active site with a medium low free energy barrier (~2.6 kcal mol -1 ), when the fluctuations of the template nt are well quenched. The simulations also show that the initial base stacking of RDV-TP with the template can be specifically stabilized by motif C-S759, S682 (near motif B) with the base, and motif G-K500 with the template backbone. Although the RDV-TP insertion can be hindered by motif F-R555/R553 interaction with the triphosphate, the ATP insertion seems to be facilitated by such interactions. The inserted RDV-TP and ATP can be further distinguished by specific sugar interaction with motif B-T687 and motif A-D623, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Modeling of toroidal torques exerted by internal kink instability in a tokamak plasma

Toroidal modeling efforts are initiated to systematically compute and compare various toroidal torques, exerted by an unstable internal kink in a tokamak plasma, using the MARS-F/K/Q suite of codes. The torques considered here include the resonant electromagnetic torque due to the Maxwell stress (the EM or JXB torque), the neoclassical toroidal viscous (NTV) torque, and the torque associated with the Reynolds stress. Numerical results show that the relative magnitude of the net resonant electromagnetic and the Reynolds stress torques increases with the equilibrium flow speed of the plasma, whilst the net NTV torque follows the opposite trend. The global flow shear sensitively affects the Reynolds stress torque, but not the electromagnetic and the NTV torques. Detailed examinations reveal dominant contributions to the Maxwell and Reynolds stress torques, in terms of the poloidal harmonic numbers of various perturbation fields, as well as their relative toroidal phasing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Robust unfolding of MeV x-ray spectra from filter stack spectrometer data

Here, we present an inversion method capable of robustly unfolding MeV x-ray spectra from filter stack spectrometer (FSS) data without requiring an a priori specification of a spectral shape or arbitrary termination of the algorithm. Our inversion method is based upon the perturbative minimization (PM) algorithm, which has previously been shown to be capable of unfolding x-ray transmission data, albeit for a limited regime in which the x-ray mass attenuation coefficient of the filter material increases monotonically with x-ray energy. Our inversion method improves upon the PM algorithm through regular smoothing of the candidate spectrum and by adding stochasticity to the search. With these additions, the inversion method does not require a physics model for an initial guess, fitting, or user-selected termination of the search. Instead, the only assumption made by the inversion method is that the x-ray spectrum should be near a smooth curve. Testing with synthetic data shows that the inversion method can successfully recover the primary large-scale features of MeV x-ray spectra, including the number of x-rays in energy bins of several-MeV widths to within 10%. Fine-scale features, however, are more difficult to recover accurately. Examples of unfolding experimental FSS data obtained at the Texas Petawatt Laser Facility and the OMEGA EP laser facility are also presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Forecasting high-dimensional spatio-temporal systems from sparse measurements

This paper introduces a new neural network architecture designed to forecast high-dimensional spatio-temporal data using only sparse measurements. The architecture uses a two-stage end-to-end framework that combines neural ordinary differential equations (NODEs) with vision transformers. Initially, our approach models the underlying dynamics of complex systems within a low-dimensional space; and then it reconstructs the corresponding high-dimensional spatial fields. Many traditional methods involve decoding high-dimensional spatial fields before modeling the dynamics, while some other methods use an encoder to transition from high-dimensional observations to a latent space for dynamic modeling. In contrast, our approach directly uses sparse measurements to model the dynamics, bypassing the need for an encoder. This direct approach simplifies the modeling process, reduces computational complexity, and enhances the efficiency and scalability of the method for large datasets. We demonstrate the effectiveness of our framework through applications to various spatio-temporal systems, including fluid flows and global weather patterns. Although sparse measurements have limitations, our experiments reveal that they are sufficient to forecast system dynamics accurately over long time horizons. Our results also indicate that the performance of our proposed method remains robust across different sensor placement strategies, with further improvements as the number of sensors increases. This robustness underscores the flexibility of our architecture, particularly in real-world scenarios where sensor data is often sparse and unevenly distributed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The effect of saturated thermal conduction on clouds in a hot plasma

We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Barless flocculent galaxies: a dynamic puzzle

We draw attention to the bright galaxies that do not show a bar in their structure but have a flocculent spiral structure. Using the THINGS and HERACLES kinematic data for four barless galaxies (NGC 2841, NGC 3512, NGC 5055, NGC 7331), we built their mass models including dark haloes. We concluded that the fraction of the dark matter does not exceed 50 per cent within the optical radii of the galaxies. This is too little to explain the lack of a bar in these galaxies. In an attempt to understand the featureless structure of these galaxies, we constructed several N-body models with an initially reduced content of dark matter. We concluded that, in addition to the low mass of the dark halo, the decisive factor that leads to a barless disc is the start from an initially unstable state. An isolated dynamically cold disc (with the Toomre parameter Q < 0.5) settled into rotational equilibrium passes through the short stage of violent instability with fragmentation and formation of stellar clumps. After that, it evolves passively and ends up with a featureless structure. We assume that the barless flocculent galaxies studied in this work may be descendants of galaxies at high redshifts with rotation curves that are consistent with the high mass fraction of baryons relative to the total dark matter halo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Inference of neutrino flavor evolution through data assimilation and neural differential equations

The evolution of neutrino flavor in dense environments such as core-collapse supernovae and binary compact object mergers constitutes an important and unsolved problem. Its solution has potential implications for the dynamics and heavy-element nucleosynthesis in these environments. In this paper, we build upon recent work to explore inference-based techniques for the estimation of model parameters and neutrino flavor evolution histories. Here, we combine data assimilation, ordinary differential equation solvers, and neural networks to craft an inference approach tailored for nonlinear dynamical systems. Using this architecture, and a simple two-neutrino-beam, two-flavor model, we compare the performances of nine different optimization algorithms and expand upon previous assessments of the efficacy of inference for tackling problems in flavor evolution. We find that employing this new architecture, together with evolutionary optimization algorithms, accurately captures flavor histories in the small-scale model and allows us to quickly explore both model parameters and initial flavor content. In future work we plan to extend these inference techniques to large numbers of neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Direct Measurement of Nuclear Dependence of Charged Current Quasielasticlike Neutrino Interactions Using MINERvA

Charged-current ν μ interactions on carbon, iron, and lead with a final state hadronic system of one or more protons with zero mesons are used to investigate the influence of the nuclear environment on quasielasticlike interactions. The transferred four-momentum squared to the target nucleus, Q 2 , is reconstructed based on the kinematics of the leading proton, and differential cross sections versus Q 2 and the cross-section ratios of iron, lead, and carbon to scintillator are measured for the first time in a single experiment. The measurements show a dependence on the atomic number. While the quasielasticlike scattering on carbon is compatible with predictions, the trends exhibited by scattering on iron and lead favor a prediction with intranuclear rescattering of hadrons accounted for by a conventional particle cascade treatment. Furthermore, these measurements help discriminate between different models of both initial state nucleons and final state interactions used in the neutrino oscillation experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Iterative Linearization for Phasor-Defined Optimal Power Dispatch

Optimal power flow (OPF) problems, which dispatch power targets to controllable generating units across a network, must generally account for non-convex constraints on power flow. Furthermore, adapting those problems so as to make them solvable with convex optimization techniques is an area of much academic and operational interest. In this paper, we present a method for solving OPF as a quadratic program by iteratively refining and re-initializing a linearized model of power flow based on the outputs of an associated nonlinear solver. The linear model on which we demonstrate this method is an adapted version of an approximation designed for use with unbalanced distribution networks. As an important benefit, the model allows for the explicit inclusion of nodal voltage phasor values in both the OPF problem's objective and its constraints, which opens the door to the idea of phasor-based control (PBC) design. We show in simulations on the IEEE 13-node test feeder that our method quickly converges to a set of phasor targets that are sufficiently precise for use in operations at the distribution level.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Biolink Model: A universal schema for knowledge graphs in clinical, biomedical, and translational science

Abstract Within clinical, biomedical, and translational science, an increasing number of projects are adopting graphs for knowledge representation. Graph‐based data models elucidate the interconnectedness among core biomedical concepts, enable data structures to be easily updated, and support intuitive queries, visualizations, and inference algorithms. However, knowledge discovery across these “knowledge graphs” (KGs) has remained difficult. Data set heterogeneity and complexity; the proliferation of ad hoc data formats; poor compliance with guidelines on findability, accessibility, interoperability, and reusability; and, in particular, the lack of a universally accepted, open‐access model for standardization across biomedical KGs has left the task of reconciling data sources to downstream consumers. Biolink Model is an open‐source data model that can be used to formalize the relationships between data structures in translational science. It incorporates object‐oriented classification and graph‐oriented features. The core of the model is a set of hierarchical, interconnected classes (or categories) and relationships between them (or predicates) representing biomedical entities such as gene, disease, chemical, anatomic structure, and phenotype. The model provides class and edge attributes and associations that guide how entities should relate to one another. Here, we highlight the need for a standardized data model for KGs, describe Biolink Model, and compare it with other models. We demonstrate the utility of Biolink Model in various initiatives, including the Biomedical Data Translator Consortium and the Monarch Initiative, and show how it has supported easier integration and interoperability of biomedical KGs, bringing together knowledge from multiple sources and helping to realize the goals of translational science.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Existence and Optimality Conditions for Risk-Averse PDE-Constrained Optimization

Uncertainty is ubiquitous in virtually all engineering applications, and, for such problems, it is inadequate to simulate the underlying physics without quantifying the uncertainty in unknown or random inputs, boundary and initial conditions, and modeling assumptions. Here in this paper, we introduce a general framework for analyzing risk-averse optimization problems constrained by partial differential equations (PDEs). In particular, we postulate conditions on the random variable objective function as well as the PDE solution that guarantee existence of minimizers. Furthermore, we derive optimality conditions and apply our results to the control of an environmental contaminant. Lastly, we introduce a new risk measure, called the conditional entropic risk, that fuses desirable properties from both the conditional value-at-risk and the entropic risk measures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

The effects of space charge on the dynamics of the booster synchrotron in the Jefferson Lab EIC

Optimization of the booster synchrotron design to operate in the extreme space-charge dominated regime is proposed. This study is motivated by the ultra-high luminosity promised by the JLEIC accelerator complex, which poses several beam dynamics and lattice design challenges for its individual components. We examine the effects of space charge on the dynamics of the booster synchrotron for the proposed JLEIC electron ion collider. This booster will inject and accumulate protons and heavy ions at an energy of 285 MeV and then engage in a process of acceleration and electron cooling to bring it to its extraction energy of 7 GeV. This would then be sent into the ion collider ring part of JLEIC. In order to examine the effects of space charge on the dynamics of this process we use the software SYNERGIA. We model both the initial coasting beam as well as proposed acceleration bunch properties to evaluate the best methods of creating the JLEIC beam.

Nissen, Edith↗
DOE OSTI2024

SRF Cavity Instability Detection with Machine Learning at CEBAF

During the operation of the Continuous Electron Beam Accelerator Facility (CEBAF), one or more unstable superconducting radio-frequency (SRF) cavities often cause beam loss trips while the unstable cavities themselves do not necessarily trip off. The present RF controls for the legacy cavities report at only 1 Hz, which is too slow to detect fast transient instabilities during these trip events. These challenges make the identification of an unstable cavity out of the hundreds installed at CEBAF a difficult and time-consuming task. To tackle these issues, a fast data acquisition system (DAQ) for the legacy SRF cavities has been developed, which records the sample at 5 kHz. A Principal Component Analysis (PCA) approach is being developed to identify anomalous SRF cavity behavior. We will discuss the present status of the DAQ system and PCA model, along with initial performance metrics. Overall, our method offers a practical solution for identifying unstable SRF cavities, contributing to increased beam availability and machine reliability.

Carpenter, A.↗
DOE OSTI2017

Complete Report on the Development of Welding Parameters for Irradiated Materials

The advanced welding facility at the Radiochemical Engineering Development Center of Oak Ridge National Laboratory, which was conceived to enable research and development of weld repair techniques for nuclear power plant life extension, is now operational. The development of the facility and its advanced welding capabilities, along with the model materials for initial welding trials, were funded jointly by the U.S. Department of Energy, Office of Nuclear Energy, Light Water Reactor Sustainability Program, the Electric Power Research Institute, Long Term Operations Program and the Welding and Repair Technology Center, with additional support from Oak Ridge National Laboratory. Welding of irradiated materials was initiated on November 17, 2017, which marked a significant step in the development of the facility and the beginning of extensive welding research and development campaigns on irradiated materials that will eventually produce validated techniques and guidelines for weld repair activities carried out to extend the operational lifetimes of nuclear power plants beyond 60 years. This report summarizes the final steps that were required to complete weld process development, initial irradiated materials welding activities, near-term plans for irradiated materials welding, and plans for post-weld analyses that will be carried out to assess the ability of the advanced welding processes to make repairs on irradiated materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Additive Manufacturing Interface Profiles for Jet Studies on the Vertical Shock Tube

Understanding the mixing between two fluids of disparate densities (i.e., variable density mixing) is important for the lab's mission. Experiments are required to validate our codes, such as RAGE, FLAG, and the BHR mix model. The role initial conditions (ICs) have on instability growth and the transition to turbulence is not completely understood and comparisons with experiments are needed. The Extreme Fluids team focuses on experiments at Los Alamos that utilize gases and low shock velocities (Mach <2 in air). These experiments have the advantage of high fidelity diagnostics capable of temporally and spatially resolved velocity and density measurements. One main difficulty is setting up a specific interface profile between the two gases. While attempts have been made to do so without the aid of a solid structure (i.e., a membrane), the results are limited in the range of possible shapes and reproducibility. Membranes provide the best way to provide a deterministic initial condition. However, fragments from the membrane can interfere with the flow and instability growth, as well as the diagnostics. In this project, we focused on exploring the materials and techniques available via modern additive manufacturing (i.e., 3D printing) to test if they could address some of the shortcomings of prior membrane designs. Our goal was to find a material that would have minimal impact on our experiment while also providing us with a means to produce complex initial conditions. This report details our work to explore the range of additive manufacturing capabilities and the jet-like structures we were able to generate using three-dimensional initial conditions.

36 MATERIALS SCIENCE↗