Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “tensor”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Two-Year Cosmology Large Angular Scale Surveyor (CLASS) Observations: 40GHz Telescope Pointing, Beam Profile, Window Function, and Polarization Performance

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over 75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale (1°~< θ >=90°) CMB polarization to constrain the tensor-to-scalar ratio at the r ∼ 0.01 level and the optical depth to last scattering to the sample variance limit. This paper presents the optical characterization of the 40 GHz telescope during its first observation era, from 2016 September to 2018 February. High signal-to-noise observations of the Moon establish the pointing and beam calibration. The telescope boresight pointing variation is <0°. 023 (<1.6% of the beam’s full width at half maximum (FWHM)). We estimate beam parameters per detector and in aggregate, as in the CMB survey maps. The aggregate beam has an FWHM of 1°.579 ± 0°.001 and a solid angle of 838 ± 6 μsr, consistent with physical optics simulations. The corresponding beam window function has a sub-percent error per multipole at ℓ < 200. An extended 90° beam map reveals no significant far sidelobes. The observed Moon polarization shows that the instrument polarization angles are consistent with the optical model and that the temperature-to-polarization leakage fraction is <10(exp −4) (95% C.L.). We find that the Moon-based results are consistent with measurements of M42, RCW 38, and Tau A from CLASS’s CMB survey data. In particular, Tau A measurements establish degree- level precision for instrument polarization angles.

Xu, Zhilei↗
NASA NTRS2019

Adjoint-Based Anisotropic Mesh Adaptation for a Stabilized Finite-Element Flow Solver

An adjoint solver is implemented in the FUN3D stabilized finite-element flow solver. The adjoint solution is used to generate anisotropic, adapted meshes to control error in scalar out- put functionals, such as lift or drag coefficient. The inviscid and turbulent adjoints are verified with a finite-difference derivative approximation and can be used in design optimization in addition to mesh adaptation. The adjoint capability represents an extension of previous research using the FUN3D stabilized finite-element capability with metric-based mesh adaptation and interpolation-based error estimates to generate highly anisotropic adapted meshes for turbulent flows. In the present work, a metric-based approach is again utilized, where the adjoint and the primal solutions both contribute to the generation of a metric tensor field that is subsequently used to produce the required anisotropic mesh for each adaptation cycle. Adaptive results are then shown for an inviscid supersonic flow over a diamond airfoil, inviscid transonic flow over the ONERA M6 wing, and viscous laminar flow over the NACA 0012 airfoil, all using drag force as the output functional. Adjoint-based adaptation is compared with a multiscale solution-based approach that controls the L norm of Mach number interpolation error to demonstrate the effectiveness and effciency of the adjoint-based adaptive mesh technology.

Balan, Aravind↗
NASA NTRS2000

Coseismic Excitation of the Earth's Polar Motion

Apart from the shaking near the epicenter that is the earthquake, a seismic event creates a permanent field of dislocation in the entire Earth. This redistribution of mass changes (slightly) the Earth's interia tensor; and the Earth's rotation will change in accordance with the conservation of angular momentum.

Earth↗
NASA NTRS2019

Efficient method for approximating nonlinear dynamics: applications to uncertainty propagation and estimation

High-order Taylor series expansions can be used to model nonlinear dynamics at the cost of integrating a large set of variational equations to obtain high-order state-transition tensors (STTs). This paper presents an innovative technique for approximating the high-order STTs that reduces significantly the computational cost by retaining only the dominant secular terms. We propagate the low-order partial derivatives of Kepler’s equation, which only requires the integration of six additional equations to extend an n-th order approximation to order (n + 1). The approximation stems from the Lindstedt-Poincare procedure and exploits the stability properties of orbital motion. Since the method makes no dynamical assumptions, it can accommodate any source of orbital perturbations. We show how the approximation of the second-order STT significantly increases the accuracy of the linear method for uncertainty propagation with only a small computational overhead. Finally, we derive a high-order approximate extended Kalman filter that implements the proposed approximation of the STT and improves the performance of linear filters. Examples of application with different perturbation sources include the heliocentric orbit of an asteroid, an orbiter around Europa, and an Earth-orbiting satellite.

Park, Ryan S↗
NASA NTRS2020

Small Aperture Telescopes for the Simons Observatory

The Simons Observatory (SO) is an upcoming cosmic microwave background(CMB) experiment located on Cerro Toco, Chile, that will map the microwave sky in temperature and polarization in six frequency bands spanning 27 to 285 GHz. SO will consist of one 6-meter Large Aperture Telescope (LAT) fielding∼30,000 detectors and an array of three 0.42-meter Small Aperture Telescopes (SATs) fielding an additional 30,000 detectors. This synergy will allow for the extremely sensitive characterization of the CMB over an-gular scales ranging from an arcmin to tens of degrees, enabling a wide range of scientific output. Here we focus on the SATs targeting degree angular scales with successive dichroic instruments observing at Mid-Frequency (MF: 93/145 GHz), Ultra-High-Frequency (UHF:225/285 GHz), and Low-Frequency (LF: 27/39 GHz). The three SATs will be able to map∼10% of the sky to a noise level of∼2 μK-arcmin when combining 93 and 145 GHz. The multiple frequency bands will allow the CMB to be separated from galactic foregrounds (primarily synchrotron and dust), with the primary science goal of characterizing the primordial tensor-to-scalar ratio, r, at a target level ofσ(r)≈0.003.

Small aperature telescope↗
NASA NTRSDate not supplied

Multi-Scale Analysis of Effective Mechanical Properties of 3D Woven Composite Materials

This work aims to compute and analyze the mechanical properties of 3D woven composite materials that share the same yarn structure but differ in their matrix porosity and fiber volume fraction. Assuming linear elastic behavior, the software PuMA, developed at NASA, was used to characterize the different composites from their constituents' data averaging their elastic properties in three main directions to obtain their effective orthotropic mechanical properties. This was performed in three steps: firstly, the different phases (yarns/matrix) were analyzed as anisotropic materials separately at the micro-scale; then, the results from these analyses were homogenized and used to obtain the anisotropic elasticity tensor at the meso-scale for the different woven materials; and finally, assuming orthotropic mechanical behavior for the unit cell the mechanical properties were obtained at the macro-scale. The 3D woven composite materials selected for this study are made of phenolic resin with porosity values between 0 and 0.6, and resin infiltrated carbon fiber with volume fractions between 0.5 and 0.9 for the yarns. Results from PuMA were validated with external FEA simulations for the non-porous 3D woven material with a fiber volume fraction of 0.9 in the yarns. Results will be included in the final paper after completing the study detailed in Section II.

Sergio Fraile Izquierdo↗
NASA NTRS2021

Organized kinetic energy backscatter in the hurricane boundary layer from radar measurements

The fluid mechanics of hurricanes strongly depends on boundary layer energetics due to the warm-core nature of the system with peak velocities located at lower levels. One barrier that has inhibited a more complete characterization of energy transfer in the boundary layer is a lack of observations that resolve large, turbulent eddies. In particular, the occurrence and structure of upscale energy transfer (backscatter) in the hurricane boundary layer as well as the effects of backscatter on the vortex intensity are unknown. The analysis presented here of very high-resolution, three-dimensional wind observations from Hurricane Rita (2005) at peak intensity reveals large regions of organized backscatter in the boundary layer associated with coherent, turbulent eddies. Strong forwardscatter is also found next to the backscatter regions due to the interaction between adjacent eddies. Two components of the stress tensor are primarily responsible for this alternating scatter structure, as shown by large correlation coefficients between the fields: the radial–vertical component (τ13) and the azimuthal–vertical component (τ23) with average correlations of 79% and 49 %, respectively. The Leonard, Reynolds and cross-term stress components are also provided. The impact of the sub-filter-scale energy transfer is estimated by computing the kinetic energy budget for the resolved-scale and eddy-scale motions. The results show that the sub-filter-scale energy transfer term is of the same order as the other terms in the eddy-scale budgets, contributing between 16% and 40% to the local time tendency with an average contribution of approximately 30 %. These results indicate that the coherent turbulent eddies can affect the vortex dynamics through wave–wave nonlinear interactions, which can subsequently influence the wave–mean flow interactions. This is the first study to examine the full sub-filter-scale

Sydney Sroka↗
NASA NTRSDate not supplied

Mechanical Erosion Modeling of TPS Materials

This work describes the development of a model that accounts for the additional total surface recession in Thermal Protection Systems materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics capability module was integrated within the Porous material Analysis Toolbox based on OpenFOAM, PATO. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver module computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model of the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.

Stress Analysis↗
NASA NTRSDate not supplied

The Influence of Adaptive Mesh Refinement on the Prediction of Vortex Interactions about a Generic Missile Airframe

The complex interaction of forebody and wing vortices significantly impacts missile aerodynamics. The formation of these vortices involves smooth regions of the geometry or geometric discontinuities like leading edges, trailing edges, tips, and corners. Regions of supersonic flow and complex shock topologies interact with boundary layers and vortices. Smooth-body separation and 3D viscous effects strain current Reynolds-averaged Navier-Stokes (RANS) techniques. The quantification and control of discretization error is critical to obtaining reliable simulation results and often turbulence model assessments are made in the presence of unquantified (and potentially large) discretization errors. Two mesh adaptation schemes are applied to steady RANS simulations. Multiscale unstructured mesh adaptation is applied to control interpolation error estimates of the Mach field, which resolves boundary layers, vortices, and shocks. A dual-mesh approach with overset communication is applied between an expert-crafted near-body unstructured mesh and an adaptive off-body Cartesian mesh refined with Q-criterion scaled by the strain tensor magnitude. A generic missile configuration is examined in a supersonic flow field to show the interaction of mesh adaptation and turbulence model. Turbulence model modifications for rotational correction and a quadratic constitutive relationship show a strong influence on adaptive mesh refinement and predicted rolling moment.

slender body aerodynamics↗
NASA NTRS2021

Preventing Anomalous Torques in Circumbinary Accretion Simulations

Numerical experiments are the primary method of studying the evolution of circumbinary disks due to the strong nonlinearities involved. Many circumbinary simulations also require the use of numerical mass sinks: source terms that prevent gas from unphysically accumulating around the simulated point masses by removing gas at a given rate. However, special care must be taken when drawing physical conclusions from such simulations to ensure that results are not biased by numerical artifacts. We demonstrate how improved sink methods reduce some of these biases, using simulations of aspect ratio 0.1 accretion disks around binaries with mass ratios between 0.1 and 1. We show that sink terms that do not reduce the angular momentum of gas relative to the accreting object (1) reduce the dependence on the sink rate of quantities such as the torque on the binary, distribution of accretion between binary components, and evolution of the binary semimajor axis; (2) reduce the degree to which the sink rate affects the structure of the accretion disks around each binary component; (3) alter the variability of accretion onto the binary, regularizing it in time. We also investigate other potential sources of systematic error, such as gravitational softening and simplifications to the viscous stress tensor. Because of the strong dependence of binary orbital evolution on both the torque and the distribution of mass between binary components, the sink method influences the orbital evolution of the binary at all mass ratios, with mass ratios below ∼0.3 most affected.

Alexander Dittmann↗
NASA NTRS2021

Cluster after 20 Years of Operations: Science Highlights and Technical Challenges

The Cluster mission was the first constellation using four identical spacecraft to study Sun-Earth connection plasma processes. Using four spacecraft in a tetrahedron shape, it could measure, for the first time, 3D quantities such as electrical currents, plasma gradients or divergence of the electron pressure tensor and 3D structures such as boundaries, surface waves or vortices. Launched in pairs in July and August 2000, on two Soyuz rockets from Baikonur, the four spacecraft have been collecting data continuously for more than 20 years. The mission faced many challenges during the years of operations as some spacecraft subsystems had a lifetime of a few years beyond the initial two-year mission. The major one was to operate without functioning batteries and to successfully pass short and long eclipses, up to 3 h long, without damaging the on-board computers and transmitters and without freezing the fuel. More than 1,000 eclipses have been successfully passed since 2010 using a specially made procedure which switches off the complete spacecraft before entering into eclipse and switches it on when the Sun is again illuminating the solar panels. During 20 years, many discoveries and science results have been published in more than 2,700 scientific papers. A few highlights are presented here, focusing on how varying the spacecraft separation was essential to achieve the science goals of the mission. The Cluster Science Data System and the Cluster archive allows public access to all science data as well as spacecraft ancillary data.

C P Escoubet↗
NASA NTRS2022

Simons Observatory: Constraining inflationary gravitational waves with multitracer B-mode delensing

We introduce and validate a delensing framework for the Simons Observatory (SO), which will be used to improve constraints on inflationary gravitational waves by reducing the lensing noise in measurements of the B modes in CMB polarization. SO will initially observe CMB by using three small aperture telescopes and one large-aperture telescope. While polarization maps from small-aperture telescopes will be used to constrain inflationary gravitational waves, the internal CMB lensing maps used to delens will be reconstructed from data from the large-aperture telescope. Since lensing maps obtained from the SO data will be noise dominated on subdegree scales, the SO lensing framework constructs a template for lensing-induced B modes by combining internal CMB lensing maps with maps of the cosmic infrared background from Planck as well as galaxy density maps from the LSST survey. We construct a likelihood for constraining the tensor-to-scalar ratio r that contains auto and cross spectra between observed B modes and the lensing B-mode template. We test our delensing analysis pipeline on map-based simulations containing survey nonidealities, but that, for this initial exploration, does not include contamination from Galactic and extragalactic foregrounds. We find that the SO survey masking and inhomogeneous and atmospheric noise have very little impact on the delensing performance, and the r constraint becomes σ(r)≈0.0015 which is close to that obtained from the idealized forecasts in the absence of the Galactic foreground and is nearly a factor of 2 tighter than without delensing. We also find that uncertainties in the external large-scale structure tracers used in our multitracer delensing pipeline lead to bias much smaller than the 1σ statistical uncertainties.

Toshiya Namikawa↗
NASA NTRS2020

TopiQAL: Topic-aware Question Answering using Scalable Domain-specific Supercomputers

We all have questions. About today's temperature, scores of our favorite baseball team, the Universe, and about vaccine for COVID-19. Life, physical, and natural scientists have been trying to find answers to various topics using scientific methods and experiments, while computer scientists have built language models as a tiny step towards automatically answering all of these questions across domains given a little bit of context. In this paper, we propose an architecture using state-of-the-art Natural Language Processing language models namely Topic Models and Bidirectional Encoder Representations from Transformers (BERT) that can transparently and automatically retrieve articles of relevance to questions across domains, and fetch answers to topical questions related to COVID-19 current and historical medical research literature. We demonstrate the benefits of using domain-specific supercomputers like Tensor Processing Units (TPUs), residing on cloud-based infrastructure, using which we could achieve significant gains in training and inference times, also with very minimal cost.

Penberthy, Scott↗
NASA NTRS2021

On Parametric Generation of Deep UV Light in Monolithic Cavities With Cubic Nonlinearity

Coherent ultra violet (UV) light has broad applications in optical atomic clocks and quantum computers. Since availability of semiconductor tunable monochromatic lasers is limited at short wavelengths, optical nonlinear processes are frequently used to produce the desirable radiation wavelength. Optical crystals with quadratic nonlinearity are promising for this purpose because of their transparency and high optical damage threshold. In this paper we discuss a possibility of generation of coherent deep UV light utilizing high quality factor (-Q) monolithic optical microcavities made of crystalline optical materials characterized with cubic nonlinearities. We show that the generation is feasible in spite of the relatively small cubic nonlinearities, due to the ultra-broad-band high optical transparency of the crystals and associated high-Q achievable in the cavities created out of the materials. We found that the birefringence as well as the tensor nature of the nonlinear susceptibility of the crystalline optical materials can simplify the phase matching of the nonlinear processes.

Matsko, A.B.↗
NASA NTRSDate not supplied

Mechanical Erosion Modeling of TPS Materials

This work describes the development of a model that accounts for the additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. A computational solid mechanics module was integrated within the Porous material Analysis Toolbox (PATO) based on OpenFOAM. The mechanical erosion was modeled in three steps: first, the implemented stress analysis solver computes the stress and the displacement fields for orthotropic materials using the wall shear stress tensor as a boundary condition; then, regions on the surface where the stress meets the failure criteria are identified; and last, the mesh is moved accordingly to remove the failed material. The outcome is a model capable to predict the total recession of the material due to surface chemistry and mechanical erosion. Results will be included in the final paper after verifying and completing the study.

Stress Analysis↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗
NASA NTRSDate not supplied

Computation of Effective Mechanical Properties and Mechanical Erosion Modeling of TPS Materials

The goal of this presentation is to provide a general overview of the multi-scale modeling formulation to determine if there is additional surface recession in Thermal Protection Systems (TPS) materials as a result of mechanical erosion due to high shear conditions during atmospheric entry. This modeling process is performed at different scales by leveraging two computational frameworks developed at NASA: the Porous Microstructure Analysis (PuMA) software, and the Porous material Analysis Toolbox based on OpenFOAM (PATO). PuMA specializes in computing effective macro-scale material properties by performing material response simulations on 3D digital micro-scale representations of porous micro-structures. The modeling of TPS materials at the micro-scale is essential to understand how they behave as part of a heat shield assembly. The first part of the presentation will detail the implementation of PuMA’s cell-centered finite volume elasticity solver, which allows the computation of macro-scale effective mechanical properties of heterogeneous and anisotropic materials such as fibrous and woven TPS composites. These homogenized mechanical properties are used by PATO’s mechanical erosion model to predict the TPS material’s recession at a larger scale. This work will also provide some examples of multi-scale analysis from the fiber level up to the unit cell. The second part of the presentation will focus on the macro-scale approach to determine if erosion at the heat shield’s surface occurs due to mechanical and thermal loads experienced during atmospheric entry. To accomplish this, a solid mechanics module was integrated within PATO enabling it to model the potential mechanical erosion in three steps: first, after obtaining the effective mechanical properties with PuMA, the implemented stress analysis solver computes the stress and the displacement fields for the TPS material using the wall shear stress tensor, computed using a CFD solver, as boundary conditions; then, regions on the surface where the stress meets the failure criteria are identified; finally, the failed material is removed and the mesh is redistributed accordingly. The outcome is a model capable of predicting the total recession in the material due to surface chemistry and mechanical erosion.

Mechanical Properties↗
NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗