Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stokes equations”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Solutions for laminar boundary layers with separation and reattachment

Numerical solutions of the laminar, incompressible boundary-layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified, and the pressure is deduced from the solution. A stream function/vorticity formulation of the boundary-layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite-difference equations, in which the flow direction is taken into account, are obtained by a successive column iteration scheme. Results are compared with Klineberg and Steger's (1974) separated boundary-layer calculations, and with Briley's (1971) solution of the Navier-Stokes equations for a separated region.

Carter, J. E.↗
NASA NTRS1976

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗
NASA NTRS1999

Use of Boundary Layer Transition Detection to Validate Full-Scale Flight Performance Predictions

Full-scale flight performance predictions can be made using CFD or a combination of CFD and analytical skin-friction predictions. However, no matter what method is used to obtain full-scale flight performance predictions knowledge of the boundary layer state is critical. The implementation of CFD codes solving the Navier-Stokes equations to obtain these predictions is still a time consuming, expensive process. In addition, to ultimately obtain accurate performance predictions the transition location must be fixed in the CFD model. An example, using the M2.4-7A geometry, of the change in Navier-Stokes solution with changes in transition and in turbulence model will be shown. Oil flow visualization using the M2.4-7A 4.0% scale model in the 14'x22' wind tunnel shows that fixing transition at 10% x/c in the CFD model best captures the flow physics of the wing flow field. A less costly method of obtaining full-scale performance predictions is the use of non-linear Euler codes or linear CFD codes, such as panel methods, combined with analytical skin-friction predictions. Again, knowledge of the boundary layer state is critical to the accurate determination of full-scale flight performance. Boundary layer transition detection has been performed at 0.3 and 0.9 Mach numbers over an extensive Reynolds number range using the 2.2% scale Reference H model in the NTF. A temperature sensitive paint system was used to determine the boundary layer state for these conditions. Data was obtained for three configurations: the baseline, undeflected flaps configuration; the transonic cruise configuration; and, the high-lift configuration. It was determined that at low Reynolds number conditions, in the 8 to 10 million Reynolds number range, the baseline configuration has extensive regions of laminar flow, in fact significantly more than analytical skin-friction methods predict. This configuration is fully turbulent at about 30 million Reynolds number for both 0.3 and 0.9, Mach numbers. Both the transonic cruise and the high-lift configurations were fully turbulent aft of the leading-edge flap hingeline at all Reynolds numbers.

Hamner, Marvine↗
NASA NTRS2005

Boundary-Layer Receptivity and Integrated Transition Prediction

The adjoint parabold stability equations (PSE) formulation is used to calculate the boundary layer receptivity to localized surface roughness and suction for compressible boundary layers. Receptivity efficiency functions predicted by the adjoint PSE approach agree well with results based on other nonparallel methods including linearized Navier-Stokes equations for both Tollmien-Schlichting waves and crossflow instability in swept wing boundary layers. The receptivity efficiency function can be regarded as the Green's function to the disturbance amplitude evolution in a nonparallel (growing) boundary layer. Given the Fourier transformed geometry factor distribution along the chordwise direction, the linear disturbance amplitude evolution for a finite size, distributed nonuniformity can be computed by evaluating the integral effects of both disturbance generation and linear amplification. The synergistic approach via the linear adjoint PSE for receptivity and nonlinear PSE for disturbance evolution downstream of the leading edge forms the basis for an integrated transition prediction tool. Eventually, such physics-based, high fidelity prediction methods could simulate the transition process from the disturbance generation through the nonlinear breakdown in a holistic manner.

Chang, Chau-Lyan↗
NASA NTRS1974

Hypersonic nonequilibrium flow past a sphere at low Reynolds numbers

Detailed flow field properties are computed in the nose region of a sphere from the Navier-Stokes equations for a range of Reynolds numbers. Comparisons are made of theoretical flow fields with those obtained in experiments and also with other analyses. In particular, the effects of various approximations on the flow field are investigated. It is shown that the Navier-Stokes model consistently predicts a smaller region of disturbance upstream of a sphere than that in experiments.

Li, C. P.↗
NASA NTRS2003

Computation of Large-Scale Structure Jet Noise Sources With Weak Nonlinear Effects Using Linear Euler

An approximate technique is presented for the prediction of the large-scale turbulent structure sound source in a supersonic jet. A linearized Euler equations code is used to solve for the flow disturbances within and near a jet with a given mean flow. Assuming a normal mode composition for the wave-like disturbances, the linear radial profiles are used in an integration of the Navier-Stokes equations. This results in a set of ordinary differential equations representing the weakly nonlinear self-interactions of the modes along with their interaction with the mean flow. Solutions are then used to correct the amplitude of the disturbances that represent the source of large-scale turbulent structure sound in the jet.

Dahl, Milo D.↗
NASA NTRS1975

Inverse solutions for laminar boundary-layer flows with separation and reattachment

Numerical solutions of the laminar, incompressible boundary layer equations are presented for flows involving separation and reattachment. Regular solutions are obtained with an inverse approach in which either the displacement thickness or the skin friction is specified; the pressure is deduced from the solution. A vorticity-stream-function formulation of the boundary layer equations is used to eliminate the unknown pressure. Solutions of the resulting finite difference equations, in which the flow direction is taken into account, are obtained by several global iteration schemes which are stable and have unconditional diagonal dominance. Results are compared with Klineberg and Steger's separated boundary layer calculations, and with Briley's solution of Navier-Stokes equations for a separated region. In addition, an approximate technique is presented in which the streamwise convection of vorticity is set equal to zero in the reversed flow region; such a technique results in a quick forward marching procedure for separated flows.

Carter, J. E.↗
NASA NTRS1977

Downstream boundary conditions for viscous flow problems

The problem of the specification of artificial outflow conditions in flow problems is studied. It is shown that for transport type equations incorrect outflow conditions will adversely affect the solution only in a small region near the outflow boundary, while for elliptic equations, e.g. those governing the streamfunction or pressure, a correct boundary specification is essential. In addition, integral outflow boundary conditions for fluid dynamical problems are considered. It is shown that such conditions are well posed, and their effect on the solutions of the Navier-Stokes equations is also considered.

Fix, G.↗
NASA NTRS1993

Analysis of the Harrier forebody/inlet design using computational techniques

Under the support of this Cooperative Agreement, computations of transonic flow past the complex forebody/inlet configuration of the AV-8B Harrier II have been performed. The actual aircraft configuration was measured and its surface and surrounding domain were defined using computational structured grids. The thin-layer Navier-Stokes equations were used to model the flow along with the Chimera embedded multi-grid technique. A fully conservative, alternating direction implicit (ADI), approximately-factored, partially flux-split algorithm was employed to perform the computation. An existing code was altered to conform with the needs of the study, and some special engine face boundary conditions were developed. The algorithm incorporated the Chimera technique and an algebraic turbulence model in order to deal with the embedded multi-grids and viscous governing equations. Comparison with experimental data has yielded good agreement for the simplifications incorporated into the analysis. The aim of the present research was to provide a methodology for the numerical solution of complex, combined external/internal flows. This is the first time-dependent Navier-Stokes solution for a geometry in which the fuselage and inlet share a wall. The results indicate the methodology used here is a viable tool for transonic aircraft modeling.

Chow, Chuen-Yen↗
NASA NTRS2001

Development of Flow and Heat Transfer Models for the Carbon Fiber Rope in Nozzle Joints of the Space Shuttle Reusable Solid Rocket Motor

Methodologies have been developed for modeling both gas dynamics and heat transfer inside the carbon fiber rope (CFR) for applications in the space shuttle reusable solid rocket motor joints. Specifically, the CFR is modeled using an equivalent rectangular duct with a cross-section area, friction factor and heat transfer coefficient such that this duct has the same amount of mass flow rate, pressure drop, and heat transfer rate as the CFR. An equation for the friction factor is derived based on the Darcy-Forschheimer law and the heat transfer coefficient is obtained from pipe flow correlations. The pressure, temperature and velocity of the gas inside the CFR are calculated using the one-dimensional Navier-Stokes equations. Various subscale tests, both cold flow and hot flow, have been carried out to validate and refine this CFR model. In particular, the following three types of testing were used: (1) cold flow in a RSRM nozzle-to-case joint geometry, (2) cold flow in a RSRM nozzle joint No. 2 geometry, and (3) hot flow in a RSRM nozzle joint environment simulator. The predicted pressure and temperature history are compared with experimental measurements. The effects of various input parameters for the model are discussed in detail.

Wang, Q.↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

An approximate calculation of the strong interaction on a transonic airfoil

This paper describes an interactive computational scheme that calculates the transonic flow field over an airfoil. The inviscid flow, the turbulent boundary layer in the weak interaction region, and the elliptic flow in the strong interaction region are computed separately, and the solutions are matched iteratively in a direct, second-order interaction scheme. The pressure distribution, obtained from a transonic, small-perturbation solution, is prescribed on the boundary layer and the Navier-Stokes equations and the ensuing displacement thicknesses are used to modify the inviscid flow. The best overall results are obtained by using a turbulence model that eliminates the eddy viscosity in the strong interaction region. Computed pressure distributions are compared with experimental data and with Navier-Stokes solutions. In weak interactions, the agreement with experimental data is very good, as it is also in strong interactions except for the pressures in the separated-flow region. There the agreement is qualitative only, and the pressure recovery is too high. The results in both cases are superior or equal to those of the Navier-Stokes solutions, and the present method is significantly faster.

Seginer, A.↗
NASA NTRS1996

Two-Dimensional Dynamic Simulation of a Continuous Foil Bearing

In this paper, the two dimensional(radial and circumferential) transient Navier-Stokes equations are used to solve the hydrodynamic problem in conjunction with the time dependent motion of the journal, and the deformable, spring supported foil. The elastic deformation of the foil and its supports are simulated by a finite element model. The time-dependent Navier-Stokes formulation is used to solve for the interaction between the fluid lubricant, the motion of the journal and the deformable foil boundary. The steady state, the quasi-transient and the full transient dynamic simulation of the foil-fluid journal interaction are examined on a comparative basis. For the steady state simulation, the fluid lubricant pressures are evaluated for a particular journal position, by means of an iterative scheme until convergence is achieved in both the fluid pressures and the corresponding foil deformation. For the quasi-transient case, the transient motion of the journal is calculated using a numerical integration scheme for the velocity and displacement of the journal. The deformation of the foil is evaluated through numerical iteration in feedback mode with the fluid film pressure generated by the journal motion until convergence at every time step is achieved. For the full transient simulation, a parallel real-time integration scheme is used to evaluate simultaneously the new journal position and the new deformed shape of the foil at each time step. The pressure of the fluid lubricant is iterated jointly with the corresponding journal position and the deformed foil geometry until convergence is achieved. A variable time-stepping Newmark-Beta integration procedure is used to evaluate the transient dynamics at each time step of the bearing.

Braun, M. Jack↗
NASA NTRS2016

Investigation of Damping Physics and CFD Tool Validation for Simulation of Baffled Tanks at Variable Slosh Amplitude

To meet the flight control damping requirement, baffles of various configurations have been devised to increase the natural viscous damping and decrease the magnitude of the slosh forces and torques. In the design of slosh baffles, the most widely used damping equation is the one derived by Miles, which is based on the experiments of Keulegan and Carpenter. This equation has been used in predicting damping of the baffled tanks in different diameters ranging from 12 to 112 inches. The analytical expression of Miles equation is easy to use, especially in the design of complex baffle system. Previous investigations revealed that some experiments had shown good agreements with the prediction method of Miles, whereas other experiments have shown significant deviations. For example, damping from Miles equation differs from experimental measurements by as much as 100 percent over a range of tank diameters from 12 to 112 inches, oscillation amplitudes from 0.1 to 1.5 baffle widths, and baffle depths of 0.3 to 0.5 tank radius. Previously, much of this difference has been attributed to experimental scatter. A systematical study is needed to understand the damping physics of baffled tanks, to identify the difference between Miles equation and experimental measurement, and to develop new semi-empirical relations to better represent the real damping physics. The approach of this study is to use CFD technology to shed light on the damping mechanisms of a baffled tank. First, a 1-D Navier-Stokes equation representing different length scales and time scales in the baffle damping physics is developed and analyzed. A well validated CFD solver, developed at NASA MSFC, Loci-STREAM-VOF, is applied to study vorticity field around the baffle and around the fluid interface to highlight the dissipation mechanisms at different slosh amplitudes. Previous measurement data are then used to validate the CFD damping results. The study found several critical parameters controlling fluid damping from a baffle: local slosh amplitude to baffle thickness (A/t), surface liquid depth to tank radius (h/R), local slosh amplitude to baffle width (A/W); and non-dimensional slosh frequency. The simulation highlights three significant damping regimes where different mechanisms dominate. The study proves that the previously found discrepancies between Miles equation and experimental measurement are not due to the measurement scatter, but rather due to different damping mechanisms at various slosh amplitudes. The limitations on the use of Miles equation are discussed based on the flow regime.

Yang, H. Q.↗
NASA NTRS1986

A continuum analysis of chemical nonequilibrium under hypersonic low-density flight conditions

Results of employing the continuum model of Navier-Stokes equations under the low-density flight conditions are presented. These results are obtained with chemical nonequilibrium and multicomponent surface slip boundary conditions. The conditions analyzed are those encountered by the nose region of the Space Shuttle Orbiter during reentry. A detailed comparison of the Navier-Stokes (NS) results is made with the viscous shock-layer (VSL) and direct simulation Monte Carlo (DSMC) predictions. With the inclusion of new surface-slip boundary conditions in NS calculations, the surface heat transfer and other flowfield quantities adjacent to the surface are predicted favorably with the DSMC calculations from 75 km to 115 km in altitude. This suggests a much wider practical range for the applicability of Navier-Stokes solutions than previously thought. This is appealing because the continuum (NS and VSL) methods are commonly used to solve the fluid flow problems and are less demanding in terms of computer resource requirements than the noncontinuum (DSMC) methods.

Gupta, R. N.↗
NASA NTRS1995

High-order essentially non-oscillatory methods for computational aeroacoustics

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. In addition, it is desirable to highly resolve discontinuities that can develop in the solutions of the Euler or Navier-Stokes equations. The class of essentially non-oscillatory (ENO) shock-capturing schemes has been designed to have both of these properties. The dual capacity of ENO schemes for high-order accuracy and non-oscillatory shock-capturing is achieved through the use of adaptive stenciling, which makes these schemes highly nonlinear. These schemes are briefly described and referenced herein. A fourth-order algorithm is then applied to the solution of an acoustic wave in a quasi-one-dimensional converging-diverging nozzle.

Casper, Jay↗
NASA NTRS2001

Hypersonic Boundary Layer Instability Over a Corner

A boundary-layer transition study over a compression corner was conducted under a hypersonic flow condition. Due to the discontinuities in boundary layer flow, the full Navier-Stokes equations were solved to simulate the development of disturbance in the boundary layer. A linear stability analysis and PSE method were used to get the initial disturbance for parallel and non-parallel flow respectively. A 2-D code was developed to solve the full Navier-stokes by using WENO(weighted essentially non-oscillating) scheme. The given numerical results show the evolution of the linear disturbance for the most amplified disturbance in supersonic and hypersonic flow over a compression ramp. The nonlinear computations also determined the minimal amplitudes necessary to cause transition at a designed location.

Balakumar, Ponnampalam↗
NASA NTRS1984

Prediction of transonic separated flows

Johnson et al. (1982) have provided a detailed comparison between a thoroughly documented transonic flow with shock-induced separations and solutions of the flow using the Navier-Stokes equations. According to this comparison, there were several deficiencies in the computations. The present investigation takes into account new experimental data which have been obtained in a larger wind tunnel with the same test model for a wider range of freestream Mach numbers. The results of new Navier-Stokes computations using more compatible boundary conditions are shown, and the effects of the turbulence model choice on predicting Mach number trends are assessed.

Horstman, C. C.↗