Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large grids”

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,099 records · Page 61

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Balakumar, P.↗

Evaluation of Model Microphysics within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is approx 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were approx 0.25 m/s too slow, while the SBU-YLIN was 0.25 to 0.5 m/s too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗

Evaluation of Model Microphysics Within Precipitation Bands of Extratropical Cyclones

Recent studies evaluating the bulk microphysical schemes (BMPs) within cloud resolving models (CRMs) have indicated large uncertainties and errors in the amount and size distributions of snow and cloud ice aloft. The snow prediction is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33‐km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Coldseason Precipitation Experiment (GCPEx) experiment, as well as a few years (15 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku-band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. A Particle Size and Velocity (PARSIVEL) disdrometer was also used to measure the surface size distribution and fall speeds of snow at SBNY. For the 15 cases at SBNY, the WSM6, Morrison (MORR), Thompson (THOM2), and Stony Brook (SBU-YLIN) BMPs were validated. A non-spherical snow assumption (THOM2 and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU-YLIN schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM2 is ~0.25 meters per second too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM2, and MORR schemes were ~0.25 meters per second too slow, while the SBU-YLIN was 0.25 to 0.5 meters per second too fast. Overall, the BMPs simulate a size distribution close to the observed for D < 4 mm in the dendritic, plates, and mixed habit periods. The model BMPs underestimate the size distribution when large aggregates were observed. For D > 6 mm in the dendrites, side planes, and mixed habit periods, the BMPs are likely not simulating enough aggregation to create a larger size distribution, although the MORR (double moment) scheme seemed to perform best. These SBNY results will be compared with some results from GCPEx for a warm frontal snow band observed at 18 February 2012.

Colle, Brian A.↗

TECHNO-ECONOMIC ANALYSIS OF A CRYOGENIC FLUX CAPACITOR FOR GRID STORAGE

The Cryogenic Flux Capacitor (CFC) is a cold, dense fluid storage core with integrated design features that afford the designer flexibility and provide new possibilities for the storage and discharge of energy. The stored energy, in this case, is represented by hydrogen physically bonded within the nanoscale pores within the aerogel composite blanket material, and the process of bonding or debonding is governed by principles of physical adsorption (physisorption) and thermodynamics. The large surface area afforded by the nanoporous aerogel (~1,000 m2/g) allows for storage densities close to, or in some cases exceeding, that of normal boiling point liquids. Its performance easily exceeds what can be achieved via ambient temperature, high-pressure gas storage for an equivalent volume. CFC storage is predicted to be easily scalable, is constructed from readily available commercial materials, lends itself to a range of pressure applications, and is geometry insensitive.

Joshua Schmitt↗

Optimized silicon solar cells for space exploration power systems

A program is described aimed at designing and fabricating improved silicon solar cells for a range of missions extending from 0.1 to 15 astronomical units (Mercury to Jupiter). For missions lying inside Earth radius (Mercury, Venus) the major cell property required was very low series resistance, allowing high curve fill factor to be maintained at the higher intensities. For the Mercury mission, the temperature of the cell had to be kept low. This was achieved by reflecting more of the incident sunlight by use of large area front contacts. For the outer missions (Mars, the Asteroid belts and Jupiter) the formation of a Schottky barrier at the back contact had to be avoided (by use of a P+ layer under the back contact) and the excess leakage current of the PN junction had to be reduced. Optimum grid patterns were also derived and used for these missions.

Iles, P. A.↗

Effect of an oscillating flow direction on leading edge heat transfer

An experimental investigation was conducted to examine the effect of a periodic variation in the angle of attack on heat transfer at the leading edge of a gas turbine blade. A circular cylinder was used as a large-scale model of the leading edge region. The cylinder was placed in a wind tunnel and was oscillated rotationally about its axis. The incident flow Reynolds number and the Strouhal number of oscillation were chosen to model an actual turbine condition. Incident turbulence levels up to 4.9 percent were produced by grids placed upstream of the cylinder. The transfer rate was measured using a mass transfer technique and heat transfer rates inferred from the results. A direct comparison of the unsteady and steady results indicate that the effect is dependent on the Strouhal number, turbulence level, and the turbulence length scale, but that the largest observed effect was only a 10 percent augmentation at the nominal stagnation position.

Marziale, M. L.↗

Structure of turbulence in the presence of uniform shear

The structure of the vorticity field in homogeneous turbulent shear flow is analyzed using a database generated by direct numerical solution of the unsteady Navier-Stokes equations with up to 128x128x128 grid points. For two Reynolds numbers, the probability distribution of the inclination angle of the vorticity vectors, two-point correlations of the velocity and vorticity fields, and the instantaneous vorticity vectors and vortex lines in three-dimensional space were examined. It is shown that homogeneous turbulent shear flow contains a large number of horseshoe vortices. These vortices are most often found in planes inclined at 45 deg to the flow direction and are formed from the roll-up of sheets of transverse vorticity. These findings and similar results obtained in turbulent channel flow, lead to the conclusion that hairpin vortices are the characteristic vortical structures in all turbulent shear flows.

Moin, P.↗

Numerical simulation of the incompressible internal flow through a tilting disk valve

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady and unsteady flow calculations are performed by solving the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Chang, I-Dee↗

Simulation of blood flow through an artificial heart

A numerical simulation of the incompressible viscous flow through a prosthetic tilting disk heart valve is presented in order to demonstrate the current capability to model unsteady flows with moving boundaries. Both steady state and unsteady flow calculations are done by solving the incompressible Navier-Stokes equations in 3-D generalized curvilinear coordinates. In order to handle the moving boundary problems, the chimera grid embedding scheme which decomposes a complex computational domain into several simple subdomains is used. An algebraic turbulence model for internal flows is incorporated to reach the physiological values of Reynolds number. Good agreement is obtained between the numerical results and experimental measurements. It is found that the tilting disk valve causes large regions of separated flow, and regions of high shear.

Kiris, Cetin↗

Application of an unsteady Navier-Stokes solver to transonic turbine design

This study presents a numerical evaluation of the performance of the first stage of a new-generation turbine design. The numerical method solves the two-dimensional Navier-Stokes equations using a system of patched grids. Three-dimensional effects of stream-tube contraction are also modeled. The study focuses on the effects of axial gap variation on the unsteady rotor-stator interactions and on stage performance. Results are presented for three different axial gaps. The results indicate that the unsteady interactions can be very large in this design. These interactions affect not only the stage efficiency but also substantially alter the time-averaged features of the flow. In particular, for the case of the smallest axial gap, it was found that there was an unsteady shock on the stator suction surface which spanned the gap region and impinged upon the moving rotor airfoils.

Rangwalla, Akil A.↗

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗

An object-oriented approach for parallel self adaptive mesh refinement on block structured grids

Self-adaptive mesh refinement dynamically matches the computational demands of a solver for partial differential equations to the activity in the application's domain. In this paper we present two C++ class libraries, P++ and AMR++, which significantly simplify the development of sophisticated adaptive mesh refinement codes on (massively) parallel distributed memory architectures. The development is based on our previous research in this area. The C++ class libraries provide abstractions to separate the issues of developing parallel adaptive mesh refinement applications into those of parallelism, abstracted by P++, and adaptive mesh refinement, abstracted by AMR++. P++ is a parallel array class library to permit efficient development of architecture independent codes for structured grid applications, and AMR++ provides support for self-adaptive mesh refinement on block-structured grids of rectangular non-overlapping blocks. Using these libraries, the application programmers' work is greatly simplified to primarily specifying the serial single grid application and obtaining the parallel and self-adaptive mesh refinement code with minimal effort. Initial results for simple singular perturbation problems solved by self-adaptive multilevel techniques (FAC, AFAC), being implemented on the basis of prototypes of the P++/AMR++ environment, are presented. Singular perturbation problems frequently arise in large applications, e.g. in the area of computational fluid dynamics. They usually have solutions with layers which require adaptive mesh refinement and fast basic solvers in order to be resolved efficiently.

Lemke, Max↗

A Comparison of Seasonal and Interannual Variability of Soil Dust Aerosols Over the Atlantic Ocean as Inferred by the Toms AI and AVHRR AOT Retrievals

The seasonal cycle and interannual variability of two estimates of soil (or 'mineral') dust aerosols are compared: Advanced Very High Resolution Radiometer (AVHRR) aerosol optical thickness (AOT) and Total Ozone Mapping Spectrometer (TOMS) aerosol index (AI), Both data sets, comprising more than a decade of global, daily images, are commonly used to evaluate aerosol transport models. The present comparison is based upon monthly averages, constructed from daily images of each data set for the period between 1984 and 1990, a period that excludes contamination from volcanic eruptions. The comparison focuses upon the Northern Hemisphere subtropical Atlantic Ocean, where soil dust aerosols make the largest contribution to the aerosol load, and are assumed to dominate the variability of each data set. While each retrieval is sensitive to a different aerosol radiative property - absorption for the TOMS AI versus reflectance for the AVHRR AOT - the seasonal cycles of dust loading implied by each retrieval are consistent, if seasonal variations in the height of the aerosol layer are taken into account when interpreting the TOMS AI. On interannual time scales, the correlation is low at most locations. It is suggested that the poor interannual correlation is at least partly a consequence of data availability. When the monthly averages are constructed using only days common to both data sets, the correlation is substantially increased: this consistency suggests that both TOMS and AVHRR accurately measure the aerosol load in any given scene. However, the two retrievals have only a few days in common per month so that these restricted monthly averages have a large uncertainty. Calculations suggest that at least 7 to 10 daily images are needed to estimate reliably the average dust load during any particular month, a threshold that is rarely satisfied by the AVHRR AOT due to the presence of clouds in the domain. By rebinning each data set onto a coarser grid, the availability of the AVHRR AOT is increased during any particular month, along with its interannual correlation with the TOMS AI The latter easily exceeds the sampling threshold due to its greater ability to infer the aerosol load in the presence of clouds. Whether the TOMS AI should be regarded as a more reliable indicator of interannual variability depends upon the extent of contamination by sub-pixel clouds.

Cakmur, R. V.↗

Three-Dimensional Effects in Modeling of Dual-Mode Scramjets

A numerical investigation of an experimental dual-mode scramjet configuration is performed. Both experimental and numerical results indicate significant upstream interaction for this case. Several computational cases are examined: these include the use of jet-to-jet symmetry and entire half-duct modeling. Grid convergence, turbulence modeling, and wall temperature effects are studied in terms of wall pressure predictions and flow-field characteristics. Wall pressure comparisons between CFD and experiment show fair agreement for the jet-to-jet case. However, further computations of the entire half-duct show the development of a large sidewall separation zone extending much further upstream than the separation zone at the duct centerline. This sidewall separation is the dominant feature in the CFD-generated flowfield but is not evident in the experimental data, resulting in a unfavorable comparison between CFD and experimental data. Current work aimed at resolving this issue and at further understanding asymmetric flow-structures in dual-mode flow-fields is discussed.

Rodriguez, C. G.↗

Discrete Adjoint-Based Design Optimization of Unsteady Turbulent Flows on Dynamic Unstructured Grids

An adjoint-based methodology for design optimization of unsteady turbulent flows on dynamic unstructured grids is described. The implementation relies on an existing unsteady three-dimensional unstructured grid solver capable of dynamic mesh simulations and discrete adjoint capabilities previously developed for steady flows. The discrete equations for the primal and adjoint systems are presented for the backward-difference family of time-integration schemes on both static and dynamic grids. The consistency of sensitivity derivatives is established via comparisons with complex-variable computations. The current work is believed to be the first verified implementation of an adjoint-based optimization methodology for the true time-dependent formulation of the Navier-Stokes equations in a practical computational code. Large-scale shape optimizations are demonstrated for turbulent flows over a tiltrotor geometry and a simulated aeroelastic motion of a fighter jet.

Nielsen, Eric J.↗

Discovery of Nine Gamma-Ray Pulsars in Fermi-Lat Data Using a New Blind Search Method

We report the discovery of nine previously unknown gamma-ray pulsars in a blind search of data from the Fermi Large Area Telescope (LAT). The pulsars were found with a novel hierarchical search method originally developed for detecting continuous gravitational waves from rapidly rotating neutron stars. Designed to find isolated pulsars spinning at up to kHz frequencies, the new method is computationally efficient, and incorporates several advances, including a metric-based gridding of the search parameter space (frequency, frequency derivative and sky location) and the use of photon probability weights. The nine pulsars have spin frequencies between 3 and 12 Hz, and characteristic ages ranging from 17 kyr to 3 Myr. Two of them, PSRs Jl803-2149 and J2111+4606, are young and energetic Galactic-plane pulsars (spin-down power above 6 x 10(exp 35) ergs per second and ages below 100 kyr). The seven remaining pulsars, PSRs J0106+4855, J010622+3749, Jl620-4927, Jl746-3239, J2028+3332,J2030+4415, J2139+4716, are older and less energetic; two of them are located at higher Galactic latitudes (|b| greater than 10 degrees). PSR J0106+4855 has the largest characteristic age (3 Myr) and the smallest surface magnetic field (2x 10(exp 11)G) of all LAT blind-search pulsars. PSR J2139+4716 has the lowest spin-down power (3 x l0(exp 33) erg per second) among all non-recycled gamma-ray pulsars ever found. Despite extensive multi-frequency observations, only PSR J0106+4855 has detectable pulsations in the radio band. The other eight pulsars belong to the increasing population of radio-quiet gamma-ray pulsars.

Celik-Tinmaz, Ozlem↗

Multi-Angle Implementation of Atmospheric Correction for MODIS (MAIAC). Part 3: Atmospheric Correction

This paper describes the atmospheric correction (AC) component of the Multi-Angle Implementation of Atmospheric Correction algorithm (MAIAC) which introduces a new way to compute parameters of the Ross-Thick Li-Sparse (RTLS) Bi-directional reflectance distribution function (BRDF), spectral surface albedo and bidirectional reflectance factors (BRF) from satellite measurements obtained by the Moderate Resolution Imaging Spectroradiometer (MODIS). MAIAC uses a time series and spatial analysis for cloud detection, aerosol retrievals and atmospheric correction. It implements a moving window of up to 16 days of MODIS data gridded to 1 km resolution in a selected projection. The RTLS parameters are computed directly by fitting the cloud-free MODIS top of atmosphere (TOA) reflectance data stored in the processing queue. The RTLS retrieval is applied when the land surface is stable or changes slowly. In case of rapid or large magnitude change (as for instance caused by disturbance), MAIAC follows the MODIS operational BRDF/albedo algorithm and uses a scaling approach where the BRDF shape is assumed stable but its magnitude is adjusted based on the latest single measurement. To assess the stability of the surface, MAIAC features a change detection algorithm which analyzes relative change of reflectance in the Red and NIR bands during the accumulation period. To adjust for the reflectance variability with the sun-observer geometry and allow comparison among different days (view geometries), the BRFs are normalized to the fixed view geometry using the RTLS model. An empirical analysis of MODIS data suggests that the RTLS inversion remains robust when the relative change of geometry-normalized reflectance stays below 15%. This first of two papers introduces the algorithm, a second, companion paper illustrates its potential by analyzing MODIS data over a tropical rainforest and assessing errors and uncertainties of MAIAC compared to conventional MODIS products.

Lyapustin, A.↗

Exploratory Climate Data Visualization and Analysis Using DV3D and UVCDAT

Earth system scientists are being inundated by an explosion of data generated by ever-increasing resolution in both global models and remote sensors. Advanced tools for accessing, analyzing, and visualizing very large and complex climate data are required to maintain rapid progress in Earth system research. To meet this need, NASA, in collaboration with the Ultra-scale Visualization Climate Data Analysis Tools (UVCOAT) consortium, is developing exploratory climate data analysis and visualization tools which provide data analysis capabilities for the Earth System Grid (ESG). This paper describes DV3D, a UV-COAT package that enables exploratory analysis of climate simulation and observation datasets. OV3D provides user-friendly interfaces for visualization and analysis of climate data at a level appropriate for scientists. It features workflow inte rfaces, interactive 40 data exploration, hyperwall and stereo visualization, automated provenance generation, and parallel task execution. DV30's integration with CDAT's climate data management system (COMS) and other climate data analysis tools provides a wide range of high performance climate data analysis operations. DV3D expands the scientists' toolbox by incorporating a suite of rich new exploratory visualization and analysis methods for addressing the complexity of climate datasets.

Maxwell, Thomas↗