Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

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 127 records · Page 7

Development of a Unified SHOC-EDMF PBL Parameterization in the NASA GEOS Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies.

SHOC-EDMF↗

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗

A Unified SHOC-MF Boundary Layer Parameterization in the NASA Goddard Earth Observing System Model

We present results from a novel extension of the boundary layer parameterization in the NASA Goddard Earth Observing System (GEOS) model. An Eddy Diffusivity-Mass Flux (EDMF) multi-plume updraft scheme was previously implemented in GEOS. Here we replace the K-profile eddy diffusivity scheme of Lock et. al. (2000) with the Simplified Higher Order Closure (SHOC; Bogenschutz and Krueger, 2013). The scheme prognoses turbulent kinetic energy, and derives buoyancy flux, cloud fraction and condensate from a trivariate double-gaussian joint PDF. The original SHOC scheme is modified in two significant respects: (i) contributions from the EDMF updrafts are included in estimates of all higher order moments used to determine the PDF; and (ii) the second and third moments of total water are made prognostic variables, with source terms associated with EDMF updrafts. The predicted higher order moments as well as updraft and environmental properties are compared with large eddy simulations in several case studies. Results indicate that many higher order moments are underestimated when the mass flux contribution is neglected.

Nathan Arnold↗

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherents

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

bondline thickness↗

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherends

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

X-ray computed tomography↗

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

Siegel, R.↗

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

Polarization radiation in the planetary atmosphere delimited by a heterogeneous diffusely reflecting surface

Spatial frequency characteristics (SFC) and the scattering functions were studied in the two cases of a uniform horizontal layer with absolutely black bottom, and an isolated layer. The mathematical model for these examples describes the horizontal heterogeneities in a light field with regard to radiation polarization in a three dimensional planar atmosphere, delimited by a heterogeneous surface with diffuse reflection. The perturbation method was used to obtain vector transfer equations which correspond to the linear and nonlinear systems of polarization radiation transfer. The boundary value tasks for the vector transfer equation that is a parametric set and one dimensional are satisfied by the SFC of the nonlinear system, and are expressed through the SFC of linear approximation. As a consequence of the developed theory, formulas were obtained for analytical calculation of albedo in solving the task of dissemination of polarization radiation in the planetary atmosphere with uniform Lambert bottom.

Strelkov, S. A.↗

Seasonal and latitudinal variations of gravity wave-driven fluctuations in OH nightglow

A model which incorporates extended OH-layer emission and gravity-wave dynamics with eddy diffusivities is set forth to analyze seasonal and latitudinal variations in the OH nightglow. The 2D analysis is aimed at determining the nature of the relationship between the oscillation in intensity about a mean intensity and the oscillation in temperature of the emission region about a mean temperature. Competing eddy diffusivities modify the local values and altitudes of maximum wave amplitude, so long-period characterization of the relationship is imprecise. Seasonal trends in the relationship are noted for periods of several hrs or less that are primarily related to seasonal changes in the mesopause undisturbed temperature. Latitudinal trends are shown to be less pronounced, and trends that exist at intermediate periods are difficult to isolate due to interference effects.

Hickey, M. P.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are a relatively new and promising tool for todays aerospace technology challenges. This paper uses two stress-omega full Reynolds stress models to evaluate challenging flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSG/LRR full second-moment Reynolds stress models have been implemented into the FUN3D (Fully Unstructured Navier-Stokes Three Dimensional) unstructured Navier-Stokes code and are evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST-V) two-equation turbulence models.

reynolds averaging↗

Modeling Radar Scatter from Icy and Young Rough Lunar Craters

For lunar orbital synthetic aperture radars, such as the Chandrayaan Mini-RF operating at S- band (13-cm) wavelength and the Lunar Reconnaissance Orbiter Mini-RF operating at S- band and X-band (3-cm) wavelengths, it is important to understand the radar backscattering characteristics of the icy and young, rough craters. Assuming a mixing model consisting of diffuse and quasi-specular scattering components, we have modeled the opposite-sense circular (OC) and same-sense circular (SC) backscattering characteristics. The specular component, consisting of only OC echoes, represents the echoes from the surface and subsurface layers that are oriented perpendicular to the radar's line-of-sight. The diffuse component, consisting of both SC and OC echoes, represents the echoes associated with either rocks or ice. Also, diffuse echoes have backscatter that is proportional to the cosine of the incidence angle. We modeled how these two (specular and diffuse) radar scattering components could be modulated by factors such as surface roughness associated with young craters. We also modeled how ice radar scattering components could be modulated by a thin regolith covering, and/or by the situation where ice occupies small patches within a larger radar pixel. We tested this modeling by examining 4 nonpolar craters and 12 polar craters using LRO Mini-RF data. Results indicate that icy and young rough craters can be distinguished based upon their SC enhancements (Alpha) and OC enhancements (Gamma). In addition, we also examined the craters that have unusual circular polarization ratios (CPRs) that likely result from a double bounce mode of scattering. Blocky fresh craters, icy craters, and craters exhibiting double bounce scattering can be separated based on the values of Alpha, Gamma, the ratio of Alpha/Gamma and the weighted sum of Alpha and Gamma.

lunar craters↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D

Full seven-equation Reynolds stress turbulence models are promising tools for today’s aerospace technology challenges. This paper examines two such models for computing challenging turbulent flows including shock-wave boundary layer interactions, separation and mixing layers. The Wilcox and the SSG/LRR full second-moment Reynolds stress models have been implemented into the FUN3D (Fully Unstructured Navier-Stokes Three Dimensional) unstructured Navier-Stokes code and were evaluated for four problems: a transonic two-dimensional diffuser, a supersonic axisymmetric compression corner, a compressible planar shear layer, and a subsonic axisymmetric jet. Simulation results are compared with experimental data and results computed using the more commonly used Spalart-Allmaras (SA) one-equation and the Menter Shear Stress Transport (SST-V) two-equation turbulence models.

reynolds averaging↗

Simulated interannual variability of the Greenland Sea deep water formation and its connection to surface forcing

A fully prognostic Arctic ice-ocean model is used to study the interannual variability of deepwater formation in the Greenland Sea Gyre based on the simulations for the Arctic ice-ocean system for the period 1955 and 1960 - 1985. The model uses monthly climatology for thermodynamic forcing components (such as air temperature and cloudiness), together with constant annual net precipitation and river runoff. The daily wind forcing is derived from analyzed sea level air pressures from the National Center for Atmospheric Research (NCAR). In summary, the model shows that the occurence of deep convection in the Greenland Sea Gyre is controlled by the extensive Fram Strait ice export and/or local wind conditions in the Greenland Sea. In the latter case the weakening of the local wind curl allows the Polar Front to move eastward. The movement of the Polar Front causes adverse ice conditions, often together with much larger than normal ice export from the Arctic, such as in 1968, which can block convection in the gyre. The density difference between upper and lower layers is investigated as an indication of water mass formation through convection, occurring as strong diffusion in the model. The model-simulated density difference between the average top 100 m and deep levels reveals that the period 1960 - 1985 had only a few distinct years with weak stratification, and, especially, the model predicts no deep convection since the nid-1970s. The common factor for the years of the weakest decrease of the model-predicted heat content of the upper 2000 m which can, to a high degree, be explained by local heat loss.

Haekkinen, Sirpa↗

High-Temperature Oxidation Behavior of Iridium-Rhenium Alloys

The life-limiting mechanism for radiation-cooled rockets made from iridium-coated rhenium (Ir/Re) is the diffusion of Re into the Ir layer and the subsequent oxidation of the resulting Ir-Re alloy from the inner surface. In a previous study, a life model for Ir/Re rockets was developed. It incorporated Ir-Re diffusion and oxidation data to predict chamber lifetimes as a function of temperature and oxygen partial pressure. Oxidation testing at 1540 deg C suggested that a 20-wt percent Re concentration at the inner wall surface should be established as the failure criterion. The present study was performed to better define Ir-oxidation behavior as a function of Re concentration and to supplement the data base for the life model. Samples ranging from pure Ir to Ir-40 wt percent Re (Ir-40Re) were tested at 1500 deg C, in two different oxygen environments. There were indications that the oxidation rate of the Ir-Re alloy increased significantly when it went from a single-phase solid solution to a two-phase mixture, as was suggested in previous work. However, because of testing anomalies in this study, there were not enough dependable oxidation data to definitively raise the Ir/Re rocket failure criterion from 20-wt percent Re to a Re concentration corresponding to entry into the two-phase region.

Reed, Brian D.↗

An Eulerian transport/transformation/removal model for SO2 and sulfate. I - Model development. II - Model calculation of SO(x) transport in the eastern United States

A three-dimensional, time dependent Eulerian atmospheric SO2 and sulfate transport/transformation/removal model is described and applied to the eastern U.S. The model was developed in anticipation of increased input to the atmospheric sulfur content by coal-burning power plants in the near future and is intended as an aid in identifying sources of SO2. The Eulerian transport model incorporates functions for chemical transformations, dry deposition, spatial topographical variations, spatial and temporal variations of mixing layer heights, the wind field, eddy diffusivities, deposition velocities, and temperature and water concentration profiles. Attention is also given to the SO2 photochemical oxidation mechanism and rates. Results from a 72-hr simulation of SO(x) transport over the eastern U.S., using actual 1974 meteorological data, illustrate the model's capability to depict interactions between emissions, transport, chemistry and removal. Concentration distributions are demonstrated to have significant spatial and temporal variations.

Carmichael, G. R.↗

Wet model of Saturn's ionosphere: Water from the rings

Current theoretical models of Saturn's ionosphere are difficult to reconcile with the ionospheric electron density profiles obtained from the Pioneer and Voyager radio occultation observations and the large diurnal variation of maximum ionospheric electron density deduced from studies of Saturn lightning discharges. A model of Saturn's ionosphere is proposed in which water plays a major role as a minor constituent present by virtue of downward diffusion from an external source. This model of the Saturn ionosphere is a classical 'F2' type layer resulting from the photodissociative production of H(+) from H2 and rapid chemical loss due to a series of charge exchange reactions with water. A planet-wide influx of about 4x10 to the 7th power molecules/sec/sq cm of water from the rings is consistent with the observed ionospheric electron densities and estimates of influx due to micrometeoride bombardment of the rings. An enhanced influx of water occurs at latitudes (-38 deg, +44 deg) magnetically connected to the inner edge of Saturn's B ring which results from an electromagnetic erosion process contributing substantially to the (local) upper atmosphere water content. Present day influx at these latitudes is possibly as large as 2x10 to the 9th power molecules/sec/sq cm.

Connerney, J. E. P.↗

Application of the pressure sensitive paint technique to steady and unsteady flow

Pressure sensitive paint is a newly-developed optical measurement technique with which one can get a continuous pressure distribution in much shorter time and lower cost than a conventional pressure tap measurement. However, most of the current pressure sensitive paint applications are restricted to steady pressure measurement at high speeds because of the small signal-to-noise ratio at low speed and a slow response to pressure changes. In the present study, three phases of work have been completed to extend the application of the pressure sensitive paint technique to low-speed testing and to investigate the applicability of the paint technique to unsteady flow. First the measurement system using a commercially available PtOEP/GP-197 pressure sensitive paint was established and applied to impinging jet measurements. An in-situ calibration using only five pressure tap data points was applied and the results showed good repeatability and good agreement with conventional pressure tap measurements on the whole painted area. The overall measurement accuracy in these experiments was found to be within 0.1 psi. The pressure sensitive paint technique was then applied to low-speed wind tunnel tests using a 60 deg delta wing model with leading edge blowing slots. The technical problems encountered in low-speed testing were resolved by using a high grade CCD camera and applying corrections to improve the measurement accuracy. Even at 35 m/s, the paint data not only agreed well with conventional pressure tap measurements but also clearly showed the suction region generated by the leading edge vortices. The vortex breakdown was also detected at alpha=30 deg. It was found that a pressure difference of 0.2 psi was required for a quantitative pressure measurement in this experiment and that temperature control or a parallel temperature measurement is necessary if thermal uniformity does not hold on the model. Finally, the pressure sensitive paint was applied to a periodically changing pressure field with a 12.8s time period. A simple first-order pole model was applied to deal with the phase lag of the paint. The unsteady pressure estimated from the time-changing pressure sensitive paint data agreed well with the pressure transducer data in regions of higher pressure and showed the possibility of extending the technique to unsteady pressure measurements. However, the model still needs further refinement based on the physics of the oxygen diffusion into the paint layer and the oxygen quenching on the paint luminescence.

Shimbo, Y.↗