Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multi-zone model”

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 19 records

An Extended Multi-Zone Model for the MCG-6-30-15 Warm Absorber

The variable warm absorber seen with ASCA in the X-ray spectrum of MCG 6-30-15 shows complex time behaviour in which the optical depth of O VIII anticorrelates with the flux whereas that of O VII is unchanging. The explanation in terms of a two zone absorber has since been challenged by BeppoSAX observations. These present a more complicated behaviour for the O VII edge. The explanation we offer for both ASCA and BeppoSAX observations requires a very simple photoionization model together with the presence of a third, intermediate, zone and a period of very low luminosity. In practice warm absorbers are likely to be extended, multi-zone regions of which only part causes directly observable absorption edges at any given time depending on the value of the luminosity.

Morales, R.↗

Heat transfer on accreting ice surfaces

Based on previous observations of glaze ice accretion, a 'Multi-Zone' model with distinct zones of different surface roughness is demonstrated. The use of surface roughness in the LEWICE ice accretion prediction code is examined. It was found that roughness is used in two ways: to determine the laminar to turbulent transition location and to calculate the turbulent heat transfer coefficient. A two zone version of the Multi-Zone model is implemented in the LEWICE code, and compared with experimental heat transfer coefficient and ice accretin results. The analysis of the boundary layer transition, surface roughness, and viscous flow field effects significantly increased the accuracy in predicting heat transfer coefficients. The Multi-Zone model was found to greatly improve the ice accretion prediction for the cases compared.

Yamaguchi, Keiko↗

Deterministic multi-zone ice accretion modeling

The study focuses on a deterministic model of the surface roughness transition behavior of glaze ice and analyzes the initial smooth/rough transition location, bead formation, and the propagation of the transition location. Based on a hypothesis that the smooth/rough transition location coincides with the laminar/turbulent boundary-layer transition location, a multizone model is implemented in the LEWICE code. In order to verify the effectiveness of the model, ice accretion predictions for simple cylinders calculated by the multizone LEWICE are compared to experimental ice shapes. The glaze ice shapes are found to be sensitive to the laminar surface roughness and bead thickness parameters controlling the transition location, while the ice shapes are found to be insensitive to the turbulent surface roughness.

Yamaguchi, K.↗

Deterministic multi-zone ice accretion modeling

The focus here is on a deterministic model of the surface roughness transition behavior of glaze ice. The initial smooth/rough transition location, bead formation, and the propagation of the transition location are analyzed. Based on the hypothesis that the smooth/rough transition location coincides with the laminar/turbulent boundary layer transition location, a multizone model is implemented in the LEWICE code. In order to verify the effectiveness of the model, ice accretion predictions for simple cylinders calculated by the multizone LEWICE are compared to experimental ice shapes. The glaze ice shapes are found to be sensitive to the laminar surface roughness and bead thickness parameters controlling the transition location, while the ice shapes are found to be insensitive to the turbulent surface roughness.

Yamaguchi, K.↗

CFD applications in pump flows

The objective of the proposed paper is to develop a computational procedure that solves incompressible Navier-Stokes equations for pump flows. The solution method is based on the pseudo-compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. As a benchmark problem, the flow through the Rocketdyne inducer is numerically simulated. A coarse grid solution is obtained with a single zone by using an algebraic turbulence model. In multi-zone fine grid computation, a one-equation Baldwin-Barth turbulence model is utilized. Numerical results are compared with experimental measurements and a good agreement is found between the two. The resulting computer code is then applied to the flow analysis inside a two-stage fuel pump impeller operating at 80 percent, 100 percent, and 120 percent of design flow.

Kiris, Cetin↗

An investigation of air transportation technology at the Massachusetts Institute of Technology, 1990-1991

Brief summaries are given of research activities at the Massachusetts Institute of Technology (MIT) under the sponsorship of the FAA/NASA Joint University Program. Topics covered include hazard assessment and cockpit presentation issues for microburst alerting systems; the situational awareness effect of automated air traffic control (ATC) datalink clearance amendments; a graphical simulation system for adaptive, automated approach spacing; an expert system for temporal planning with application to runway configuration management; deterministic multi-zone ice accretion modeling; alert generation and cockpit presentation for an integrated microburst alerting system; and passive infrared ice detection for helicopter applications.

Simpson, Robert W.↗

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗

Computational/experimental analysis of three low sonic boom configurations with design modifications

The Euler code, designated AIRPLANE, which uses an unstructured tetrahedral mesh was used to compute near-field sonic boom pressure signatures on three modern low sonic boom configurations: the Mach 2, Mach 3, and Haglund models. The TEAM code which uses a multi-zoned structured grid was used to calculate pressure signatures for the Mach 2 model. The computational pressure signatures for the Mach 2 and Mach 3 models are compared with recent experimental data. The computed pressure signatures were extracted at distances less than one body length below the configuration and extrapolated to the experimental distance. The Mach 2 model was found to have larger overpressures off-ground-track than on-ground-track in both computational and experimental results. The correlations with the experiment were acceptable where the signatures were not contaminated by instrumentation and model-support hardware. AIRPLANE was used to study selected modifications to improve the overpressures of the Mach 2 model.

Cliff, Susan E.↗

Development of a computational testbed for numerical simulation of combustion instability

A synergistic hierarchy of analytical and computational fluid dynamic techniques is used to analyze three-dimensional combustion instabilities in liquid rocket engines. A mixed finite difference/spectral procedure is employed to study the effects of a distributed vaporization zone on standing and spinning instability modes within the chamber. Droplet atomization and vaporization are treated by a variety of classical models found in the literature. A multi-zone, linearized analytical solution is used to validate the accuracy of the numerical simulations at small amplitudes for a distributed vaporization region. This comparison indicates excellent amplitude and phase agreement under both stable and unstable operating conditions when amplitudes are small and proper grid resolution is used. As amplitudes get larger, expected nonlinearities are observed. The effect of liquid droplet temperature fluctuations was found to be of critical importance in driving the instabilities of the combustion chamber.

Grenda, Jeffrey↗

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

Hybrid MPI+OpenMP Programming of an Overset CFD Solver and Performance Investigations

This report describes a two level parallelization of a Computational Fluid Dynamic (CFD) solver with multi-zone overset structured grids. The approach is based on a hybrid MPI+OpenMP programming model suitable for shared memory and clusters of shared memory machines. The performance investigations of the hybrid application on an SGI Origin2000 (O2K) machine is reported using medium and large scale test problems.

Djomehri, M. Jahed↗

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.↗

Multi-Zone Liquid Thrust Chamber Performance Code with Domain Decomposition for Parallel Processing

Computational Fluid Dynamics (CFD) has considerably evolved in the last decade. There are many computer programs that can perform computations on viscous internal or external flows with chemical reactions. CFD has become a commonly used tool in the design and analysis of gas turbines, ramjet combustors, turbo-machinery, inlet ducts, rocket engines, jet interaction, missile, and ramjet nozzles. One of the problems of interest to NASA has always been the performance prediction for rocket and air-breathing engines. Due to the complexity of flow in these engines it is necessary to resolve the flowfield into a fine mesh to capture quantities like turbulence and heat transfer. However, calculation on a high-resolution grid is associated with a prohibitively increasing computational time that can downgrade the value of the CFD for practical engineering calculations. The Liquid Thrust Chamber Performance (LTCP) code was developed for NASA/MSFC (Marshall Space Flight Center) to perform liquid rocket engine performance calculations. This code is a 2D/axisymmetric full Navier-Stokes (NS) solver with fully coupled finite rate chemistry and Eulerian treatment of liquid fuel and/or oxidizer droplets. One of the advantages of this code has been the resemblance of its input file to the JANNAF (Joint Army Navy NASA Air Force Interagency Propulsion Committee) standard TDK code, and its automatic grid generation for JANNAF defined combustion chamber wall geometry. These options minimize the learning effort for TDK users, and make the code a good candidate for performing engineering calculations. Although the LTCP code was developed for liquid rocket engines, it is a general-purpose code and has been used for solving many engineering problems. However, the single zone formulation of the LTCP has limited the code to be applicable to problems with complex geometry. Furthermore, the computational time becomes prohibitively large for high-resolution problems with chemistry, two-equation turbulence model, and two-phase flow. To overcome these limitations, the LTCP code is rewritten to include the multi-zone capability with domain decomposition that makes it suitable for parallel processing, i.e., enabling the code to run every zone or sub-domain on a separate processor. This can reduce the run time by a factor of 6 to 8, depending on the problem.

Homayun K. Navaz↗

Numerical simulation of rotorcraft

The objective of the research is to develop and validate accurate, user-oriented viscous CFD codes (with inviscid options) for three-dimensional, unsteady aerodynamic flows about arbitrary rotorcraft configurations. Unsteady, three-dimensional Euler and Navier-Stokes codes are developed, adapted, and extended to rotor-body combinations. Flow solvers are coupled with zonal grid topologies, including rotating and nonrotating blocks. Special grid clustering and wave-fitting techniques were developed to capture low-level radiating acoustic waves. Significant progress was made in computing the propagation of acoustic waves due to the interaction of a concentrated vortex and a helicopter airfoil. The need for higher-order schemes was firmly established in relatively inexpensive two-dimensional calculations. In three dimensions, the number of grid points required to capture the low-level acoustic waves becomes very large, so that large supercomputer memory becomes essential. Good agreement was obtained between the numerical results obtained with a thin-layer Navier-Stokes code and experimental data from a model rotor. In addition, several nonrotating configurations that are sometimes proposed to simulate rotor blade tips in conventional wind tunnels were examined, and the complex flow around the radical tip shape of the world's fastest helicopter is under investigation. These studies demonstrate the flexibility and power of CFD to gain physical insight, study novel ideas, and examine various possibilities that might be difficult or impossible to set up in physical experiments. As a prelude to studies of rotor-body aerodynamic interactions, a preliminary grid topology and moving-interface strategy were developed. A new Euler/Navier-Stokes code using these techniques computes the vortical wake directly, rather than modeling it, as in most previous rotorcraft studies. Several hover cases were run for conventional and advanced-geometry blades. Numerical schemes using multi-zones and/or adaptive grids appear to be necessary to simulate the complex vortical flows in rotor wakes.

Mccroskey, William J.↗

The multi-zone calculation of turbomachinery flows. II - The multi-zone calculation of the turbulent, two-specie flow through the SSME HPFTP first and second stage cavities

A multi-zone Navier-Stokes methodology to calculate the two-specie flow through the first and second stage cooling cavities of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) is developed. A simplified two-component fluid formulation is used to model the interaction of coolant and hot gas. Johnston's secant approximation is used to define an appropriate near wall velocity for use in a three-dimensional law of the wall. The basic Navier-Stokes algorithm used is a finite-volume, predictor-corrector algorithm which uses a pressure correction technique. A multi-zone method is used to partition each cavity into easily handled subdomains. The results show that coolant flow is pumped up the turbine wheel for both cavities, creating a region of large temperature gradients on the turbine shank.

Williams, M.↗

Analysis of surface roughness generation in aircraft ice accretion

Patterns of roughness evolution have been studied analysis of high magnification video observations of accreting ice surfaces provided by the NASA Lewis Research Center. Three distinct patterns of surface roughness generation have been identified within the parametric regions studied. They include: Rime, Multi-Zone Glaze, and Uniform Glaze. Under most icing conditions, a brief period of transient rime ice growth was observed caused by heat conduction into the body. The resulting thin rime layer explains previously observed insensitivity of some ice accretions to substrate insensitivity of some ice accretions to substrate surface chemistry and may provide justification for simplifying assumptions in ice accretion sailing and modeling effects.

Hansman, R. J., Jr.↗

CFD Parametric Study of Consortium Impeller

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design impact in a productive manner. The main goal of this study is to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A Finite Difference Navier-Stokes flow solver, FDNS, which includes the extended k-epsilon turbulence model and appropriate moving interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. A second-order central difference scheme plus adaptive dissipation terms was employed in the FDNS code, along with a predictor plus multi-corrector pressure-based solution procedure. The multi-zone, multi-block capability allows the FDNS code to efficiently solve flow fields with complicated geometry. The FDNS code has been benchmarked by analyzing the pump consortium inducer, and it provided satisfactory results. In the present study, a CFD parametric study of the pump consortium impeller was conducted using the FDNS code. The pump consortium impeller, with partial blades, is a new design concept of the advanced rocket engines. The parametric study was to analyze the baseline design of the consortium impeller and its modification which utilizes TANDEM blades. In the present study, the TANDEM blade configuration of the consortium impeller considers cut full blades for about one quarter chord length from the leading edge and clocks the leading edge portion with an angle of 7.5 or 22.5 degrees. The purpose of the present study is to investigate the effect and trend of the TANDEM blade modification and provide the result as a design guideline. A 3-D flow analysis, with a 103 x 23 x 30 mesh grid system and with the inlet flow conditions measured by Rocketdyne, was performed for the baseline consortium impeller. The numerical result shows that the mass flow rate splits through various blade passages are relatively uniform. Due to the complexity of blade geometries, the TANDEM blade configurations were analyzed with the multi-zone grid structure. Both the 7.5 deg- and the 22.5 deg-clocking TANDEM blade cases utilized a 80K mesh system. The numerical result of two TANDEM blade modifications indicates the efficiency and the head are worse than those of the baseline case due to larger flow distortion. The gap between the TANDEM blade and the full blade allows the flow passes through and heavily loads the pressure side of the partial blade such that flow reversal occurs near the suction side of the splitter. The flow split at the exit of impeller blades is very non-uniform for TANDEM blade cases, and this will greatly induce the side load on the diffuser. Therefore, the TANDEM blade modification in the present CFD analysis does not improve the performance of the consortium impeller.

Cheng, Gary C.↗

Advanced Multi-Phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

It is known that the simulations of solid rocket motor internal flow field with AL-based propellants require complex multi-phase turbulent flow model. The objective of this study is to develop an advanced particulate multi-phase flow model which includes the effects of particle dynamics, chemical reaction and hot gas flow turbulence. The inclusion of particle agglomeration, particle/gas reaction and mass transfer, particle collision, coalescence and breakup mechanisms in modeling the particle dynamics will allow the proposed model to realistically simulate the flowfield inside a solid rocket motor. The Finite Difference Navier-Stokes numerical code FDNS is used to simulate the steady-state multi-phase particulate flow field for a 3-zone 2-D axisymmetric ASRM model and a 6-zone 3-D ASRM model at launch conditions. The 2-D model includes aft-end cavity and submerged nozzle. The 3-D model represents the whole ASRM geometry, including additional grain port area in the gas cavity and two inhibitors. FDNS is a pressure based finite difference Navier-Stokes flow solver with time-accurate adaptive second-order upwind schemes, standard and extended k-epsilon models with compressibility corrections, multi zone body-fitted formulations, and turbulence particle interaction model. Eulerian/Lagrangian multi-phase solution method is applied for multi-zone mesh. To simulate the chemical reaction, penalty function corrected efficient finite-rate chemistry integration method is used in FDNS. For the AL particle combustion rate, the Hermsen correlation is employed. To simulate the turbulent dispersion of particles, the Gaussian probability distribution with standard deviation equal to (2k/3)(exp 1/2) is used for the random turbulent velocity components. The computational results reveal that the flow field near the juncture of aft-end cavity and the submerged nozzle is very complex. The effects of the turbulent particles affect the flow field significantly and provide better prediction of the ASRM performance. The multi-phase flow analysis using the FDNS code in the present research can be used as a design tool for solid rocket motor applications.

Liaw, Paul↗