Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal & statistical models”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Application of Land Surface Data Assimilation to Simulations of Sea Breeze Circulations

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite- observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The sea/land breeze is a well-documented mesoscale circulation that affects many coastal areas of the world including the northern Gulf Coast of the United States. The focus of this paper is to examine how the satellite assimilation technique impacts the simulation of a sea breeze circulation observed along the Mississippi/Alabama coast in the spring of 2001. The technique is implemented within the PSUNCAR MM5 V3-5 and applied at spatial resolutions of 12- and 4-km. It is recognized that even 4-km grid spacing is too coarse to explicitly resolve the detailed, mesoscale structure of sea breezes. Nevertheless, the model can forecast certain characteristics of the observed sea breeze including a thermally direct circulation that results from differential low-level heating across the land-sea interface. Our intent is to determine the sensitivity of the circulation to the differential land surface forcing produced via the assimilation of GOES skin temperature tendencies. Results will be quantified through statistical verification techniques.

Mackaro, Scott↗
NASA NTRS2006

An Analytical Thermal Model for Autonomous Soaring Research

A viewgraph presentation describing an analytical thermal model used to enable research on autonomous soaring for a small UAV aircraft is given. The topics include: 1) Purpose; 2) Approach; 3) SURFRAD Data; 4) Convective Layer Thickness; 5) Surface Heat Budget; 6) Surface Virtual Potential Temperature Flux; 7) Convective Scaling Velocity; 8) Other Calculations; 9) Yearly trends; 10) Scale Factors; 11) Scale Factor Test Matrix; 12) Statistical Model; 13) Updraft Strength Calculation; 14) Updraft Diameter; 15) Updraft Shape; 16) Smoothed Updraft Shape; 17) Updraft Spacing; 18) Environment Sink; 19) Updraft Lifespan; 20) Autonomous Soaring Research; 21) Planned Flight Test; and 22) Mixing Ratio.

Allen, Michael↗
DOE OSTI2026

Database and deep-learning scalability of anharmonic phonon properties by automated brute-force first-principles calculations

Understanding the anharmonic phonon properties of crystal compounds—such as phonon lifetimes and thermal conductivities—is essential for investigating and optimizing their thermal transport behaviors. These properties also impact optical, electronic, and magnetic characteristics through interactions between phonons and other quasiparticles and fields. In this study, we develop an automated first-principles workflow to calculate anharmonic phonon properties and build a comprehensive database encompassing more than 6500 inorganic compounds. Utilizing this dataset, we train a graph neural network model to predict thermal conductivity values and spectra from structural parameters, demonstrating a scaling law in which prediction accuracy improves with increasing training data size. High-throughput screening with the model enables the identification of materials exhibiting extreme thermal conductivities—both high and low. The resulting database offers valuable insights into the anharmonic behavior of phonons, thereby accelerating the design and development of advanced functional materials.

Ohnishi, Masato [University of Tokyo (Japan); Inst↗
NASA NTRS1984

The simulation of transient statistics of the Southern Hemispheric circulation by the GLAS seasonal cycle model: Preliminary results

The general circulation of the Southern Hemisphere is quite different from that of the Northern Hemisphere in many important ways. These include the barotropic nature of the stationary waves and the presence of a strong barotropic component to the mean zonal wind, the lack of a strong seasonal dependence of the transient eddies, and the dominant role played by eddies with periods less than 10 days compared to longer period fluctuations. Such differences attest to the importance of the altered nature of the orographic and thermal land-sea forcings in the Southern Hemisphere compared to the Northern Hemisphere. Some of the important features of the Southern Hemisphere circulation as simulated by the GLAS Seasonal Cycle Model (SCM) are presented. The geographical patterns of local variability and their seasonal shifts in the SCM are discussed and compared to observations.

Straus, D. M.↗
DOE OSTI2020

Issue Summary of INL Phase IV Transient Results for IAEA CRP on HTGR UAM Benchmark

This report details the Parallel and Highly Innovative Simulation for Idaho National Laboratory (INL) Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the transient core exercises defined for Phase IV of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase III models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the Uncertainty/Sensitivity Assessment (U/SA) of the prismatic modular high-temperature gas cooled reactor (MHTGR)-350 design, and specifically for Exercises IV-1 and IV-2 of the benchmark: the Control Rod Withdrawal (CRW) and Pressurised Loss of Cooling (PLOFC) events. The statistical U/SA methodology is implemented and demonstrated using the RAVEN code, based on perturbed cross-section libraries obtained from the SCALE/Sampler sequence. Uncertainties in nuclear data (cross-sections and the average number of neutrons produced per fission, 235U[¯v ]) lead to standard deviations (uncertainties of one s) of approximately 0.5% in the core eigenvalues of the MHTGR-350 and core models. For the coupled neutronics/thermal fluid model, local power density uncertainties up to 3.6% were observed in the colder regions of the core, while the local maximum fuel temperature uncertainties reached 1.5% for the models that included thermal fluid uncertainties. The addition of thermal fluid uncertainties dominated the impacts of nuclear data uncertainties in all cases. The main contributors to uncertainties in the power density and fuel temperatures during the transients were uncertainties in the reactor operating conditions (total power, inlet mass flow rate and inlet gas temperature). Variations in the bypass flows did not have significant impact on any of the output variables. For the nuclear data uncertainties it was found that the 235U(¯v ) / 235U(¯v ) covariance produced the largest sensitivities in terms of its impact on the eigenvalue and peak reactor power. It was also observed that the impact of any nuclear data uncertainties on the maximum fuel temperature was much less significant that the impact on eigenvalue and power. Another important finding was that although the use of eight or more energy groups is recommended for best-estimate HTGR simulation, two-group models produced acceptable uncertainty and sensitivity results for most FOMs. Since the statistical U/SA methodology is computationally expensive, and most transient solver requirements will scale directly with the number of energy groups, two energy groups could be used by HTGR developers during the early stages of design when larger uncertainty margins can be tolerated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
NASA NTRS2010

A Comparison of Tropical Storm (TS) and Non-TS Gust Factors for Assessing Peak Wind Probabilities at the Eastern Range

Peak wind speed is an important forecast element to ensure the safety of personnel and flight hardware at Kennedy Space Center (KSC) and the Cape Canaveral Air Force Station (CCAFS) in East-Central Florida. The 45th Weather Squadron (45 WS), the organization that issues forecasts for the KSC/CCAFS area, finds that peak winds are more difficult to forecast than mean winds. This difficulty motivated the 45 WS to request two independent studies. The first (Merceret 2009) was the development of a reliable model for gust factors (GF) relating the peak to the mean wind speed in tropical storms (TS). The second (Lambert et al. 2008) was a climatological study of non-TS cool season (October-April) mean and peak wind speeds by the Applied Meteorology Unit (AMU; Bauman et al. 2004) without the use of GF. Both studies presented their statistics as functions of mean wind speed and height. Most of the few comparisons of TS and non-TS GF in the literature suggest that non-TS GF at a given height and mean wind speed are smaller than the corresponding TS GF. The investigation reported here converted the non-TS peak wind statistics calculated by the AMU to the equivalent GF statistics and compared them with the previous TS GF results. The advantage of this effort over all previously reported studies of its kind is that the TS and non-TS data were taken from the same towers in the same locations. This eliminates differing surface attributes, including roughness length and thermal properties, as a major source of variance in the comparison. The goal of this study is two-fold: to determine the relationship between the non-TS and TS GF and their standard deviations (GFSD) and to determine if models similar to those developed for TS data in Merceret (2009) could be developed for the non-TS environment. The results are consistent with the literature, but include much more detailed, quantitative information on the nature of the relationship between TS and non-TS GF and GFSD as a function of height and mean wind speed.

Merceret, Francis J.↗
DOE OSTI2024

Robust finite-temperature many-body scarring on a quantum computer

Mechanisms for suppressing thermalization in disorder-free many-body systems, such as Hilbert space fragmentation and quantum many-body scars, have recently attracted much interest in foundations of quantum statistical physics and potential quantum information processing applications. However, their sensitivity to realistic effects such as finite temperature remains largely unexplored. Here, we have utilized IBM's Kolkata quantum processor to demonstrate an unexpected robustness of quantum many-body scars at finite temperatures when the system is prepared in a thermal Gibbs ensemble. We identify such robustness in the PXP model, which describes quantum many-body scars in experimental systems of Rydberg atom arrays and ultracold atoms in tilted Bose-Hubbard optical lattices. By contrast, other theoretical models which host exact quantum many-body scars are found to lack such robustness and their scarring properties quickly decay with temperature. Our study sheds light on the important differences between scarred models in terms of their algebraic structures, which impacts their resilience to finite temperature. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Statistical evaluation of thermal advection and stratification effects in scatterometer observations

The effects of thermal advection and atmospheric stratification are statistically evaluated using Seasat scatterometer observations as a data base. The results indicate that, whenever the surface winds or wind stress are related to the atmospheric pressure field, the appropriate stratification and baroclinic corrections should be applied. Without such corrections, errors of 15-25 percent are likely to arise in the surface fluxes computed from model low-level winds or pressure measurements.

Levy, Gad↗
DOE OSTI2020

Comparison of fixed charge and polarizable models for predicting the structural, thermodynamic, and transport properties of molten alkali chlorides

Results from extensive molecular dynamics simulations of molten LiCl, NaCl, KCl, and RbCl over a wide range of temperatures are reported. Comparison is made between the “Polarizable Ion Model” (PIM) and the non-polarizable “Rigid Ion Model” (RIM). Densities, self-diffusivities, shear viscosities, ionic conductivities, and thermal conductivities are computed and compared with experimental data. In addition, radial distribution functions are computed from ab initio molecular dynamics simulations and compared with the two sets of classical simulations as well as experimental data. Here, the two classical models perform reasonably well at capturing structural and dynamic properties of the four molten alkali chlorides, both qualitatively and often quantitatively. With the singular exception of liquid density, for which the PIM is more accurate than the RIM, there are few clear trends to suggest that one model is more accurate than the other for the four alkali halide systems studied here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

HIGH-FIDELITY SIMULATION OF SOOT FORMATION AND THERMAL RADIATION IN A LABORATORY-SCALE RICH-QUENCH-LEAN BURNER

High-fidelity simulations of a swirl-stabilized turbulent spray flame in a laboratory-scale aero-combustor have been performed to evaluate the predictability of state-of-the-art models in capturing soot formation. The simulations employ a complex chemical mechanism developed for Jet-A with PAH chemistry, coupled with the Hybrid Method of Moments (HMOM) soot model, and a Lagrangian dilute spray model for the fuel injection. Two simulations are performed to compare the results when thermal radiation is neglected or included in the solution with a mean spectral model. Modeling closures for the soot differential diffusion effects in mixture fraction space, as well as turbulence-radiation interaction are also evaluated using the data generated by the simulations. Given the degree of complexity of the simulation, the results showed good agreement with experimental measurements of the spatial distribution of the soot volume fraction ensemble average. A closer agreement with the experiment is observed when thermal radiation is included in the solution. Thermal radiation is observed to reduce the flame temperature and increase the flame intermittency, denoted by the increase in the temperature standard deviation in mixture fraction space. The reduction in temperature also leads to a reduction in PAH production and soot volume fraction. Turbulence is observed to have different effects on radiative emission depending on the mixture fraction. Turbulent scalar fluctuations significantly enhance radiative emission in fuel lean mixtures and can also play a role for fuel rich conditions. The statistical description of the turbulence-radiation interaction, previously proposed in the literature, was observed to correctly reproduce the high-fidelity results. Model coefficients were provided for swirl-stabilized flames. The soot differential diffusion model, previously proposed in the literature, based on the residual between the exact term and its model approximation, was also evaluated. The residual correction term further improved the agreement with exact differential diffusion term evaluated with the high-fidelity simulation data in mixture fraction space. The results suggest that the effective turbulent Lewis number can be equal to unity in simulations of turbulent non-premixed recirculating flames.

Soriano, Bruno [Sandia National Laboratories (SNL)↗
DOE OSTI2023

Spin polarization induced by magnetic field and the relativistic Barnett effect

First, I study the analogy between the magnetization of a material and the spin polarization of particles in a fluid. Using the relativistic version of the Barnett effect, i.e. the magnetization of a material induced by mechanical rotation, the spin polarization induced by thermal vorticity is obtained within a purely classical model, where spin is treated as an intrinsic magnetic moment and rotation is included as a non-inertial effect. I argue that since spin polarization induced by thermal vorticity can be obtained in a classical theory, it can not be dominated by quantum anomalies. Second, the spin polarization induced by magnetic field is obtained for a fluid at local thermal equilibrium using statistical quantum field theory. The obtained formula is valid beyond the weak field approximation and when contributions from the non-homogeneity of the magnetic field are small. The exact form of spin polarization is studied for a free Dirac field at global equilibrium, and, like magnetic susceptibility, it oscillates according to the de Haas - van Alphen effect. Finally, I briefly review how magnetic field contributes to the difference between the spin polarization of Λ and $\barΛ$ observed in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ensembles of Metastable States Govern Heterogeneous Catalysis on Dynamic Interfaces

Heterogeneous catalysis is at the heart of chemical industry. Being able to tune and design efficient catalysts for processes of interest is of an utmost importance, and for this, the molecular-level understanding of heterogeneous catalysts is the first step, and indeed a prime focus of the modern catalysis research. For a long time, a single most thermodynamically stable structure of the catalytic interface attained in reaction conditions had been envisioned as the reactive phase. However, some catalytic interfaces continue to undergo structural dynamics in the steady state, triggered by high temperatures, pressures, and binding and changing reagents. Among particularly dynamic interfaces are such widely-used catalysts as crystalline and amorphous surfaced supporting (sub-)nano metallic clusters. Recently, it became clear that this dynamic fluxionality causes the supported clusters to populate many distinct structural and stoichiometric states in catalytic conditions. Hence, the catalytic interface should be viewed as an evolving statistical ensemble of many (not one) structures. As such, every member in the ensemble contributes to the properties of the catalyst differently, and in proportion to its probability to be populated. This new notion flips the established paradigm and calls for new theory, modeling approaches, operando measurements, and updated design strategies. The statistical ensemble nature of surface-supported sub-nano cluster catalysts can be exemplified by oxide-supported and adsorbate-covered Pt, Pd, Cu, CuPd clusters, catalytic toward oxidative and non-oxidative dehydrogenation. They have access to a variety of 3D and quasi-2D shapes. The compositions of their thermal ensembles are dependent on the cluster size, leading to size-specific catalytic activities and the famous “every atom counts” phenomenon. The support and adsorbates affect catalyst structures, and state of the reacting species causes the ensemble to change in every reaction intermediate. The most stable member of the ensemble dominates the thermodynamic properties of the corresponding intermediate, whereas the kinetics can be determined by more active but less populated metastable catalyst states, and that suggests that many earlier studies might have overlooked the actual active sites. Both effects depend on the relative time-scales of catalyst restructuring and reaction dynamics. The catalyst may routinely operate off-equilibrium. Ensemble phenomena lead to surprising exceptions from established rules of catalysis, such as scaling relations, and Arrhenius behavior. Catalyst deactivation is also an ensemble property, and its extent of mitigation can be predicted through the new paradigm. These findings were enabled by advances in theory, such as global optimization and subsequent utilization of multiple local minima, and pathways sampling, as well as operando catalyst characterization. The fact that the per-site and per-species resolution is needed for the description and predicting of catalyst properties gives theory the central role in catalysis research, as most experiments provide ensemble-average information and cannot detect the crucial minority species that may be responsible for the catalytic activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The STAT7 Code for Statistical Propagation of Uncertainties in Steady State Thermal Hydraulics Analysis of Plate Fueled Reactors

The STAT7 software was developed to perform steady-state thermal hydraulic analyses. Application of the software is for non-power research and test reactors, including conversion to low-enriched uranium fuel of U.S. High-Performance Research Reactors such as MITR-II. Since it can be necessary to repeat analysis during fuel reloading, STAT7 accommodates flexibility in analyzing many realistic aspects of reactor fuel management. STAT7 uses a Monte Carlo approach to model common fabrication parameters and other key reactor analysis uncertainties required for research and test reactor thermal hydraulic analyses. These safety calculations are ultimately intended to protect against high fuel plate temperatures due to critical heat flux or departure from nucleate boiling or onset of flow instability; but additional margin is obtained by basing the limiting safety settings on avoiding onset of nucleate boiling. STAT7 can simultaneously analyze all of the axial nodes of all of the fuel plates and all of the coolant channels for one lateral stripe of a fuel element. The stripes run the length of the fuel, from the bottom to the top. Power splits are calculated for each axial node of each plate to determine how much of the power goes out each face of the plate. By running STAT7 multiple times, full core analysis can be performed by analyzing the margin to onset of nucleate boiling and onset of flow instability for each axial node of each stripe of each plate of each element in the core.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

Validation of chemistry models employed in a particle simulation method

The chemistry models employed in a statistical particle simulation method, as implemented in the Intel iPSC/860 multiprocessor computer, are validated and applied. Chemical relaxation of five-species air in these reservoirs involves 34 simultaneous dissociation, recombination, and atomic-exchange reactions. The reaction rates employed in the analytic solutions are obtained from Arrhenius experimental correlations as functions of temperature for adiabatic gas reservoirs in thermal equilibrium. Favorable agreement with the analytic solutions validates the simulation when applied to relaxation of O2 toward equilibrium in reservoirs dominated by dissociation and recombination, respectively, and when applied to relaxation of air in the temperature range 5000 to 30,000 K. A flow of O2 over a circular cylinder at high Mach number is simulated to demonstrate application of the method to multidimensional reactive flows.

Haas, Brian L.↗
NASA NTRS1998

Proposed Mars Surveyor Landing Sites in Northern Meridiani Sinus, Southern Elysium Planitia, and Argyre Planitia

Our objective is to propose two landing sites that the Mars Surveyor 2001 Lander and Athena Rover could go to on Mars that should meet the safety requirements of the spacecraft landing system and optimize surface operations (chiefly driven by power and communications requirements). An additional site within Argyre Planitia, initially proposed by Parker to the Mars Surveyor Landing Site program, is also proposed for potential consideration for post-2001 missions to Mars, as it is well outside the current latitude limits for the Athena Rover. All three sites are designed to be situated as close to a diversity of geologic units within a few kilometers of the landing site so that diversity can be placed in a geologic context. This objective is very different from the Mars Pathfinder requirement to land at a site with a maximum chance for containing a diversity of rocks within a few tens of meters of the lander. That requirement was driven by the Sojourner mobility limit of a few tens of meters. It can be argued that the Athena project, with its much larger mobility capability, might actually want to avoid such a site, because placing collected samples in geologic context would be difficult. While it has been argued, both before and after the Mars Pathfinder landing, that the provenance for local blocks may be determined by orbiter spectra, primarily from the MGS TES instrument, our ability to do so has yet to be demonstrated. Indeed, several months after conclusion of the Pathfinder mission, we have yet to reach a consensus on the composition of local materials. Our primary data set for selecting a landing site within the latitude and elevation constraints of the 2001 mission is the Viking Orbiter image archive. The site must be selected to place the landing ellipse so as to avoid obvious hazards, such as steep slopes, large or numerous craters, or abundant large knobs. For this purpose, we chose a resolution limit of better than 50 m/pixel. This necessarily excludes from the present study images from current and future orbiter spacecraft, until such data does become readily available. Within each proposed region, it may be possible to identify additional sites once these data become available. Second, the fine-component thermal inertia data, should be greater than about 5 or 6 cgs Units (10(exp -3) cal/sq cm s(exp -0.5)/K). Low thermal inertias imply dusty environments, which could pose a mobility hazard. Similarly, the albedo of the site should not be particularly high, which would also suggest dusty surfaces. Low albedos are preferred, as they often coincide with low Viking red:violet ratios and indicate less dusty surfaces. Next, the Modeled Block Abundance should also not be too high or too low. Based on the Viking Lander and Mars Pathfinder experiences, percentages of blocks should be on the order of 5-25%. Too many blocks could pose a hazard to the landing and mobility. Too few blocks could also indicate a dusty surface. Primary Landing Site: Northern Meridiani Sinus (Proposed by T. J. Parker and K., S. Edgett) Vital Statistics: (1) Latitude, Longitude: 0-3 N, 350-2 W. *Elevation (Viking): about0.5-1.5 Ian. (2) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbits 709A and 410B). Possible stereo coverage in region where two orbits overlap (probably small parallax angle, as these orbits are not listed in NASA Contractor Report 3501) (3) Albedo: about .18 -.26 (4) Block Abundance: 5-26% (5)Fine-Component Thermal Inertia: 5-9 cgs units This region consists of bright deposits similar to those described by Edgett et al, that also lie within a prominent dark albedo region. These deposits are flat-lying, to such a degree that they ramp against topography rather than draping over it. This led Edgett and Parker to suggest that they may be subaqueous sediments, possibly lacustrine or marine evaporites, laid down sometime from the late Noachian to middle Hesperian (age determination pending crater counts). A contact between this material and elevated, dissected highlands to the south was identified , and is described by Edgett et al. Our desire in proposing this landing site is to sample the edge of this deposit where it has been exposed through etching, presumably eolian deflation (the deposit, though in the highlands, is itself only lightly to moderately cratered). This should enable access to in situ stratigraphy. The actual landing site will be selected where slopes are not expected to be steep, such that the rover itself should be able to traverse them and sample layered materials on the way, either up or down the slope. Perhaps due to uncertainties at this time as to the friability or meter-scale roughness of the deposit, it might make sense to place the landing ellipse on the exhumed highland surface adjacent to the deflated margin of the deposit and plan on driving to the deposit rather than landing on it and driving downslope. This should also enable imaging the margin for evidence of layering should it prove too difficult to climb. A target ellipse on the highland surface should also allow Athena access to ancient Noachian highland materials, particularly if placed near crater ejecta or an inlier of knobby material. Secondary Landing Site: Southern Elysium Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 1.5-3.5 S, 195-198 W. (2) Elevation (Viking): -1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbit 725). Possible stereo coverage between images from beginning and end of orbit that overlap (probably small parallax angle) (4) Albedo: about .27-.28 (5) Block Abundance: 4-7% (6) Fine-Component Thermal Inertia: about 3 cgs units This region consists of eroded knobby material, probably of Noachian age, though much of the crater population has been destroyed, that is onlapped at a sharp contact by an extensive plains unit in southern Elysium Planitia that is Amazonian in age. The plains materials have been attributed to unusually low-viscosity flood lavas from fissures south of the Elysium volcanic rise, or to lacustrine materials associated with a large, Amazonian lake at the source of Marte Vallis. Parker and Schenk presented evidence in support of the latter interpretation, though they attributed the putative shore morphology to an embayment of a northern plains ocean into the southern Elysium region. Detailed examination of the margin of the deposit, showing erosion, not simply burial, of small crater rims and fluidized ejecta blankets, also points to lacustrine or marine sedimentation rather than volcanic plains burial. The plains surface exhibits a "crusty" appearance that many researchers have attributed to pressure ridges in lava flows. In a lacustrine context, they also resemble pressure ridges in desiccated evaporite deposits and salt-rimmed pools (now dry) similar in scale and morphology to spectacular, hundred meter-scale pool rims in alkaline Lake Natron, East African Rift. The eroded highland margin surface adjacent to these plains appears to be fairly smooth, even at 15 m/pixel. Isolated knob inliers are scattered from a few kilometers to several tens of "kilometers apart. Heights of the knobs have not been measured yet but, based on experience with similar features in the Pathfinder landing ellipse, are probably typically on the order of several tens of meters high and smaller, though some of the largest knobs in the region are probably up to a few hundred meters high. Two craters larger than a kilometer in diameter, with fluidized deposits, lie nearby the proposed landing site. Very high-resolution images from MOC should help to determine whether a landing site navigable by the Athena rover could be placed in this region. The space between knobs and craters is large enough to enable placement of a target landing ellipse between them but still provide access to one or more of them and to the margin of the Elysium plains material. Post-2001 Mars Surveyor Landing Site: Argyre Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 55-56 S, 41-43 W. (2) Elevation (Viking): 1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 40 m/pixel images (orbits 567B, 568B, and 569B). Excellent stereo coverage with large parallax angles over the entire landing site region, and much of central and southern Argyre. (4) Albedo: about .23-.24 (5) Block Abundance: No data (6) Fine-Component Thermal Inertia: No data The floors of both the Argyre and Hellas basins contain etched layered materials that are probably thick accumulations of channel or lacustrine sediments. The deposits in Hellas are much more eroded than those in Argyre, and Hellas lacks a channel outlet. Argyre is unique in that Uzboi Vallis flowed out of the basin, requiring overflow of a standing body of water within Argyre. This makes it the largest impact basin on Mars with channels both draining into it and flowing out from it. Hellas' channels may be catastrophic flood channels, whereas Argyre was fed by modest-scale valley networks, though the outlet at Uzboi Vallis was a catastrophic flood Highland craters and basins of this kind should be high-priority landing targets for missions intended to focus on the search for either prebiotic organic materials or even simple fossil microorganisms. Basins with internally-draining valley networks should be preferred over flood channels, as they could have provided the long-term influx of water favorable to the origin of life. (Catastrophic floods are not conducive to fossil preservation, due to their very short durations and high transportation energies). They also afford an opportunity to study the evolution of the planet's climate and volatiles during the period of time between the late Noachian and early Hesperian, when a drastic change from a proposed early warm, wet climate to one more closely resembling the modern environment is thought to have occurred. Large basin

Parker, T. J.↗
DOE OSTI2022

Supersaturation Variability from Scalar Mixing: Evaluation of a New Subgrid-Scale Model Using Direct Numerical Simulations of Turbulent Rayleigh–Bénard Convection

We report that supersaturation fluctuations in the atmosphere are critical for cloud processes. A nonlinear dependence on two scalars - water vapor and temperature - leads to different behavior than single scalars in turbulent convection. For modeling such multiscalar processes at subgrid scales (SGS) in large-eddy simulations (LES) or convection-permitting models, a new SGS scheme is implemented in CM1 that solves equations for SGS water vapor and temperature fluctuations and their covariance. The SGS model is evaluated using benchmark direct-numerical simulations (DNS) of turbulent Rayleigh–Bénard convection with water vapor as in the Michigan Tech Pi Cloud Chamber. This idealized setup allows thorough evaluation of the SGS model without complications from other atmospheric processes. DNS results compare favorably with measurements from the chamber. Results from LES using the new SGS model compare well with DNS, including profiles of water vapor and temperature variances, their covariance, and supersaturation variance. SGS supersaturation fluctuations scale appropriately with changes to the LES grid spacing, with the magnitude of SGS fluctuations decreasing relative to those at the resolved scale as the grid spacing is decreased. Sensitivities of covariance and supersaturation statistics to changes in water vapor flux relative to thermal flux are also investigated by modifying the sidewall conditions. Relative changes in water vapor flux substantially decrease the covariance and increase supersaturation fluctuations even away from boundaries.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Computer Design Tool for Ceramic Receivers

Gen3 solar receivers and other components will experience a combination of high temperatures and high stress caused by thermal stress and internal pressure. Under these conditions metallic components, even those manufactured from nickel-based superalloys, have poor reliability. Engineering ceramics could be a solution, as these materials have excellent high temperature strength. However, accurately assessing the reliability of a ceramic component operating in these conditions requires an entirely different approach compared to metallic materials. This paper describes the implementation of statistical models for evaluating the reliability of high temperature ceramic components in srlife – an open-source software package for estimating the life of high temperature concentrating solar power equipment. This new capability allows users to make fair comparisons between competing metallic and ceramic component designs and to accurately assess the plant efficiency and economics of ceramic solar receivers and other components.

Messner, Mark (ORCID:0000000200404385)↗