Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Cryogenic Material Characterization”

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.

1,489 records · Page 15

Research and Development Related to the Commercialization of Perovskite Radiation Detectors (CRADA Final Report)

This work will focus on the development of perovskite radiation detectors fabricated with nanocrystal, thick polycrystalline films, or single crystal materials. This will harness NREL’s instrumentation and expertise in a variety of research areas including but not limited to metal halide perovskites, device characterization, device packaging, accelerated lifetime testing, and ink chemistry development.

14 SOLAR ENERGY

CatTestHub: A benchmarking database of experimental heterogeneous catalysis for evaluating advanced materials

The ability to quantitatively compare newly evolving catalytic materials and technologies is hindered by the widespread availability of catalytic data collected in a consistent manner. While certain catalytic chemistries have been widely studied across decades of scientific research, quantitative comparisons based on literature information is hindered by variability in reaction conditions, types of reported data, and reporting procedures. Here, we present CatTestHub, an open-access database dedicated to benchmarking experimental heterogeneous catalysis data. Combining systematically reported catalytic activity data for selected probe chemistries, with relevant material characterization and reactor configuration information, the database provides a collection of catalytic benchmarks for distinct classes of active site functionality. Through key choices in data access, availability, and traceability, CatTestHub seeks to balance the fundamental information needs of chemical catalysis and the FAIR data design principles. Details of the database architecture and the means through which to navigate it are presented, highlighting examples of catalytic insights readily drawn from the available benchmarking data. In its current iteration, CatTestHub spans over 250 unique experimental data points, collected over 24 solid catalysts, that facilitated the turnover of 3 distinct catalytic chemistries. Here, a roadmap is presented through which to expand the open-access platform that serves as a community wide benchmark, primarily through continuous addition of kinetic information on select catalytic systems by members of the heterogeneous catalysis community at large.

Benchmark

Tracing Phase Transformation and Lattice Evolution in a TRIP Sheet Steel under High-Temperature Annealing by Real-Time In Situ Neutron Diffraction

Real-time in situ neutron diffraction was used to characterize the crystal structure evolution in a transformation-induced plasticity (TRIP) sheet steel during annealing up to 1000 °C and then cooling to 60 °C. Based on the results of full-pattern Rietveld refinement, critical temperature regions were determined in which the transformations of retained austenite to ferrite and ferrite to high-temperature austenite during heating and the transformation of austenite to ferrite during cooling occurred, respectively. The phase-specific lattice variation with temperature was further analyzed to comprehensively understand the role of carbon diffusion in accordance with phase transformation, which also shed light on the determination of internal stress in retained austenite. These results prove the technique of real-time in situ neutron diffraction as a powerful tool for heat treatment design of novel metallic materials.

Yu, Dunji [ORNL] (ORCID:0000000189467851)

Observation of Extraordinary Vibration Scatterings Induced by Strong Anharmonicity in Lead‐Free Halide Double Perovskites

Abstract Lead‐free halide double perovskites provide a promising solution for the long‐standing issues of lead‐containing halide perovskites, i.e., the toxicity of Pb and the low stability under ambient conditions and high‐intensity illumination. Their light‐to‐electricity or thermal‐to‐electricity conversion is strongly determined by the dynamics of the corresponding lattice vibrations. Here, the measurement of lattice dynamics is presented in a prototypical lead‐free halide double perovskite(Cs 2 NaInCl 6 ). The quantitative measurements and first‐principles calculations show that the scatterings among lattice vibrations at room temperature are at the timescale of ≈1 ps, which stems from the extraordinarily strong anharmonicity in Cs 2 NaInCl 6 . Further the degree of anharmonicity of each type of atom is quantitatively characterized in the Cs 2 NaInCl 6 single crystal, which stems from the interatomic forces, and demonstrate that this strong anharmonicity is synergistically contributed by the bond hierarchy, the tilting of the NaCl 6 and InCl 6 octahedral units, and the rattling of Cs + ions. Consequently, the crystalline Cs 2 NaInCl 6 possesses an ultralow thermal conductivity of ≈0.43 W mK −1 at room temperature, and a weak temperature dependence ofT −0.41 . These findings uncovered the underlying mechanisms behind the dynamics of lattice vibrations in double perovskites, which can largely benefit the design of optoelectronics and thermoelectrics based on halide double perovskites.

Chemistry

Advances in Design Capabilities for Planetary Missions from the NASA Entry Systems Modeling and Instrumentation Portfolio

The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.

Justin Haskins

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO

Influence of gradation on failure mode of saturated sand at particle scale

Extensive research has been reported in the literature to characterize the failure mode of dry sand using various experimental techniques such as surface optical imaging, photo-elastic materials, three-dimensional (3D) computed tomography (CT), and 3D synchrotron micro-computed tomography (SMT). However, there is a limited literature about the behavior of saturated sand. This paper presents the results of axisymmetric triaxial compression (ATC) experiments that were conducted on saturated sand specimens. The behavior of specimens composed of a uniform sand with grain size between US sieves #40 and #50 is compared to specimens conducted on the same sand that has a wider gradation. 3D SMT technique was used to acquire 3D scans while shearing the specimens to probe localized events that are completely missed or misinterpreted when analyzing ATC measurements based on global standard measurements. The results show a higher effective principal stress ratio (EPSR) for the non-uniform specimen and a thicker shear band when compared to uniform specimen.

Elnur, Mohammed [University of Tennessee]

Uncovering Hidden Entanglement in Twin Beams

Proper characterization of quantum correlations in multimode optical quantum states is critical for applications in quantum information science. However, standard entanglement measurements can lead to incomplete state reconstruction and characterization. Here, we implement a resonator-based detection system that reveals entanglement between sideband modes of twin beams, achieving full tomography and retrieving often ignored quantum correlations. Unlike standard spectral measurements such as homodyne detection, resonator detection can independently address the sidebands of each beam, thereby accessing these hidden correlations. Additionally, we show how phase shifts between the carrier and the sideband modes of the involved fields redistribute information and modify the observation of entanglement for different witnesses. The ability of the resonant detection to independently address sideband modes of entangled states can contribute to enhancing the capacity for secure communication and quantum networking protocols.

Rincon Celis, Raul [University of Sao Paulo, Brazi

Well-Defined Iron Sites in Crystalline Carbon Nitride

Carbon nitride materials can be hosts for transition metal sites, but Mössbauer studies on iron complexes in carbon nitrides have always shown a mixture of environments and oxidation states. Here we describe the synthesis and characterization of a crystalline carbon nitride with stoichiometric iron sites that all have the same environment. The material (formula C 6 N 9 H 2 Fe 0.4 Li 1.2 Cl, abbreviated PTI/FeCl 2 ) is derived from reacting poly(triazine imide)·LiCl (PTI/LiCl) with a low-melting FeCl 2 /KCl flux, followed by anaerobic rinsing with methanol. X-ray diffraction, X-ray absorption and Mössbauer spectroscopies, and SQUID magnetometry indicate that there are tetrahedral high-spin iron(II) sites throughout the material, all having the same geometry. As a result, the material is active for electrocatalytic nitrate reduction to ammonia, with a production rate of ca. 0.1 mmol cm –2 h –1 and Faradaic efficiency of ca. 80% at −0.80 V vs RHE.

Anions

Modeling Study of Hatch Spacing’s Effect on Grain Morphology in Repairing Damaged Metal Parts with Welding and Hatch Spacing’s Potential Use in Lunar Exploration

The welding process is a potential way of repairing a damaged metal component, especially cavity damage caused by a harsh environment like the lunar environment, which is characterized by large temperature differences and reduced gravity. The adjustment of welding parameter (e.g., hatch spacing) can improve production efficiency in the repair process. Seen from the microstructural level, hatch spacing sensitivity affects the metallic grain evolution and morphology in the welding process, which can further influence a repaired part’s mechanical properties; however, the study of hatch spacing’s effect on microstructure is challenging. Traditional experimental procedures are costly and time-consuming, and any change in hatch spacing value needs roll-back of experimental procedure. A modeling study can address the above challenges in experimental observation. In this research, a modeling approach based on the Kinetic Monte Carlo (KMC) Potts theory was used to simulate grain evolution and morphology with three hatch spacings. Through quantifying and analyzing the predicted grain morphologies, the effect of hatch spacing on microstructure in a welding-fabricated part was investigated. The predicted grain morphologies were validated with an EBSD image of welding microstructure, which has been published before. The primary grain morphologies were columnar grains with a small amount of fine equiaxed grains formed in the scanning path centerline. When increasing the hatch spacing, the columnar grains become larger and more lengthy, while the effect of hatch spacing on the equiaxed grains is not obvious.

Welding for repairing

Properties and Potential of Two (Ni,Pt)Ti Alloys for Use as High-Temperature Actuator Materials

The microstructure, transformation temperatures, basic tensile properties, shape memory behavior, and work output for two (Ni,Ti)Pt high-temperature shape memory alloys have been characterized. One was a Ni 30 Pt 20 Ti 50 alloy (referred to as 20Pt) with transformation temperatures above 230 °C and the other was a Ni 20 Pt 30 Ti 50 alloy (30Pt) with transformation temperatures above 530 °C. Both materials displayed shape memory behavior and were capable of 100% (no-load) strain recovery for strain levels up to their fracture limit (3-4%) when deformed at room temperature. For the 20Pt alloy, the tensile strength, modulus, and ductility dramatically increased when the material was tested just above the austenite finish (A f ) temperature. For the 30Pt alloy, a similar change in yield behavior at temperatures above the A f was not observed. In this case the strength of the austenite phase was at best comparable and generally much weaker than the martensite phase. A ductility minimum was also observed just below the A s temperature in this alloy. As a result of these differences in tensile behavior, the two alloys performed completely different when thermally cycled under constant load. The 20Pt alloy behaved similar to conventional binary NiTi alloys with work output due to the martensite-to-austenite transformation initially increasing with applied stress. The maximum work output measured in the 20Pt alloy was nearly 9 J/cm 3 and was limited by the tensile ductility of the material. In contrast, the martensite-to-austenite transformation in the 30Pt alloy was not capable of performing work against any bias load. The reason for this behavior was traced back to its basic mechanical properties, where the yield strength of the austenite phase was similar to or lower than that of the martensite phase, depending on temperature. Hence, the recovery or transformation strain for the 30Pt alloy under load was essentially zero, resulting in zero work output.

High-Temperature Shape Memory Alloy

The Growth of Zeolites A, X and Mordenite in Space

Zeolites are a class of crystalline aluminosilicate materials that form the backbone of the chemical process industry worldwide. They are used primarily as adsorbents and catalysts and support to a significant extent the positive balance of trade realized by the chemical industry in the United States (around $19 billion in 1991). The magnitude of their efforts can be appreciated when one realizes that since their introduction as 'cracking catalysts' in the early 1960's, they have saved the equivalent of 60 percent of the total oil production from Alaska's North Slope. Thus the performance of zeolite catalysts can have a profound effect on the U.S. economy. It is estimated that a 1 percent increase in yield of the gasoline fraction per barrel of oil would represent a savings of 22 million barrels of crude oil per year, representing a reduction of $400 million in the United States' balance of payments. Thus any activity that results in improvement in zeolite catalyst performance is of significant scientific and industrial interest. In addition, due to their 'stability,' uniformity, and, within limits, their 'engineerable' structures, zeolites are being tested as potential adsorbents to purify gases and liquids at the parts-per-billion levels needed in today's electronic, biomedical, and biotechnology industries and for the environment. Other exotic applications, such as host materials for quantum-confined semiconductor atomic arrays, are also being investigated. Because of the importance of this class of material, extensive efforts have been made to characterize their structures and to understand their nucleation and growth mechanisms, so as to be able to custom-make zeolites for a desired application. To date, both the nucleation mechanics and chemistry (such as what are the 'key' nutrients) are, as yet, still unknown for many, if not all, systems. The problem is compounded because there is usually a 'gel' phase present that is assumed to control the degree of supersaturation, and this gel undergoes a continuous 'polymerization' type reaction during nucleation and growth. Generally, for structure characterization and diffusion studies, which are useful in evaluating zeolites for improving yield in petroleum refining as well as for many of the proposed new applications (e.g., catalytic membranes, molecular electronics, chemical sensors) large zeolites (greater than 100 to 1000 times normal size) with minimum lattice defects are desired. Presently, the lack of understanding of zeolite nucleation and growth precludes the custom design of zeolites for these or other uses. It was hypothesized that the microgravity levels achieved in an orbiting spacecraft could help to isolate the possible effects of natural convection (which affects defect formation) and minimize sedimentation, which occurs since zeolites are twice as dense as the solution from which they are formed. This was expected to promote larger crystals by allowing growing crystals a longer residence time in a high-concentration nutrient field. Thus it was hypothesized that the microgravity environment of Earth orbit would allow the growth of large, more defect-free zeolite crystals in high yield.

A Sacco, Jr

3D Woven Mid-Density Carbon Phenolic (3MDCP) Full-Scale Heatshield Development for Mars Sample Return (MSR) Earth Entry System (EES)

The Mars Sample Return (MSR) mission will be returning samples of Martian soil to Earth. The samples will be returned in the Earth Entry System (EES). The EES capsule is protected by a new thermal protection system; the 3-Dimensionally Woven, Mid-Density, Carbon Phenolic (3MDCP) material which is derived from the Heatshield for Extreme Entry Environment Technology (HEEET). The baseline 3MDCP design is a single piece thermal protection system that avoids the manufacturing and certification challenges associated with a tiled configuration. A 3MDCP heatshield begins as a flat woven preform that is formed to the final heatshield shape and then infused with phenolic resin. A NASA Ames lead team, in collaboration with Spirit Textiles (formerly TEAM Inc.), Fiber Materials Inc. (FMI, a Spirit AeroSystems Company), and Kratos SRE, has been developing and demonstrating the processes to fabricate the 3MDCP heatshield and characterize the resulting properties. The process of forming the flat woven preform into the final heatshield shape, a sphere-cone geometry with a 52.5° cone angle, involves local but substantial movement of the yarns in the weave. This results in continuous fibers across the single piece heatshield, albeit with property variations between different regions of the heatshield. This presentation will provide a high-level status of 3MDCP development for the MSR EES heatshield. This will include an overview of the manufacturing processes with an emphasis on the impact of forming on fiber orientation, material properties, and performance. Results comparing preliminary formed and flat materials will be shown. Additionally, the presentation will layout the broader plan for property testing and tie-in to aerothermal performance.

Peter Edward Marshall

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Helmholtz Motor: A Novel Electric Machine That Enables High Temperature Superconducting Armatures

Electrified aircraft are being developed to increase the efficiency and reduce the cost of operating subsonic transport aircraft. Achieving a substantial impact necessitates focusing on single- and twin-aisle aircraft which use propulsion systems with >20 MW ratings. For these aircraft, multi-MW superconducting electric machines are being developed due to their high specific power and high efficiency. Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Second generation high temperature superconductors (HTS) are the ideal conductor for the field winding (typically the rotor) of superconducting machines, and they have attractive features for the armature winding (typically the stator) of these machines. However, their use in armature windings has been severely limited due to excessive AC losses in conventional electric machines. This presentation will present the working principle, design, and a performance trade study of a novel electric machine – termed a Helmholtz machine – that enables HTS armatures by significantly reducing AC losses. The working principle of the patent pending Helmholtz motor is to optimize the rotor-produced magnetic field (i.e., the field winding) so that it produces a magnetic field almost purely in the plane of the HTS. This is accomplished using 2 rotors that contain a set of matched Helmholtz coil pairs (or matched permanent magnet Halbach arrays). The armature (stator) HTS coils are positioned between the rotors and oriented to minimize the out of plane magnetic field. The armature HTS coils are also designed to minimize the out of plane component of their self field. The presentation will include the lessons learned in the sizing and design of this motor, including equations for how to match the magnet arrays. A performance trade study will be presented for a partially superconducting, radial flux version of the motor. The trade study involved an optimization of a motor design code based on analytical calculations (electromagnetic, AC loss, thermofluid, mechanical, and motor sizing) that determine the motor’s total efficiency and total specific power. The design code will be described in the presentation, with an emphasis on the calculation of magnetic fields, AC loss, and temperature distribution.

Superconducting electric machines

Crystal Growth Scale-Up and Stability of RbSr 2 X 5 :Eu Scintillators

The discovery and development of new scintillation materials support national security needs. In these applications, large volumes of scintillator crystals are needed to achieve efficient screening for contraband. Therefore, one important step in the discovery of new scintillator compositions is testing their feasibility for scale-up and their stability. In this work, high-quality Ø22 mm crystals of two new scintillators RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu were grown via the Vertical Bridgman method, and their scintillation properties were characterized. Here, the Ø22 mm RbSr 2 I 5 :Eu crystals could be grown with fast translation rates up to 3.5 mm/h. Both RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu had high scintillation performance, including light yields of 46,000 and 61,000 ph/MeV, respectively, for Ø22 × 35 mm crystals. Additionally, properties related to physical stability were investigated, including the coefficients of thermal expansion via high-temperature X-ray diffraction (HTXRD) as well as moisture sensitivity. HTXRD confirmed the absence of solid–solid phase transitions and showed that RbSr 2 Br 5 had minimal thermal expansion anisotropy compared to RbSr 2 I 5 and some other inorganic metal halide scintillators, which favors the growth of large-sized crystals.

36 MATERIALS SCIENCE

Surface deposition of nanometer scale carbon structures on Bi 2 Sr 2 CaCu 2 O 8+δ using a low-cost scanning electron microscope

In this work, we investigate optically active nanostructures on Bi 2 Sr 2 CaCu 2 O 8+δ (BSCCO). The nanostructures are constructed from carbon nanocrystals formed using the electron beam of a scanning electron microscope to locally decompose residual hydrocarbons at the BSCCO surface. Atomic force microscopy is used to characterize the dimensions of these structures as a function of beam exposure time. This technique has been used to induce localized intercalation to form carbon nanoparticles up to tens of nanometers into a variety of 2D materials. Cross-sectional scanning transmission electron microscopy in BSCCO shows the presence of these carbon deposits on the surface, but there is no evidence they penetrate into the BSCCO substrate. Therefore, this technique is not suitable for localized intercalation in BSCCO, and it does have promise as a cost-effective approach to pattern nanometer scale etch stops.

2D materials

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation