Search NASA⌕ Search

SEARCH · Search NASA

Results for “STAGNATION REGION”

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 109 records · Page 6

Adaptive grid methods for RLV environment assessment and nozzle analysis

Rapid access to highly accurate data about complex configurations is needed for multi-disciplinary optimization and design. In order to efficiently meet these requirements a closer coupling between the analysis algorithms and the discretization process is needed. In some cases, such as free surface, temporally varying geometries, and fluid structure interaction, the need is unavoidable. In other cases the need is to rapidly generate and modify high quality grids. Techniques such as unstructured and/or solution-adaptive methods can be used to speed the grid generation process and to automatically cluster mesh points in regions of interest. Global features of the flow can be significantly affected by isolated regions of inadequately resolved flow. These regions may not exhibit high gradients and can be difficult to detect. Thus excessive resolution in certain regions does not necessarily increase the accuracy of the overall solution. Several approaches have been employed for both structured and unstructured grid adaption. The most widely used involve grid point redistribution, local grid point enrichment/derefinement or local modification of the actual flow solver. However, the success of any one of these methods ultimately depends on the feature detection algorithm used to determine solution domain regions which require a fine mesh for their accurate representation. Typically, weight functions are constructed to mimic the local truncation error and may require substantial user input. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. These weight functions can then be used to construct blending functions for algebraic redistribution, interpolation functions for unstructured grid generation, forcing functions to attract/repel points in an elliptic system, or to trigger local refinement, based upon application of an equidistribution principle. The popularity of solution-adaptive techniques is growing in tandem with unstructured methods. The difficultly of precisely controlling mesh densities and orientations with current unstructured grid generation systems has driven the use of solution-adaptive meshing. Use of derivatives of density or pressure are widely used for construction of such weight functions, and have been proven very successful for inviscid flows with shocks. However, less success has been realized for flowfields with viscous layers, vortices or shocks of disparate strength. It is difficult to maintain the appropriate mesh point spacing in the various regions which require a fine spacing for adequate resolution. Mesh points often migrate from important regions due to refinement of dominant features. An example of this is the well know tendency of adaptive methods to increase the resolution of shocks in the flowfield around airfoils, but in the incorrect location due to inadequate resolution of the stagnation region. This problem has been the motivation for this research.

Thornburg, Hugh J.↗

Modeling of surface roughness effects on glaze ice accretion

A series of experimental investigations focused on studying the cause and effect of roughness on accreting glaze ice surfaces were conducted. Detailed microvideo observations were made of glaze ice accretions on 1 to 4 inch diameter cylinders in three icing wind tunnels (the Data Products of New England six inch test facility, the NASA Lewis Icing Research Tunnel, and the B. F. Goodrich Ice Protection Research Facility). Infrared thermal video recordings were made of accreting ice surfaces in the Goodrich facility. Distinct zones of surface water behavior were observed; a smooth wet zone in the stagnation region with a uniform water film; a rough zone where surface tension effects caused coalescence of surface water into stationary beads; a horn zone where roughness elements grow into horn shapes; a runback zone where surface water ran back as rivulets; and a dry zone where rime feathers formed. The location of the transition from the smooth to the rough zone was found to migrate with time towards the stagnation point. The behavior of the transition appeared to be controlled by boundary layer transition and bead formation mechanisms at the interface between the smooth and rough zones. Regions of wet ice growth and enhanced heat transfer were clearly visible in the infrared video recordings of glaze ice surfaces. A simple multi-zone modification to the current glaze ice accretion model was proposed to include spatial variability in surface roughness.

Hansman, R. John, Jr.↗

Calculation of radiation from argon shock layers

The accuracy of calculations of the radiation emissions from argon plasmas produced by the shock layers over blunt bodies is assessed. The existing theoretical and experimental spectroscopic data on argon are collated. A set of such data is selected for use in the radiative transfer calculations. Calculations are performed for the stagnation regions of the shock layers over laboratory-sized models using these data, and the results are compared with the existing experimental results obtained in a shock-tube. Through this comparison and a parametric study it is shown that radiative heat fluxes at the stagnation point in an argon environment can be calculated within an uncertainty of about 15%. It is shown also that radiative heat fluxes of the order of 100 kW/sq cm can be produced in the existing laboratory facilities.

Park, C.↗

Analysis of PICA-NuSil at the Hypersonic Materials Environmental Test System (HyMETS)

Phenolic Impregnated Carbon Ablator (PICA) gained heritage during the stardust mission and is the baseline material for missions to Mars. The phenolic phase of PICA is friable, and therefore the surface is treated with a siloxane resin (NuSil) to mitigate the shedding of particulate matter. Historically, material response models have not accounted for the presence of this coating. Therefore, a test campaign was executed at the Hypersonic Materials Environmental Test System (HyMETS) to investigate the performance of NuSil under multiple heating rates and gas composition. Sphere-cone models were instrumented with thermocouples to measure the in-depth temperature response, and surface temperatures were monitored with a pyrometer. Furthermore, emission spectrometers were used to collect spectra from species emitting at the surface of the test article and the post-shock stagnation region. Post-test analysis suggests that the coating decomposes according to multiple stages. Inspection of silicon emission data reveals a rapid rise and decline in intensity within the first couple of seconds of testing, which is attributed to pyrolysis of NuSil. Pyrolysis is followed by a stagnation in silicon emission and surface temperature which are attributed to the formation of a thin layer of oxidation-resistant silicon oxycarbide. Moreover, the presence of the silicon oxycarbide layer significantly reduces the in-depth thermal response when compared to the baseline material. Eventually, the silicon oxycarbide layer decomposes via high-temperature reactions, which expose the underlying char layer to a highly reactive flow.

PICA-NuSil↗

Investigation of surface water behavior during glaze ice accretion

A series of experimental investigations that focused on isolating the primary factors that control the behavior of unfrozen surface water during glaze ice accretion were conducted. Detailed microvideo observations were made of glaze ice accretions on 2.54 cm diam cylinders in a closed-loop refrigerated wind tunnel. Distinct zones of surface water behavior were observed; a smooth wet zone in the stagnation region with a uniform water film, a rough zone where surface tension effects caused coalescence of surface water into stationary beads, and a zone where surface water ran back as rivulets. The location of the transition from the smooth to the rough zone was found to migrate towards the stagnation point with time. Comparative tests were conducted to study the effect of the substrate thermal and roughness properties on ice accretion. The importance of surface water behavior was evaluated by the addition of a surface tension reducing agent to the icing tunnel water supply, which significantly altered the accreted glaze ice shape. Measurements were made to determine the contact angle behavior of water droplets on ice. A simple multizone modification to current glaze ice accretion models was proposed to include the observed surface roughness behavior.

Hansman, R. John, Jr.↗

The effect of a turbulent wake on the stagnation point. I - Skin friction results

The response of a boundary layer in the stagnation region of a two-dimensional body to fluctuations in the freestream is examined. The analysis is restricted to laminar incompressible flow. The assumed form of the velocity distribution at the edge of the boundary layer represents both a pulsation of the incoming flow, and an oscillation of the stagnation point streamline. Both features are essential in accurately representing the effect which freestream spatial and temporal nonuniformities have upon the unsteady boundary layer. Finally, a simple model is proposed which relates the characteristic parameters in a turbulent wake to the unsteady boundary-layer edge velocity. Numerical results are presented for both an arbitrary two-dimensional geometry and a circular cylinder.

Wilson, Dennis E.↗

Motion and Heating During Atmosphere Reentry of Space Vehicles

The results of an analysis of the motion and heating during atmospheric reentry of manned space vehicles has shown the following: 1. Flight-corridor depths which allow reentry in a single pass decrease rapidly as the reentry speed increases if the maximum deceleration is limited to 10 g. 2. Use of aerodynamic lift can result in a three-to five fold increase in corridor depth over that available to a ballistic vehicle for the same deceleration limits. 3. Use of aerodynamic lift to widen these reentry corridors causes a heating penalty which becomes severe for values of the lift-drag ratio greater than unity for constant lift-drag entry. 4. In the region of most intense convective heating the inviscid flow is generally in chemical equilibrium but the boundary-layer flows are out of equilibrium. Heating rates for the nonequilibrium boundary layer are generally lower than for the corresponding equilibrium case. 5. Radiative heating from the hot gas trapped between the shock wave and the body stagnation region may be as severe as the convective heating and unfortunately occurs at approximately the same time in the flight.

Wong, Thomas J.↗

The effect of electron temperature and impact ionization on Martian return AOTV flowfields

Various electron impact ionization models in conjunction with a quasi-equilibrium electron temperature model have been investigated and applied to the stagnation region of a hypothetical 2.3 m nose radius Martian return AOTV. For the conditions considered, U = 12 km/sec at 80 km, both multi-temperature inviscid and viscous results indicate that a two-step ionization impact model predicts ionization distances in agreement with experimental data, that nonequilibrium chemistry and radiation effects are important throughout the stagnation zone, and that the quasi-equilibrium electron temperature model is reasonable. Also, using a non-grey emission-absorption radiation step model, it is shown that nonequilibrium causes a reduction in radiative heating from that predicted for equilibrium conditions and that compared to an adiabatic wall a cool wall (1650 deg K) results in a 28 to 45 percent reduction in radiative heating due to absorption near the wall.

Carlson, Leland A.↗

Characterization of Ice Roughness Variations in Scaled Glaze Icing Conditions

Because of the significant influence of surface tension in governing the stability and breakdown of the liquid film in flooded stagnation regions of airfoils exposed to glaze icing conditions, the Weber number is expected to be a significant parameter governing the formation and evolution of ice roughness. To investigate the influence of the Weber number on roughness formation, 53.3-cm (21-in.) and 182.9-cm (72-in.) NACA 0012 airfoils were exposed to flow conditions with essentially the same Weber number and varying stagnation collection efficiency to illuminate similarities of the ice roughness created on the different airfoils. The airfoils were exposed to icing conditions in the Icing Research Tunnel (IRT) at the NASA Glenn Research Center. Following exposure to the icing event, the airfoils were then scanned using a ROMER Absolute Arm scanning system. The resulting point clouds were then analyzed using the self-organizing map approach of McClain and Kreeger (2013) to determine the spatial roughness variations along the surfaces of the iced airfoils. The roughness characteristics on each airfoil were then compared using the relative geometries of the airfoil. The results indicate that features of the ice shape and roughness such as glaze-ice plateau limits and maximum airfoil roughness were captured well by Weber number and collection efficiency scaling of glaze icing conditions. However, secondary ice roughness features relating the instability and waviness of the liquid film on the glaze-ice plateau surface are scaled based on physics that were not captured by the local collection efficiency variations.

Roughness Formation↗

Nonequilibrium viscous flow over Jovian entry probes at high altitudes

The viscous chemical nonequilibrium flow around a Jovian entry body is investigated at high altitudes using two different methods. First method is only for the stagnation region and integrates the full Navier-Stokes equations from the body surface to the freestream. The second method uses viscous shock layer equations between the body surface and the shock. Due to low Reynolds numbers, both methods use surface slip boundary conditions and the second method also uses shock slip boundary conditions. The results of the two methods are compared at the stagnation point. It is found that the entire shock layer is under chemical nonequilibrium at higher altitudes and that the slip boundary conditions are important at these altitudes.

Kumar, A.↗

A numerical study of hypersonic stagnation heat transfer predictions at a coordinate singularity

The problem of grid induced errors associated with a coordinate singularity on heating predictions in the stagnation region of a three-dimensional body in hypersonic flow is examined. The test problem is for Mach 10 flow over an Aeroassist Flight Experiment configuration. This configuration is composed of an elliptic nose, a raked elliptic cone, and a circular shoulder. Irregularities in the heating predictions in the vicinity of the coordinate singularity, located at the axis of the elliptic nose near the stagnation point, are examined with respect to grid refinement and grid restructuring. The algorithm is derived using a finite-volume formulation. An upwind-biased total-variation diminishing scheme is employed for the inviscid flux contribution, and central differences are used for the viscous terms.

Grasso, Francesco↗

Radiative Heating in MSL Entry: Comparison of Flight Heating Discrepancy to Ground Test and Predictive Models

During the recent entry of the Mars Science Laboratory (MSL), the heat shield was equipped with thermocouple stacks to measure in-depth heating of the thermal protection system (TPS). When only convective heating was considered, the derived heat flux from gauges in the stagnation region was found to be underpredicted by as much as 17 W/sq cm, which is significant compared to the peak heating of 32 W/sq cm. In order to quantify the contribution of radiative heating phenomena to the discrepancy, ground tests and predictive simulations that replicated the MSL entry trajectory were performed. An analysis is carried through to assess the quality of the radiation model and the impact to stagnation line heating. The impact is shown to be significant, but does not fully explain the heating discrepancy.

MSL↗

Stagnation Point Heat Transfer with Gas Injection Cooling

The present paper deals with an experimental study of the stagnation-point heat transfer to a cooled copper surface with gas injection under subsonic conditions. Test were made with a probe that combined a steady-state water-cooled calorimeter that allows the capability to study convective blockage and to perform heat transfer measurements in presence of gas injection in the stagnation region. The copper probe was pierced by 52 holes, representing 2.4% of the total probe surface. The 1.2 MW high enthalpy plasma wind tunnel was operated at anode powers between 130 and 230 kW and a static pressures from 35 hPa up to 200 hPa. Air, carbon dioxide and argon were injected in the mass flow range 0-0.4 g/s in the boundary layer developed around the 50 mm diameter probe. The measured stagnation-point heat transfer rates are reported and discussed.

Vancrayenest, B.↗

Design fabrication, testing, and delivery of shuttle heat pipe leading edge test modules. Volume 1: Executive summary

The construction of two test modules is presented for a feasibility demonstration of a concept for reusable space shuttle wing leading edge surfaces. In this leading edge concept high temperature heat pipes were incorporated into the structure to cool the stagnation region, allowing the use of super-alloys in place of refractory metal, ablator protected, or carbon-carbon structures. The program included the analysis and design of the heat pipes, their integration into the test module structure, heat pipe development testing, construction of the test modules and a facility adapter, and formulation of recommended testing conditions. The results of the heat pipe and leading edge module thermal analyses indicate the test modules will meet the design goal; reducing the leading edge temperature at the stagnation line from 1315 C (2400 F) to less than 1010 C (1850 F). The development tests demonstrated that the module assembly could be brazed with active heat pipes, as was borne out by the subsequent successful brazing of both modules with active heat pipes loaded with sodium.

Source record↗

Radiative Heating for MSL Entry: Verification of Simulations from Ground Test to Flight Data

The heat shield of the Mars Science Laboratory (MSL) was equipped with thermocouplestacks to measure in-depth heating of the thermal protection system during atmosphericentry. The heat load derived from the thermocouples in the stagnation region was found tobe 33 lower than corresponding post-flight predictions of convective heating alone. It washypothesized that this difference could be attributed to radiation from the shock-heated gas,a mechanism not considered in post-flight analyses of flow fields. In order to test thehypothesis and quantify the contribution of shock-layer radiation to total surface heating,ground tests and simulations (both flow and radiation) were performed at several pointsalong the best-estimated entry trajectory of MSL. The present paper provides anassessment of the quality of the radiation model and its impact to stagnation point heating.Although the impact of radiative heating is shown to be significant, it only accounts for 43of the discrepancy. Additional factors behind the remaining discrepancy are discussed.

Aeroheating↗

Experimental determination of the thermal, turbulent, and rotational ion motion and magnetic field profiles in imploding plasmas

A tutorial is presented on advances in spectroscopic diagnostic methods developed for measuring key plasma properties in pulsed-power systems such as Z-pinches, magnetized-plasma compression devices, ion and electron diodes, and plasma switches. The parameters measured include the true ion temperature in Z-pinch implosions, which led to a discovery that much of the ion kinetic energy at stagnation is stored in hydrodynamic rather than in thermal motion. This observation contributed a new important insight into the understanding of the ion thermalization at stagnation and stimulated further investigations of turbulence at stagnation, discussed here too. The second part of this tutorial is devoted to the development of measurements for magnetic-field distributions in Z-pinches and in other pulsed-power systems, as well as their use in studying the plasma dynamics, resistivity, and pressure and energy balance. The latter study raises intriguing questions on the implosion process. In particular, in Z-pinches, the current during stagnation was found to largely flow at relatively large radii, outside the stagnation region. The magnetic-field measurements also enable investigations into the compression of a pre-magnetized cylindrical plasma that uncover striking phenomena related to the current flow, where the current was found to redistribute toward the outer regions during the implosion. Observation of the rotation of the magnetized plasma is also discussed. Lastly, experimental and theoretical investigations of a non-diffusive fast penetration of magnetic field into a low-density plasma, including its effect on the plasma dynamics, are described

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Granular Flow in Novel Octet Shape–Based Lattice Frame Material

Moving packed-bed heat exchangers in concentrated solar power (CSP) plants involves heat transfer between heated falling particles and supercritical carbon dioxide. The overall effective thermal conductivity of the moving packed bed and particle-side channel contact resistances are still the bottlenecks in achieving the desirable thermal transport levels. To this end, a novel moving packed bed heat exchanger consisting of an Octet lattice packed between the walls of the particle-side channel is proposed in this study. Granular flow analysis in Octet lattice moving packed bed heat exchanger (OLHX) was conducted through experiments and discrete element method (DEM)-based numerical simulations. The experimental images clearly demonstrated stagnation regions upstream of lattice fibers, void regions downstream of the fiber junctions, and wavy-type unobstructed flow on the lateral sides of the fibers. DEM simulations were successful in capturing all these critical flow phenomena. Larger flow velocities were observed on the lateral sides of the fibers in the simulations. Also, when the particles in the silo were emptied, the final images showed an accumulation of particles on the inter-fiber as well as fiber–channel wall junctions. Moreover, the fiber connections resulted in some regions devoid of particle contact on the channel endwall, which means that these regions would suffer from poor thermal exchange. Finally, the overall mass flowrate increased with increasing porosity for a fixed particle diameter.

14 SOLAR ENERGY↗

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗