SEARCH · Engineering Papers

Results for “Algebraic models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Navier-Stokes calculation of solid-propellant rocket motor internal flowfields

A comprehensive numerical analysis has been carried out to study the detailed physical and chemical processes involved in the combustion of homogeneous propellant in a rocket motor. The formulation is based on the time-dependent full Navier-Stokes equations, with special attention devoted to the chemical reactions in both gas and condensed phases. The turbulence closure is achieved using both the Baldwin-Lomax algebraic model and a modified k-epsilon two-equation scheme with a low Reynolds number and near-wall treatment. The effects of variable thermodynamic and transport properties are also included. The system of governing equations are solved using a multi-stage Runge-Kutta shceme with the source terms treated implicitly. Preliminary results clearly demonstrate the presence of various combustion regimes in the vicinity of propellant surface. The effects of propellant combustion on the motor internal flowfields are investigated in detail.

Hsieh, Kwang-Chung↗

NASA NTRS1988

Navier-Stokes calculations of scramjet-afterbody flowfields

Scramjet simulant-gas exhaust flows are computed using two-dimensional Navier-Stokes (NS) equations. The NS equations are solved using an implicit, upwind, finite-volume scheme, and the Reynolds stresses are modeled algebraically. The solutions are compared with experimental data where the freestream and the nozzle exhaust fluids are air. The NS equations for the flow, where the scramjet exhaust simulant gas is a mixture of Freon-12 and argon, are solved by an explicit, finite difference scheme. The results are compared with experimental data, and certain discrepancies are explained.

Baysal, O.↗

NASA NTRS1990

Numerical simulation of an F-16A at angle of attack

The transonic flowfield around an F-16A fighter configuration at a moderate incidence angle is simulated by solving the Navier-Stokes equations on a single-block grid. The numerical solution matches experimental freestream conditions with a mach number of 0.85, 16 degrees angle of attack, and a characteristic Reynolds number of 12.75 million. MacCormack's explicit algorithm is used in conjuction with a local time step and consecutive mesh refinement procedure to accelerate numerical convergence. The Baldwin-Lomax algebraic model provides turbulent closure. Computed surface pressure distributions and the aircraft lift coefficient compare favorably with wind tunnel data. The drag coefficient in the simulation overpredicts the experimental value by 8 percent.

Huband, G. W.↗

NASA NTRS1991

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q-omega low Reynolds number two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Gorla, Rama S. R.↗

NASA NTRS1991

An efficient solution technique for shockwave-boundary layer interactions with flow separation and slot suction effects

An efficient method for computing two-dimensional compressible Navier-Stokes flow fields is presented. The solution algorithm is a fully-implicit approximate factorization technique based on an unsymmetric line Gauss-Seidel splitting of the equation system Jacobian matrix. Convergence characteristics are improved by the addition of acceleration techniques based on Shamanskii's method for nonlinear equations and Broyden's quasi-Newton update. Characteristic-based differencing of the equations is provided by means of Van Leer's flux vector splitting. In this investigation, emphasis is placed on the fast and accurate computation of shock-wave-boundary layer interactions with and without slot suction effects. In the latter context, a set of numerical boundary conditions for simulating the transpiration flow in an open slot is devised. Both laminar and turbulent cases are considered, with turbulent closure provided by a modified Cebeci-Smith algebraic model. Comparisons with computational and experimental data sets are presented for a variety of interactions, and a fully-coupled simulation of a plenum chamber/inlet flowfield with shock interaction and suction is also shown and discussed.

Edwards, Jack R.↗

NASA NTRS1992

Overview of the relevant CFD work at Thiokol Corporation

An in-house developed proprietary advanced computational fluid dynamics code called SHARP (Trademark) is a primary tool for many flow simulations and design analyses. The SHARP code is a time dependent, two dimensional (2-D) axisymmetric numerical solution technique for the compressible Navier-Stokes equations. The solution technique in SHARP uses a vectorizable implicit, second order accurate in time and space, finite volume scheme based on an upwind flux-difference splitting of a Roe-type approximated Riemann solver, Van Leer's flux vector splitting, and a fourth order artificial dissipation scheme with a preconditioning to accelerate the flow solution. Turbulence is simulated by an algebraic model, and ultimately the kappa-epsilon model. Some other capabilities of the code are 2-D two-phase Lagrangian particle tracking and cell blockages. Extensive development and testing has been conducted on the 3-D version of the code with flow, combustion, and turbulence interactions. The emphasis here is on the specific applications of SHARP in Solid Rocket Motor design. Information is given in viewgraph form.

Chwalowski, Pawel↗

NASA NTRS1992

Navier-Stokes analysis of turbomachinery blade external heat transfer

The two-dimensional, compressible, thin-layer Navier-Stokes and energy equations were solved numerically to obtain heat transfer rates on turbomachinery blades. The Baldwin-Lomax algebraic model and the q - omega low Reynolds number, two-equation model were used for modeling of turbulence. For the numerical solution of the governing equations a four-stage Runge-Kutta solver was employed. The turbulence model equations were solved using an implicit scheme. Numerical solutions are presented for two-dimensional flow within two vane cascades. The heat transfer results and the pressure distributions were compared with published experimental data. The agreement between the numerical calculations and the experimental values were found to be generally favorable. The position of transition from laminar to turbulent flow was also predicted accurately.

Ameri, A. A.↗

NASA NTRS1992

Computation of vortex wake flows and control of their effects on trailing wings

The near-vortex-wake flow of a large aspect-ratio rectangular wing is accurately computed by using the thin-layer and full Navier-Stokes (NS) equations. The chordwise section of the wing is a NACA-0012 airfoil and its tip is round. The computations have been carried out on a fine C-O grid using an implicit, upwind, flux-difference splitting, finite-volume scheme. The thin-layer NS results have been obtained with and without flux limiters, and the full NS results have been obtained without flux limiters. Flow transition from laminar to turbulent is mimicked by turning-on the Baldwin-Lomax algebraic model at an experimentally prescribed chord-station location of 0.05. Comparison of computed results and experimental data shows that the full NS results give the best resolution of the near-vortex-wake flow. Next, the strength of the wing-tip vortex has been reduced substantially without reducing the lift coefficient by using flow-injection from a slot along a portion of the wing tip. The flow injection is directed in the wing plane at 45 deg with the wing-tip chord.

Wong, Tin-Chee↗

NASA NTRS1987

Research on complex turbulent flows at the UW (overall effort)

Development of a pressure-strain model, an algebraic stress model, and wall functions appropriate for flows with spanwise variations in the local wall shear stress are accomplished. Furthermore, a hot-wire measurement technique was also developed for determining the local mean velocity and Reynolds stresses in a complex flow. Experiments were performed on supersonic and subsonic turbulent flow in a square duct, flow about a strut-endwall, flow within a transition duct, and on co-flowing annular jets with swirl. All results are presented in a viewgraph format.

Gessner, F. B.↗

NASA NTRS1993

Turbulence modelling of flow fields in thrust chambers

Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.

Chen, C. P.↗

NASA NTRS1992

Variable-complexity aerodynamic-structural design of a high-speed civil transport wing

A variable-complexity strategy of combining simple and detailed analysis methods is presented for the design optimization of a high-speed civil transport (HSCT) wing. Two sets of results are shown: the aerodynamic design of the wing using algebraic weight equations for structural considerations, and optimization results of the internal wing structure for a fixed wing configuration. We show example results indicating that using simple analysis methods alone for the calculation of a critical constraint can allow an optimizer to exploit weaknesses in the analysis. The structural optimization results provide a valuable check for the weight equations used in the aerodynamic design. In addition, these results confirm the need for using simple, algebraic models in conjunction with more detailed analysis methods. A strategy of interlaced aerodynanic-structural design is proposed.

Hutchison, M. G.↗

NASA NTRS1995

Heat transfer predictions for two turbine nozzle geometries at high Reynolds and Mach numbers

Predictions of turbine vane and endwall heat transfer and pressure distributions are compared with experimental measurements for two vane geometries. The differences in geometries were due to differences in the hub profile, and both geometries were derived from the design of a high rim speed turbine (HRST). The experiments were conducted in the Isentropic Light Piston Facility (ILPF) at Pyestock at a Reynolds number of 5.3 x 10(exp 6), a Mach number of 1.2, and a wall-to-gas temperature ratio of 0.66. Predictions are given for two different steady-state three-dimensional Navier-Stokes computational analyses. C-type meshes were used, and algebraic models were employed to calculate the turbulent eddy viscosity. The effects of different turbulence modeling assumptions on the predicted results are examined. Comparisons are also given between predicted and measured total pressure distributions behind the vane. The combination of realistic engine geometries and flow conditions proved to be quite demanding in terms of the convergence of the CFD solutions. An appropriate method of grid generation, which resulted in consistently converged CFD solutions, was identified.

Boyle, R. J.↗

NASA NTRS1996

Application of Hybrid Optimization-Expert System for Optimal Power Management on Board Space Power Station

The space power system has two sources of energy: photo-voltaic blankets and batteries. The optimal power management problem on-board has two broad operations: off-line power scheduling to determine the load allocation schedule of the next several hours based on the forecast of load and solar power availability. The nature of this study puts less emphasis on speed requirement for computation and more importance on the optimality of the solution. The second category problem, on-line power rescheduling, is needed in the event of occurrence of a contingency to optimally reschedule the loads to minimize the 'unused' or 'wasted' energy while keeping the priority on certain type of load and minimum disturbance of the original optimal schedule determined in the first-stage off-line study. The computational performance of the on-line 'rescheduler' is an important criterion and plays a critical role in the selection of the appropriate tool. The Howard University Center for Energy Systems and Control has developed a hybrid optimization-expert systems based power management program. The pre-scheduler has been developed using a non-linear multi-objective optimization technique called the Outer Approximation method and implemented using the General Algebraic Modeling System (GAMS). The optimization model has the capability of dealing with multiple conflicting objectives viz. maximizing energy utilization, minimizing the variation of load over a day, etc. and incorporates several complex interaction between the loads in a space system. The rescheduling is performed using an expert system developed in PROLOG which utilizes a rule-base for reallocation of the loads in an emergency condition viz. shortage of power due to solar array failure, increase of base load, addition of new activity, repetition of old activity etc. Both the modules handle decision making on battery charging and discharging and allocation of loads over a time-horizon of a day divided into intervals of 10 minutes. The models have been extensively tested using a case study for the Space Station Freedom and the results for the case study will be presented. Several future enhancements of the pre-scheduler and the 'rescheduler' have been outlined which include graphic analyzer for the on-line module, incorporating probabilistic considerations, including spatial location of the loads and the connectivity using a direct current (DC) load flow model.

Momoh, James↗

NASA NTRS1981

Numerical Simulation of Three-Dimensional Boattail Afterbody Flowfields

The thin shear-layer approximations of the three-dimensional, compressible Navier-Stokes equations are solved for subsonic, transonic, and supersonic now over axisymmetric boattail bodies at moderate angles of attack. The plume is simulated by a solid body configuration identical to those used In experimental tests. An implicit algorithm of second-order accuracy is used to solve the equations on the ILLIAC 4 computer. The turbulence is expressed by an algebraic model applicable to three-dimensional flowfields with moderate separation. The formulation used is attractive in its independence of boundary-layer parameters. Such a simple model, however, is incapable of supporting detailed quantitative descriptions of complex shear flows. Never-the-less, good qualitative comparisons are found with three different sets of experimental date. Quantitative improvement will depend on improved turbulence transport descriptions.

Deiwert, G. S.↗

NASA NTRS1995

Analysis of Aerothermal Heating Over Reusable Single-Stage to Orbit Vehicles

Using the General Aerodynamic Simulation Program (GASP) with finite-rate surface catalysis and radiative equilibrium surface temperature conditions, a computational study has been conducted to predict three-dimensional aero-convective heating over the proposed reusable single-stage to orbit vehicles. The solutions for three basic classes of usable single-stage-to-orbit vehicles, including wing body, lifting body and vertical lander, are presented and compared to provide information for vehicle design. The free stream conditions used in this study occur at an altitude of 57.5 km, ambient density of 4.129e-4 kg/cubic m, and velocity of 4,849 m/s. The convective heating over the vehicle surface in both laminar and turbulent flow regimes were computed. The Baldwin-Lomax algebraic model with a modified van Driest damping factor was applied in the turbulent flow calculations. A bifurcation model which guarantees the conservation of mass diffusive fluxes was adopted to account for the multi-species mass diffusion. The effects of turbulence and surface catalysis on the surface heating are discussed in detail The solutions indicate the degree of oxygen dissociation over the leeward side is sensitive to the shape of vehicle body. Predicted convective heating over the entire surface of lifting body is strongly affected by the presence of surface catalysis. For wing body and vertical lander, the surface catalysis effect does not appear to be important in the far downstream region, and turbulence is the dominant effect on the heating prediction in this region.

Chen, Y. K.↗

NASA NTRS2006

General Flow-Solver Code for Turbomachinery Applications

Phantom is a computer code intended primarily for real-fluid turbomachinery problems. It is based on Corsair, an ideal-gas turbomachinery code, developed by the same authors, which evolved from the ROTOR codes from NASA Ames. Phantom is applicable to real and ideal fluids, both compressible and incompressible, flowing at subsonic, transonic, and supersonic speeds. It utilizes structured, overset, O- and H-type zonal grids to discretize flow fields and represent relative motions of components. Values on grid boundaries are updated at each time step by bilinear interpolation from adjacent grids. Inviscid fluxes are calculated to third-order spatial accuracy using Roe s scheme. Viscous fluxes are calculated using second-order-accurate central differences. The code is second-order accurate in time. Turbulence is represented by a modified Baldwin-Lomax algebraic model. The code offers two options for determining properties of fluids: One is based on equations of state, thermodynamic departure functions, and corresponding state principles. The other, which is more efficient, is based on splines generated from tables of properties of real fluids. Phantom currently contains fluid-property routines for water, hydrogen, oxygen, nitrogen, kerosene, methane, and carbon monoxide as well as ideal gases.

Dorney, Daniel↗

NASA NTRSDate not supplied

Unsteady CFD Simulations of a Compression Corner Geometry Using Wall-Modeled LES Methods in Loci/CHEM

Several wall-modeled large eddy simulation (WMLES) methods are tested by simulating an unsteady Mach 2.0 compression corner geometry in Mississippi State’s Loci/CHEM solver. This study is conducted to evaluate the usage and requirements of these WMLES methods for applications regarding fluctuating pressure environments on launch vehicles with computational fluid dynamics (CFD). Two hybrid Reynolds-averaged Navier-Stokes (RANS)-large eddy simulation (LES) methods, Dynamic Hybrid RANS-LES (DHRL) and Improved Delayed Detached Eddy Simulation (IDDES), and one wall-stress-model, the Algebraic Wall Model for Wall-Modeled LES (AWMLES), are tested on varying grid and timestep refinement levels. These grid and timestep sizes are chosen to test the minimum requirements for successfully running these WMLES methods. The simulations are evaluated based on turbulent boundary layer properties in the developed boundary layer as well as unsteady quantities relating to fluctuating pressure environments in the region of the compression corner. The DHRL method shows good agreement with the comparison wind tunnel data and shows good grid and timestep convergence. The results from the IDDES and WMLES simulations show good agreement for several quantities with some discrepancies regarding others. The results presented in this paper will be used to inform further studies in predicting unsteady environments on higher-complexity geometries.

Bryson Frank↗

NASA NTRS1996

A Galilean Invariant Explicit Algebraic Reynolds Stress Model for Curved Flows

A Galilean invariant weak-equilbrium hypothesis that is sensitive to streamline curvature is proposed. The hypothesis leads to an algebraic Reynolds stress model for curved flows that is fully explicit and self-consistent. The model is tested in curved homogeneous shear flow: the agreement is excellent with Reynolds stress closure model and adequate with available experimental data.

Girimaji, Sharath↗