Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial refinement”

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 91 records · Page 5

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Using Raman scattering to study reacting gas flow over a catalyst

l present a series of computational methods for analyzing and correcting Raman scattering from a reacting gas flow. My work focused on the laser diagnostic measurements of formaldehyde production over a silver catalyst, specifically the Raman spectroscopy component. This technique offers a non-intrusive, in situ way to measure the coupling between the gas phase and surface activity of a catalytic reaction at relevant industrial conditions. This report outlines the steps followed for correcting the non-idealities in the data due to both the detection optics and scattering mechanics. By implementing such refinement steps, l can generate an accurate 2D spatial map of the temperature profile, number densities for major species, and relative consumption of reactants throughout the reactor cell.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Analysis of Published DAS Studies for Application to SPE Phase III

Distributed Acoustic Sensing (DAS) is an emerging technology capable of recording the acoustic wavefield at unprecedented spatial resolution. However, this new tool requires significant refinements before it becomes operational for explosion monitoring objectives. Recent studies have shown significant development of array processing with DAS data. In this report we explore three such array processing methods including DAS strain-rate data versus geophone measured ground motion, beamforming for event parameters, and machine learning based denoising. Prior to applying these algorithms to the Source Physics Experiment Phase II and Phase III data, we validate these methods through replication analysis.

47 OTHER INSTRUMENTATION↗

Parallel, adaptive finite element methods for conservation laws

We construct parallel finite element methods for the solution of hyperbolic conservation laws in one and two dimensions. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. A posteriori estimates of spatial errors are obtained by a p-refinement technique using superconvergence at Radau points. The resulting method is of high order and may be parallelized efficiently on MIMD computers. We compare results using different limiting schemes and demonstrate parallel efficiency through computations on an NCUBE/2 hypercube. We also present results using adaptive h- and p-refinement to reduce the computational cost of the method.

Biswas, Rupak↗

Global Flowfield About the V-22 Tiltrotor Aircraft

This final report includes five publications that resulted from the studies of the global flowfield about the V-22 Tiltrotor Aircraft. The first of the five is 'The Chimera Method of Simulation for Unsteady Three-Dimensional Viscous Flow', as presented in 'Computational Fluid Dynamics Review 1995.' The remaining papers, all presented at AIAA conferences, are 'Unsteady Simulation of the Viscous Flow About a V-22 Rotor and Wing in Hover', 'An Efficient Means of Adaptive Refinement Within Systems of Overset Grids', 'On the Spatial and Temporal Accuracy of Overset Grid Methods for MOving Body Problems', and 'Moving Body Overset Grid Methods for Complete Aircraft Tiltrotor Simulations.'

Meakin, Robert L.↗

Integrating Ideas for International Data Collaborations Through The Committee on Earth Observation Satellites (CEOS) International Directory Network (IDN)

The capabilities of the International Directory Network's (IDN) version MD9.5, along with a new version of the metadata authoring tool, "docBUILDER", will be presented during the Technology and Services Subgroup session of the Working Group on Information Systems and Services (WGISS). Feedback provided through the international community has proven instrumental in positively influencing the direction of the IDN s development. The international community was instrumental in encouraging support for using the IS0 international character set that is now available through the directory. Supporting metadata descriptions in additional languages encourages extended use of the IDN. Temporal and spatial attributes often prove pivotal in the search for data. Prior to the new software release, the IDN s geospatial and temporal searches suffered from browser incompatibilities and often resulted in unreliable performance for users attempting to initiate a spatial search using a map based on aging Java applet technology. The IDN now offers an integrated Google map and date search that replaces that technology. In addition, one of the most defining characteristics in the search for data relates to the temporal and spatial resolution of the data. The ability to refine the search for data sets meeting defined resolution requirements is now possible. Data set authors are encouraged to indicate the precise resolution values for their data sets and subsequently bin these into one of the pre-selected resolution ranges. New metadata authoring tools have been well received. In response to requests for a standalone metadata authoring tool, a new shareable software package called "docBUILDER solo" will soon be released to the public. This tool permits researchers to document their data during experiments and observational periods in the field. interoperability has been enhanced through the use of the Open Archives Initiative s (OAI) Protocol for Metadata Harvesting (PMH). Harvesting of XML content through OAI-MPH has been successfully tested with several organizations. The protocol appears to be a prime candidate for sharing metadata throughout the international community. Data services for visualizing and analyzing data have become valuable assets in facilitating the use of data. Data providers are offering many of their data-related services through the directory. The IDN plans to develop a service-based architecture to further promote the use of web services. During the IDN Task Team session, ideas for further enhancements will be discussed.

Olsen, Lola M.↗

The Gravity Field, Orientation, and Ephemeris of Mercury from MESSENGER Observations After Three Years in Orbit

We have analyzed three years of radio tracking data from the MESSENGER spacecraft in orbit around Mercury and determined the gravity field, planetary orientation, and ephemeris of the innermost planet. With improvements in spatial coverage, force modeling, and data weighting, we refined an earlier global gravity field both in quality and resolution, and we present here a spherical harmonic solution to degree and order 50. In this field, termed HgM005, uncertainties in low-degree coefficients are reduced by an order of magnitude relative to the earlier global field, and we obtained a preliminary value of the tidal Love number k(sub 2) of 0.451+/-0.014. We also estimated Mercury's pole position, and we obtained an obliquity value of 2.06 +/- 0.16 arcmin, in good agreement with analysis of Earth-based radar observations. From our updated rotation period (58.646146 +/- 0.000011 days) and Mercury ephemeris, we verified experimentally the planet's 3: 2 spin-orbit resonance to greater accuracy than previously possible. We present a detailed analysis of the HgM005 covariance matrix, and we describe some near-circular frozen orbits around Mercury that could be advantageous for future exploration.

MESSENGER↗

Behavioral simulations and experimental evaluations of stress induced spatial nonuniformity of dynamic bulk modulus in additive friction stir deposited AA 6061

Additive friction stir deposition (AFSD) is a rapidly developing solid-state additive manufacturing (AM) technology that enables the fabrication of customized products on a relatively large scale. The unique forge-liked process offers high quality products with refined equiaxed microstructure and exceptional mechanical performance. However, in AFSD, the spatially distributed mechanical properties within a deposited layer have not been established so far, due to the lack of a suitable large-scale mapping method. In the present work, an ultrasonic dynamic-elasticity mapping method was employed to map the distribution of physical properties in a AFSD 6061 Al. Strong spatial asymmetries in elastographies were mainly due to the residual stresses. The origin of the asymmetries was investigated through thermo-mechanical simulations of the AFSD process. It was observed that majority of the asymmetries were associated with the mechanical residual stresses instead of the symmetric thermal residual stresses during linear deposition of the AFSD process. However, the thermal residual stresses can be manipulated to become asymmetric with a 90° turn in the linear deposition path, due to the difference between the levels of cooling along the edges of the inner and outer corners. Hence, the mechanical and thermo-mechanical residual stresses can be additive or subtractive depending on the direction of the turns. Here, such numerically predicted behaviors were also experimentally observed on the AFSD fabricated samples by ultrasonic elastography. Furthermore, additional preliminary studies on multi-layer AFSD fabrications show the asymmetric mechanical stresses generated during the deposition process on different layers generate constructive or destructive interference which results to spatially variable dynamic bulk modulus depending on the deposition orientations of layers.

36 MATERIALS SCIENCE↗

GPS Altimetry

The advent of satellite altimetry has greatly improved our ability to observe global ocean circulation. However, the swath of a single, nadir-viewing satellite altimeter is only a few km and the track spacing is several hundred km to resolve the two-dimensional structure of ocean eddies. Our goal is to increase spatial and temporal coverage by monitoring Global Positioning System (GPS) signals reflected from the ocean. A constellation of spacecraft would each carry a GPS receiver capable of recording 8 reflections simultaneously. The reflections are well distributed in azimuth and elevation and can be tracked continuously while the satellite is in view, and another is then acquired, as illustrated below. The diagram depicts a new approach at altimetry measurements where ocean surface reflected GPS signals are simultaneously tracked and processed in a GPS flight receiver in space. The reflected GPS signals from the ocean must be compared precisely with the direct GPS signals in order to infer the characteristics of the ocean from the combined data set. Understanding the features and accuracy of GPS altimetry measurement is crucial to establishing its suitability for oceanography. Preliminary work has enabled us to theoretically model the signal output of the correlator for a variety of system parameters such as wind speed (sea roughness), receiver height, incidence angle, receiver range and Doppler filter bandwidth and antenna gain. Expected signal-to-noise ratio has been estimated from which we have inferred, to a first approximation, the basic receiver gain requirements for a space-based altimeter and the expected range raw error. In 1998, work on a different task led to the extraction of the first reflected GPS signal observed from a spaceborne receiver during the 1995 Space Transportation System-68 (STS-68) Shuttle Radar Laboratory-2 (SRL-2) high resolution synthetic aperture radar mission. Good comparisons with our signal models have been obtained. Having established that only modest signal-to-noise ratios are obtainable unless very high gain antennas are used, we expect that single measurements might not provide an estimate of sea state parameters as accurate as that obtainable with traditional remote sensing instruments, if costs are limited. Therefore spatial and temporal averaging of many measurements is required. Since the receiving satellite tracks do not repeat, measurements in a given area will be used to refine the local solution as a function of time, and will define the spatial resolution. In order to make this instrument viable in space, we need to detect and process many scattered signals. In 1998 we made a first step at understanding what configurations of antenna gain and orientations capture the largest number of viable signals to be used in the subsequent spatial/temporal averaging process. Additional information is contained in the original.

Lichten, Stephen M.↗

New modelling capabilities in IDT

This work concerns the enhanced modelling capabilities of the discrete ordinates transport solver IDT. The novelties introduced allow for modelling unstructured geometries composed by a collection of X/Y segments and circles, and the use of reciprocity and conservation relations reduce the memory imprint as well as the computational cost of the method. IDT decomposes geometries in modular Cartesian patterns, which are the so-called Heterogeneous Cartesian Cells (HCCs), containing a chunk of the original unstructured geometries. Each HCC can be then discretized by superimposing a XY grid to refine locally the HCC. Unlike the most popular MOC, IDT performs the spatial sweeping by directional collision probabilities instead of trajectories. The sources and interface angular fluxes are expanded up to linear order. The accuracy of ray-tracing, the memory imprints together with the novel mesh refinement capabilities have been verified. A first set of preliminary results on PWR lattice problems will be presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

AMReX: Block-structured adaptive mesh refinement for multiphysics applications

Block-structured adaptive mesh refinement (AMR) provides the basis for the temporal and spatial discretization strategy for a number of Exascale Computing Project applications in the areas of accelerator design, additive manufacturing, astrophysics, combustion, cosmology, multiphase flow, and wind plant modeling. AMReX is a software framework that provides a unified infrastructure with the functionality needed for these and other AMR applications to be able to effectively and efficiently utilize machines from laptops to exascale architectures. AMR reduces the computational cost and memory footprint compared to a uniform mesh while preserving accurate descriptions of different physical processes in complex multiphysics algorithms. AMReX supports algorithms that solve systems of partial differential equations in simple or complex geometries and those that use particles and/or particle–mesh operations to represent component physical processes. In this article, we will discuss the core elements of the AMReX framework such as data containers and iterators as well as several specialized operations to meet the needs of the application projects. In addition, we will highlight the strategy that the AMReX team is pursuing to achieve highly performant code across a range of accelerator-based architectures for a variety of different applications.

Zhang, Weiqun↗

Four Years of Atmospheric Boundary Layer Height Retrievals Using COSMIC-2 Satellite Data

This work aimed to study the atmospheric boundary layer height (ABLH) from COSMIC-2 refractivity data, endeavoring to refine existing ABLH detection algorithms and scrutinize the resulting spatial and seasonal distributions. Through validation analyses involving different ground-based methodologies (involving data from lidar, ceilometer, microwave radiometers, and radiosondes), the optimal ABLH determination relied on identifying the lowest refractivity gradient negative peak with a magnitude at least $τ$% times the minimum refractivity gradient magnitude, where $τ$ is a fitting parameter representing the minimum peak strength relative to the absolute minimum refractivity gradient. Different $τ$ values were derived accounting for the moment of the day (daytime, nighttime, or sunrise/sunset) and the underlying surface (land or sea). Results show discernible relations between ABLH and various features, notably, the land cover and latitude. On average, ABLH is higher over oceans (≈1.5 km), but extreme values (maximums > 2.5 km, and minimums < 1 km) are reached over intertropical lands. Variability is generally subtle over oceans, whereas seasonality and daily evolution are pronounced over continents, with higher ABLHs during daytime and local wintertime (summertime) in intertropical (middle) latitudes.

54 ENVIRONMENTAL SCIENCES↗

High-speed compressible flow and other advection-dominated problems of fluid dynamics

Finite element methods are described for modeling high speed compressible flows with strong advection, problems important to aerodynamics. The situations are characterized by high pressure and temperature gradients, transients and the appearance of discontinuities, factors which require mesh refinement during computations. Techniques are developed for temporal and spatial discretization of a model problem. Several observations are made regarding the explicit and implicit features of the calculations, the use of the Lax-Wendroff scheme to produce a mass-matrix for obtaining accurate results for transients, methods of performing stability analyses, and simplification techniques. Examples are provided of solving the nonlinear shallow-water equations and describing compressible flows, particularly transonic flows. Domain splitting is defined for improving the calculations at each time step and in different parts of the flow regime while simultaneously advancing the calculations towards a solution.

Zienkiewicz, O. C.↗

Remote sensing information sciences research group: Browse in the EOS era

The problem of science data browse was examined. Given the tremendous data volumes that are planned for future space missions, particularly the Earth Observing System in the late 1990's, the need for access to large spatial databases must be understood. Work was continued to refine the concept of data browse. Further, software was developed to provide a testbed of the concepts, both to locate possibly interesting data, as well as view a small portion of the data. Build II was placed on a minicomputer and a PC in the laboratory, and provided accounts for use in the testbed. Consideration of the testbed software as an element of in-house data management plans was begun.

Estes, John E.↗

Twenty Years of Progress on Global Ocean Tides: The Impact of Satellite Altimetry

At the dawn of the era of high-precision altimetry, before the launch of TOPEX/Poseidon, ocean tides were properly viewed as a source of noise--tidal variations in ocean height would represent a very substantial fraction of what the altimeter measures, and would have to be accurately predicted and subtracted if altimetry were to achieve its potential for ocean and climate studies. But to the extent that the altimetry could be severely contaminated by tides, it also represented an unprecedented global-scale tidal data set. These new data, together with research stimulated by the need for accurate tidal corrections, led to a renaissance in tidal studies in the oceanographic community. In this paper we review contributions of altimetry to tidal science over the past 20 years, emphasizing recent progress. Mapping of tides has now been extended from the early focus on major constituents in the open ocean to include minor constituents, (e.g., long-period tides; non-linear tides in shelf waters, and in the open ocean), and into shallow and coastal waters. Global and spatially local estimates of tidal energy balance have been refined, and the role of internal tide conversion in dissipating barotropic tidal energy is now well established through modeling, altimetry, and in situ observations. However, energy budgets for internal tides, and the role of tidal dissipation in vertical ocean mixing remain controversial topics. Altimetry may contribute to resolving some of these important questions through improved mapping of low-mode internal tides. This area has advanced significantly in recent years, with several global maps now available, and progress on constraining temporally incoherent components. For the future, new applications of altimetry (e.g., in the coastal ocean, where barotropic tidal models remain inadequate), and new mission concepts (studies of the submesoscale with SWOT, which will require correction for internal tides) may bring us full circle, again pushing further development of tidal models as corrections.

Egbert, Gary↗

A robust fourth-order finite-difference discretization for the strongly anisotropic transport equation in magnetized plasmas

We propose a second-order temporally implicit, fourth-order-accurate spatial discretization scheme for the strongly anisotropic heat transport equation characteristic of hot, fusion-grade plasmas. Following Du Toit et al. (2018), the scheme transforms mixed-derivative diffusion fluxes (which are responsible for the lack of a discrete maximum principle) into nonlinear advective fluxes, amenable to nonlinear-solver-friendly monotonicity-preserving limiters. The scheme enables accurate multi-dimensional heat transport simulations with up to seven orders of magnitude of heat-transport-coefficient anisotropies with low cross-field numerical error pollution and excellent algorithmic performance, with the number of linear iterations scaling very weakly with grid resolution and grid anisotropy, and scaling with the square-root of the implicit timestep. We propose a multigrid preconditioning strategy based on a lower-order approximation that renders the scheme efficient and scalable under grid refinement. Several numerical tests are presented that display the expected spatial convergence rates and strong algorithmic performance, including fully nonlinear magnetohydrodynamics simulations of kink instabilities in a Bennett pinch in 2D helical geometry and of ITER in 3D toroidal geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GW-PINN: A deep learning algorithm for solving groundwater flow equations

Machine learning methods provide new perspective for more convenient and efficient prediction of groundwater flow. In this study, a deep learning method “GW-PINN” without labeled data for solving groundwater flow equations with wells was proposed. GW-PINN takes the physics inform neural network (PINN) as the backbone and uses either the hard or soft constraint in the loss function for training. A locally refined sampling strategy (LRS) is adopted to generate the consistent spatial sampling points for problems with strong hydraulic head change, and then combined with an appropriate temporal sampling scheme to obtain the final spatial-temporal sampling points. A snowball-style two-stage training strategy by dividing the temporal domain into two subdomains is designed to decrease the sampling points. Five cases were designed to test the training performance of GW-PINN under different sampling strategies and two constraints. The predicted results of GW-PINN were compared with MODFLOW and the analytical solution. The results demonstrate that GW-PINN possesses strong ability in capturing the hydraulic head change for both confined and un-confined aquifers. The hard constraint owns more robust learning ability than the soft constraint. The LRS strategy can generate more accurate results with much fewer sampling points than traditional sampling strategies, and the snowball-style two-stage training strategy is significantly efficient for problems with the drastic change of hydraulic head. Additionally, the application of GW-PINN as a surrogate model for parameterized groundwater flow equations is illustrated. This study provides an option tool for efficient groundwater flow simulation, especially for those with local refinements are needed.

54 ENVIRONMENTAL SCIENCES↗

Laser terminal attitude determination via autonomous star tracking

Freespace laser communication systems require an initial spatial acquisition process to coalign the cooperating terminal line of sites (typically to sub-microrad accuracy). The complexity and time required for this process can grow dramatically with the size of the host platform(s) initial uncertainty field-of-view. The feasibility of significantly refining the laser terminal attitude knowledge prior to the initiation of the spatial acquisition process via a preliminary star tracking operation is addressed. It is shown that currently planned laser terminal designs are capable of supporting this star tracking function with minimal hardware augmentation. An analysis of stellar spectral irradiance is presented which enables the computation of the tracking noise equivalent angle versus stellar class and visual magnitude.

Chapman, William↗