Search NASA⌕ Search

SEARCH · Search NASA

Results for “wall modeling”

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 1,009 records · Page 56

Development of an algebraic turbulence model for analysis of propulsion flows

A simple turbulence model that will be applicable to propulsion flows having both wall bounded and unbounded regions was developed and installed within the PARC Navier-Stokes code by linking two existing algebraic turbulence models. The first is the Modified Mixing Length (MML) model which is optimized for wall bounded flows. The second is the Thomas model, the standard algebraic turbulence model in PARC which has been used to calculate both bounded and unbounded turbulent flows but was optimized for the latter. This paper discusses both models and the method employed to link them into one model (referred to as the MMLT model). The PARC code with the MMLT model was applied to two dimensional turbulent flows over a flat plate and over a backward facing step to validate and optimize the model and to compare its predictions to those obtained with the three turbulence models already available in PARC.

Georgiadis, N. J.↗

The wall of reconnection-driven magnetohydrodynamic turbulence in a large solar flare

LaRosa and Moore (1993) recently proposed that the bulk dissipation of magnetic field that is required for the electron energization in the explosive phase of solar flares occurs in a 'fat current sheet', a wall of cascading magnetohydrodynamic (MHD) turbulence sustained by highly disordered driven reconnection of opposing magnetic fields impacting at a turbulent boundary layer. Here, we use the well-observed great two-ribbon eruptive flare of 1984 April 24/25 to assess the feasibility of both (1) the standard model for the overall three-dimensional form and action of the magnetic field and (2) the turbulent reconnection wall within it. We find (1) that the morphology of this flare closely matched that of the standard model; (2) the preflare sheared core field had enough nonpotential magnetic energy to power the flare; (3) the model turbulent wall required to achieve the flare's peak dissipative power easily fit within the overall span of the flaring magnetic field; (4) this wall was thick enough to have turbulent eddies large enough (diameters approximately 10(exp 8 cm) to produce the approximately ergs energy release fragments typically observed in the explosive phase of flares; (5) the aspect ratio (thickness/vertical extent) of the turbulent reconnection wall was in the 0.1-1 range expected by (Parker 1973). We therefore conclude that the viability of our version of the standard model (i.e., having the magnetic field dissipation occur in our turbulent reconnection wall) is well confirmed by this typical great two-ribbon eruptive flare.

Moore, R. L.↗

Assessing the hygrothermal performance of bio-based materials in building wall systems

Building envelope systems are crucial in regulating thermal and moisture exchange between interior and exterior environments, accounting for approximately 28 % of building energy consumption in the United States with walls being the primary contributors. Improper selection of building envelope materials can lead to moisture-related issues, reduced resilience, and compromised durability. Hygrothermal performance assessment is a key factor in efficient building design. As such, improving the energy and hygrothermal performance of opaque wall materials, through careful assessment of material choices, is essential to enhancing building resilience, lowering energy costs, and improving occupant comfort. As the building industry seeks new strategies to reduce material energy intensity, bio-based materials emerge as a promising solution. However, their long-term hygrothermal performance in building envelope systems remains underexplored. To fill this gap, this study evaluates the hygrothermal behavior of 13 bio-based materials in residential wall systems across three U.S. climate zones. Laboratory experiments were performed to measure material properties such as density, thermal conductivity, moisture transmission, and sorption isotherms. These data were integrated into the WUFI® simulation tool to assess wall hygrothermal performance in Houston, Baltimore, and Chicago. A three-phase modeling approach was used: (1) baseline residential walls with oriented strand board (OSB) and gypsum board; (2) replacing OSB with bio-based materials; and (3) replacing drywall with bio-based materials. Results showed that the evaluated bio-based materials maintained acceptable moisture thresholds of ≤ 16 % across all climates, confirming their viability as an alternative for current sheathing materials. Furthermore, this study provides a foundation for future research and innovation in material science on the use of certain bio-based materials in high-performance, low energy use residential construction. Ultimately, providing critical data, offering a database of bio-based material properties, and supplying a simulation-based approach will help designers make informed decisions for future efficient building practices.

Bio-based materials↗

On the Fidelity of Computational Models for the Flow of Milled Loblolly Pine: A Benchmark Study on Continuum-Mechanics Models and Discrete-Particle Models

The upstream of bioenergy industry has suffered from unreliable operations of granular biomass feedstocks in handling equipment. Computational modeling, including continuum-mechanics models and discrete-particle models, offers insightful understandings and predictive capabilities on the flow of milled biomass and can assist equipment design and optimization. This paper presents a benchmark study on the fidelity of the continuum and discrete modeling approaches for predicting granular biomass flow. We first introduce the constitutive law of the continuum-mechanics model and the contact law of the coarse-grained discrete-particle model, with model parameters calibrated against laboratory characterization tests of the milled loblolly pine. Three classical granular material flow systems (i.e., a lab-scale rotating drum, a pilot-scale hopper, and a full-scale inclined plane) are then simulated using the two models with the same initial and boundary conditions as the physical experiments. The close agreement of the numerical predictions with the experimental measurements on the hopper mass flow rate, the hopper critical outlet width, the material stopping thickness on the inclined plane, and the dynamic angle of repose, clearly indicates that the two methods can capture the critical flow behavior of granular biomass. The qualitative comparison shows that the continuum-mechanics model outperforms in parameterization of materials and wall friction, and large-scale systems, while the discrete-particle model is more preferred for discontinuous flow systems at smaller scales. Industry stakeholders can use these findings as guidance for choosing appropriate numerical tools to model biomass material flow in part of the optimization of material handling equipment in biorefineries.

Jin, Wencheng↗

Modeling the Restraint of Liquid Jets by Surface Tension in Microgravity

An axisymmetric phase field model is developed and used to model surface tension forces on liquid jets in microgravity. The previous work in this area is reviewed and a baseline drop tower experiment selected 'for model comparison. A mathematical model is developed which includes a free surface. a symmetric centerline and wall boundaries with given contact angles. The model is solved numerically with a compact fourth order stencil on a equally spaced axisymmetric grid. After grid convergence studies, a grid is selected and all drop tower tests modeled. Agreement was assessed by comparing predicted and measured free surface rise. Trend wise agreement is good but agreement in magnitude is only fair. Suspected sources of disagreement are suspected to be lack of a turbulence model and the existence of slosh baffles in the experiment which were not included in the model.

Chato, David J.↗

Plasmaspheric Erosion via Plasmasphere Coupling to Ring Current Plasmas: EUV Observations and Modeling

During a geomagnetic storm on 24 May 2000, the IMAGE Extreme Ultraviolet (EUV) camera observed a plasmaspheric density trough in the evening sector at L-values inside the plasmapause. Forward modeling of this feature has indicated that plasmaspheric densities beyond the outer wall of the trough are well below model expectations. This diminished plasma condition suggests the presence of an erosion process due to the interaction of the plasmasphere with ring current plasmas. We present an overview of EUV, energetic neutral atom (ENA), and Far Ultraviolet (FUV) camera observations associated with the plasmaspheric density trough of 24 May 2000, as well as forward modeling evidence of the lie existence of a plasmaspheric erosion process during this period. FUV proton aurora image analysis, convolution of ENA observations, and ring current modeling are then presented in an effort to associate the observed erosion with coupling between the plasmasphere and ring-current plasmas.

Adrian, M. L.↗

Heating distribution comparison between asymmetric and symmetric blunt cones

An experiment was performed to compare the heating distribution between symmetric and asymmetric large-angle blunt cones, with cone angles of 100, 120, and 140 deg. These hot-wall data were obtained from models made from typical thermal protection insulation for proposed aeroassisted orbital transfer vehicles. Experimental data are compared with predictions using a boundary-layer integral matrix procedure with kinetics to determine how well the heating distribution over an asymmetric cone could be approximated using axisymmetric solutions for a cone and spherical segment. In addition, a relationship between the stagnation-point heat-transfer rate and the bow-shock standoff distance for these cones is discussed. The heat-distribution data from the symmetric and asymmetric cones were very similar. Numerical results compared well with the measured wall temperatures at the stagnation point but slightly underpredicted them over the conical portion of the models.

Stewart, D. A.↗

Ground testing for the no-vent fill of cryogenic tanks - Results of tests for a 71 cubic foot tank

NASA Lewis Research has been investigating the no-vent till method, since it is a promising approach to transfer liquid while handling the problems of low-g venting. This paper reports the results of a test series for filling a 71 cu ft tank with liquid hydrogen without venting. Twenty two tests were conducted, ten with a bottom orifice as the inlet and 12 with a spray bar. Parameters investigated included inlet saturation pressures of approximately 5, 15, and 25 psia; transfer pressures of 20, 30, and 45 psia; and various starting wall temperatures. Of the tests, only the one run at the highest wall temperature (238 R) failed to fill the tank. Test results are compared to a thermodynamic equilibrium model. Overall model-data agreement was good except for the tendency of the model to overshoot during the initial wall cool-down of the higher starting wall temperature fills.

Chato, David J.↗

Ground testing for the no-vent fill of cryogenic tanks: Results of tests for a 71 cubic foot tank

NASA Lewis Research has been investigating the no-vent fill method, since it is a promising approach to transfer liquid while handling the problems of low-g venting. This paper reports the results of a test series for filling a 71 cu ft tank with liquid hydrogen without venting. 22 tests were conducted, 10 with a bottom orifice as the inlet and 12 with a spray bar. Parameters investigated included inlet saturation pressures of approximately 5, 15, and 25 psia, transfer pressures of 20, 30, and 45 psia, and various starting wall temperatures. Of the tests, only the one run at the highest wall temperature (238 R) failed to fill the tank. Test results are compared to a thermodynamic equilibrium model. Overall model-data agreement was good except for the tendency of the model to overshoot during the initial wall cool down of the higher starting wall temperature fills.

Chato, David J.↗

Modeling of turbulence and transition

The first objective is to evaluate current two-equation and second order closure turbulence models using available direct numerical simulations and experiments, and to identify the models which represent the state of the art in turbulence modeling. The second objective is to study the near-wall behavior of turbulence, and to develop reliable models for an engineering calculation of turbulence and transition. The third objective is to develop a two-scale model for compressible turbulence.

Shih, Tsan-Hsing↗

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗

11 Foot Unitary Plan Tunnel Facility Optical Improvement Large Window Analysis

The test section of the 11 by 11-foot Unitary Plan Transonic Wind Tunnel (11-foot UPWT) may receive an upgrade of larger optical windows on both the North and South sides. These new larger windows will provide better access for optical imaging of test article flow phenomena including surface and off body flow characteristics. The installation of these new larger windows will likely produce a change to the aerodynamic characteristics of the flow in the Test Section. In an effort understand the effect of this change, a computational model was employed to predict the flows through the slotted walls, in the test section and around the model before and after the tunnel modification. This report documents the solid CAD model that was created and the inviscid computational analysis that was completed as a preliminary estimate of the effect of the changes.

computer aided design↗

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↗

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↗

Czochralski growth of crystals - Simple models for growth rate and interface shape

A simple model for the crystal growth by the Czochralski (CZ) process has been proposed based on semiquantitative arguments. The model provides empirical relationships for the dependence of the pulling rate and the interface shape on the important process variables such as crystal radius, crucible temperature, height of the melt level, and the height of the exposed portion of the crucible wall. The parameters of the model can be evaluated by matching the results obtained from a detailed mathematical model of the CZ process or from extensive experimental data. The model has, therefore, the potential application for determining the best process conditions and for on-line control and optimization of the crystal puller to grow crystals with constant diameter and nearly planar interface.

Srivastava, R. K.↗

Predicting equilibrium states with Reynolds stress closures in channel flow and homogeneous shear flow

Turbulent channel flow and homogeneous shear flow have served as basic building block flows for the testing and calibration of Reynolds stress models. A direct theoretical connection is made between homogeneous shear flow in equilibrium and the log-layer of fully-developed turbulent channel flow. It is shown that if a second-order closure model is calibrated to yield good equilibrium values for homogeneous shear flow it will also yield good results for the log-layer of channel flow provided that the Rotta coefficient is not too far removed from one. Most of the commonly used second-order closure models introduce an ad hoc wall reflection term in order to mask deficient predictions for the log-layer of channel flow that arise either from an inaccurate calibration of homogeneous shear flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations are presented to demonstrate this point which has important implications for turbulence modeling.

Abid, R.↗

Predicting equilibrium states with Reynolds stress closures in channel flow and homogeneous shear flow

Turbulent channel flow and homogeneous shear flow have served as basic building block flows for the testing and calibration of Reynolds stress models. A direct theoretical connection is made between homogeneous shear flow in equilibrium and the log-layer of fully-developed turbulent channel flow. It is shown that if a second-order closure model is calibrated to yield good equilibrium values for homogeneous shear flow it will also yield good results for the log-layer of channel flow provided that the Rotta coefficient is not too far removed from one. Most of the commonly used second-order closure models introduce an ad hoc wall reflection term in order to mask deficient predictions for the log-layer of channel flow that arise either from an inaccurate calibration of homogeneous shear flow or from the use of a Rotta coefficient that is too large. Illustrative model calculations are presented to demonstrate this point which has important implications for turbulence modeling.

Abid, R.↗

Wall interference tests of a CAST 10-2/DOA 2 airfoil in an adaptive-wall test section

A wind-tunnel investigation of a CAST 10-2/DOA 2 airfoil model has been conducted in the adaptive-wall test section of the Langley 0.3-Meter Transonic Cryogenic Tunnel (TCT) and in the National Aeronautical Establishment High Reynolds Number Two-Dimensional Test Facility. The primary goal of the tests was to assess two different wall-interference correction techniques: adaptive test-section walls and classical analytical corrections. Tests were conducted over a Mach number range from 0.3 to 0.8 and over a chord Reynolds number range from 6 million to 70 million. The airfoil aerodynamic characteristics from the tests in the 0.3-m TCT have been corrected for wall interference by the movement of the adaptive walls. No additional corrections for any residual interference have been applied to the data, to allow comparison with the classically corrected data from the same model in the conventional National Aeronautical Establishment facility. The data are presented graphically in this report as integrated force-and-moment coefficients and chordwise pressure distributions.

Mineck, Raymond E.↗