Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical analysis”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Geomechanical Simulation of Big Hill Strategic Petroleum Reserve (Model Calibration)

A finite element numerical analysis model has been constructed that consists of a realistic mesh capturing the geometries of Big Hill (BH) Strategic Petroleum Reserve (SPR) site using the multi-mechanism deformation (M-D) salt constitutive model and including data taken daily of the wellhead pressure and level of the oil-brine interface. The salt creep rate is not uniform in the salt dome, and creep test data for BH salt is limited. Therefore, a model calibration is necessary to simulate the geomechanical behavior of the salt dome. Cavern volumetric closures of SPR caverns calculated from sonar survey reports are used for the field baseline measurement. The structure factor, A 2 , and transient strain limit factor, K o , in the M-D constitutive model are used for model calibration. An A 2 value obtained experimentally from the BH salt and K o value of WIPP salt are used as the baseline values. To adjust the magnitude of A 2 and K o , multiplication factors A2F and KOF are defined, respectively. The A2F and KOF values of the salt dome and salt drawdown layer of elements surrounding each SPR cavern have been determined through a number of back fitting analyses. The trendlines of the predictions and sonar data match up well for BH 101, 103, 104, 106, 110, 111, and 113. The prediction curves are close to the sonar data for BH 102 and 114. However, the prediction curves for BH 105, 107, 108, 109, and 112 are not close to the sonar data. An inconsistency was found in the sonar data, i.e. the sonar measurements of cavern volumes increase with time, during some periods for BH 101, 104, 106, 107, and 112. A follow-up report in 2019 will provide a resolution for these issues.

02 PETROLEUM↗
DOE OSTI2019

Assessment of the Available Drawdowns for Oil Storage Caverns at the Big Hill SPR Site - Cavern Integrity

This report updates the estimated values of the baseline available drawdowns for the caverns at the Big Hill storage facility, and an updated table listing the available drawdowns. A new finite element numerical analysis model was constructed that consists of a realistic mesh capturing the sonar-measured geometries of Big Hill SPR site and used the daily data of actual wellhead pressures and oil-brine interfaces. The number of available drawdowns for each of the Big Hill SPR caverns is estimated using the new model. All caverns are predicted to have five available drawdowns remaining from a geomechanical perspective. BC-101, 105, and 110 have a region of concern at the floor edge and/or sloping floor, where tensile and dilatant stresses are predicted to occur during each workover. The tensile state is predicted to occur because of the geometries of the edge and floor. Therefore, geomechanical examination for three caverns would be recommended after a drawdown leach. The well integrity of each cavern is not investigated in this report. The estimate of the number of baseline available drawdowns for the Big Hill caverns in this report will be incorporated in future assessments of the available drawdowns for all the SPR caverns. The estimates for the number of baseline available drawdowns are subject to change in the future as the knowledge of physical phenomena at the sites, and the further development of the models of geomechanical behavior at the sites, evolve over time.

02 PETROLEUM↗
DOE OSTI2022

Computational Physics Overview [Slides]

Computational physics is an important part of the overall investment in National Security Science at LANL. Computational physics is the study and implementation of numerical analysis to solve problems in physics for which a quantitative theory already exists. Historically, computational physics was the first application of modern computers in science. There are three key elements to computational physics: mathematical models of physical phenomena and conservation equations, computer codes that implement these models, and computer platforms that execute the code instructions and manipulate the data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Irradiation of MiniFuel Targets Bearing TRISO Fuel Compacts (Status Report)

Irradiation testing of MiniFuel compacts bearing tristructural isotropic (TRISO) fuel particles was performed at Oak Ridge National Laboratory (ORNL) to support the development of Kairos Power’s (KP’s) fluoride salt–cooled high-temperature reactor concept. The fuel compacts were fabricated with TRISO fuel particles of different types—including low-enriched uranium oxide, uranium carbide (LEUCO), natural uranium oxide, uranium carbide (NUCO), and low-enriched uranium dioxide (LEUO 2 )—and inserted into MiniFuel irradiation targets. Five targets were assembled and inserted in the High Flux Isotope Reactor (HFIR) for four cycles. The data collected post-irradiation will provide experimental input to validate TRISO fuel performance models for high particle power operations. This report summarizes the completion of the HFIR irradiation, the as-irradiated numerical analysis, and the post-irradiation work performed to date. This work was performed under the Nuclear Science User Facility program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Uncertainties in Local Intense Precipitation Flood Modeling

In the Task 1 report for the U.S. Nuclear Regulatory Commission (NRC) funded Local Intense Precipitation (LIP) Probabilistic Flood Hazard Assessment (PFHA) Pilot Study, the Pacific Northwest National Laboratory (PNNL) reviewed hydrologic and hydraulic modeling approaches for estimation of LIP flooding and their implementation in readily available simulation software packages (Prasad and Yuan 2020). NRC’s interest is in developing guidance for performing LIP PFHA at nuclear power plant (NPP) sites. This report describes the work done under Task 2 of the study that focused on describing uncertainties associated with LIP flood modeling. Two main types of uncertainties exist in LIP flood modeling: (1) aleatory uncertainties that arise from the inherent natural variability of the hydrometeorological system and (2) epistemic uncertainties that arise from the analysts’ incomplete knowledge of the hydrometeorological and site configuration. Principal sources for aleatory variability in LIP flood simulations at NPP sites include precipitation (i.e., magnitude, duration, and temporal distribution of LIP events), initial conditions (e.g., soil moisture content, stormwater drainage discharge, surface storage), boundary conditions (e.g., upstream discharge, downstream water surface elevations), and long-term temporal trends (e.g., climate change). This report describes approaches used for estimating the aleatory variability in precipitation including precipitation-frequency analysis, numerical weather prediction, and stochastic weather generation. Data sources that can be used to estimate aleatory variability in initial and boundary conditions are also described. Approaches to include effects of long-term trends like those from climate change into estimation of aleatory variability are summarized. Sources of epistemic uncertainty in LIP flood simulations at NPP sites include process representation (e.g., multiple approaches to represent runoff generation, stormwater drainage, and hydraulic routing), site configurations (e.g., site layout, flow features, status of temporary flood protection, blockage of drains), model resolution, and long-term temporal trends (e.g., known/planned site alterations, land-use changes at and in the vicinity of the site). The report describes alternative process representations (methods and models), particularly those implemented in the LIP flood simulation software packages reviewed in Task 1 report, and lists the associated model parameters. Approaches for estimating model parameters when surface and subsurface water exchanges occur are also described.

42 ENGINEERING↗
DOE OSTI2022

High-Power and Flexible Low-Frequency AC Transmission for Renewable Integration

This project develops a steady-state, dynamic, and production cost models of a multi-terminal, high-power, flexible low-frequency ac (LF-HVac) transmission system with grid-forming converters. Through both numerical analysis and EMTP simulation, the project goal is to demonstrate the feasibility of such an innovative and transformative LF-HVac system to support bulk integration of large-scale, cheap, and steady renewable resources from remote locations such as offshore wind farms to meet recent clean energy goals from DOE. The proposed adoption of low-frequency ac transmission aims to improve power transfer capability for remote large-scale renewable integration. In addition, the control for power electronics converters provides grid-support functions to improve the flexibility and resilience the existing 60-Hz HVac system. Future projects on this topic can be supported by both Energy Efficiency and Renewable Energy (DOEEE) and Office of Electricity Delivery & Energy Reliability (DOEOE) offices.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Loop Thermosyphon Enhanced Solar Collector

A two-phase loop thermosyphon solar collector was developed to transfer thermal energy efficiently and passively from a concentrated solar collector to a thermal desalination process. The two-phase flow in the loop thermosyphon is driven by a difference in gravitational pressure head between the liquid return line and the two-phase riser located between the evaporator and condenser. A numerical model incorporating mass, energy, and pressure drop balances was developed to predict system performance and use as a design tool. Three iterations of loop thermosyphon solar collectors were designed, fabricated, and tested to validate the numerical analysis and demonstrate system integration. Additionally, a novel solar receiver design was evaluated based on volumetric absorption of the concentrated sunlight by an optically absorptive two-phase working fluid. Life tests were performed to evaluate the working fluid thermal stability and resistance to ultraviolet degradation. A glass receiver tube and associated glass-to-metal junctions were developed. Finally, a full-scale loop thermosyphon was integrated with a parabolic trough solar collector and tested on-sun.

14 SOLAR ENERGY↗
DOE OSTI2024

Acceptance criteria for in situ surveillance of MSR materials based on thermally-loaded mechanical test articles

This report describes practices and acceptance test procedures for designing, running, and maintaining a material surveillance program in a future operating molten salt reactor. The programs described here rely on test data from passively actuated mechanical test articles inserted into critical regions of the reactor and periodically removed for out-of-reactor testing. The report defines definite acceptance procedures, based on the results of these tests, to determine whether a component can continue to operate accounting for the accumulation of environmentally-assisted mechanical damage in the component materials to date, and extrapolated out through the next inspection period. Additionally, the report describes work on a software tool implementing many of the surveillance methods and procedures described here and progress on simplified methods for inferring damage accumulation in the test articles, based on out-of-reactor thermal cycling, that do not rely on sophisticated numerical analysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elastic Changepoint Detection for Globally-indexed Functional Time Series Data with Climate Applications

Changepoint detection is a vital tool in the application of climate data analysis. Numerous types of climate observation data are most properly represented by functional time series, implying a need for accurate changepoint detection methods applicable to functional time series data. Such data taken at a global scale often contain both spatial heterogeneity and dependence as well as phase (time) misalignment. In this report, we present methods which can detect spatially-dependent changepoints while allowing different estimates of change time and change strength depending on location. Additionally, we provide extensions to this spatially-predicted model which controls for phase variability among observations. Our methods provide the ability to detect a single change, or control for epidemic changes (where a “return-to-normal” change is more likely to be detected than the initial change). We showcase results analyzing the June 1991 eruption of Mt. Pinatubo, where our methods demonstrate the ability to accurately detect both single and epidemic changepoints even in the presence of strong seasonal variability. We find that our spatially-predicted model improves the detection of relevant changepoints versus methods which do not take spatial information into account, and we find that controlling for phase variability helps to control the false discovery rate during the detection process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

6G Sparse Transmitarray

This powerpoint presents a sparse transmitarray design for 6G technology. The proposed approach leverages aperiodic arrays with large inter-element spacings, offering grating lobe-free operation, enhanced gain, and narrow beamwidths. The transmitarray topology, designed for sub-THz and THz frequencies, simplifies feed-network complexity, and minimizes losses. Numerical analysis and full wave simulation results show this approach is feasible.

Christian, Thomas Edward↗
DOE OSTI2024

Evaluation of Additively Manufactured Monolithic SiC and SiC Ceramic Matrix Composites for Concentrating Solar Receiver Applications

Supported by an award from the Solar Technology Office, US Department of Energy, GE Aerospace Research in collaboration with Heliogen Holdings Inc, is engaged in the development of ultra-High Operating Temperature SiC-matrix Solar Thermal Air Receiver (HOTSSTAR). We report the results of our study of the effective heat transfer characteristic of several candidate structure motifs, or feature geometries, made of SiC via additive manufacturing in a simulative laboratory test. The SiC structure motifs studied include different permutations of three-dimensional periodic lattices and defined shapes. The solar-thermal simulating laboratory test setup is constructed using a 4kW CO2 laser system with beam shaping optics to apply radiative heating power on one face of 2”-diameter cylindrical feature specimens representing the structure motifs of interest for receiver element design. Using the test setup, simulative test conditions representative of a concentrated solar flux of up to ~2000 suns could be achieved in the lab tests under varying air flow through the test structure. A numerical analysis scheme is developed to extract an effective or compound heat transfer coefficient representative of the test structure under steady-state heat flow conditions.

14 SOLAR ENERGY↗
DOE OSTI2024

Evaluation of Additively Manufactured Monolithic SiC and SiC-Matrix Ceramic Matrix Composites for Concentrating Solar Receiver Applications - Fabrication and Testing of Receiver Design Feature Specimens in a Simulating Lab Test via Laser Heating

Increasing operating temperatures of solar receivers is paramount to the efficiency of concentrated solar thermal (CST) and solar power (CSP) systems. Owing to its high temperature stability combined with excellent thermal and optical properties, SiC has been the material of choice for application in high-temperature solar receivers. We report the results of our study of the effective heat transfer characteristics of several candidate SiC structure motifs, or feature geometries, which are fabricated via additive manufacturing. The SiC structure motifs studied include different permutations of three-dimensional periodic lattices and defined shapes. A solar-thermal simulating laboratory test setup is constructed using a 4kW CO2 laser system with beam shaping optics to apply concurrent radiative heating power on one face of 2”-diameter cylindrical feature specimens, representing the structure motifs of interest for receiver element design, while flowing through the sample as heat transfer fluid. Using the test setup, simulative test conditions representative of a concentrated solar flux of up to ~2000 suns could be achieved in the lab tests under varying air flow through the test structure. A simple 1D numerical analysis scheme is developed to extract an effective or compound heat transfer coefficient representative of the test structure under steady-state heat flow conditions. The test results and their use to guide the selection and optimization of SiC material and structure motifs for the receiver element design fabrication are discussed.

14 SOLAR ENERGY↗
DOE OSTI2021

A Novel Feature of Valence Quark Distributions in Hadrons

Examining the evolution of the maximum of valence quark distribution, qV, weighted by Bjorken x, h(x,t)≡xqV(x,t), it is observed that h(x,t) at the peak becomes a one-parameter function; h(xp,t)=Φ(xp(t)), where xp is the position of the peak, t=logQ2, and Q2 is the resolution scale. This observation is used to derive a new model-independent relation which connects the partial derivative of the valence parton distribution functions (PDFs) in xp to the quantum chromodynamics (QCD) evolution equation through the xp derivative of the logarithm of the function Φ(xp(t)). A numerical analysis of this relation using empirical PDFs results in an observation of the exponential form of the Φ(xp(t))=h(xp,t)=CeDxp(t) for leading to next-to-next leading order approximations of PDFs for the range of Q2, covering four orders in magnitude. The exponent, D, of the observed “height-position” correlation function converges with the increase in the order of approximation. This result holds for all the PDF sets considered. A similar relation is observed also for the pion valence quark distribution, indicating that the obtained relation may be universal for any non-singlet partonic distribution. The observed “height-position” correlation is used also to indicate that no finite number of exchanges can describe the analytic behavior of the valence quark distribution at the position of the peak at fixed Q2.

Leon, Christopher↗
DOE OSTI2017

Selective Mode Focusing in a Plate of Arbitrary Shape Applying Time Reversal Mirrors

In this study, a time reversal mirror is used to remotely focus symmetric or antisymmetric modes in a plate of arbitrary shape without the need of precise knowledge about material properties and geometry. The addition or subtraction of the forward motions recorded by two laser beams located on both sides of the plate allows, respectively, to focus a symmetric or an antisymmetric mode. The concept is validated using experimental and numerical analysis on an aluminum plate of complex machined geometry which exhibits various thicknesses as well as a bi-materials zone. Finally, the limitations and possible ways to overcome them are then presented.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Gas Atomization of Molten Metal: Part I. Numerical Modeling Conception

This numerical analysis study entails creating and assessing a model that is capable of simulating molten metal droplets and the production of metal powder during the Gas Atomization (GA) method. The essential goal of this research aims to gather more information on simulating the process of creating metal powder. The model structure and perspective was built through the application of governing equations and aspects that utilized factors such as gas dynamics, droplet dynamics, energy balance, heat transfer, fluid mechanics and thermodynamics that were proposed from previous studies. The model is very simple and can be broken down into having a set of inputs to produce outputs. The inputs are the processing parameters such as the initial temperature of the metal alloy, the gas pressure and the size of the droplets. Additional inputs include the selection of the metal alloy and the atomization gas and factoring in their properties. The outputs can be designated by the velocity and thermal profiles of the droplet and gas. These profiles illustrate the speed of both as well as the rate of temperature change or cooling rate of the droplets. Here, the main focus is the temperature change and finding the right parameters to ensure that the metal powder is efficiently produced. Once the model was conceptualized and finalized, it was employed to verify the results of other previous studies.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Fully Kinetic Large-scale Simulations of the Collisionless Magnetorotational Instability

We present two-dimensional particle-in-cell simulations of the fully kinetic collisionless magnetorotational instability (MRI) in weakly magnetized (high β) pair plasma. The central result of this numerical analysis is the emergence of a self-induced turbulent regime in the saturation state of the collisionless MRI, which can only be captured for large enough simulation domains. One of the underlying mechanisms for the development of this turbulent state is the drift-kink instability (DKI) of the current sheets resulting from the nonlinear evolution of the channel modes. The onset of the DKI can only be observed for simulation domain sizes exceeding several linear MRI wavelengths. The DKI and ensuing magnetic reconnection activate the turbulent motion of the plasma in the late stage of the nonlinear evolution of the MRI. At steady-state, the magnetic energy has an MHD-like spectrum with a slope of k –5/3 for kρ < 1 and k –3 for sub-Larmor scale (kρ > 1). We also examine the role of the collisionless MRI and associated magnetic reconnection in the development of pressure anisotropy. We study the stability of the system due to this pressure anisotropy, observing the development of mirror instability during the early-stage of the MRI. We further discuss the importance of magnetic reconnection for particle acceleration during the turbulence regime. Furthermore, consistent with reconnection studies, we show that at late times the kinetic energy presents a characteristic slope of ϵ –2 in the high-energy region.

79 ASTRONOMY AND ASTROPHYSICS↗