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 703 records · Page 39

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

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↗

Integrated Geomechanics and Geophysics in Induced Seismicity: Mechanisms and Monitoring

Quantifying in-situ subsurface stresses and predicting fracture development are critical to reducing risks of induced seismicity and improving modern energy activities in the subsurface. In this work, we developed a novel integration of controlled mechanical failure experiments coupled with microCT imaging, acoustic sensing, modeling of fracture initiation and propagation, and machine learning for event detections and waveform characterization. Through additive manufacturing (3D printing), we were able to produce bassanite-gypsum rock samples with repeatable physical, geochemical and structural properties. With these "geoarchitected" rock, we provided the role of mineral texture orientation on fracture surface roughness. The impact of poroelastic coupling on induced seismicity has been systematically investigated to improve mechanistic understanding of post shut-in surge of induced seismicity. This research will set the groundwork for characterizing seismic waveforms by using multiphysics and machine learning approaches and improve the detection of low-magnitude seismic events leading to the discovery of hidden fault/fracture systems.

42 ENGINEERING↗

Screening Level Ecological Exposure Assessment Work Plan

This work plan is part of an effort to complete an updated evaluation of the potential for adverse effects from miscellaneous treatment unit emissions to potential ecological receptors within and surrounding Site 200. It will address potential pathways between air emissions from the Lawrence Livermore National Laboratory’s (LLNL) Decontamination and Waste Treatment Facility (DWTF), single emission stack, and ecological receptors (i.e., plants and wildlife) present onsite and within a one-mile radius from the DWTF stack. This work plan outlines the procedures by which data collected as part of the previous investigations will be evaluated to assess risks to ecological receptors. It presents the problem formulation, initial conceptual site models, methods for characterizing exposure, toxicity, and risks that will be used if a more detailed “baseline” ecological risk assessment (BERA) is required.

54 ENVIRONMENTAL SCIENCES↗

Enabling Real-time Scattering Data Analysis with Scalable Optimization [Slides]

Diffraction experiments produce datasets with rich multidimensional physics information such as microstructure, equations of state, crystal structure, elastoplastic properties, and other key inputs to LANL mission-essential multiphysics models. This information is typically extracted through a process called Rietveld refinement, which involves selecting appropriate models of the instrument, crystal structure, and microstructure, identifying suitable starting values, and then fitting often hundreds of model parameters using a sequence of empirical parameter turnon/off sequences within a non-global gradient-based optimization. Extensive user expertise is required to properly setup a refinement, identify appropriate models, and select initial parameter values close to truth, such that the refinement will yield parameter values that are optimally predictive. This is a very tedious manual process performed far after the beamline campaign has ended. As facilities have become capable of generating larger volumes of data, the limitation in throughput due to Rietveld refinement has led to a dramatic increase in unanalyzed data as opposed to an intended increase in new science. In our FY22 TED, we demonstrated an integrated toolset providing near real-time automated Rietveld analysis. If this toolset can be optimized to provide automated Rietveld analysis in real-time, this could alleviate the bottleneck in unanalyzed diffraction data, aid in decision-making during experiments, and increase efficiency of the facility.

74 ATOMIC AND MOLECULAR PHYSICS↗

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.↗