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 541 records · Page 30

Three-dimensional inviscid analysis of radial-turbine flow and a limited comparison with experimental data

The three-dimensional inviscid DENTON code is used to analyze flow through a radial-inflow turbine rotor. Experimental data from the rotor are compared with analytical results obtained by using the code. The experimental data available for comparison are the radial distributions of circumferentially averaged values of absolute flow angle and total pressure downstream of the rotor exit. The computed rotor-exit flow angles are generally underturned relative to the experimental values, which reflect the boundary-layer separation at the trailing edge and the development of wakes downstream of the rotor. The experimental rotor is designed for a higher-than-optimum work factor of 1.126 resulting in a nonoptimum positive incidence and causing a region of rapid flow adjustment and large velocity gradients. For this experimental rotor, the computed radial distribution of rotor-exit to turbine-inlet total pressure ratios are underpredicted due to the errors in the finite-difference approximations in the regions of rapid flow adjustment, and due to using the relatively coarser grids in the middle of the blade region where the flow passage is highly three-dimensional. Additional results obtained from the three-dimensional inviscid computation are also presented, but without comparison due to the lack of experimental data. These include quasi-secondary velocity vectors on cross-channel surfaces, velocity components on the meridional and blade-to-blade surfaces, and blade surface loading diagrams. Computed results show the evolution of a passage vortex and large streamline deviations from the computational streamwise grid lines. Experience gained from applying the code to a radial turbine geometry is also discussed.

Choo, Y. K.↗

Temporal Large-Eddy Simulation

In 1999, Stolz and Adams unveiled a subgrid-scale model for LES based upon approximately inverting (defiltering) the spatial grid-filter operator and termed .the approximate deconvolution model (ADM). Subsequently, the utility and accuracy of the ADM were demonstrated in a posteriori analyses of flows as diverse as incompressible plane-channel flow and supersonic compression-ramp flow. In a prelude to the current paper, a parameterized temporal ADM (TADM) was developed and demonstrated in both a priori and a posteriori analyses for forced, viscous Burger's flow. The development of a time-filtered variant of the ADM was motivated-primarily by the desire for a unifying theoretical and computational context to encompass direct numerical simulation (DNS), large-eddy simulation (LES), and Reynolds averaged Navier-Stokes simulation (RANS). The resultant methodology was termed temporal LES (TLES). To permit exploration of the parameter space, however, previous analyses of the TADM were restricted to Burger's flow, and it has remained to demonstrate the TADM and TLES methodology for three-dimensional flow. For several reasons, plane-channel flow presents an ideal test case for the TADM. Among these reasons, channel flow is anisotropic, yet it lends itself to highly efficient and accurate spectral numerical methods. Moreover, channel-flow has been investigated extensively by DNS, and a highly accurate data base of Moser et.al. exists. In the present paper, we develop a fully anisotropic TADM model and demonstrate its utility in simulating incompressible plane-channel flow at nominal values of Re(sub tau) = 180 and Re(sub tau) = 590 by the TLES method. The TADM model is shown to perform nearly as well as the ADM at equivalent resolution, thereby establishing TLES as a viable alternative to LES. Moreover, as the current model is suboptimal is some respects, there is considerable room to improve TLES.

Pruett, C. D.↗

Using MERRA Gridded Innovations for Quantifying Uncertainties in Analysis Fields and Diagnosing Observing System Inhomogeneities

MERRA is a NASA reanalysis for the satellite era using a major new version of the Goddard Earth Observing System Data Assimilation System Version 5 (GEOS-5). The project focuses on historical analyses of the hydrological cycle on a broad range of weather and climate time scales and places the NASA EOS suite of observations in a climate context. The characterization of uncertainty in reanalysis fields is a commonly requested feature by users of such data. While intercomparison with reference data sets is common practice for ascertaining the realism of the datasets, such studies typically are restricted to long term climatological statistics and seldom provide state dependent measures of the uncertainties involved. In principle, variational data assimilation algorithms have the ability of producing error estimates for the analysis variables (typically surface pressure, winds, temperature, moisture and ozone) consistent with the assumed background and observation error statistics. However, these "perceived error estimates" are expensive to obtain and are limited by the somewhat simplistic errors assumed in the algorithm. The observation minus forecast residuals (innovations) by-product of any assimilation system constitutes a powerful tool for estimating the systematic and random errors in the analysis fields. Unfortunately, such data is usually not readily available with reanalysis products, often requiring the tedious decoding of large datasets and not so-user friendly file formats. With MERRA we have introduced a gridded version of the observations/innovations used in the assimilation process, using the same grid and data formats as the regular datasets. Such dataset empowers the user with the ability of conveniently performing observing system related analysis and error estimates. The scope of this dataset will be briefly described. We will present a systematic analysis of MERRA innovation time series for the conventional observing system, including maximum-likelihood estimates of background and observation errors, as well as global bias estimates. Starting with the joint PDF of innovations and analysis increments at observation locations we propose a technique for diagnosing bias among the observing systems, and document how these contextual biases have evolved during the satellite era covered by MERRA.

da Silva, Arlindo↗

A Mixed Finite Volume Element Method for Flow Calculations in Porous Media

A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Jim E Jones↗

Global Ozone Loss Following Extreme Solar Proton Storms Based on the July 2012 Coronal Mass Ejection

Large solar coronal mass ejections pose a threat in the near-Earth space. As a cause of extreme periods of space weather, they can damage satellite-based communications and create geomagnetically induced currents in power and energy grids. Further, the solar wind energetic particles can reduce the protecting layer of atmospheric ozone and pose a threat to life on Earth. The large coronal mass ejection (CME) of July 2012, although directed away from the Earth, is often highlighted as a prime example of a potentially devastating super storm. Here we show, based on proton fluxes recorded by the instruments aboard the STEREO-A satellite, that the atmospheric response to the July 2012 event would have been comparable to those of the largest solar proton events of the satellite era. Significant impact on total ozone outside polar regions would require a much larger event, similar to those recorded in historical proxy data sets. Such an extreme event would cause long-term ozone reduction all the way to the equator and increase the size, duration, and depth of the Antarctic ozone hole. The impact would be comparable to predicted drastic and sudden ozone reduction from major volcanic eruptions, regional nuclear conflicts, or long-term stratospheric geoengineering.

ozone depletion↗

Space Shuttle filament wound case compressive strength study. II - Analysis

In order to explain the results of the testing conducted in Part I and relate the subscale test specimen performance to case performance during prelaunch loading, a comprehensive analysis program was conducted. Finite element models of both full-scale case and subscale test coupons were made. Validated axisymmetric models were used for the case analyses. Modeling of coupons was accomplished using a generalized plane strain finite element model in most cases, although some three-dimensional finite element modeling was used to investigate the importance of three-dimensional effects on the coupon performance. The large number of plies and complexity of the test coupons used limited the finite element grids to a through-the-thickness discretization of one element per ply. Delaminations observed in the helical ply dropoff region were modeled using frictionless sliders. Comparisons of predicted strains and displacements were made with the test data from both cases and coupons to validate the finite element models. Based on observations of the test results and finite element predictions, a failure criterion for both coupons and cases was postulated and then validated by comparing performance predictions to test results. Finally, the validated failure criterion was used to predict the performance of the first two flight cases under prelaunch loads.

Messick, M. J.↗

Operator induced multigrid algorithms using semirefinement

A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two and three dimensional model problems are presented, together with a two level analysis explaining these results.

Decker, Naomi↗

Operator induced multigrid algorithms using semirefinement

A variant of multigrid, based on zebra relaxation, and a new family of restriction/prolongation operators is described. Using zebra relaxation in combination with an operator-induced prolongation leads to fast convergence, since the coarse grid can correct all error components. The resulting algorithms are not only fast, but are also robust, in the sense that the convergence rate is insensitive to the mesh aspect ratio. This is true even though line relaxation is performed in only one direction. Multigrid becomes a direct method if an operator-induced prolongation is used, together with the induced coarse grid operators. Unfortunately, this approach leads to stencils which double in size on each coarser grid. The use of an implicit three point restriction can be used to factor these large stencils, in order to retain the usual five or nine point stencils, while still achieving fast convergence. This algorithm achieves a V-cycle convergence rate of 0.03 on Poisson's equation, using 1.5 zebra sweeps per level, while the convergence rate improves to 0.003 if optimal nine point stencils are used. Numerical results for two- and three-dimensional model problems are presented, together with a two level analysis explaining these results.

Decker, Naomi Henderson↗

Aerodynamic shape optimization using sensitivity analysis on third-order Euler equations

Two major advancements of an aerodynamic optimization method with two design variables are described. The first is represnted by an improved flow prediction method that uses the third-order accurate discretization of the Euler equations. This method makes it possible to predict the flowfield of a perturbed shape which generates shocks and other large gradients without intermediate CFD analysis. The second advancement is the use of every surface grid point as a design variable in the aerodynamic shape optimization problem. The improved algorithm is demonstrated by optimizing the ramp shape of a scramjet-afterbody configuration for maximum axial thrust. It is concluded that the improved flow-field prediction method eliminates most of the flow analysis and a priori guessing of all possible shapes from which the optimum is to be selected.

Baysal, Oktay↗

A Numerical Method for Direct Simulation of Turbulence in Complex Geometries

The ultimate goal of this work is to study the flow inside a channel with riblets on one of the two walls. The method has been tested for two dimensional flows in the presence of bodies with a geometrical singularity and for three dimensional flows inside domains described by Cartesian coordinates. The results have been compared with previous numerical simulations and with experimental results. The cases considered are: (1) the growth of Orr-Sommerfield waves in plane Poiseulle flow; (2) the flow over a backward facing step; (3) the flow past a wedge; and (4) the flow inside a narrow channel. Finally, the case of a channel with two large riblets on a wall has been simulated. In this case a limited number of grid points is sufficient, and in spite of the slow convergence for the pressure solver, one is able to obtain solutions with a reasonable amount of computer time. At present solutions with very fine grids in all three directions can not be obtained, due to lack of a fast pressure solver for general curvilinear coordinates.

P Orlandi↗

Surface photometry of disk galaxies at 60 and 100 microns - Radial gradients in dust temperature in optically thin disks

Surface, aperture, and grid photometry are performed at 60 and 100 microns for a sample of 28 nearby, large angular size galaxies. The spatial resolution of the data is 0.5-2.0 kpc, which allows for a coarse investigation of predicted radial gradients in physical dust temperature. A model is derived which predicts a drop of 3 K for every optical scale length, and this model agrees well with the data. In addition, for a subsample of 16 galaxies the radial profile at 60 microns is directly compared with that in the B or I band. The model predicts that optically thick disks should have scale lengths a factor of two smaller at 60 microns compared to B. In the optically thin case, the scale lengths are about equal except for the effect of the gradient in dust temperature. Virtually all of the data support the prediction of the optically thin case. Finally, 10 galaxies in the sample have sufficient angular sizes to allow for a coarse 2D spatial mapping of the 60/100 microns flux ratio. These maps clearly reveal that, when present, small amounts of hot dust dominate the emission at 60 microns and effectively mask the existence of a much larger amount of cold dust.

Bothun, Gregory D.↗

Reduced Navier Stokes Relaxation Procedures for Internal Flows

In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.

Rubin, Stanley G.↗

Latency Hiding in Dynamic Partitioning and Load Balancing of Grid Computing Applications

The Information Power Grid (IPG) concept developed by NASA is aimed to provide a metacomputing platform for large-scale distributed computations, by hiding the intricacies of highly heterogeneous environment and yet maintaining adequate security. In this paper, we propose a latency-tolerant partitioning scheme that dynamically balances processor workloads on the.IPG, and minimizes data movement and runtime communication. By simulating an unsteady adaptive mesh application on a wide area network, we study the performance of our load balancer under the Globus environment. The number of IPG nodes, the number of processors per node, and the interconnected speeds are parameterized to derive conditions under which the IPG would be suitable for parallel distributed processing of such applications. Experimental results demonstrate that effective solution are achieved when the IPG nodes are connected by a high-speed asynchronous interconnection network.

Das, Sajal K.↗

Preservation Environments

The long-term preservation of digital entities requires mechanisms to manage the authenticity of massive data collections that are written to archival storage systems. Preservation environments impose authenticity constraints and manage the evolution of the storage system technology by building infrastructure independent solutions. This seeming paradox, the need for large archives, while avoiding dependence upon vendor specific solutions, is resolved through use of data grid technology. Data grids provide the storage repository abstractions that make it possible to migrate collections between vendor specific products, while ensuring the authenticity of the archived data. Data grids provide the software infrastructure that interfaces vendor-specific storage archives to preservation environments.

Moore, Reagan W.↗

CAA for Jet Noise Physics: Issues and Recent Progress

Dr. Mankbadi summarized recent CAA results. Examples of the effect of various boundary condition schemes on the computed acoustic field, for a point source in a uniform flow, were shown. Solutions showing the impact of inflow excitations on the result were also shown. Results from a large eddy simulation, using a fourth-order MacCormack scheme with a Smagorinsky sub-grid turbulence model, were shown for a Mach 2.1 unheated jet. The results showed that the results were free from spurious modes. Results were shown for a Mach 1.4 jet using LES in the near field and the Kirchhoff method for the far field. Predicted flow field characteristics were shown to be in good agreement with data and predicted far field directivities were shown to be in qualitative agree with experimental measurements.

Mankbadi, Reda↗

CAA for Jet Noise Physics

Dr. Mankbadi summarized recent CAA results. Examples of the effect of various boundary condition schemes on the computed acoustic field, for a point source in a uniform flow, were shown. Solutions showing the impact of inflow excitations on the result were also shown. Results from a large eddy simulation, using a fourth-order MacCormack scheme with a Smagorinsky sub-grid turbulence model, were shown for a Mach 2.1 unheated jet. The results showed that the results were free from spurious modes. Results were shown for a Mach 1.4 jet using LES in the near field and the Kirchhoff method for the far field. Predicted flow field characteristics were shown to be in good agreement with data and predicted far field directivities were shown to be in qualitative agree with experimental measurements.

Mankbadi, Reda↗

Performance Analysis of a Hybrid Overset Multi-Block Application on Multiple Architectures

This paper presents a detailed performance analysis of a multi-block overset grid compu- tational fluid dynamics app!ication on multiple state-of-the-art computer architectures. The application is implemented using a hybrid MPI+OpenMP programming paradigm that exploits both coarse and fine-grain parallelism; the former via MPI message passing and the latter via OpenMP directives. The hybrid model also extends the applicability of multi-block programs to large clusters of SNIP nodes by overcoming the restriction that the number of processors be less than the number of grid blocks. A key kernel of the application, namely the LU-SGS linear solver, had to be modified to enhance the performance of the hybrid approach on the target machines. Investigations were conducted on cacheless Cray SX6 vector processors, cache-based IBM Power3 and Power4 architectures, and single system image SGI Origin3000 platforms. Overall results for complex vortex dynamics simulations demonstrate that the SX6 achieves the highest performance and outperforms the RISC-based architectures; however, the best scaling performance was achieved on the Power3.

Djomehri, M. Jahed↗

Lake Hoare, Antarctica: sedimentation through a thick perennial ice cover

Lake Hoare in the Dry Valleys of Antarctica is covered with a perennial ice cover more than 3 m thick, yet there is a complex record of sedimentation and of growth of microbial mats on the lake bottom. Rough topography on the ice covering the lake surface traps sand that is transported by the wind. In late summer, vertical conduits form by melting and fracturing, making the ice permeable to both liquid water and gases. Cross-sections of the ice cover show that sand is able to penetrate into and apparently through it by descending through these conduits. This is the primary sedimentation mechanism in the lake. Sediment traps retrieved from the lake bottom indicate that rates of deposition can vary by large amounts over lateral scales as small as 1 m. This conclusion is supported by cores taken in a 3 x 3 grid with a spacing of 1.5 m. Despite the close spacing of the cores, the poor stratigraphic correlation that is observed indicates substantial lateral variability in sedimentation rate. Apparently, sand descends into the lake from discrete, highly localized sources in the ice that may in some cases deposit a large amount of sand into the lake in a very short time. In some locations on the lake bottom, distinctive sand mounds have been formed by this process. They are primary sedimentary structures and appear unique to the perennially ice-covered lacustrine environment. In some locations they are tens of centimetres high and gently rounded with stable slopes; in others they reach approximately 1 m in height and have a conical shape with slopes at angle of repose. A simple formation model suggests that these differences can be explained by local variations in water depth and sedimentation rate. Rapid colonization of fresh sand surfaces by microbial mats composed of cyanobacteria, eukaryotic algae, and heterotrophic bacteria produces a complex intercalation of organic and sandy layers that are a distinctive form of modern stromatolites.

NASA Discipline Exobiology↗