Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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,243 records · Page 69

1DVAR Analysis of Temperature and Humidity using GPS Radio Occulation Data

The Global Positioning System (GPS) enables positioning with a very small receiver. The signals transmitted by the GPS satellites are sensitive to the atmosphere and can be used to perform soundings with the radio occultation technique (e.g., Kursinski et al., 1997). The GPS signal can be converted to refractivity N via the Abel transform. The refractivity can then be related to atmospheric pressure P, temperature T, and water vapor partial pressure P (sub w) the GPS measurement, (between 0.5 and 1.5 km), its self-calibration, and it's nearly all-weather capabilities make it a good candidate for use in data assimilation systems (DAS) and numerical weather prediction (NWP). In order to demonstrate its usefulness in a DAS or NWP system, a first step is to assess its impact oil the analysis. A one-dimensional variational off-line analysis (1DVAR), meaning the data are not assimilated 'In the 3D DAS, constitutes a starting approach to which further enhancements can be made. The chosen observable to be analyzed in this study is the refractivity. One way to extract temperature (humidity) from the refractivity, is to assume a humidity (temperature) profile. One variable may then be retrieved without any a priori information. The 1DVAR approach used here resolves the ambiguity problem raised in the interpretation of these data. It enables retrieving these two atmospheric variables at a reasonable computing cost.

Poli, Paul↗

Model of Atmospheric Links on Optical Communications from High Altitude

Optical communication links have the potential to solve many of the problems of current radio and microwave links to satellites and high-altitude aircraft. The higher frequency involved in optical systems allows for significantly greater signal bandwidth, and thus information transfer rate, in excess of 10 Gbps, and the highly directional nature of laser-based signals eliminates the need for frequency-division multiplexing seen in radio and microwave links today. The atmosphere, however, distorts an optical signal differently than a microwave signal. While the ionosphere is one of the most significant sources of noise and distortion in a microwave or radio signal, the lower atmosphere affects an optical signal more significantly. Refractive index fluctuations, primarily caused by changes in atmospheric temperature and density, distort the incoming signal in both deterministic and nondeterministic ways. Additionally, suspended particles, such as those in haze or rain, further corrupt the transmitted signal. To model many of the atmospheric effects on the propagating beam, we use simulations based on the beam-propagation method. This method, developed both for simulation of signals in waveguides and propagation in atmospheric turbulence, separates the propagation into a diffraction and refraction problem. The diffraction step is an exact solution, within the limits of numerical precision, to the problem of propagation in free space, and the refraction step models the refractive index variances over a segment of the propagation path. By applying refraction for a segment of the propagation path, then diffracting over that same segment, this method forms a good approximation to true propagation through the atmospheric medium. Iterating over small segments of the total propagation path gives a good approximation to the problem of propagation over the entire path. Parameters in this model, such as initial beam profile and atmospheric constants, are easily modified in a simulation such as this, which allows for the rapid analysis of different propagation scenarios. Therefore, this method allows the development of a near-optimal system design for a wide range of situations, typical of what would be seen in different atmospheric conditions over a receiving ground station. A simulation framework based upon this model was developed in FORTRAN, and for moderate grid sizes and propagation distances these simulations are computable in reasonable time on a standard workstation. This presentation will discuss results thus far.

Subich, Christopher↗

Numerical Study of the Turbulent Prandtl Number in Supersonic Plane-Channel Flow – the Effect of Thermal Boundary Conditions

Direct numerical simulation is used to investigate the turbulent Prandtl number Prt above cold (isothermal) and hot (adiabatic) walls in a family of supersonic channel flows. A range of Reynolds numbers and mean temperature/density variations is considered. The value of Prt away from the wall approaches 0.85 above both the isothermal and adiabatic walls, with the width of the Prt = 0.85 region increasing with wall-variable-based Reynolds number Reτw (the latter strongly affected by the thermal boundary condition). The variable, near-wall Prt profiles from both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗ wall scaling proposed by Huang et al., with the y∗-dependence weakly parameterized by Reτw.

Turbulence↗

The Crack-contact and the Free End Problem for a Strip Under Residual Stress

The plane problem for an infinite strip with two edge cracks under a given state of residual stress is considered. The residual stress is compressive near and at the surfaces and tensile in the interior of the strip. If the crack is deep enough to penetrate into the tensile zone, then the problem is one of crack-contact problem in which the depth of the contact area is an unknown which depends on the crack depth and the residual stress profile. The problem has applications to the static fatigue of glass plates and is solved for three typical residual stress profiles. In the limiting case of the crack crossing the entire plate thickness, the problem becomes a stress-free end problem for a semi-infinite strip under a given residual stress state away from the end. This is a typical stress diffusion problem in which decay behavior of the residual stress near and the nature of the normal displacement at the end of the semi-infinite strip are of special interest. For two typical residual stress states the solution is obtained, and some numerical results are given.

Bakioglu, M.↗

Parameterization of radiative flux profiles within layer clouds

This paper presents simple formulations for the vertical structure of solar and infrared radiative fluxes within a layer cloud overlying a partially reflective surface. The parameterized profile shapes are analytic, with governing parameters derived from bulk radiative properties and more physically based radiative transfer models. The bulk cloud, subcloud layer, and surface radiative properties are assumed to be known. The parameterizations are based on exponential functions of height, with decay scales related to cloud liquid water content. Although the results presented here are based on very simple assumptions about the cloud structure in the vertical, the method used is applicable to more general cases as well as to various other analytic and/or numerical radiative transfer calculations.

Hanson, Howard P.↗

Evolution of the solar nebula. I - Nonaxisymmetric structure during nebula formation

Numerical solutions of the equations of hydrodynamics, gravitation, and radiative transfer in three spatial dimensions are used to model the formation and time evolution of the early solar nebula in order to learn whether or not gravitational torques between nonaxisymmetric structures in the solar nebula can transport angular momentum rapidly enough to produce nebula clearing on astronomically indicated (10 to the 5 to 10 to the 7 yr) time scales. The models involve solutions for the collapse of spherical clouds with assumed initial density and rotation profiles onto protosuns of variable mass. Most of the models assume uniform initial density and rotation, and have variations in the initial parameters of cloud mass, cloud rotation rate, and protosun mass which are chosen to simulate a range of possible phases of early solar nebula evolution. The models show little tendency for directly forming small numbers of giant gaseous protoplanets through gaseous gravitational instability.

Boss, Alan P.↗

Computational Study of Chaotic and Ordered Solutions of the Kuramoto-Sivashinsky Equation

We report the results of extensive numerical experiments on the Kuramoto-Sivashinsky equation in the strongly chaotic regime as the viscosity parameter is decreased and increasingly more linearly unstable modes enter the dynamics. General initial conditions are used and evolving states do not assume odd-parity. A large number of numerical experiments are employed in order to obtain quantitative characteristics of the dynamics. We report on different routes to chaos and provide numerical evidence and construction of strange attractors with self-similar characteristics. As the 'viscosity' parameter decreases the dynamics becomes increasingly more complicated and chaotic. In particular it is found that regular behavior in the form of steady state or steady state traveling waves is supported amidst the time-dependent and irregular motions. We show that multimodal steady states emerge and are supported on decreasing windows in parameter space. In addition we invoke a self-similarity property of the equation, to show that these profiles are obtainable from global fixed point attractors of the Kuramoto-Sivashinsky equation at much larger values of the viscosity.

Smyrlis, Yiorgos S.↗

HST Studies of the Chromospheres, Wind, and Mass-Loss Rates of Cool Giant and Supergiant Stars

UV spectra of K-M giant and supergiant stars and of carbon stars have been acquired with the Goddard High Resolution Spectrograph (GHRS) on the Hubble Space Telescope (HST). These spectra have been used to measure chromospheric flow and turbulent velocities, study the acceleration of their stellar winds, acquire constraints on their outer atmospheric structure, and enable estimates of their mass-loss rates. Results from our observations of the giant stars Gamma Dra (K5 III hybrid), Alpha Tau (K5 III), Gamma Cru (M3.4 III), Mu Gem (M3 IIIab), and 30 Her (MG III), the supergiants Alpha Ori (M2 Iab) and Lambda Vel (K5 Ib), and the carbon stars TX Psc (NO; C6,2) and TW Hor (NO; C7,2) will be summarized and compared. The high resolution and wavelength accuracy of these data have allowed the direct measurement of the acceleration of the stellar winds in the chromospheres of several of these stars (from initial velocities of 3-9 km/s to upper velocities of 15-25 km/s) and of the chromospheric macroturbulence (-25-35 km/s). The high signal-to-noise and large dynamic range of these spectra have allowed the detection and identification of numerous new emission features, including weak C IV emission indicative of hot transition-region plasma in the non-coronal giant Alpha Tau, many new fluorescent lines of Fe II, and the first detection of fluorescent molecular hydrogen emission and of Ca II recombination lines in the UV spectrum of a giant star. The UV spectrum of two carbon stars have been studied with unprecedented resolution and reveal extraordinarily complicated Mg II lines nearly smothered by circumstellar absorptions. Finally, comparison of synthetic UV emission line profiles computed with the Lamers et al. (1987) Sobolev with Exact Integration (SEI) code with observations of chromospheric emission lines overlain with wind absorption features provides estimates of the mass-loss rates for four of these stars.

Carpenter, Kenneth G.↗

Mars Climate Orbiter's Investigation of the Atmosphere and Polar Caps

The Mars Climate Orbiter (MCO) is now on its way to Mars. It carries an atmospheric sounder whose observations will provide a continuous, global data set on weather and climate for a full Martian year. This paper describes the observation strategy and anticipated results from the Pressure Modulator Infrared Radiometer (PMIRR). PMIRR will measure vertical profiles of atmospheric infrared radiance in the 7 to 50 micron wavelength region extending from the surface of Mars to 80-km altitude. The observations have a vertical resolution of 5 km, or one-half the atmospheric scale height. From these radiance profiles we will retrieve profiles of atmospheric temperature, pressure, and the amounts of dust, condensates and water vapor. In addition, PMIRR will measure the radiative balance of the polar regions of Mars in an effort to better understand the short-term climate variability of the planet. The information obtained with PMIRR on MCO will be complementary to data obtained by the Thermal Emission Spectrometer (TES) and Radio Science (RS) experiments on the Mars Global Surveyor. A major emphasis of our research will be the assimilation of PMIRR data into numerical models of the Martian atmosphere. Assimilation schemes, of which several are currently in development, will permit the extension of measurements to spatial and temporal scales and to phenomena (e.g. winds) not observed directly by PMIRR.

McCleese, D. J.↗

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.↗

A surface temperature and moisture parameterization for use in mesoscale numerical models

A modified multi-level soil moisture and surface temperature model is presented for use as in defining lower boundary conditions in mesoscale weather models. Account is taken of the hydraulic and thermal diffusion properties of the soil, their variations with soil type, and the mixing ratio at the surface. Techniques are defined for integrating the surface input into the multi-level scheme. Sample simulation runs were performed with the modified model and the original model defined by Pielke, et al. (1977, 1981). The models were applied to regional weather forecasting over soils composed of sand and clay loam. The new form of the model avoided iterations necessary in the earlier version of the model and achieved convergence at reasonable profiles for surface temperature and moisture in regions where the earlier version of the model failed.

Tremback, C. J.↗

Numerical Study of Wave Propagation in a Non-Uniform Flow

The propagation of acoustic waves originating from cylindrical and spherical pulses, in a non-uniform mean flow, and in the presence of a reflecting wall is investigated by Hardin and Pope approach using compact approximation of spatial derivatives. The 2-D and 3-D stagnation flows and a flow around a cylinder are taken as prototypes of real world flows with strong gradients of mean pressure and velocity. The intensity and directivity of acoustic wave patterns appear to be quite different from the benchmark solutions obtained in a static environment for the same geometry. The physical reasons for amplification and weakening of sound are discussed in terms of dynamics of wave profile and redistribution of acoustic energy and its potential and kinetic components. For an acoustic wave in the flow around a cylinder, the observed mean acoustic pressure is approximately doubled (upstream pulse position) and halved (downstream pulse position) in comparison with the sound propagation in static ambient conditions.

Povitsky, Alex↗

Preliminary study of a wall stabilized constructed arc

An iterative, implicit, finite-difference numerical technique is described which is suitable for obtaining solutions to the governing equations for a gas flowing in an axially symmetric constricted-arc heater. The method is shown to provide adequate solutions for three cases of simple pipe flows found in the literature, and for flow in a constricted-arc heater by direct comparison with experimental data. The comparison with arc-heater data includes static pressure, arc voltage, and wall heat flux, all as a function of axial location, and a radial temperature profile at one axial station. The arc-heater data were taken with air as the test gas at a heater inlet pressure approximately 0.40 atm and two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 584 amperes.

Graves, R. A.↗

Infrared cooling rate calculations in operational general circulation models - Comparisons with benchmark computations

The performance of several parameterized models is described with respect to numerical prediction and climate research at GFDL, NCAR, and GISS. The radiation codes of the models were compared to benchmark calculations and other codes for the intercomparison of radiation codes in climate models (ICRCCM). Cooling rates and fluxes calculated from the models are examined in terms of their application to established general circulation models (GCMs) from the three research institutions. The newest radiation parameterization techniques show the most significant agreement with the benchmark line-by-line (LBL) results. The LBL cooling rates correspond to cooling rate profiles from the models, but the parameterization of the water vapor continuum demonstrates uncertain results. These uncertainties affect the understanding of some lower tropospheric cooling, and therefore more accurate parameterization of the water vapor continuum, as well as the weaker absorption bands of CO2 and O3 is recommended.

Kiehl, J. T.↗

Aeroloads and secondary flows in a transonic mixed flow turbine stage

A numerical simulation of a transonic mixed flow turbine stage has been carried out using an average passage Navier-Stokes analysis. The mixed flow turbine stage considered here consists of a transonic nozzle vane and a highly loaded rotor. The simulation was run at the design pressure ratio and is assessed by comparing results with those of an established throughflow design system. The three-dimensional aerodynamic loads are studied as well as the development and migration of secondary flows and their contribution to the total pressure loss. The numerical results indicate that strong passage vortices develop in the nozzle vane, mix out quickly, and have little impact on the rotor flow. The rotor is highly loaded near the leading edge. Within the rotor passage, strong spanwise flows and other secondary flows exist along with the tip leakage vortex. The rotor exit loss distribution is similar in character to that found in radial inflow turbines. The secondary flows and non-uniform work extraction also tend to significantly redistribute a non-uniform inlet total temperature profile by the exit of the stage.

Kirtley, K. R.↗

CALIPSO Satellite Lidar Identification Of Elevated Dust Over Australia Compared With Air Quality Model PM60 Forecasts

Global measurements of the vertical distribution of clouds and aerosols have been recorded by the lidar on board the CALIPSO (Cloud Aerosol Lidar Infrared Pathfinder Satellite Observations) satellite since June 2006. Such extensive, height-resolved measurements provide a rare and valuable opportunity for developing, testing and validating various atmospheric models, including global climate, numerical weather prediction, chemical transport and air quality models. Here we report on the initial results of an investigation into the performance of the Australian Air Quality Forecast System (AAQFS) model in forecasting the distribution of elevated dust over the Australian region. The model forecasts of PM60 dust distribution are compared with the CALIPSO lidar Vertical Feature Mask (VFM) data product. The VFM classifies contiguous atmospheric regions of enhanced backscatter as either cloud or aerosols. Aerosols are further classified into six subtypes. By comparing forecast PM60 concentration profiles to the spatial distribution of dust reported in the CALIPSO VFM, we can assess the model s ability to predict the occurrence and the vertical and horizontal extents of dust events within the study area.

Young, Stuart A.↗

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The planned microgravity experiment Growth and Extinction Limit (GEL) is designed to study the material flammability and extinction of thermally thick solid with an emphasis to quantify the effect of conductive heat loss into the solid interior by measuring the in-depth temperature profile. To obtain a more detailed understanding of the extinction phenomena, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is expanded to account for the excessive heat loss into the solid beyond the steady burning state. This is important because for thick solids, this is important because the burning period can have a long and slow solid transient with a quasi-steady gas -phase. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper and extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The classical U-shaped flammability boundary in the oxygen vs. stretch rate plot now becomes a series of U-shaped boundaries as a function of Ψ. For a given oxygen ambient, extinction boundary is a function of Ψ. The boundary consists of a blowoff branch and a radiative quenching branch. This is similar to the steady burning flammability boundary (Ψ =0) in the oxygen vs stretch rate plot. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, and reaction rate are presented along the extinction boundaries. In particular, we note that the non-dimensional gas reaction rate is the ratio of diffusion reaction time to the reaction diffusion time, so it has the classical Damhkoler number character. By examining the maximum non-dimensional reaction rates, we observe that along the extinction boundaries, they are approximate constant across the entire range of stretch rates (from quenching to blow-off). These critical values are of the order unity but may vary in different oxygen ambient. The implication will be discussed in the presentation together with the near-limit flame structures.

Chengyao Li↗

A global model of the earth's ionosphere for use in space applications

A general expression is derived for the F layer electron density profile as a function of latitude and longitude for that part of the earth which is in direct sunlight including dawn and dusk. Furthermore, the derived model is extended to encompass the night-time ionosphere. The expressions allow determination by standard means of the range correction for arbitrary ray path directions. It is also shown that the naive application of the Chapman ionospheric model entails range correction errors which for low elevation angles (less than 20 deg) and large solar zenith angles (40 deg) cannot be tolerated. Numerical calculations are displayed showing the dependence of the range correction on the pertinent parameters.

Von Roos, O. H.↗