Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - BOUNDARY LAYER - PRANDTL”

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 199 records · Page 11

Euler/Navier-Stokes calculations of transonic flow past fixed- and rotary-wing aircraft configurations

Computational fluid dynamics has an increasingly important role in the design and analysis of aircraft as computer hardware becomes faster and algorithms become more efficient. Progress is being made in two directions: more complex and realistic configurations are being treated and algorithms based on higher approximations to the complete Navier-Stokes equations are being developed. The literature indicates that linear panel methods can model detailed, realistic aircraft geometries in flow regimes where this approximation is valid. As algorithms including higher approximations to the Navier-Stokes equations are developed, computer resource requirements increase rapidly. Generation of suitable grids become more difficult and the number of grid points required to resolve flow features of interest increases. Recently, the development of large vector computers has enabled researchers to attempt more complex geometries with Euler and Navier-Stokes algorithms. The results of calculations for transonic flow about a typical transport and fighter wing-body configuration using thin layer Navier-Stokes equations are described along with flow about helicopter rotor blades using both Euler/Navier-Stokes equations.

Deese, J. E.↗

Effect of subjectively enhancing the initial vorticity field in a short-range numerical forecast

The paper tests a simple method to incorporate subsynoptic-scale information into the initial conditions of a forecast model. The amplitude, orientation and vertical structure of the vorticity field associated with the subsynoptic-scale disturbance are estimated subjectively; this is done using a combination of climatological and physical models together with all available data. The initial values of geopotential are then obtained in the vicinity of the subsynoptic-scale disturbance, as over the rest of the domain, through the balance equation. Two sets of initial conditions are tested in 12-hr forecasts using a six-layer version of a primitive equation forecast model (Anthes and Warner, 1978). The analysis schemes consisted of (1) a careful subjective analysis of radiosonde data with a high-resolution objective interpolation, and (2) the same preceding analysis but including data derived from the enhanced vorticity fields in the vicinity of the subsynoptic-scale disturbance. The results are discussed in some detail.

Fried, G. D.↗

Numerical study of a delta planform with multiple jets in ground effect

The flow past a 60-deg delta wing equipped with two thrust-reverser jets near the inboard trailing edge has been analyzed by numerical solution of the 3D thin-layer Navier-Stokes equations. An implicit, partially flux-split, approximately-factored Navier-Stokes solver coupled with a multiple grid embedding scheme has been adapted to this problem. Studies of the impact of numerical parameters (e.g., grid refinement and dissipation levels), and flow-field parameters such as the height of the delta wing above the ground plane and the jet size on the solution, were performed. Results of these numerical studies indicate some challenges in the accurate resolution of complex 3D free shear layers and jets. Nevertheless, flow features such as jet deformation and ground vortex formation observed in experimental flow visualizations are captured. Further, comparisons with experimental data confirm the ability to simulate the loss of wing-borne lift, commonly referred to 'suckdown, as the delta planform flies at slow speeds in close proximity to the ground. Detailed analysis of the numerical results has also given additional insight into the structure of the ground vortex and the mechanisms of lift loss.

Chawla, K.↗

Light transmission and reflection in perennially ice-covered Lake Hoare, Antarctica

We have investigated the transmission and albedo of the perennial ice cover on Lake Hoare, Antarctica. Our database consists of year-round measurements of the photosynthetically active radiation (400-700 nm) under the ice, measurements of the spatial variation of the under-ice light in midsummer, and spectrally resolved measurements from 400 to 700 nm of the albedo and transmission of the ice cover in early (November) and in midsummer (January). Our results show that the transmission decreases in the first part of summer, dropping by a factor of approximately 4 from November to January. We suggest that this is due to heating in the upper layers of the ice cover and the formation of Tyndall figures. Later in the summer when a significant liquid water fraction occurs within the ice cover, the transmission increases. In the fall when the ice cover freezes solid the transmission drops markedly. The spectrally resolved measurements from 400 to 700 nm show that approximately 2-5% of the incident light in this spectral region penetrates the 3.5-m thick ice cover. We have analyzed the spectral data using a two-stream scattering solution to the radiative transfer equation with three vertical layers in the ice cover. A surficial glaze of scattering ice 1 cm thick overlies a layer of sandy, bubbly ice about a meter thick, and below this is a thick layer of sand-free ice with bubbles. We find that the ice cover is virtually opaque at wavelengths longer than 800 nm. The net transmission of solar energy is approximately 2%. Significant changes in the thickness of the ice cover have been reported at Lake Hoare. These are due primarily to changes in the thickness of the bottom layer only. Because this layer is relatively clear, the effect on the transmission through the ice cover from these changes is less than would be predicted assuming a homogeneous ice cover.

Mckay, C. P.↗

Modified radiative transfer theory for a two-layer anisotropic random medium

The modified radiative transfer (MRT) equations which describe propagation and scattering of the electromagnetic field intensity in a layered anisotropic random medium are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The extinction matrices and scattering phase matrices for both mean and incoherent field intensities are obtained. Backscattering enhancement is observed due to the spatial coherent effect between upward and downward propagating waves. It also occurs for the half-space case, because of coupling between ordinary and extraordinary waves. The depolarization effect is predicted in the first-order renormalization approximation to the MRT equation.

Lee, J. K.↗

Reformulation of Oxide Growth Equations for Oxidation of Silicon Bond Coat in Environmental Barrier Coating Systems

The formulation of Deal and Grove for the oxidation of silicon and the growth of a silica layer is revisited to include the effects of a protective coating layer. The linear-parabolic equation of the form x_o^2+Ax_o=Bt is still applicable, but the expressions for the constants A and B must be modified to account for the coating. The modified expression for A is a linear function of the coating thickness. The expression for B is modified by the relative solubility of the outer coating layer. The modified expressions are applied to assess recently published experimental results for the oxidation of a silicon bond coat in an environmental barrier coating (EBC). Conclusions are drawn regarding the diffusivity and solubility of oxidants in the EBC top coating.

Sullivan, Roy M.↗

Modified radiative transfer theory for a two-layer random medium

Modified radiative transfer (MRT) equations appropriate for electromagnetic wave propagation in bounded random media are derived from the Bethe-Salpeter equation with the ladder approximation and the Dyson equation with the nonlinear approximation. The MRT equations are then solved with the first-order renormalization approximation to obtain analytical results for the backscattering cross sections of a two-layer random medium with arbitrary three-dimensional correlation functions. The coherent effects of the MRT theory are illustrated and comparisons are made with backscattering cross sections obtained with the first Born approximation to the wave equation.

Zuniga, M. A.↗

Theoretical study of multiple equilibria in simple axisymmetric tropical circulations

The possibility that the asymmetric part of the atmospheric circulation can possess multiple equilibrium states is examined using a two-layer axisymmetric model involving balance equations on an equatorial beta plane. Mountains are excluded from consideration and a Newtonian cooling formulation represents thermal forcing. A temperature maximum at 25 deg N is selected to simulate summer conditions in the Northern Hemisphere. Steady-state solutions obtained are investigated for stability with regard to first and second y-mode perturbations. A single stable mode is found, together with two other quasi-stable states. Attention is given to numerically modeling multiple equilibria in symmetric circulations, and one steady-state is determined for the two-layer model. A model employing primitive equations with the Boussinesq approximation is also examined, and it also furnishes only one steady state. The reasons for the lack of multiple steady-states as derived by the models are discussed.

Goswami, B. N.↗

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗

Universal reduced basis for the calibration of covariant energy density functionals

The reduced basis method is used to construct a “universal” basis of Dirac orbitals that may be applicable throughout the nuclear chart to calibrate covariant energy density functionals. Relative to the successful development of a reduced basis emulator for the nonrelativistic Schrödinger equation, the Dirac equation adds an extra layer of complexity due to the existence of negative energy states, which complicates building an efficient reduced basis. However, once this problem is mitigated, the resulting reduced basis is able to accurately and efficiently reproduce the high-fidelity model at a fraction of the computational cost. We are confident that the resulting reduced basis will serve as a foundational element in developing rapid and accurate emulators. In turn, these emulators will play a critical role in the Bayesian optimization of covariant energy density functionals.

Bayesian methods↗

Mode coupling in partial reflections from the ionosphere at vertical incidence.

The effect of mode coupling has been calculated for radio waves partially reflected from the lower ionosphere at vertical incidence. Calculations were made for parameters appropriate to the partial reflection experiments (Belrose, 1970). The four reflection coefficients (two direct and two coupled) are calculated for a plane discontinuity using a matrix perturbation analysis. The ratios of coupled to direct reflection coefficients are calculated for a general one-dimensional reflecting layer using the coupled equations of Clemmow and Heading (1954). At the highest electron density considered, 10,000 per cu cm, about 10% of the reflected extraordinary wave amplitude might be due to coupling from the incident ordinary wave. The effect is smaller at lower electron density.

Connolly, D. J.↗

A second-order closure analysis of turbulent diffusion flames

A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.

Varma, A. K.↗

Acoustic resonances and sound scattering by a shear layer

The reflection and transmission characteristics of plane waves scattered by a finite-thickness shear layer having a linear velocity profile and bounded by two otherwise uniform parallel flows is examined using the pressure perturbation equation solutions in the shear layer shown previously to be in terms of Whittaker M functions. In addition to the angle of plane wave incidence and the relative Mach number of the flows bounding the shear layer, it is found that the scattering properties of the shear layer depend crucially upon a parameter tau in such a manner that the case tau approaching zero characterizes the long wavelength properties of the layer and the case tau approaching infinity characterizes the short wavelength properties of the layer. In contrast to the region of ordinary reflection in the cases where the corresponding vortex sheet does not have a Brewster angle, the values of the reflection coefficient up to tau of 2 follow those of the vortex sheet; for the case for which the corresponding vortex sheet has a Brewster angle, the magnitude of the reflection coefficient may be sensitive to even small changes in tau in certain cases.

Koutsoyannis, S. P.↗

A numerical simulation of hypersonic viscous flow over arbitrary geometries at angle of attack

An implicit conservative, noniterative, finite-difference algorithm that predicts the supersonic, laminar or turbulent viscous flow about arbitrary geometries at large angles of attack is presented. The three-dimensional parabolized form of the thin-layer Navier-Stokes equations are written in generalized coordinates. These equations are solved using the delta form of the Beam-Warming implicit algorithm. Flow field simulations have been obtained for a blunt biconic with windward and leeward cuts and an X-24C lifting body for both laminar and turbulent flow at various Mach numbers and angles of attack. When compared with experiment or with previous theories, these computational predictions show good agreement.

Chaussee, D. S.↗

A numerical simulation of two-dimensional inlet flow fields

Inlet flow fields for airbreathing missiles are calculated by the adaptation of a two-dimensional computational method developed for the flow around airfoils. A supersonic free stream is assumed to allow the forebody calculation to be uncoupled from the inlet calculation. The inlet calculation employs an implicit, time-marching finite-difference procedure to solve the thin-layer Navier-Stokes equations formulated in body-fitted coordinates. Because the method can be used for a flow field with both subsonic and supersonic regions, it is applicable to subcritical as well as supercritical inlet operation.

Biringen, S.↗