SEARCH · Engineering Papers

Results for “load uncertainty”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis – Simulated Wind

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wind turbine. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen . While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the simulated wind energy profiles, NLR used OpenFAST to simulate a 3.4-MW International Energy Agency (IEA) reference wind turbine. The hour-long wind energy profiles varied over wind turbulence intensity (Class A or Class C) and average wind speed (5, 7, or 9 m/s). To match the power limits of the 1.25-MW electrolyzer and 3.4-MW IEA wind turbine most effectively and to maximize the efficiency of hydrogen production at a given average wind speed, the profiles were sometimes scaled by two times. This means that, in some cases, the experimental setup assumed two 1.25-MW electrolyzers were coupled with the wind turbine, representing a total maximum electrolysis load of 2.5 MW. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wind turbine electrolysis experiment and is formatted as follows: {technology}-{average wind speed}-{turbulence class}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “windIEA3.4-5ms-C_2-400.zip” represents the hour-long experiment using the IEA 3.4-MW turbine, subjected to an average wind speed of 5 m/s and Class C wind turbulence, and connected to two 1.25-MW electrolyzers with the power supply set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wind turbine power. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all simulated wind experiments combined into one dataset labeled "combined_wind_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis .

08 HYDROGEN↗

DOE OSTI2020

Pumped Storage Hydropower (PSH) FAST Commissioning Prize Technical Analysis

The US energy landscape has undergone major changes over the past 10 years and will continue to see significant changes in future decades as the power grid increases its reliance on variable renewable energy resources. Because of the inherent variability of these resources, renewable energy growth may require additional energy storage capacity to provide flexible load-following capabilities and other grid services that can quickly adjust to changes in energy demand and generation. Pumped storage hydropower (PSH)—one such energy storage technology—uses pumps to convey water from a lower reservoir to an upper reservoir for energy storage and releases water back to the lower reservoir via a powerhouse for hydropower generation. PSH facility pump and generation cycling often follows economic and energy demand conditions. Across the United States, 43 PSH facilities are in operation and 55 projects are in various permitting or licensing stages. Altogether, the 43 operational projects provide the wide majority (95%) of utility-scale electricity storage in the United States. These facilities also provide significant power and nonpower grid benefits, including large-scale electrical system reserve capacity, grid reliability support, and electricity supply-demand balancing through quick-response capabilities and operational flexibility. PSH systems can accomplish these at a scale (e.g., size) and cost that makes these systems highly attractive from a technical standpoint. Although these research concepts are still in their infancy, they demonstrate promising potential as future PSH energy storage technologies. Although PSH has many advantages, development in the United States has effectively stalled since the 1990s, partially because of the magnitude of project costs and financing interest during development and construction, the length of time from project investment until project revenue begins, permitting challenges, construction risks, competition from other storage technologies (e.g., batteries, hydrogen storage), and electricity market evolution and uncertainty. In short, the time, cost, and risk associated with modern PSH development have resulted in limited growth in the United States recently, despite the growing energy storage demand stemming from increased wind and solar power deployment. Technology innovation is needed to help reduce PSH commissioning time, cost, and risk, particularly during the post-licensing phase of project development. To address challenges facing the PSH industry and to improve PSH commissioning timelines, the US Department of Energy (DOE) Water Power Technologies Office (WPTO) initiated the PSH Furthering Advancements to Shorten Time to (FAST) Commissioning Prize project.

13 HYDRO ENERGY↗

NASA NTRS2018

A13H-2546: Intercomparison of biomass burning aerosol properties from in-situ and remote-sensing instruments in ORACLES-2016

The magnitude of the cumulative effects of aerosols, both directly and on cloud properties, remains the biggest source of uncertainty in our understanding of anthropogenic radiative forcing on the climate. The magnitude of aerosol effects in any given region will depend on not just the aerosol concentration, but also on the characteristics of the aerosol under consideration. The NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) campaign was conducted in 2016-2018 over three deployments in the Southeast Atlantic Ocean (SEA), with the goal of making observations of seasonal biomass burning (BB) smoke plumes overlying persistent stratocumulus cloud decks. During the first ORACLES deployment (September 2016), two NASA aircraft (a P-3 and an ER-2) were flown with a suite of aerosol, cloud, radiation, and meteorological instruments for remote-sensing and in-situ observations. In ORACLES, data collected by eight different instrument teams are used to derive aerosol properties over this region by independent methods. Here we present a comparison of results from the different instrument retrievals with a focus on the measures of aerosol absorption. We focus on specific cases with near-coincident measurements from in-situ (PSAP or PTI for aerosol absorption combined with a nephelometer for aerosol scattering) and remote sensing (4STAR, AirMSPI, RSP, and SSFR) instruments from coordinated flights between the two aircraft. We use the single scattering albedo (SSA), a measure of the relative absorption versus total extinction by aerosols, as the main comparable property between the various remote sensing and in-situ methods. For the remote sensing instruments, we see general good agreement in the total aerosol loading (aerosol optical depth, AOD), though see some differences in the absorbing component (AAOD). We report values for the retrieved aerosol parameters measured in the ORACLES flights and comment upon the observed variability between different instruments. Finally, we discuss these ORACLES results in the context of previous measurements of BB aerosol over the SEA.

Kristina Pistone↗

NASA NTRS2014

Retrieval, Inter-Comparison, and Validation of Above-Cloud Aerosol Optical Depth from A-train Sensors

Absorbing aerosols produced from biomass burning and dust outbreaks are often found to overlay lower level cloud decks and pose greater potentials of exerting positive radiative effects (warming) whose magnitude directly depends on the aerosol loading above cloud, optical properties of clouds and aerosols, and cloud fraction. Recent development of a 'color ratio' (CR) algorithm applied to observations made by the Aura/OMI and Aqua/MODIS constitutes a major breakthrough and has provided unprecedented maps of above-cloud aerosol optical depth (ACAOD). The CR technique employs reflectance measurements at TOA in two channels (354 and 388 nm for OMI; 470 and 860 nm for MODIS) to retrieve ACAOD in near-UV and visible regions and aerosol-corrected cloud optical depth, simultaneously. An inter-satellite comparison of ACAOD retrieved from NASA's A-train sensors reveals a good level of agreement between the passive sensors over the homogeneous cloud fields. Direct measurements of ACA such as carried out by the NASA Ames Airborne Tracking Sunphotometer (AATS) and Spectrometer for Sky-Scanning, Sun-Tracking Atmospheric Research (4STAR) can be of immense help in validating ACA retrievals. We validate the ACA optical depth retrieved using the CR method applied to the MODIS cloudy-sky reflectance against the airborne AATS and 4STAR measurements. A thorough search of the historic AATS-4STAR database collected during different field campaigns revealed five events where biomass burning, dust, and wildfire-emitted aerosols were found to overlay lower level cloud decks observed during SAFARI-2000, ACE-ASIA 2001, and SEAC4RS- 2013, respectively. The co-located satellite-airborne measurements revealed a good agreement (RMSE less than 0.1 for AOD at 500 nm) with most matchups falling within the estimated uncertainties in the MODIS retrievals. An extensive validation of satellite-based ACA retrievals requires equivalent field measurements particularly over the regions where ACA are often observed from satellites, i.e., south-eastern Atlantic Ocean, tropical Atlantic Ocean, northern Arabian Sea, South-East and North-East Asia.

validation↗

NASA NTRS2020

Design and Analysis of Buckling-Critical Large-Scale Sandwich Composite Cylindrical Test Articles

It has long been established in the literature that the buckling response of thin-shell structures can be very sensitive to the presence of small geometric and loading imperfections. The Shell Buckling Knockdown Factor Project (SBKF) was established by the NASA Engineering and Safety Center (NESC) to develop analysis-based shell buckling design recommendations for stiffened-metallic and composite launch-vehicle shell structures. Large-scale buckling tests were used to validate the modeling and analysis methods applied in developing these analysis-based recommendations. Herein, the test article design methodology for 8-ft-diameter, honeycomb-core sandwich composite cylinder validation tests is discussed and cylinder designs are presented. In this methodology, first, the sandwich composite design space was defined using several nondimensional parameters, and the desired test article design space was determined by examining the designs of launch-vehicle cylinder structures. Essentially all test article designs within certain design parameters were generated and then downselected based on simple closed-form failure calculations and the nondimensional design-space parameters. Four of these designs that spanned a significant portion of the design space of interest and had global buckling as the first predicted failure mode were selected and subjected to higher-fidelity finite element analyses (FEAs): shell-element-based analyses, axisymmetric-element-based analyses, and global-local analyses. The analysis flow discussed in this report supported the design objective. As the analysis flow progressed, designs were downselected so the fidelity of the analysis methods, and consequently their computational cost and accuracy, was increased. The selection of the FEA types created an analysis framework where particular methods complemented each other and reduced the uncertainty of the predicted test article responses. The analysis results are illustrated using several designs when the computationally expeditious closed-form analysis stage is discussed. Once this stage is complete, the higher-fidelity FEA types are illustrated using one selected detailed test article design. Both perfect and imperfect test article geometries were considered.

Buckling↗

DOE OSTI2021

Fracture toughness evaluations for spent nuclear fuel dry storage canister welds and spent nuclear fuel clad-pellet structures

Integrity of spent nuclear fuel (SNF) interim storage canisters is very important to the safety of the back-end nuclear fuel cycle. Stress corrosion cracking (SCC) potential of interim storage canister has been considered as a high priority. Because no post-weld heat treatment was required for forming these canisters, the high tensile residual stress existed within these canister welds. This can change the fracture resistance capacity significantly as well as increase SCC potential. Due to relative thin shell thickness of a canister weldment, the spiral notch torsion test (SNTT) method was used to estimate the canister weldment fracture toughness. SNTT was developed to measure the intrinsic fracture toughness of structural materials using small specimens. The SNTT method has been applied to a wide variety of structural materials, such as low-alloy steels, stainless steel, aluminum alloy, ceramics, concrete, and composites. The SNTT system operates by applying pure torsion to cylindrical specimens with a notch line that spirals around the specimen at a 45° pitch. In order to carry out pure torsion load mode, biaxial tension/torsion tester was developed accordingly to perform SNTT protocol. Moreover, applying fracture mechanics approach to SNF system reliability investigation is warranted due to the inherent flaws and hydride structures existed in a SNF system after nuclear reactor operation. However, none of the existing fracture toughness data deal with fuel cladding specific geometry or spent fuel material conditions, such as cladding structure with the pellet-inserts and the associated pellet clad mechanical interactions induced mixed-mode damage mechanisms. Thus, the development of an intrinsic fracture mechanics approach that is suitable for SNF materials is needed. Furthermore, due to thin wall and small dimension of clad tubing structure, the SNTT method was used to estimate the clad tubing structure fracture toughness. Fracture testing were performed on the received stainless steel canister weldment, most SNTT weld samples fracture initiation sites are at heat-affected zone (HAZ) regions. The estimated fracture toughness J Q ’ for the baseline SS304 steel is at 283 kJ/m². The estimated JQ’ for the SS304/308 weld and baseline metals are 148 kJ/m 2 and 459 kJ/m 2 , respectively. Out of cell fracture testing for spent fuel structure were carried out on the surrogate rods made of Zr-4 clad and alumina inserts, the estimated fracture toughness values for baseline Zr-4 cladding with alumina-pellet inserts are: (1) For SNTT samples with a short or medium crack length, between 5.4-mm and 8-mm, the estimate J IQ upon fracture initiation for the baseline Zr-4 cladding is at 50 kJ/m 2 with 2-sigma uncertainty of 3.26 kJ/m 2 , and the associated K IQ is at 67.46 MPa$\sqrt{m}$; and (2) For SNTT samples with a long crack length, around 13-mm, the crack initialization is deviated from that of the Mode-I tensile fracture and appears to be a mixed-mode fracture of Mode I - tensile stress and Mode III - out of plane shear stress; the estimated J MQ is at 18.9 kJ/m 2 , the associated K MQ is at 41.4 MPa$\sqrt{m}$.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2000

Experimental and Numerical Characterization of a Cylindrical Blackbody Cavity

During hypersonic flight, high temperatures and high heat fluxes are generated on the surfaces of vehicles. The Flight Loads Laboratory (FLL) at Dryden Flight Research Center (DFRC) is equipped with a calibration furnace, capable of calibrating heat flux gages up to 1100kW per square meters, and temperature sensors up to 2600 C. One heating configuration of the calibration furnace is a cylindrical blackbody cavity. Throughout the blackbody there are temperature gradients due to various boundary conditions. These boundary conditions include resistance heating, radiant heat transfer, and conduction to water-cooled electrodes. Also, an inert gas is purged through the graphite blackbody to prevent it from oxidizing. Consequently, the various modes of heat transfer present during operation of the blackbody cavity must be well understood in order to produce accurate heat flux gage and temperature sensor calibrations for use in ground testing or flight testing of hypersonic vehicles. The first step towards understanding the heat transfer in the blackbody cavity was to perform experiments at 1100 C, with and without outer surface insulation, while taking detailed temperature measurements inside the blackbody cavity. Steady state thermal models of the blackbody cavity were then developed. These models included detailed thermal analysis using commercial thermal analysis software. Conduction, radiation, and convection were considered in the thermal models for two cases: one with the outside of the blackbody cavity insulated and the second without insulation. This paper describes the experimental and numerical efforts used to characterize the steady state operation of the blackbody cavity. It describes the analysis of the test measurements, the boundary conditions used in the numerical models, and how the models were calibrated to fit the experimental data. Effects of various uncertainties, such as material properties, and convection are discussed.Initial thermal models predicted temperatures in the deepest part of the blackbody cavity within 7 C of the measured value and produced trends comparable to the experimental data, throughout the models. Adjustment of the boundary conditions, which were included in the thermal models, produced good agreement with measured temperatures. Free and forced convection of the purge gas inside the blackbody was found to be insignificant.

Abdelmessih, Amanie N.↗

NASA NTRS1994

Design of a Flush Airdata System (FADS) for the Hypersonic Air Launched Option (HALO) Vehicle

This paper presents a design study for a pressure based Flush airdata system (FADS) on the Hypersonic Air Launched Option (HALO) Vehicle. The analysis will demonstrate the feasibility of using a pressure based airdata system for the HALO and provide measurement uncertainty estimates along a candidate trajectory. The HALO is a conceived as a man-rated vehicle to be air launched from an SR-71 platform and is proposed as a testbed for an airbreathing hydrogen scramjet. A feasibility study has been performed and indicates that the proposed trajectory is possible with minimal modifications to the existing SR71 vehicle. The mission consists of launching the HALO off the top of an SR-71 at Mach 3 and 80,000 ft. A rocket motor is then used to accelerate the vehicle to the test condition. After the scramjet test is completed the vehicle will glide to a lakebed runway landing. This option provides reusability of the vehicle and scramjet engine. The HALO design will also allow for various scramjet engine and flowpath designs to be flight tested. For the HALO flights, measurements of freestream airdata are considered to be a mission critical to perform gain scheduling and trajectory optimization. One approach taken to obtaining airdata involves measurement of certain parameters such as external atmospheric winds, temperature, etc to estimate the airdata quantities. This study takes an alternate approach. Here the feasibility of obtaining airdata using a pressure-based flush airdata system (FADS) methods is assessed. The analysis, although it is performed using the HALO configuration and trajectory, is generally applicable to other hypersonic vehicles. The method to be presented offers the distinct advantage of inferring total pressure, Mach number, and flow incidence angles, without stagnating the freestream flow. This approach allows for airdata measurements to be made using blunt surfaces and significantly diminishes the heating load at the sensor. In the FADS concept a matrix of flush ports is placed in the vicinity of the aircraft nose, and the airdata are inferred indirectly from the measured pressures.

Whitmore, Stephen A.↗

NASA NTRS2008

Magnetic Compensation for Second-Order Doppler Shift in LITS

The uncertainty in the frequency of a linear-ion-trap frequency standard (LITS) can be reduced substantially by use of a very small magnetic inhomogeneity tailored to compensate for the residual second-order Doppler shift. An effect associated with the relativistic time dilatation, one cause of the second-order Doppler shift, is ion motion that is attributable to the trapping radio-frequency (RF)electromagnetic field used to trap ions. The second-order Doppler shift is reduced by using a multi-pole trap; however it is still the largest source of systematic frequency shift in the latest generation of LITSs, which are among the most stable clocks in the world. The present compensation scheme reduces the frequency instability of the affected LITS to about a tenth of its previous value. The basic principles of prior generation LITSs were discussed in several prior NASA Tech Briefs articles. Below are recapitulated only those items of basic information necessary to place the present development in context. A LITS includes a microwave local oscillator, the frequency of which is stabilized by comparison with the frequency of the ground state hyperfine transition of 199Hg+ ions. The comparison involves a combination of optical and microwave excitation and interrogation of the ions in a linear ion trap in the presence of a nominally uniform magnetic field. In the current version of the LITS, there are two connected traps (see figure): (1) a quadrupole trap wherein the optical excitation and measurement take place and (2) a 12-pole trap (denoted the resonance trap), wherein the microwave interrogation takes place. The ions are initially loaded into the quadrupole trap and are thereafter shuttled between the two traps. Shuttling ions into the resonance trap allows sensitive microwave interrogation to take place well away from loading interference. The axial magnetic field for the resonance trap is generated by an electric current in a finely wound wire coil surrounded by magnetic shields. In the quadrupole and 12-pole traps, the potentials are produced by RF voltages applied to even numbers (4 and 12, respectively) of parallel rods equally spaced around a circle. The polarity of the voltage on each rod is opposite that of the voltage on the adjacent rod. As a result, the amplitude of the RF trapping field is zero along the centerline and increases, with radius, to a maximum value near the rods.

Burt, Eric↗

NASA NTRS2003

Impact of Climatic Variability on Atmospheric Mass Distribution and GRACE-Derived Gravity Fields

During the period we calculated the atmospheric data sets related to its mass and angular momentum distribution. For mass, we determined the various harmonics from the NCEP-NCAR reanalysis, especially the low-order harmonics that are useful in studying the gravitation distribution as will be determined from the GRACE mission. Atmospheric mass is also related to the atmospheric loading on the solid Earth; we cooperated with scientists who needed the atmospheric mass information for understanding its contributions to the overall loading, necessary for vertical and horizontal coordinate estimation. We calculated atmospheric angular momentum from the NCEP-NCAR reanalyses and 4 operational meteorological centers, based on the motion (wind) terms and the mass (surface pressure) terms. These are associated with motions of the planet, including its axial component causing changes in the length of day, more related to the winds, and the equatorial component related to motions of the pole, more related to the mass. Tasks related to the ocean mass and angular momentum were added to the project as well. For these we have noted the ocean impact on motions of the pole as well as the torque mechanisms that relate the transfer of angular momentum between oceans and solid earth. The activities of the project may be summarized in the following first manuscript written in December 2002, for a symposium that Dr. Salstein attended on Geodynamics. We have continued to assess ocean angular momentum (OAM) quantities derived from bottom pressure and velocity fields estimated with our finite-difference barotropic (single layer) model. Three years of output (1993-95) from a run without any data constraints was compared to output from a corresponding run that was constrained by altimeter data using a Kalman filter and smoother scheme. Respective OAM time series were combined with corresponding atmospheric series and compared to observed polar motion. The constrained OAM series provided slightly better variance reduction than the unconstrained series. Analysis provided a check on the estimation scheme and pointed to further work to improve the determination of OAM using this method. A significant effort was also devoted to quantifying effects of uncertainties in high frequency winds on the mean and seasonal momentum exchange between atmosphere and oceans.

Salstein, David A.↗

DOE OSTI2020

Using the Carbon Capture Simulation Initiative (CCSI) Tool to Design the Experiments in the Parametric Campaign of a Novel Compact Absorber for Carbon Capture

Gas absorption towers with structured packing and solvent have been used for Carbon Dioxide (CO 2 ) Capture for about many decades. To overcome process limitations and practical disadvantages for CO 2 capture from the stationary emitter (e.g. NG and coal power plant), many new designs have been proposed and explored at the various scales in the last decade with aim of either low energy penalty or low capital cost. To reduce the size of the absorption tower and hence the total cost of CO 2 capture, the University of Kentucky Center for Applied Energy Research Center (UK CAER) has designed and built a novel CO 2 capture absorption tower or Compact Absorber, integrated into an existing large-bench scale CO 2 capture unit. The Compact Absorber has three sections. The top of the column is a fogging section where the solvent is sprayed through a nozzle producing droplets flowing downward in a co-current fashion with the flue gas. The center of the column is a frothing section where the solvent and flue gas flow through regenerative frothing screens designed by Industrial Climate Solutions, Inc. The bottom of the column is a typical structured packing section were the flue gas and solvent flow in a counter-current fashion. The parametric campaign will be conducted in order to optimize the operating parameters for CO 2 capture including liquid/gas ratio, lean loading, and temperature, liquid residence time. A simulated flue gas with 14% CO 2 will be used along with a UK CAER developed proprietary solvent. The 100-hour parametric campaign is designed using a statistical approach of the Sequential Design of Experiments (sDOE). sDOE is one of the CCSI tools that provides an adaptive statistical approach for designing future experiments based on the results of previous experiments. Application of a typical DOE provides the user with the minimum number of experiments required to get the same data, but sDOE allows the user to make an informed choice of experiments based on the results of previous experiments. The complete absorption column has been constructed and has been partially commissioned. Initial data has been collected by operating using the fogging section and the frothing section. The fogging section produces solvent droplets of about 100 μm sauter mean diameter and as small as 25 μm using a hydraulic nozzle by BETE. The frothing section produces bubbles of about 5mm with high mixing of solvent promoting the higher mass transfer from gas to liquid. The absorber reaches the capture efficiency of about 50% with only two sections in operation. Based on the current results, it can be deduced that increasing the solvent feed temperature and including the packed section for absorption the capture efficiency will increase further. Initial data will be collected using all three sections of the absorber and will be used for sDOE. Non-Uniform Space Filling model of sDOE will be used to prioritize the input conditions resulting into maximum capture efficiency. sDOE is performed using the platform called Framework Optimization, Quantification of Uncertainty, and Surrogates (FOQUS). The method and results demonstrating the progress of the parametric campaign from the initial set of experiments to the final stage of obtaining optimized parameters using sDOE tool will be presented in detail.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗

DOE OSTI2023

Subcomponent Validation of Composite Joints for the Marine Energy Advanced Materials Project

The Marine Energy Advanced Materials project is an ongoing multi-year, multi-lab project with the main goals of addressing barriers and uncertainties facing marine energy developers in adopting advanced materials for structural applications. NREL's goals of the project were to address subcomponents testing needs for marine energy materials, to improve understanding of design allowables at the full-scale and provide near net-scale static and fatigue data of composite subcomponents using materials applicable to the marine energy industry. In the long term, the test method development and data generated would be used to inform standards development. This report outlines perhaps one of the largest-scale studies conducted with regards to saltwater conditioning of various composite material subcomponents and their subsequent structural validation, specifically directed at the marine renewable energy industry. A variety of fiberglass composite panels with epoxy and vinyl ester epoxy resin systems were manufactured at Montana State University, which were then used to manufacture an array of different types of subcomponent test specimens at the National Renewable Energy Laboratory's Flatirons Campus. These subcomponents were in the form of T-bolt and double-ended-insert specimens, which were intended to represent bonded and mechanical bolted connections for thick composite laminates, metal-metal and composite lap shear specimens to evaluate adhesion of constituent materials, and adhesive beam-shear specimens as part of an effort to better evaluate the characteristics of thick adhesive bondlines. Overall, the materials used were fiberglass reinforced epoxy and vinyl ester matrix composites, epoxy and methacrylate adhesives, and 316 and 2507 stainless steels. Specimens were then conditioned in salt water at various temperatures and for various periods of time at Florida Atlantic University and Pacific Northwest National Laboratory. All specimens were then mechanically characterized and validated using various test methods under static and fatigue loading conditions at NREL's Structural Technology Laboratory. Throughout the conditioning and mechanical validation process, valuable experience was gained, which will help guide future test method development for marine energy materials. In many instances, the results indicated similar observations as to what had been observed during previous coupon scale characterization efforts that provided a vital understanding of the scale up process. However in some instances, unexpected phenomena were observed, such as interactions between the adhesives and 316 steel. Furthermore, some materials exhibited significant degradation due to the saltwater conditioning. Ultimately, this report provides a detailed summary of the specimens that were designed, the subcomponent test methods that were developed, and the results that were generated, which will serve as important guidance for marine renewable energy developers and researchers for future structural designs and validation.

16 TIDAL AND WAVE POWER↗

NASA NTRS2010

Orion Entry Performance-Based Center-of-Gravity Box

The Orion capsule is designed both for Low Earth Orbit missions to the ISS and for missions to the moon. For ISS class missions, the capsule will use an Apollo-style direct entry. For lunar return missions, depending on the timing of the mission, the capsule could perform a direct entry or a skip entry of up to 4800 n.mi. in order to land in the coastal waters of California. The physics of atmospheric re-entry determine the capability of the Orion vehicle. For a given vehicle mass and shape, physics tells us that the driving parameters for an entry vehicle are the hypersonic lift-to-drag ratio (L/D) and the flight path angle at entry interface (gamma(sub EI)). The design of the Orion atmospheric re-entry must meet constraints during both nominal and dispersed flight conditions on landing accuracy, heating rate, total heat load, sensed acceleration, and proper disposal of the Service Module. These constraints define an entry corridor in the space of L/D-gamma(sub EI); if the vehicle falls within this corridor, then all constraints are met. The gamma(sub EI) dimension of the corridor can be further constrained by the gloads experienced during emergency entries. Thus, the entry performance for the Orion vehicle can be described completely by the L/D. Bounds on the hypersonic L/D necessary to achieve all the mission requirements can be defined for the given entry corridor. Landing accuracy performance drives the lower limit on L/D. In order to achieve the desired landing accuracy, a minimum L/D must be ensured. The design of the Thermal Protection System (TPS) drives the upper limit on L/D. A higher L/D can drive mass into the design of the TPS. Conversely, once the TPS is designed, the L/D must be ensured to stay below a certain limit in order for the TPS to stay within its design envelop. The L/D must stay within its upper and lower bounds during dispersed flight conditions. L/D is a function of both the aerodynamics and the center-of-gravity (CG) of the vehicle. The aerodynamics of the vehicle are determined by Computational Fluid Mechanics (CFD) and wind tunnel tests. However, the aerodynamics are not known precisely. Instead, an aerodynamic database has been developed where the aerodynamic coefficients are known to fall within a probabilistic band defined by upper and lower bounds. It is expected that the probabilistic band will shrink after the first missions are flown and real-world data is collected. Until that time, the Orion must be designed to the current aerodynamic database. Thus, for a given aerodynamic database with given uncertainties, the allowable range in L/D can be mapped to an allowable box for the CG location. The CG box is used to set requirements on the dispersions allowed for vehicle packaging and cargo storage. As the aerodynamic uncertainties decrease, the size of the CG box can increase. This paper discusses the technique used to map the minimum and maximum L/D bounds set by the entry performance requirements to the allowable dispersions in CG while accounting for aerodynamic uncertainties. The L/D is defined as the ratio of the lift force to the drag force. It is equivalent to the ratio of lift coefficient (C(sub L)) over drag coefficient (C(sub D)). C(sub L) and C(sub D) are functions of Mach number (M) and angle of attack (alpha). A Mach number of 25 is used as a measuring point of the hypersonic L/D. Variations in C(sub L), C(sub D) and alpha cause variations in L/D. Equation (1) shows the three contributions to the variation in L/D.

Rea, Jeremy R.↗

DOE OSTI2022

Small field campaign: aerosol – ice formation closure pilot study. Final Report

Prediction of atmospheric ice formation from aerosol particles by heterogeneous nucleation represents one of the grand challenges in atmospheric science. Our insufficient predictive understanding of primary ice formation is the reason that climate models typically do not include heterogeneous ice nucleation with subsequent effects on climate uncertainty. Mixed-phase clouds, where supercooled water droplets and ice crystals coexist play globally an important role regulating climate. This is especially the case for the Arctic region that experiences the greatest warming due to climate change compared to other regions in the world. Immersion freezing initiated by ice-nucleating particles (INPs) in supercooled water droplets is recognized as the dominant primary ice formation pathway in mixed-phase cloud regimes. For this reason, it is crucial to evaluate our capability to predict immersion freezing for a given ambient aerosol population. The goal of this project is conducting a field-based pilot study at the U.S. DOE Atmospheric Radiation Measurement (ARM) user facility at Southern Great Plains (SGP) to evaluate our capability to predict the number concentration of aerosol particles that serve as INPs in the immersion freezing mode. Successful prediction of INP number concentrations is also termed “closure”. This field-observational approach represents a first-of-its kind attempt of an aerosol–ice formation closure study (AEROICESTUDY). Very few closure studies related to INPs have been conducted, and to our knowledge, none using robust size-resolved ambient aerosol composition measurements as a starting point. Achievement of aerosol–ice formation closure relies on our ability to characterize the ambient aerosol population with respect to particles size and composition and to determine INP number concentrations for specified freezing temperatures. This requires numerous online and offline instrumentation resulting in this pilot field campaign being a multi-institutional and community-collaborative effort. We chose the ARM SGP megasite for this first aerosol-ice formation closure pilot study due to its significant measurement capabilities available to obtain detailed physical characterization of the local aerosol population including size distribution, mass loading, and chemical composition of non-refractory aerosol particles. The overall objective of this project is to identify ice nucleation parameterizations that produce the most robust predictions of INP numbers and thus are best suited to be included in cloud and climate models. This objective includes the following goals for this field and laboratory-based project: i) What are the crucial aerosol physicochemical property measurements needed to accurately guide ice nucleation representations in models and long-term INP measurements? ii) What level of parameter details needs to be known to achieve aerosol–ice formation closure? iii) What are the leading causes for climate model bias in INP predictions? We found that the advances in our understanding of immersion freezing garnered over the last 20 years allowed us to yield partial and full closures of atmospheric immersion freezing from ambient aerosol particles. When the aerosol population is physicochemically complex and parameterizations for representative INP types are not yet available, we still struggle to accurately predict INP number concentrations. This project clearly demonstrates that with more laboratory and field measurements that are accompanied by particle composition analysis, the necessary datasets to achieve aerosol–ice formation closure for various locations will emerge, thus providing a robust foundation for guiding the representation of INPs in cloud and climate models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

NASA NTRS2024

Development of a Data Fusion Methodology for Lineload Aerodynamic Databases for a Launch Vehicle during Liftoff and Transition

The need for databases for the distributed loading on launch vehicles during the early portion of flight necessitates the use of expensive computational flows in regimes where wake effects dominate. While also being expensive, this is a regime that computational tools tend to historically have problems simulating accurately. To help tackle this problem, a method of data fusion to combine computational results to wind tunnel derived force and moment data is developed. Using this method, significant reduction in computational costs and increases in confidence of the final product is possible and has been used to generate several databases for the Space Launch System (SLS) at NASA. While the full details of database generation are not part of this work, the crucial method at its core is developed here. Two SLS geometries are used throughout the work to demonstrate the techniques. These are two of the larger geometries and represent both planned crewed missions to the Moon as well as potential cargo missions to deep space. The method uses principal component analysis (PCA) to generate a reduced ordered model (ROM) to help fill in the full parameter space. Other similar techniques are explored, but were not found to have a significant result on the predictions of the ROM. Because the full number of components are kept to generate the model, this lack of difference is expected. This method is then extended to ensure that predicted surfaces match trusted force and moment data derived from wind tunnel testing. This extension is done by setting up a constrained optimization problem in order to minimize the deviation from the surface resolved computational data while still integrating to the desired values. When generating the constrained optimization problem, a weighting factor to balance these competing needs is introduced. The work compares previously introduced weighting terms from similar work to the proposed terms and shows that the previously used terms do not have as desirable behavior in this flow regime. This method is then expanded by developing a technique to incorporate uncertainty quantification into the developed data fusion methodology. This expansion takes a two pronged approach. One examines transferring the uncertainties in the force and moment database and characterizes how those adjustments change the predicted lineloads. The second looks at model form error and looks how rebuilding the model using slightly different data changes the predictions. These two terms are then combined in order to create an uncertainty model that takes both effects into account. The limitations of the proposed methods is then discussed as well as possible techniques to address these shortcomings.

Launch Vehicles↗

NASA NTRS2020

SWIRP (Submm-Wave and Long Wave InfraRed Polarimeter); Development and Characterization of a Sub-Mm Polarimeter for Ice Cloud Investigations

A major source of uncertainty in climate models is the presence, shape and distribution of ice particles in the uppermost layers of the clouds. The effects of this component are poorly constrained, turning ice particles into an almost-free variable in many climate models.NASA-GSFC is developing a new instrument aimed at measuring the size and shape of ice particles. The instrument consists of two sub-mm polarimeters (at 220 and 670 GHz) coupled with a long-wave infrared polarimeter at 10 micron. Each polarimeter has identical V-pol and H-pol channels; the axes of polarization are defined geometrically by the orientation of the waveguide elements, and the purity has been measured in the lab. The instrument is configured as a conical scanner, suitable for deployment as a payload on a small satellite or on a high-altitude sub-orbital platform. From a 400 km orbit, the instrument has a 3dB spatial resolution of 20 (10) km at 220 (670) GHz and a swath of 600 km over 180 degrees of view.The BAPTA (Bearing And Power Transfer Assembly) carries heritage from the SSMIS design, now in its 22nd year of on-orbit operation, but with a much reduced SWaP (Size Weight and Power) footprint, suitable for a small satellite.The main components of the instrument have been fabricated and are undergoing final testing prior to their integration as a single unit. The sub-mm channels have dedicated secondary reflectors which illuminate a shared primary reflector. The receiving units are placed behind the focal point of the optical arrangement, so that all beams equally illuminate the primary reflector and are almost co-located on the ground (within a single 220 GHz footprint). Primary and secondary beam patterns have been measured and verified to match the as-designed expectations. A Zytex (TM) window is deployed to protect the secondary reflectors and the feed horns from debris and other contaminants, and to reduce the heat load from the active (hot) IR calibration unit. The insertion loss of Zytex has been measured and is accounted in the calibration equation of the sub-mm channels.The radiometric performance of the sub-mm receivers has been characterized in the lab and under operational conditions of temperature and pressure.This paper discusses the design constraints on the sub-mm components, details of the scientific goals and their flowdown, and describes the characterization of the polarimeters. Options to optimize the layout and distribution of the masses within the assembly, with the goal of making the instrument even more compact and fully-compatible with cubesat-class satellites will be presented.

De Amici, G.↗