Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model benchmark”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,549 records · Page 86

Validation Assessment for the Soil Moisture Active Passive (SMAP) Level 4 Carbon (L4_C) Data Product Version 5

The post-launch Cal/Val phase of the SMAP mission is guided by two primary objectives for each science product team: 1) to calibrate, verify, and improve the performance of the science algorithms, and 2) validate accuracies of the science data products as specified in the SMAP Level-1 mission science requirements. Algorithm science and product maintenance activities during the SMAP extended mission phase have also involved periodic algorithm calibration and product refinements to maintain or enhance product consistency and performance as well as science utility. This report provides an assessment of the latest (Version 5) SMAP Level 4 Carbon (L4_C) product. The L4_C Version 5 (v5) global record now spans more than six years (March 2015 – present) of SMAP operations and has benefited from five major reprocessing updates to the operational product. These reprocessing events and L4_C product release updates have incorporated various algorithm refinements and calibration adjustments to account for similar refinements to the upstream GEOS land model assimilation system, SMAP brightness temperatures, and MODIS vegetation inputs used for L4_C processing. The SMAP L4_C algorithms utilize a terrestrial carbon flux model informed by daily surface and root zone soil moisture information contributed from the SMAP Level 4 Soil Moisture (L4_SM) product along with optical remote sensing-based (e.g. MODIS-based) land cover and canopy fractional photosynthetic active radiation (fPAR), and other ancillary biophysical data. The carbon flux model estimates global daily net ecosystem CO2 exchange (NEE) and the component carbon fluxes, namely, vegetation gross primary production (GPP) and soil heterotrophic respiration (Rh). Other L4_C product elements include surface (~0-5 cm depth) soil organic carbon (SOC) stocks and associated environmental constraints to these processes, including soil moisture-related controls on GPP and ecosystem respiration (Kimball et al. 2014, Jones et al. 2017). The L4_C product addresses SMAP carbon cycle science objectives by: 1) providing a direct link between terrestrial carbon fluxes and underlying freeze/thaw and soil moisture-related constraints to these processes, 2) documenting primary connections between terrestrial water, energy and carbon cycles, and 3) improving understanding of terrestrial carbon sink activity. The SMAP L4_C algorithms and operational product are mature and at a CEOS Validation Stage 4 level (Jackson et al. 2012) based on extensive validation of the multi-year record against a diverse array of independent benchmarks, well characterized global performance, and systematic refinements gained from five major reprocessing events. There are no Level-1 mission science requirements for the L4_C product; however, self-imposed requirements have been established focusing on NEE as the primary product field for validation, and on demonstrating L4_C accuracy and success in meeting product science requirements (Jackson et al. 2012). The other L4_C product fields also have strong utility for carbon science applications (e.g., Liu et al. 2019, Endsley et al. 2020); however, analysis of these other fields is considered secondary relative to primary validation activities focusing on NEE. The L4_C targeted accuracy requirements are to meet or exceed a mean unbiased root-mean-square error (ubRMSE, or standard deviation of the error) for NEE of 1.6 g C m-2 d-1 and 30 g C m-2 yr-1, emphasizing northern (≥45°N) boreal and arctic ecosystems; this accuracy is similar to that of tower eddy covariance measurement-based observations (Baldocchi 2008). Methods used for the latest v5 L4_C product performance and validation assessment have been established from the SMAP Cal/Val plan and previous studies (Jackson et al. 2012, Jones et al. 2017) and include: 1) consistency evaluations of the product fields against earlier product releases (version 4 or earlier); 2) comparisons of daily carbon flux estimates with independent tower eddy 4covariance measurement-based daily carbon (CO2) flux observations from core tower validation sites (CVS); and 3) consistency checks against other global carbon products, including soil carbon inventory records, global GPP records derived from tower observation upscaling methods, and satellite-based observations of canopy solar induced chlorophyll fluorescence (SIF) as a surrogate for GPP. Metrics used to evaluate relative agreement between L4_C product fields and observational benchmarks include correlation (r-value), RMSE differences, bias and model sensitivity diagnostics. Following these validation criteria, the present report provides a validation assessment of the latest L4_C product release (v5). Detailed descriptions of the L4_C algorithm and additional global product accuracy and performance results are given elsewhere (Jones et al. 2017, Endsley et al. 2020). The v5 L4_C product replaces earlier product versions and continues to show: (i) accuracy and performance levels meeting or exceeding SMAP L4_C science requirements; (ii) improvement over the previous product version (version 4); and (iii) suitability for a diversity of science applications. Example L4_C applications from the recent literature include clarifying environmental trends and controls on the northern terrestrial carbon sink (Liu et al. 2019), diagnosing drought-related impacts on ecosystem productivity (Li et al. 2020), and regional monitoring of cropland conditions for projecting annual yields (Wurster et al. 2020). 2EXPECTED L4_C ALGORITHM AND PRODUCTPERFORMANCE The L4_C algorithm performance, including variance and uncertainty estimates of model outputs, was determined during the mission pre-launch phase through spatially explicit model sensitivity studies using available model inputs similar to those currently being used for operational production and evaluating the resulting model simulations over the observed range of northern (≥45 °N) and global conditions (Kimball et al. 2012, Entekhabi et al. 2014). The L4_C algorithm options were also evaluated during the mission prelaunch phase, including deriving canopy fPAR from lower order NDVI (Normalized Difference Vegetation Index) inputs in lieu of using MODIS (MOD15) fPAR; and including an explicit model representation of boreal fire disturbance recovery impacts. These results indicated that the L4_C accuracy requirements (i.e., NEE ubRMSE ≤ 30 g C m-2 yr-1or ≤ 1.6 g C m-2 d-1) could be met from the baseline algorithms over more than 82% and 89% of global and northern vegetated land areas, respectively (Yi et al. 2013, Kimball et al. 2014). The global L4_C algorithm error budget for NEE derived during the mission prelaunch phase indicated that the estimated NEE ubRMSE uncertainty is proportional to GPP and is therefore larger in higher biomass productivity areas, including forests and croplands (Kimball et al. 2014). Likewise, NEE ubRMSE uncertainty is expected to be lower in less-productive areas, including grasslands and shrublands. Expected model NEE ubRMSE levels were also generally within targeted accuracy levels for characteristically less-productive boreal and Arctic biomes, even though relative model error as a proportion of total productivity (NEE RMSE / GPP) may be large in these areas. The estimated NEE uncertainty was lower than expected in some warmer tropical high biomass productivity areas (e.g. Amazon rainforest) because of reduced low temperature and moisture constraints to the L4_C respiration calculations so that the bulk of model uncertainty is contributed by GPP in these areas. Model NEE uncertainty in the African Congo was estimated to

SMAP↗

A numerical evaluation of the ambient air temperature in the Electron-Ion Collider tunnel

The Electron-Ion Collider (EIC) is a next-generation collider-accelerator that may require consistent operating temperature conditions for the beams within the accelerator tunnels to maintain stable operation. Variations in ambient temperature within the tunnel can cause thermal expansion of beampipe and component supports and can negatively affect the tunnel equipment, impacting the stability of the beamline. Modifications will be made to the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory (BNL) to create the EIC, which necessitates a temperature model that addresses these modifications. To approach this problem, the consistency of temperature changes in different tunnel sections was first evaluated by plotting RHIC tunnel temperature data at various times of the day and year. From this data, a tunnel section was selected and a 2D temperature model was created for RHIC, EIC, and EIC with added cooling configurations. Soil temperature data was analyzed to determine the maximum, average, and mode soil temperatures, which were used as boundary conditions in different temperature scenarios. Computational fluid dynamics modeling was used to create 2D temperature profiles for the configurations. From this model, the predicted temperatures indicate that further analysis is required to validate the boundary conditions and benchmark the current conditions to allow the prediction of the tunnel ambient conditions at EIC. This research can be used as a preliminary model to create an EIC tunnel cooling system that will increase the operational stability of the EIC. As a result of my work this summer, I have become familiar with computational fluid dynamics, including creating fluid dynamic simulations using ANSYS Fluent and related software. I have also learned about the project process required for planning large-scale engineering projects.

43 PARTICLE ACCELERATORS↗

Inferring Plant Acclimation and Improving Model Generalizability With Differentiable Physics‐Informed Machine Learning of Photosynthesis

Net photosynthesis (A N ) is a key component of the global carbon cycle influencing climate feedback over decadal scales. Although plant acclimation to environmental changes can modify A N , traditional vegetation models in Earth system models (ESMs) often rely on plant functional type (PFT)-specific parameterizations or simplified acclimation assumptions limiting generalizability across time, space, and PFTs. In this study, we developed a differentiable photosynthesis model to learn the environmental dependencies of V c,max25 (maximum carboxylation rate at 25°C, representing photosynthetic capacity), as this genre of hybrid physics-informed machine learning can seamlessly train neural networks and process-based equations together. Compared to PFT-specific parameterization of V c,max25 , learning the environment dependencies of key photosynthetic parameters improved model spatiotemporal generalizability. Applying environmental acclimation to V c,max25 led to substantial variations in global mean A N indicating the need to address acclimation in ESMs. The model effectively captured multivariate observations (V c,max25 , A N , and stomatal conductance (g s )) simultaneously with multivariate constraints, improving generalization across space and PFTs. It also learned sensible acclimation relationships of V c,max25 to different environmental conditions. The model explained more than 54%, 57%, and 62% of the variance of A N , g s , and V c,max25 , respectively, presenting a first global-scale spatial test benchmark of A N and g s . These results highlight the potential for differentiable modeling to enhance process-based modules in ESMs and effectively leverage information from large, multivariate data sets.

54 ENVIRONMENTAL SCIENCES↗

Did You Win the GPU Cloud Lottery? Benchmarking from TFLOPS to Tokens/$

Cloud GPUs are commonly assumed to deliver consistent performance for a given GPU model. This assumption does not always hold: cloud providers employ diverse system configurations and virtualization mechanisms, and GPUs themselves exhibit non-negligible manufacturing variability (the silicon lottery). In this work, we present a large-scale measurement study of GPU performance variability across 11 cloud providers, covering over 3,500 physical GPUs and 6,800 benchmark runs. Our hierarchical analysis shows that while execution-level variation stays below 9%, performance varies by up to 38% across devices and providers for the same GPU model. Regression analysis indicates that driver- and OS-related software factors contribute less than 1% of the variance; instead, silicon lottery effects dominate observed performance variation, and cloud providers further amplify them through persistent, systematic second-order effects.

Slynko, Platon [Silicon Data, New York, USA] (ORCI↗

Speedup of UEDGE Parameter Scans Using Machine-Learning Optimized OpenMP Parallelization and a Continuation Solver

This article presents the OpenMP parallelization of the preconditioning Jacobian assembly and right‐hand side residual evaluation in UEDGE. A continuation algorithm, utilizing the internal NKSOL implicit Jacobian‐Free Newton‐Krylov solver to efficiently scan physical parameters, is also presented. The implemented parallelization reduces the computational time for a benchmark scan run on 32 threads by compared to the serial version when using trained random forest regression models to identify the optimal decomposition of the system of equations. Random forest regression models applied to the UEDGE time‐dependent and continuation solver algorithms did not yield meaningful improvement in computational performance. A benchmark DIII‐D gas injection rate scan in the 0.35–0.75 kA interval, performed on a test cluster using the parallelized code and continuation solver, produced 1066 steady‐state solutions with a 22 s average wall‐clock computational time per steady‐state solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Measurements and modeling of type-I and type-II ELMs heat flux to the DIII-D divertor

Type-I and type-II edge-localized-modes (ELMs) heat flux profiles measured at the DIII-D divertor feature a peak in the vicinity of the strike-point and a plateau in the scrape-off-layer (SOL), which extends to the first wall. The plateau is present in attached and detached divertors and it is found to originate with plasma bursts upstream in the SOL. The integrated ELM heat flux is distributed at ~65% in the peak and ~35% in this plateau. The parallel loss model, currently used at ITER to predict power loads to the walls, is benchmarked using these results in the primary and secondary divertors with unprecedented constraints using experimental input data for ELM size, radial velocity, energy, electron temperature and density, heat flux footprints and number of filaments. The model can reproduce the experimental near-SOL peak within ~20%, but cannot match the SOL plateau. Employing a two-component approach for the ELM radial velocity, as guided by intermittent data, the full radial heat flux profile can be well matched. The ELM-averaged radial velocity at the separatrix, which explains profile widening, increases from ~0.2 km s –1 in attached to ~0.8 km s –1 in detached scenarios, as the ELM filaments' path becomes electrically disconnected from the sheath at the target. The results presented here indicate filaments fragmentation as a possible mechanism for ELM transport to the far-SOL and provide evidence on the beneficial role of detachment to mitigate ELM flux in the divertor far-SOL. However, these findings imply that wall regions far from the strike points in future machines should be designed to withstand significant heat flux, even for small-ELM regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Continued Verification of MOOSE Structural Mechanics Tools for Modeling Core Bowing Phenomena in Fast Reactors

Under the U.S. Department of Energy Office of Nuclear Energy’s Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism whereby increased power (which leads to temperature and flux gradients) influences the core to bow into less reactive configurations when the restraint system is properly designed. The phenomenon includes a complex interplay of radiation transport, duct temperature calculations involving fluid flow and heat transfer, and thermo-mechanical responses to the induced temperature and flux gradients. Structural material properties are also important to determining inelastic response to longer term flux gradients which cause irradiation creep and swelling. While core bowing provides a strong negative reactivity feedback when the restraint system is designed properly, it also results in additional forces between assemblies which increase the loads required to extricate them during refueling or control rod movement. Therefore, the restraint system must be designed with these tradeoffs in mind. The first stage of the work, which commenced in FY21 and continues through FY22, assesses thermo-mechanical modeling tools for producing core bowing predictions consistent with conventional tools. The Multiphysics Object Oriented Simulation Environment (MOOSE) Tensor Mechanics and Contact Modules are employed. This status report describes work on additional thermo-mechanical benchmark verification problems with increased complexity from the examples demonstrated in FY21. Several benchmark verification examples were selected from the IAEA verification and validation report. These examples involve clusters of ducts representative of a sector of a hexagonal reactor core which bow into each other and cause contact and load pad elevations, as well as single ducts subjected to irradiation fields undergoing swelling and subsequent bowing. The MOOSE-based results were compared to both IAEA benchmark participants’ results, analytic equations as available, and NUBOW-3D, a beam model code developed by Argonne National Laboratory. In every case, the MOOSE results agreed with other simulations results, providing additional verification basis of the tools for this particular physics application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Autonomous Synthesis of Metastable Materials Using a Modular Mixed-Flow Reactor

Understanding and controlling atomic-level processes at solid-liquid interfaces is key to advancing technologies in energy storage, carbon capture, critical element recovery, and materials synthesis. Many of these processes are dominated by the formation of short-lived intermediate precipitates that determine the final properties of synthesized materials. However, studying these intermediates is challenging due to their sensitivity and the reliance on trial-and-error methods. To address this, we developed an automated variable-volume mixed-flow reactor (MFR) to optimize metastable material synthesis and investigate rapid kinetic processes. This state-of-the-art MFR system, paired with an automated modeling framework, enables efficient synthesis and real-time analysis of transient phases. Benchmarking with advanced capabilities, such as wide-/small-angle X-ray scattering, allows us to resolve fast nucleation and growth dynamics that were previously inaccessible. By combining automation, ML-guided optimization, and tailored kinetic modeling, this approach provides a robust platform for improving material design and achieving precise control over solid-liquid reactions.

36 MATERIALS SCIENCE↗

When do short-range atomistic machine-learning models fall short?

In this work, we explore the role of long-range interactions in atomistic machine-learning models by analyzing the effects on fitting accuracy, isolated cluster properties, and bulk thermodynamic properties. Such models have become increasingly popular in molecular simulations given their ability to learn highly complex and multi-dimensional interactions within a local environment; however, many of them fundamentally lack a description of explicit long-range interactions. In order to provide a well-defined benchmark system with precisely known pairwise interactions, we chose as the reference model a flexible version of the Extended Simple Point Charge (SPC/E) water model. Our analysis shows that while local representations are sufficient for predictions of the condensed liquid phase, the short-range nature of machine-learning models falls short in representing cluster and vapor phase properties. These findings provide an improved understanding of the role of long-range interactions in machine learning models and the regimes where they are necessary.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

IPRT Polarized Radiative Transfer Model Intercomparison Project-Phase A

The polarization state of electromagnetic radiation scattered by atmospheric particles such as aerosols, cloud droplets, or ice crystals contains much more information about the optical and microphysical properties than the total intensity alone. For this reason an increasing number of polarimetric observations are performed from space, from the ground and from aircraft. Polarized radiative transfer models are required to interpret and analyse these measurements and to develop retrieval algorithms exploiting polarimetric observations. In the last years a large number of new codes have been developed, mostly for specific applications. Benchmark results are available for specific cases, but not for more sophisticated scenarios including polarized surface reflection and multi-layer atmospheres. The International Polarized Radiative Transfer (IPRT) working group of the International Radiation Commission (IRC) has initiated a model intercomparison project in order to fill this gap. This paper presents the results of the first phase A of the IPRT project which includes ten test cases, from simple setups with only one layer and Rayleigh scattering to rather sophisticated setups with a cloud embedded in a standard atmosphere above an ocean surface. All scenarios in the first phase A of the intercomparison project are for a one-dimensional plane-parallel model geometry. The commonly established benchmark results are available at the IPRT website

radiative transfer↗

Toward a Climate OSSE Framework for Satellite Mission Design

The rich history of observing system simulation experiments (OSSEs) does not yet include a well-established framework for using climate models. The need for a climate OSSE is triggered by the need to quantify the value of a particular measurement for reducing the uncertainty in climate predictions, which differ from numerical weather predictions in that they depend on future atmospheric composition rather than the current state of the weather. However, both weather and climate modeling communities share a need for motivating major observing system investments. Here, we outline a new framework for climate OSSEs that leverages the use of machine learning to calibrate climate model physics against existing satellite data. We demonstrate its application using NASA’s GISS-E3 model to objectively quantify the value of potential future improvements in spaceborne measurements of Earth’s planetary boundary layer. A mature climate OSSE framework should be able to quantitatively compare the ability of proposed observing system architectures to answer a climate-related question, thus offering added value throughout the mission design process, which is subject to increasingly rapid advances in instrument and satellite technology. Technical considerations include selection of observational benchmarks and climate projection metrics, approaches to pinpoint the sources of model physics uncertainty that dominate uncertainty in projections, and the use of instrument simulators. Community and policy-making considerations include the potential to interface with an established culture of model intercomparison projects and a growing need to economically assess the value-driven efficiency of social spending on Earth observations.

54 ENVIRONMENTAL SCIENCES↗

Radiative Forcing by Well-Mixed Greenhouse Gases: Estimates from Climate Models in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4)

The radiative effects from increased concentrations of well-mixed greenhouse gases (WMGHGs) represent the most significant and best understood anthropogenic forcing of the climate system. The most comprehensive tools for simulating past and future climates influenced by WMGHGs are fully coupled atmosphere-ocean general circulation models (AOGCMs). Because of the importance of WMGHGs as forcing agents it is essential that AOGCMs compute the radiative forcing by these gases as accurately as possible. We present the results of a radiative transfer model intercomparison between the forcings computed by the radiative parameterizations of AOGCMs and by benchmark line-by-line (LBL) codes. The comparison is focused on forcing by CO2, CH4, N2O, CFC-11, CFC-12, and the increased H2O expected in warmer climates. The models included in the intercomparison include several LBL codes and most of the global models submitted to the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4). In general, the LBL models are in excellent agreement with each other. However, in many cases, there are substantial discrepancies among the AOGCMs and between the AOGCMs and LBL codes. In some cases this is because the AOGCMs neglect particular absorbers, in particular the near-infrared effects of CH4 and N2O, while in others it is due to the methods for modeling the radiative processes. The biases in the AOGCM forcings are generally largest at the surface level. We quantify these differences and discuss the implications for interpreting variations in forcing and response across the multimodel ensemble of AOGCM simulations assembled for the IPCC AR4.

Collins, W. D.↗

Land carbon models underestimate the severity and duration of drought’s impact on plant productivity

The ability to accurately predict ecosystem drought response and recovery is necessary to produce reliable forecasts of land carbon uptake and future climate. Using a suite of models from the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP), we assessed modeled net primary productivity (NPP) response to, and recovery from, drought events against a benchmark derived from tree ring observations between 1948 and 2008 across forested regions of the US and Europe. We find short lag times (0–6 months) between climate anomalies and modeled NPP response. Although models accurately simulate the direction of drought legacy effects (i.e. NPP decreases), projected effects are approximately four times shorter and four times weaker than observations suggest. This discrepancy between observed and simulated vegetation recovery from drought reveals a potential critical model deficiency. Since productivity is a crucial component of the land carbon balance, models that underestimate drought recovery time could overestimate predictions of future land carbon sink strength and, consequently, underestimate forecasts of atmospheric CO2.

carbon cycle↗

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient thin film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better turbulence models. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity and also the mean temperature and normalized root mean square (rms) temperatures in a single injector cooling flow arrangement. In addition to flowfield measurements, thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures and Mach numbers, across a 30.48 cm long plate equipped with a single injector cooling hole. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 15 axial stations were collected in the same flows. The velocity and temperature data were then compared against Wind-US CFD code predictions for the same flow conditions. The results of this and planned follow-on studies will support NASA's development and assessment of turbulence models for heated flows.

Raman Scattering↗

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient thin film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better turbulence models. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity and also the mean temperature and normalized root mean square (rms) temperatures in a single injector cooling flow arrangement. In addition to flowfield measurements, thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures and Mach numbers, across a 30.48cm long plate equipped with a single injector cooling hole. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 15 axial stations were collected in the same flows. The velocity and temperature data were then compared against Wind-US CFD code predictions for the same flow conditions. The results of this and planned follow-on studies will support NASA's development and assessment of turbulence models for heated flows.

Wernet, Mark P.↗

PIV and Rotational Raman-Based Temperature Measurements for CFD Validation in a Single Injector Cooling Flow

Film cooling is used in a wide variety of engineering applications for protection of surfaces from hot or combusting gases. The design of more efficient thin film cooling geometries/configurations could be facilitated by an ability to accurately model and predict the effectiveness of current designs using computational fluid dynamics (CFD) code predictions. Hence, a benchmark set of flow field property data were obtained for use in assessing current CFD capabilities and for development of better turbulence models. Both Particle Image Velocimetry (PIV) and spontaneous rotational Raman scattering (SRS) spectroscopy were used to acquire high quality, spatially-resolved measurements of the mean velocity, turbulence intensity and also the mean temperature and normalized root mean square (rms) temperatures in a single injector cooling flow arrangement. In addition to flowfield measurements, thermocouple measurements on the plate surface enabled estimates of the film effectiveness. Raman spectra in air were obtained across a matrix of radial and axial locations downstream from a 68.07 mm square nozzle blowing heated air over a range of temperatures and Mach numbers, across a 30.48cm long plate equipped with a single injector cooling hole. In addition, both centerline streamwise 2-component PIV and cross-stream 3-component Stereo PIV data at 15 axial stations were collected in the same flows. The velocity and temperature data were then compared against Wind-US CFD code predictions for the same flow conditions. The results of this and planned follow-on studies will support NASA's development and assessment of turbulence models for heated flows.

Turbulence Modeling↗

Atmospheric Boundary-Layer and Flutter Computations Using CFD Model of the Transonic Dynamics Tunnel

This paper presents two Computational Fluid Dynamics (CFD) models of the flow in the NASA Langley Research Center Transonic Dynamics Tunnel (TDT). The TDT is a continuous-flow, closed circuit wind tunnel with a 16- by 16-foot slotted test section with cropped corners. The tunnel was originally built as the 19-ft Pressure Tunnel in 1938, but it was converted to the current transonic tunnel in the 1950s, with capabilities to use either air or heavy gas as the test medium. The first computational model describes the generation of an atmospheric-boundary-layer (ABL) profile inside the tunnel. An ABL, which includes both a wind profile and turbulence content, is one of the aerodynamic characteristics affecting the occurrence of wind-induced oscillations for a launch vehicle sitting on a pad. The challenging part of this analysis was modeling the turbulent flow inside the tunnel. This is due to the special ABL-generating hardware that was installed at the entrance of the TDT test section in order to change the downstream velocity profile and to introduce velocity fluctuations into the flow. The second CFD model builds on the computational aeroelastic results that were generated in support of the second Aeroelastic Prediction Workshop (AePW) for the NASA Benchmark Supercritical Wing (BSCW) configuration. During the AePW, the wing-only configuration (classical free-air model) was analyzed. In the current study, the flutter computations were conducted on the configuration as it was mounted in the TDT during the experiment. This includes the wing attached to the splitter plate that was attached to the wind-tunnel walls. The preliminary results show that the wind-tunnel walls marginally affect flutter prediction.

Pawel Chwalowski↗

Supporting hierarchical soil biogeochemical modeling: version 2 of the Biogeochemical Transport and Reaction model (BeTR-v2)

Abstract. Reliable soil biogeochemical modeling is a prerequisite for credible projections of climate change and associated ecosystem feedbacks. This recognition has called for frameworks that can support flexible and efficient development and application of new or alternative soil biogeochemical modules in Earth system models (ESMs). The the Biogeochemical Transport and Reaction model version 1 (BeTR-v1) code (i.e., CLM4-BeTR) is one such framework designed to accelerate the development and integration of new soil biogeochemistry formulations into ESMs and to analyze structural uncertainty in ESM simulations. With a generic reactive transport capability, BeTR-v1 can represent multiphase (e.g., gaseous, aqueous, and solid), multi-tracer (e.g., nitrate and organic carbon), and multi-organism (e.g., plants, bacteria, and fungi) dynamics. Here, we describe the new version, Biogeochemical Transport and Reaction model version 2 (BeTR-v2), which adopts more robust numerical solvers for multiphase diffusion and advection and coupling between biogeochemical reactions and improves code modularization over BeTR-v1. BeTR-v2 better supports different mathematical formulations in a hierarchical manner by allowing the resultant model be run for a single topsoil layer or a vertically resolved soil column, and it allows the model to be fully coupled with the land component of the Energy Exascale Earth System Model (E3SM). We demonstrate the capability of BeTR-v2 with benchmark cases and example soil biogeochemical (BGC) implementations. By taking advantage of BeTR-v2's generic structure integrated in E3SM, we then found that calibration could not resolve biases introduced by different numerical coupling strategies of plant–soil biogeochemistry. These results highlight the importance of numerically robust implementation of soil biogeochemistry and coupling with hydrology, thermal dynamics, and plants – capabilities that the open-source BeTR-v2 provides. We contend that Earth system models should strive to minimize this uncertainty by applying better numerical solvers.

58 GEOSCIENCES↗