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,387 records · Page 77

The decay of isotropic turbulence in a rapidly rotating frame

A direct numerical simulation of the decay of initially isotropic turbulence in a rapidly rotating frame was conducted. This 128 x 128 x 128 simulation was completed for a Reynolds number Re sub lambda = 15.3 and a Rossby number Ro sub lambda = 0.07 based on the initial turbulent kinetic energy and Taylor microscale. The numerical results indicate that the turbulence remains essentially isotropic during the major part of the decay (i.e., beyond the point where the turbulent kinetic energy has decayed to less than 10 percent of its initial value). The rapid rotation has the primary effect of shutting off the energy transfer so that the turbulence dissipation (and hence the rate of decay of the turbulent kinetic energy) is substantially reduced. Consequently, the anisotropy tensor remains essentially unchanged while the energy spectrum undergoes a nearly linear viscous decay (the same results that are predicted by Rapid Distortion Theory which is only formally valid for much shorter elapsed times. Surprisingly, no Taylor-Proudman reorganization of the flow to a two-dimensional state is observed. The implications that these results have on turbulence modeling are discussed briefly along with prospective future research.

Speziale, C. G.↗

Shock breakout in SN 1987A

Detailed models for the UV burst from SN 1987A are calculated using a two-temperature radiation hydrodynamics code. The sequence of events which occur during shock breakout is described, including the formation of a radiative precursor, the disappearance of the shock, and the growth of a thin, dense shell above the photosphere which gives rise to a new viscous, isothermal shock. Continuum spectra are calculated under the assumption that the emergent radiation is a dilute blackbody with a color temperature equal to the radiation temperature at the thermalization depth. This color temperature is 2-3 times higher than the effective temperature. Light-travel-time corrections have been included in the light curves and the spectra. Quantities most relevant to accurate modeling of the fluorescing ring around SN 1987A are highlighted.

Ensman, Lisa↗

Analysis of Aerospike Plume Induced Base-Heating Environment

Computational analysis is conducted to study the effect of an aerospike engine plume on X-33 base-heating environment during ascent flight. To properly account for the effect of forebody and aftbody flowfield such as shocks and to allow for potential plume-induced flow-separation, thermo-flowfield of trajectory points is computed. The computational methodology is based on a three-dimensional finite-difference, viscous flow, chemically reacting, pressure-base computational fluid dynamics formulation, and a three-dimensional, finite-volume, spectral-line based weighted-sum-of-gray-gases radiation absorption model computational heat transfer formulation. The predicted convective and radiative base-heat fluxes are presented.

Wang, Ten-See↗

Unsteady Simulation of the Viscous Flow About a V-22 Rotor and Wing in Hover

Results of an unsteady thin-layer Navier-Stokes simulation of a 0.658-scale V-22 rotor and wing configuration in hover are presented. All geometric components of the flapped-wing and rotor test rig, including rotor blades, are accurately modeled. Rotor motion and rotor/airframe interference effects are simulated directly using moving body overset grid methods. Tiltrotor hover aerodynamics are visualized via unsteady particle trace images. Wing download predictive ability is demonstrated. Simulation results are compared with experimental data.

Meakin, Robert L.↗

Interfacial solitary waves in a two-fluid medium

Solitary waves are generated on the interface between two immiscible liquids with free upper surface; their behavior is generally consistent with that predicted by the Korteweg-de Vries equation. Inviscid theory expanded to include capillary effects predicts unchanged wave speed and narrowed wave shape to first order in amplitude regardless of upper boundary condition; but measured speeds are down 8% and narrowing is substantially greater at very small amplitudes. The speed difference appears consistent with an interfacial viscous boundary layer. The critical depth ratio separating the elevation and depression modes of the stable gravity solitary wave agrees with prediction in the inviscid limit. The design of a programmed, variable flux wave maker for efficiently generating interfacial waves is discussed.

Walker, L. R.↗

Particle motion in atmospheric boundary layers of Mars and Earth

To study the eolian mechanics of saltating particles, both an experimental investigation of the flow field around a model crater in an atmospheric boundary layer wind tunnel and numerical solutions of the two- and three-dimensional equations of motion of a single particle under the influence of a turbulent boundary layer were conducted. Two-dimensional particle motion was calculated for flow near the surfaces of both Earth and Mars. For the case of Earth both a turbulent boundary layer with a viscous sublayer and one without were calculated. For the case of Mars it was only necessary to calculate turbulent boundary layer flow with a laminar sublayer because of the low values of friction Reynolds number; however, it was necessary to include the effects of slip flow on a particle caused by the rarefied Martian atmosphere. In the equations of motion the lift force functions were developed to act on a single particle only in the laminar sublayer or a corresponding small region of high shear near the surface for a fully turbulent boundary layer. The lift force functions were developed from the analytical work by Saffman concerning the lift force acting on a particle in simple shear flow.

White, B. R.↗

The stability of motion of satellites with cavities partially filled with liquid

The stability and time dependent motion of a spinning satellite, simulated by a rigid body with a cavity partially filled with liquid is examined. The problem formulation, consisting of the boundary-value problem for the liquid and moment equations for the entire system is presented. Because of large Reynold's numbers involved, viscosity effects are negligible everywhere except for a thin boundary layer near the wetted surface. Using a boundary-layer analysis, the effect of the boundary layer is replaced by modified boundary conditions for the liquid. The solution of the differential equations for the inviscid problem is solved in closed form. A semi-analytical numerical solution of the inviscid equations subject to the viscous boundary condition has proved unsucessful.

Nayfeh, A. H.↗

Electron paramagnetic resonance studies of slowly tumbling vanadyl spin probes in nematic liquid crystals

An analysis of EPR line shapes by the method of Polnaszek, Bruno, and Freed is made for slowly tumbling vanadyl spin probes in viscous nematic liquid crystals. The use of typical vanadyl complexes as spin probes for nematic liquid crystals is shown to simplify the theoretical analysis and the subsequent interpretation. Rotational correlation times tau and orientational ordering parameters S sub Z where slow tumbling effects are expected to be observed in vanadyl EPR spectra are indicated in a plot. Analysis of the inertial effects on the probe reorientation, which are induced by slowly fluctuating torque components of the local solvent structure, yield quantitative values for tau and S sub Z. The weakly ordered probe VOAA is in the slow tumbling region and displays these inertial effects throughout the nematic range of BEPC and Phase V. VOAA exhibits different reorientation behavior near the isotropic-nematic transition temperature than that displayed far below this transition temperature.

Bruno, G. V.↗

Application of an implicit technique to the shock-layer flow around general bodies

The Euler and the parabolic Navier-Stokes equations are solved by an iterative implicit scheme along the body axis for angles of incidence up to 30 deg. A complete solution, including both the inviscid and viscous flows, can be obtained at stations normal to the axis. The technique is developed from a comparatively simple linearization procedure and has an option to iterate between body stations for higher accuracy. Reliable procedures have been introduced to account for the effect of the axial pressure gradient and to adjust the step increment for a given convergence requirement. A blunt-nosed conical and a Shuttle-Orbiter-like configuration were studied extensively and are compared with available laminar boundary-layer solutions and experimental data.

Li, C. P.↗

Measurements of the boundary layer and near wake of a supercritical airfoil at cruise conditions

The flow behavior within the upper-surface boundary layer and near wake of a supercritical airfoil operating at cruise conditions is discussed. Experimental results obtained from wind-tunnel tests are presented which provide a more detailed description of the flow in these regions than was previously available. Mean streamwise velocity profiles measured by pitot-pressure-probe and laser-velocimeter techniques were found to be in excellent agreement. Other mean-flow properties obtained by the laser-velocimeter technique were the local flow angles in the viscous layers and the static pressures at the edges of the boundary layer and wake. The data set also includes measurements of the turbulence intensity and turbulent Reynolds-stress distributions as obtained by the laser-velocimeter technique. To assess the effects of the shock wave, a less extensive set of measurements was realized at a subcritical test condition. The two test conditions (Mach number at free-stream conditions = 0.72, airfoil section lift coefficient = 0.76 and Mach number of free-stream conditions = 0.5, airfoil section lift coefficient = 0.75) provide a good test for state-of-the-art prediction methods because the upper-surface-boundary layer is separated just upstream of the trailing edge in both cases.

Johnson, D. A.↗

Direct simulation of stagnation streamline flow for hypersonic reentry

The direct simulation Monte Carlo (DSMC) method is applied to simulate one-dimensional flow along the stagnation streamline. The freestream conditions considered are those encountered by the nose region of the Space Shuttle Orbiter during the hypersonic reentry. The range of altitudes (75 to 92 km) considered in the present calculations covers continuum to the less rarefied portion of the transition flow regime. The calculations account for thermal as well as the chemical nonequilibrium effects. The attention is focused on the flow structure along the stagnation streamline for different specified shock locations at a given altitude. The effect of shock location on the stagnation-point gas composition and heat transfer is analyzed for specific altitude; consequently, these results are appropriate for hemispherical bodies with different nose radii. Finally, comparison is made between the present DSMC results and the viscous shock-layer calculations, which show good agreement at the lowest altitude.

Dogra, Virendra K.↗

A patched-grid algorithm for complex configurations directed towards the F-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 flowfield is shown.

Thomas, James L.↗

On vortical flow and sound generation

The generation of an acoustic field by an unsteady viscous vortical flow is investigated analytically, with a focus on the case where the reference scale length of the flow is much smaller than that of the acoustic field. The expansion schemes, governing equations, and far-field behavior are described; the effect of pressure and density fluctuations on the next-order far-field solution is explored; and the method of matched asymptotics is applied to construct the leading-order solution of the induced acoustic field. Particular attention is given to the sound generated by turbulence and by slender vortex filaments.

Ting, LU↗

Application of a patched-grid algorithm to the F/A-18 forebody-leading-edge extension configuration

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/A-18 forebody-Ieading-edge extension 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 flowfield is shown.

Leading edge extension↗

Transonic swirling nozzle flow

A numerical model of viscous transonic swirling flow in axisymmetric nozzles is developed. MacCormack's implicit Gauss-Seidel method is applied to the thin-layer Navier-Stokes equations in transformed coordinates. Numerical results are compared with experimental data to validate the method. The effect of swirl and viscosity on nozzle performance are demonstrated by examining wall pressures, Mach contours, and integral parameters.

Keith, Theo G., Jr.↗

In-flight surface-flow measurements on a subsonic transport high-lift flap system

As part of a multiphased program for subsonic transport high-lift flight research, flight tests were conducted on the Transport Systems Research Vehicle (B737-100 aircraft) at the NASA Langley Research Center, to obtain detailed flow characteristics of the high-lift flap system for correlation with computational and wind-tunnel investigations. Pressure distributions, skin friction, and flow-visualization measurements were made on a triple-slotted flap system for a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.36). Experimental test results are given for representative flap settings indicating flow separation on the fore-flap element for the largest flap deflection. Comparisons of the in-flight flow measurements were made with predictions from available viscous multielement computational methods modified with simple-sweep theory. Computational results overpredicted the experimentally measured pressures, particularly in the case involving separation of the fore lap, indicating the need for better modeling of confluent boundary layers and three-dimensional sweep effects.

Yip, Long P.↗

Agglomeration multigrid for viscous turbulent flows

Agglomeration multigrid, which has been demonstrated as an efficient and automatic technique for the solution of the Euler equations on unstructured meshes, is extended to viscous turbulent flows. For diffusion terms, coarse grid discretizations are not possible, and more accurate grid transfer operators are required as well. A Galerkin coarse grid operator construction and an implicit prolongation operator are proposed. Their suitability is evaluated by examining their effect on the solution of Laplace's equation. The resulting strategy is employed to solve the Reynolds-averaged Navier-Stokes equations for aerodynamic flows. Convergence rates comparable to those obtained by a previously developed non-nested mesh multigrid approach are demonstrated, and suggestions for further improvements are given.

Mavriplis, D. J.↗

A New Approach to Measure Contact Angle and Evaporation Rate with Flow Visualization in a Sessile Drop

The contact angle and the spreading process of sessile droplet are very crucial in many technological processes, such as painting and coating, material processing, film-cooling applications, lubrication, and boiling. Additionally, as it is well known that the surface free energy of polymers cannot be directly, measured for their elastic and viscous restraints. The measurements of liquid contact angle on the polymer surfaces become extremely important to evaluate the surface free energy of polymers through indirect methods linked with the contact angle data. Due to the occurrence of liquid evaporation is inevitable, the effects of evaporation on the contact angle and the spreading become very important for more complete understanding of these processes. It is of interest to note that evaporation can induce Marangoni-Benard convection in sessile drops. However, the impacts of the inside convection on the wetting and spreading processes are not clear. The experimental methods used by previous investigators cannot simultaneously measure the spreading process and visualize the convection inside. Based on the laser shadowgraphic system used by the present author, a very simple optical procedure has been developed to measure the contact angle, the spreading speed, the evaporation rate, and to visualize inside convection of a sessile drop simultaneously. Two CCD cameras were used to synchronously record the real-time diameter of the sessile drop, which is essential for determination of both spreading speed and evaporation rate, and the shadowgraphic image magnified by the sessile drop acting as a thin plano-convex lens. From the shadowgraph, the inside convection of the drop can be observed if any and the image outer diameter, which linked to the drop profile, can be measured. Simple equations have been derived to calculate the drop profile, including the instantaneous contact angle, height, and volume of the sessile drop, as well as the evaporation rate. The influence of the inside convection on the wetting and spreading processes can be figured out through comparison of the drop profiles with and without inside convection when the sessile drop is placed at different evaporation conditions.

Zhang, Nengli↗