Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature Design Methodology”

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 451 records · Page 25

Developing a Nuclear Quality Assurance Compliant Design Methodology for Neutronic Analysis of Xe-100 Design

The primary objective of this work is to develop a design methodology compliant with nuclear quality assurance standards for the Xe-100 neutronic design verification studies. To achieve this, a Monte Carlo model of the Xe-100 reactor was constructed using the exclusion principle, transformation technique, and universe-based level specification following Idaho National Laboratory (INL) NQA level-1 compliant standards and an NQA-1 compliant version of MCNP6. The model encompasses the entire reactor core structures, including the upper plenum, core region, and lower plenum sections, along with all sub-components. The active core section was represented using the spectral regions, each comprising a particular fuel composition and temperature averaged over the considered zone, calculated by X-energy using Very Superior Old Programs (VSOP). Additionally, a component-wise temperature map was implemented into the model, not only for the core region but also for the structural components. Temperature-dependent cross-section libraries, along with thermal scattering law libraries, generated using INL NQA-1 compliant version of NJOY21, were utilized for each isotope in the burnt fuel and the structural materials. Furthermore, the volume of each modeled component was estimated using a stochastic approach with the ray tracing method in MCNP and criticality calculations were performed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Predicting large-scale pool fire dynamics using an unsteady flamelet- and large-eddy simulation-based model suite

A low-Mach, unstructured, large-eddy-simulation-based, unsteady flamelet approach with a generalized heat loss combustion methodology (including soot generation and consumption mechanisms) is deployed to support a large-scale, quiescent, 5-m JP-8 pool fire validation study. The quiescent pool fire validation study deploys solution sensitivity procedures, i.e., the effect of mesh and time step refinement on capturing key fire dynamics such as fingering and puffing, as mesh resolutions approach O(1) cm. A novel design-order, discrete-ordinate-method discretization methodology is established by use of an analytical thermal/participating media radiation solution on both low-order hexahedral and tetrahedral mesh topologies in addition to quadratic hexahedral elements. Here, the coupling between heat losses and the flamelet thermochemical state is achieved by augmenting the unsteady flamelet equation set with a heat loss source term. Soot and radiation source terms are determined using flamelet approaches for the full range of heat losses experienced in fire applications including radiative extinction. The proposed modeling and simulation paradigm are validated using pool surface radiative heat flux, maximum centerline temperature location, and puffing frequency data, all of which are predicted within 10% accuracy. Simulations demonstrate that under-resolved meshes predict an overly conservative radiative heat flux magnitude with improved comparisons as compared to a previously deployed hybrid Reynolds-averaged Navier–Stokes/eddy dissipation concept-based methodology.

42 ENGINEERING↗

Experimental Studies of the Heat Transfer to RBCC Rocket Nozzles for CFD Application to Design Methodologies

Rocket thrusters for Rocket Based Combined Cycle (RBCC) engines typically operate with hydrogen/oxygen propellants in a very compact space. Packaging considerations lead to designs with either axisymmetric or two-dimensional throat sections. Nozzles tend to be either two- or three-dimensional. Heat transfer characteristics, particularly in the throat, where the peak heat flux occurs, are not well understood. Heat transfer predictions for these small thrusters have been made with one-dimensional analysis such as the Bartz equation or scaling of test data from much larger thrusters. The current work addresses this issue with an experimental program that examines the heat transfer characteristics of a gaseous oxygen (GO2)/gaseous hydrogen (GH2) two-dimensional compact rocket thruster. The experiments involved measuring the axial wall temperature profile in the nozzle region of a water-cooled gaseous oxygen/gaseous hydrogen rocket thruster at a pressure of 3.45 MPa. The wall temperature measurements in the thruster nozzle in concert with Bartz's correlation are utilized in a one-dimensional model to obtain axial profiles of nozzle wall heat flux.

Santoro, Robert J.↗

Spin-phonon relaxation from a universal ab initio density-matrix approach

Designing new quantum materials with long-lived electron spin states urgently requires a general theoretical formalism and computational technique to reliably predict intrinsic spin relaxation times. We present a new, accurate and universal first-principles methodology based on Lindbladian dynamics of density matrices to calculate spin-phonon relaxation time of solids with arbitrary spin mixing and crystal symmetry. This method describes contributions of Elliott-Yafet and D’yakonov-Perel’ mechanisms to spin relaxation for systems with and without inversion symmetry on an equal footing. We show that intrinsic spin and momentum relaxation times both decrease with increasing temperature; however, for the D’yakonov-Perel’ mechanism, spin relaxation time varies inversely with extrinsic scattering time. We predict large anisotropy of spin lifetime in transition metal dichalcogenides. The excellent agreement with experiments for a broad range of materials underscores the predictive capability of our method for properties critical to quantum information science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of an Advanced Multiphysics Simulation Capability for Radiant's Microreactor Design

Argonne National Laboratory and Idaho National Laboratory, through a Department of Energy Gateway for Accelerated Innovation in Nuclear Voucher, supported key analysis needs of Radiant related to (i) air jacket thermal fluid performance, (ii) evaluation of decay heat source terms defining air jacket technical requirements, and (iii) assessment of modeling methodologies employed for core analysis. All of these activities center on numerical simulation of various aspects of Kaleidos using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, the Cardinal multiphysics application, the OpenMC Monte Carlo code, and the Nek5000 computational fluid dynamics (CFD) code. This project builds upon an earlier Nuclear Energy Advanced Modeling and Simulation (NEAMS) Thermal-Hydraulic (T/H) Center of Excellence (CoE) project focused on initial demonstration of Cardinal multiphysics simulation of High Temperature Gas Reactors (HTGRs) and now focuses on Radiant’s Kaleidos concept.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Accelerated fission rate irradiation design, pre-irradiation characterization, and adaptation of conventional PIE methods for U-10Mo and U-17Mo

Metallic U alloys have high U density and thermal conductivity and thus have been explored since the beginning of nuclear power research. Alloys of U with modest amounts of Mo, such as U-10 wt % Mo (U-10Mo), are of particular interest because the γ-U crystal structure in this alloying addition shows prolonged stability in reactor service. Historically, radiation data on U-10Mo fuels were collected in Na fast reactors or lower temperature research reactor conditions, but little is known about irradiation behavior, particularly swelling and creep, at irradiation temperatures between 250 and 500°C. This work discusses the methodology and pre-irradiation characterization results from a U-Mo irradiation campaign performed in the High Flux Isotope Reactor at Oak Ridge National Laboratory. U-10Mo and U-17Mo samples irradiations are being completed at temperatures ranging from 250 to 500°C to three targeted fission densities between 2 × 10 20 and 1.5 × 10 21 fissions per cubic centimeter. Swelling measurement of the specimen sizes studied here required development and assessment of new methods for volume determination before and after irradiation. Laser profilometry and X-ray computation tomography (XCT) were used to provide preirradiation characterization of samples to determine the error and applicability of each to determine swelling following irradiation. These outcomes are contextualized through use of BISON simulations performed to assess the predicted expansion of U-Mo fuels subjected to the irradiation conditions of this work. Use of existing BISON fuel performance models predicted a maximum of 7% swelling under the irradiation conditions of this study. Pre-irradiation characterization revealed the as-cast U-Mo fuel samples were uniformly large-grained fully cubic U crystals with small U-C/N bearing precipitates and pores distributed throughout. Samples were found to contain a bulk porosity between .4 and 3% because of the casting process. Local porosity in areas far from large, interconnected pores was found by Slice-and-View to be under .2%. Nanometer-sized precipitates rich in C and N were identified in all samples, likely because of impurities during the fabrication process. Dendritic bands were also observed throughout the samples. These bands were characterized by variable Mo content that deviated from the overall Mo content by 2–3 wt %. No other microstructural features were correlated to these bands. Mechanical properties were found to be slightly strengthened compared to literature reports of bulk U-Mo fuels due to the nano-scale precipitates throughout the sample.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Trends in Performance and Characteristics of Ultra-Stable Oscillators for Deep Space Radio Science Experiments

Telecommunication systems of spacecraft on deep space missions also function as instruments for Radio Science experiments. Radio scientists utilize the telecommunication links between spacecraft and Earth to examine very small changes in the phase/frequency, amplitude, and/or polarization of radio signals to investigate a host of physical phenomena in the solar system. Several missions augmented the radio communication system with an Ultra-Stable Oscillator (USO) in order to provide a highly stable reference signal for oneway downlink. This configuration is used in order to enable better investigations of the atmospheres of the planets occulting the line-of-sight to the spacecraft; one-way communication was required and the transponders' built-in auxiliary oscillators were neither sufficiently stable nor spectrally pure for the occultation experiments. Since Radio Science instrumentation is distributed between the spacecraft and the ground stations, the Deep Space Network (DSN) is also equipped to function as a world-class instrument for Radio Science research. For a detailed account of Radio Science experiments, methodology, key discoveries, and the DSN's historical contribution to the field, see Asmar and Renzetti (1993). The tools of Radio Science can be and have also been utilized in addressing several mission engineering challenges; e.g., characterization of spacecraft nutation and anomalous motion, antenna calibrations, and communications during surface landing phases. Since the first quartz USO was flown on Voyager, the technology has advanced significantly, affording future missions higher sensitivity in reconstructing the temperature pressure profiles of the atmospheres under study as well as other physical phenomena of interest to Radio Science. This paper surveys the trends in stability and spectral purity performance, design characteristics including size and mass, as well as cost and history of these clocks in space.

Asmar, Sami↗

Development of an AQUA Based Near-Surface Parameter Retrieval

The production of a satellite based turbulent surface flux product relies critically upon the near-surface input parameters. Development of retrieval algorithms for the necessary near-surface variables of wind speed, specific humidity, air temperature, and sea surface temperature has proceeded relatively independent of each another until recently. The use of a neural network approach using Special Sensor Microwave/Imager (SSM/I) data in conjunction with a first guess sea surface temperature has led to successful retrieval of all parameters simultaneously. However, SSM/I frequencies lack inherent sensitivity to the sea surface temperature (SST). Recent studies have found improved air temperature and humidity retrievals can be obtained via inclusion of microwave sounding channels weighted in the lower troposphere. The inclusion of SSM/I-like frequencies as well as SST-sensitive microwave channels on AMSR-E along with AMSU-A sounding data onboard the AQUA platform provides an unique opportunity. That is the ability to provide near-simultaneous (in space and time) measurements allowing the retrieval of all the near-surface variables, including SST. This study shows results of a new algorithm designed to take advantage of the unique sampling ability of AQUA based sensors. Results from a neural network based methodology will be shown as compared to in-situ based observations of near-surface variables. Implications for creation of an AQUA based turbulent surface product are also discussed.

Roberts, Brent↗

System Modeling of Metabolic Heat Regenerated Temperature Swing Adsorption (MTSA) Subassembly for Prototype Design

This paper describes modeling methods for the three core components of a Metabolic heat regenerated Temperature Swing Adsorption (MTSA) subassembly: the sorbent bed, a sublimation (cooling) heat exchanger (SHX), and a condensing icing (warming) heat exchanger (CIHX). The primary function of the MTSA, removing carbon dioxide from a ventilation loop, is performed via the sorbent bed. The CIHX is used to heat the sorbent bed for desorption and to remove moisture from the ventilation loop while the SHX is alternately employed to cool the sorbent bed via sublimation of a spray of water at low pressure to prepare the reconditioned bed for the next cycle. This paper describes a system level model of the MTSA as developed in Thermal Desktop and SINDA/FLUINT including assumptions on geometry and physical phenomena, modeling methodology and relevant pa ra mete rizatio ns. Several areas of particular modeling interest are discussed. In the sorbent bed, capture of the translating CO2 saturation front and associated local energy and mass balance in both adsorbing and desorbing modes is covered. The CIHX poses particular challenges for modeling in SINDA/FLUINT as accounting for solids states in fluid submodels are not a native capability. Methods for capturing phase change and latent heat of ice as well as the transport properties across a layer of low density accreted frost are developed. This extended modeling capacity is applicable to temperatures greater than 258 K. To extend applicability to the minimum device temperature of 235 K, a method for a mapped transformation of temperatures from below the limit temperatures to some value above is given along with descriptions for associated material property transformations and the resulting impacts to total heat and mass transfer. Similar considerations are shown for the SHX along with assumptions for flow mechanics and resulting model methods for sublimation in a flow.

Bower, Chad↗

Advanced Materials for the Lunar Surface: Multiscale Computational Design of Refractory Alloys and Carbides

Emerging operational environments, such as the lunar surface, present novel challenges for NASA and drive the need for advanced materials in applications like fission surface power systems. To address these demands, computational materials science is rapidly evolving to augment or replace costly and hazardous empirical testing. Although materials selection at NASA remains predominantly experimentally driven, advanced simulation methodologies are being steadily integrated into the engineering lifecycle. This work details the application of multiscale simulation techniques—including first-principles calculations, CALPHAD, dislocation dynamics, and molecular dynamics—at NASA's Ames Research Center to evaluate advanced materials for extreme environments. First, we present contributions to the Space Nuclear Propulsion Project. Be-cause propellant channel coatings in nuclear thermal rockets must withstand high-pressure, high-temperature hydro-gen, optimizing these materials is critical. First-principles calculations were employed to establish a rigorous quantitative and qualitative understanding of the behavior of the refractory carbides ZrC, NbC, and their mixtures in high-enthalpy hydrogen environments. This necessitated the generation of high-fidelity thermodynamic models for both stoichiometric and carbon-depleted carbides, both with and without the presence of hydrogen. Furthermore, we highlight efforts under the Refractory Alloy Additive Manufacturing Build Optimization (RAAMBO) project, where existing and novel alloy compositions were assessed for additive manufacturing printability and subsequent performance in applications such as heat pipes and rocket nozzle extensions. This was accomplished through a comprehensive multiscale simulation framework that bridged the gap from the nanometer to the millimeter scale. Across both initiatives, rigorous validation against empirical data was prioritized. By systematically employing a verified and validated computational frame-work, we demonstrate how simulation effectively supports multidisciplinary engineering efforts, builds project-wide confidence, and drives critical materials development.

computational materials↗

Multi-Fidelity Scheme for Accelerating 4D Finite Element Analysis for Mircoreactors

Next-generation microreactors are currently being designed to be operated terrestrial and extraterrestrial for remote surface power production. These systems will provide an alternative source of carbon-free energy that is versatile and can be utilized for various applications. These applications of microreactors have prompted the usage of high-fidelity unstructured finite element (FE) based approaches to provide time-dependent solutions for multiple physics fields. Computing these high-fidelity (4-D) solutions for multiple design iterations, physics fields, and transient events requires an immense computational cost. These high-fidelity solutions are computational expensive and the cost can be reduced through the implementation of less accurate low-fidelity solutions. Unlike the current fleet of commercial nuclear reactors, these next-gen systems present challenges due to the material and physical limitations required. Due to these constraints, the brute force technique of parameterizing important system characteristics determines whether a design meets the project objective. A system such as a nuclear reactor could have thousands of design parameters that affect system performance. In order to analyze the entire parameter space of such a complex system would require millions of CPU hours and countless design iterations. This costly approach is not practical due to regulatory and budget limitations. In this proposal, an approach that utilizes hybrid high-fidelity and low-fidelity FE models to reduce the computational cost of evaluating these applications in 4-D will be presented. The accuracy of the high-fidelity model and the computational efficiency of the low-fidelity model are taken advantage of to produce a solution that closely resembles the full order high-fidelity solution. By utilizing an unconverged coarse FE mesh, operation limits such as temperature, structural loading, etc., can be evaluated in an accelerated fashion and then can be spatially interpolated onto a finer mesh. The resulting error arising from the coarse mesh can be mitigated with a discrepancy function that actively quantifies and corrects the error in the coarse solution. This discrepancy function can be periodically updated with high-fidelity calculations across the temporal domain thus requiring less iterations on the fine mesh. As a result, the computational cost can be reduced for the evaluation and design iteration of microreactors. The proposal for this research contains three sections: Section 2 provides a literature review, Section 3 outlines the methodology for the proposed multi-fidelity scheme, and Section 4 displays preliminary results of the proposed research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Interlaminar shear test method development for long term durability testing of composites

The high speed civil transport is a commercial aircraft that is expected to carry 300 passengers at Mach 2.4 over a range of more than 6000 nautical miles. With the existing commercial structural material technology (i.e., aluminum) the performance characteristics of the high speed civil transport would not be realized. Therefore there has been a concerted effort in the development of light weight materials capable of withstanding elevated temperatures for long duration. Thermoplastic composite materials are such candidate materials and the understanding of how these materials perform over the long term under harsh environments is essential to safe and effective design. The matrix dominated properties of thermoplastic composites are most affected by both time and temperature. There is currently an effort to perform short term testing to predict long term behavior of in-plane mechanical properties E22 (transverse modulus of elasticity) and G12 (shear modulus). Out-of-plane properties such as E33, G13, and G23 are inherently more difficult to characterize. This is especially true for the out-of-plane shear modulus G23 and hence there is no existing acceptable standard test method. Since G23 is the most matrix dominated property, it is essential that a test method be developed. A shear test methodology is developed to do just that. The test method, called the double notched specimen, along with the previously developed shear gage was tested at room temperature. Mechanical testing confirmed the attributes of the methodology. A finite element parametric study was conducted for specimen optimization. Moire interferometry, a high sensitivity laser optical method, was used for full-field analysis of the specimen. From this work, material parameters will be determined and thus enable the prediction of long term material behavior of laminates subjected to general loading states.

Ifju, Peter G.↗

ElectroCat: Platinum-Group-Metal-Free Engineered Framework Nano-Structure Catalyst (Final Technical Report)

Traditional synthesis methodologies that produce durable, highly active PGM-free electrocatalysts typically consist of a top-down approach utilizing a nitrogen-coordinating precursor, a metal precursor, and a support material. Nitrogen-coordinated transition metal complexes typically are absorbed onto a support material, which is then activated via high-temperature pyrolysis. Selection of the support and the nitrogen coordinating precursor has long been the focal point of efforts to maximize the number of accessible, catalytically active sites. In this project, GWE in collaboration with SRNL and NU will develop durable, highly active, low-cost, PGM-free electrocatalysts for polymer electrolyte membrane fuel cells designed and synthesized using a unique, bottom-up, rationally designed approach. Electrocatalysts are produced via polymerization reaction followed by a high-temperature pyrolysis step. Additionally, this project focuses on catalyst site and ORR mechanistic modeling to better understand the structure-property relationship.

08 HYDROGEN↗

A microprocessor application to a strapdown laser gyro navigator

The replacement of analog circuit control loops for laser gyros (path length control, cross axis temperature compensation loops, dither servo and current regulators) with digital filters residing in microcomputers is addressed. In addition to the control loops, a discussion is given on applying the microprocessor hardware to compensation for coning and skulling motion where simple algorithms are processed at high speeds to compensate component output data (digital pulses) for linear and angular vibration motions. Highlights are given on the methodology and system approaches used in replacing differential equations describing the analog system in terms of the mechanized difference equations of the microprocessor. Standard one for one frequency domain techniques are employed in replacing analog transfer functions by their transform counterparts. Direct digital design techniques are also discussed along with their associated benefits. Time and memory loading analyses are also summarized, as well as signal and microprocessor architecture. Trade offs in algorithm, mechanization, time/memory loading, accuracy, and microprocessor architecture are also given.

Giardina, C.↗

An Overview of Inter-comparison Methodologies for Terra and Aqua MODIS Calibration

With increasing efforts on data fusion and long-term climate data records (CDR) using observations made by multiple sensors, on the same or different platforms, their cross-calibration and validation work has become more and more important. The uncertainties of the climate models and data records depend not only on the calibration quality of individual sensors but also on their calibration consistency. This paper provides an overview of methodologies used by the MODIS Characterization Support Team (MCST) at NASA GSFC for the inter-comparison studies of Terra and Aqua MODIS on-orbit calibration. Each MODIS was built with a set of onboard calibrators (OBC) that include a blackbody (BB), a solar diffuser (SD), and a solar diffuser stability monitor (SDSM). The BB is primarily used for the thermal emissive bands (TEB) calibration and the SD/SDSM system for the reflective solar bands (RSB) calibration. Although the-instrument design and calibration approach are nearly identical for both Terra and Aqua MODIS and they all went through an extensive and similar pre-launch calibration and characterization process, still their on-orbit calibration consistency needs to be carefully examined and validated as many science products have been generated from observations made by both instruments. Methodologies discussed in this paper include inter-comparison studies using the Moon, a third sensor, and ground targets. Our results show that Terra and Aqua reflective solar bands and thermal emissive bands have been calibrated consistently with excellent long-term stabilities. For the 11 and 12 micrometers sea surface temperature (SST) bands, the calibration difference of Terra and Aqua MODIS is less than 0.2K.

Xiong, X.↗

Thermal Modeling and Testing of High-Temperature Refractory Ceramic Insulation Felts

Heat transfer in high-temperature, high-porosity, flexible refractory ceramic fibrous insulation felts is investigated. Heat transfer in these insulation materials consists of combined gas conduction, solid conduction, and radiation modes, with the precise theoretical modeling of the latter two modes being formidable. A semi-empirical model that requires inverse methods and steady-state thermal test data to infer some of the required model parameters is further developed in this study, and applied to five insulation materials for temperatures between 300 K and 1900 K. The steady-state thermal test setup at NASA Langley Research Center with recent modifications to increase its testing capability to 1900 K is discussed. Design considerations to ensure one-dimensional heat transfer in the test setup are described. Test data and corresponding thermal models for alumina and zirconia-based fibrous insulation felts are presented. Furthermore, test data and thermal models on two fibrous insulation samples containing additives to further suppress either radiation or gas conduction modes of heat transfer are presented. Previously published alumina-based insulation data are also re-processed using the updated modeling methodology. The significance of various heat transfer modes in typical insulation samples is discussed and used to provide general guidance on optimum insulation layups.

Kamran Daryabeigi↗

Instrumentation and Methodology Development for Mars Mission

The Mars environment comprises a dry, cold and low air pressure atmosphere with low gravity (0.38g) and high resistivity soil. The global dust storms that cover a large portion of Mars were observed often from Earth. This environment provides an idea condition for triboelectric charging. The extremely dry conditions on the Martian surface have raised concerns that electrostatic charge buildup will not be dissipated easily. If triboelectrically generated charge cannot be dissipated or avoided, then dust will accumulate on charged surfaces and electrostatic discharge may cause hazards for future exploration missions. The low surface temperature on Mars helps to prolong the charge decay on the dust particles and soil. To better understand the physics of Martian charged dust particles is essential to future Mars missions. We research and design two sensors, velocity/charge sensor and PZT momentum sensors, to detect the velocity distribution, charge distribution and mass distribution of Martian charged dust particles. These sensors are fabricated at NASA Kenney Space Center, Electromagnetic Physics Testbed. The sensors will be tested and calibrated for simulated Mars atmosphere condition with JSC MARS-1 Martian Regolith simulant in this NASA laboratory.

Chen, Yuan-Liang Albert↗