Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscous effects”

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,405 records · Page 78

A High Angle of Attack Inviscid Shuttle Orbiter Computation

An accurate description of the aerothermal environment is required to minimize the weight of the Thermal Protection System required on the leeside of winged reentry vehicles. The inability of ground-test facilities to reproduce the high enthalpy, separated flow present during reentry flight conditions, coupled with the prohibitive expense of flight tests, leads to the use of an analytical method - namely Computational Fluid Dynamics (CFD) - to describe the flow. While the ultimate goal of this work is to accurately predict the leeside flow and its associated thermal environment, an essential and reasonable first step towards that goal is a comparison of pressure predictions by a code with wind-tunnel data. Until such CFD pressure predictions agree with wind-tunnel test cases, there is little hope of accurately predicting the thermal environment at flight conditions. Thus, the objective of this study is to compare the pressure distributions predicted by inviscid, perfect gas CFD to Shuttle Orbiter wind-tunnel data and to address any significant issues encountered during the computation. While flight data is available for the Shuttle Orbiter, a wind-tunnel case is chosen for this study to allow a tractable problem for preliminary investigation. A wind-tunnel case allows the perfect gas assumption for the flow chemistry. This provides a significant computational savings over a several species finite-rate chemistry model which would be necessary if high-temperature effects present at flight conditions were to be included. In addition, by concentrating on the surface pressures, the analysis need only consider inviscid flow for general evaluation of the code capability. This further reduces the computational expense due to the absence of viscous terms and the associated decrease in the number of points required for the computational grid. Previous computational efforts (such as STEIN and HALIS) have been directed toward the windward surface quantities, primarily due to restrictions in treating either the winged geometry or its associated subsonic regions at high angle of attack. The code used for this study, the LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) code of Gnoffo, represents a state-of-the-art code for computing the flow over complex configurations at hypersonic speeds. In the study, the LAURA code is applied to a wind-tunnel condition to initiate the assessment of the code's ability to predict the flow over a relatively complex hypersonic vehicle at high angles of attack. This presentation is used to highlight the pertinent results of a more detailed investigation of the inviscid calculation over the Shuttle Orbiter with the LAURA code.

William L. Kleb↗

A Patched-Grid Algorithm for Complex Configurations Directed Towards the F/A-18 Aircraft

A patched-grid algorithm for the analysis of complex configurations with an implicit, upwind-biased Navier-Stokes solver is presented. Results from both a spatial-flux and a time-flux conservation approach to patching across zonal boundaries are presented. A generalized coordinate transformation with a biquadratic geometric element is used at the zonal interface in order to treat highly stretched viscous grids and arbitrarily-shaped zonal boundaries. Applications are made to the F-18 forebody-strake configuration at subsonic, high-alpha conditions. Computed surface flow patterns compare well with ground-based and flight-test results; the large effect of Reynolds number on the forebody flow-field is shown.

High alpha research vehicle↗

Adjoint-Based Anisotropic Mesh Adaptation for a Stabilized Finite-Element Flow Solver

An adjoint solver is implemented in the FUN3D stabilized finite-element flow solver. The adjoint solution is used to generate anisotropic, adapted meshes to control error in scalar out- put functionals, such as lift or drag coefficient. The inviscid and turbulent adjoints are verified with a finite-difference derivative approximation and can be used in design optimization in addition to mesh adaptation. The adjoint capability represents an extension of previous research using the FUN3D stabilized finite-element capability with metric-based mesh adaptation and interpolation-based error estimates to generate highly anisotropic adapted meshes for turbulent flows. In the present work, a metric-based approach is again utilized, where the adjoint and the primal solutions both contribute to the generation of a metric tensor field that is subsequently used to produce the required anisotropic mesh for each adaptation cycle. Adaptive results are then shown for an inviscid supersonic flow over a diamond airfoil, inviscid transonic flow over the ONERA M6 wing, and viscous laminar flow over the NACA 0012 airfoil, all using drag force as the output functional. Adjoint-based adaptation is compared with a multiscale solution-based approach that controls the L norm of Mach number interpolation error to demonstrate the effectiveness and effciency of the adjoint-based adaptive mesh technology.

Balan, Aravind↗

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Upgraded viscous flow analysis of multi-element airfoils

A description of an improved version of the NASA/Lockheed multi-element airfoil analysis computer program is presented. The improvements include several major modifications of the aerodynamic model as well as substantial changes of the computer code. The modifications of the aerodynamic model comprise the representation of the boundary layer and wake displacement effects with an equivalent source distribution, the prediction of wake parameters with Green's lag-entrainment method, the calculation of turbulent boundary layer separation with the method of Nash and Hicks, the estimation of the onset of confluent boundary layer separation with a modified form of Goradia's method, and the prediction of profile drag with the formula of Squire and Young. The modifications of the computer program for which the structured approach to computer software development was employed are also described. Important aspects of the structured program development such as the functional decomposition of the aerodynamic theory and its numerical implementation, the analysis of the data flow within the code, and the application of a pseudo code are discussed.

Brune, G. W.↗

Significance of shock and body slip conditions on Jovian entry heating

The influence of the body and shock slip conditions on the heating of a Jovian entry body is investigated. The flow in the shock layer is considered to be axisymmetric, steady, laminar, viscous, and in chemical equilibrium. Realistic thermophysical and step-function spectral models are employed and results are obtained by implicit finite-difference and iteractive procedures. The freestream conditions correspond to a typical Jovian entry trajectory point. The results indicate that the effect of the slip conditions is significant when the altitudes are higher than 225 km and that the contribution of a radiative heat-flux term in the energy equation should not be neglected at any altitude.

Tiwari, S. N.↗

Flow analysis of nozzle installations with strong airplane flow interactions

A numerical procedure has been developed to calculate the flow fields resulting from the viscous-inviscid interactions that occur when a strong jet exhaust and aircraft flow field coupling exists. The approach used in the current procedure is to divide the interaction region into zones which are either predominantly viscous or inviscid. The flow in the inviscid zone, which surrounds most of the aircraft, is calculated using an existing linearized potential flow code. The viscous flow zone, which encompasses the jet plume, is modeled using a parabolized Navier-Stokes code. The key feature of the present procedure is the coupling of the zonal solutions such that sufficient information is transferred between the zones to preserve the effects of the interactions. The zonal boundaries overlap with the boundary conditions being the information link between zones. An iteraction scheme iterates the coupled analysis until convergence has been obtained. The procedure has been successfully used for several test cases for which the computed results are presented.

Roberts, D. W.↗

Analytical modeling of circuit aerodynamics in the new NASA Lewis wind tunnel

Rehabilitation and extention of the capability of the altitude wind tunnel (AWT) was analyzed. The analytical modeling program involves the use of advanced axisymmetric and three dimensional viscous analyses to compute the flow through the various AWT components. Results for the analytical modeling of the high speed leg aerodynamics are presented; these include: an evaluation of the flow quality at the entrance to the test section, an investigation of the effects of test section bleed for different model blockages, and an examination of three dimensional effects in the diffuser due to reentry flow and due to the change in cross sectional shape of the exhaust scoop.

Towne, C. E.↗

Analytical modeling of circuit aerodynamics in the new NASA Lewis Altitude Wind Tunnel

Rehabilitation and extention of the capability of the altitude wind tunnel (AWT) was analyzed. The analytical modelling program involves the use of advanced axisymmetric and three dimensional viscous analyses to compute the flow through the various AWT components. Results for the analytical modelling of the high speed leg aerodynamics are presented; these include: an evaluation of the flow quality at the entrance to the test section, an investigation of the effects of test section bleed for different model blockages, and an examination of three dimensional effects in the diffuser due to reentry flow and due to the change in cross sectional shape of the exhaust scoop.

Towne, C. E.↗

An asymptotic investigation of the stationary modes of instability of the boundary layer on a rotating disc

An investigation of high Reynolds number stationary instabilities in the boundary layer on a rotating disc is given. It is shown that in addition to the inviscid mode at high Reynolds numbers, there is a stationary short wavelength mode. This mode has its structure fixed by a balance between viscous and Coriolis forces and cannot be described by an inviscid theory. The asymptotic structure of the wavenumber and orientation of this mode is obtained. A similar analysis is given for the inviscid mode, the expansion procedure used is capable of taking nonparallel effects into account in a self consistent manner. The results are compared to numerical calculations and experimental observations.

Hall, P.↗

An asymptotic investigation of the stationary modes of instability of the boundary layer on a rotating disc

An investigation of high Reynolds number stationary instabilities in the boundary layer on a rotating disc is given. It is shown that in addition to the inviscid mode at high Reynolds numbers, there is a stationary short wavelength mode. This mode has its structure fixed by a balance between viscous and Coriolis forces and cannot be described by an inviscid theory. The asymptotic structure of the wavenumber and orientation of this mode is obtained. A similar analysis is given for the inviscid mode, the expansion procedure used is capable of taking nonparallel effects into account in a self consistent manner. The results are compared to numerical calculations and experimental observations.

Hall, P.↗

An explicit Runge-Kutta method for turbulent reacting flow calculations

The paper presents a numerical method for the solution of the conservation equations governing steady, reacting, turbulent viscous flow in two-dimensional geometries, in both Cartesian and axisymmetric coordinates. These equations are written in Favre-averaged form and closed with a first order model. A two-equation K-epsilon model, where low Reynolds number and compressibility effects are included, and a modified eddy-break up model are used to simulate fluid mechanics turbulence, chemistry and turbulence-combustion interaction. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multi-stage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. The method is applied to the computation of both diffusion and premixed turbulent reacting flows. The computed temperature distributions compare favorably with experimental data.

Boretti, A. A.↗

Several families of viscous optimized waveriders - A review of waverider research at the University of Maryland

Theoretical and experimental research conducted by the Hypersonic Group at the University of Maryland in the area of viscous optimized waveriders is reviewed. In particular, the motivation for this class of waveriders is discussed, and results obtained to date are presented. Various waverider families generated from different flows and designed for different applications are examined, including waveriders using chemically reacting flow effects and very high altitude applications.

Anderson, John D., Jr.↗

Departure solutions of the unsteady thin-layer and full Navier-Stokes equations solved using streamline curvature based iteration techniques

The development of a thorough understanding of the mechanisms for vortex eruptions from viscous layers, which are believed to be associated with phenomena such as dynamic stall onset and transition, is crucial if accurate models of such phenomena are to be formulated. The development of such models may, in turn, allow for the possibility that such effects could be accounted for during the design of various aerodynamic devices such as wings, helicopter rotors, and turbomachinery blading and thus lead to designs which are stall free or stall resistant and which have better stall-recovery properties. The present investigation is being conducted as part of an effort to develop analytical and numerical tools which can be used to help improve our understanding of the vortex-eruption mechanism at high Reynolds numbers. The addition of the normal-momentum equation to the classical unsteady boundary-layer equations is crucial according to recent asymptotic analyses of the vortex-eruption problem and is a key feature of the analyses being developed by the present authors. The purpose of this paper is as follows: to describe departure solution behavior observed when using unsteady, streamline-curvature based solution procedures in which nontrivial transverse pressure gradient effects are included; and to show that special treatment of the time-derivative of the normal velocity is needed to eliminate the ill-posed solution behavior, which is observed when small spatial and temporal step sizes are used.

Barnett, M.↗

Parsec-Scale Accretion and Winds Irradiated by a Quasar

We present numerical simulations of properties of a parsec-scale torus exposed to illumination by the central black hole in an active galactic nucleus (AGN). Our physical model allows to investigate the balance between the formation of winds and accretion simultaneously. Radiation-driven winds are allowed by taking into account radiation pressure due to UV and IR radiation along with X-ray heating and dust sublimation. Accretion is allowed through angular momentum transport and the solution of the equations of radiative, viscous radiation hydrodynamics. Our methods adopt flux-limited diffusion radiation hydrodynamics for the dusty, infrared pressure driven part of the flow, along with X-ray heating and cooling. Angular momentum transport in the accreting part of the flow is modeled using effective viscosity. Our results demonstrate that radiation pressure on dust can play an important role in shaping AGN obscuration. For example, when the luminosity illuminating the torus exceeds L greater than 0.01 L(sub Edd), where L(sub Edd) is the Eddington luminosity, we find no episodes of sustained disk accretion because radiation pressure does not allow a disk to form. Despite the absence of the disk accretion, the flow of gas to smaller radii still proceeds at a rate 10(exp -4)-10(exp -1)M dot yr(exp -1) through the capturing of the gas from the hot evaporative flow, thus providing a mechanism to deliver gas from a radiation-pressure dominated torus to the inner accretion disk. As L L(sub edd) increases, larger radiation input leads to larger torus aspect ratios and increased obscuration of the central black hole. We also find the important role of the X-ray heated gas in shaping the obscuring torus.

accretion disks - galaxies: active - galaxies: Sey↗

Parameter identification methods for improving structural dynamic models

There is an increasing need to develop Parameter Identification methods for improving structural dynamic models, based on the inability of engineers to produce mathematical models which correlate with experimental data. This research explores the efficiency of combining Component Mode Synthesis (substructuring) methods with Parameter Identification procedures in order to improve analytical modeling of structural components and their connections. Improvements are computed in terms of physical stiffness and damping parameters in order that the physical characteristics of the model can be better understood. Connections involving both viscous and friction damping are investigated. Substructuring methods are utilized to reduce the complexity of the identification problem. Component and inter-component structural connection properties are evaluated and identified independently, thus simplifying the identification problem. It is shown that modal test data is effective for identifying modeling problems associated with structural components, and for determining the stiffness and damping properties of intercomponent connections. In general, Parameter Identification is improved when greater quantities of experimental data are available.

Lawrence, Charles↗

Cosmic-ray viscosity

The transport equation for cosmic rays scattered by magnetic field irregularities carried in a rarefied conducting fluid has been reexamined. To lowest order in the ratio, U/w, of flow speed to random particle speed, the analysis gives the standard equation first derived by Parker (1965), but additional terms are found in the next order of this ratio. One new term, which reflects viscous damping of fluid motions by the energetic-particle gas, describes both a change in the mean particle momentum and a spreading around the mean. Other new terms, which derive from accelerations of the fluid, describe inertial drift and energy changes. Although these effects are small, they are potentially important at shocks, because they are proportional to the square of velocity derivatives. In addition, they can be significant for the case of pure velocity shear, in which the adiabatic energy change is zero.

Earl, J. A.↗

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↗