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 703 records · Page 39

Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers

The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long.

Kinney, R. B.↗

Hypersonic turbulent wall boundary layer computations

The Baldwin-Lomax algebraic turbulence model was modified for hypersonic flow conditions. Two coefficients in the outer layer eddy viscosity model were determined as functions of Mach number and temperature ratio. By matching the solutions from the Baldwin-Lomax model to those from the Cebeci-Smith model for a flat plate at hypersonic speed, the new values of the coefficient were obtained. The results show that the values of C sub cp and C sub kleb are functions of both Mach number and wall temperature ratio. The C sub cp and C sub kleb variations with Mach number and wall temperature were used for the calculations of both a 4 deg wedge flow at Mach 18 and an axisymmetric Mach 20 nozzle flow. The Navier-Stokes equations with thin layer approximation were solved for the above hypersonic flow conditions and the results were compared with existing experimental data. The agreement between the numerical solutions and the existing experimental data were good. The modified Baldwin-Lomax model thus is useful in the computations of hypersonic flows.

Kim, S. C.↗

Direct numerical study of leading-edge contamination

Instability, turbulence, and relaminarization in the attachment-line region of swept and unswept cylindrical bodies are studied by numerical solution of the full Navier-Stokes equations. The flow is simulated over a strip containing the attachment-line and treated as homogeneous in the spanwise direction; the disturbances decay exponentially upstream. Transpiration through the wall may be prescribed. The new method, which admits completely general disturbance, agrees with published linear-stability results, which were limited to an apparently restrictive form of disturbance. Fully developed turbulent solutions with sweep are generated and compare well with the experiment. The turbulence is subcritical (except for blowing), resulting in large hysteresis loops. By lowering the sweep Reynolds number, or increasing the suction, the turbulent flow is made to relaminarize. The relaminarization Reynolds number is much less sensitive to suction than the linear-stability Reynolds number. Extensive attempts to detect the postulated nonlinear instability of the unswept flow failed, suggesting that this flow is linearly and nonlinearly stable.

Spalart, P. R.↗

CFD code evaluation for internal flow modeling

Research on the computational fluid dynamics (CFD) code evaluation with emphasis on supercomputing in reacting flows is discussed. Advantages of unstructured grids, multigrids, adaptive methods, improved flow solvers, vector processing, parallel processing, and reduction of memory requirements are discussed. As examples, researchers include applications of supercomputing to reacting flow Navier-Stokes equations including shock waves and turbulence and combustion instability problems associated with solid and liquid propellants. Evaluation of codes developed by other organizations are not included. Instead, the basic criteria for accuracy and efficiency have been established, and some applications on rocket combustion have been made. Research toward an ultimate goal, the most accurate and efficient CFD code, is in progress and will continue for years to come.

Chung, T. J.↗

High alpha aerodynamic control by tangential fuselage blowing

This work explores the effect of tangential blowing on the vortical structures that develop around a tangent-ogive cylinder configuration at high angle of attack. A lateral force results if blowing is applied asymmetrically. The study is conducted numerically by solving the three-dimensional, compressible-flow Navier-Stokes equations. The computation was done for a Reynolds number of 52,000, Mach number of 0.2, blowing momentum coefficients of 0.0, 0.1, 0.2, and 0.4, and angle of attack of 10, 30, and 45 deg. Only asymmetrical blowing was considered. Surface streamlines, helicity contours and pressure distributions are obtained for each case.

Font, G. I.↗

Computer Code Generates Homotopic Grids

HOMAR is computer code using homotopic procedure to produce two-dimensional grids in cross-sectional planes, which grids then stacked to produce quasi-three-dimensional grid systems for aerospace configurations. Program produces grids for use in both Euler and Navier-Stokes computation of flows. Written in FORTRAN 77.

Moitra, Anutosh↗

Comparison Of Two Viscous-Flow Computer Codes

Two viscous-flow computer codes compared by applying them to five test cases of steady-state transonic viscous flows about transonic airfoils. Two codes were: FLOMG, which solves Navier-Stokes differential equations of flow by implementing explicit, Runge-Kutta, finite-volume, multigrid numerical-integration algorithm; and ARC2D, which implements implicit, finite-difference, approximate-factorization, eigenvector-diagonalization numerical-integration algorithm.

Maksymiuk, C. M.↗

Parallelization of a Parabolized Navier-Stokes Solver with a Design Optimizer

The design of future supersonic aircraft, such as the High Speed Civil Transport (HSCT), will rely heavily on computational methods for aircraft design and the prediction of the complex aerodynamics encountered in flight. Parabolized Navier-Stokes (PNS) equation flow solvers are recognized as efficient and accurate computational tools for the solution of supersonic and hypersonic flow-fields, while design optimizers have the potential to be valuable tools within the overall design process. Presently, however, the execution of the flow solver in conjunction with a design optimizer presents a computationally intensive and formidable problem. To meet the challenges and increasing demand for multidisciplinary numerical tools which are faster, more robust and provide greater functionality, alternative strategies are explored to increase computational throughput by coupling a design optimizer and flow solver in a parallel processing environment. To address this problem the parallel processing of a PNS flow solver with a nonlinear constraint design optimizer is investigated as an alternative computational method.

Pallis, J. M.↗

Direct Numerical Simulation of a Temporally Evolving Incompressible Plane Wake: Effect of Initial Conditions on Evolution and Topology

Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Sondergaard, R.↗

Overview of High-Fidelity Modeling Activities in the Numerical Propulsion System Simulations (NPSS) Project

A high-fidelity simulation of a commercial turbofan engine has been created as part of the Numerical Propulsion System Simulation Project. The high-fidelity computer simulation utilizes computer models that were developed at NASA Glenn Research Center in cooperation with turbofan engine manufacturers. The average-passage (APNASA) Navier-Stokes based viscous flow computer code is used to simulate the 3D flow in the compressors and turbines of the advanced commercial turbofan engine. The 3D National Combustion Code (NCC) is used to simulate the flow and chemistry in the advanced aircraft combustor. The APNASA turbomachinery code and the NCC combustor code exchange boundary conditions at the interface planes at the combustor inlet and exit. This computer simulation technique can evaluate engine performance at steady operating conditions. The 3D flow models provide detailed knowledge of the airflow within the fan and compressor, the high and low pressure turbines, and the flow and chemistry within the combustor. The models simulate the performance of the engine at operating conditions that include sea level takeoff and the altitude cruise condition.

Veres, Joseph P.↗

Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.

Swanson, Roy C., Jr.↗

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

A computer code has been developed to solve the full two-dimensional Navier-Stokes equations in a scramjet inlet. The analysis uses a numerical coordinate transformation which generates a set of boundary fitted curvilinear coordinates. The explicit finite-difference algorithm of MacCormack is used to solve the governing equations. A two-layer eddy viscosity model is used for the turbulent flow. The code can analyze both inviscid and viscous flows with multiple struts in the flow field. Detailed results are presented for two model problems and two scramjet inlets with one and two struts. The application of the two-dimensional analysis in the preliminary design of the actual scramjet inlet is briefly discussed.

Kumar, A.↗

Analysis of the Scramjet inlet flow field using two-dimensional Navier-Stokes equations

A computer code was developed to solve the full two dimensional Navier-Stokes equations in a scramjet inlet. The analysis uses a numerical coordinate transformation which generates a set of boundary-fitted curvilinear coordinates. The explicit finite difference algorithm of MacCormack is used to solve the governing equations. A two-layer eddy viscosity model is used for the turbulent flow. The code can analyze both inviscid and viscous flows with multiple struts in the flow field. Detailed results are presented for two model problems and two scramjet inlets with one and two struts. The application of the two dimensional analysis in the preliminary design of the actual scramjet inlet is briefly discussed.

Kumar, A.↗

A sharp immersed method for 2D flow-body interactions using the vorticity-velocity Navier-Stokes equations

Immersed methods discretize boundary conditions for complex geometries on background Cartesian grids. Here, this makes such methods especially suitable for two-way coupled flow-body problems, where the body mechanics are partially driven by hydrodynamic forces. However, for the vorticity-velocity form of the Navier-Stokes equations, existing immersed geometry discretizations for two-way coupled problems only achieve first order spatial accuracy near solid boundaries. Here we introduce a sharp-interface approach based on the immersed interface method to handle the one- and two-way coupling between an incompressible flow and one or more rigid bodies using the 2D vorticity-velocity Navier-Stokes equations. Our main contributions are three-fold. First, we develop and analyze a moving boundary treatment for sharp immersed methods that can be applied to PDEs with implicitly defined boundary conditions, such as those commonly imposed on the vorticity field. Second, we develop a two-way coupling methodology for the vorticity-velocity Navier-Stokes equations based on control-volume momentum balance that does not require the pressure field. Third, we show through extensive testing and validation that our resulting flow-body solver reaches second-order accuracy for most practical scenarios, and provides significant efficiency benefits compared to a representative first-order approach.

42 ENGINEERING↗

Computation of non-equilibrium flow downstream of a plasma torch

Numerical solutions of the Navier-Stokes equations for compressible reacting flow in an axisymmetric geometry are presented for a nitrogen plasma torch with both thermal and chemical nonequilibrium. The Navier-Stokes equations are solved using a new axisymmetric finite element/finite volume formulation in which the convective flux is treated by an extension of the approximate Riemann solver due to Osher. The numerical scheme is validated by comparison with a previous solution of the same problem using a different computational scheme. Results obtained using two different models of nonthermal dissociation rates are compared to experimental data.

Harle, Christophe↗

Numerical simulation of F-18 fuselage forebody flows at high angles of attack

Fine-grid Navier-Stokes solutions were obtained for flow over the fuselage forebody and wing leading edge extension of the F/A-18 High Alpha Research Vehicle at large incidence. The resulting flows are complex, and exhibit cross flow separation from the sides of the forebody and from the leading edge extension. A well-defined vortex pattern is observed in the leeward-side flow. Results obtained for laminar flow show good agreement with flow visualizations obtained in ground-based experiments. Further, turbulent flows computed at high Reynolds-number flight-test conditions show good agreement with surface and off-surface visualizations obtained in flight.

Schiff, Lewis B.↗

Computations of supersonic viscous flow over a finite-width plate

Finite difference methods were applied to solve the parabolic Navier-Stokes equations for the flow over a finite width plate at 0 deg and 10 deg angles of attack. The methods were developed on the basis of the operator factorization concept resulting in the split of a three dimensional equation into successive two dimensional equations. Backward and centered implicit factorization schemes, were used and their results were compared. Available numerical solutions and experimental data obtained at low Reynolds number conditions were also used for comparison. The backward implicit method provides a more successful solution, which ranges from the merged layer to the strong interaction regimes. Detailed structures were revealed of the shear layer around and behind the side edge.

Li, C. P.↗

CFD applications in pump flows

The objective of the proposed paper is to develop a computational procedure that solves incompressible Navier-Stokes equations for pump flows. The solution method is based on the pseudo-compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. As a benchmark problem, the flow through the Rocketdyne inducer is numerically simulated. A coarse grid solution is obtained with a single zone by using an algebraic turbulence model. In multi-zone fine grid computation, a one-equation Baldwin-Barth turbulence model is utilized. Numerical results are compared with experimental measurements and a good agreement is found between the two. The resulting computer code is then applied to the flow analysis inside a two-stage fuel pump impeller operating at 80 percent, 100 percent, and 120 percent of design flow.

Kiris, Cetin↗