Search NASA⌕ Search

SEARCH · Search NASA

Results for “AXISYMMETRIC FLOW”

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 199 records · Page 11

Modeling Turbulent Combustion for Variable Prandtl and Schmidt Number

This report consists of two abstracts submitted for possible presentation at the AIAA Aerospace Science Meeting to be held in January 2005. Since the submittal of these abstracts we are continuing refinement of the model coefficients derived for the case of a variable Turbulent Prandtl number. The test cases being investigated are a Mach 9.2 flow over a degree ramp and a Mach 8.2 3-D calculation of crossing shocks. We have developed an axisymmetric code for treating axisymmetric flows. In addition the variable Schmidt number formulation was incorporated in the code and we are in the process of determining the model constants.

Hassan, H. A.↗

Automatic Boundary-Layer Adaptation of Structured Grids in VULCAN-CFD

In supersonic and hypersonic flow computations, well-resolved boundary layers are essential for accurate quantification of surface heating and transition prediction, particularly via linear stability analysis. Grid design for hypersonic flows with shocks, boundary-layer separation, and/or complex mean flow features incorporating spanwise/azimuthal inhomogeneities is a difficult issue. In comparison to a manual grid adaptation procedure, an autonomous grid adaptation technique offers significant improvements in computing time and solution quality. The VULCAN-CFD solver already includes a validated procedure for automatic adaption of structured grids to the bow shock. The present focus is on implementing an automatic boundary-layer adaptation capability in VULCAN-CFD that adapts structured, multiblock grids to both the bow shock and the boundary layer at the same time. The boundary-layer adaptation algorithm allows the user to specify the number of cells within the boundary layer, along with the input parameters used for detecting the edge of the boundary layer, namely, the variable used in the edge detection criterion,the edge detection method, the detection direction, and the relaxation factor used during the morphing of the grid. The algorithm automatically distributes grid points along the wall-normal direction to achieve a smooth variation in grid spacing from the edge of the boundary layer to a"junction" location within the outer part of the grid. Illustrative results are presented for three different high-speed configurations: the two-dimensional flow over a cylinder at Mach 17.6 and unit Reynolds number of Re=0.38x10^6 m^-1, the axisymmetric flow over a cone-cylinder-flare model at Mach 6.0 and Re = 10.5×10^6 m^-1, and the three-dimensional flow over a blunt, 7-degree half-angle cone at 5-degree angle of attack in a Mach 9.79 flow with Re = 17.1×10^6 m^-1. The automated boundary-layer adaptation is shown to provide an adequate grid topology that is aligned with the bow shock in the outer part of the grid and also resolves the viscous boundary-layer region close to the surface.

boundary layer transition↗

Deceleration of a supersonic flow behind a curved shock wave with isentropic precompression

Three-dimensional supersonic flows of an ideal fluid in the neighborhood of bodies formed by being cut out along the streamlines of an axisymmetric flow are investigated. The flow consists of a region of isoentropic compression and a region of vortex flow. An exact solution with variable entropy is used to describe the flow in the vortex region. In the continuous flow region an approximate solution is constructed by expanding the solution in a series in a small parameter. The effect of the shape of the excision and the vorticity of the flow on compression of the jet and and the total pressure loss coefficient is studied.

Dulov, V. G.↗

Large Scale Flows through the Solar Cycle

Large scale flows within the solar convection zone are the primary drivers of the Sun's magnetic activity cycle. Differential rotation amplifies the magnetic field and converts poloidal fields into toroidal fields. Poleward meridional flow near the surface carries magnetic flux that reverses the magnetic poles and converts toroidal fields into poloidal fields. The deeper, equatorward meridional flow carries magnetic flux toward the equator where it must reconnect with oppositely directed fields in the other hemisphere. These axisymmetric flows are themselves driven by large scale convective motions. Given these intimate connections between the large scale flows and solar activity, it would be surprising if there weren't solar cycle variations in the flow characteristics. Some variations, namely the torsional oscillations, are well established. Other variations, namely changes in the meridional flow and in the convective motions themselves, are more controversial. In this presentation I will describe the observed characteristics of the large scale flows and discuss possible variations associated with the solar cycle.

Hathaway, D. H.↗

Off-design computer code for calculating the aerodynamic performance of axial-flow fans and compressors

An off-design axial-flow compressor code is presented and is available from COSMIC for predicting the aerodynamic performance maps of fans and compressors. Steady axisymmetric flow is assumed and the aerodynamic solution reduces to solving the two-dimensional flow field in the meridional plane. A streamline curvature method is used for calculating this flow-field outside the blade rows. This code allows for bleed flows and the first five stators can be reset for each rotational speed, capabilities which are necessary for large multistage compressors. The accuracy of the off-design performance predictions depend upon the validity of the flow loss and deviation correlation models. These empirical correlations for the flow loss and deviation are used to model the real flow effects and the off-design code will compute through small reverse flow regions. The input to this off-design code is fully described and a user's example case for a two-stage fan is included with complete input and output data sets. Also, a comparison of the off-design code predictions with experimental data is included which generally shows good agreement.

Schmidt, James F.↗

How Large Scale Flows in the Solar Convection Zone may Influence Solar Activity

Large scale flows within the solar convection zone are the primary drivers of the Sun s magnetic activity cycle. Differential rotation can amplify the magnetic field and convert poloidal fields into toroidal fields. Poleward meridional flow near the surface can carry magnetic flux that reverses the magnetic poles and can convert toroidal fields into poloidal fields. The deeper, equatorward meridional flow can carry magnetic flux toward the equator where it can reconnect with oppositely directed fields in the other hemisphere. These axisymmetric flows are themselves driven by large scale convective motions. The effects of the Sun s rotation on convection produce velocity correlations that can maintain the differential rotation and meridional circulation. These convective motions can influence solar activity themselves by shaping the large-scale magnetic field pattern. While considerable theoretical advances have been made toward understanding these large scale flows, outstanding problems in matching theory to observations still remain.

Hathaway, D. H.↗

Abel Inversion of Deflectometric Measurements in Dynamic Flows

We present an Abel-inversion algorithm to reconstruct mean and rms refractive-index profiles from spatially resolved statistical measurements of the beam-deflection angle in time-dependent, axisymmetric flows. An oscillating gas-jet diffusion flame was investigated as a test case for applying the algorithm. Experimental data were obtained across the whole field by a rainbow schlieren apparatus. Results show that simultaneous multipoint measurements are necessary to reconstruct the rms refractive index accurately.

Agrawal, Ajay K.↗

A numerical method for predicting hypersonic flowfields

The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to react chemically and reach states in thermal nonequilibrium. In this paper, a new procedure based on Gauss-Seidel line relaxation is shown to solve the equations of hypersonic flow fields containing finite reaction rate chemistry and thermal nonequilibrium. The method requires a few hundred time steps and small computer times for axisymmetric flows about simple body shapes. The extension to more complex two-dimensional body geometries appears straightforward.

Maccormack, Robert W.↗

A Search for Giant Convection Cells on the Sun

Giant convection cells (slow, long-lived cellular flows in the Sun's convection zone with typical diameters greater than about 100,000 km) have been the subject of many searches over the last 30 years. If such flows exist, they should play a key role in structuring the Sun's large scale magnetic field and in driving the large scale axisymmetric flows: the differential rotation and meridional circulation. Detailed observations of the flows in these cells may also allow us to better predict future magnetic field configurations and the solar activity associated with them. Line-of-sight velocity data from the Michelson Doppler Interferometer on the ESA/NASA Solar and Heliospheric Observatory provides us with new opportunities to search for giant cells. This data is free of any atmospheric distortion and has been obtained continuously without any day/night gaps for more than two months at a time. These two-month datasets are important because giant cells are expected to have lifetimes somewhat longer that the Sun's 27 day rotation period. Any reappearance of a flow pattern after 27 days would be an important confirmation of the existence of these cells. The approach taken in this search is to separate the giant cell velocity signal from the other, stronger velocity signals by using a spherical harmonic representation of the spatial structures and a fourier decomposition of the temporal behavior. Any giant cell signal should be characterized by low spatial wavenumbers with temporal frequencies appropriate to the solar rotation of these patterns.

Hathaway, David H.↗

An integral analysis of transonic normal shock wave/turbulent boundary layer interactions in internal flow

An approximate integral viscous-inviscid interaction method is presented for calculating the development of a turbulent boundary layer subjected to a normal shock wave induced adverse pressure gradient in an internal axisymmetric flow. The inflow conditions and the downstream pressure are provided for the computation. In the supersonic region of shock pressure rise, the Prandtl-Meyer function is used to couple the viscous and inviscid flows. An analytical model for the coupling process is postulated and appropriate equations are defined. Downstream of the sonic point, one-dimensional inviscid flow is assumed for coupling with the viscous flow. The turbulent boundary layer is calculated using Green's integral lag-entrainment method. Comparisons of the solutions with the experimental data are made for interactions which are unseparated, near separation and separated. For comparison purposes, solutions to the time-dependent, mass-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model are also shown. The computed results from the integral method show good agreement with experimental data for unseparated interactions and reasonable agreement with the trend of the viscous effects when the interaction becomes increasingly separated.

Om, D.↗

Flow in a differentially rotated cylindrical drop at moderate Reynolds number

Galerkin finite-element approximations are combined with computer-implemented perturbation methods for tracking families of solutions to calculate the steady axisymmetric flows in a differentially rotated cylindrical drop as a function of Reynolds number Re, drop aspect ratio and the rotation ratio between the two end disks. The flows for Reynolds numbers below 100 are primarily viscous and reasonably described by an asymptotic analysis. When the disks are exactly counter-rotated, multiple steady flows are calculated that bifurcate to higher values of Re from the expected solution with two identical secondary cells stacked symmetrically about the axial midplane. The new flows have two cells of different size and are stable beyond the critical value Re sub c. The slope of the locus of Re sub c for drops with aspect ratio up to 3 disagrees with the result for two disks of infinite radius computed assuming the similarity form of the velocity field. Changing the rotation ratio for exact counter-rotation ruptures the junction of the multiple flow fields into two separated flow families.

Harriott, G. M.↗

Computer studies of baroclinic flow

Programs necessary for computing the transition curve on the regime diagram for the atmospheric general circulation experiment (AGOE) were completed and used to determine the regime diagram for the rotating annulus and some axisymmetric flows for one possible AGOE configuration. The effect of geometrical constraints on the size of eddies developing from a basic state is being examined. In AGOE, the geometric constraint should be the width of the shear zone or the baroclinic zone. Linear and nonlinear models are to be used to examine both barotropic and baroclinic flows. The results should help explain the scale selection mechanism of baroclinic eddies in the atmosphere experimental models such as AGOE, and the multiple vortex phenomenon in tornadoes.

Gall, R.↗

Selection principles and pattern formation in fluid mechanics and nonlinear shell theory

Research accomplishments are summarized and publications generated under the contract are listed. The general purpose of the research was to investigate various symmetry breaking problems in fluid mechanics by the use of structure parameters and selection principles. Although all of the nonlinear problems studied involved systems of partial differential equations, many of these problems led to the study of a single nonlinear operator equation of the form F(w, lambda, gamma) = 0, (w is an element of H), (lambda is an element of R1), (gamma is an element of R1). Instead of varying only the load parameter lambda, as is often done in the study of such equations, one of the main ideas used was to vary the structure parameter gamma in such a way that stable solutions were obtained. In this way one determines detailed stability results by making use of the structure of the model equations and the known physical parameters of the problem. The approach was carried out successfully for Benard-type convection problems, Taylor-like problems for short cylinders, rotating Couette-Poiseuille channel flows, and plane Couette flows. The main focus of the research was on wave theory of vortex breakdown in a tube. A number of preliminary results for inviscid axisymmetric flows were obtained.

Sather, Duane P.↗

Hollow conical jet models for SS 433 - A paradigm lost?

A precessing jet such as that in SS 433 may be approximated as an axisymmetric flow, if the precession time is short by comparison to the propagation time. A series of simulations has been conducted for precessing jets using an axisymmetric finite-difference hydrodynamics code. Examinations are made of hollow cylindrical jets, which lack the complication of a growing interior volume, and conical jets, which model the behavior of a precessing jet propagating on the surface of its precession cone.

Kochanek, Christopher S.↗

Computational analysis of flow in 3D propulsive transition ducts

A numerical analysis of fully three dimensional, statistically steady flows in propulsive transition ducts being considered for use in future aircraft of higher maneuverability is investigated. The purpose of the transition duct is to convert axisymmetric flow from conventional propulsion systems to that of a rectangular geometry of high aspect ratio. In an optimal design, the transition duct would be of minimal length in order to reduce the weight penalty, while the geometrical change would be gradual enough to avoid detrimental flow perturbations. Recent experiments conducted at the Propulsion Aerodynamics Branch have indicated that thrust losses in ducts of superelliptic cross-section can be surprisingly low, even if flow separation occurs near the divergent walls. In order to address the objective of developing a rational design procedure for optimal transition ducts, it is necessary to have available a reliable computational tool for the analysis of flows achieved in a sequence of configurations. Current CFD efforts involving complicated geometries usually must contend with two separate but interactive aspects: namely, grid generation and flow solution. The first two avenues of the present investigation were comprised of suitable grid generation for a class of transition ducts of superelliptic cross-section, and the subsequent application of the flow solver PAB3D to this geometry. The code, PAB3D, was developed as a comprehensive tool for the solution of both internal and external high speed flows. The third avenue of investigation has involved analytical formulations to aid in the understanding of the nature of duct flows, and also to provide a basis of comparison for subsequent numerical solutions. Numerical results to date include the generation of two preliminary grid systems for duct flows, and the initial application of PAB3D to the corresponding geometries, which are of the class tested experimentally.

Sepri, Paavo↗

Natural convection between concentric spheres

A finite-difference solution for steady natural convective flow in a concentric spherical annulus with isothermal walls has been obtained. The stream function-vorticity formulation of the equations of motion for the unsteady axisymmetric flow is used; interest lying in the final steady solution. Forward differences are used for the time derivatives and second-order central differences for the space derivatives. The alternating direction implicit method is used for solution of the discretization equations. Local one-dimensional grid adaptation is used to resolve the steep gradients in some regions of the flow at large Rayleigh numbers. The break-up into multi-cellular flow is found at high Rayleigh numbers for air and water, and at significantly low Rayleigh numbers for liquid metals. Excellent agreement with previous experimental and numerical data is obtained.

Garg, Vijay K.↗

Latitudinal Transport of Angular Momentum by Cellular Flows Observed with MDI

We have analyzed Doppler velocity images from the MDI instrument on SOHO to determine the latitudinal transport of angular momentum by the cellular photospheric flows. Doppler velocity images from 60-days in May to July of 1996 were processed to remove the p-mode oscillations, the convective blue shift, the axisymmetric flows, and any instrumental artifacts. The remaining cellular flows were examined for evidence of latitudinal angular momentum transport. Small cells show no evidence of any such transport. Cells the size of supergranules (30,000 km in diameter) show strong evidence for a poleward transport of angular momentum. This would be expected if supergranules are influenced by the Coriolis force, and if the cells are elongated in an east-west direction. We find good evidence for just such an east-west elongation of the supergranules. This elongation may be the result of differential rotation shearing the cellular structures. Data simulations of this effect support the conclusion that elongated supergranules transport angular momentum from the equator toward the poles, Cells somewhat larger than supergranules do not show evidence for this poleward transport. Further analysis of the data is planned to determine if the direction of angular momentum transport reverses for even larger cellular structures. The Sun's rapidly rotating equator must be maintained by such transport somewhere within the convection zone.

Hathaway, David H.↗

Numerical investigations in the backflow region of a vacuum plume

Four tasks were completed in this period and results were published in AIAA papers. First, a Boltzmann-2D code, was developed and applied to compute MSFC-A2 nozzle/plume flow field. It solved the two-dimensional Boltzmann-BGK equation using the Finite Difference Discrete Ordinate (FDDO) numerical technique. The code was validated by experimental data for one-dimensional shock structure predictions, paper 95-2056. Successful results for nozzle/plume flow simulation using the developed Boltzmann-2D code were presented at the 1995 AIAA Aerospace Science Conference, paper 95-0627. Second, a computer code solving two-dimensional Burnett equations was developed and applied to low-density nozzle flow field calculation. Results were also published at the 1994 AIAA Thermophysics Conference, paper 94-2055. Third, the developed two-dimensional Burnett code was extended to compute axisymmetric flow field inside MSFC-A2 nozzle, paper 95-2008. The computed nozzle exit conditions are used as input data for Direct Simulation Monte Carlo (DSMC) plume calculation. Fourth, a DSMC code was modified to compute the exhausted plume near the nozzle exit and in the backflow region.

Liaw, Goang-Shin↗