Search NASA⌕ Search

SEARCH · Search NASA

Results for “Overlapping grid”

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 163 records · Page 9

Plane Smoothers for Multiblock Grids: Computational Aspects

Standard multigrid methods are not well suited for problems with anisotropic discrete operators, which can occur, for example, on grids that are stretched in order to resolve a boundary layer. One of the most efficient approaches to yield robust methods is the combination of standard coarsening with alternating-direction plane relaxation in the three dimensions. However, this approach may be difficult to implement in codes with multiblock structured grids because there may be no natural definition of global lines or planes. This inherent obstacle limits the range of an implicit smoother to only the portion of the computational domain in the current block. This report studies in detail, both numerically and analytically, the behavior of blockwise plane smoothers in order to provide guidance to engineers who use block-structured grids. The results obtained so far show alternating-direction plane smoothers to be very robust, even on multiblock grids. In common computational fluid dynamics multiblock simulations, where the number of subdomains crossed by the line of a strong anisotropy is low (up to four), textbook multigrid convergence rates can be obtained with a small overlap of cells between neighboring blocks.

Llorente, Ignacio M.↗

Assessment of Numerical and Modeling Errors of RANS based Transition Models for Low-Reynolds Numbers 2-D Flows

In this paper we report the outcome of selected workshops organized as part of the NATO Applied Vehicle Technology (AVT)-313 activity Incompressible Laminar-to-Turbulent Flow Transition Study that focused on assessing the numerical and modeling accuracy of the γ−Reθ and γ transition models coupled to the k−ω Shear-Stress Transport (SST) two-equation eddy-viscosity model. Three different test cases involving nominally 2D flow configurations were selected: flow over a flat plate with two different levels of turbulence intensity at the inlet; flow around the Eppler 387 foil at a Reynolds number of 3×10^5 and angles of attack of 1 deg. and 7 deg. flow around the NACA 0015 foil at a Reynolds number of 1.8×10^5 and angles of attack of 5 deg. and10 deg. The flat plate flow conditions correspond to natural and by-pass transition, whereas the other two test cases include laminar separation bubbles that lead to separation-induced transition. For each test case, the selected quantities of interest include both integral and local flow quantities. Geometrically similar grids with a wide range of grid refinement ratios were generated for each of the test cases to allow the estimation of numerical uncertainties for all quantities of interest selected for this study. Several RANS flow solvers were used, employing common grids with the same boundary conditions and mathematical models. Therefore, it is possible to analyze the consistency of the results, i.e., to check if the intervals defined by the different numerical solutions with their respective uncertainties overlap with each other. Modeling errors can also be addressed for the selected flow quantities that have experimental data available. However, the experimental information available in these cases is not sufficient to guarantee that experiments and simulations are performed with the same settings. Nonetheless, the available experimental data is sufficient to guarantee that modeling errors are significantly reduced with the use of the transition models when compared to simulations performed using only the k−ω SST model.

CFD Modeling↗

A Novel Method for Estimating Shortwave Direct Radiative Effect of Above-Cloud Aerosols Using CALIOP and MODIS Data

This paper describes an efficient and unique method for computing the shortwave direct radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds using CALIOP and MODIS data. It addresses the overlap of aerosol and cloud rigorously by utilizing the joint histogram of cloud optical depth and cloud top pressure while also accounting for subgrid-scale variations of aerosols. The method is computationally efficient because of its use of grid-level cloud and aerosol statistics, instead of pixel-level products, and a pre-computed look-up table based on radiative transfer calculations. We verify that for smoke over the southeast Atlantic Ocean the method yields a seasonal mean instantaneous (approximately 1:30PM local time) shortwave DRE of above cloud aerosol (ACA) that generally agrees with more rigorous pixel-level computation within 4 percent. We also estimate the impact of potential CALIOP aerosol optical depth (AOD) retrieval bias of ACA on DRE. We find that the regional and seasonal mean instantaneous DRE of ACA over southeast Atlantic Ocean would increase, from the original value of 6.4 W m(-2) based on operational CALIOP AOD to 9.6 W m(-2) if CALIOP AOD retrieval are biased low by a factor of 1.5 (Meyer et al., 2013) and further to 30.9 W m(-2) if CALIOP AOD retrieval are biased low by a factor of 5 as suggested in (Jethva et al., 2014). In contrast, the instantaneous ACA radiative forcing efficiency (RFE) remains relatively invariant in all cases at about 53 W m(-2) AOD(-1), suggesting a near linear relation between the instantaneous RFE and AOD. We also compute the annual mean instantaneous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) over global oceans based on 4 years of CALIOP and MODIS data. We find that the variability of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is mainly driven by the optical depth of the underlying clouds. While we demonstrate our method using CALIOP and MODIS data, it can also be extended to other satellite data sets, as well as climate model outputs.

CALIOP and MODIS↗

A De-Centralized Scheduling and Load Balancing Algorithm for Heterogeneous Grid Environments

In the past two decades, numerous scheduling and load balancing techniques have been proposed for locally distributed multiprocessor systems. However, they all suffer from significant deficiencies when extended to a Grid environment: some use a centralized approach that renders the algorithm unscalable, while others assume the overhead involved in searching for appropriate resources to be negligible. Furthermore, classical scheduling algorithms do not consider a Grid node to be N-resource rich and merely work towards maximizing the utilization of one of the resources. In this paper we propose a new scheduling and load balancing algorithm for a generalized Grid model of N-resource nodes that not only takes into account the node and network heterogeneity, but also considers the overhead involved in coordinating among the nodes. Our algorithm is de-centralized, scalable, and overlaps the node coordination time of the actual processing of ready jobs, thus saving valuable clock cycles needed for making decisions. The proposed algorithm is studied by conducting simulations using the Message Passing Interface (MPI) paradigm.

Arora, Manish↗

ISCCP Cloud Properties Associated with Standard Cloud Types Identified in Individual Surface Observations

Individual surface weather observations from land stations and ships are compared with individual cloud retrievals of the International Satellite Cloud Climatology Project (ISCCP), Stage C1, for an 8-year period (1983-1991) to relate cloud optical thicknesses and cloud-top pressures obtained from satellite data to the standard cloud types reported in visual observations from the surface. Each surface report is matched to the corresponding ISCCP-C1 report for the time of observation for the 280x280-km grid-box containing that observation. Classes of the surface reports are identified in which a particular cloud type was reported present, either alone or in combination with other clouds. For each class, cloud amounts from both surface and C1 data, base heights from surface data, and the frequency-distributions of cloud-top pressure (p(sub c) and optical thickness (tau) from C1 data are averaged over 15-degree latitude zones, for land and ocean separately, for 3-month seasons. The frequency distribution of p(sub c) and tau is plotted for each of the surface-defined cloud types occurring both alone and with other clouds. The average cloud-top pressures within a grid-box do not always correspond well with values expected for a reported cloud type, particularly for the higher clouds Ci, Ac, and Cb. In many cases this is because the satellites also detect clouds within the grid-box that are outside the field of view of the surface observer. The highest average cloud tops are found for the most extensive cloud type, Ns, averaging 7 km globally and reaching 9 km in the ITCZ. Ns also has the greatest average retrieved optical thickness, tau approximately equal 20. Cumulonimbus clouds may actually attain far greater heights and depths, but do not fill the grid-box. The tau-p(sub c) distributions show features that distinguish the high, middle, and low clouds reported by the surface observers. However, the distribution patterns for the individual low cloud types (Cu, Sc, St) occurring alone overlap to such an extent that it is not possible to distinguish these cloud types from each other on the basis of tau-p(sub c) values alone. Other cloud types whose tau-p(sub c) distributions are indistinguishable are Cb, Ns, and thick As. However, the tau-p(sub c) distribution patterns for the different low cloud types are nevertheless distinguishable when all occurrences of a low cloud type are included, indicating that the different low types differ in their probabilities of co-occurrence with middle and high clouds.

Hahn, Carole J.↗

Application of FUN3D and CFL3D to the Third Workshop on CFD Uncertainty Analysis

Two Reynolds-averaged Navier-Stokes computer codes - one unstructured and one structured - are applied to two workshop cases (for the 3rd Workshop on CFD Uncertainty Analysis, held at Instituto Superior Tecnico, Lisbon, in October 2008) for the purpose of uncertainty analysis. The Spalart-Allmaras turbulence model is employed. The first case uses the method of manufactured solution and is intended as a verification case. In other words, the CFD solution is expected to approach the exact solution as the grid is refined. The second case is a validation case (comparison against experiment), for which modeling errors inherent in the turbulence model and errors/uncertainty in the experiment may prevent close agreement. The results from the two computer codes are also compared. This exercise verifies that the codes are consistent both with the exact manufactured solution and with each other. In terms of order property, both codes behave as expected for the manufactured solution. For the backward facing step, CFD uncertainty on the finest grid is computed and is generally very low for both codes (whose results are nearly identical). Agreement with experiment is good at some locations for particular variables, but there are also many areas where the CFD and experimental uncertainties do not overlap.

Rumsey, C. L.↗

Parallel volume ray-casting for unstructured-grid data on distributed-memory architectures

As computing technology continues to advance, computational modeling of scientific and engineering problems produces data of increasing complexity: large in size and unstructured in shape. Volume visualization of such data is a challenging problem. This paper proposes a distributed parallel solution that makes ray-casting volume rendering of unstructured-grid data practical. Both the data and the rendering process are distributed among processors. At each processor, ray-casting of local data is performed independent of the other processors. The global image composing processes, which require inter-processor communication, are overlapped with the local ray-casting processes to achieve maximum parallel efficiency. This algorithm differs from previous ones in four ways: it is completely distributed, less view-dependent, reasonably scalable, and flexible. Without using dynamic load balancing, test results on the Intel Paragon using from two to 128 processors show, on average, about 60% parallel efficiency.

Ma, Kwan-Liu↗

Rectangular beam (5 X 40 cm multipole ion source)

A 5 x 40 cm rectangular-beam ion source was designed and fabricated. A multipole field configuration was used to facilitate design of the modular rectangular chamber, while a three-grid ion optics system was used for increased ion current densities. For the multipole chamber, a magnetic integral of 0.000056 Tesla-m was used to contain the primary electrons. This integral value was reduced from the initial design value, with the reduction found necessary for discharge stability. The final value of magnetic integral resulted in discharge losses at typical operating conditions which ranged from 600 to 1000 eV/ion, in good agreement with the design value of 800 eV/ion. The beam current density at the ion optics was limited to about 3.2 mA/sq cm at 500 eV and to about 3.5 mA/sq cm at 1000 ev. The effects of nonuniform ion current, dimension tolerance, and grid thermal warping were considered. The use of multiple rectangular-beam ion sources to process wider areas than would be possible with a single source (approx. 40 cm) was also studied. Beam profiles were surveyed at a variety of operating conditions and the results of various amounts of beam overlap calculated.

Haynes, C. M.↗

A zonal flow analysis method for two-dimensional airfoils

A closed-loop, overlapped, velocity-coupling procedure has been utilized to combine a two-dimensional potential flow panel code and a Navier-Stokes code. The fully coupled, two-zone code has been used to compute the flow past a NACA 0012 airfoil. For this case, the zonal method has shown that the grid domain size can be reduced to 0.14 chord lengths with less than 1 percent loss in accuracy. Further, the required computation time is reduced by a factor of approximately four.

Summa, J. Michael↗

Advanced 30/20 GHz multiple beam antenna for future communications satellites

A multiple offset reflector system is being used to provide multiple fixed and scanning beams at 30 GHz by means of a polarization grid diplexer to separate two orthogonally polarized feed assemblies. Each feed assembly consists of multiple fixed and scanning beam feeds. The technology developments required for successful operation of this proposed system are presented. These technologies include: (1) an offset Cassegrain system for low sidelobe wide angle scan beams; (2) a diplexed beam forming network for overlapping low sidelobe beams; (3) a fast switching network for scanning beams; and (4) fabrication of precision feed components and large offset reflectors.

Chen, C. C.↗

The role of shocks in NGC 6240

The case of the infrared-bright system NGC 6240 encapsulates in many ways the major issues in understanding processes in infrared-bright systems in general. Its morphology strongly suggests that this is a merger in progress, with energetic phenomena being triggered as evidenced by strong radio emission and optical emission lines as well as the strong far-infrared output. Evidence has been produced supporting various schemes for producing the energy being radiated by this object. Early work on the optical spectrum suggested a Seyfert 2 classification, but the ionization level and spatial extent of the emission implied a more distributed source of energy. IR line studies found evidence for young supergiants, and for immense amounts of H2 excited by low-velocity shocks, which might contribute to a large enhancement in star-formation rate. Finally, the radio-continuum morphology shows compact sources possibly associated with active nuclei in both remnants of the merger. This is one of the nearest IR-bright galaxies objects at high luminosity (Cz = 7500 km s(exp -1), so that spatially resolved studies can yield further insight into which of these processes are most important - star formation, nuclear activity, or shock excitation. The velocity field, linewidth distribution, and ionization structure of the gas in NGC 6240 was mapped with the Dense-Pak fiber-optic array on the Kitt Peak 4-m telescope. In all, 135 spectra were obtained through an array of 2 inch apertures, with considerable overlap between many of the adjacent apertures produced by offsetting the telescope between the exposures. As the strongest features with lambda less than 1.1 ma, H alpha and the adjacent (NII) lines were observed, along with the neighboring (O I) and (S II) features. Maps of velocity and linewidth were constructed by interpolating measurements from these spectra onto a regular grid, and show irregular patterns with large velocity gradients and, particularly, an extended region of large local (shear?) velocity width.

Keel, William C.↗

Finite element computation on nearest neighbor connected machines

Research aimed at faster, more cost effective parallel machines and algorithms for improving designer productivity with finite element computations is discussed. A set of 8 boards, containing 4 nearest neighbor connected arrays of commercially available floating point chips and substantial memory, are inserted into a commercially available machine. One-tenth Mflop (64 bit operation) processors provide an 89% efficiency when solving the equations arising in a finite element problem for a single variable regular grid of size 40 by 40 by 40. This is approximately 15 to 20 times faster than a much more expensive machine such as a VAX 11/780 used in double precision. The efficiency falls off as faster or more processors are envisaged because communication times become dominant. A novel successive overrelaxation algorithm which uses cyclic reduction in order to permit data transfer and computation to overlap in time is proposed.

Mcaulay, A. D.↗

A novel method of calculating far-field patterns of large aperture antennas

A method is described for calculation of the radiation pattern of large aperture antennas. A piece-wise linear approximation of the aperture field using overlapping pyramidal basis functions allows the radiation pattern of an aperture antenna to be calculated as though it were a two-dimensional array. The calculation of radiation pattern data versus theta and phi, suitable for 3-D or contour plot algorithms, is achieved by locating the array in the yz-plane and performing a summation over the aperture field data sampled on a square grid. A FORTRAN subroutine is provided for performing radiation pattern calculations. Numerical results are included to demonstrate the accuracy and convergence of the method. These numerical results indicate that typical accuracies of + or - 0.1 dB for Directivity, + or - dB for the 1st Sidelobe Level, and + - 2dB for the 2nd Sidelobe Level can be obtained with an aperture grid of 45x45 points and requires approximately 0.02 seconds CPU time per far-field data point on a VAX 11/750 with a floating point accelerator.

Bailey, M. C.↗

The European Russia Drought Atlas (1400-2016 CE)

We present the European Russia Drought Atlas (ERDA) that covers the East European Plain to the Ural Mountains from 1400–2016 CE. Like the Old World Drought Atlas (OWDA) for the Euro-Mediterranean region, the ERDA is a one-half degree gridded reconstruction of summer Palmer Drought Severity Indices estimated from a network of annual tree-ring chronologies. Ensemble point-by-point regression is used to generate the ERDA with the identical protocols used for developing the OWDA. Split calibration/validation tests of the ERDA indicate that it has significant skill over most of its domain and is much more skillful than the OWDA where they overlap in the western part of ERDA domain. Comparisons to historical droughts over European Russia additionally support the ERDA’s overall validity. The ERDA has been spatially smoothed and infilled using a local regression method to yield a spatially complete drought atlas back to 1400 CE. EOF analysis indicates that there are three principal modes of hydroclimatic variability in the ERDA. After Varimax rotation, these modes correlate significantly with independent climate data sets extending back to the late nineteenth century in a physically interpretable way and relate to atmospheric circulation dynamics of droughts and heatwaves over European Russia based on more recent instrumental data.

Cook, Edward R.↗

A 5 x 40 cm rectangular-beam multipole ion source

A rectangular ion source particularly suited for the continuous sputter processing of materials over a wide area is discussed. A multipole magnetic field configuration was used to design an ion source with a 5 x 40 cm beam area, while a three-grid ion optics system was used to maximize ion current density at the design ion energy of 500 eV. An average extracted current density of about 4 mA/sq cm could be obtained from 500 eV Ar ions. The difference between the experimental performance and the design value of 6 mA/sq cm is attributed to grid misalignment due to thermal expansion. The discharge losses at typical operating conditions ranged from about 600 to 1000 eV/ion, in reasonable agreement with the design value of 800 eV/ion. The use of multiple rectangular-beam ion sources to process wider areas than would be possible with a single source was also studied, and the most uniform coverage was found to be obtainable with a 0 to 2 cm overlap.

Robinson, R. S.↗

Mach number validation of a new zonal CFD method (ZAP2D) for airfoil simulations

A closed-loop overlapped velocity coupling procedure has been utilized to combine a two-dimensional potential-flow panel code and a Navier-Stokes code. The fully coupled two-zone code (ZAP2D) has been used to compute the flow past a NACA 0012 airfoil at Mach numbers ranging from 0.3 to 0.84 near the two-dimensional airfoil C(lmax) point for a Reynolds number of 3 million. For these cases, the grid domain size can be reduced to 3 chord lengths with less than 3-percent loss in accuracy for freestream Mach numbers through 0.8. Earlier validation work with ZAP2D has demonstrated a reduction in the required Navier-Stokes computation time by a factor of 4 for subsonic Mach numbers. For this more challenging condition of high lift and Mach number, the saving in CPU time is reduced to a factor of 2.

Strash, Daniel J.↗

A Tool for Defining Tow-Steered Laminates for Finite Element Grids

Tailoring of composite laminates is traditionally performed by changing the orientation of straight fibers in one or more plies. Modern automated fiber placement machines facilitate placement of bundles of curved fibers (tows) in a process called tow-steering, but additional variables must be used to define the shapes of tow-steered fiber paths. In this paper, a Python-based tool, called Automated Tool for Steered COmposite Optimizable Laminates (ATSCOOL), is presented for modeling steered tows as chains of circular arcs. The geometry for defining tow path shapes using four input variables is described. The formulae for determining gaps and overlaps between clusters of consecutive tows, called courses, and for determining the number of steered courses needed to cover a rectangular panel are presented. Graphical representation of courses and thickness distribution in an example panel using the ATSCOOL software are presented. Finally, an example of a finite element analysis performed using the property information output from ATSCOOL is shown.

finite element analysis↗

The NASA/MSFC Global Reference Atmospheric Model: 1999 Version (GRAM-99)

The latest version of Global Reference Atmospheric Model (GRAM-99) is presented and discussed. GRAM-99 uses either (binary) Global Upper Air Climatic Atlas (GUACA) or (ASCII) Global Gridded Upper Air Statistics (GGUAS) CD-ROM data sets, for 0-27 km altitudes. As with earlier versions, GRAM-99 provides complete geographical and altitude coverage for each month of the year. GRAM-99 uses a specially-developed data set, based on Middle Atmosphere Program (MAP) data, for 20-120 km altitudes, and NASA's 1999 version Marshall Engineering Thermosphere (MET-99) model for heights above 90 km. Fairing techniques assure smooth transition in overlap height ranges (20-27 km and 90-120 km). GRAM-99 includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He and H). A variable-scale perturbation model provides both large-scale (wave) and small-scale (stochastic) deviations from mean values for thermodynamic variables and horizontal and vertical wind components. The small-scale perturbation model includes improvements in representing intermittency ("patchiness"). A major new feature is an option to substitute Range Reference Atmosphere (RRA) data for conventional GRAM climatology when a trajectory passes sufficiently near any RRA site. A complete user's guide for running the program, plus sample input and output, is provided. An example is provided for how to incorporate GRAM-99 as subroutines in other programs (e.g., trajectory codes).

Justus, C. G.↗