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,531 records · Page 85

Development and application of computational aerothermodynamics flowfield computer codes

Multiple nozzle plume flow field is computed with a 3-D, Navier-Stokes solver. Numerical simulation is performed with a flux-split, two-factor, time asymptotic viscous flow solver of Ying and Steger. The two factor splitting provides a stable 3-D solution procedure under ideal-gas assumptions. An ad-hoc acceleration procedure that shows promise in improving the convergence rate by a factor of three for steady state problems is utilized. Computed solutions to generic problems at various altitude and flight conditions show flow field complexity and three-dimensional effects due to multiple nozzle jet interactions. Viscous, ideal gas solutions for the symmetric nozzle are compared with other numerical solutions.

Venkatapathy, Ethiraj↗

3-D plume flow computations with an upwind solver

In the present study multi-nozzle plume flows are computed with a three-dimensional, Navier-Stokes solver. Numerical simulations are performed with the flux-split, two-factor, time symptotic, viscous flow solver. The two factor splitting provides a stable three-dimensional solution procedure under ideal-gas assumptions. Viscous, ideal-gas solutions for a thin lip symmetrical nozzle are compared with experimental and numerical solutions. Computed solutions to axisymmetric and three-dimensional, multi-nozzle problems at various altitudes and flight conditions demonstrate flow field complexity and three-dimensional effects.

Venkatapathy, Ethiraj↗

Numerical solutions of 2-D multi-stage rotor/stator unsteady flow interactions

The Rai method of single-stage rotor/stator flow interaction is extended to handle multistage configurations. In this study, a two-dimensional Navier-Stokes multi-zone approach was used to investigate unsteady flow interactions within two multistage axial turbines. The governing equations are solved by an iterative, factored, implicit finite-difference, upwind algorithm. Numerical accuracy is checked by investigating the effect of time step size, the effect of subiteration in the Newton-Raphson technique, and the effect of full viscous versus thin-layer approximation. Computer results compared well with experimental data. Unsteady flow interactions, wake cutting, and the associated evolution of vortical entities are discussed.

Yang, R.-J.↗

Space Propulsion Design and Analysis

This software provides an improved methodology for predicting launcher base pressure and heat loads for RSRM (Reusable Solid Rocket Motor) launchers by accounting for complex anisotropic stress/strains and variable turbulent Prandtl and Schmidt numbers. A "building block" approach to turbulence model development, and validation has been applied for improved missile/launcher base region analysis. Modifications to existing kappa - epsilon turbulence models and application of scalar variance models are incorporated into a RANS-based method for aeropropulsive flow modeling, directly related to base flow methodology. (RANS stands for Reynolds-averaged Navier-Stokes.) The models are applied in a RANS solver framework and can improve analysis of other complex flow fields. The enhanced models provide a more accurate predictive capability for improving the design and analysis of RSRM launcher configuration. The kappa - epsilon model enhancements have been shown to improve the capability for predicting turbulence effects in base blow environments. The scalar variance models have been assessed over a wide range of flow configurations to improve prediction of turbulent scalar mixing.

Sinha, Neeraj↗

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approximate factorization for incompressible flow

For computational solution of the incompressible Navier-Stokes equations, the approximate factorization (AF) algorithm is used to solve the vectorized momentum equation in delta form based on the pressure calculated in the previous time step. The newly calculated velocities are substituted into the pressure equation (obtained from a linear combination of the continuity and momentum equation), which is then solved by means of line SOR. Computational results are presented for the NACA 66 sub 3 018 airfoil at Reynolds numbers of 1000 and 40,000 and attack angles of 0 and 6 degrees. Comparison with wind tunnel data for Re = 40,000 indicates good qualitative agreement between measured and calculated pressure distributions. Quantitative agreement is only fair, however, with the calculations somewhat displaced from the measurements. Furthermore, the computed velocity profiles are unrealistically thick around the airfoil, due to the excessive amount of artificial viscosity needed for stability. Based on the performance of the algorithm with regard to stability, it is concluded that AF/SOR is suitable for calculations at Reynolds numbers less than 10,000. Speedwise, the method is faster than point SOR by at least a factor of two.

Bernard, R. S.↗

Separated transonic airfoil flow calculations with a nonequilibrium turbulence model

Navier-Stokes transonic airfoil calculations based on a recently developed nonequilibrium, turbulence closure model are presented for a supercritical airfoil section at transonic cruise conditions and for a conventional airfoil section at shock-induced stall conditions. Comparisons with experimental data are presented which show that this nonequilibrium closure model performs significantly better than the popular Baldwin-Lomax and Cebeci-Smith equilibrium algebraic models when there is boundary-layer separation that results from the inviscid-viscous interactions.

King, L. S.↗

A space-marching method for the computation of viscous internal flows

A space-marching method has been developed to compute 3-D viscous flows in internal geometries. The Navier-Stokes equations have been posed as an initial-value problem by neglecting the effects of streamwise diffusion and treating the streamwise pressure gradient as a known source term. The fully coupled system of equations has been solved by a noniterative algorithm at each streamwise step of the computation. A low Mach number formulation of the equations has been used to compute incompressible flow fields. A computer program has been written to implement all aspects of the space-marching algorithm. The program is modular and is easily adapted to the widely varying geometries of internal flows. The space-marching algorithm has been tested by computing simple flows with known analytical solutions. The method has been used to predict complex 3-D turbulent flows. The algorithm is stable and very economical. A single sweep of the flow field by the space-marching method is approximately equivalent to one time-step of the time-marching method.

Govindan, T. R.↗

Simulation of vortex breakdown in tubes

The breakdown of a vortex flow in a tube is studied for a slightly diverging tube by means of a numerical solution of the Navier-Stokes equations for a quasi-stationary, incompressible, laminar flow. Numerical results are compared to experiments of Faler and Leibovich. The numerical data display features similar to the experimental results concerning breakdown conditions, location, and structure.

Menne, Stefan↗

Simulation of a three-dimensional vortex breakdown

The breakdown of a vortex flow in a tube is studied for a slightly diverging tube by means of a numerical solution of the Navier-Stokes equations for a quasi-stationary, incompressible, laminar flow. Numerical results are compared to experiments of Faler and Leibovich. The numerical data display features similar to the experimental results concerning breakdown conditions, location, and structure.

Liu, C. H.↗

Corrugated and Composite Nozzle-Inlets for Thrust and Noise Benefits

The following research results are based on development of an approach previously proposed and investigated in for optimum nozzle design to obtain maximum thrust. The design was denoted a Telescope nozzle. A Telescope nozzle contains one or several internal designs, which are inserted at certain locations into a divergent conical or planar main nozzle near its exit. Such a design provides additional thrust augmentation over 20% by comparison with the optimum single nozzle of equivalent lateral area, What is more, experimental acoustic tests have discovered an essential noise reduction due to application of Telescope nozzles. In this paper, some additional theoretical results are presented for Telescope nozzles and a similar approach is applied for aero-performance improvement of a supersonic inlet. Numerical simulations were conducted for supersonic flow into the divergent portion of a 2D or axisymmetric nozzle with several plane or conical designs as well as into a 2D or axisymmetric supersonic inlet with a forebody. The Kryko-Godunov marching numerical scheme for inviscid supersonic flows was used. Several cases were tested using the NASA CFL3d and IM/MSU Russian codes based on the full Navier-Stokes equations. Numerical simulations were conducted for non reacting flows (both codes) as well as for real high temperature gas flows with non-equilibrium chemical reactions (the latter code). In general, these simulations have confirmed essential benefits of Telescope design applications in propulsion system. Some preliminary numerical simulations of several typical inlet designs were conducted with the goal of inlet design optimization for maneuvering flight conditions.

Gilinsky, M.↗

Numerical solution of the time-dependent compressible Navier-Stokes equations in inlet regions

The results of a study to determine the effects of compressibility on the viscous flow through channels that have straight, parallel walls are presented. Two channel configurations are considered, the flow between two semi-infinite flat plates with uniform flow prescribed at the inlet plane and a cascade of semi-infinite flat plates with uniform flow introduced upstream. The flow field is modeled by using the time dependent, compressible Navier-Stokes equations. Time dependent solutions are obtained by using an explicit finite difference technique which advances the pressure on near field subsonic boundaries such that accurate steady state solutions are obtained. Steady state results at Reynolds number 20 and 150 are presented for Mach numbers between 0.09 and 0.36 and compared with the incompressible solutions of previous studies.

Olson, L. E.↗

The 3-D viscous flow CFD analysis of the propeller effect on an advanced ducted propeller subsonic inlet

A time marching Navier-Stokes code called PARC3D was used to study the 3-D viscous flow associated with an advanced ducted propeller (ADP) subsonic inlet at take-off operating conditions. At a free stream Mach number of 0.2, experimental data for the inlet-with-propeller test model indicated that the airflow was attached on the cowl windward lip at an angle of attack of 25 degrees became unstable at 29 degrees, and separated at 30 degrees. An experimental study with a similar inlet and with no propeller (through-flow) indicated that flow separation occurred at an angle of attack a few degrees below the value observed when the inlet was tested with the propeller. This tends to indicate that the propeller exerts a favorable effect on the inlet performance. During the through-flow experiment a stationary blockage device was used to successfully simulate the propeller effect on the inlet flow field at angles of attack. In the present numerical study, this flow blockage was modeled via a PARC3D computational boundary condition (BC) called the screen BC. The principle formulation of this BC was based on the one-and-half dimension actuator disk theory. This screen BC was applied at the inlet propeller face station of the computational grid. Numerical results were obtained with and without the screen BC. The application of the screen BC in this numerical study provided results which are similar to the results of past experimental efforts in which either the blockage device or the propeller was used.

Iek, Chanthy↗

Navier-Stokes Computations On Zonal Grids

Report describes numerical simulations of transonic flow of air about isolated wing and about wing and fuselage (not tail assembly) of F-16A aircraft. Fast, diagonalized Beam-Warming numerical-integration algorithm for solution of Euler and Navier-Stokes equations applied via Transonic Navier-Stokes computer program, using computational grids divided into zones of coarse and finer meshes.

Flores, Jolen↗

The curled wake model: a three-dimensional and extremely fast steady-state wake solver for wind plant flows

Abstract. Wind turbine wake models typically require approximations, such as wake superposition and deflection models, to accurately describe wake physics. However, capturing the phenomena of interest, such as the curled wake and interaction of multiple wakes, in wind power plant flows comes with an increased computational cost. To address this, we propose a new hybrid method that uses analytical solutions with an approximate form of the Reynolds-averaged Navier–Stokes equations to solve the time-averaged flow over a wind plant. We compare results from the solver to supervisory control and data acquisition data from the Lillgrund wind plant obtaining wake model predictions which are generally within 1 standard deviation of the mean power data. We perform simulations of flow over the Columbia River Gorge to demonstrate the capabilities of the model in complex terrain. We also apply the solver to a case with wake steering, which agreed well with large-eddy simulations. This new solver reduces the time – and therefore the related cost – it takes to simulate a steady-state wind plant flow (on the order of seconds using one core). Because the model is computationally efficient, it can also be used for different applications including wake steering for wind power plants and layout optimization.

17 WIND ENERGY↗

Large Eddy/Reynolds-Averaged Navier-Stokes Simulations of CUBRC Base Heating Experiments

ven with great advances in computational techniques and computing power during recent decades, the modeling of unsteady separated flows, such as those encountered in the wake of a re-entry vehicle, continues to be one of the most challenging problems in CFD. Of most interest to the aerothermodynamics community is accurately predicting transient heating loads on the base of a blunt body, which would result in reduced uncertainties and safety margins when designing a re-entry vehicle. However, the prediction of heat transfer can vary widely depending on the turbulence model employed. Therefore, selecting a turbulence model which realistically captures as much of the flow physics as possible will result in improved results. Reynolds Averaged Navier Stokes (RANS) models have become increasingly popular due to their good performance with attached flows, and the relatively quick turnaround time to obtain results. However, RANS methods cannot accurately simulate unsteady separated wake flows, and running direct numerical simulation (DNS) on such complex flows is currently too computationally expensive. Large Eddy Simulation (LES) techniques allow for the computation of the large eddies, which contain most of the Reynolds stress, while modeling the smaller (subgrid) eddies. This results in models which are more computationally expensive than RANS methods, but not as prohibitive as DNS. By complimenting an LES approach with a RANS model, a hybrid LES/RANS method resolves the larger turbulent scales away from surfaces with LES, and switches to a RANS model inside boundary layers. As pointed out by Bertin et al., this type of hybrid approach has shown a lot of promise for predicting turbulent flows, but work is needed to verify that these models work well in hypersonic flows. The very limited amounts of flight and experimental data available presents an additional challenge for researchers. Recently, a joint study by NASA and CUBRC has focused on collecting heat transfer data on the backshell of a scaled model of the Orion Multi-Purpose Crew Vehicle (MPCV). Heat augmentation effects due to the presence of cavities and RCS jet firings were also investigated. The high quality data produced by this effort presents a new set of data which can be used to assess the performance of CFD methods. In this work, a hybrid LES/RANS model developed at North Carolina State University (NCSU) is used to simulate several runs from these experiments, and evaluate the performance of high fidelity methods as compared to more typical RANS models. .

Salazar, Giovanni↗

Computation of the steady viscous flow over a tri-element 'augmentor wing' airfoil

The augmentor wing consists of a main airfoil with a slotted trailing edge for blowing, and two smaller aft airfoils which shroud the jet. This configuration has been modeled for numerical simulation by a novel discretization procedure which generates four separate grids: three surface-oriented airfoil grids and one outer free-stream grid. Grid lines and slopes are continuous across boundaries, so grid overlap at common boundaries provides boundary information without interpolation. A two-dimensional unsteady thin-layer Navier-Stokes code is used to calculate the flow for the no-blowing case at freestream Mach number = 0.7, Re = 12,600.000, and angles-of-incidence = 1.05 deg. Qualitative agreement with experimental data indicates the utility of this procedure in the analysis of multi-element configurations.

Lasinski, T. A.↗

Numerical simulation of viscous supersonic flow over a generic fighter configuration

A procedure is presented, as well as some results, to calculate the flow over a generic fighter configuration. A parabolized marching Navier-Stokes code is used to obtain the solution over a wing-canopy body. The flow conditions simulate supersonic cruise with a freestream Mach number of 2.169 and angles of attack of 4 and 10 deg. The body surface is considered to be an adiabatic wall and the flow is assumed to be turbulent for the given Reynolds number.

Chaussee, D. S.↗