Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Average 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

Analytical skin friction and heat transfer formula for compressible internal flows

An analytic, closed-form friction formula for turbulent, internal, compressible, fully developed flow was derived by extending the incompressible law-of-the-wall relation to compressible cases. The model is capable of analyzing heat transfer as a function of constant surface temperatures and surface roughness as well as analyzing adiabatic conditions. The formula reduces to Prandtl's law of friction for adiabatic, smooth, axisymmetric flow. In addition, the formula reduces to the Colebrook equation for incompressible, adiabatic, axisymmetric flow with various roughnesses. Comparisons with available experiments show that the model averages roughly 12.5 percent error for adiabatic flow and 18.5 percent error for flow involving heat transfer.

Dechant, Lawrence J.↗

Estimating the urban bias of surface shelter temperatures using upper-air and satellite data. Part 2: Estimation of the urban bias

A methodology is presented for estimating the urban bias of surface shelter temperatures due to the effect of the urban heat island. Multiple regression techniques were used to predict surface shelter temperatures based on the time period 1986-89 using upper-air data from the European Centre for Medium-Range Weather Forecasts (ECMWF) to represent the background climate, site-specific data to represent the local landscape, and satellite-derived data -- the normalized difference vegetation index (NDVI) and the Defense Meteorological Satellite Program (DMSP) nighttime brightness data -- to represent the urban and rural landscape. Local NDVI and DMSP values were calculated for each station using the mean NDVI and DMSP values from a 3 km x 3 km area centered over the given station. Regional NDVI and DMSP values were calculated to represent a typical rural value for each station using the mean NDVI and DMSP values from a 1 deg x 1 deg latitude-longitude area in which the given station was located. Models for the United States were then developed for monthly maximum, mean, and minimum temperatures using data from over 1000 stations in the U.S. Cooperative (COOP) Network and for monthly mean temperatures with data from over 1150 stations in the Global Historical Climate Network (GHCN). Local biases, or the differences between the model predictions using the observed NDVI and DMSP values, and the predictions using the background regional values were calculated and compared with the results of other research. The local or urban bias of U.S. temperatures, as derived from all U.S. stations (urban and rural) used in the models, averaged near 0.40 C for monthly minimum temperatures, near 0.25 C for monthly mean temperatures, and near 0.10 C for monthly maximum temperatures. The biases of monthly minimum temperatures for individual stations ranged from near -1.1 C for rural stations to 2.4 C for stations from the largest urban areas. The results of this study indicate minimal problems for global application once global NDVI and DMSP data become available.

Epperson, David L.↗

Banded Electron Structure Formation in the Inner Magnetosphere

Banded electron structures in energy-time spectrograms have been observed in the inner magnetosphere concurrent with a sudden relaxation of geomagnetic activity. In this study, the formation of these banded structures is considered with a global, bounce-averaged model of electron transport, and it is concluded that this structure is a natural occurrence when plasma sheet electrons are captured on closed drift paths near the Earth. These bands do not appear unless there is capture of plasma sheet electrons; convection along open drift paths making open pass around the Earth do not have time to develop this feature. The separation of high-energy bands from the injection population due to the preferential advection of the gradient-curvature drift creates spikes in the energy distribution, which overlap to form a series of bands in the energy spectrograms. The lowest band is the bulk of the injected population in the sub-key energy range. Using the Kp history for an observed banded structure event, a cloud of plasma sheet electrons is captured and the development of their distribution function is examined and discussed.

Liemohn, M. W.↗

Photochemistry of Saturn's Atmosphere: Hydrocarbon Chemistry and Comparisons with ISO Observations - 1

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Moses, Julianne I.↗

Hydrocarbon photochemistry in the upper atmosphere of Jupiter

The hydrocarbon photochemistry in the upper atmosphere of Jupiter is investigated using a one-dimensional, photochemical-diffusive, and diurnally averaged model. The important chemical cycles and pathways among the major species are outlined and a standard model for the North Equatorial Belt region is examined in detail. It is found that several traditionally dominant chemical pathways among the C and C2 species are replaced in importance by cycles involving C-C4 species. The pressure and altitude profiles of mixing ratios for several observable hydrocarbon species are compared with available ultraviolet- and infrared-derived abundances. The results of sensitivity studies on the standard model with respect to variations in eddy diffusion profile, solar flux, atomic hydrogen influx, latitude, temperature, and important chemical reaction rates are presented. Measured and calculated airglow emissions of He at 584 angstroms and H at 1216 angstroms are also used to provide some constraints on the range of model parameters. The relevance of the model results to the upcoming Galileo mission is briefly discussed. The model is subject to considerable improvement; there is a great need for laboratory measurements of basic reaction rates and photodissociation quantum yields, even for such simple species as methylacetylene and allene. Until such laboratory measurements exist there will be considerable uncertainty in the understanding of the C3 and higher hydrocarbons in the atmospheres of the jovian planets.

NASA Discipline Number 52-20↗

Dynamics of the Lunar Spin Axis

The evolution of the lunar spin axis is studied. Prior work has assumed that the inclination of the lunar orbit is constant and that the node regresses uniformly. This work takes into account the nonconstant inclination and nonuniform regression of the node as determined from averaged models of the motion of the lunar orbit. The resulting dynamics is considerably more rich, exhibiting additional resonances, period doubling and tripling, and chaos.

Wisdom, Jack↗

Hill Problem Analytical Theory to the Order Four. Application to the Computation of Frozen Orbits around Planetary Satellites

Frozen orbits of the Hill problem are determined in the double averaged problem, where short and long period terms are removed by means of Lie transforms. The computation of initial conditions of corresponding quasi periodic solutions in the non-averaged problem is straightforward for the perturbation method used provides the explicit equations of the transformation that connects the averaged and non-averaged models. A fourth order analytical theory reveals necessary for the accurate computation of quasi periodic, frozen orbits.

Lara, Martin↗

Transient Climate Effects of Large Impacts on Titan

Titan's thick atmosphere and volatile-rich surface cause it to respond to big impacts in a somewhat Earth-like manner. Here we construct a simple globally-averaged model that tracks the flow of energy through the environment in the weeks, years, and millenia after a big comet strikes Titan. The model Titan is endowed with 1.4 bars of N2 and 0.07 bars of CH4, methane lakes, a water ice crust, and enough methane underground to saturate the regolith to the surface. We find that a nominal Menrva impact is big enough to raise the surface temperature by approx. 80 K and to double the amount of methane in the atmosphere. The extra methane drizzles out of the atmosphere over hundreds of years. An upper-limit Menrva is just big enough to raise the surface to water's melting point. The putative Hotei impact (a possible 800-1200 km diameter basin, Soderblom et al., 2009) is big enough to raise the surface temperature to 350-400 K. Water rain must fall and global meltwaters might range between 50 m to more than a kilometer deep, depending on the details. Global meltwater oceans do not last more than a few decades or centuries at most, but are interesting to consider given Titan's organic wealth. Significant near-surface clathrate formation is possible as Titan cools but faces major kinetic barriers.

Titan↗

Big Impacts and Transient Oceans on Titan

We have studied the thermal consequences of very big impacts on Titan [1]. Titan's thick atmosphere and volatile-rich surface cause it to respond to big impacts in a somewhat Earth-like manner. Here we construct a simple globally-averaged model that tracks the flow of energy through the environment in the weeks, years, and millenia after a big comet strikes Titan. The model Titan is endowed with 1.4 bars of N2 and 0.07 bars of CH4, methane lakes, a water ice crust, and enough methane underground to saturate the regolith to the surface. We assume that half of the impact energy is immediately available to the atmosphere and surface while the other half is buried at the site of the crater and is unavailable on time scales of interest. The atmosphere and surface are treated as isothermal. We make the simplifying assumptions that the crust is everywhere as methane saturated as it was at the Huygens landing site, that the concentration of methane in the regolith is the same as it is at the surface, and that the crust is made of water ice. Heat flow into and out of the crust is approximated by step-functions. If the impact is great enough, ice melts. The meltwater oceans cool to the atmosphere conductively through an ice lid while at the base melting their way into the interior, driven down in part through Rayleigh-Taylor instabilities between the dense water and the warm ice. Topography, CO2, and hydrocarbons other than methane are ignored. Methane and ethane clathrate hydrates are discussed quantitatively but not fully incorporated into the model.

methane↗

Acceleration and Velocity Sensing from Measured Strain

A simple approach for computing acceleration and velocity of a structure from the strain is proposed in this study. First, deflection and slope of the structure are computed from the strain using a two-step theory. Frequencies of the structure are computed from the time histories of strain using a parameter estimation technique together with an autoregressive moving average model. From deflection, slope, and frequencies of the structure, acceleration and velocity of the structure can be obtained using the proposed approach. Simple harmonic motion is assumed for the acceleration computations, and the central difference equation with a linear autoregressive model is used for the computations of velocity. A cantilevered rectangular wing model is used to validate the simple approach. Quality of the computed deflection, acceleration, and velocity values are independent of the number of fibers. The central difference equation with a linear autoregressive model proposed in this study follows the target response with reasonable accuracy. Therefore, the handicap of the backward difference equation, phase shift, is successfully overcome.

shape sensing↗

Acceleration and Velocity Sensing from Measured Strain

A simple approach for computing acceleration and velocity of a structure from the strain is proposed in this study. First, deflection and slope of the structure are computed from the strain using a two-step theory. Frequencies of the structure are computed from the time histories of strain using a parameter estimation technique together with an Autoregressive Moving Average model. From deflection, slope, and frequencies of the structure, acceleration and velocity of the structure can be obtained using the proposed approach. shape sensing, fiber optic strain sensor, system equivalent reduction and expansion process.

shape sensing↗

Unsteady Aerodynamic Force Sensing from Measured Strain

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm. A cantilevered rectangular wing built and tested at the NASA Langley Research Center (Hampton, Virginia, USA) in 1959 is used to validate the simple approach. Unsteady aerodynamic forces as well as wing deflections, velocities, accelerations, and strains are computed using the CFL3D computational fluid dynamics (CFD) code and an MSC/NASTRAN code (MSC Software Corporation, Newport Beach, California, USA), and these CFL3D-based results are assumed as measured quantities. Based on the measured strains, wing deflections, velocities, accelerations, and aerodynamic forces are computed using the proposed approach. These computed deflections, velocities, accelerations, and unsteady aerodynamic forces are compared with the CFL3D/NASTRAN-based results. In general, computed aerodynamic forces based on the lifting surface theory in subsonic speeds are in good agreement with the target aerodynamic forces generated using CFL3D code with the Euler equation. Excellent aeroelastic responses are obtained even with unsteady strain data under the signal to noise ratio of -9.8dB. The deflections, velocities, and accelerations at each sensor location are independent of structural and aerodynamic models. Therefore, the distributed strain data together with the current proposed approaches can be used as distributed deflection, velocity, and acceleration sensors. This research demonstrates the feasibility of obtaining induced drag and lift forces through the use of distributed sensor technology with measured strain data. An active induced drag control system thus can be designed using the two computed aerodynamic forces, induced drag and lift, to improve the fuel efficiency of an aircraft. Interpolation elements between structural finite element grids and the CFD grids and centroids are successfully incorporated with the unsteady aeroelastic computation scheme. The most critical technology for the success of the proposed approach is the robust on-line parameter estimator, since the least-squares curve fitting method depends heavily on aeroelastic system frequencies and damping factors.

shape sensing↗

Scale Dependence of Cirrus Horizontal Heterogeneity Effects on TOA Measurements: MODIS Brightness Temperatures in the Thermal Infrared - Part I

This paper presents a study on the impact of cirrus cloud heterogeneities on MODIS simulated thermal infrared (TIR) brightness temperatures (BTs) at the top of the atmosphere (TOA) as a function of spatial resolution from 50 meters to 10 kilometers. A realistic 3-D (three-dimensional) cirrus field is generated by the 3DCLOUD model (average optical thickness of 1.4, cloudtop and base altitudes at 10 and 12 kilometers, respectively, consisting of aggregate column crystals of D (sub eff) equals 20 microns), and 3-D thermal infrared radiative transfer (RT) is simulated with the 3DMCPOL (3-D Monte Carlo Polarized) code. According to previous studies, differences between 3-D BT computed from a heterogenous pixel and 1-D (one-dimensional) RT computed from a homogeneous pixel are considered dependent at nadir on two effects: (i) the optical thickness horizontal heterogeneity leading to the plane-parallel homogeneous bias (PPHB); and the (ii) horizontal radiative transport (HRT) leading to the independent pixel approximation error (IPAE). A single but realistic cirrus case is simulated and, as expected, the PPHB mainly impacts the low-spatial resolution results (above approximately 250 meters), with averaged values of up to 5-7 K (thousand), while the IPAE mainly impacts the high-spatial resolution results (below approximately 250 meters) with average values of up to 1-2 K (thousand). A sensitivity study has been performed in order to extend these results to various cirrus optical thicknesses and heterogeneities by sampling the cirrus in several ranges of parameters. For four optical thickness classes and four optical heterogeneity classes, we have found that, for nadir observations, the spatial resolution at which the combination of PPHB and HRT effects is the smallest, falls between 100 and 250 meters. These spatial resolutions thus appear to be the best choice to retrieve cirrus optical properties with the smallest cloud heterogeneity-related total bias in the thermal infrared. For off-nadir observations, the average total effect is increased and the minimum is shifted to coarser spatial resolutions.

TOA↗

Unsteady Aerodynamic Force Sensing from Strain Data

A simple approach for computing unsteady aerodynamic forces from simulated measured strain data is proposed in this study. First, the deflection and slope of the structure are computed from the unsteady strain using the two-step approach. Velocities and accelerations of the structure are computed using the autoregressive moving average model, on-line parameter estimator, low-pass filter, and a least-squares curve fitting method together with analytical derivatives with respect to time. Finally, aerodynamic forces over the wing are computed using modal aerodynamic influence coefficient matrices, a rational function approximation, and a time-marching algorithm.

aerodynamic load sensing↗

Development of a two-dimensional zonally averaged statistical-dynamical model. III - The parameterization of the eddy fluxes of heat and moisture

A number of perpetual January simulations are carried out with a two-dimensional zonally averaged model employing various parameterizations of the eddy fluxes of heat (potential temperature) and moisture. The parameterizations are evaluated by comparing these results with the eddy fluxes calculated in a parallel simulation using a three-dimensional general circulation model with zonally symmetric forcing. The three-dimensional model's performance in turn is evaluated by comparing its results using realistic (nonsymmetric) boundary conditions with observations. Branscome's parameterization of the meridional eddy flux of heat and Leovy's parameterization of the meridional eddy flux of moisture simulate the seasonal and latitudinal variations of these fluxes reasonably well, while somewhat underestimating their magnitudes. New parameterizations of the vertical eddy fluxes are developed that take into account the enhancement of the eddy mixing slope in a growing baroclinic wave due to condensation, and also the effect of eddy fluctuations in relative humidity. The new parameterizations, when tested in the two-dimensional model, simulate the seasonal, latitudinal, and vertical variations of the vertical eddy fluxes quite well, when compared with the three-dimensional model, and only underestimate the magnitude of the fluxes by 10 to 20 percent.

Stone, Peter H.↗

Models of space averaged energetics of plates

The analysis of high frequency vibrations in plates is of particular interest in the study of structure borne noise in aircrafts. The current methods of analysis are either too expensive (finite element method) or may have a confidence band wider than desirable (Statistical Energy Analysis). An alternative technique to model the space and time averaged response of structural acoustics problems with enough detail to include all significant mechanisms of energy generation, transmission, and absorption is highly desirable. The focus of this paper is the development of a set of equations which govern the space and time averaged energy density in plates. To solve this equation, a new type of boundary value problem must be treated in terms of energy density variables using energy and intensity boundary conditions. A computer simulation verification study of the energy governing equation is performed. A finite element formulation of the new equations is also implemented and several test cases are analyzed and compared to analytical solutions.

Bouthier, O. M.↗

Interpretation and modeling of the averaged equations for a fluid-solid flow

A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Shen, Hayley H.↗

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗