Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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,603 records · Page 89

Self-consistent theory of the Darrieus–Landau and Rayleigh–Taylor instabilities with self-generated magnetic fields

The Rayleigh–Taylor (RT) and Darrieus–Landau (DL) instabilities are studied in an inertial confinement fusion context within the framework of small critical-to-shell density ratio D R and weak acceleration regime, i.e., large Froude number Fr. The quasi-isobaric analysis in Sanz et al. is completed with the inclusion of non-isobaric and self-generated magnetic-field effects. The analysis is restricted to perturbation wavelengths k~1 larger than the conduction length scale at the ablation front, yet its validity ranges from wavelengths shorter and larger than the conduction layer width (distance between the ablation front and the critical surface). The use of a sharp boundary model leads to a single analytical expression of the dispersion relation encompassing both instabilities. The two new effects come into play by modifying the perturbed mass and momentum fluxes at the ablation front. The momentum flux (perturbed pressure at the spike) is the predominant stabilizing mechanism in the RT instability (overpressure) and the driving mechanism in the DL instability (underpressure). The non-isobaric effects notably modify the scaling laws in the DL limit, leading to an underpressure scaling as ~k –11/15 rather than ~k –2/5 obtained in the quasi-isobaric model. The magnetic fields are generated due to misalignment between pressure and density gradients (Biermann battery effect). They affect the hydrodynamics by bending the heat flux lines. Within the framework of this paper, they enhance ablation, resulting in a stabilizing effect that peaks for perturbation wavelengths comparable to the conduction layer width. Here, the combination of parameters D R Fr 2/3 defines the region of predominance of each instability in the dispersion relation. It is proven that the DL region falls outside of the parameter range in inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TNet: A Model-Constrained Tikhonov Network Approach for Inverse Problems

Deep learning (DL), in particular deep neural networks, by default is purely data-driven and in general does not require physics. This is the strength of DL but also one of its key limitations when applied to science and engineering problems in which underlying physical properties—such as stability, conservation, and positivity—and accuracy are required. DL methods in their original forms are often not capable of respecting the underlying mathematical models or achieving desired accuracy even in big-data regimes. On the other hand, many data-driven science and engineering problems, such as inverse problems, typically have limited experimental or observational data, and DL would overfit the data in this case. Leveraging information encoded in the underlying mathematical models, we argue, not only compensates for missing information in low data regimes but also provides opportunities to equip DL methods with the underlying physics, hence promoting better generalization. This paper develops a model-constrained DL approach and its variant TNet—a Tikhonov neural network—which are capable of learning not only information hidden in the training data but also in the underlying mathematical models to solve inverse problems governed by partial differential equations in low data regimes. We provide the constructions and some theoretical results for the proposed approaches for both linear and nonlinear inverse problems. Since TNet is designed to learn inverse solutions with Tikhonov regularization, it is interpretable: in fact it recovers Tikhonov solutions for linear cases while potentially approximating Tikhonov solutions for nonlinear inverse problems. We also prove that data randomization can enhance not only the smoothness of the networks but also their generalizations. Comprehensive numerical results confirm the theoretical findings and show that with even as little as 1 training data sample for one-dimensional (1D) deconvolution, 5 for an inverse 2D heat conductivity problem, 100 for inverse initial conditions for a time-dependent 2D Burgers’s equation, and 50 for inverse initial conditions for 2D Navier–Stokes equations, TNet solutions can be as accurate as Tikhonov solutions while being several orders of magnitude faster. Furthermore, this is possible owing to the model-constrained term, replications, and randomization.

97 MATHEMATICS AND COMPUTING↗

DEM Modeling and Validation of Pebble Bed Packing Using Chrono::GPU

Accurate prediction of pebble packing structure is important for pebble bed reactors because the spatial distribution of void fraction directly affects coolant flow, pressure drop, heat transfer, and neutronic behavior. However, experimentally validated DEM studies that directly evaluate local void-fraction structure in reactor-relevant pebble beds remain limited. In this work, the pebble bed experiment conducted at Missouri University of Science and Technology is simulated using the graphics processing unit (GPU)-based discrete element method (DEM) code Chrono::GPU. The study focuses on evaluating the ability of Chrono::GPU to reproduce the packing arrangement and void-fraction distribution of a randomly packed spherical pebble bed. The DEM results are first verified against established radial void-fraction correlations, including the Mueller and Vortmeyer-Schuster models, to assess the predicted bulk porosity, near-wall behavior, and oscillatory packing structure. The simulation is then verified against reference DEM data and validated against gamma-ray computed tomography (CT) experimental data at three axial locations. The Chrono::GPU results reproduce the main features of the experimental packing, including the high void fraction near the wall, the first near-wall trough, and the damped oscillatory radial profile caused by wall-induced ordering. Quantitative comparison with DEM data and the CT-based radial profiles shows good agreement, with mean absolute errors on the order of 0.07 and root-mean-square errors below 0.09 for the averaged profiles. These results demonstrate that Chrono::GPU can accurately capture the void-fraction structure of spherical pebble beds and provides a reliable DEM framework for future pebble bed reactor packing, recycling, and thermal-hydraulic studies.

97 - MATHEMATICS AND COMPUTING↗

Temperature distribution in a laser-heated diamond anvil cell as described by finite element analysis

Finite element analysis (FEA) is a powerful tool for numerically solving partial differential equations over complex geometries and is thus useful for analyzing heat transport in laser-heated diamond anvil cell (LHDAC) experiments. Our models expand on previously published simulations by calculating the volume-averaged temperatures of both the sample and insulation/pressure media under steady-state heating to determine the thermal pressure of the hot sample. Our goal is to produce an accurate relationship between the measured surface temperature of the absorbing sample and the temperature of the transparent insulating media, which is used to determine thermal pressure but susceptible to steep temperature gradients. We find that in doing so, our FEA models of temperature within the pressure/insulation media can differ from simplified estimates of temperature gradients by more than a factor of 2. We also explore temperature-dependent and temperature-independent thermal conductivity models and find that the volume-averaged temperatures differ by up to a factor of 1.3, forcing the predicted thermal pressures determined to also differ by up to a factor of 1.5 at a temperature of 2000 K at 50 GPa for neon. Higher temperatures exacerbate this difference. We also find that unintentional asymmetric sample insertion and sample heating, which are common in LHDAC experiments, do not have a first-order effect on volume-averaged temperatures. The FEA models, available in both Python and FlexPDE, are versatile across different sample geometries, materials, and heat source laser shapes.

Farah, Frederick↗

Multiscale Modeling of Microstructural Evolution in Fused-Coating Additive Manufacturing

A quantitative phase field model with adaptive-computational-domain-expansion (ACDE) method was developed to simulate the microstructural evolution in fused-coating additive manufacturing process. A heat transfer and fluid flow model was established to investigate the solidification behavior of the molten pool in a single deposited layer. The solidification parameters such as temperature gradient and solidification rate were calculated from the macroscopic model and then fed into the phase field model to get the microstructure at different positions within a single deposited layer. Columnar dendrites were formed at the bottom of the deposited layer, and equiaxed dendrites appeared in the upper part of the layer. The ACDE method is used to improve the simulation efficiency. The metallographic experiments and energy-dispersive spectrometer analysis for composition were conducted on the fused-coating specimen, and the results were in good accordance with the simulations.

36 MATERIALS SCIENCE↗

An Online Tool for Preliminary Design and Techno-Economic Analysis of District Geothermal Heating and Cooling Systems

District geothermal heating and cooling systems (DGHCS) have significant benefits for reducing energy consumption as well as building- and grid-level peak electric demand. Currently, no publicly available tools are available to effectively design and conduct techno-economic analysis of DGHCS. GeoWISE was originally developed for preliminary design and techno-economic analysis of geothermal heating and cooling systems in an individual commercial or residential building. This paper introduces recent upgrades of GeoWISE that allow users to design and conduct techno-economic analysis of DGHCS. Several new features are implemented in GeoWISE to allow selection and specification of multiple new or existing buildings. A database of information for over 125 million existing U.S. buildings was used in GeoWISE that allows users easily locate existing buildings of interest based on street addresses, and optionally edit information of the buildings (e.g., footprint, vintage, principal functions, number of floors, window-to-wall ratio). Unique energy simulation models of the selected buildings are then automatically created using the Automatic Building Energy Modeling (AutoBEM) and EnergyPlus simulations are performed to predict thermal loads of the buildings. A simplified DGHCS is then designed and simulated to predict its energy use. A central borehole heat exchanger (BHE) of the DGHCS is sized using the RowWise algorithm of GHEDesigner to meet the thermal loads within user-specified land areas for installing BHE. The upgraded GeoWISE reports the needed capacity of heating and cooling equipment in each building, design of the central BHE, energy consumption reduction, and energy cost saving resulting from using DGHCS compared with conventional HVAC systems. A case study is showcased using the upgraded GeoWISE to design and conduct techno-economic analysis of a simplified DGHCS.

Prem Anand Jayaprabha, Jyothis Anand [ORNL] (ORCID↗

Modeling of Aerobrake Ballute Stagnation Point Temperature and Heat Transfer to Inflation Gas

A trailing Ballute drag device concept for spacecraft aerocapture is considered. A thermal model for calculation of the Ballute membrane temperature and the inflation gas temperature is developed. An algorithm capturing the most salient features of the concept is implemented. In conjunction with the thermal model, trajectory calculations for two candidate missions, Titan Explorer and Neptune Orbiter missions, are used to estimate the stagnation point temperature and the inflation gas temperature. Radiation from both sides of the membrane at the stagnation point and conduction to the inflating gas is included. The results showed that the radiation from the membrane and to a much lesser extent conduction to the inflating gas, are likely to be the controlling heat transfer mechanisms and that the increase in gas temperature due to aerodynamic heating is of secondary importance.

Bahrami, Parviz A.↗

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗

Uncertainty Determination for Aeroheating in Uranus and Saturn Probe Entries by the Monte Carlo Method

The 2013-2022 Decaedal survey for planetary exploration has identified probe missions to Uranus and Saturn as high priorities. This work endeavors to examine the uncertainty for determining aeroheating in such entry environments. Representative entry trajectories are constructed using the TRAJ software. Flowfields at selected points on the trajectories are then computed using the Data Parallel Line Relaxation (DPLR) Computational Fluid Dynamics Code. A Monte Carlo study is performed on the DPLR input parameters to determine the uncertainty in the predicted aeroheating, and correlation coefficients are examined to identify which input parameters show the most influence on the uncertainty. A review of the present best practices for input parameters (e.g. transport coefficient and vibrational relaxation time) is also conducted. It is found that the 2(sigma) - uncertainty for heating on Uranus entry is no more than 2.1%, assuming an equilibrium catalytic wall, with the uncertainty being determined primarily by diffusion and H(sub 2) recombination rate within the boundary layer. However, if the wall is assumed to be partially or non-catalytic, this uncertainty may increase to as large as 18%. The catalytic wall model can contribute over 3x change in heat flux and a 20% variation in film coefficient. Therefore, coupled material response/fluid dynamic models are recommended for this problem. It was also found that much of this variability is artificially suppressed when a constant Schmidt number approach is implemented. Because the boundary layer is reacting, it is necessary to employ self-consistent effective binary diffusion to obtain a correct thermal transport solution. For Saturn entries, the 2(sigma) - uncertainty for convective heating was less than 3.7%. The major uncertainty driver was dependent on shock temperature/velocity, changing from boundary layer thermal conductivity to diffusivity and then to shock layer ionization rate as velocity increases. While radiative heating for Uranus entry was negligible, the nominal solution for Saturn computed up to 20% radiative heating at the highest velocity examined. The radiative heating followed a non-normal distribution, with up to a 3x variation in magnitude. This uncertainty is driven by the H(sub 2) dissociation rate, as H(sub 2) that persists in the hot non-equilibrium zone contributes significantly to radiation.

Palmer, Grant↗

The core and the magnetic field of Uranus

The existence of a magnetic field on Uranus is suggested by the tentative observation of 0.5 MHz pulses from the planet which are analogous to the Jovian decametric radiation. Analysis of a recent three-layer model leads to the conclusion that the electrical conductivity, physical state, and thermal gradients are such that a thermally driven hydromagnetic dynamo is probably operating in the core. An investigation is conducted to see whether the considered structure is suitable for the generation of a field. It is found that the existence of a magnetic field on Uranus can be understood as being a natural result of gravitational differentiation in a liquid, metallic core. The presence of a heat source in the core implies that the central temperature is approximately 8400 K rather than 7000 K in the original model

Torbett, M.↗

Augmentation of Stagnation Region Heat Transfer Due to Turbulence From a DLN Can Combustor

Heat transfer measurements have been made in the stagnation region of a flat plate with a circular leading edge. Electrically heated aluminum strips placed symmetrically about the leading edge stagnation region were used to measure spanwise averaged heat transfer coefficients. The maximum Reynolds number obtained, based on leading edge diameter, was about 100,000. The model was immersed in the flow field downstream of an approximately half scale model of a can-type combustor from a low NO(x), ground based power-generating turbine. The tests were conducted with room temperature air; no fuel was added. Room air flowed into the combustor through six vane type fuel/air swirlers. The combustor can contained no dilution holes. The fuel/air swirlers all swirled the incoming airflow in a counter clockwise direction (facing downstream). A 5-hole probe flow field survey in the plane of the model stagnation point showed the flow was one big vortex with flow angles up to 36' at the outer edges of the rectangular test section. Hot wire measurements showed test section flow had very high levels of turbulence, around 28.5 percent, and had a relatively large axial-length scale-to-leading edge diameter ratio of 0.5. X-wire measurements showed the turbulence to be nearly isotropic. Stagnation heat transfer augmentation over laminar levels was around 77 percent and was about 14 percent higher than predicted by a previously developed correlation for isotropic grid generated turbulence.

VanFossen, G. James↗

Sensitivity analysis of an automated fault detection algorithm for residential air-conditioning systems

The state of the art of fault detection and diagnosis (FDD) for residential air-conditioning systems is expensive and not yet amenable to widespread implementation. FDD for homes can significantly reduce utility costs, and increase the lifespan of the equipment. The cost barriers currently, however, make FDD for homes economically unviable for large scale implementation. In prior work, we offered a solution to reduce FDD costs by proposing an automated fault detection algorithm to serve as a screening step before more expensive FDD tests can be conducted. The algorithm uses only the home thermostat and local weather information to identify thermodynamic parameters and detect high-impact air-conditioning faults, including those that occur during equipment installation. We had tested the algorithm on a single EnergyPlus™ model of a home in Orlando, Florida. The thermodynamic parameter identification process is highly nonconvex involving several local optimal solutions. In this paper we propose a novel method to select the best model for fault detection from among the list of local optimal solutions to make the algorithm more robust to homes of different construction, without which the fault detection process would be infeasible. Another unique contribution of the paper is implementing the solution on real-world data. We also bring the algorithm closer to market by testing it on real-world data. We implement the algorithm on data obtained from experiments conducted by the Florida Solar Energy Center (FSEC) on a laboratory home equipped with a heat pump where faults were intentionally added for a period of seven months. The algorithm successfully detected an undercharge fault with 70.6% accuracy, concurrent duct leakage and undercharge faults with 85.2% accuracy, and duct leakage faults with 69.1% accuracy. A sensitivity analysis is also performed on EnergyPlus models of nine types of homes that vary in construction to demonstrate the robustness of the algorithm. Finally, the algorithm achieves an average accuracy of 71% for no-fault condition, 77% for 40% undercharge fault, and 76% for duct-leak fault.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Characterization of the indoor far-field aerosol transmission in a model commercial office building

As businesses gradually reopen and employees return to work, the potential spread of SARS-CoV-2 and its variants through airborne transmission via the heating, ventilation, and air-conditioning (HVAC) systems of commercial building raises concerns. Since the general practice in commercial buildings is to use low-efficiency air filters and given that indoor air is generally recycled, the degree to which cross-zone aerosol transmission occurs is of interest. To quantify the cross-zone aerosol transmission, experiments were conducted using a synthetic test aerosol in the five zones on the first floor of a model commercial office building at the Oak Ridge National Laboratory. Because the synthetic aerosol was tagged with fluorescent salt, the aerosol generated from the source zone can be distinguished from the background aerosols due to its unique fluorescent signal. Data from cross-zone campaigns showed that submicron-aerosol transmission was higher than the micron aerosols. In campaigns with doors closed, the submicron aerosol transmission was less than 16% and less than 11% for micron aerosol transmission. Finally, opening the interior doors that connecting different zones can significantly enhance the aerosol transmission for zones at the close proximity to the source, but has less impact on those farther away.

60 APPLIED LIFE SCIENCES↗

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗

Why Do Oceanic Nonlinearities Contribute Only Weakly to Extreme El Niño Events?

Abstract Extreme El Niño events have outsized global impacts and control the El Niño Southern Oscillation (ENSO) warm/cold phases asymmetries. Yet, a consensus regarding the relative contributions of atmospheric and oceanic nonlinearities to their genesis remains elusive. Here, we isolate the contribution of oceanic nonlinearities by conducting paired experiments forced with opposite wind stress anomalies in an oceanic general circulation model, which realistically simulates extreme El Niño events and oceanic nonlinearities thought to contribute to ENSO skewness (Tropical Instability Waves (TIWs), Nonlinear Dynamical Heating (NDH)). Our findings indicate a weak contribution of oceanic nonlinearities to extreme El Niño events in the eastern Pacific, owing to compensatory effects between lateral (NDH and TIWs) and vertical processes. These results hold across different vertical mixing schemes and modifications of the upper‐ocean heat budget mixed layer criterion. Our study reinforces previous research underscoring the pivotal role of atmospheric nonlinearities in shaping extreme El Niño events.

Geology↗

Device-Level Thermal Management of Gallium Oxide Field-Effect Transistors

The ultra-wide bandgap (~4.8 eV) and melt-grown substrate availability of β-Ga 2 O 3 gives promise to the development of next generation power electronic devices with dramatically improved size, weight, power, and efficiency over current state-of-the-art wide bandgap devices based on 4H-SiC and GaN. Also, with recent advancements made in GHz frequency RF applications, the potential for monolithic or heterogenous integration of RF and power switches has attracted researchers’ attention. Yet, it is expected that Ga 2 O 3 devices will suffer from self-heating due to the poor thermal conductivity of the material. Thermoreflectance thermal imaging and infrared thermography were used to understand the thermal characteristics of a MOSFET fabricated via homo-epitaxy. A 3D coupled electro-thermal model was constructed based on the electrical and thermal characterization results. The device model shows that a homo-epitaxial device suffers from an unacceptable junction temperature rise of ~1500°C under a targeted power density of 10 W/mm indicating the importance of employing device level thermal managements to individual Ga 2 O 3 transistors. The effectiveness of various active and passive cooling solutions was tested to achieve a goal of reducing the device operating temperature below 200°C at a power density of 10 W/mm. Results show that flip-chip hetero-integration is a viable option to enhance both the steady-state and transient thermal characteristics of Ga 2 O 3 devices without sacrificing the intrinsic advantage of high-quality native substrates. Also, it is not an active thermal management solution that entails peripherals requiring additional size and cost implications.

42 ENGINEERING↗

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗

Global 1 km land surface parameters for kilometer-scale Earth system modeling

Earth system models (ESMs) are progressively advancing towards the kilometer scale (“k-scale”). However, the surface parameters for land surface models (LSMs) within ESMs running at the k-scale are typically derived from coarse-resolution and outdated datasets. This study aims to develop a new set of global land surface parameters with a resolution of 1 km for multiple years from 2001 to 2020, utilizing the latest and most accurate available datasets. Specifically, the datasets consist of parameters related to land use and land cover, vegetation, soil, and topography. Differences between the newly developed 1 km land surface parameters and conventional parameters emphasize their potential for higher accuracy due to the incorporation of the most advanced and latest data sources. To demonstrate the capability of these new parameters, we conducted 1 km resolution simulations using the E3SM Land Model version 2 (ELM2) over the contiguous United States. Our results demonstrate that land surface parameters contribute to significant spatial heterogeneity in ELM2 simulations of soil moisture, latent heat, emitted longwave radiation, and absorbed shortwave radiation. On average, about 31 % to 54 % of spatial information is lost by upscaling the 1 km ELM2 simulations to a 12 km resolution. Using eXplainable Machine Learning (XML) methods, the influential factors driving the spatial variability and spatial information loss of ELM2 simulations were identified, highlighting the substantial impact of the spatial variability and information loss of various land surface parameters, as well as the mean climate conditions. The comparison against four benchmark datasets indicates that ELM generally performs well in simulating soil moisture and surface energy fluxes. The new land surface parameters are tailored to meet the emerging needs of k-scale LSM and ESM modeling with significant implications for advancing our understanding of water, carbon, and energy cycles under global change. The 1 km land surface parameters are publicly available at https://doi.org/10.5281/zenodo.10815170 (Li et al., 2024).

54 ENVIRONMENTAL SCIENCES↗