Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mesh refinement”

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.

604 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

Three-Dimensional Aeroelastic and Aerothermoelastic Behavior in Hypersonic Flow

The aeroelastic and aerothermoelastic behavior of three-dimensional configurations in hypersonic flow regime are studied. The aeroelastic behavior of a low aspect ratio wing, representative of a fin or control surface on a generic hypersonic vehicle, is examined using third order piston theory, Euler and Navier-Stokes aerodynamics. The sensitivity of the aeroelastic behavior generated using Euler and Navier-Stokes aerodynamics to parameters governing temporal accuracy is also examined. Also, a refined aerothermoelastic model, which incorporates the heat transfer between the fluid and structure using CFD generated aerodynamic heating, is used to examine the aerothermoelastic behavior of the low aspect ratio wing in the hypersonic regime. Finally, the hypersonic aeroelastic behavior of a generic hypersonic vehicle with a lifting-body type fuselage and canted fins is studied using piston theory and Euler aerodynamics for the range of 2.5 less than or equal to M less than or equal to 28, at altitudes ranging from 10,000 feet to 80,000 feet. This analysis includes a study on optimal mesh selection for use with Euler aerodynamics. In addition to the aeroelastic and aerothermoelastic results presented, three time domain flutter identification techniques are compared, namely the moving block approach, the least squares curve fitting method, and a system identification technique using an Auto-Regressive model of the aeroelastic system. In general, the three methods agree well. The system identification technique, however, provided quick damping and frequency estimations with minimal response record length, and therefore o ers significant reductions in computational cost. In the present case, the computational cost was reduced by 75%. The aeroelastic and aerothermoelastic results presented illustrate the applicability of the CFL3D code for the hypersonic flight regime.

McNamara, Jack J.↗
DOE OSTI2025

Refining NEAMS MARVEL Reactor Model Accuracy

This report outlines the progress of Idaho National Laboratory (INL) in improving the fidelity of the multiphysics model of the Microreactor Applications Research Validation and Evaluation (MARVEL) reactor developed in FY24. The multiphysics model firstly developed in FY24 leverages three single-physics models coupled via the MOOSE’s MultiApps and Transfer systems. The first single-physics model, that functions as the driver application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). The second single-physics model uses BISON to compute the solid temperature in the reactor. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. In this report we detail the following improvement and enhancements: • An updated mesh leveraging the latest feature in the reactor module (see Section 3.1.).

42 - ENGINEERING↗
NASA NTRS1990

Applications of an adaptive unstructured solution algorithm to the analysis of high speed flows

An upwind cell-centered scheme for the solution of steady laminar viscous high-speed flows is implemented on unstructured two-dimensional meshes. The first-order implementation employs Roe's (1981) approximate Riemann solver, and a higher-order extension is produced by using linear reconstruction with limiting. The procedure is applied to the solution of inviscid subsonic flow over an airfoil, inviscid supersonic flow past a cylinder, and viscous hypersonic flow past a double ellipse. A detailed study is then made of a hypersonic laminar viscous flow on a 24-deg compression corner. It is shown that good agreement is achieved with previous predictions using finite-difference and finite-volume schemes. However, these predictions do not agree with experimental observations. With refinement of the structured grid at the leading edge, good agreement with experimental observations for the distributions of wall pressure, heating rate and skin friction is obtained.

Thareja, R. R.↗
NASA NTRS1993

An upwind multigrid method for solving viscous flows on unstructured triangular meshes

A multigrid algorithm is combined with an upwind scheme for solving the two dimensional Reynolds averaged Navier-Stokes equations on triangular meshes resulting in an efficient, accurate code for solving complex flows around multiple bodies. The relaxation scheme uses a backward-Euler time difference and relaxes the resulting linear system using a red-black procedure. Roe's flux-splitting scheme is used to discretize convective and pressure terms, while a central difference is used for the diffusive terms. The multigrid scheme is demonstrated for several flows around single and multi-element airfoils, including inviscid, laminar, and turbulent flows. The results show an appreciable speed up of the scheme for inviscid and laminar flows, and dramatic increases in efficiency for turbulent cases, especially those on increasingly refined grids.

Bonhaus, Daryl Lawrence↗
NASA NTRS1998

A Domain-Decomposed Multi-Level Method for Adaptively Refined Cartesian Grids with Embedded Boundaries

The work presents a new method for on-the-fly domain decomposition technique for mapping grids and solution algorithms to parallel machines, and is applicable to both shared-memory and message-passing architectures. It will be demonstrated on the Cray T3E, HP Exemplar, and SGI Origin 2000. Computing time has been secured on all these platforms. The decomposition technique is an outgrowth of techniques used in computational physics for simulations of N-body problems and the event horizons of black holes, and has not been previously used by the CFD community. Since the technique offers on-the-fly partitioning, it offers a substantial increase in flexibility for computing in heterogeneous environments, where the number of available processors may not be known at the time of job submission. In addition, since it is dynamic it permits the job to be repartitioned without global communication in cases where additional processors become available after the simulation has begun, or in cases where dynamic mesh adaptation changes the mesh size during the course of a simulation. The platform for this partitioning strategy is a completely new Cartesian Euler solver tarcreted at parallel machines which may be used in conjunction with Ames' "Cart3D" arbitrary geometry simulation package.

Aftosmis, M. J.↗
DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗
DOE OSTI2026

SAVY-4000 Finite-Element Drop Test Analysis

PFE Auxiliary Systems conducted drop testing on SAVY-4000 containers to evaluate structural response under 12-foot drop conditions. In support of that effort, a finite-element modeling capability was developed to simulate drop response across multiple container sizes and impact orientations. The purpose of this work was to provide a consistent analysis framework that could support interpretation of testing, compare response trends across multiple configurations, and generate quantities of interest for later comparison with experimental data. More broadly, the analysis and testing were intended to assess whether the containers continued to perform their primary function after a 12-foot drop, namely maintaining structural integrity and containment of the contents. The modeling approach combined an implicit preload analysis with an explicit drop simulation so that each drop event began from a mechanically realistic assembled condition, including compression of the silicone O-ring. Separate models were developed for 2-quart, 5-quart, 12-quart, and 10-gallon containers. The results were evaluated in terms of strain-gauge response, collar-lid gap behavior, and accumulated plastic strain. In addition, parametric studies were performed on the 2-quart container to assess sensitivity to O-ring stiffness, friction, canister thickness, geometry tolerance, and mesh density. The simulations showed that predicted drop responses depended strongly on both container size and drop orientation. Gap metrics identified cases in which the predicted collar-lid opening exceeded the nominal O-ring cross-section threshold, while plastic strain metrics identified localized regions of elevated permanent deformation. Parametric studies showed that the predicted response was especially sensitive to the assumed O-ring stiffness and contact friction, while the geometry tolerance study produced smaller changes in the cases examined. The main value of this work was that it established a repeatable modeling and simulation workflow to support drop-test implementation, evaluate effects of future configuration changes, and understand modeling assumptions that most influenced predicted response. At the current stage, the results were viewed as preliminary model predictions rather than validated predictions. The next step would be to compare drop-test data to the model so that predictive values of the workflow could be refined and used with greater confidence to assess whether the containers maintained structural integrity and containment of the contents after a 12-foot drop.

42 ENGINEERING↗
NASA NTRS2007

The Use of Tribocharging in the Electrostatic Beneficiation of Lunar Simulant

Any future lunar base and habitat must be constructed from strong dense materials in order to provide for thermal and radiation protection. Lunar soil may meet this need. Lunar regolith has high concentrations of aluminum, silicon, calcium, iron, sodium, and titanium oxides. Refinement or enrichment of specific minerals in the soil before it is chemically processed may be more desirable as it would reduce the size and energy requirements required to produce the virgin material and it may significantly reduce the process' complexity. Also, investigations into the potential production of breathable oxygen from oxidized mineral components are a major research initiative by NASA. In this study. the objective was to investigate the use of tribocharging to charge lunar simulants and pass them through a parallel plate separator to enrich different mineral fractions. This technique takes advantage of the high Lunar vacuum in which much higher voltages can be used on the separation plates than in air. Additionally, the Lunar g1avity, only being 1/6 that of Earth, allows the particles more separation time between the plates and therefore enhances separation. For the separation studies, two lunar stimulants were used. The first simulant was created in-house, labeled KSC-1. using commercially supplied (sieved to 325 mesh) materials, and was composed of 40 wt. % feldspar ((Na,K,Ca)AlSi3O8;SiO2), 40 wt. % olivine ((Mg,Fe)2SiO4), 10 wt. % ilmenite (FeTiO3). and 10 wt. % spodumene (LiAlSi2O6) (pyroxene). The advantage of the in-house mixture is that the composition can he varied to simulate different soil compositions from different areas on the moon. This simulant was used to show proof-of-concept using the designed separator in air. The second stimulant was JSC-1. used for the vacuum experiments. JSC-1 is principally basalts, containing phases of plagioclase. pyroxene. olivine, and ilmenite. The JSC-1 was sieved to provide a 50-75 micron size range to correlate with the mean grain size found on the moon's surface [1]. Four different materials were investigated for the triboelectrification process; aluminum, copper. stainless steel, and PTFE. These materials were selected because they offer a wide variation in work functions (aluminum 4.28 eV, copper 4.65 eV. stainless steel 5.04 eV, and PTFE 5.75 eV). The difference between the work function of each material and the simulant influences the charge obtained by the grains. Each simulant was analyzed before and after separation using X-ray Photoelectron Spectroscopy (XPS) to determine mineral surface composition. In addition. Raman spectroscopy was performed on the JSC-1 before and after separation in vacuum to determine the mineral composition. Charge-to-mass (Q/M) measurements were performed using a fluidizing bed in air and passing the simulant through a static mixer of a particular material and collecting it in a Faraday pail grounded through an electrometer. To measure the Q/M in vacuum, a special device was constructed consisting of a heater/shaker cup that fed into a solid block of material (either PTFE, copper, or aluminum) in which a channel composed of a "zig-zag" series of inclines greater than 50 degrees has been cut. The voltage to the vibrating motor can be varied to control the amount of simulant passing through the channel. Figure I shows the Q/M measurements for JSC-1 tribocharged using the static mixers and the incline plane chargers in air, and the incline plane chargers in vacuum.

Trigwell, S.↗
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

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↗