Development and validation of a slurry model for chemical hydrogen storage in fuel cell vehicle applications
Not Available
SEARCH · Engineering Papers
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.
Not Available
Flow visualization, Particle Image Velocimetry (PIV) and dissolved oxygen charging experiments are carried out in a small-scale cylindrical air-water bubble column at different superficial gas velocities. Void fraction profile, total gas holdup, liquid phase velocity field and volumetric mass transfer coefficient are obtained. Details of each experiment are given, with a new procedure for obtaining a 2D void fraction profile from flow visualization by image processing and time averaging which could provide useful information for CFD model validation. The experimental results are analyzed and compared with simulations where good agreement is found for both local and global flow patterns. The averaged governing equations and constitutive relations used in the simulation are discussed briefly, with emphasis on the proper comparison between simulations and experiments based on the derivation of the two-fluid model. Finally, by conducting three sets of experiments related to phase distribution, velocity field and mass transfer in the same apparatus, a validated CFD model for multiphase mass transfer is established and can be used in the future design of a multi-phase reactor.
Here our paper addresses the problem of validating models in the absence of observed events, in the area of Weapons of Mass Destruction terrorism risk assessment. We address that problem with a broadened definition of “Validation,” based on “backing up” to the reason why modelers and decision makers seek validation, and from that basis re-define validation as testing how well the model can advise decision makers in terrorism risk management decisions. We develop that into two conditions: Validation must be based on cues available in the observable world; and it must focus on what can be done to affect that observable world, i.e. risk management. That in turn leads to two foci: 1.) the risk generating process, 2.) best use of available data. Based on our experience with nine WMD terrorism risk assessment models, we then describe three best use of available data pitfalls: SME confidence bias, lack of SME cross-referencing, and problematic initiation rates. Those two foci and three pitfalls provide a basis from which we define validation in this context in terms of four tests -- Does the model: … capture initiation? … capture the sequence of events by which attack scenarios unfold? … consider unanticipated scenarios? … consider alternative causal chains? Finally, we corroborate our approach against three key validation tests from the DOD literature: Is the model a correct representation of the simuland? To what degree are the model results comparable to the real world? Over what range of inputs are the model results useful?
Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.
Due to increased pollutant loads and water use from coastal development and population growth, occurrences of low-dissolved oxygen and "hypoxic zones" have increased. Reports of fish kills and water quality impairment are also becoming more frequent in many coastal waters. Water quality managers and regulatory agencies rely on numerical modeling tools to quantify the relative contributions of anthropogenic and "natural" pollutant loads (nutrients and biochemical oxygen demand) on dissolved oxygen levels and use the results for remedial activities and source control. The ability to conduct seasonlong simulations with sufficient nearshore resolution is therefore a key requirement. Mesh flexibility and the ability to increase site specific resolution without disturbing the larger domain setup and calibration are critical. The objective of this effort was to develop a robust biogeochemical model suitable for simulation of water quality dynamics including dissolved oxygen in complex coastal environments with multiple tidal channels, tidal flats, and density-driven circulation using unstructured-grid formulation. This paper presents an offline unstructured biogeochemical model that uses the Finite Volume Coastal Ocean Model (FVCOM) discretization of the study domain and the corresponding hydrodynamic solution to drive biogeochemical kinetics based on a water quality model CE-QUAL-ICM. In this paper, the linkage between selected hydrodynamic and water quality models is subjected to several scalar transport and biogeochemical module tests (plume transport and dilution, BOD/DO sag, and phytoplankton/nutrients reaction), and results are compared to their analytical solutions as part of model validation. A preliminary application of the biogeochemical model with a year-long simulation of Hood Canal basin in Puget Sound, USA, is presented as an example and a test of the tool in a real estuary setting. The model reproduced the dynamics and seasonal variations in the biogeochemical state variables and was used to test short-term wind-driven dynamics that could influence dissolved oxygen concentrations in Hood Canal.
Desirable outcomes for geologic carbon storage include maximizing storage efficiency, preserving injectivity, and avoiding unwanted consequences such as caprock or wellbore leakage or induced seismicity during and post injection. Here, to achieve these outcomes, three control measures are evident including pore pressure, injectate chemistry, and knowledge and prudent use of geologic heterogeneity. Field, experimental, and modeling examples are presented that demonstrate controllable GCS via these three measures. Observed changes in reservoir response accompanying CO 2 injection at the Cranfield (Mississippi, USA) site, along with lab testing, show potential for use of injectate chemistry as a means to alter fracture permeability (with concomitant improvements for sweep and storage efficiency). Further control of reservoir sweep attends brine extraction from reservoirs, with benefit for pressure control, mitigation of reservoir and wellbore damage, and water use. State-of-the-art validated models predict the extent of damage and deformation associated with pore pressure hazards in reservoirs, timing and location of networks of fractures, and development of localized leakage pathways. Experimentally validated geomechanics models show where wellbore failure is likely to occur during injection, and efficiency of repair methods. Use of heterogeneity as a control measure includes where best to inject, and where to avoid attempts at storage. Lastly, an example is use of waste zones or leaky seals to both reduce pore pressure hazards and enhance residual CO 2 trapping.
Abstract Solidification processing of structural alloys can take place over an extremely wide range of solid–liquid interface velocities spanning six orders of magnitude, from the low-velocity constitutional supercooling limit of microns/s to the high-velocity absolute stability limit of m/s. In between these two limits, the solid–liquid interface is morphologically unstable and typically forms cellular-dendritic microstructures, but also other microstructures that remain elusive. Rapid developments in additive manufacturing have renewed the interest in modeling the high-velocity range, where approximate analytical theories provide limited predictions. In this article, we discuss recent advances in phase-field modeling of rapid solidification of metallic alloys, including a brief description of state-of-the-art experiments used for model validation. We describe how phase-field models can cope with the dual challenge of carrying out simulations on experimentally relevant length- and time scales and incorporating nonequilibrium effects at the solid–liquid interface that become dominant at rapid rates. We present selected results, illustrating how phase-field simulations have yielded unprecedented insights into high-velocity interface dynamics, shedding new light on both the absolute stability limit and the formation of banded microstructures that are a hallmark of rapid alloy solidification near this limit. We also discuss state-of-the-art experiments used to validate those insights. Graphical abstract
The phase 2 effort consisted of a series of ground and flight test measurements to obtain data for validation of the Statistical Energy Analysis (SEA) model. Included in the gound tests were various transfer function measurements between vibratory and acoustic subsystems, vibration and acoustic decay rate measurements, and coherent source measurements. The bulk of these, the vibration transfer functions, were used for SEA model validation, while the others provided information for characterization of damping and reverberation time of the subsystems. The flight test program included measurements of cabin and cockpit sound pressure level, frame and panel vibration level, and vibration levels at the main transmission attachment locations. Comparisons between measured and predicted subsystem excitation levels from both ground and flight testing were evaluated. The ground test data show good correlation with predictions of vibration levels throughout the cabin overhead for all excitations. The flight test results also indicate excellent correlation of inflight sound pressure measurements to sound pressure levels predicted by the SEA model, where the average aircraft speech interference level is predicted within 0.2 dB.
A recent expert elicitation showed that model validation remains one of the largest barriers for commercial wind farm control deployment. The Gaussian-shaped wake deficit model has grown in popularity in wind farm field experiments, yet its validation for larger farms and throughout annual operation remains limited. This article addresses this scientific gap, providing a model comparison of the Gaussian wind farm model with historical data of three offshore wind farms. The energy ratio is used to quantify the model’s accuracy. We assume a fixed turbulence intensity of $I_∞$ = 6% and a standard deviation on the inflow wind direction of $σ_{wd}$ = 3° in our Gaussian model. First, we demonstrate the non-uniqueness issue of $I_∞$ and $σ_{wd}$, which display a waterbed effect when considering the energy ratios. Second, we show excellent agreement between the Gaussian model and historical data for most wind directions in the Offshore Windpark Egmond aan Zee (OWEZ) and Westermost Rough wind farms (36 and 35 wind turbines, respectively) and wind turbines on the outer edges of the Anholt wind farm (110 turbines). Turbines centrally positioned in the Anholt wind farm show larger model discrepancies, likely due to deep-array effects that are not captured in the model. A second source of discrepancy is hypothesized to be inflow heterogeneity. In future work, the Gaussian wind farm model will be adapted to address those weaknesses.
Here, this study evaluates the ability of the Community Atmospheric Model version 5 (CAM5) to reproduce low clouds observed by the Atmospheric Radiation Measurement (ARM) cloud radar at Manus Island of the tropical western Pacific during the Years of Tropical Convection. Here low clouds are defined as clouds with their tops below the freezing level and bases within the boundary layer. Low-cloud statistics in CAM5 simulations and ARM observations are compared in terms of their general occurrence, mean vertical profiles, fraction of precipitating versus nonprecipitating events, diurnal cycle, and monthly time series. Other types of clouds are included to put the comparison in a broader context. The comparison shows that the model overproduces total clouds and their precipitation fraction but underestimates low clouds in general. The model, however, produces excessive low clouds in a thin layer between 954 and 930 hPa, which coincides with excessive humidity near the top of the mixed layer. This suggests that the erroneously excessive low clouds stem from parameterization of both cloud and turbulence mixing. The model also fails to produce the observed diurnal cycle in low clouds, not exclusively due to the model coarse grid spacing that does not resolve Manus Island. Lastly, this study demonstrates the utility of ARM long-term cloud observations in the tropical western Pacific in verifying low clouds simulated by global climate models, illustrates issues of using ARM observations in model validation, and provides an example of severe model biases in producing observed low clouds in the tropical western Pacific.
Many important aspects of the dynamic thermo-mechanical response of materials occur at the mesoscale, i.e. a physical scale of interactions smaller than what can be adequately described by homogenous behaviors, yet larger than the scale of the atomic lattice. Concurrent advancements in computational power, continuum theory, and experimental diagnostics are enabling unprecedented understanding of such interactions. However, we cannot develop a sufficient level of confidence in such mesoscale capability until the constitutive description of the underlying constituents is reliably representative of their actual physical behavior. Therefore, there is a strong need to combine experimental, modeling, and data-science techniques to validate models of the thermomechanical response of individual single crystals. One experimental diagnostic with high potential impact to shock physics and materials science is in-situ x-ray diffraction. This paper is primarily focused on simulation of x-ray diffraction in shock physics, but with an aim toward quantifying parametric uncertainty of simulation models. Here, we develop and demonstrate a data-science and model-driven approach to constrain the parameterization of continuum models of crystal lattice deformation associated with the shock response of crystalline materials. The framework is built around the connection between continuum hydrodynamic simulations of lattice deformation and a new Bragg diffraction simulation code, BarberShop. The dynamic deformation of a crystal lattice is modeled using the DiscoFlux model within an arbitrary Lagrangian-Eulerian hydrodynamic code, FLAG. These detailed continuum simulations of lattice deformation can be computationally slow, thus a statistical model is used to emulate the evolution of lattice deformation fields in time and across the considered model parameter space. Emulated lattice deformation fields can then be generated rapidly for any combination of physics model parameters. In turn, these fields can be fed into BarberShop to realize a rapid prediction of Bragg diffraction patterns associated with particular values of physics model parameters. The framework enables parameterization of the single crystal model to obtain Bragg diffraction patterns that most closely resemble a corresponding measurement. Furthermore, the framework naturally provides sensitivities of the lattice deformation to the physics parameters. We highlight the utility of this framework through the application to a synthetic closed-loop inverse problem leading to the parameterization of a single crystal material model. As a model problem, we consider the dynamic response of the energetic molecular crystal, cyclotrimethylenetrinitramine (or RDX), under dynamic compression induced by simulated flyer plate impact experiments.
Turbulence, low-temperature chemistry, and their interactions in the form of turbulent cool flames are critical to understanding and improving advanced engines. Design of such engines requires tractable simulations which in turn necessitate turbulent combustion models that can account for cool flames. While manifold-based turbulent combustion models are an attractive option for hot flames, their applicability to cool flames is not yet fully understood. This is partially due to the lack of turbulent cool flame experiments, which has made model validation difficult. This work addresses these points with a combined experimental and computational investigation of turbulent nonpremixed dimethyl ether cool flames. First, a Co-flow Axisymmetric Reactor-Assisted Turbulent (CARAT) burner is developed and tested. Turbulent cool flames are studied using the formaldehyde planar laser-induced fluorescence (PLIF), acetone PLIF, and planar Rayleigh scattering techniques. The acetone PLIF signals are converted into mixture fraction values, and measurements of time-averaged temperature are derived from the Rayleigh scattering signals by taking advantage of the similarities between cool flames and unburned mixtures. Second, a Direct Numerical Simulation (DNS) of the CARAT burner is conducted. The flame is shown to be sensitive to thermal boundary conditions such that the stabilization method is dependent on a 10 K change in inlet temperature. Comparisons of first- and second-order statistics of temperature, mixture fraction, and formaldehyde between the DNS and experiment show good agreement. The validity of manifold-based turbulent combustion models in turbulent cool flames is then explored by first studying the effective Lewis numbers in the flame using a local differential diffusion parameter analysis. Turbulent cool flames are shown to have effective unity Lewis number transport even at a much lower Reynolds number than typically required for hot flames. Good agreement is shown for temperature and formaldehyde between DNS conditional means and one-dimensional nonpremixed flame solutions. Additionally, the motion of the cool flame in mixture fraction space is shown to be well described by one-dimensional nonpremixed flame solutions. These points indicate that the fundamental physics of turbulent cool flames are analogous to turbulent hot flames, implying that the same reduced-order manifold-based modeling approaches utilized for turbulent hot flames can be utilized for turbulent cool flames.
Radioisotope thermoelectric generators (RTGs) have been utilized in the USA to power satellites and space exploration equipment for over half a century. In our preceding paper, Dustin and Borrelli (2021), we assessed the performance of 241 Am as an alternative to 238 Pu by modeling the specific decay-heat and gamma spectra of several potential radioisotopes in SCALE-ORIGEN, and comparing these with 238 Pu An analysis of heat output, heat stability, safety, availability and required mass showed that of the isotopes reviewed, 241 Am had the greatest potential for replacing 238 Pu, with additional design considerations. In our current work, we have assessed radiation transport for two configurations of an 241 Am powered Multi-Mission Radioisotope Thermoelectric Generator (MMRTG), due to the potential health effects for operators, fabricators, and others with the potential to be exposed to neutrons and lower energy gamma rays. We have modeled the use of americium as a fuel source replacement for plutonium in current MMRTG designs. The assessment has been performed using MCNP6, in terms of neutron and gamma emissions that might interfere with sensitive equipment and power system production. The radiation transport behavior in terms of dose rate of the 241 Am fueled models is compared with that of a 238 Pu fueled model. Finally, additional work was performed to validate modeling simplifications with regards to thermocouple material configurations.
A wave resource characterization was performed for the southern coast of Alaska based on a 32-year hindcast covering the period from 1979 to 2010. The characterization closely followed International Electrotechnical Commission (IEC) Technical Specifications. An unstructured-grid Simulating WAves Nearshore (SWAN) model, which had an approximate spatial resolution of 300 m within 30 km from the nearest shoreline, was developed. Additionally, extensive model validation and error characterization was performed. The model was found to perform well with an average absolute percent error of 8.6% in significant wave height, averaged over 18 buoys. Statistics for the six IEC wave resource parameters were calculated and aggregated at 20 km from shore to quantify the incident wave power and its variability at a regional scale. The southern coast of the Aleutian Archipelago was found to have the most available wave energy in the region. A nearshore resource assessment was performed by evaluating resource hotspots located 1 km from shore. More than 900 nearshore stations had an Optimum Hotspot Identifier value of 5 (kW/m) at diverse water depths, thereby positioning Alaska as a promising location for wave energy development.
Recent developments in machine learning interatomic potentials (MLIPs) have empowered even nonexperts in machine learning to train MLIPs for accelerating materials simulations. However, reproducibility and independent evaluation of presented MLIP results is hindered by a lack of clear standards in current literature. In this Perspective, we aim to provide guidance on best practices for documenting MLIP use while walking the reader through the development and deployment of MLIPs including hardware and software requirements, generating training data, training models, validating predictions, and MLIP inference. We also suggest useful plotting practices and analyses to validate and boost confidence in the deployed models. Finally, we provide a step-by-step checklist for practitioners to use directly before publication to standardize the information to be reported. Altogether, we hope that our work will encourage the reliable and reproducible use of these MLIPs, which will accelerate their ability to make a positive impact in various disciplines including materials science, chemistry, and biology, among others.
Wind turbines are typically deployed to provide energy, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL), it is demonstrated that advanced wind turbine controls can be employed to support higher contributions of wind, and to demonstrate ways that wind can play a role in supporting grid stability in islanded or grid-connected configurations. This paper documents models of various subsystem comprising a portion of NREL's Flatirons campus that will be used in three subsequent reports to demonstrate capabilities of advanced wind turbine controls. The series of reports will detail advanced capabilities of distributed wind turbines to provide support to isolated grids, distribution grids, and microgrids. We developed models to simulate a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a diesel generator (2 MW) and various types of loads (critical, dynamic). The model of the subsystems in MATLAB/Simulink are validated with available data from real-world components on NREL's Flatirons Campus. These validated models can be configured for various studies including four MIRACL use cases: 1) isolated grids, 2) microgrids, and 3) behind-the-meter, and 4) front-of-the-meter wind turbine deployments.
In an isolated grid, wind turbines are typically deployed to provide energy to maximize energy production, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL) this paper demonstrates, through desktop simulations, advanced wind turbine controls that can be employed to support higher contributions of wind in isolated grids, and to demonstrate ways that wind can play a role in supporting stability of an isolated grid. This isolated grid used in these desktop simulations is comprised of a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a simulated diesel generator (2 MW) and various types of loads (critical, dynamic). We developed a model of the subsystems in MATLAB/Simulink and validated them with available data from real-world components on NREL's Flatirons Campus. These validated models are then configured for various case studies. We compare the output of the desktop simulation with a baseline case with the diesel generator. Active and reactive power control of the wind turbine can help improve frequency and voltage responses in the isolated grid, respectively. By utilizing a small integrated battery energy storage system in the DC-link of the wind turbine, we also demonstrate that wind turbines can help blackstart a critical load comparable to its rated power and support other renewables (e.g. solar PV) come online and pick up an additional load. This report illustrates some of these reliability and resilience services a wind turbine can provide in an isolated grid.
The Stable Salt Reactor (SSR) combines the proven structural design of light water reactor fuel assemblies with the inherent safety and fuel-cycle advantages of molten salt technology. In its fast reactor configuration, the SSR utilizes recycled nuclear waste as fuel, sealed within narrow salt-filled fuel pins and cooled by a surrounding liquid salt coolant. Reliable transfer of heat from the molten fuel salt through the cladding to the external coolant is essential for both reactor safety and performance. This work investigates conjugate heat transfer (CHT) in the SSR’s salt-filled fuel pins using NekRS, a high-fidelity spectral element computational fluid dynamics (CFD) solver. The analyses capture internal natural convection within the molten fuel salt and external forced convection in the coolant, under steady-state and transient operating conditions. Parametric studies evaluate how variations in reactor power and coolant flow rate influence heat transfer distributions and system response. The high-fidelity CFD results are time-averaged and post-processed for direct comparison with moderate-fidelity Reynolds-averaged Navier–Stokes (RANS) models, and for the development of reduced-order models within the SAM system code. These validated models support fast-running safety analyses of normal and off-normal transients, improving predictive capability for key safety margins. By integrating advanced CFD with system-level safety tools, this study strengthens the modeling framework for SSR design, reduces uncertainty in molten salt CHT simulations, and accelerates the engineering and licensing of next-generation nuclear reactors.