Parabolized Navier-Stokes analysis of three-dimensional supersonic and subsonic jet mixing problems
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg., (2) Re =1000, alpha = 15 deg., (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).
The results of an investigation of the effects of far field boundary conditions on the solution of the three dimensional Euler equations governing the flow field of a high speed single rotation propeller are presented. The results show that the solutions obtained with the nonreflecting boundary conditions are in good agreement with experimental data. The specification of nonreflecting boundary conditions is effective in reducing the dependence of the solution on the location of the far field boundary. Details of the flow field within the blade passage and the tip vortex are presented. The dependence of the computed power coefficient on the blade passage and the tip vortex are presented. The dependence of the computed power coefficient on the blade setting angle is examined.
Chapman chemistry has been used in a two-dimensional model to simulate ozone balance phenomenology. The similarity between regions of ozone production and loss calculated using Chapman chemistry and those computed using LIMS and SAMS data with a photochemical equilibrium model indicate that such simplified chemistry is useful in studying gross features in stratospheric ozone balance. Net ozone production or loss rates are brought about by departures from the photochemical equilibrium (PCE) condition. If transport drives ozone above its PCE condition, then photochemical loss dominates production. If transport drives ozone below its PCE condition, then photochemical production dominates loss. Gross features of ozone loss/production (L/P) inferred for the real atmosphere from data are also simulated using only eddy diffusion. This indicates that one must be careful in assigning a transport scheme for a two-dimensional model that mimics only behavior of the observed ozone L/P.
An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies subjected to mixed mode loading is presented. The total and product integrals consist of the sum of an area or domain integral and line integrals on the crack faces. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all the problems analyzed.
Certain Z-mode wave emissions in the earth's magnetosphere have been identified using the wave spectra and polarization measurements of the DE 1 satellite. Although such emissions accompany the aurora, and thus presumably originate from the evening-sector auroral zone, they are found to occur over much wider ranges of latitude and longitude. Since the predicted cyclotron maser emission at the cyclotron frequency could not have produced waves which travel such great distances, as shown by three-dimensional ray tracing, it is proposed instead that these emissions must originate from lower altitudes within the auroral zone and probably from near the plasma frequency inside the auroral plasma cavity.
Explore the source record for details and available documents.
The present paper describes a transient translating elastoplastic thermomechanical finite element model to study two-dimensional frictional rolling contact. Frictional two-dimensional contact is simulated by repeatedly translating a nonuniform thermomechanical distribution across the surface of an elastoplastic half space. The half space is represented by a two-dimensional finite element mesh with appropriate boundaries. Calculations are for an elastic-perfectly plastic material and the selected thermophysical properties are assumed to be temperature independent. The paper presents temperature variations, stress and plastic strain distributions and deformations. Residual tensile stresses are observed. The magnitude and depth of these stresses depends on (1) the temperature gradients and (2) the magnitudes of the normal and tangential tractions.
A hybrid numerical technique is presented for a characterization of the scattering and radiation properties of three-dimensional cavity arrays recessed in a ground plane. The technique combines the finite element and boundary integral methods and invokes Floquet's representation to formulate a system of equations for the fields at the apertures and those inside the cavities. The system is solved via the conjugate gradient method in conjunction with the Fast Fourier Transform (FFT) thus achieving an O(N) storage requirement. By virtue of the finite element method, the proposed technique is applicable to periodic arrays comprised of cavities having arbitrary shape and filled with inhomogeneous dielectrics. Several numerical results are presented, along with new measured data, which demonstrate the validity, efficiency, and capability of the technique.
In this paper, the boundary element method (BEM) is applied to two-dimensional acoustic ducts with sudden area changes. The mathematical basis of BEM along with the technique of handling singularities is presented in a systematic manner. The computer code is verified by comparison with the exact solutions and a number of illustrative examples are included.
This study presents a numerical method for solving the 3D Navier-Stokes equations for unsteady, viscous flow through multiple turbomachinery blade rows. The method solves the fully 3D Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms, while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flowfield of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flowfields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.
A numerical method for solving the three-dimensional, Navier-Stokes equations for unsteady, viscous flow and heat transfer through multiple turbomachinery blade rows is presented. The method solves the fully three-dimensional Navier-Stokes equations with an implicit scheme which is based on a control volume approach. A two-equation turbulence model with a low Reynolds number modification is employed. A third-order accurate upwinding scheme is used to approximate convection terms while a second order accurate central difference scheme is used for the discretization of viscous terms. A second-order accurate scheme is employed for the temporal discretization. The numerical method is applied to study the unsteady flow and heat transfer field of the High Pressure Fuel side Turbo-Pump (HPFTP) of the Space Shuttle Main Engine (SSME). The stage calculation is performed by coupling the stator and the rotor flow fields at each time step through an over-laid grid. Numerical results for the complete geometry with the vane trailing edge cutback are presented and compared with the available experimental data.
The renormalization group (RNG) method is used to study the physics of two-dimensional (2D) magnetohydrodynamic (MHD) turbulence. It is shown that, for a turbulent magnetofluid in two dimensions, no RNG transformation fixed point exists on account of the coexistence of energy transfer to small scales and mean-square magnetic flux transfer to large scales. The absence of a fixed point renders the RNG method incapable of describing the 2D MHD system. A similar conclusion is reached for 2D hydrodynamics, where enstrophy flows to small scales and energy to large scales. These analyses suggest that the applicability of the RNG method to turbulent systems is intrinsically limited, especially in the case of systems with dual-direction transfer.
Full 3D Navier-Stokes computational results are presented for compressible flows within nondiffusing and diffusing S-ducts. The present study provides an understanding of the performance characteristics of typical S-ducts with attached and separated flows and provides a frame of reference for future computational fluid dynamic studies of internal flows with strong secondary flows. The predicted results, which were obtained using both H- and O-grids, are compared with the experimental wall pressures, static and total pressure fields, and velocity vectors. In addition, computed boundary-layer thickness, velocity profiles in wall coordinates and skin friction values are presented. The inviscid contributions to the secondary flows are quantified. The S-duct entrance Mach number was 0.6, and the Reynolds number was 1.76 x 10 exp 6 based on the upstream duct diameter.
A method has been developed for making photoelastic analyses of three-dimensional stress systems by utilizing the polarization phenomena associated with the scattering of light. By this method, the maximum shear and the directions of the three principal stresses at any point within a model can be determined, and the two principal stresses at a free-bounding surface can be separately evaluated. Polarized light is projected into the model through a slit so that it illuminates a plane section. The light is continuously analyzed along its path by scattering and the state of stress in the illuminated section is obtained. By means of a series of such sections, the entire stress field may be explored. The method was used to analyze the stress system of a simple beam in bending. The results were found to be in good agreement with those expected from elementary theory.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.