SEARCH · Engineering Papers

Results for “Subgrid scale”

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

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2024

Accurate Assessment of Land-Atmosphere Coupling in Climate Models Requires High Frequency Data Output

Land-atmosphere (L-A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centred climate anomalies. Local L-A coupling (LoCo) metrics capture relevant L-A processes, highlighting the impact of soil and vegetation states on surface flux partitioning, and the impact of surface fluxes on boundary layer (BL) growth, development, and entrainment of air above the BL. A primary goal of the Climate Process Team on Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L-A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land-atmosphere interactions span time scales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability of behavioural regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyse, and better understand L-A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behaviour. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant time scales, (ii) enable demonstration of inter-model spread in L-A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

Kirsten L. Findell↗

DOE OSTI2021

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1998

Influence of Soil Heterogeneity on Mesoscale Land Surface Fluxes During Washita '92

The influence of soil heterogeneity on the partitioning of mesoscale land surface energy fluxes at diurnal time scales is investigated over a 10(exp 6) sq km domain centered on the Little Washita Basin, Oklahoma, for the period June 10 - 18, 1992. The sensitivity study is carried out using MM5/PLACE, the Penn State/NCAR MM5 model enhanced with the Parameterization for Land-Atmosphere-Cloud Exchange or PLACE. PLACE is a one-dimensional land surface model possessing detailed plant and soil water physics algorithms, multiple soil layers, and the capacity to model subgrid heterogeneity. A series of 12-hour simulations were conducted with identical atmospheric initialization and land surface characterization but with different initial soil moisture and texture. A comparison then was made of the simulated land surface energy flux fields, the partitioning of net radiation into latent and sensible heat, and the soil moisture fields. Results indicate that heterogeneity in both soil moisture and texture affects the spatial distribution and partitioning of mesoscale energy balance. Spatial averaging results in an overprediction of latent heat flux, and an underestimation of sensible heat flux. In addition to the primary focus on the partitioning of the land surface energy, the modeling effort provided an opportunity to examine the issue of initializing the soil moisture fields for coupled three-dimensional models. For the present case, the initial soil moisture and temperature were determined from off-line modeling using PLACE at each grid box, driven with a combination of observed and assimilated data fields.

Jasinski, Michael F.↗

NASA NTRS2015

Mapping Biophysical Parameters for Land Surface Modeling over the Continental US Using MODIS and Landsat

In terms of the space cities occupy, urbanization appears as a minor land transformation. However, it permanently modifies land's ecological functions, altering its carbon, energy, and water fluxes. It is therefore necessary to develop a land cover characterization at fine spatial and temporal scales to capture urbanization's effects on surface fluxes. We develop a series of biophysical vegetation parameters such as the fraction of photosynthetically active radiation, leaf area index, vegetation greenness fraction, and roughness length over the continental US using MODIS and Landsat products for 2001. A 13-class land cover map was developed at a climate modeling grid (CMG) merging the 500mMODIS land cover and the 30m impervious surface area from the National Land Cover Database. The landscape subgrid heterogeneity was preserved using fractions of each class from the 500 m and 30 m into the CMG. Biophysical parameters were computed using the 8-day composite Normalized Difference Vegetation Index produced by the North American Carbon Program. In addition to urban impact assessments, this dataset is useful for the computation of surface fluxes in land, vegetation, and urban models and is expected to be widely used in different land cover and land use change applications.

Biophysical↗

DOE OSTI2023

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks. Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Bridging the Gap between Cosmological Simulations with Graph Neural Networks and Domain Adaptation

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross dataset tasks (up to 28% better relative error and up to almost an order of magnitude better χ 2 ). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks (up to $28\%$ better relative error and up to almost an order of magnitude better $\chi^2$). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window

ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.

Sharbaf, Z. (ORCID:0009000450545946)↗

DOE OSTI2024

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

NASA NTRS2024

Regimes of Cloud Vertical Structure From Active Observations

We introduce a new type of cloud class, which we call “active cloud regime” (ACR), owing to its provenance from active (lidar and cloud radar) spaceborne cloud observations. ACRs intend to provide a climatological description based on cloud vertical structure (CVS) of the most prevalent monthly CVS mixtures encountered at large spatial scales of ~400 km. ACRs are thus a way to create a gridded dataset of a vertically-resolved cloud mask that can facilitate joint analysis with other gridded datasets. The detailed 2D cloud mask comes from the 2B-CLDCLASS-LIDAR CloudSat dataset fusing CALIPSO (lidar) and CloudSat (cloud radar) cloud detections. We show that the global classification of cloudiness under the ACR framework provides valuable insights on how the world’s dominant cloud systems regulate the two major components of atmospheric energetics, precipitation and radiative cooling. NASA’s GEOS model allows us to demonstrate the feasibility of applying the ACR concept in Earth System Models that have the capability to produce subgrid cloudiness obeying pre-specified vertical overlap rules. Comparison of observed and simulated ACRs provides thus another means to assess the realism of modeled clouds.

cloud vertical structure↗

DOE OSTI2021

A Bimodal Diagnostic Cloud Fraction Parameterization. Part I: Motivating Analysis and Scheme Description

Cloud fraction parameterizations are beneficial to regional, convection-permitting numerical weather prediction. For its operational regional midlatitude forecasts, the Met Office uses a diagnostic cloud fraction scheme that relies on a unimodal, symmetric subgrid saturation-departure distribution. This scheme has been shown before to underestimate cloud cover and hence an empirically based bias correction is used operationally to improve performance. This first of a series of two papers proposes a new diagnostic cloud scheme as a more physically based alternative to the operational bias correction. The new cloud scheme identifies entrainment zones associated with strong temperature inversions. Additionally, for model grid boxes located in this entrainment zone, collocated moist and dry Gaussian modes are used to represent the subgrid conditions. The mean and width of the Gaussian modes, inferred from the turbulent characteristics, are then used to diagnose cloud water content and cloud fraction. It is shown that the new scheme diagnoses enhanced cloud cover for a given gridbox mean humidity, similar to the current operational approach. It does so, however, in a physically meaningful way. Using observed aircraft data and ground-based retrievals over the southern Great Plains in the United States, it is shown that the new scheme improves the relation between cloud fraction, relative humidity, and liquid water content. An emergent property of the scheme is its ability to infer skewed and bimodal distributions from the large-scale state that qualitatively compare well against observations. A detailed evaluation and resolution sensitivity study will follow in Part II.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A Novel Method for Diagnosing Land–Atmosphere Coupling Sensitivity in a Single-Column Model

The response of boundary layer properties and cloudiness to changes in surface evaporative fraction (EF) is investigated in a single-column model to quantify the locally coupled impact of subgrid surface variations on the atmosphere during summer. Sensitive coupling days are defined when the model atmosphere exhibits large variations across a range of EFs centered on the analyzed value. Coupling sensitivity exists as both positive feedback (cloudiness increases with EF) and negative feedback (clouds increase with decreasing EF) regimes. The positive regime manifests in shallow convection situations, which are capped by a strengthened inversion and subsidence, restricting the vertical extent of convection to just above the boundary layer. Surfaces with larger EF (greater surface latent heat flux) can inject more moisture into the vertically confined system, lowering the cloud base and an increasing cloud liquid water path (LWP). Negative feedback regimes tend to manifest when large-scale deep convection, such as from mesoscale convective systems and fronts, is advected through the domain, where convection strengthens over surfaces with a lower EF (greater surface sensible heat flux). The invigoration of these systems by the land surface leads to an increase in LWP through strengthened updrafts and stronger coupling between the boundary layer and the free atmosphere. These results apply in the absence of heterogeneity-induced mesoscale circulations, providing a one-dimensional dynamical perspective on the effect of surface heterogeneity. Finally, this study provides a framework of intermediate complexity, lying between parcel theory and high-resolution coupled land–atmosphere modeling, and therefore isolates the relevant first-order processes in land–atmosphere interactions.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1999

Statistics of Cirrus Horizontal Inhomogeneity in the Southern Great Plains

Variability of cloud properties on scales smaller than that of a GCM grid is potentially important both for realistic parameterizations of microphysical processes and for the prediction of the large-scale radiative effects of clouds, have suggested that a simple model of marine stratocumulus variability, based on the assumption of Gaussian variation statistics of cloud depth, can explain the liquid water path histogram shapes observed in Landsat data. In advance of ARM SGP MMCR ice water path climatologies, we have examined aircraft ice water content statistics for cirrus clouds observed over Coffeyville, Kansas during FIRE 2. We find similar associations of histogram shape and cloud cover for these clouds, and we show that a simple modification of the model for cirrus combined with observed mean cloud depths, their standard deviations, and ambient thermodynamic conditions predicts both the histogram shape and cirrus cloud cover fairly well. This suggests that subgrid variability of cloud properties may be similar for vastly different cloud types, and that a universal parameterization of the effects of subgrid variability in GCMs as a function of only a few parameters may be a realistic goal.

Smith, S. A.↗

DOE OSTI2023

Region and cloud regime dependence of parametric sensitivity in E3SM atmosphere model

Abstract The Department of Energy (DOE)’s Energy Exascale Earth System Model (E3SM), including its atmosphere model (EAM), has many relatively new features. In a previous study we conducted a systematic parametric sensitivity analysis for EAM based on short, perturbed parameter ensemble (PPE) simulations, mainly focusing on global mean climate features and metrics. While parameter values in global climate models are generally invariant in space and time, model response to parameters perturbation may vary by regions and climate regimes, which motivates the need to better understand the EAM model behaviors and physics at regional scale and process level. In this study, using the same set of PPE simulations and a similar sensitivity analysis framework, we identify parameters that cause largest sensitivities over different regions and compare model responses in fast atmospheric processes to the parameters across different cloud regimes for several important cloud-related fidelity metrics. We find that cloud forcing has opposite response to some parameters over mid-latitude vs. tropical land. We also analyze how the parametric sensitivity varies as stratocumulus transitions to shallow convection and to deep convection over ocean. Low cloud forcing and shortwave cloud forcing in the subtropical eastern Pacific are most sensitive to the parameters controlling the width of the probability density function (PDF) of the subgrid vertical velocity ( w’ ) ( gamma ) and the damping of the w’ skewness ( c8 ) near the coast but become more sensitive to the parameter affecting the damping of the w’ variance ( c1 ) further offshore. Detailed interpretation of the spatial dependence of parametric sensitivity is provided. We also investigate how the parametric sensitivity evolves with prediction duration. This study improves our process-level understanding of cloud physics and parameterization and provides insights for developing more advanced regime-aware parameterization schemes in global climate model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Evaluation of Autoconversion Representation in E3SMv2 Using an Ensemble of Large-Eddy Simulations of Low-Level Warm Clouds

In numerical atmospheric models that treat cloud and rain droplet populations as separate condensate categories, precipitation initiation in warm clouds is often represented by an autoconversion rate (Au), which is the rate of formation of new rain droplets through the collisions of cloud droplets. Being a function of the cloud droplet size distribution (DSD), the local Au is commonly parameterized as a function of DSD moments: cloud droplet number (n c ) and mass (q c ) concentrations. When applied in a large-scale model, the grid-mean Au must also include a correction, or enhancement factor, to account for the horizontal variability of the cloud properties across the model grid. In this study, we evaluate the Au representation in the Energy Exascale Earth System Model version 2 (E3SMv2) climate model using large-eddy simulations (LES), which explicitly resolve cloud droplet spectra, and therefore the local Au, as well as its spatial variability. The analysis of an ensemble of warm low-level cloud cases shows that the E3SMv2 formulation represents the Au reasonably well compared to the horizontally averaged explicitly computed rate from LES. The agreement, however, comes from a combination of an underestimated E3SM-tuned local Au rate and an overestimated subgrid cloud variability enhancement factor. The latter bias is traced to neglecting the horizontal variability of n c and its co-variability with q c in parameterizing the grid-mean Au.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Near‐Surface Hydrology and Soil Properties Drive Heterogeneity in Permafrost Distribution, Vegetation Dynamics, and Carbon Cycling in a Sub‐Arctic Watershed

Abstract Discontinuous permafrost environments exhibit strong spatial heterogeneity at scales too small to be driven by weather forcing or captured by Earth System Models. Here we analyze effects of observed spatial heterogeneity in soil and vegetation properties, hydrology, and thermal dynamics on ecosystem carbon dynamics in a watershed on the Seward Peninsula in Alaska. We apply a Morris global sensitivity analysis to a process‐rich, successfully tested terrestrial ecosystem model (TEM), ecosys , varying soil properties, boundary conditions, and weather forcing. We show that landscape heterogeneity strongly impacts soil temperatures and vegetation composition. Snow depth, O‐horizon thickness, and near‐surface water content, which vary at scales of O(m), control the soil thermal regime more than an air temperature gradient corresponding to a 140 km north–south distance. High shrub productivity is simulated only in talik (perennially unfrozen) soils with high nitrogen availability. Through these effects on plant and permafrost dynamics, landscape heterogeneity impacts ecosystem productivity. Simulations with near‐surface taliks have higher microbial respiration (by 78.0 gC m −2 yr −1 ) and higher net primary productivity (by 104.9 gC m −2 yr −1 ) compared to runs with near‐surface permafrost, and simulations with high shrub productivity have outlying values of net carbon uptake. We explored the prediction uncertainty associated with ignoring observed landscape heterogeneity, and found that watershed net carbon uptake is 60% larger when heterogeneity is accounted for. Our results highlight the complexity inherent in discontinuous permafrost environments and demonstrate that missing representation of subgrid heterogeneity in TEMs could bias predictions of high‐latitude carbon budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Adaptive Multirate Mass Transfer (aMMT) Model: A New Approach to Upscale Regional‐Scale Transport Under Transient Flow Conditions

Abstract The long‐term evaluation of regional‐scale groundwater quality needs efficient upscaling methods for transient flow. Upscaling techniques, such as the Multirate Mass Transfer (MRMT) method with constant upscaling parameters, have been used for transport with steady‐state flow, yet the upscaling parameters (i.e., rate coefficients) may be time dependent. This study proposed and validated an adaptive MRMT (aMMT) method by allowing the mass transfer coefficients in MRMT to change with the flow field. Advective‐dispersive contaminant transport simulated in a 3‐D heterogeneous medium was used as a reference solution. Equivalent transport under homogeneous flow conditions was evaluated by applying the MRMT and aMMT models for upscaling. The relationship between mass transfer coefficients and flow rates was fitted under steady‐state flow driven by various hydraulic gradients. A power law relationship was obtained, which was then used to update the mass transfer coefficients in each stress period under transient flow conditions in the aMMT method. Results indicated that for advection‐dominated transport, both the MRMT and aMMT methods can upscale the anomalous transport dynamics affected by subgrid heterogeneity under transient flow conditions. Whereas for diffusion‐dominated systems, the MRMT model failed to capture the tails of tracer breakthrough curves after the boundary condition changed, but the results from the aMMT model were significantly improved. However, if the overall flow direction changed, both MRMT and aMMT failed to represent the breakthrough curve tail generated by the heterogeneous system. The results point toward a promising path for upscaling transport in complex aquifers with transient flow.

Guo, Zhilin↗

DOE OSTI2020

Laser propagation in a subcritical foam: Subgrid model

Here, we present a subgrid model for laser propagation in a subcritical foam. Our model describes the expansion of laser-irradiated foam elements that are below the resolution of the simulation grid and predicts the plasma conditions that result from burning down the foam. Our model can be included as a module within a larger multiphysics code, and we have implemented it within the code pF3D, which is used for simulating a laser-plasma interaction. The model predicts a reduced propagation velocity for a laser through a subcritical foam compared to simulating that foam as a homogeneous gas. This is attributed to the laser energy that goes into burning down the foam microstructure. We compare our model against experimental data by simulating a 2 mg/cc SiO 2 foam shot performed at the Janus laser facility at the Lawrence Livermore National Laboratory. pF3D simulations with the foam model predict hot ion temperatures. This leads to a reduction in the level of stimulated Brillouin scattering (SBS), bringing the simulated level of SBS into agreement with the data. Intensity fluctuations at the foam front due to laser speckles and refraction result in ion temperature fluctuations when the foam burns down. These drive long-lived electron density fluctuations on scales that are large compared to the pore size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗