Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes 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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Development of the Glenn-Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

heat transfer↗
NASA NTRS1985

Numerical studies of unsteady transonic flow over oscillating airfoil

A finite difference solution to the Navier-Stokes equations combined with a time varying grid generation technique was used to compute unsteady transonic flow over an oscillating airfoil. These computations were compared with experimental data (obtained at Ames Research Center) which form part of the AGARD standard configuration for aeroelastic analysis. A variety of approximations to the full Navier-Stokes equations was used to determine the effect of frequency, shock wave motion, flow separation, and airfoil geometry on unsteady pressures and overall air loads. Good agreement is shown between experiment and theory with the limiting factor being the lack of a reliable turbulence model for high Reynolds number, unsteady transonic flows.

Chyu, W. J.↗
NASA NTRS1997

Enhanced Multiobjective Optimization Technique for Comprehensive Aerospace Design

A multidisciplinary design optimization procedure which couples formal multiobjectives based techniques and complex analysis procedures (such as computational fluid dynamics (CFD) codes) developed. The procedure has been demonstrated on a specific high speed flow application involving aerodynamics and acoustics (sonic boom minimization). In order to account for multiple design objectives arising from complex performance requirements, multiobjective formulation techniques are used to formulate the optimization problem. Techniques to enhance the existing Kreisselmeier-Steinhauser (K-S) function multiobjective formulation approach have been developed. The K-S function procedure used in the proposed work transforms a constrained multiple objective functions problem into an unconstrained problem which then is solved using the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm. Weight factors are introduced during the transformation process to each objective function. This enhanced procedure will provide the designer the capability to emphasize specific design objectives during the optimization process. The demonstration of the procedure utilizes a computational Fluid dynamics (CFD) code which solves the three-dimensional parabolized Navier-Stokes (PNS) equations for the flow field along with an appropriate sonic boom evaluation procedure thus introducing both aerodynamic performance as well as sonic boom as the design objectives to be optimized simultaneously. Sensitivity analysis is performed using a discrete differentiation approach. An approximation technique has been used within the optimizer to improve the overall computational efficiency of the procedure in order to make it suitable for design applications in an industrial setting.

Chattopadhyay, Aditi↗
NASA NTRS1985

Finite-amplitude steady waves in plane viscous shear flows

Computations of two-dimensional solutions of the Navier-Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.

Milinazzo, F. A.↗
NASA NTRS1980

Numerical simulation of three-dimensional boattail afterbody flow fields

The thin shear layer approximations of the three-dimensional, compressible Navier-Stokes equations are solved for subsonic, transonic, and supersonic flow over axisymmetric boattail bodies at moderate angles of attack. The plume is modeled 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 IV computer. The turbulence is expressed by an algebraic model applicable to three-dimensional flow fields with moderate separation. The computed results compare favorably with three different sets of experimental data reported by Reubush, Shrewsbury, and Benek, respectively

Deiwert, G. S.↗
NASA NTRS1992

A unified approach for numerical simulation of viscous compressible and incompressible flows over adiabatic and isothermal walls

A new formulation (including the choice of variables, their non-dimensionalization, and the form of the artificial viscosity) is proposed for the numerical solution of the full Navier-Stokes equations for compressible and incompressible flows with heat transfer. With the present approach, the same code can be used for constant as well as variable density flows. The changes of the density due to pressure and temperature variations are identified and it is shown that the low Mach number approximation is a special case. At zero Mach number, the density changes due to the temperature variation are accounted for, mainly through a body force term in the momentum equation. It is also shown that the Boussinesq approximation of the buoyancy effects in an incompressible flow is a special case. To demonstrate the new capability, three examples are tested. Flows in driven cavities with adiabatic and isothermal walls are simulated with the same code as well as incompressible and supersonic flows over a wall with and without a groove. Finally, viscous flow simulations of an oblique shock reflection from a flat plate are shown to be in good agreement with the solutions available in literature.

Hafez, M.↗
NASA NTRS1983

Numerical studies of laminar and turbulent drag reduction, part 2

The flow over wave shaped surfaces is studied using a Navier Stokes solver. Detailed comparisons with theoretical results are presented, including the stability of a laminar flow over wavy surfaces. Drag characteristics of nonplanar surfaces are predicted using the Navier-Stokes solver. The secondary instabilities of wall bounded and free shear flows are also discussed.

Balasubramanian, R.↗
NASA NTRS1988

Large-scale turbulence simulation on the Navier-Stokes Computer

Analytical techniques for turbulent air or water flow over vehicles and for the development of active control methods are described, with a focus on simulations being implemented on the Navier-Stokes Computer (NSC) at Princeton University. Consideration is given to the general NSC architecture (the central control, memory, and processing units; the switching network; and the prototype NSC), the algorithm for large-scale MHD turbulence (Orszag and Tang, 1979), and the pipeline configurations. Extensive diagrams and flow charts are provided.

Hayder, M. E.↗
NASA NTRS1975

Development of an efficient computer code to solve the time-dependent Navier-Stokes equations

A research effort was conducted with the goal of reducing computer time of a Navier Stokes Computer Code for prediction of viscous flow fields about lifting bodies. A two-dimensional, time-dependent, laminar, transonic computer code (STOKES) was modified to incorporate a non-uniform timestep procedure. The non-uniform time-step requires updating of a zone only as often as required by its own stability criteria or that of its immediate neighbors. In the uniform timestep scheme each zone is updated as often as required by the least stable zone of the finite difference mesh. Because of less frequent update of program variables it was expected that the nonuniform timestep would result in a reduction of execution time by a factor of five to ten. Available funding was exhausted prior to successful demonstration of the benefits to be derived from the non-uniform time-step method.

Harp, J. L., Jr.↗
NASA NTRS2012

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗
NASA NTRS1978

Equations of motion and two-equation turbulence model for plane or axisymmetric turbulent flows in body-oriented orthogonal curvilinear coordinates and mass-averaged dependent variables

The full Navier-Stokes time-dependent, compressible, turbulent, mean-flow equations in mass-averaged variables for plane or axisymmetric flow are presented. The equations are derived in a body-oriented, orthogonal, curvilinear coordinate system. Turbulence is modelled by a system of two equations for mass-averaged turbulent kinetic energy and dissipation rate proposed. These equations are rederived and some new features are discussed. A system of second order boundary layer equations is then derived which includes the effects of longitudinal curvature and the normal pressure gradient. The Wilcox and Chambers approach is used in considering effects of streamline curvature on turbulence phenomena in turbulent boundary layer type flows. Their two-equation turbulence model with curvature terms are rederived for the cases considered in the present report. The derived system equations serves as a basis for an investigation of problems where streamline curvature is of the order of the characteristic length in the longitudinal direction.

Sislian, J. P.↗
NASA NTRS1991

CFD validation of a supersonic laminar flow control concept

A three-dimensional Navier-Stokes code is used in conjunction with a linear compressible stability analysis code to develop a numerical procedure for prediction of laminar flow transition. The procedure is applied to a modified F-16XL fighter with a laminar flow control glove at supersonic speed. Details of boundary layer stability analysis indicate that, computationally, laminar flow could be realized on the highly swept wing in the absence of the leading edge attachment-line contamination. Effects of the three-dimensionality of the flow were shown to be important in the boundary layer stability analysis. The numerically predicted surface pressures compare favorably with the flight test data.

Woan, Chung-Jin↗
NASA NTRS1989

Numerical investigation of hypersonic exhaust plume/afterbody flow fields

An upwind, implicit Navier-Stokes computer program has been applied to hypersonic exhaust plume/afterbody flowfields. The sensitivity of gross thrust to operating conditions has been assessed through parametric variations. Comparison of the numerical results with available experimental data shows good agreement in all cases investigated. Results show that, for moderately underexpanded jets, the afterbody force varies linearly with the nozzle exit pressure. Exhaust gases with low isentropic exponents (gamma) were found to contribute up to 25 percent more afterbody force than high-gamma exhaust gases. Modifying the nozzle geometry influenced the exhaust plume development, which had a significant effect on the afterbody force. Grid density, while important to resolving the initial plume/afterbody interaction, had only a minor impact on the resultant afterbody force.

Edwards, T. A.↗
NASA NTRS1992

Segmented multigrid domain decomposition procedure for incompressible viscous flows

In this paper, the Navier-Stokes (NS) equations are approximated with a reduced NS system, that represents the lowest-order terms in an asymptotic Re expansion. This system allows for simplified boundary conditions, more generality in the location of the outflow boundary, and insures mass conservation in all subdomain grid interfaces, and also at the outflow boundary. The higher-order NS diffusion terms are included through a deferred corrector, in chosen subdomains, when required.

Srinivasan, Kumar↗
NASA NTRS1997

Effective Control of Computationally Simulated Wing Rock in Subsonic Flow

The unsteady compressible, full Navier-Stokes (NS) equations and the Euler equations of rigid-body dynamics are sequentially solved to simulate the delta wing rock phenomenon. The NS equations are solved time accurately, using the implicit, upwind, Roe flux-difference splitting, finite-volume scheme. The rigid-body dynamics equations are solved using a four-stage Runge-Kutta scheme. Once the wing reaches the limit-cycle response, an active control model using a mass injection system is applied from the wing surface to suppress the limit-cycle oscillation. The active control model is based on state feedback and the control law is established using pole placement techniques. The control law is based on the feedback of two states: the roll-angle and roll velocity. The primary model of the computational applications consists of a 80 deg swept, sharp edged, delta wing at 30 deg angle of attack in a freestream of Mach number 0.1 and Reynolds number of 0.4 x 10(exp 6). With a limit-cycle roll amplitude of 41.1 deg, the control model is applied, and the results show that within one and one half cycles of oscillation, the wing roll amplitude and velocity are brought to zero.

Kandil, Osama A.↗
NASA NTRS1986

Navier-Stokes computations for circulation controlled airfoils

Navier-Stokes computations of subsonic to transonic flow past airfoils with augmented lift due to rearward jet blowing over a curved trailing edge are presented. The approach uses a spiral grid topology. Solutions are obtained using a Navier-Stokes code which employs an implicit finite difference method, an algebraic turbulence model, and developments which improve stability, convergence, and accuracy. Results are compared against experiments for no jet blowing and moderate jet pressures and demonstrate the capability to compute these complicated flows.

Pulliam, T. H.↗
NASA NTRS1987

Navier-Stokes computations for circulation control airfoils

Navier-Stokes computations of subsonic to transonic flow past airfoils with augmented lift due to rearward jet blowing over a curved trailing edge are presented. The approach uses a spiral grid topology. Solutions are obtained using a Navier-Stokes code which employs an implicit finite difference method, an algebraic turbulence model, and developments which improve stability, convergence, and accuracy. Results are compared against experiments for no jet blowing and moderate jet pressures and demonstrate the capability to compute these complicated flows.

Pulliam, Thomas H.↗