Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model validation”

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 559 records · Page 31

Two-Fluid and Discrete Element Modeling of a Parallel Plate Fluidized Bed Heat Exchanger for Concentrating Solar Power

A novel high-temperature particle solar receiver is developed using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a fluidized bed at the lower section, leaving the system from the bottom at a constant flowrate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built by collaborators at the Colorado School of Mines and their data are used for model validation. In this experimental setup, near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using the discrete element method and a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A new continuum conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2 K.

concentrating solar power↗

Validation of Reanalysis-Based Offshore Wind Resource Characterization Using Lidar Buoy Observations

The offshore wind industry in the U.S. is gaining strong momentum to achieve sustainable energy goals, and the need for observations to provide resource characterization and model validation is greater than ever. Pacific Northwest National Laboratory (PNNL) operates two lidar buoys for the U.S. Department of Energy (DOE) in order to collect hub height wind data and associated meteorological and oceanographic information near the surface in areas of interest for offshore wind development. This work evaluates the performance of commonly used reanalysis products and spatial approximation techniques using lidar buoy observations off the coast of New Jersey and Virginia, USA. Reanalysis products are essential tools in order to set performance expectations and quantify the wind resource variability at a given site. Long-term accurate observations at typical wind turbine hub-heights have been lacking at offshore locations. Using wind speed observations from both lidar buoy deployments, biases and degrees of correspondence for the Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), and the analysis system of the Rapid Refresh (RAP) are examined both at hub height and near surface. Results provide insights on the performance and uncertainty of using reanalysis products for long-term wind resource characterization.

lidar buoy, model validation, offshore wind energy↗

Cookoff experiments of a melt cast explosive (Comp-B3)

Validated models of melt cast explosives exposed to accidental fires are essential for safety analysis. In the current work, we provide several experiments that can be used to develop and validate cookoff models of melt cast explosives such as Comp-B3 composed of 60:40 wt% RDX:TNT. We present several vented and sealed experiments from 2.5 mg to 4.2 kg of Comp-B3 in several configurations. We measured pressure, spatial temperature, and ignition time. Some experiments included borescope images obtained during both vented and sealed decomposition. We observed the TNT melt, the suspension of RDX particles in the melt, bubble formation caused by RDX decomposition, and bubble-induced mixing of the suspension. The RDX suspension did not completely dissolve, even as temperatures approached ignition. Our results contrast with published measurements of RDX solubility in hot TNT that suggest RDX would be completely dissolved at these high temperatures. These different observations are attributed to sample purity. We did not observe significant movement of the two-phase mixture until decomposition gases formed bubbles. Bubble generation was inhibited in our sealed experiments and suppressed mixing.

42 ENGINEERING↗

Simulation studies of the interaction of laser radiation with additively manufactured foams

The interaction of laser radiation with foams of various porosities and low densities has been the subject of several numerical and experimental studies. In all cases, the modeling of low-Z under-dense foams as uniform gases of equivalent average density using standard radiation-hydrodynamics codes has resulted in heat-front velocities that are considerably faster than those observed experimentally. It has been theoretically conjectured that this difference may be attributed to the breakdown of the foam’s morphology, leading to a dynamics of filament expansion where the ion and electron energy partitions are significantly different from those calculated using the uniform gas model. We found that 3D computer simulations employing a disconnected representation of the foam’s microstructure which allowed for the dynamics of foam element heating, expansion, and stagnation largely supported the theoretical picture. Simulations using this model for laser experiments on under-dense 2 mg cc -1 SiO 2 aerogel foams reproduced the experimental data fairly well. We used the validated model in simulations of low-density structured foam-like materials (produced via additive manufacturing) with a variety of morphologies. We found that the log-pile configurations were consistent with the analytical propagation model of Gus’kov et al. Further validation of the model was obtained by simulating experiments performed at the Jupiter Laser Facility using the log-pile and octet-truss foam morphologies. Simulations of the foam-laser interaction using a wave propagation code showed that the microstructure was able to enhance stimulated Brillouin scattering (SBS) by concentrating the light energy into density holes. In turn, this promotes laser filamentation, reducing SBS and bringing the predicted values closer to the experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Phenomenological Model for the Transport and Plugging of Aerosol through Stress Corrosion Cracks

Recent experimental evidence from Sandia (2018) demonstrated partial microchannel plugging due to particulate deposition. A direct consequence is that the accuracy of consequence assessments can be improved by accounting for the leak path deposition of aerosol in the source term. The filtration effect of microchannels, such as cracks, is relevant to the source term assessments. An important side effect of aerosol deposition in leak paths could be the plugging of the leak path. This report summarizes current progress (as of end of FY19) on the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for prediction of aerosol transport, deposition, and plugging in leak paths that are similar to stress corrosion cracks while accounting for potential plugging formation. This effort includes (i) development of a numerical model to analyze the various deposition processes in leak paths and provide quantitative estimates of penetration factors, as well as an understanding of the variables that affect them, (ii) confirmation of model validity with theoretical and experimental data, and (iii) parametric analysis for different scenarios including various particle sizes, pressures, and crack dimensions. Such a model could be used to improve the accuracy of consequence assessments and reduce the uncertainty of radiological consequence predictive analyses by taking the filtering effect of leak path aerosol deposition and plugging into account in the source term. This report includes: (i) an overview of main stress corrosion cracking characteristics, (ii) an overview of the developed numerical model to analyze the various deposition processes in leak paths, to provide quantitative estimates of penetration factors, and to gain an understanding of the variables that affect them, (iii) a summary of validation results and verification of the model’s validity with recent experimental results from Sandia and others, and (iv) a model application example.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Foundational Open Source Solar System Modeling Through Improvement and Validation of the System Advisor Model and PVWatts (FY19-FY21 Final Technical Report)

Accelerating intelligent deployment of solar energy technologies demands accurate system modeling every step of the way. From project development to policy research, grid integration studies to development of novel technologies, industry and researchers alike need validated, transparent, easy-to-use, extensible, cutting-edge, and accurate models of both the performance and financing of solar systems. The System Advisor Model (SAM) and PVWatts tools provide a platform to fill that need. The overarching goal of this set of software tools is to enable accurate PV system modeling across the industry, and our usage metrics indicate that we continue to succeed in that endeavor, with a user starting SAM every 2 minutes globally, and over 17 million PVWatts hits per month. This project leveraged DOEs past investment in the SAM and PVWatts platforms to continue to provide valuable and extensible PV, battery, and financial modeling resources to the larger solar community. We pursued multiple avenues in parallel: software maintenance and technical support that are foundational to the continued usability of the SAM and PVWatts platforms; platform and PV model improvements and stakeholder engagement activities that are core to the continued relevance of the platforms; and open source activities to foster the continued creation of a vibrant open-source community around the SAM and PVWatts tools, which opens up exciting new opportunities for industry interaction.

14 SOLAR ENERGY↗

Experiments and CFD simulation of mass transfer and hydrodynamics in a cylindrical bubble column

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.

42 ENGINEERING↗

Validation in the Absence of Observed Events

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?

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Analysis of soot formation from aviation fuels in laminar counterflow flames

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.

Counterflow↗

Heterogeneity, pore pressure, and injectate chemistry: Control measures for geologic carbon storage

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.

58 GEOSCIENCES↗

Microstructure development during rapid alloy solidification

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

36 MATERIALS SCIENCE↗

Comparison of the Gaussian Wind Farm Model with Historical Data of Three Offshore Wind Farms

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.

17 WIND ENERGY↗

Low-cloud characteristics over the tropical western Pacific from ARM observations and CAM5 simulations

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.

54 ENVIRONMENTAL SCIENCES↗

Simulation and Emulation of X-Ray Diffraction from Dynamic Compression Experiments

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.

36 MATERIALS SCIENCE↗

Turbulent nonpremixed cool flames: Experimental measurements, Direct Numerical Simulation, and manifold-based combustion modeling

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of Am-241 as an alternative fuel source in a radioisotope thermoelectric generator

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.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Wave resource characterization at regional and nearshore scales for the U.S. Alaska coast based on a 32-year high-resolution hindcast

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.

54 ENVIRONMENTAL SCIENCES↗