Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regular 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗
NASA NTRS2002

Sudden ionospheric delay decorrelation and its impact on WAAS

In the absence of selective availability, the ionosphere represents the largest source of positioning error for single-frequency users of the GPS. In differential GPS systems such as the Wide Angle Augmentation System, vertical ionospheric delays are modeled at regularly-spaced intervals in latitude and longitude. The broadcast bound on the error at each of these points is designated the grid ionospheric vertical erorr (GIVE). A critical integrity requirement of WAAS is that the broadcast GIVE bounds residual error with a very high degree of confidence.

ionosphere WAAS GPS↗
NASA NTRS2000

Use of a GCM to Explore Sampling Issues in Connection with Satellite Remote Sensing of the Earth Radiation Budget

Collocated in time and space, top-of-the-atmosphere measurements of the Earth radiation budget (ERB) and cloudiness from passive scanning radiometers, and lidar- and radar-in-space measurements of multilayered cloud systems, are the required combination to improve our understanding of the role of clouds and radiation in climate. Experiments to fly multiple satellites "in formation" to measure simultaneously the radiative and optical properties of overlapping cloud systems are being designed. Because satellites carrying ERB experiments and satellites carrying lidars- or radars-in space have different orbital characteristics, the number of simultaneous measurements of radiation and clouds is reduced relative to the number of measurements made by each satellite independently. Monthly averaged coincident observations of radiation and cloudiness are biased when compared against more frequently sampled observations due, in particular, to the undersampling of their diurnal cycle, Using the Colorado State University General Circulation Model (CSU GCM), the goal of this study is to measure the impact of using simultaneous observations from the Earth Observing System (EOS) platform and companion satellites flying lidars or radars on monthly averaged diagnostics of longwave radiation, cloudiness, and its cloud optical properties. To do so, the hourly varying geographical distributions of coincident locations between the afternoon EOS (EOS-PM) orbit and the orbit of the ICESAT satellite set to fly at the altitude of 600 km, and between the EOS PM orbit and the orbits of the PICASSO satellite proposed to fly at the altitudes of 485 km (PICA485) or 705 km (PICA705), are simulated in the CSU GCM for a 60-month time period starting at the idealistic July 1, 2001, launch date. Monthly averaged diagnostics of the top-of-the-atmosphere, atmospheric, and surface longwave radiation budgets and clouds accumulated over grid boxes corresponding to satellite overpasses are compared against monthly averaged diagnostics obtained from hourly samplings over the entire globe. Results show that differences between irregularly (satellite) and regularly (true) sampled diagnostics of the longwave net radiative budgets are the greatest at the surface and the smallest in the atmosphere and at the top-of-the-atmosphere, under both cloud-free and cloudy conditions. In contrast, differences between the satellite and the true diagnostics of the longwave cloud radiative forcings are the largest in the atmosphere and at the top-of-the-atmosphere, and the smallest at the surface. A poorer diurnal sampling of the surface temperature in the satellite simulations relative to the true simulation contributes a major part to sampling biases in the longwave net radiative budgets, while a poorer diurnal sampling of cloudiness and its optical properties directly affects diagnostics of the longwave cloud radiative forcings. A factor of 8 difference in the number of satellite overpasses between PICA705 and PICA485 and ICESAT leads to a systematic factor of 3 difference in the spatial standard deviations of all radiative and cloudiness diagnostics.

Fowler, Laura D.↗
NASA NTRS2010

Improving Global Mass Flux Solutions from Gravity Recovery and Climate Experiment (GRACE) Through Forward Modeling and Continuous Time Correlation

We describe Earth's mass flux from April 2003 through November 2008 by deriving a time series of mas cons on a global 2deg x 2deg equal-area grid at 10 day intervals. We estimate the mass flux directly from K band range rate (KBRR) data provided by the Gravity Recovery and Climate Experiment (GRACE) mission. Using regularized least squares, we take into account the underlying process dynamics through continuous space and time-correlated constraints. In addition, we place the mascon approach in the context of other filtering techniques, showing its equivalence to anisotropic, nonsymmetric filtering, least squares collocation, and Kalman smoothing. We produce mascon time series from KBRR data that have and have not been corrected (forward modeled) for hydrological processes and fmd that the former produce superior results in oceanic areas by minimizing signal leakage from strong sources on land. By exploiting the structure of the spatiotemporal constraints, we are able to use a much more efficient (in storage and computation) inversion algorithm based upon the conjugate gradient method. This allows us to apply continuous rather than piecewise continuous time-correlated constraints, which we show via global maps and comparisons with ocean-bottom pressure gauges, to produce time series with reduced random variance and full systematic signal. Finally, we present a preferred global model, a hybrid whose oceanic portions are derived using forward modeling of hydrology but whose land portions are not, and thus represent a pure GRACE-derived signal.

Sabaka, T. J.↗
NASA NTRS2019

Daily evaluation of 26 precipitation datasets using Stage-IV gauge-radar data for the CONUS

New precipitation (P) datasets are released regularly, following innovations in weather forecasting models, satellite retrieval methods, and multi-source merging techniques. Using the conterminous US as a case study, we evaluated the performance of 26 gridded (sub-)daily P datasets to obtain insight into the merit of these innovations. The evaluation was performed at a daily timescale for the period 2008–2017 using the Kling–Gupta efficiency (KGE), a performance metric combining correlation, bias, and variability. As a reference, we used the high-resolution (4 km) Stage-IV gauge-radar P dataset. Among the three KGE components, the P datasets performed worst overall in terms of correlation (related to event identification). In terms of improving KGE scores for these datasets, improved P totals (affecting the bias score) and improved distribution of P intensity (affecting the variability score) are of secondary importance. Among the 11 gauge-corrected P datasets, the best overall performance was obtained by MSWEP V2.2, underscoring the importance of applying daily gauge corrections and accounting for gauge reporting times. Several uncorrected P datasets outperformed gauge-corrected ones. Among the 15 uncorrected P datasets, the best performance was obtained by the ERA5-HRES fourth-generation reanalysis, reflecting the significant advances in earth system modeling during the last decade. The (re)analyses generally performed better in winter than in summer, while the opposite was the case for the satellite-based datasets. IMERGHH V05 performed substantially better than TMPA-3B42RT V7, attributable to the many improvements implemented in the IMERG satellite P retrieval algorithm. IMERGHH V05 outperformed ERA5-HRES in regions dominated by convective storms, while the opposite was observed in regions of complex terrain. The ERA5-EDA ensemble average exhibited higher correlations than the ERA5-HRES deterministic run, highlighting the value of ensemble modeling. The WRF regional convection-permitting climate model showed considerably more accurate P totals over the mountainous west and performed best among the uncorrected datasets in terms of variability, suggesting there is merit in using high-resolution models to obtain climatological P statistics. Our findings provide some guidance to choose the most suitable P dataset for a particular application.

Hylke E. Beck↗
NASA NTRS1978

Finite element concepts in computational aerodynamics

Finite element theory was employed to establish an implicit numerical solution algorithm for the time averaged unsteady Navier-Stokes equations. Both the multidimensional and a time-split form of the algorithm were considered, the latter of particular interest for problem specification on a regular mesh. A Newton matrix iteration procedure is outlined for solving the resultant nonlinear algebraic equation systems. Multidimensional discretization procedures are discussed with emphasis on automated generation of specific nonuniform solution grids and accounting of curved surfaces. The time-split algorithm was evaluated with regards to accuracy and convergence properties for hyperbolic equations on rectangular coordinates. An overall assessment of the viability of the finite element concept for computational aerodynamics is made.

Baker, A. J.↗
NASA NTRS1985

Three-dimensional algorithms for grid restructuring in Free-Lagrangian calculations

Grid restructuring algorithms which lower the price of three-dimensional Free-Lagrange calculations are presented. The algorithms are first given for the case of planar triangulated surfaces embedded in and spanning a three-dimensional region. The tetrahedra generated by this technique form a Delaunay mesh if the interplane spacing is comparable to the resolution within the planes. The algorithm can therefore be used for efficient determinations of Voronoi connections for initial grids. Modifications of the algorithm for the case of closely spaced surfaces are demonstrated in the context of restructuring algorithms which can accommodate colliding surfaces. Then, the restriction to planar surfaces is removed and regular surfaces are examined. The basic algorithm is the same, with an additional operation to project the vertices of one surface onto another. Finally, vertices on the surface are allowed to migrate anywhere in space.

Fritts, M.↗
NASA NTRS1975

A grid of model atmospheres for metal-deficient giant stars. I

A grid of flux-constant model atmospheres for stars with effective temperatures between 3750 and 6000 K, log g between 0.75 and 3.0, and (A/H) between -3.0 and 0.0 has been constructed. The line absorption is approximated by opacity distribution functions. Metal lines and molecular lines, including those from the infrared bands of CO and CN, are taken into account. The variation of the structure of the models with metal abundance and microturbulence parameter is found to be quite regular and not very drastic. The surface cooling produced by CO is important for all the models with a maximum temperature of 5000 K, while CN mainly causes a backwarming effect but is not very important for solar CNO abundances. The effects of convection, estimated by using the mixing-length approximation, are important only for the coldest models. The models compare very well with models from overlapping regions of other grids. A solar model consistent with the grid models is found to agree satisfactorily with empirical solar model atmospheres. The careful use of scaled solar models for stars with (A/H) approximately equal to zero is justified.

Gustafsson, B.↗
NASA NTRS2020

The Evolution of a Circumplanetary Disc with a Dead Zone

We investigate whether the regular Galilean satellites could have formed in the dead zone of a circumplanetary disc. A deadzone is a region of weak turbulence in which the magnetorotational instability is suppressed, potentially an ideal environmentfor satellite formation. With the grid-based hydrodynamic codeFAR GO3D, we examine the evolution of a circumplanetary discmodel with a dead zone. Material accumulates in the dead zone of the disc leading to a higher total mass and but a similartemperature profile compared to a fully turbulent disc model. The tidal torque increases the rate of mass transport through thedead zone leading to a steady-state disc with a dead zone that does not undergo accretion outbursts. We explore a range of disc,dead zone, and mass inflow parameters and find that the maximum mass of the disc is around 0.001MJ. Since the total solidmass of such a disc is much lower, we find that there is not sufficient material in the disc forin situformation of the Galileansatellites and that external supplement is required.

Cheng Chen↗
NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗
NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2017

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗
NASA NTRS2021

Evaluation of Global Surface Water Temperature Data Sets for use in Passive Remote Sensing of Soil Moisture

Inland open water bodies often pose a systematic error source in the passive remote sensing retrievals of soil moisture. Water temperature is a necessary variable used to compute water emissions that is required to be subtracted from satellite observation to yield actual emissions from the land portion, which in turn generates accurate soil moisture retrievals. Therefore, overestimation of soil moisture can often be corrected using concurrent water temperature data in the overall mitigation procedure. In recent years, several data sets of lake water temperature have become available, but their specifications and accuracy have rarely been investigated in the context of passive soil moisture remote sensing on a global scale. For this reason, three lake temperature products were evaluated against in-situ measurements from 2007 to 2011. The data sets include the lake surface water temperature (LSWT) from Global Observatory of Lake Responses to Environmental Change (GloboLakes), the Copernicus Global Land Operations Cryosphere and Water (C-GLOPS), as well as the lake mix-layer temperature (LMLT) from the European Centers for Medium-Range Weather Forecast (ECMWF) ERA5 Land Reanalysis. GloboLakes, C-GLOPS, and ERA5 Land have overall comparable performance with Pearson correlations (R) of 0.87, 0.92 and 0.88 in comparison with in-situ measurements. LSWT products exhibit negative median biases of −0.27 K (GloboLakes) and −0.31 K (C-GLOPS), whereas the median bias of LMLT is 1.56 K. When mapped from their respective native resolutions to a common 9 km Equal-Area Scalable Earth (EASE) Grid 2.0 projection, similar relative performance was observed. LMLT and LSWT data are closer in performance over the 9 km grid cells that exhibit a small range of lake cover fractions (0.05–0.5). Despite comparable relative performance, ERA5 Land shows great advantages in spatial coverage and temporal resolution. In summary, an integrated evaluation on data accuracy, long-term availability, global coverage, temporal resolution, and regular forward processing with modest data latency led us to conclude that LMLT from the ERA5 Land Reanalysis product represents the most optimal path for use in the development of a long-term soil moisture product.

Runze Zhang↗
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS1995

Quantitative characterization of the small-scale fracture patterns on the plains of Venus

The objectives of this research project were to (1) compile a comprehensive database of the occurrence of regularly spaced kilometer scale lineations on the volcanic plains of Venus in an effort to verify the effectiveness of the shear-lag model developed by Banerdt and Sammis (1992), and (2) develop a model for the formation of irregular kilometer scale lineations such as typified in the gridded plains region of Guinevere Planitia. Attached to this report is the paper 'A Tectonic Model for the Formation of the Gridded Plains on Guinevere Planitia, Venus, and Implications for the Elastic Thickness of the Lithosphere'.

Sammis, Charles G.↗
NASA NTRS1990

Highly parallel sparse Cholesky factorization

Several fine grained parallel algorithms were developed and compared to compute the Cholesky factorization of a sparse matrix. The experimental implementations are on the Connection Machine, a distributed memory SIMD machine whose programming model conceptually supplies one processor per data element. In contrast to special purpose algorithms in which the matrix structure conforms to the connection structure of the machine, the focus is on matrices with arbitrary sparsity structure. The most promising algorithm is one whose inner loop performs several dense factorizations simultaneously on a 2-D grid of processors. Virtually any massively parallel dense factorization algorithm can be used as the key subroutine. The sparse code attains execution rates comparable to those of the dense subroutine. Although at present architectural limitations prevent the dense factorization from realizing its potential efficiency, it is concluded that a regular data parallel architecture can be used efficiently to solve arbitrarily structured sparse problems. A performance model is also presented and it is used to analyze the algorithms.

Gilbert, John R.↗
NASA NTRS1980

Effect of grid system on finite element calculation

Detailed parametric studies of the effect of grid system on finite element calculation for potential flows were made. These studies led to the formulation of a design criteria for optimum mesh system and the development of two methods to generate the optimum mesh system. The guidelines for optimum mesh system are: (1) the mesh structure should be regular; (2) the element should be as regular and equilateral as possible; (3) the distribution of size of element should be consistent with that of flow variables to insure maximum uniformity in error distribution; (4) for non-Dirichlet boundary conditions, smaller boundary elements or higher order interpolation functions should be used; and (5) the mesh should accommodate the boundary geometry as accurately as possible. The results of the parametric studies are presented.

Lee, K. D.↗
NASA NTRS1990

Fresh look at floating shock fitting

A fast implicit upwind procedure for the two-dimensional Euler equations is described that allows accurate computations of shocked flows on nonadapted meshes. Away from shocks, the second-order accurate upwinding is based on the split-coefficient-matrix (SCM) method. In the presence of shocks, the difference stencils are modified using a floating shock fitting technique. Rapid convergence to steady-state solutions is attained with a diagonalized approximate factorization (AF) algorithm. Results are presented for Riemann's problem, for a regular shock reflection at an inviscid wall, for supersonic flow past a cylinder, and for a transonic airfoil. All computed shocks are ideally sharp and in excellent agreement with other numerical results or 'exact' solutions. Most importantly, this has been accomplished on unusually crude meshes without any attempt to align grid lines with shock fronts or to cluster grid lines around shocks.

Hartwich, PETER-M.↗