Search NASA⌕ Search

SEARCH · Search NASA

Results for “permeability models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Under Pressure: The Artemis I Heatshield Char Loss Investigation and Artemis II’s Successful Entry

Following the successful Artemis I skip reentry on December 11, 2022, unexpected char liberation from the Orion heat- shield prompted formation of an Anomaly Response Team to determine root cause and establish flight rationale and corrective actions for subsequent missions. Through a comprehensive investigation including detailed hardware analysis, modeling & simulation, and ground testing, the team determined that the low permeability Avcoat experienced extreme internal gas pressure buildup during the skip entry that could not adequately outgas, leading to crack formation and char liberation. Key discoveries included significant performance differences between permeable and impermeable heatshield regions and successful replication of the anomaly through ground testing. Based on extensive ground testing and analysis, Artemis II flew a modified trajectory without the skip entry to minimize crack-inducing conditions, and successfully splashed down on April 10, 2026 with significantly reduced char loss. For Artemis III and beyond, a corrective action was implemented to produce a more permeable version of Avcoat that will substantially reduce internal pressure accumulation during entry.

Artemis II↗

Theory of global thermoremanent magnetization of planetary lithospheres in dipole fields of internal origin

A model is presented for the global thermoremanent magnetization of spherical lithospheres which cool in the presence of central dipole fields. Reversals and intensity variations of the field are incorporated in this model, which is applicable to bodies whose interiors have remained above the Curie point throughout their evolution. The model demonstrates that even considering Runcorn's (1975) magnetostatics theorems for spherical shells, a nonzero magnetic permeability and a finite cooling rate in the lithosphere permit the acquisition of a sizable global remanent dipole moment, which would be detectable by external measurements after the magnetizing field has disappeared. Preliminary application of this model to Mercury, Venus, and Mars suggests that only the combination of a nonreversing ancient source field with a surface value near 1 Oe plus a sizable concentration (about 1% by volume) of ferromagnetic material in their crusts could produce remanent planetary dipole fields as large as those measured by spacecraft. On the other hand, if ancient reversing dynamos existed in these planets, it is unlikely that large planetary-scale fields like those observed at Mercury could be due to remanence in their crusts, irrespective of their composition.

Srnka, L. J.↗

A groundwater convection model for Rio Grande rift geothermal resources

It has been proposed that forced convection, driven by normal groundwater flow through the interconnected basins of the Rio Grande rift is the primary source mechanism for the numerous geothermal anomalies along the rift. A test of this concept using an analytical model indicates that significant forced convection must occur in the basins even if permeabilities are as low as 50-200 millidarcies at a depth of 2 km. Where groundwater flow is constricted at the discharge areas of the basins forced convection can locally increase the gradient to a level where free convection also occurs, generating surface heat flow anomalies 5-15 times background. A compilation of groundwater data for the rift basins shows a strong correlation between constrictions in groundwater flow and hot springs and geothermal anomalies, giving strong circumstantial support to the convection model.

Morgan, P.↗

Iron abundance and magnetic permeability of the moon

A set of simultaneous data from the Apollo 12 lunar surface magnetometer and the Explorer 35 Ames magnetometer are used to construct a whole-moon hysteresis curve, from which a new value of global lunar permeability is determined to be mu = 1.012 + or - 0.006. The corresponding global induced dipole moment is 2.1 x 10 to the 18th power gauss-cucm for typical inducing fields of .1000 gauss in the lunar environment. From the permeability measurement, lunar free iron abundance is determined to be 2.5 + or - 2.0 wt. %. Total iron abundance is calculated for two assumed compositional models of the lunar interior: a free iron/orthopyroxene lunar composition and a free iron/olivine composition. The overall lunar total iron abundance is determined to be 9.0 + or - 4.7 wt. %. Other lunar models with a small iron core and with a shallow iron-rich layer are discussed in light of the measured global permeability. Effects on permeability and iron content calculations due to a possible lunar ionosphere are also considered.

Parkin, C. W.↗

Aeroacoustic Study of a Subscale Large Civil Transport (STAR) Model – Part 1: Simulations

Results from a computational aeroacoustic study of the Subsonic Transport Aeroacoustic Research (STAR) model are presented. The model, a 26%-scale semispan Boeing 777-200 aircraft, was an inceptive configuration in a comprehensive simulation campaign to accurately predict the airframe noise of a full-scale, large civil transport during landing. The lengthy process of obtaining a high-fidelity, CFD-ready, vetted, digital representation of the model is described in detail. The simulations were performed with the STAR model installed in the NASA Ames 40- by 80-foot full-scale wind tunnel to enable direct comparisons with aerodynamic and acoustic measurements of the model, presented separately in Part 2 of the study. Two versions of the wind tunnel test section were considered. The entire test section was used to simulate the aerodynamic behavior of the model; an open-wall representation of the section was adopted during the acoustic simulations to eliminate wall reflections. To preserve the aerodynamic equivalency of the model for the two simulated tunnel configurations, the open-wall computations were executed at an angle-of-attack one degree higher than that of the corresponding closed-wall tunnel. The effect on aerodynamic behavior of several minor geometry alterations necessary to con-form the simulated to the tested model was assessed – the modifications considerably altered the strength and extent of a flow separation zone on the inboard segment of the flap. Also evaluated was the impact of spatial resolution on aerodynamic and acoustic characteristics of the model – in general, medium resolution was sufficient to establish convergence trends for averaged surface pressure behavior and proper localization of noise sources. After the aeroacoustic characteristics of the baseline configuration were predicted, a toboggan-shaped fairing designed to reduce main landing gear noise was evaluated computationally. Far-field noise spectra for the baseline and toboggan-equipped STAR model were calculated via a Ffowcs-Williams and Hawkings integral approach, with flow quantities on a permeable data surface enclosing the source regions used as input. Results from the simulations indicated that, while the fairing reduced gear noise in the near-field, it had a negligible effect on the far-field noise signature of the model. This trend agreed with previously published experimental results obtained during sub-scale, isolated gear tests and full-scale flight tests

airframe noise↗

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

Inhomogeneous diffusion-limited aggregation

It is demonstrated here that inhomogeneous diffusion-limited aggregation (DLA) model can be used to simulate viscous fingering in a medium with inhomogeneous permeability and homogeneous porosity. The medium consists of a pipe-pore square-lattice network in which all pores have equal volume and the pipes have negligible volume. It is shown that fluctuations in a DLA-based growth process may be tuned by noise reduction, and that fluctuations in the velocity of the moving interface are multiplicative in form.

Selinger, Robin Blumberg↗

Iron abundance and magnetic permeability of the moon

A larger set of simultaneous data from the Apollo 12 lunar surface magnetometer and the Explorer 35 Ames magnetometer are used to construct a whole-moon hysteresis curve, from which a new value of global lunar permeability is determined to be mu = 1.012 + or - 0.006. The corresponding global induced dipole moment is 2.1 x 10 to the 18th power gauss-cu cm for typical inducing fields of .0001 gauss in the lunar environment. From the permeability measurement, lunar free iron abundance is determined to be 2.5 + or - 2.0 wt. %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior: a free iron/orthopyroxene lunar composition and a free iron/olivine composition. The overall lunar total iron abundance is determined to be 9.0 + or - 4.7 wt. %. Other lunar models with a small iron core and with a shallow iron-rich layer are discussed in light of the measured global permeability.

Parkin, C. W.↗

Iron abundance and magnetic permeability of the moon

A larger set of simultaneous data from the Apollo 12 lunar surface magnetometer and the Explorer 35 Ames magnetometer are used to construct a whole-moon hysteresis curve, from which a new value of global lunar permeability is determined to be mu = 1.012 + or - 0.006. The corresponding global induced dipole moment is 2.1 times ten to the eighteenth power gauss-cu cm for typical inducing fields of one ten-thousandth gauss in the lunar environment. From the permeability measurement, lunar free iron abundance is determined to be 2.5 + or - 2.0 wt %. Total iron abundance (sum of iron in the ferromagnetic and paramagnetic states) is calculated for two assumed compositional models of the lunar interior: a free iron/orthopyroxene lunar composition and a free iron/olivine composition. The overall lunar total iron abundance is determined to be 9.0 + or - 4.7 wt %. Other lunar models with a small iron core and with a shallow iron-rich layer are discussed in light of the measured global permeability. Effects on permeability and iron content calculations due to a possible lunar ionosphere are also considered.

Parkin, C. W.↗

RFI and SCRIMP Model Development and Verification

Vacuum-Assisted Resin Transfer Molding (VARTM) processes are becoming promising technologies in the manufacturing of primary composite structures in the aircraft industry as well as infrastructure. A great deal of work still needs to be done on efforts to reduce the costly trial-and-error methods of VARTM processing that are currently in practice today. A computer simulation model of the VARTM process would provide a cost-effective tool in the manufacturing of composites utilizing this technique. Therefore, the objective of this research was to modify an existing three-dimensional, Resin Film Infusion (RFI)/Resin Transfer Molding (RTM) model to include VARTM simulation capabilities and to verify this model with the fabrication of aircraft structural composites. An additional objective was to use the VARTM model as a process analysis tool, where this tool would enable the user to configure the best process for manufacturing quality composites. Experimental verification of the model was performed by processing several flat composite panels. The parameters verified included flow front patterns and infiltration times. The flow front patterns were determined to be qualitatively accurate, while the simulated infiltration times over predicted experimental times by 8 to 10%. Capillary and gravitational forces were incorporated into the existing RFI/RTM model in order to simulate VARTM processing physics more accurately. The theoretical capillary pressure showed the capability to reduce the simulated infiltration times by as great as 6%. The gravity, on the other hand, was found to be negligible for all cases. Finally, the VARTM model was used as a process analysis tool. This enabled the user to determine such important process constraints as the location and type of injection ports and the permeability and location of the high-permeable media. A process for a three-stiffener composite panel was proposed. This configuration evolved from the variation of the process constraints in the modeling of several different composite panels. The configuration was proposed by considering such factors as: infiltration time, the number of vacuum ports, and possible areas of void entrapment.

Loos, Alfred C.↗

An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators

The objective of this paper is to define empirical parameters for an initial thermal non-equilibrium porous-media model for use in Computational Fluid Dynamics (CFD) codes for simulation of Stirling regenerators. The two codes currently used at Glenn Research Center for Stirling modeling are Fluent and CFD-ACE. The codes porous-media models are equilibrium models, which assume solid matrix and fluid are in thermal equilibrium. This is believed to be a poor assumption for Stirling regenerators; Stirling 1-D regenerator models, used in Stirling design, use non-equilibrium regenerator models and suggest regenerator matrix and gas average temperatures can differ by several degrees at a given axial location and time during the cycle. Experimentally based information was used to define: hydrodynamic dispersion, permeability, inertial coefficient, fluid effective thermal conductivity, and fluid-solid heat transfer coefficient. Solid effective thermal conductivity was also estimated. Determination of model parameters was based on planned use in a CFD model of Infinia's Stirling Technology Demonstration Converter (TDC), which uses a random-fiber regenerator matrix. Emphasis is on use of available data to define empirical parameters needed in a thermal non-equilibrium porous media model for Stirling regenerator simulation. Such a model has not yet been implemented by the authors or their associates.

Tew, Roy C.↗

On the global TRM of the lunar lithosphere

A global thermoremanent magnetization (TRM) theory is developed which takes into account the direction and magnitude of a hypothetical source dipole field, the distribution of magnetic permeability and the coefficient of TRM and the cooling rate in the lithosphere, and models the low-field TRM characteristics of the material. When applied to the moon in conjunction with lunar magnetic parameter ranges now available, the theory gives a permanent lunar dipole moment which is uncertain to a factor of about 10 to the 5th power.

Srnka, L. J.↗

Optimum dimensions of power solenoids for magnetic suspension

Design optimization of power solenoids for controllable and stabilizable magnetic suspensions with force compensation in a wind tunnel is shown. It is assumed that the model of a levitating body is a sphere of ferromagnetic material with constant magnetic permeability. This sphere, with a radius much smaller than its distance from the solenoid above, is to be maintained in position on the solenoid axis by balance of the vertical electromagnetic force and the force of gravitation. The necessary vertical (axial) force generated by the solenoid is expressed as a function of relevant system dimensions, solenoid design parameters, and physical properties of the body. Three families of curves are obtained which depict the solenoid power for a given force as a function of the solenoid length with either outside radius or inside radius as a variable parameter and as a function of the outside radius with inside radius as a variable parameter. The curves indicate the optimum solenoid length and outside radius, for minimum power, corresponding to a given outside radius and inside radius, respectively.

Kaznacheyev, B. A.↗

Processing and Characterization of PETI Composites Fabricated by High Temperature VARTM

The use of composites as primary structures on aerospace vehicles has increased dramatically over the past decade. As these advanced structures increase in size and complexity, their production costs have grown significantly. A major contributor to these manufacturing costs is the requirement of elevated processing pressures, during the thermal cure, to create fully consolidated composites. For certain composite parts, high temperature vacuum assisted resin transfer molding (HT-VARTM) can offer reduced fabrication costs compared to conventional autoclave techniques. The process has been successfully used with phenylethynyl terminated imide (PETI) resins developed by NASA LaRC. In the current study, two PETI resins, LARC(TradeMark) PETI-330 and LARC(TradeMark) PETI-8 have been used to make test specimens using HT-VARTM. Based on previous work at NASA LaRC, larger panels with a quasi-isotropic lay-up were fabricated. The resultant composite specimens exhibited void contents <3% by volume depending on the type of carbon fabric preform used. Mechanical properties of the panels were determined at both room and elevated temperatures. These included open-hole compressive (OHC) and short beam shear (SBS) properties. Limited process modeling efforts were carried out including infusion times, composite panel size limitations and fabric permeability characterization. Work has also been carried out to develop new PETI based resins specifically geared towards HT-VARTM. The results of this work are presented herein.

Ghose, Sayata↗

Estimates for the Aerodynamic Coefficients of Ringsail and Disk-Gap-Band Parachutes Operating on Mars

The creation of flight dynamics simulations is often part of the design, development, and flight operations of missions that include entry, descent, and landing (EDL) in their concept of operation. Such flight dynamics simulations include models for many of the EDL sub-systems. For EDL systems that include a parachute, models of the parachute’s aerodynamic characteristics must be created for inclusion in the flight dynamics simulations. Models are presented for the aerodynamic coefficients (i.e., C(D),C(Tot),C(T),C(N), and C(m,SLCP)) of Supersonic Ringsail (SSRS) and Disk-Gap-Band (DGB) parachutes as functions of total porosity, λ(T), Mach number, M, and total angle of attack, λ(T) (when necessary). The source aerodynamic coefficients data used for creating these models were obtained from a wind tunnel test of subscale parachutes. In this wind tunnel test, subscale parachutes of both parachute types were fabricated from two different fabrics with very different permeabilities. By varying the fabric permeability, while maintaining the parachute geometry constant, it was possible to vary λ(T). The fabric permeability test data necessary for the calculation of λ(T) were obtained from samples of the same fabrics used to fabricate the subscale parachutes. Although the models for the aerodynamic coefficients are simple polynomial functions of λ(T) and M, they are capable of producing good reproductions of the source data. The (λ(T), M) domains over which these models are applicable are clearly defined. The models and their domains are applicable to flight operations on Mars. The models presented in this report are not unique – other models could be created from the source data. To make the creation of alternative numerical models, the aerodynamic coefficients source data are included in this report.

Cruz, Juan R.↗

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↗