Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space 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 559 records · Page 31

File Specification for GEOS-5 FP (Forward Processing)

The GEOS-5 FP Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The assimilation is performed at a horizontal resolution of 0.3125-degree longitude by 0.25- degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational FP, using the most recent version of the GEOS-5 assimilation system. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary. Documentation about the current access methods for products described in this document can be found on the GMAO products page: http://gmao.gsfc.nasa.gov/products/.

GSI↗

File Specification for GEOS-5 FP-IT (Forward Processing for Instrument Teams)

The GEOS-5 FP-IT Atmospheric Data Assimilation System (GEOS-5 ADAS) uses an analysis developed jointly with NOAA's National Centers for Environmental Prediction (NCEP), which allows the Global Modeling and Assimilation Office (GMAO) to take advantage of the developments at NCEP and the Joint Center for Satellite Data Assimilation (JCSDA). The GEOS-5 AGCM uses the finite-volume dynamics (Lin, 2004) integrated with various physics packages (e.g, Bacmeister et al., 2006), under the Earth System Modeling Framework (ESMF) including the Catchment Land Surface Model (CLSM) (e.g., Koster et al., 2000). The GSI analysis is a three-dimensional variational (3DVar) analysis applied in grid-point space to facilitate the implementation of anisotropic, inhomogeneous covariances (e.g., Wu et al., 2002; Derber et al., 2003). The GSI implementation for GEOS-5 FP-IT incorporates a set of recursive filters that produce approximately Gaussian smoothing kernels and isotropic correlation functions. The GEOS-5 ADAS is documented in Rienecker et al. (2008). More recent updates to the model are presented in Molod et al. (2011). The GEOS-5 system actively assimilates roughly 2 × 10(exp 6) observations for each analysis, including about 7.5 × 10(exp 5) AIRS radiance data. The input stream is roughly twice this volume, but because of the large volume, the data are thinned commensurate with the analysis grid to reduce the computational burden. Data are also rejected from the analysis through quality control procedures designed to detect, for example, the presence of cloud. To minimize the spurious periodic perturbations of the analysis, GEOS-5 FP-IT uses the Incremental Analysis Update (IAU) technique developed by Bloom et al. (1996). More details of this procedure are given in Appendix A. The analysis is performed at a horizontal resolution of 0.625-degree longitude by 0.5-degree latitude and at 72 levels, extending to 0.01 hPa. All products are generated at the native resolution of the horizontal grid. The majority of data products are time-averaged, but four instantaneous products are also available. Hourly data intervals are used for two-dimensional products, while 3-hourly intervals are used for three-dimensional products. These may be on the model's native 72-layer vertical grid or at 42 pressure surfaces extending to 0.1 hPa. This document describes the gridded output files produced by the GMAO near real-time operational GEOS-5 FP-IT processing in support of the EOS instrument teams. Additional details about variables listed in this file specification can be found in a separate document, the GEOS-5 File Specification Variable Definition Glossary.

3D Var↗

Adaptive finite differencing in high accuracy electronic structure calculations

Abstract A multi-order Adaptive Finite Differencing (AFD) method is developed for the kinetic energy operator in real-space, grid-based electronic structure codes. It uses atomic pseudo orbitals produced by the corresponding pseudopotential codes to optimize the standard finite difference (SFD) operators for improved precision. Results are presented for a variety of test systems and Bravais lattice types, including the well-known Δ test for 71 elements in the periodic table, the Mott insulator NiO, and borax decahydrate, which contains covalent, ionic, and hydrogen bonds. The tests show that an 8th-order AFD operator leads to the same average Δ value as that achieved by plane-wave codes and is typically far more accurate and has a much lower computational cost than a 12th-order SFD operator. The scalability of real-space electronic calculations is demonstrated for a 2016-atom NiO cell, for which the computational time decreases nearly linearly when scaled from 18 to 144 CPU-GPU nodes.

Briggs, E. L. (ORCID:0000000343983492)↗

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗

Self-Driving Microscopy for AI/ML-Enabled Physics Discovery and Materials Optimization

Materials are the bedrock of economy and foundation for all real-world technologies. The viability of space travel, grid energy storage, solar to fuels conversion, methane removal, and photovoltaic energy solutions hinge on the discovery and optimization of novel materials and rapid scaling toward manufacturing. The last 20 years have seen an exponential growth in the theoretical predictive capability for crystalline materials and small molecules. However, it is only in the last five years that we have seen the rapid expansion of high-throughput synthesis enabled by laboratory robotics and microfluidics, as well as a resurgence of combinatorial synthesis (Abolhasani and Kumacheva 2023; Epps and Abolhasani 2021; Jiang et al. 2022; Rajan 2008; Soldatov et al. 2021; Szymanski et al. 2023). Combinatorial synthesis, microfluidics, and ultimately dip-pen megalibraries have demonstrated the ability to “write” multicomponent nanomaterials at high throughput scale, generating millions of material examples in the 3D, 4D, and 5D composition spaces (Chen et al. 2016, 2019; Jibril et al. 2022).

36 MATERIALS SCIENCE↗

Energy Beam Highways Through the Skies

The emergence of Energy Beam Flight Transportation Systems could dramatically change the way we travel in the 21st Century. A framework for formulating 'Highways of Light' and the top level architectures that invoke radically new Space Power Grid infrastructure, are introduced. Basically, such flight systems, hereafter called Lightcraft, would employ off-board energy beam sources (either laser or microwave) to energize on-board dependent 'motors' -- instead of the traditional autonomous 'engines' with their on-board energy sources (e.g., chemical fuels). Extreme reductions in vehicle dry mass appear feasible with the use of off-board power and a high degree of on-board artificial intelligence. Such vehicles may no longer need airports for refueling (since they require no propellant), and could possibly pick up travelers at their homes -- before motoring over to one of many local boost stations, for the flight out. With off-board power, hyper-energetic acceleration performance and boost-glide trajectories become feasible. Hypersonic MS airbreathing propulsion can enable boosts up to twice escape velocity, which will cut trip times to the moon down to 5.5 hours. The predominant technological, environmental and social factors that will result from such transportation systems will be stressed. This presentation first introduces the remote source siting options for the space power system infrastructure, and then provides three representative laser/microwave Lightcraft options (derived from historical Case Studies): i.e., 'Acorn', 'Toy Top', and 'Disc.' Next the gamut of combined-cycle engine options developed for these Lightcraft are examined -- to illuminate the 'emerging technologies' that must be harnessed to produce flight hardware. Needed proof-of concept experiments are identified, along with the Macro-Level Issues that can springboard these revolutionary concepts into hardware reality.

Myrabo, Leik N.↗

A three-dimensional space marching algorithm for the solution of the Euler equations on unstructured grids

A three-dimensional space marching algorithm using an unstructured discretization is proposed. The governing equations and the spatial discretization of the domain and the interpolation from one arbitrary unstructured mesh to another are discussed, and the grid generation technique is described. Solutions for several different geometries are presented, and the results are compared with existing numerical techniques and experimental results.

Mcgrory, W. D.↗

Numerical computation of complex multi-body Navier-Stokes flows with applications for the integrated Space Shuttle launch vehicle

An enhanced grid system for the Space Shuttle Orbiter was built by integrating CAD definitions from several sources and then generating the surface and volume grids. The new grid system contains geometric components not modeled previously plus significant enhancements on geometry that has been modeled in the old grid system. The new orbiter grids were then integrated with new grids for the rest of the launch vehicle. Enhancements were made to the hyperbolic grid generator HYPGEN and new tools for grid projection, manipulation, and modification, Cartesian box grid and far field grid generation and post-processing of flow solver data were developed.

Chan, William M.↗

A combined integrating and differentiating matrix formulation for boundary value problems on rectangular domains

Integrating and differentiating matrices allow the numerical integration and differential of functions whose values are known at points of a discrete grid. Previous derivations of these matrices were restricted to one dimensional grids or to rectangular grids with uniform spacing in at least one direction. Integrating and differentiating matrices were developed for grids with nonuniform spacing in both directions. The use of these matrices as operators to reformulate boundary value problems on rectangular domains as matrix problems for a finite dimensional solution vector is considered. The method requires nonuniform grids which include near boundary points. An eigenvalue problem for the transverse vibrations of a simply supported rectangular plate is solved to illustrate the method.

Lakin, W. D.↗

A combined integrating- and differentiating-matrix formulation for boundary-value problems on rectangular domains

Integrating and differentiating matrices allow the numerical integration and differential of functions whose values are known at points of a discrete grid. Previous derivations of these matrices were restricted to one dimensional grids or to rectangular grids with uniform spacing in at least one direction. Integrating and differentiating matrices were developed for grids with nonuniform spacing in both directions. The use of these matrices as operators to reformulate boundary value problems on rectangular domains as matrix problems for a finite dimensional solution vector is considered. The method requires nonuniform grids which include near boundary points. An eigenvalue problem for the transverse vibrations of a simply supported rectangular plate is solved to illustrate the method.

Lakin, W. D.↗

A Flexible Field Mapping System for Accelerator Magnets

Magnetic field mapping is a fundamental magnetic measurement method that typically uses Hall and NMR sensors. In magnet measurement facilities, such systems are likely used in various configurations suitable for a specific task at hand. To address this diversity, the authors developed a flexible field mapping system capable of being configured and tailored to each particular measurement case. Further, the system needs to address the variability introduced by differences in sensors and their readout systems, probe positioning systems, power supply systems, and required mapping geometry (mapping space and grid, measurement steps and sequences). Although the discussed field mapping systems range from a self-propelled multi-sensor mapper of a large detector magnet to a single 3D Hall sensor system to scan a small permanent magnet, they were all built with the same core mapping system. The variability present in field mapping systems, the measurement system architecture addressing this variability, as well as examples of several field mapping systems built in this architecture are presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TRT Compton Scattering Data: Comparisons between LLNL and LANL: Pointwise

This memo is one of a pair of reports comparing Compton scattering data in the context of thermal radiation transport (TRT). This document focuses on Pointwise (PW) data while the other report focuses on Multigroup (MG) data. These complementary analyses, via PW comparisons, show if the same physics are used and if fundamental bugs are present, and, via MG comparisons, quantify macroscopic Compton effects such as mean amplification factors. No averaging of Compton data over photon energy is done for the PW comparisons, while lab-dependent numerical resolution choices become relevant for the MG comparisons. This document plots and quantifies differences in PW Compton scattering data between LANL and LLNL for four disparate temperatures on a logarithmically-spaced energy grid. PW comparisons means we use the received LLNL data as-is, without any averaging or integration over energy. In addition to LANL-to-LLNL comparisons, we investigate the impact of numerical choices on LANL PW data. For the purposes of this comparison, we neglect induced scattering.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Flexible Field Mapping System for Accelerator Magnets

Magnetic field mapping is a fundamental magnetic measurement method that typically uses Hall and NMR sensors. In magnet measurement facilities, such systems are likely used in various configurations suitable for a specific task at hand. To address this diversity, the authors developed a flexible field mapping system capable of being configured and tailored to each particular measurement case. The system needs to address the variability introduced by differences in sensors and their readout systems, probe positioning systems, power supply systems and in required mapping geometry (mapping space and grid, measurement steps and sequences). Although the discussed field mapping systems range from a self-propelled, multi-sensor mapper of a large detector magnet to a single 3D Hall sensor system to scan a small permanent magnet, they were all built with the same core mapping system. The variability present in field mapping systems and the measurement system architecture addressing this variability, as well as examples of several field mapping systems built in this architecture are presented.

43 PARTICLE ACCELERATORS↗

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗

Influence of nonconservative differencing on transonic streamline shapes

A computer program recently developed by South and Brandt (1976) which contained the Murman (1973) conservative finite-difference scheme is easily modified to use the Garabedian and Korn (1971) nonconservative finite difference scheme. This program solves the transonic small disturbance equation for only symmetric flow, but incorporates several iterative solution techniques. Results are presented for the case where the equally spaced computational grid extended to infinity in both the streamwise and normal directions. Streamline shapes are obtained along several grid lines by a streamwise integration of the normal component of the perturbation velocity. Comparison cases are run for a 10% thick parabolic arc airfoil at zero incidence for freestream Mach numbers of 0, 0.70, 0.84, and 0.95. It is shown that the use of a nonconservative finite-difference scheme in transonic flow calculations destroys the global mass balance when shocks are present. This lack of mass balance may prove to be more crucial in the case of an unconfined external flow.

Newman, P. A.↗

Apparatus and method for destructive removal of particles contained in flowing fluid

An apparatus and method for destructively removing particles from a flowing gas containing the particles is described. In the specific embodiments disclosed the apparatus is adapted to remove carbon particles from diesel engine exhaust products. The exhaust products are directed to a predetermined location where they are rapidly vaporized and combine with oxygen in the exhaust products to form carbon dioxide. Vaporization in one embodiment is effected by a discharge grid located within an exhaust conduit, the grid being chosen so that alternate conductors defining the grid are spaced apart a distance approximately 125 times the mean diameter of the particles to be removed. A voltage differential of approximately 690 volts is applied across adjacent conductors.

Yang, L. C.↗

Grid generation in three dimensions by Poisson equations with control of cell size and skewness at boundary surfaces

An algorithm for generating computational grids about arbitrary three-dimensional bodies is developed. The elliptic partial differential equation (PDE) approach developed by Steger and Sorenson and used in the NASA computer program GRAPE is extended from two to three dimensions. Forcing functions which are found automatically by the algorithm give the user the ability to control mesh cell size and skewness at boundary surfaces. This algorithm, as is typical of PDE grid generators, gives smooth grid lines and spacing in the interior of the grid. The method is applied to a rectilinear wind-tunnel case and to two body shapes in spherical coordinates.

Sorenson, R. L.↗

Fourier analysis of the SOR iteration

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. It is shown that for the standard model problem of Poisson's equation on a rectangle, the optimal omega and corresponding convergence rate can be rigorously obtained by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives insight into the relation between convergence rates of several variants.

Leveque, R. J.↗