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 721 records · Page 40

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

Helium recovery system at IB3a

The increasing need for optimal and sustainable use of cryogenic resources to support Fermilab’s scientific mission has highlighted the necessity of improving the Laboratory’s helium management practices. An assessment of cryogenic test facilities identified the Technical Division’s Industrial Building 3a (IB3A) as a key site requiring upgrades to integrate a helium recovery system. The IB3A facility is essential for characterizing and testing superconductors, cables, and coils for various R&D projects, including the US High-Luminosity LHC Accelerator Upgrade Project (AUP), Mu2e, and other external collaborations. Currently, the facility relies on 500 L helium Dewars and vents the vaporized helium directly into the atmosphere, leading to significant helium loss. Given the non-renewable nature of helium, recovering and reusing this resource is critical for the sustainability of Fermilab’s operations. To address this challenge, a project has been initiated to connect IB3A to an existing helium purification station and refrigeration system located in another building via a dedicated pipeline pass over the roof of several buildings. This solution will enable the efficient capture of vented helium, its transfer to the purification station, and subsequent liquefaction for reuse in future operations. The project includes a detailed design phase, specifying the pipeline route, flow control mechanisms, and integration with the existing cryogenic infrastructure, followed by phased implementation and commissioning. By implementing this pipeline connection and upgrading IB3A, Fermilab aims to significantly reduce helium waste, lower operational costs, and align with its commitment to sustainability. This initiative provides a model for resource-efficient cryogenic operations and reinforces the Laboratory’s capacity to support its science mission for the long term.

Porwisiak, D. [Fermilab]↗

Data Structure Alchemy

In an increasingly more data-driven world, the project set out to uncover the first principles of data-structure design, chart the immense design space they form, and build automation that can synthesize an optimal structure, or even a whole storage engine, for any given workload, hardware platform, and cost target. Data structures are at the center of every computational system and are directly responsible for its performance. Two core technical thrusts were defined: 1) Mapping design spaces for key data-centric abstractions (filters, hash functions, storage-engine layouts, neural-network topologies, blockchain protocols, image layouts, etc.). 2) Developing search & synthesis algorithms, initially analytical cost models, later neural-guided bi-level optimisers that navigate sextillions of candidate designs in seconds and materialise the best one as ready‐to-run code. This report distills the key insights, accomplishments, and impact.

97 MATHEMATICS AND COMPUTING↗

Adaptable Standards for Discovery, Access, and Usability of Oak Ridge National Laboratory’s Data Portals and Catalogs

Oak Ridge National Laboratory (ORNL) is leveraging its established capabilities and subject matter expertise in data curation, governance, management, national security, and risk assessment and mitigation to support the US Department of Energy (DOE) Grid Modernization Initiative. Using standards modeled by the National Institute of Standards and Technology (NIST), the Data Curation Network (DCN), the Oak Ridge Leadership Computing Facility (OLCF), and other leading organizations in the fields of energy research, high-performance computing, and national and homeland security, ORNL seeks to provide a federated approach to research data discovery, use, and interoperability.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

MAT 332: Recyclable-By-Design Carbon Fiber Composites for Battery Enclosures

This work aims to fabricate a small-scale battery module with the baseline RBD-CFRC and conduct mechanical crush test to develop an initial thermal runaway model. Test the adhesive properties of the RBD-CFRCs using their resin constituent, ultrasonic welding, and other resins employed today. Recommend a future direction for resin development to enable the welding of RBDCFRCs to dissimilar materials. Time and Cost Permitting. Iterate on resin development for flame resistance and adhesive properties, demonstrating how thermal runaway can be reduced or joining properties can be improved.

33 ADVANCED PROPULSION SYSTEMS↗

What does the matter created in high multiplicity proton-nucleus collisions teach us about the 3-D structure of the proton?

The study of multiparticle correlations and collectivity in the hot and dense matter created in collisions of heavy-ions as well as those of smaller systems such as proton--heavy-ion collisions has progressed to the point where detailed knowledge of the three-dimensional structure of the proton is needed to confront experimental data. We discuss results from a simple proof-of-principle initial state parton model which reproduces many features of the data on proton--heavy-ion collisions that are often ascribed to hydrodynamic flow. We outline how this model can be further improved with particular emphasis on missing elements in our understanding of the dynamical spatial and momentum structure of the proton.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Classification Analysis of Southwest Pacific Tropical Cyclone Intensity Changes Prior to Landfall

This study evaluates the ability of a random forest classifier to identify tropical cyclone (TC) intensification or weakening prior to landfall over the western region of the Southwest Pacific Ocean (SWPO) basin. For both Australia mainland and SWPO island cases, when a TC first crosses land after spending ≥24 h over the ocean, the closest hour prior to the intersection is considered as the landfall hour. If the maximum wind speed (V max ) at the landfall hour increased or remained the same from the 24-h mark prior to landfall, the TC is labeled as intensifying and if the V max at the landfall hour decreases, the TC is labeled as weakening. Geophysical and aerosol variables closest to the 24 h before landfall hour were collected for each sample. The random forest model with leave-one-out cross validation and the random oversampling example technique was identified as the best-performing classifier for both mainland and island cases. The model identified longitude, initial intensity, and sea skin temperature as the most important variables for the mainland and island landfall classification decisions. Incorrectly classified cases from the test data were analyzed by sorting the cases by their initial intensity hour, landfall hour, monthly distribution, and 24-h intensity changes. TC intensity changes near land strongly impact coastal preparations such as wind damage and flood damage mitigations; hence, this study will contribute to improve identifying and prioritizing prediction of important variables contributing to TC intensity change before landfall.

54 ENVIRONMENTAL SCIENCES↗

A Data-Driven Workflow Approach to Optimization of Fracture Spacing in Multi-Fractured Shale Oil Wells

A data-driven workflow approach is presented in this study for optimizing fracture spacing of multifractured horizontal wells (MFHW) in shale oil reservoirs. The workflow employs a simple well productivity model for the initial design of hydraulic fracturing well completions. This provides a transparent approach to the identification of key fracturing parameters affecting well productivity. The workflow uses transient pressure or production data to identify fracture interference. This offers a reliable and cost-effective means for assessment of well production potential in terms of optimization of fracture spacing in the MFHW. Result of a field case study indicated that three wells were drilled in an area with dense natural fractures, and the fracture spacing of MFHW in this area was short enough to effectively drain the stimulated reservoir volume (SRV), while the other three wells were drilled in an area with less natural fractures, and the fracture spacing of MFHW in this area could be shortened to double well productivity.

58 GEOSCIENCES↗

Calibration and Localization of Optically Pumped Magnetometers Using Electromagnetic Coils

In this paper, we propose a method to estimate the position, orientation, and gain of a magnetic field sensor using a set of (large) electromagnetic coils. We apply the method for calibrating an array of optically pumped magnetometers (OPMs) for magnetoencephalography (MEG). We first measure the magnetic fields of the coils at multiple known positions using a well-calibrated triaxial magnetometer, and model these discreetly sampled fields using vector spherical harmonics (VSH) functions. We then localize and calibrate an OPM by minimizing the sum of squared errors between the model signals and the OPM responses to the coil fields. We show that by using homogeneous and first-order gradient fields, the OPM sensor parameters (gain, position, and orientation) can be obtained from a set of linear equations with pseudo-inverses of two matrices. The currents that should be applied to the coils for approximating these low-order field components can be determined based on the VSH models. Computationally simple initial estimates of the OPM sensor parameters follow. As a first test of the method, we placed a fluxgate magnetometer at multiple positions and estimated the RMS position, orientation, and gain errors of the method to be 1.0 mm, 0.2°, and 0.8%, respectively. Lastly, we calibrated a 48-channel OPM array. The accuracy of the OPM calibration was tested by using the OPM array to localize magnetic dipoles in a phantom, which resulted in an average dipole position error of 3.3 mm. The results demonstrate the feasibility of using electromagnetic coils to calibrate and localize OPMs for MEG.

47 OTHER INSTRUMENTATION↗

Convective-core Overshoot and Suppression of Oscillations: Constraints from Red Giants in NGC 6811

Using data from the NASA spacecraft Kepler , we study solar-like oscillations in red giant stars in the open cluster NGC 6811. We determine oscillation frequencies, frequency separations, period spacings of mixed modes, and mode visibilities for eight cluster giants. The oscillation parameters show that these stars are helium-core-burning red giants. The eight stars form two groups with very different oscillation power spectra; the four stars with the lowest Δ ν values display rich sets of mixed l = 1 modes, while this is not the case for the four stars with higher Δ ν . For the four stars with lowest Δ ν , we determine the asymptotic period spacing of the mixed modes, Δ P , which together with the masses we derive for all eight stars suggest that they belong to the so-called secondary clump. Based on the global oscillation parameters, we present initial theoretical stellar modeling that indicates that we can constrain convective-core overshoot on the main sequence and in the helium-burning phase for these ∼2 M {sub ⊙} stars. Finally, our results indicate less mode suppression than predicted by recent theories for magnetic suppression of certain oscillation modes in red giants.

79 ASTRONOMY AND ASTROPHYSICS↗

Massive Star Cluster Formation and Destruction in Luminous Infrared Galaxies in GOALS

We present the results of a Hubble Space Telescope ACS/HRC FUV, ACS/WFC optical study into the cluster populations of a sample of 22 Luminous Infrared Galaxies in the Great Observatories All-Sky LIRG Survey. Through integrated broadband photometry, we have derived ages and masses for a total of 484 star clusters contained within these systems. This allows us to examine the properties of star clusters found in the extreme environments of LIRGs relative to lower luminosity star-forming galaxies in the local universe. We find that by adopting a Bruzual and Charlot simple stellar population model and Salpeter initial mass function, the age distribution of the clusters declines as dN/dτ=τ{sup −0.9+/−0.3}, consistent with the age distribution derived for the Antennae Galaxies, and interpreted as evidence for rapid cluster disruption occurring in the strong tidal fields of merging galaxies. The large number of 10{sup 6} M{sub ⊙} young clusters identified in the sample also suggests that LIRGs are capable of producing more high-mass clusters than what is observed to date in any lower luminosity star-forming galaxy in the local universe. The observed cluster mass distribution of dN/dM=M{sup −1.95+/−0.11} is consistent with the canonical −2 power law used to describe the underlying initial cluster mass function (ICMF) for a wide range of galactic environments. We interpret this as evidence against mass-dependent cluster disruption, which would flatten the observed CMF relative to the underlying ICMF distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

Planet Formation by Gas-assisted Accretion of Small Solids

We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗

Quantifying vertical wind shear effects in shallow cumulus clouds over Amazonia

Abstract. This study analyses and quantifies the effects of vertical wind shear (VWS) on the properties of shallow cumulus cloud fields over Central Amazonia. We perform idealised simulations with high resolution (50 m horizontally and 20 m vertically) using the Dutch Atmospheric Large-Eddy Simulation (DALES) model, changing the initial conditions and large-scale forcing of VWS. The resulting cloud field is analysed by applying a cloud tracking algorithm to generate Lagrangian datasets of the life cycle of individual clouds as well as their time-varying core and margin dimensions. The reference run has no wind speed or directional shear and represents a typical day in the local dry season. Numerical experiments with moderate and high wind speed shear are simulated by adding linear increases in the wind speed of 1.2 and 2.4 m s−1 km−1, respectively. Three additional runs are made by adding 90∘ of wind rotation between the surface and the top of the domain (5 km) on top of the three wind speed shear conditions. We find that clouds developing in a sheared environment have horizontal equivalent diameter increased by up to 100 m on average, but the cloud depth is reduced. Our quantification shows that VWS tends to increase the size of the cloud cores but reduces their relative area, volume, and mass fractions compared to the overall cloud dimensions. The addition of 2.4 m s−1 km−1 of VWS decreases the relative core area by about 0.03 (about 10 % of the overall average) and its volume and mass ratios by about 0.05 (10 %–25 % in relative terms). Relevant for the cloud transport properties is that the updraught speed and the liquid water content are lower within the cores, and consequently so is the upward mass flux. All quantifications of mean cloud properties point to the inhibition of convective strength by VWS, therefore hampering the shallow-to-deep transition. However, open questions still remain given that the individually deepest clouds were simulated under high environmental shear, even though they occur in small numbers. This could indicate other indirect effects of VWS that have opposite effects on cloud development if found to be significant in the future.

54 ENVIRONMENTAL SCIENCES↗

Seasonal forecasting skill for the High Mountain Asia region in the Goddard Earth Observing System

Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorologically relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks to months is an area of active research and development. NASA's Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1–3-month lead times in the HMA region, including a portion of the Indian subcontinent, during the retrospective forecast period, 1981–2016. To assess forecast skill, we evaluate 2 m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to the similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from an anomaly correlation of R anom =0.18 for precipitation to R anom =0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are evaluated against independent observations; results for the 1-month forecast skill range from R anom =0.13 for snow water equivalent to R anom =0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system's ability to forecast HMA hydrometeorology.

54 ENVIRONMENTAL SCIENCES↗