Search NASASearch

SEARCH · Search NASA

Results for “superconducting properties and materials”

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.

3,255 records · Page 14

Stellar Population Properties in the Stellar Streams around SPRC047

Abstract We have investigated the properties (e.g., age, metallicity) of the stellar populations of a ringlike tidal stellar stream (or streams) around the edge-on galaxy SPRC047 (z= 0.031) using spectral energy distribution (SED) fits to integrated broadband aperture flux densities. We used visual images in six different bands and Spitzer/IRAC 3.6μm data. We have attempted to derive best-fit stellar population parameters (metallicity, age) in three noncontiguous segments of the stream. Due to the very low surface brightness of the stream, we have performed a deconvolution with a Richardson–Lucy–type algorithm of the low spatial resolution 3.6μm IRAC image, thereby reducing the effect of the point-spread function aliasedemissionfrom the bright edge-on central galaxy at the locations of our three stream segments. Our SED fits that used several different star formation (SF) history priors, from an exponentially decaying SF burst to continuous SF, indicate that the age–metallicity–dust degeneracy is not resolved, most likely because of inadequate wavelength coverage and low signal-to-noise ratios of the low surface brightness features. We also discuss how future deep visual–near-infrared observations, combined with absolute flux calibration uncertainties at or below the 1% level, complemented by equally well absolute flux-calibrated observations in ultraviolet and mid-infrared bands, would improve the accuracy of broadband SED fitting results for low surface brightness targets, such as stellar streams around nearby galaxies that are not resolved into stars.

Astronomy & Astrophysics

Development of large-scale, 3D Printed high temperature ceramic material

Through this collaborative effort, a new 3D printing platform for ceramics called laser-induced slip casting (LIS) was explored and developed for improving the processing and manufacture of silicon carbide (SiC), a high temperature ceramic material. This method prints layers of ceramic slips/slurries with subsequent selective-laser heating to dry each layer to build a 3D structure. The completed and dried part is then sintered. This manufacturing process is inherently lower cost than alternative ceramic printing methods such as binder jet technology or stereolithography when it comes to the feedstock material but is more expensive than robocasting or direct ink writing. However, it has the potential to make more controlled parts with less defects compared to robocasting. The largest cost is the heating source for the printer. With this technique, there is potential to make large ceramic parts in a near-net shape. Further, there is no known commercial manufacturing of 3D printed ceramics that creates a large range of different high-density ceramics at large scale. The goals and outcomes for the development of the new printing technology were: 1) assessing the technology with alumina by characterizing coupons, 2) printing and sintering of silicon carbide (SiC), and 3) characterization and properties testing of materials printed and a scale-up of SiC part(s). Through this collaborative project, it was sought to provide the best solution to achieving highly dense and near-net shaped ceramic parts at large scale and reduced cost. The goal was to produce high density sintered materials for applications in defense, energy generating systems, armor, and wear parts using 3D printing methods. Variables including powder particle sizes, dispersant molecular weights (MWs), binders, printing parameters, and post processing were explored as well as the characterization of the physical properties (add specifics here).

36 MATERIALS SCIENCE

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence

Phosphorylation toggles the SARS-CoV-2 nucleocapsid protein between two membrane-associated condensate states

Abstract The Nucleocapsid protein (N) of SARS-CoV-2 plays a critical role in the viral lifecycle by regulating RNA replication and by packaging the viral genome. N and RNA phase separate to form condensates that may be important for these functions. Both functions occur at membrane surfaces, but how N toggles between these two membrane-associated functional states is unclear. Here, we reveal that phosphorylation switches how N condensates interact with membranes, in part by modulating condensate material properties. Our studies also show that phosphorylation alters N’s interaction with viral membrane proteins. We gain mechanistic insight through structural analysis and molecular simulations, which suggest phosphorylation induces a conformational change in N that softens condensate material properties. Together, our findings identify membrane association as a key feature of N condensates and provide mechanistic insights into the regulatory role of phosphorylation. Understanding this mechanism suggests potential therapeutic targets for COVID infection.

Science & Technology - Other Topics

Catalytic Borylation of Poly(vinyl chloride) Produces Adhesive Materials

Postpolymerization functionalization of polymers can create new applications for existing materials, while retaining their most favorable, intrinsic properties. Polyvinyl chloride (PVC) is a widely used, commodity polymer that is particularly challenging to modify. We report a copper-catalyzed protocol that replaces a small fraction of the C-Cl bonds in PVC with C-B bonds to boronic esters. The reaction occurs with an inexpensive catalyst comprising copper(II) chloride and an NHC ligand derived from a common ionic liquid that is distinct from the N-heterocyclic carbene (NHC) used for the borylation of small alkyl halides. The resulting materials adhere strongly to common surfaces, such as glass and metals, even more strongly than do commercial glues.

D’Angelo, Kyan A

Predicting Band-Gap of Inorganic Materials Using Neuromorphic Graph Learning

Predicting properties of inorganic materials is a heavily researched topic, with several new prediction approaches emerging as competitors. One such competitor is graph neural networks, which leverage the structure of the graph to aid in the prediction process. In this work, we propose integration of neuromorphic computation into the graph neural network pipeline. We call this approach Neuromorphic Graph Learning (NGL). We utilize the NGL approach to leverage evolutionary algorithms and a novel Spike Pipeline for Raster Analysis (SPIRE) for the prediction of band gap in inorganic materials.

Mulet, Ian [University of Tennessee (UT)]

Thermoplastic Matrix Composite Design for Cryotanks Using Multiscale Modeling and Bayesian Optimization

Designing lightweight, robust cryogenic storage tanks is critical for future launch vehicles, in-space propellant storage, and hydrogen powered aircraft. This work presents a multiscale modeling and Bayesian optimization framework for the design of thermoplastic matrix composite cryotanks. Molecular dynamics simulations are first used to determine temperature-dependent constituent properties for candidate thermoplastic matrices, which are homogenized to the lamina scale using NASA’s Multiscale Analysis Tool (NASMAT). These lamina properties, in combination with laminate family generation rules, are evaluated in HyperX structural optimization software to identify stacking sequences that meet all cryogenic load requirements. A Bayesian optimization framework is applied, with HyperX in the loop (via the HyperX API) to efficiently search across material and laminate design variables, yielding an optimized cryotank configuration with significant reductions in design cycle time compared to exhaustive search approaches.

thermoplastics

Mechanisms of Alkali Ionic Transport in Amorphous Oxyhalides Solid State Conductors

Amorphous oxyhalides have attracted significant attention due to their relatively high ionic conductivity (1 mS cm –1 ), excellent chemical stability, mechanical softness, and facile synthesis routes via standard solid‐state reactions. These materials exhibit an ionic conductivity that is almost independent of the underlying chemistry, in stark contrast to what occurs in crystalline conductors. In this work, we employ machine learning interatomic potentials to construct large‐scale molecular dynamics trajectories encompassing hundreds of nanoseconds to obtain statistically converged transport properties. We find that the amorphous state consists of chain fragments of metal‐anion tetrahedra of various lengths. By analyzing the residence time of alkali cations migrating around tetrahedrally‐coordinated metals, we find that oxygen anions limit alkali diffusion. By computing the full Einstein expression of the ionic conductivity, we demonstrate that the alkali transference number of these materials is strongly influenced by distinct‐particles correlations, while alkali transport is dictated by uncorrelated self‐diffusion. By extending this analysis to chemical compositions AMX 2.5 O 0.75 , spanning different alkaline (A = Li, Na, K), metallic (M = Al, Ga, In), and halogen (X = Cl, Br, I) species, we clarify why the diffusion properties of these materials remain largely insensitive to variations in atomic isovalent chemistry.

amorphous materials

Real-space visualization of a defect-mediated charge density wave transition

Here, we study the coupled charge density wave (CDW) and insulator-to-metal transitions in the 2D quantum material 1T-TaS 2 . By applying in situ cryogenic 4D scanning transmission electron microscopy with in situ electrical resistance measurements, we directly visualize the CDW transition and establish that the transition is mediated by basal dislocations (stacking solitons). We find that dislocations can both nucleate and pin the transition and locally alter the transition temperature T c by nearly ~75 K. This finding was enabled by the application of unsupervised machine learning to cluster five-dimensional, terabyte scale datasets, which demonstrate a one-to-one correlation between resistance—a global property—and local CDW domain-dislocation dynamics, thereby linking the material microstructure to device properties. This work represents a major step toward defect-engineering of quantum materials, which will become increasingly important as we aim to utilize such materials in real devices.

4D-STEM

Electron-Deficient Phenazines: Synthesis, Structures, Redox, and Optoelectronic Properties in the Gas Phase, in Solution, and in the Solid State

High-temperature gas-phase reactions of 1,4-C4F8I2 with phenazine (PHNZ) afforded new thermally- and air-stable derivatives with cyclo-C4F8 substituents, viz. 1,2-PHNZ(C4F8), 2,3-PHNZ(C4F8), 1,2:6,7-PHNZ(C4F8)2, 1,2:7,8-PHNZ(C4F8)2, and 1,2:8,9-PHNZ(C4F8)2. Reactions in the presence of Cu powder resulted in reductive-defluorination/aromatization of some of the cyclo-C4F8 substituents to produce 2,3-PHNZ(C4F4), 2,3-PHNZ(C4HF3), and 1,2:6,7-PHNZ(C4F8)(C4F4), which are formally tri- or tetrafluoro derivatives of benzo[a]phenazine. Single-crystal X-ray diffraction structures of these compounds demonstrate the planarity of 2,3-PHNZ(C4F4) and 2,3-PHNZ(C4HF3) and ordered p-stacking in 1,2-PHNZ(C4HF3), 1,2-PHNZ(C4F4), 1,2-PHNZ(C4F8), and 2,3-PHNZ(C4F8). Low-temperature gas-phase photoelectron spectra of the PHNZ(C4F8)1,2 compounds show that they are strong electron acceptors, with the three PHNZ(C4F8)2 isomers having electron affinities exceeding that of C60. Cyclic voltammetry in acetonitrile show that the five PHNZ(C4F8)1,2 compounds exhibit reversible one-electron reductions with large positive shifts in their reduction potentials relative to unsubstituted PHNZ. Spectroelectrochemical absorption and EPR spectra of PHNZ(C4F8)1,2- • anion radicals and extensive DFT calculations revealed the significant role that the different substitution patterns have on determining the optoelectronic properties of the PHNZ(C4F8)1,2 derivatives and PHNZ(C4F8)1,2- • anion radicals. Time-resolved microwave conductivity measurements and photoluminescence spectra of blends of the PHNZ(C4F8)1,2 derivatives with the donor poly(3-hexylthiophene) (P3HT) were recorded to assess the possible use of the new compounds in optoelectronic devices.

Balser, Sebastian

Expanding the horizon of bio-naphtha beyond gasoline blend: property characterization and conversion opportunity assessment through technoeconomic and life-cycle analyses

Bio-naphtha, a common by-product of biorefineries, is expected to experience substantial growth in supply due to increasing demands for renewable diesel and synthetic aviation fuel (SAF). However, demand for bio-naphtha itself as a gasoline blendstock is limited because of the electrification of light-duty vehicles. This work investigated valorization opportunities for bio-naphtha from catalytic fast pyrolysis, hydrothermal liquefaction, Fischer–Tropsch synthesis, and hydrotreated esters and fatty acids pathways. These opportunities include producing polymer-grade olefin via steam cracking, SAF via steam cracking followed by olefin oligomerization, and renewable aromatics benzene, toluene, and xylene (BTX) and hydrogen via catalytic reforming. Process models were developed in Aspen Plus V14 and Aspen HYSYS V14 to calculate the mass and energy balances for each conversion step. Technoeconomic assessment and life-cycle analysis were conducted to evaluate the minimum fuel/product selling price, conversion cost, and life-cycle CO2 equivalent (CO2e) emission reduction. Technoeconomic assessment results suggest a minimum fuel/product selling price as low as $1.9/kg of olefins, $6.30/gal of gasoline-equivalent SAF, and $1.2/kg of BTX without any incentives. For all pathways, these prices are dominated by bio-naphtha feedstock costs, which account for at least 76% of the total cost. Compared with petroleum baselines, bio-naphtha-derived SAF, olefins, and BTX can achieve significant CO2e emission reductions from the use of renewable carbon resources. The results of life-cycle analysis and subsequent technoeconomic assessment, incorporating carbon credits, indicate the economic viability of using bio-naphtha for polymer-grade olefin and BTX production, with product costs comparable to market prices.

Xu, Shuang

Lubricating Properties of Ceramic-Bonded Calcium Fluoride Coatings on Nickel-Base Alloys From 75° to 1900° F

The endurance life and the friction coefficient of ceramic-bonded calcium fluoride (CaF2) coatings on nickel-base alloys were determined at temperatures from 75° to 1900° F. The specimen configuration consisted of a hemispherical rider (3/16-in. rad.) sliding against the flat surface of a rotating disk. Increasing the ambient temperature (up to 1500° F) or the sliding velocity generally reduced the friction coefficient and improved coating life. Base-metal selection was critical above 1500° F. For instance, cast Inconel sliding against coated Inconel X was lubricated effectively to 1500° F, but at 1600° F severe blistering of the coatings occurred. However, good lubrication and adherence were obtained for René 41 sliding against coated René 41 at temperatures up to 1900° F; no blisters developed, coating wear life was fairly good, and the rider wear rate was significantly lower than for the unlubricated metals . Friction coefficients were 0.12 at 1500° F, 0.15 at 1700° F, and 0.17 at 1800° and 1900° F. Because of its ready availability, Inconel X appears to be the preferred substrate alloy for applications in which the temperature does not exceed 1500° F. René 41 would have to be used in applications involving higher temperatures. Improved coating life was derived by either preoxidizing the substrate metals prior to the coating application or by applying a very thin (<0.0002 in.) burnished and sintered overlay to the surface of the coating. Preoxidation did not affect the friction coefficient. The overlay generally resulted i n a higher friction coefficient than that obtained without the overlay. The combination of both modifications resulted in longer coating life and in friction coefficients intermediate between those obtained with either modification alone.

Nickel alloy

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen

Evaluating Crystallinity in Thermoplastic Composites for Aerospace Applications

Polymer matrix composites (PMCs) offer many benefits for the aerospace industry due to their potential for weight reduction when compared to metal or ceramic based materials. Most PMCs currently in flight use thermoset matrices, however, thermoplastic resins are being explored as alternatives due to their ability to be remelted, which is of particular interest due to the potential for in-situ repair and faster production. Most thermoplastic resins are semicrystalline polymers. The properties of semicrystalline thermoplastics are directly influenced by their crystallinity, which can vary due to many factors including thermal treatments, environmental conditions, and mechanical deformation. Within thermoplastic PMC parts, crystallinity gradients can arise due to variations in part geometry, across part thicknesses, and along bonded joints. Monitoring the crystallinity of thermoplastic composites is key to ensuring these materials meet the high demands required for aerospace. Several different analytical techniques exist that can be used to characterize the bulk crystallinity of thermoplastic materials. However, many existing methods lack the specificity required to identify the subtle variations in crystallinity that may play a significant role in the performance and durability of PMC parts. Because of this, a significant amount of work is still required to fully characterize and understand the crystallinity profiles of thermoplastic PMCs and the resulting impact to material properties. This talk discusses the use of multiple techniques such as Differential Scanning Calorimetry, Polarized Light Optical Microscopy, and Fourier-Transform Infrared Spectroscopy to characterize the crystallinity in carbon fiber/thermoplastic composites. Samples of different crystallinity profiles were manufactured using various cooling procedures. This work aims to provide the fundamental data necessary to understand the effects of crystallinity on thermoplastic PMCs, which is key to advancing their use in aerospace applications.

Thermoplastics

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

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