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 775 records · Page 43

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.↗

Turbulence Model Validation Through Joint Experimental /Computational Studies of Separated Flow Over A Three-Dimensional Tapered Bump: Part I - Experimental Investigation

The results of a three-year experimental campaign aimed at comprehensively documenting the separated flow over a three-dimensional bump are presented with the purpose of generating a benchmark experimental database useful in validating computational fluid dynamics (CFD) flow simulations. The bump model geometry was designed to provide well-defined and repeatable smooth-body flow separation conditions that were suitable for both experiments and simulations. The bump had a Gaussian streamwise profile with a constant height equal to 8.5% of its width over the central 60% of the test section width. The remaining 40% were outboard spanwise portions that gradually taper to zero using an error function profile to minimize side-wall boundary layer effects. The model was immersed in a turbulent boundary layer that was developed on a suspended flat plate in the Notre Dame Mach 0.6 Wind Tunnel. In order to document the effect of the incoming boundary layer thickness on the flow separation, the bump model could be located at two streamwise positions. The mean velocity and turbulence intensity of the wind tunnel freestream flow field and approaching turbulent boundary layer were fully documented. The measurements of the flow separation region included surface visualization, wall shear stress using oil-film interferometry, mean and dynamic surface pressure, and planar and stereoscopic particle image velocimetry. The experiments were conducted over a range of Mach numbers from 0.05 to 0.2 corresponding to a range of Reynolds numbers based on the test section spanwise dimension (0.914m) of 1.0 × 10 6 ≤ Re L = U ∞ L/ν ≤ 4.0 × 10 6 . The bulk of the results are presented for the higher Mach number conditions of 0.1 and 0.2 with Re L = 2.0 × 10 6 and 4.0 × 10 6 , respectively. Extensive uncertainty analysis of the data was performed. The data is archived in the NASA Langley Turbulence Modeling Resource website at https://turbmodels.larc.nasa.gov/Other_exp_Data/speedbump_sep_exp.html. In addition to the experiments, a computational effort was made in parallel by the CFD group at Boeing Research & Technology highlighting the usefulness of the data set, which is outlined in the accompanying CFD report.

Flow separation↗

Vibrations and structureborne noise in space station

Theoretical models were developed capable of predicting structural response and noise transmission to random point mechanical loads. Fiber reinforced composite and aluminum materials were considered. Cylindrical shells and circular plates were taken as typical representatives of structural components for space station habitability modules. Analytical formulations include double wall and single wall constructions. Pressurized and unpressurized models were considered. Parametric studies were conducted to determine the effect on structural response and noise transmission due to fiber orientation, point load location, damping in the core and the main load carrying structure, pressurization, interior acoustic absorption, etc. These analytical models could serve as preliminary tools for assessing noise related problems, for space station applications.

Vaicaitis, R.↗

Semi-span model testing in the National Transonic Facility

A semi-span testing technique has been proposed for the NASA Langley Research Center's National Transonic Facility (NTF). Semi-span testing has several advantages including (1) larger model size, giving increased Reynolds number capability; (2) improved model fidelity, allowing ease of flap and slat positioning which ultimately improves data quality; and (3) reduced construction costs compared with a full-span model. In addition, the increased model size inherently allows for increased model strength, reducing aeroelastic effects at the high dynamic pressure levels necessary to simulate flight Reynolds numbers. The Energy Efficient Transport (EET) full-span model has been modified to become the EET semi-span model. The full-span EET model was tested extensively at both NASA LRC and NASA Ames Research Center. The available full-span data will be useful in validating the semi-span test strategy in the NTF. In spite of the advantages discussed above, the use of a semi-span model does introduce additional challenges which must be addressed in the testing procedure. To minimize the influence of the sidewall boundary layer on the flow over the semi-span model, the model must be off-set from the sidewall. The objective is to remove the semi-span model from the sidewall boundary layer by use of a stand-off geometry. When this is done however, the symmetry along the centerline of the full-span model is lost when the semi-span model is mounted on the wind tunnel sidewall. In addition, the large semi-span model will impose a significant pressure loading on the sidewall boundary layer, which may cause separation. Even under flow conditions where the sidewall boundary layer remains attached, the sidewall boundary layer may adversely effect the flow over the semi-span model. Also, the increased model size and sidewall mounting requires a modified wall correction strategy. With these issues in mind, the semi-span model has been well instrumented with surface pressure taps to obtain data on the expected complex flow field in the near wall region. This status report summarizes the progress to date on developing the semi-span geometry definition suitable for generating structured grids for the computational research. In addition, the progress on evaluating three state-of-the-art Navier-Stokes codes is presented.

Chokani, Ndaona↗

Effect of a surface-to-gap temperature discontinuity on the heat transfer to reusable surface insulation tile gaps

An experimental investigation is presented that was performed to determine the effect of a surface-to-gap wall temperature discontinuity on the heat transfer within space shuttle, reusable surface insulation, tile gaps submerged in a thick turbulent boundary layer. Heat-transfer measurements were obtained on a flat-plate, single-gap model submerged in a turbulent tunnel wall boundary layer at a nominal free-stream Mach number of 10.3 and free-stream Reynolds numbers per meter of 1.5 million, 3.3 million and 7.8 million. Surface-to-gap wall temperature discontinuities of varying degree were created by heating the surface of the model upstream of the instrumented gap. The sweep angle of the gap was varied between 0 deg and 60 deg; gap width and depth were held constant. A surface-to-gap wall temperature discontinuity (surface temperature greater than gap wall temperature) results in increased heat transfer to the near-surface portion of the gap, as compared with the heat transfer under isothermal conditions, while decreasing the heat transfer to the deeper portions of the gap. The nondimensionalized heat transfer to the near-surface portion of the gap is shown to decrease with increasing Reynolds number; in the deeper portion of the gap, the heat transfer increases with Reynolds number.

Throckmorton, D. A.↗