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 523 records · Page 29

Entrainment into a stratocumulus layer with distributed radiative cooling

It is shown that the radiative cooling of a cloud layer strongly influences the turbulent flux profiles and the entrainment rate, and that the radiative cooling should be modeled as acting inside the turbulent layer. Numerical experiments demonstrate that a cloud-topped mixed-layer model, similar to that of Lilly (1968), is quite sensitive to delta (p) sub R, the depth of the radiatively cooled layer near cloud top. As delta (p) sub R increases, the model's sensitivity to the entrainment assumption is markedly heightened; for large delta (p) sub R the cloud top and cloud base rise as the entrainment parameter k is increased, while for small delta (p) sub R an increase in k has almost no effect. The model is most sensitive to delta (p) sub R for the cold-water, strong-divergence regime of greatest interest.

Randall, D. A.↗

A comparison of transport algorithms for premixed, laminar steady state flames

The effects of different methods of approximating multispecies transport phenomena in models of premixed, laminar, steady state flames were studied. Five approximation methods that span a wide range of computational complexity were developed. Identical data for individual species properties were used for each method. Each approximation method is employed in the numerical solution of a set of five H2-02-N2 flames. For each flame the computed species and temperature profiles, as well as the computed flame speeds, are found to be very nearly independent of the approximation method used. This does not indicate that transport phenomena are unimportant, but rather that the selection of the input values for the individual species transport properties is more important than the selection of the method used to approximate the multispecies transport. Based on these results, a sixth approximation method was developed that is computationally efficient and provides results extremely close to the most sophisticated and precise method used.

Coffee, T. P.↗

Viscous, resistive MHD stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one dimensional magnetohydrodynamic (MHD) quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds like numbers. Marginal stability curves, growth rates versus Reynolds like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result which appears general is that instability was found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Closed coronal structures. V - Gasdynamic models of flaring loops and comparison with SMM observations

A time-dependent one-dimensional code incorporating energy, momentum and mass conservation equations, and taking the entire solar atmospheric structure into account, is used to investigate the hydrodynamic response of confined magnetic structures to strong heating perturbations. Model calculation results are compared with flare observations which include the light curves of spectral lines formed over a wide range of coronal flare temperatures, as well as determinations of Doppler shifts for the high temperature plasma. It is shown that the numerical simulation predictions are in good overall agreement with the observed flare coronal plasma evolution, correctly reproducing the temporal profile of X-ray spectral lines and their relative intensities. The predicted upflow velocities support the interpretation of the blueshifts as due to evaporation of chromospheric material.

Peres, G.↗

Algorithm Estimates Microwave Water-Vapor Delay

Accuracy equals or exceeds conventional linear algorithms. "Profile" algorithm improved algorithm using water-vapor-radiometer data to produce estimates of microwave delays caused by water vapor in troposphere. Does not require site-specific and weather-dependent empirical parameters other than standard meteorological data, latitude, and altitude for use in conjunction with published standard atmospheric data. Basic premise of profile algorithm, wet-path delay approximated closely by solution to simplified version of nonlinear delay problem and generated numerically from each radiometer observation and simultaneous meteorological data.

Robinson, Steven E.↗

AFE base flow computations

Hypersonic wake flows behind the Aeroassist Flight Experiment (AFE) geometry are analyzed using two Navier-Stokes flow solvers. Many of the AFE wake features observed in ballistic-range shadowgraphs are simulated using a simple, two-dimensional semicylinder geometry at moderate angles of attack. At free-stream conditions corresponding to a Hypersonic Free Flight Facility (HFFF) AFE experiment, the three-dimensional base flow for the AFE geometry is computed using an ideal-gas, Navier-Stokes solver. The computed results agree reasonably well with the shadowgraphs taken at the HFFF. An ideal-gas and a nonequilibrium Navier-Stokes solver have been coupled and applied to the complete flow around the AFE vehicle at the free-stream conditions corresponding to a nomial trajectory point. Limitations of the coupled ideal-gas and nonequilibrium solution are discussed. The nonequilibrium base flow solution is analyzed for the wake radiation and the radiation profiles along various lines of sight are compared. Finally, the wake unsteadiness is predicted using experimental correlations and the numerical solutions. An adaptive grid code, SAGE, has been used in all the simulations to enhance the solution accuracy. The grid adaptation is found to be necessary in obtaining base flow solutions with accurate flow features.

Venkatapathy, Ethiraj↗

Time-dependent cosmic ray shocks with injection - A progress report

The study presents preliminary results of a two-fluid hydrodynamical simulation of cosmic-ray-modified shock structure by extending the steady-state calculations of Zank et al. (1991) to the time-dependent case. The present simulations model diffusive shock acceleration of the energetic particle component at planar shocks. A new model for the injection of cosmic ray particles is presented. The quadratic artificial viscosity algorithm is found to be numerically unstable and therefore unusable for strong cosmic ray shocks. There is visible ringing and overshoot in the various profiles. A large enhancement in the fluid density appears downstream of the subshock. This 'overcompression' results from the upstream fluid being compressed in the precursor, then traversing a strong shock with compression ratio near the adiabatic limit of four.

Donohue, D. J.↗

Solar convection zone dynamics and rotation

Helioseismology reveals that there appears to be a discrepancy between the differential rotation profile with the radius and latitude deduced from the inversion of the observed frequency splitting of p-modes and the rotation profile anticipated from various global simulations of convection in rotating shells. The problem may be caused by deductions in numerical simulations of convection which was considered near laminar flows. A high performance computing offers paths for studying the properties of such astrophysical turbulence. A range of approaches to study the basic dynamics of convection is reviewed, along with its ability to redistribute angular momentum in rotating systems. The results of recent three-dimensional simulations of turbulent compressible convection constrained by the effects of rotation are described for the cases in spherical shells and local area f-planes. The turbulence possesses intense vortex tubes with intricate interactions and instabilities. The theoretical studies reveal that the transition to turbulent states has contributed to significant changes in flow symmetries and transports, and that such a turbulent convection can drive substantial mean flows, which are distinctly different from those in which the convection is dominantly laminar.

Toomre, Juri↗

Molecular Basis of Mechano-Signal Transduction in Vascular Endothelial Cells

Simulated microgravity studies using a random positioning machine (RPM). One RPM machine has been built for us by Fokker Science in Netherland. Using the device, we have developed an in vitro system to examine the effect of simulated microgravity on osteoblastic bone cells. Using this system, we have carried out gene chip studies to determine the gene expression profiles of osteoblasts cultured under simulated microgravity conditions in comparison to static controls. From this study, we have identified numerous genes, some of which are expected ones inducing bone loss, but many of which are unexpected and unknown. These findings are being prepared for publications.

Jo, Hanjoong↗

Scale resolving simulations of the NASA Juncture Flow Model using the LAVA solver

Two distinct simulation methodologies: Delayed Detached Eddy Simulation (DDES) and stress based Wall-modelled Large Eddy Simulation (WMLES) are evaluated using structured overset curvilinear grids for the NASA juncture flow model. While both methodologies are shown to mitigate the primary shortcomings of steady state Reynolds Averaged Navier Stokes(RANS) simulations, several unresolved aspects are identified. Strong sensitivity to RANS-typegrid refinement is observed in the DDES with a substantial deterioration of the solution quality with increasing spatial resolution associated with deficiencies in the shielding function. Mean profiles for attached boundary layers on the fuselage show spurious inflections suggesting modelled stress depletion on finer grids. Lower numerical dissipation in terms of spatial discretization and time step size is seen to improve the solution quality, and the advantage of using RC and QCR2000 corrections in the underlying RANS closure is demonstrated for DDES on RANS-type meshes. Equilibrium Wall-Modelled LES used grids that resolved the tripping dots over the fuselage nose and the wing leading edge, consistent with the experiment setup. These simulations resulted in a cost-competitive approach compared to DDES on RANS-type grids. Although the agreement between WMLES predictions of first and second order single point statistics with experimental measurements is promising, a fundamental shortcoming is noted in terms of an overshoot in streamwise momentum in corner regions of the wing-fuselage juncture. This excess momentum subsequently delays the onset of separation,thereby resulting in an underprediction in length of the separation bubble. Some quantitative sensitivity to numerical discretization is observed; lowering of numerical dissipation shows better agreement with the experiment.

T^3↗

Improved Impact of Atmospheric Infrared Sounder (AIRS) Radiance Assimilation in Numerical Weather Prediction

Improvements to global and regional numerical weather prediction (NWP) have been demonstrated through assimilation of data from NASA s Atmospheric Infrared Sounder (AIRS). Current operational data assimilation systems use AIRS radiances, but impact on regional forecasts has been much smaller than for global forecasts. Retrieved profiles from AIRS contain much of the information that is contained in the radiances and may be able to reveal reasons for this reduced impact. Assimilating AIRS retrieved profiles in an identical analysis configuration to the radiances, tracking the quantity and quality of the assimilated data in each technique, and examining analysis increments and forecast impact from each data type can yield clues as to the reasons for the reduced impact. By doing this with regional scale models individual synoptic features (and the impact of AIRS on these features) can be more easily tracked. This project examines the assimilation of hyperspectral sounder data used in operational numerical weather prediction by comparing operational techniques used for AIRS radiances and research techniques used for AIRS retrieved profiles. Parallel versions of a configuration of the Weather Research and Forecasting (WRF) model with Gridpoint Statistical Interpolation (GSI) that mimics the analysis methodology, domain, and observational datasets for the regional North American Mesoscale (NAM) model run at the National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center (EMC) are run to examine the impact of each type of AIRS data set. The first configuration will assimilate the AIRS radiance data along with other conventional and satellite data using techniques implemented within the operational system; the second configuration will assimilate AIRS retrieved profiles instead of AIRS radiances in the same manner. Preliminary results of this study will be presented and focus on the analysis impact of the radiances and profiles for selected cases.

Zavodsky, Bradley↗

Assimilation of Microwave Observations in the Rainbands of Tropical Cyclones

We propose a novel Bayesian Monte Carlo Integration (BMCI) technique to retrieve the profiles of temperature, water vapor, and cloud liquid/ice water content from microwave cloudy measurements in the presence of tropical cyclones (TC). These retrievals then can either be directly used by meteorologists to analyze the structure of TCs or be assimilated into numerical models to provide accurate initial conditions for the NWP (Numerical Weather Prediction) models. The BMCI technique is applied to the data from the Advanced Technology Microwave Sounder (ATMS) onboard Suomi National Polar-orbiting Partnership (NPP) and Global Precipitation Measurement (GPM) Microwave Imager (GMI). The retrieved profiles are then assimilated into Hurricane WRF (Weather Research and Forecasting) using the GSI (Gridpoint Statistical Interpolation) data assimilation system.

Moradi, Isaac↗

Direct Numerical Simulation of Transitional and Turbulent Flow Over a Heated Flat Plate Using Finite-Difference Schemes

The work in this report was conducted at NASA Ames Research Center during the period from August 1993 to January 1995 deals with the direct numerical simulation of transitional and turbulent flow at low Mach numbers using high-order-accurate finite-difference techniques. A computation of transition to turbulence of the spatially-evolving boundary layer on a heated flat plate in the presence of relatively high freestream turbulence was performed. The geometry and flow conditions were chosen to match earlier experiments. The development of the momentum and thermal boundary layers was documented. Velocity and temperature profiles, as well as distributions of skin friction, surface heat transfer rate, Reynolds shear stress, and turbulent heat flux were shown to compare well with experiment. The numerical method used here can be applied to complex geometries in a straightforward manner.

Madavan, Nateri K.↗

Correlations Between Topex Altimeter and Temperature Profile Data in the Equatorial Pacific

Observations of sea-surface heights from satellite altimeters have been a major contributor to monitoring the evolution of the recent major El Nino event. Assimilating altimeter data directly into numerical models only takes advantage of the surface signature provided by the data. However, that surface signature is usually indicative of processes occurring at depth, especially in the equatorial Pacific. In order for assimilation schemes to make maximum use of surface data, it is helpful to have knowledge of how to best extend that data in the vertical to account for that variability. Toward that end, the vertical correlation structure between satellite-observed sea-surface heights and in situ temperature measurements is examined using TOPEX altimeter and TOGA TAO profile data for the period 1993-1998. This time period includes several states of the tropical Pacific, including a perturbed state from 1993-1994, a quiescent state encompassing 1995-1996, and the major ENSO event of 1997-1998. Analyses for each of these periods, as well as the total period will be presented. In addition to analyses at specific depth levels, analyses for particular density surfaces will also be presented.

Adamec, David↗

Calculation of turbulent boundary layers with heat transfer and pressure gradient utilizing a compressibility transformation. Part 2: Constant property turbulent boundary layer flow with simultaneous mass transfer and pressure gradient

An analysis of the incompressible turbulent boundary layer, developing under the combined effects of mass transfer and pressure gradient, is presented in this paper. A strip-integral method is employed whereby two of the three governing equations are obtained by integrating the combined momentum and continuity equation to 50 percent and 100 percent, respectively, of the boundary-layer height. The latter equation is the usual momentum-integral equation; the former equation requires specification of shear. Accordingly, Clauser's equilibrium eddy-viscosity law is assumed valid at this point. The third and final equation is obtained by specifying that Stevenson's velocity profiles apply throughout the domain of interest, from which a skin-friction law can be derived. Comparisons of the numerical results with the experiments of McQuaid, which include combined effects of variable pressure gradient and mass transfer, show good agreement.

Boccio, J.↗

Low-frequency flute instabilities of a bounded plasma column.

Derivation of exact solutions for unstable waves (called flute waves) which occur in a radially bounded plasma column at frequencies below the ion cyclotron frequency. Both analytical and numerical solutions are presented for the m = 1 and m = 2 azimuthal modes for a variety of radial electric field profiles. It is shown that the behavior of the flute waves can depend sensitively on the radial extent of the plasma column. Moreover, it is found that the m = 1 mode and the m = 2 mode do not respond in the same way to changes in the radial boundary position or in the electric field profile.

Rognlien, T. D.↗

The photochemistry of ammonia in the Jovian atmosphere

A reinvestigation of the ammonia photochemistry in the Jovian lower atmosphere above the cloud tops is conducted in connection with certain questions concerning the model proposed by Strobel (1973). Chemical reactions and their rate coefficients are considered along with the model atmosphere, solar fluxes, the continuity equations, and numerical details. Differences between the newmodel and Strobel's model are related to more moderate ammonia depletion relations. The height profiles of H, NH2, N2H3, N2H4, and NH3 concentrations are presented in a graph.

Prasad, S. S.↗

Computer program documentation D1FLTD to drive SINDA boundary nodes: User's guide

The thermal model correlation process begins when measured thermocouple data is available from the orbital flight tests of the shuttle. For this effort, it is necessary to convert some of the system improved numerical differencing analyzer (SINDA) diffusion or arithmetic nodes to boundary nodes and then drive these boundary nodes to the temperature profile of a flight measurement. An efficient way to provide this capability within the SINDA and OFT software systems is to provide a new SINDA routine, D1FLTD, for use in VARIABLES 1 of SINDA, to access the processed (word-addressable) orbital data reduction center flight data and store the appropriate measurement temperature in the desired SINDA temperature location. The ODRC flight data that is to be used for driving the boundary nodes must be assigned a logical unit number and must reside on a word-addressable file. The user must also provide two SINDA constants for the word positions of the first and last words of the temperature record for each measurement identifier (MID), i.e. each call to D1FLTD, used in the model. D1FLTD is then called from the VARIABLES 1 block to obtain the SINDA boundary node temperature for any MID on the file at any time point.

Damico, S. J.↗