Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average model”

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 253 records · Page 14

Open-source simulation platform for air source heat pump integrated with thermal energy storage

Here, this article introduces a modular simulation platform for assessing thermal energy storage (TES) integrated with air source heat pumps (ASHP). The Python platform is an open-source library that includes classes for modeling air-air and air-water heat pumps, TES devices, and the heating load of residential buildings. To validate the ASHP model, the study utilized experimental data obtained from a commercial heat pump evaluated at the National Renewable Energy Laboratory (NREL). The results indicate a mean deviation of 0.7% for COP across the operating range, with a maximum relative deviation of 12.6%. In terms of system heating capacity, the model had an average deviation of 4.3% compared to experimental results, with a maximum deviation of 8.2%. Three classes were implemented for modeling distinct types of TES devices: a generic TES based on energy balances, a sensible isothermal water tank, and a stratified water tank. Details of the mathematical models are provided, along with their respective strengths and limitations. An example is provided showcasing the integration of a residential 10 kWh thermal storage unit with an ASHP operating in Denver, CO. The comparison of two different discharge criteria for the TES unit highlights the importance of control strategies in the system performance.

25 ENERGY STORAGE↗

Revisiting Power Systems Time-Domain Simulation Methods and Models

The changing nature of power systems dynamics is challenging present practices related to modeling and study of system-level dynamic behavior. While developing new techniques and models to handle the new modeling requirements, it is also critical to review some of the terminology used to describe existing simulation approaches and the embedded assumptions. This article provides a first-principles review of the simplifications and transformations commonly used in the formulation of time-domain simulation models. It introduces a taxonomy and classification of time-domain simulation models depending on their frequency bandwidth, network representation, and software availability. Furthermore, it focuses on the fundamental aspects of averaging techniques, and model reduction approaches that result in modeling choices, and discusses the associated challenges and opportunities of applying these methods in systems with large shares of Inverter Based Resources (IBRs). The article concludes with an illustrative simulation that compares the trajectories of an IBR-dominated system.

behavioral sciences↗

Real-time ensemble microalgae growth forecasting with data assimilation

Accurate short-range (e.g., 7-day) microalgae growth forecasts will be beneficial for both production and harvesting of microalgae. This study developed an operational microalgae growth forecasting system with ensemble data assimilation (DA). The forecasting system was validated against observed Monoraphidium minutum 26B-AM growth in two outdoor pond cultures located in Mesa, Arizona, U.S. We first examined the relative roles of uncertainty in the meteorological forecast and initial conditions (i.e., algal concentration at the time of forecast) in the microalgae 7-day forecast and found initial conditions dominated the microalgae forecasting skill, suggesting the importance of implementing DA to improve initial condition characterization. To correct the systematic bias in biomass simulations, we developed a particle filter with bias estimation (PFBE) DA method to estimate biases and correct the model forecast. We found the DA forecasting system could improve the 7-day microalgae forecasting skill by about 85% on average compared to model forecasts without DA. These results suggest the potential accuracy of biomass growth forecasts may be sufficient to inform real-time operational decisions, such as harvesting planning, for commercial-scale microalgae production.

59 BASIC BIOLOGICAL SCIENCES↗

A statistical model for the product energy distribution in reactions leading to prompt dissociation

Direct dynamics calculations have been performed for three reactions: C 3 H 8 + H → i-C 3 H 7 + H 2 , C 3 H 8 + H → n-C 3 H 7 + H 2 , and C 2 H 3 + O 2 → HCO + CH 2 O. The fraction of the population for the radical products that promptly dissociates is computed. Here the results for C 3 H 8 + H are qualitatively similar to previous results for C 3 H 8 + OH, but the new results exhibit a slightly higher branching fraction for prompt dissociation products, owing to the fact that a greater fraction of the internal energy in the transition state ends up in the radical. For C 2 H 3 + O 2 → HCO + CH 2 O, the fraction of HCO that promptly dissociates is in excess of 99%. Consequently, the main product for C 2 H 3 + O 2 at lower temperatures should be written as H + CO + CH 2 O and not HCO + CH 2 O. These results are then compared with four previous systems: CH 2 O + H → HCO + H 2 , CH 2 O + OH → HCO + H 2 O, C 3 H 8 + OH → i-C 3 H 7 + H 2 O, and C 3 H 8 + OH → n-C 3 H 7 + H 2 O. Based upon these seven system, several statistical models are presented. The goal of these statistical models is to predict the fraction of the transition state energy that ends up in the rovibrationally excited radical. On average, these statistical models provide an excellent prediction of product energy distribution. Consequently, these models can be used instead of costly trajectory simulations for predicting prompt radical dissociation for larger species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tools Assessing Performance (FY2022 Quarter 4 Report)

We were finally successful in acquiring TTU’s 200-meter tower data for the duration of the Isolated Building Wake Experiment (IBWE). This data will provide an independent measurement of the undisturbed wind profile including heights far above those that are typical of meteorological sensor networks. These measurements will characterize the mean and turbulence throughout the atmospheric surface layer, which will be particularly useful for validation of highfidelity computational fluid dynamics models. This dataset can be shared freely among the TAP team but not outside of it as TTU sells this data to the wind power industry. The journal article overview has been submitted to LANL’s publication review and will be submitted to the Journal of Wind Engineering and Industrial Aerodynamics once the review is complete. We are nearing completion on another two articles. The first discusses analyses of the evolution of the measurements of mean wind speed along the centerline of the building as a function of atmospheric stability. Rather than using individual time steps, all of the wind speed measurements are normalized using the upwind measurement at 2 building heights (H) and filtered by wind direction and atmospheric stability. The second discusses the development of the diffusive wake model and uses the averaged profiles from the first article to validate the behavior of the model. The portion of the section below on the diffusive wake model discussing comparisons with the averaged profiles, will be included in the validation section of this paper.

17 WIND ENERGY↗

Explaining word embeddings with perfect fidelity: a case study in predicting research impact

The best-performing approaches for scholarly document quality prediction are based on embedding models. In addition to their performance when used in classifiers, embedding models can also provide predictions even for words that were not contained in the labelled training data for the classification model, which is important in the context of the ever-evolving research terminology. Although model-agnostic explanation methods, such as Local interpretable model-agnostic explanations, can be applied to explain machine learning classifiers trained on embedding models, these produce results with questionable correspondence to the model. We introduce a new feature importance method, Self-Model Entities Rated (SMER), for logistic regression-based classification models trained on word embeddings. We show that SMER has theoretically perfect fidelity with the explained model, as the average of logits of SMER scores for individual words (SMER explanation) exactly corresponds to the logit of the prediction of the explained model. Quantitative and qualitative evaluation is performed through five diverse experiments conducted on 50,000 research articles (papers) from the CORD-19 corpus. In conclusion, through an AOPC curve analysis, we experimentally demonstrate that SMER produces better explanations than LIME, SHAP and global tree surrogates.

Coarse-grained models↗

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗

Modelled land use and land cover change emissions – a spatio-temporal comparison of different approaches

Abstract. Quantifying the net carbon flux from land use and land cover changes (fLULCC) is critical for understanding the global carbon cycle and, hence, to support climate change mitigation. However, large-scale fLULCC is not directly measurable and has to be inferred from models instead, such as semi-empirical bookkeeping models and process-based dynamic global vegetation models (DGVMs). By definition, fLULCC estimates are not directly comparable between these two different model types. As an important example, DGVM-based fLULCC in the annual global carbon budgets is estimated under transient environmental forcing and includes the so-called loss of additional sink capacity (LASC). The LASC results from the impact of environmental changes on land carbon storage potential of managed land compared to potential vegetation and accumulates over time, which is not captured in bookkeeping models. The fLULCC from transient DGVM simulations, thus, strongly depends on the timing of land use and land cover changes mainly because LASC accumulation is cut off at the end of the simulated period. To estimate the LASC, the fLULCC from pre-industrial DGVM simulations, which is independent of changing environmental conditions, can be used. Additionally, DGVMs using constant present-day environmental forcing enable an approximation of bookkeeping estimates. Here, we analyse these three DGVM-derived fLULCC estimations (under transient, pre-industrial, and present-day forcing) for 12 models within 18 regions and quantify their differences as well as climate- and CO2-induced components and compare them to bookkeeping estimates. Averaged across the models, we find a global fLULCC (under transient conditions) of 2.0±0.6 PgC yr−1 for 2009–2018, of which ∼40 % are attributable to the LASC (0.8±0.3 PgC yr−1). From 1850 onward, the fLULCC accumulated to 189±56 PgC with 40±15 PgC from the LASC. Around 1960, the accumulating nature of the LASC causes global transient fLULCC estimates to exceed estimates under present-day conditions, despite generally increased carbon stocks in the latter. Regional hotspots of high cumulative and annual LASC values are found in the USA, China, Brazil, equatorial Africa, and Southeast Asia, mainly due to deforestation for cropland. Distinct negative LASC estimates in Europe (early reforestation) and from 2000 onward in the Ukraine (recultivation of post-Soviet abandoned agricultural land), indicate that fLULCC estimates in these regions are lower in transient DGVM compared to bookkeeping approaches. Our study unravels the strong dependence of fLULCC estimates on the time a certain land use and land cover change event happened to occur and on the chosen time period for the forcing of environmental conditions in the underlying simulations. We argue for an approach that provides an accounting of the fLULCC that is more robust against these choices, for example by estimating a mean DGVM ensemble fLULCC and LASC for a defined reference period and homogeneous environmental changes (CO2 only).

54 ENVIRONMENTAL SCIENCES↗

Accelerating the discovery of battery electrode materials through data mining and deep learning models

The availability of crystalline materials databases allows for building accurate machine learning (ML) models that can accelerate the exploration of materials chemical space for energy storage applications. In this work, we screen all inorganic materials included in the Materials Project and AFLOW databases as potential metal-ion battery electrodes. We develop an efficient protocol to mine and screen raw data in current databases and provide a new database of electrode materials by considering pairs of charged and discharged electrodes. This effort leads to a new database with over 190,000 instances, in contrast to the original battery database which contains about 5000. The expanded battery data set is then used to build regression-based deep neural network models for predicting average voltages and percentage volume changes upon charging and discharging, which present improvements of at least 28% for target properties with respect to previous models, and are now able to predict anode electrodes (low voltage region) as well as electrodes that will not work in electrochemical cells (negative voltages), overcoming the challenges identified in previous ML models for battery electrodes. Additionally, a further screening of the expanded database itself allowed us to identify 35 novel electrode candidates with excellent battery performance metrics.

25 ENERGY STORAGE↗

High-fidelity forced convection simulations of the University of Wisconsin–Madison air-cooled reactor cavity cooling system

Among the various types of advanced nuclear technologies that exist, high-temperature gas-cooled reactors (HTGRs) stand out for their co-generation capabilities and exceptional passive safety systems. HTGRs rely on a reactor cavity cooling system (RCCS) to remove decay heat from the reactor pressure vessel (RPV) during the course of long-term transients. Several studies have focused on experimental and numerical modeling of the RCCS; however, the literature reveals a notable gap when it comes to high-fidelity simulations such as direct numerical simulations (DNS) or large eddy simulation (LES) models of the RCCS. An LES model was developed for forced convection inside the RCCS so as to obtain a high-fidelity representation of the RCCS at a scale never before attempted. Furthermore, the LES model was compared against three different Reynolds-averaged Navier–Stokes (RANS) models (Realizable k − ϵ, k − ω SST, and Reynolds stress model [RSM]). This comparison verified the robustness and accuracy of the turbulence modeling approach of the RANS models for the RCCS. Finally, based on the comparison between the LES and the RANS models, the k − ω SST RANS model of the experimental facility was developed and compared against the available experimental data.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Epidural anesthesia needle guidance by forward-view endoscopic optical coherence tomography and deep learning

Epidural anesthesia requires injection of anesthetic into the epidural space in the spine. Accurate placement of the epidural needle is a major challenge. To address this, we developed a forward-view endoscopic optical coherence tomography (OCT) system for real-time imaging of the tissue in front of the needle tip during the puncture. We tested this OCT system in porcine backbones and developed a set of deep learning models to automatically process the imaging data for needle localization. A series of binary classification models were developed to recognize the five layers of the backbone, including fat, interspinous ligament, ligamentum flavum, epidural space, and spinal cord. The classification models provided an average classification accuracy of 96.65%. During puncture, it is important to maintain a safe distance between the needle tip and the dura mater. Regression models were developed to estimate that distance based on the OCT imaging data. Based on the Inception architecture, our models achieved a mean absolute percentage error of 3.05% ± 0.55%. Overall, our results validated the technical feasibility of using this novel imaging strategy to automatically recognize different tissue structures and measure the distances ahead of the needle tip during the epidural needle placement.

60 APPLIED LIFE SCIENCES↗

Large Eddy Simulation of Externally Induced Ingress about an Axial Seal by Stator Vanes

Turbine inlet temperatures in advanced gas turbines could be as high as 2000 °C. To prevent ingress of this hot gas into the wheelspace between the stator and rotor disks, whose metals can only handle temperatures up to 850 °C, rim seals and sealing flows are used. This study examines the abilities of large eddy simulation (LES) based on the WALE subgrid model and Reynolds-averaged Navier–Stokes (RANS) based on the SST model in predicting ingress in a rotor–stator configuration with vanes but no blades, a configuration with experimental data for validation. Results were obtained for an operating condition, where the ratio of the external Reynolds number to the rotational Reynolds number is 0.538. At this operating condition, both LES and RANS were found to correctly predict the coefficient of pressure, C p , located downstream of the vanes and upstream of the seal, but only LES was able to correctly predict the sealing effectiveness. This shows C p by itself is inadequate in quantifying externally induced ingress. RANS was unable to predict the sealing effectiveness because it significantly under predicted the pressure drop in the hot gas path along the axial direction, especially about the seal region. This affected the pressure difference across the seal in the radial direction, which ultimately drives ingress.

42 ENGINEERING↗

Natural variability contributes to model–satellite differences in tropical tropospheric warming

Significance Climate models have, on average, simulated substantially more tropical tropospheric warming than satellite data, with few simulations matching observations. It has been suggested that this discrepancy arises because climate models are overly sensitive to greenhouse gas increases. Tropical tropospheric temperature trends from a large ensemble of simulations performed with a single climate model span a wide range that is solely due to natural climate variability. A subset of these simulations have warming rates in accord with satellite observations. Simulations with diminished tropical tropospheric warming due to climate variability exhibit characteristic patterns of surface warming that are similar to the observed record. Our results indicate that multidecadal variability can explain current model–observational differences in the rate of tropical tropospheric warming.

58 GEOSCIENCES↗

Super-resolution and segmentation deep learning for breast cancer histopathology image analysis

Traditionally, a high-performance microscope with a large numerical aperture is required to acquire high-resolution images. However, the images’ size is typically tremendous. Therefore, they are not conveniently managed and transferred across a computer network or stored in a limited computer storage system. As a result, image compression is commonly used to reduce image size resulting in poor image resolution. Here, we demonstrate custom convolution neural networks (CNNs) for both super-resolution image enhancement from low-resolution images and characterization of both cells and nuclei from hematoxylin and eosin (H&E) stained breast cancer histopathological images by using a combination of generator and discriminator networks so-called super-resolution generative adversarial network-based on aggregated residual transformation (SRGAN-ResNeXt) to facilitate cancer diagnosis in low resource settings. The results provide high enhancement in image quality where the peak signal-to-noise ratio and structural similarity of our network results are over 30 dB and 0.93, respectively. The derived performance is superior to the results obtained from both the bicubic interpolation and the well-known SRGAN deep-learning methods. In addition, another custom CNN is used to perform image segmentation from the generated high-resolution breast cancer images derived with our model with an average Intersection over Union of 0.869 and an average dice similarity coefficient of 0.893 for the H&E image segmentation results. Finally, we propose the jointly trained SRGAN-ResNeXt and Inception U-net Models, which applied the weights from the individually trained SRGAN-ResNeXt and inception U-net models as the pre-trained weights for transfer learning. The jointly trained model’s results are progressively improved and promising. We anticipate these custom CNNs can help resolve the inaccessibility of advanced microscopes or whole slide imaging (WSI) systems to acquire high-resolution images from low-performance microscopes located in remote-constraint settings.

60 APPLIED LIFE SCIENCES↗

A novel approach to partitioning evapotranspiration into evaporation and transpiration in flooded ecosystems

Reliable partitioning of micrometeorologically measured evapotranspiration (ET) into evaporation (E) and transpiration (T) would greatly enhance our understanding of the water cycle and its response to climate change related shifts in local-to-regional climate conditions and rising global levels of vapor pressure deficit (VPD). While some methods on ET partitioning have been developed, their underlying assumptions make them difficult to apply more generally, especially in sites with large contributions of E. Here, we report a novel ET partitioning method using artificial neural networks (ANNs) in combination with a range of environmental input variables to predict daytime E from nighttime ET measurements. The study uses eddy covariance data from four restored wetlands in the Sacramento-San Joaquin Delta, California, USA, as well as leaf-level T data for validation. The four wetlands vary in their vegetation make-up and structure, representing a range of ET conditions. The ANNs were built with increasing complexity by adding the input variable that resulted in the next highest average value of model testing R 2 across all sites. The order of variable inclusion (and importance) was: VPD > gap-filled sensible heat flux (H_gf) > air temperature (T air ) > friction velocity (u * ) > other variables. The model using VPD, H_gf, T air , and u * showed the best performance during validation with independent data and had a mean testing R 2 value of 0.853 (averaged across all sites, range from 0.728 to 0.910). In comparison to other methods, our ANN method generated T/ET partitioning results which were more consistent with CO 2 exchange data especially for more heterogeneous sites with large E contributions. Further, our method improves the understanding of T/ET partitioning. While it may be particularly suited to flooded ecosystems, it can also improve T/ET partitioning in other systems, increasing our knowledge of the global water cycle and ecosystem functioning.

59 BASIC BIOLOGICAL SCIENCES↗

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗

Comparing Designed Training Sets to Optimize Multivariate Regression Models for Pr, Nd, and Nitric Acid Using Spectrophotometry

Chemometric regression models were developed for the quantification of praseodymium (Pr, 0–1000 µg/mL), neodymium (Nd, 0–1000 µg/mL), and nitric acid (HNO 3 , 0.1–5 M) using spectrophotometry. Designed calibration sets were composed of 20 samples each: 10 model points and 10 lack-of-fit (LOF) points. The D-optimal designs effectively minimized the number of samples required to build models, and each design resulted in similar prediction performance, suggesting that statistical design of experiments can provide a reliable framework for selecting training set samples in three-variable systems. Partial least squares regression (PLSR) models were validated against a one-factor-at-a-time validation set composed of 125 samples (three variables, five levels). The top PLS-1 models resulted in average percent root mean square error of prediction error values of 3.5%, 1.7%, and 1.2% for Pr(III), Nd(III), and HNO 3 , respectively. Power set augmentations of the model and LOF samples were investigated to optimize the number of training set samples. PLSR models built using just required model points (10) had similar predictive capabilities as models including the LOF points (20) but with fewer samples. The number of validation samples was also varied systematically to learn how many samples are needed to validate regression models. This work addresses long-standing questions in the field of chemometrics to help make this approach amenable to the near-real-time quantification of hazardous species in remote settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗