Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solvers”

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 1,513 records · Page 84

CFD analysis of a hydrogen fueled ramjet engine at Mach 3.44

The full Navier-Stokes solver, RPLUS, has been used to compute the flow field of a ramjet engine in the low supersonic operating regime. Computations have been made for the case where hydrogen is transversely injected and ignited using nine species and eighteen reactions to model the combustion process. The flow fields for two methods of hydrogen injection have been computed and compared to the available experimental data. The purpose of these computations was to assess the mixing effectiveness of the injected hydrogen with the primary supersonic stream. Results indicate significant ignition delay resulting from the poor mixing dynamics in the combustor section.

Duncan, Beverly↗

Use of the PARC code for generic NASP nozzles operating at off-design transonic conditions

A three-dimensional Navier-Stokes solver, PARC3D, has been investigated for applications to the wind tunnel testing of high-speed nozzle configurations, particularly with extended expansion ramps, operating at transonic off-design conditions. Numerous cases have been run and compared with existing wind tunnel data. The code has been shown to be an efficient and robust tool in the analysis of complex nozzle configurations whose flows are substantially determined by three-dimensional effects. Although the method requires that a grid, generated independently of the code, be supplied to it, the robust nature of the code does not require that this grid be excessively detailed, nor does it seem to require that the grid be fine-tuned for individual cases. The directness and flexibility of the boundary condition specifications allow a variety of configurations and flow conditions to be handled very efficiently.

Yaros, Steven F.↗

Periodic blade loads of a high speed propeller at small angle of attack

A code has been developed to predict the periodic aerodynamic loads of an advanced turboprop propeller. The analytical formulation accounts for flow three-dimensionality and flow periodicity due to the propeller inclination. The flow past the blade sections is computed using a thin layer Navier-Stokes solver. An iterative procedure is used to account for the induced axial and rotational velocities. The viscous periodic results are obtained for an eight-bladed Hamilton Standard SR-7L advanced propeller at a cruise Mach number of 0.813 and 35,000 ft. altitude. The results are shown for flow field quantities and performance parameters during the blade passage in the plane of rotation illustrating the periodic nature of blade flow separation and shocks. The time averaged coefficients of thrust and power are computed and compared with available flight test data. The results obtained show excellent agreement at cruise conditions for small nacelle angles of attack.

Takallu, M. A.↗

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗

3-D CFD in a day - The laser digitizer project

The computation of airflow over complex configurations requires a complete description of the geometry. This can be obtained from CAD data, from blueprints, or from actual models. In any case, the time required is currently estimated at 4 to 6 months. It is proposed to shorten this time by a factor of 10 to 100 through the use of automated software, a fast, highly parallel computer and a three-dimensional laser digitizer. This device can provide (x,y,z) coordinates of surface points at rates exceeding 14,500/sec. Thus, it is possible to digitize an entire model in a few minutes. The accuracy of measurement on a flat white surface is better than 0.005 inches. Higher accuracy is available at higher cost. This work discusses the challenges which remain to be addressed. In particular, the surface point data need to be converted into a surface description, the surface description needs to be made into a surface grid, and the surface grid used to make a volume grid for the flow solver. Algorithms are kept in place or in mind for all of these problems. Integration of the more mature flow solution and visualization algorithms then allows generation of solution graphics directly from a wind tunnel model.

Merriam, Marshal↗

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗

Three-Dimensional Navier-Stokes Simulation of Space Shuttle Main Propulsion 17-inch Disconnect Valves

A steady incompressible three-dimensional viscous flow analysis has been conducted for the Space Shuttle external tank/orbiter propellant feed line disconnect flapper valves with upstream elbows. The Navier-Stokes code, INS3D, is modified to handle interior obstacles and a simple turbulence model. The flow solver is tested for stability and convergence in the presence of interior flappers. An under-relaxation scheme has been incorporated to improve the solution stability. Important flow characteristics such as secondary flows, recirculation, vortex and wake regions, and separated flows are observed. Computed values for forces, moments, and pressure drop are in satisfactory agreement with water flow test data covering a maximum tube Reynolds number of 3.5 million. The predicted hydrodynamical stability of the flappers correlates well with the measurements.

Kandula, M.↗

Iterative solution techniques in boundary element analysis

Iterative techniques for the solution of the algebraic equations associated with the direct boundary element analysis (BEA) method are discussed. Continuum structural response analysis problems are considered, employing single- and multizone boundary element models with and without zone condensation. The impact on convergence rate and computer resource requirements associated with the sparse and blocked matrices, resulting in multizone BEA, is studied. Both conjugate gradient and generalized minimum residual preconditioned iterative solvers are applied for these problems and the performance of these algorithms is reported. Included is a quantification of the impact of the preconditioning utilized to render the boundary element matrices solvable by the respective iterative methods in a time competitive with direct methods.

Kane, J. H.↗

An approximate viscous shock layer approach to calculating hypersonic flows about blunt-nosed bodies

An approximate axisymmetric method has been developed which can reliably calculate fully viscous hypersonic flows over blunt-nosed bodies. By substituting Maslen's second order pressure expression for the normal momentum equation, a simplified form of the viscous shock layer (VSL) equations is obtained. This approach can solve both the subsonic and supersonic regions of the shock layer without a starting solution for the shock shape. Since the method is fully viscous, the problems associated with coupling a boundary-layer solution with an inviscid-layer solution are avoided. This procedure is significantly faster than the parabolized Navier-Stokes (PNS) or VSL solvers and would be useful in a preliminary design environment. Problems associated with a previously developed approximate VSL technique are addressed. Surface heat transfer and pressure predictions are comparable to both VSL results and experimental data. The present technique generates its own shock shape as part of its solution, and therefore could be used to provide more accurate initial shock shapes for higher-order procedures which require starting solutions.

Cheatwood, F. MCN.↗

Two-species variable gamma formulation for a forward-firing rocket nozzle

The object of the paper is to develop an algorithm for an axisymmetric case to study the effects of variable gamma on the flowfield of a two-species mixture. This algorithm forms the basis for extension to three-dimensional SSME application where the Shuttle executes a return-to-launch-site maneuver. The PARC2DR code was used as the foundation for the Navier-Stokes solver, and the species continuity equation for the exhaust species was coupled in an explicit fashion to the Navier-Stokes equation via the variable gamma. The species continuity equation could not be solved with explicit boundary conditions; instead, implicit boundary conditions were implemented.

Saladino, Anthony J.↗

Development of a nonequilibrium radiative heating prediction method for coupled flowfield solutions

A method for predicting radiation in thermochemical nonequilibrium flows has been developed. A summary of the development is presented, with emphasis on its unique features. The optimization of the method to provide efficient solutions allowing coupling to flowfield solvers is discussed. Representative results are obtained and compared to previous solutions and flight data. The method is found to be a useful tool for studies of nonequilibrium flowfields.

Hartung, Lin C.↗

The effects of self-generated and applied magnetic fields on the computation of flow over a Mars return aerobrake

A CFD technique is developed to calculate the electromagnetic phenomena simultaneously with the fluid flow in the shock layer over an axisymmetric blunt body in a thermal-equilibrium chemical-nonequilibrium environment. The flowfield is solved using an explicit time-marching, first-order spatially accurate scheme. The electromagnetic phenomena are coupled to the real-gas flow solver through an iterative procedure. The electromagnetic terms introduce a strong stiffness, which was overcome by using significantly smaller time steps for the electromagnetic conservation equation. The technique is applied in calculating the flow over a Mars return aerobrake vehicle entering the Earth's atmosphere. For the case where no external field is applied, the electromagnetic effects have little impact on the flowfield.

Palmer, Grant↗

Assessment of meanflow solutions for instability analysis of transitioning flows

Meanflow solutions of 3D supersonic flow past a cone at incidence and a swept leading edge wing have been obtained by thre methods, viz., boundary-layer, parabolized Navier-Stokes, and thin shear-layer Navier-Stokes solvers. The smoothness and accuracy of the solution profiles are compared with a view to applying the meanflow solution to boundary-layer stability analysis.

Iyer, Venkit↗

Interactive design of hypersonic waverider geometries

The paper deals with an inverse design code utilizing the method of oscillating cones; the code integrated into an interactive graphics software package allows manipulation of both the exit-plane shock profile and leading edge of the vehicle. Another interactive feature of the system is the ability to vary freestream conditions and reevaluate the governing conditions. The development of the oscillating cones is shown on five classes each of which is chosen to demonstrate an aspect of improved design flexibility over previous studies. Results are evaluated using a robust flow solver, insuring that the shock shapes specified in the design process are recovered. It is pointed out that the expanded range of waverider geometries that may be generated using the oscillating cones technique may provide insight into visually oriented optimization parameters such as volumetric efficiency and practical planform.

Center, K. B.↗

Flowfield computations and comparison with Shuttle aerodynamic data

An in-house developed flow solver, E3D, has been applied to investigate the flow field around the Shuttle Orbiter over an angle-of-attack range of 0 deg to 60 deg at Mach 3.5 and 10.0. The 3D Euler equations are integrated by means of a time-marching finite-volume shock-capturing method, based on cell-centered and upwind evaluation of the cell face fluxes with van Leer or Roe's flux-splitting techniques. A comparison of CFD, wind tunnel, and aeronautical data base is presented. Solutions are in very good agreement with available experimental data. Through the comparisons, it is demonstrated that the E3D code is capable of evaluating total aerodynamic performance.

Ma, E. C.↗

Low thrust viscous nozzle flow field predictions

The low thrust viscous nozzle flows have been studied numerically by using the steady-state Reynolds-averaged Navier-Stokes equations in conservative form. A new Lower-Upper factorization scheme in the axisymmetric coordinate system has been developed and is analyzed in this paper. This new Navier-Stokes solver has shown that the numerical analog is accurate and stable for the convergent and divergent nozzle flowfield calculation. The range of the chamber pressure is (1.85-29.4)psi with the chamber temperature equal to 600 F. The results are compared with the data collected in the NASA/MSFC vacuum chamber, and they are in very reasonable agreement.

Chou, L. C.↗

Prediction of unsteady aerodynamic loads in cascades using the linearized Euler equations on deforming grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.

Hall, Kenneth C.↗

A study of hydrogen diffusion flames using PDF turbulence model

The application of probability density function (pdf) turbulence models is addressed in this work. For the purpose of accurate prediction of turbulent combustion, an algorithm that combines a conventional CFD flow solver with the Monte Carlo simulation of the pdf evolution equation has been developed. The algorithm has been validated using experimental data for a heated turbulent plane jet. The study of H2-F2 diffusion flames has been carried out using this algorithm. Numerical results compared favorably with experimental data. The computuations show that the flame center shifts as the equivalence ratio changes, and that for the same equivalence ratio, similarity solutions for flames exist.

Hsu, Andrew T.↗