Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two dimensional 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 901 records · Page 50

Numerical analysis of the scramjet-inlet flow field by using two-dimensional Navier-Stokes equations

A computer code was developed to solve the full two dimensional Navier-Stokes equations in a supersonic combustion ramjet (scramjet) inlet. In order to be able to consider a general inlet geometry with embedded bodies, a numerical coordinate transformation is used which generates a set of boundary-fitted curvilinear coordinates. The explicit finite difference algorithm of MacCormack is used to solve the governing equations. An algebraic, two-layer eddy-viscosity model is used for the turbulent flow. The code can analyze both inviscid and viscous flows with no strut, one strut, or multiple struts in the flow field. The application of the two dimensional analysis in the preliminary parametric design studies of a scramjet inlet is discussed. Detailed results are presented for one model problem and for several actual scramjet-inlet configurations.

Kumar, A.↗

On the prediction of swirling flowfields found in axisymmetric combustor geometries

The paper reports research restricted to steady turbulence flow in axisymmetric geometries under low speed and nonreacting conditions. Numerical computations are performed for a basic two-dimensional axisymmetrical flow field similar to that found in a conventional gas turbine combustor. Calculations include a stairstep boundary representation of the expansion flow, a conventional k-epsilon turbulence model and realistic accomodation of swirl effects. A preliminary evaluation of the accuracy of computed flowfields is accomplished by comparisons with flow visualizations using neutrally-buoyant helium-filled soap bubbles as tracer particles. Comparisons of calculated results show good agreement, and it is found that a problem in swirling flows is the accuracy with which the sizes and shapes of the recirculation zones may be predicted, which may be attributed to the quality of the turbulence model.

Rhode, D. L.↗

SSME fuelside preburner two-dimensional analysis

The flow field within the fuelside preburner of the Space Shuttle Main Engine is calculated using a reacting flow code (REACT2D). Inlet and modeling parameters involved in the numerical calculation are systematically varied to establish the sensitivity of the calculated exit temperature profile. It is found that differences in the inlet equivalence ratio have a large effect on the turbine inlet temperature profile. A variety of preburner inlet modeling changes such as inlet turbulence level, modeling the gases as burned, unburned, premixed, or unmixed, are shown to have a smaller effect on the calculated turbine inlet temperature profile. Also, the form of finite differencing used is shown to have an effect on the temperature profile.

Vanoverbeke, T. J.↗

Aerodynamics of engine-airframe interaction

The development of the variational approach for the solution of inviscid aerodynamic problems using solution adaptive grids is discussed. The formulation of a new, directional weighting, functional has been shown to have desirable properties. The scheme has been applied to compute the transonic flow past two-dimensional airfoils using the Euler equations of inviscid, compressible flow. Transonic flows in quasi-one-dimensional nozzles and over the two dimensional airfoils are solved on the various solution-adaptive-grids to demonstrate the applicability of the proposed directional-concentration functional and the grid adaptation process from the stand point of improving the solution accuracy and demonstrating the overall convergence.

Caughey, David A.↗

Coherent structures in transitional and turbulent free shear flows

The development of a quantitative understanding of large-scale coherent structures in shear flows has led to the recognition of their evolutionary features and 'rules' of nonlinear interaction among each other. These determinations, in conjunction with the study of fine-grained turbulence and mean motion, constitute a general view of hydrodynamical instabilities that has resulted in useful concepts for the achievement of shear-flow control. The dynamical role of longitudinal structures in free shear flows is presently exemplified by the axisymmetric and helical modes in a round jet, which resemble the two-dimensional and spanwise-periodic three-dimensional modes in a two-dimensional mean flow.

Liu, J. T. C.↗

A theoretical evaluation of aluminum gel propellant two-phase flow losses on vehicle performance

A one-dimensional model of a hydrocarbon/Al/O2(gaseous) fueled rocket combustion chamber was developed to study secondary atomization effects on propellant combustion. This chamber model was coupled with a two dimensional, two-phase flow nozzle code to estimate the two-phase flow losses associated with solid combustion products. Results indicate that moderate secondary atomization significantly reduces propellant burnout distance and Al2O3 particle size; however, secondary atomization provides only moderate decreases in two-phase flow induced I(sub sp) losses. Despite these two-phase flow losses, a simple mission study indicates that aluminum gel propellants may permit a greater maximum payload than the hydrocarbon/O2 bi-propellant combination for a vehicle of fixed propellant volume. Secondary atomization was also found to reduce radiation losses from the solid combustion products to the chamber walls, primarily through reductions in propellant burnout distance.

Mueller, Donn C.↗

Numerical simulation of flow through inlets/diffusers with terminal shocks

Two-dimensional Euler and Navier-Stokes solutions of the flow through three inlet/diffuser configurations with terminal shock systems are reported. Calculations without bleed indicate that the terminal shock location is very sensitive to the outflow back pressure. For cases where there are no available experimental results, it becomes difficult to estimate the back pressure that will result in a terminal shock. Estimates based on quasi-one-dimensional analysis are not found adequate for complex two-dimensional flows. It is found that since the flow downstream of the terminal shock is subsonic, and what happens at the outflow boundary affects the flow inside the inlet, enough of the subsonic diffuser must be modeled to accurately predict the terminal shock region. The diffuser portion should be fairly long with the outflow boundary occurring in a region of more or less uniform flow to be able to prescribe a uniform back pressure. The third configuration studied was investigated with and without incorporating bleed in the code. It is found that the use of bleed stabilizes the shock location and allows solutions which without bleed result in unstarting of the inlet. Comparisons are made with available experimental data.

Talcott, N. A., Jr.↗

Continuation in a parameter - Experience with viscous and free surface flows

The results of modifications in continuation methods applied to obtain solutions to the Navier-Stokes systems of equations for incompressible, two-dimensional, steady flows are reported. It is shown that parameter continuation permits prediction of accurate, initial estimates for iterative processing of nonlinear finite difference and finite element equations of motions. The new parameter steps are derived from values of the preceding parameter steps. The accuracy of the estimates is ensured with appropriate choices of the step size. The continuation predictor/iterative corrector is demonstrated to trace the branches of parameter space along which steady flow states are found, and techniques are available for tracing multiply branching paths. The techniques are applied to solving the Navier-Stokes equations for flow through a rotating square channel, the formation of a falling liquid curtain, and gyrostatic equilibria of rotating cylindrical drops.

Kheshgi, H. S.↗

Viscous flow in simple curved gaps. I - An asymptotic theory. II - Viscous stress and shape function

The present asymptotic theory for generalized incompressible two-dimensional steady flow in curved channels has been constructed in the limit when gas thickness approaches zero with its lateral dimensions fixed; successive asymptotic solution terms are analytically generated by quadratures. In the second part of this work, the curvature of the gap treated is arbitrary. It is established that each term in the series solution of velocity and pressure is the product of a scale factor and a universal shape functions. Various interaction modes between the volume rate-of-flow, curvature, and its variations, are identified and quantitatively characterized.

Fan, D.-N.↗

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗

Dependence of Enstrophy Transport and Mixed Mass on Dimensionality and Initial Conditions in the Richtmyer–Meshkov Instability Induced Flows

Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.

42 ENGINEERING↗

Thermonuclear runaways in nova outbursts. 2: Effect of strong, instantaneous, local fluctuations

In an attempt to understand the manner in which nova outbursts are initiated on the surface of a white dwarf, we investigate the effects fluctuations have on the evolution of a thermonuclear runaway. Fluctuations in temperature density, or the composition of material in the burning shell may arise due to the chaotic flow field generated by convection when it occurs, or by the accretion process itself. With the aid of two-dimensional reactive flow calculations, we consider cases where a strong fluctutation in temperature arises during the early, quiescent accretion phase or during the later, more dynamic, explosion phase. In all cases we find that an instantaneous, local temperature fluctuation causes the affected material to become Rayleigh-Taylor unstable. The rapid rise and subsequent expansion of matter immediately cools the hot blob, which prevents the lateral propagation of burning. This suggests that local temperature fluctuations do not play a significant role in directly initiating the runaway, especially during the early stages. However, they may provide an efficient mechanism of mixing core material into the envelope (thereby pre-enriching the fuel for subsequent episodes of explosive hydrogen burning) and of mixing substantial amounts of the radioactive nucleus N-13 into the surface layers, making novae potential gamma-ray sources. This suggests that it is the global not the local, evolution of the core-envelope interface to high temperatures which dominates the development of the runaway. We also present a possible new scenario for the initiation of nova outbursts based on our results.

Shankar, Anurag↗

Numerical investigations in three-dimensional internal flows

The present study is a preliminary investigation into the behavior of the flow within a 28 degree total geometric turning angle hypothetical Mach 10 inlet as calculated with the full three-dimensional Navier-Stokes equations. Comparison between the two-dimensional and three-dimensional solutions have been made. The overall compression is not significantly different between the two-dimensional and center plane three dimensional solutions. Approximately one-half to two-thirds of the inlet flow at the exit of the inlet behave nominally two-dimensionally. On the other hand, flow field non-uniformities in the three-dimensional solution indicate the potential significance of the sidewall boundary layer flows ingested into the inlet. The tailoring of the geometry at the inlet shoulder and on the cowl obtained in the two-dimensional parametric design study have also proved to be effective at controlling the boundary layer behavior in the three-dimensional code. The three-dimensional inlet solution remained started indicating that the two-dimensional design had a sufficient margin to allow for three-dimensional flow field effects. Although confidence is being gained in the use of SCRAM3D (three-dimensional full Navier-Stokes code) as applied to similar flow fields, the actual effects of the three-dimensional flow fields associated with sidewalls and wind tunnel installations can require verification with ground-based experiments.

Rose, William C.↗

Recent computation of viscous effects in transonic flow

The paper describes some improvements to a method of solving the time-dependent Reynolds-averaged Navier-Stokes equations for two-dimensional compressible flows for describing the coupling between viscous and inviscid regions in transonic flows. Mesh patching is implemented at the sublayer-boundary layer interfaces and the boundary layer-inviscid flow interface. Mesh stretching is geometric and is small in regions of large gradients and large in regions of small gradients, thus preserving second-order accuracy. Mesh lines are arranged so as to accommodate MacCormack's rapid solver for hyperbolic systems and a wider class of turbulent transport models. Results of computations of transonic flows over lifting airfoils are compared with experimental data and other solutions.

Deiwert, G. S.↗

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗

Problems in fluid dynamics

A scheme was developed for the parametric differentiation and integration of gas dynamics equations. A numerical integration of the gas dynamics equations is necessarily performed for a specific set of parameter values. The linear variational equations are obtained by differentiating the exact equations with respect to each of the relevant parameters. The resulting matrix of flow quantities is referred to as the Jacobi matrix. The subsequent procedure is then straightforward. The method was tested for two dimensional supersonic flow past an airfoil, with airfoil thickness, camber, and angle of attack varied. This approach has great potential value for rapidly assessing the effect of design changes. The other focus of the work was on problems in fluid stability, bifurcations, and turbulence.

Sirovich, Lawrence↗

The measurement of boundary layers on a compressor blade in cascade. I - A unique experimental facility

A unique cascade facility is described which makes it possible to use LDV for measuring blade boundary layer profiles. Because of the need for a laser access window, the facility cannot rely on continuous blade pack suction to achieve stable flow. Instead, a combination of tailboards and strong suction upstream of the blade pack is used to create a two-dimensional periodic flow field at a chord Reynolds number of 500,000. The distribution of the upstream suction is controlled through a complex baffling system. Details of the inlet and outlet flow are documented using five-hole pressure probes, blade static-pressure measurements, and sublimation flow visualization. Experimental data are presented on the flow field on the suction surface of the blade, the pressure surface, and near-wake flow fields.

Deutsch, Steven↗

Numerical magetohydrodynamics in astronphysics: Algorithm and tests for one-dimensional flow`

We describe a numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference scheme on an Eulerian grid, called the total variation diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. We also describe a non-linear Riemann solver for ideal MHD, which includes rarefractions as well as shocks. The numerical code and the Riemann solver have been used to test each other. Extensive tests encompassing all the possible ideal MHD structures with planar symmetries (i.e., one-dimensional flows) are presented. These include those for which the field structure is two dimensional (i.e., those flows often called '1 + 1/2 dimensional') as well as those for which the magnetic field plane rotates (i.e., those flows often called '1 + 1/2 + 1/2 dimensional'). Results indicate that the code can resolve strong fast, slow, and magnetosonic shocks within two to four cells, but more cells are required if shocks become weak. With proper steepening, we could resolve rotational discontinuities within three to five cells. However, without successful implementation of steepening, contact discontinuities are resolved with approximately 10 cells and tangential discountinuities are resolved with approximately 15 cells. Out tests confirm that slow compound structures with tow-dimensional magnetic fields are composed of intermediate shocks (so-called 2-4 intermediate shocks) followed by slow rarefaction waves. Finally, tests demostrate that in two-dimensional magnetohydrodynamics, fast compound structures, which are composed of intermediate shocks (so-called 1-3 intermediate shocks) preceeded by fast rarefaction waves, are also possible.

Ryu, Dongsu↗