Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal estimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed.

Rodriguez, G.↗
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗
NASA NTRS1987

X-ray observations of spiral galaxies. II - Images and spectral parameters of 13 galaxies

The results of a spatial and spectral analysis of a sample of spiral galaxies observed in X-rays with the Einstein satellite indicate that the X-ray emission is generally extended and complex. In some cases, discrete, very bright X-ray sources are detected, some of these coincident with the nuclear regions. The spectral fit of the IPC data to a thermal bremsstrahlung spectrum with a low-energy absorption cutoff shows that the average spectrum could be rather hard with kT greater than 2 keV, and with low-energy cutoff consistent with line-of-sight absorption. Results on NGC 4631, NGC 6946, and the nuclei of IC 342 and NGC 1313 are discussed and compared with previous observations in different energy ranges.

Fabbiano, G.↗
NASA NTRS1989

Photon counting detectors for Fabry-Perot interferometers

Sealed channel plate photomultipliers with multiple discrete anodes for use as photon counting detectors in the image plane of Fabry-Perot interferometers are described. The influence of design and construction on performance of completed devices is discussed. Effects on spatial resolution, lifetime, and counting efficiency are described. It is shown that devices can be optimized for particular applications. The results should be generally applicable to resistive anode and wedge and strip anode types of sealed detectors.

Darlington, E. H.↗
DOE OSTI2025

Spatial mapping of dissolved methane using an in situ sensor in Puget Sound

Release of methane, as gas bubbles or in the dissolved phase, from the seafloor has been observed in coastal waters (< 200 m) and deep ocean basins (> 1000 m). Methane dissolution within the water column affects the geochemistry of the surrounding water, leading to localized oxygen loss and potential escape to the atmosphere, particularly from shallower sites. Traditional methods for detecting and quantifying dissolved methane rely on collecting discrete water samples for ship- or land-based ex situ analysis and post processing. Here, we report on the use of a reduced response time, in situ methane sensor, the Sensor for Aqueous Gases in the Environment (SAGE), for detecting and quantifying dissolved methane concentrations in a wide range of seafloor environments. During a Fall 2022 research cruise on the R/V Thomas G. Thompson in Puget Sound, SAGE was integrated onto a towed conductivity/temperature/depth rosette and deep-sea camera system with live-stream 1 Hz telemetry and used to spatially map the concentration of methane approximately 1 m above the seafloor. The site had been previously identified as an active methane plume field characterized by gas bubbles, fluid venting, and a faulted seabed. The widespread background dissolved concentration of methane measured by SAGE was 83 nM, and a range of 78–670 nM was observed throughout the survey. The results highlight the capacity of SAGE to map the spatial and temporal variability of dissolved methane concentrations in situ and to identify and localize sites of variable methane emissions from the seafloor.

Padilla, Alexandra M. [Woods Hole Oceanographic In↗
DOE OSTI2021

On a nonlocal Cahn–Hilliard model permitting sharp interfaces

A nonlocal Cahn–Hilliard model with a non-smooth potential of double-well obstacle type that promotes sharp interfaces in the solution is presented. To capture long-range interactions between particles, a nonlocal Ginzburg–Landau energy functional is defined which recovers the classical (local) model as the extent of nonlocal interactions vanish. In contrast to the local Cahn–Hilliard problem that always leads to diffuse interfaces, the proposed nonlocal model can lead to a strict separation into pure phases of the substance. In this work, the lack of smoothness of the potential is essential to guarantee the aforementioned sharp-interface property. Mathematically, this introduces additional inequality constraints that, in a weak formulation, lead to a coupled system of variational inequalities which at each time instance can be restated as a constrained optimization problem. We prove the well-posedness and regularity of the semi-discrete and continuous in time weak solutions, and derive the conditions under which pure phases are admitted. Moreover, we develop discretizations of the problem based on finite element methods and implicit–explicit time-stepping methods that can be realized efficiently. Finally, we illustrate our theoretical findings through several numerical experiments in one and two spatial dimensions that highlight the differences in features of local and nonlocal solutions and also the sharp interface properties of the nonlocal model.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

A source map for dust jets observed in the coma of Comet P/Halley

In this paper, the large-scale enhancements of about 50 s duration in the flux of submicron dust particles from Comet P/Halley recorded by the Dust Counter and Mass Analyzer carried aboard the two Vega spaceprobes are traced back to their source departure points on the nucleus. Evidence is presented that such enhancements which correspond to regions of enhanced particle density in the coma with spatial widths of about 4000 km, are superimposed upon a background flux distribution which is well described by a simple fountain model. Thus, dust production on the sunlit hemisphere of the nucleus is characterized by discrete jet sources superimposed upon relatively low-level isotropic emission. The jet sources at the nucleus responsible for the major enhancements in recorded particle impact rates are narrow and discrete, with each jet accelerating about 10 to the -13th g particles to a narrow range of terminal velocities which decrease with increasing particle mass.

Rabinowitz, D. L.↗
NASA NTRS1991

Time-dependent aerodynamic analysis of ducted and unducted propfans at angle of attack

A three-dimensional unsteady aerodynamic analysis is described for predicting the time-dependent flow about ducted and unducted propfans operating at angle of attack. Although the freestream is assumed to be uniform, the flow relative to the rotating blades varies with circumferential position, resulting in an inherent unsteadiness due to the nonaxial inflow. The time-dependent Euler equations are solved utilizing a Runge-Kutta time-stepping scheme. The analysis is based on a finite-volume discretization employing a multiple-block grid network. To permit the use of large calculation time steps, an implicit residual smoothing scheme previously tested for unsteady flow calculations in two dimensions is extended to three spatial dimensions. For unducted propfans, a single H-type grid block is used for each blade passage to determine the time-periodic flowfield. For ducted propfans (ultra-high bypass fans) a body-centered C-type grid is wrapped about the cowl to improve the accuracy of the analysis in the high gradient flow region near the cowl leading edge. Numerical results are compared with available data for both ducted and unducted propfans operating at angle of attack.

Hall, Edward J.↗
NASA NTRS1993

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: A careful study of the boundary error

The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Carpenter, Mark H.↗
DOE OSTI2025

Introducing the SLICE Method for estimating pebble-bed reactor inventories at equilibrium operation with SCALE

This paper introduces the SCALE Leap-In method for Cores at Equilibrium (SLICE) for estimating pebble-bed reactor equilibrium core isotopic inventories using capabilities in the SCALE code system, requiring only a small computational cluster and a few days of computation. This method uses an iterative approach that relies on (1) a surrogate spectrum model that captures spatial and time-dependent spectral conditions, (2) a multi-pass model that captures the pebble’s evolving nuclide inventory as a function of location and time in the core, and (3) a full-core model that captures the core’s spatial neutron flux distribution. The SLICE approach is applied to a generic fluoride salt–cooled high-temperature reactor, demonstrating fuel inventory convergence through nuclide concentration inspection across iterations and comparisons for core realizations with varying discretizations. Results agree within ~5% with another state-of-the-art code, with differences attributed to input parameter or modeling assumption variations in the equilibrium generation methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A low-latitude, high-elevation 10Be production-rate calibration from the Rwenzori Mountains, Uganda

Cosmogenic nuclide surface-exposure dating has become a key tool in geomorphology. Robust application of the technique relies on determination of local nuclide production rates using independently dated calibration sites. However, there are few existing calibration sites in certain regions, particularly in low-latitude, high-elevation regions where cosmogenic nuclide production is theoretically more sensitive to changes in Earth's magnetic field and atmospheric conditions. We present a new cosmogenic beryllium-10 (10Be) production-rate calibration from the equatorial Rwenzori Mountains in Uganda. The calibration is based on radiocarbon dating of basal lake sediments (∼21.3 kyr BP) from Lake Mahoma (2900 m asl) and measured 10Be concentrations from boulders on the crest of a moraine which dams the lake. Our results indicate a time-integrated 10Be production rate of 15.83 ± 0.81 at/g/yr at Lake Mahoma. We then use these calibration data and two public online calculators to calculate ages for moraines in a separate Rwenzori valley for which there exists a minimum-limiting radiocarbon age constraint (∼12 kyr BP) on surficial glacial sediments. This enables us to assess the performance of different parameters used to scale production rates spatially and temporally (e.g., scaling frameworks, geomagnetic field reconstructions, atmospheric models) for surface-exposure age calculation in the region. Our results highlight the sensitivity of low-latitude cosmogenic nuclide surface-exposure ages to discrete calculation parameters, such as scaling framework or reconstructions of Earth's magnetic field, and the utility of determining production rates at sites where limiting age data exists for other nearby landforms.

Anderson, N [Dartmouth College]↗
NASA NTRS1982

Satellite observations of suprathermal electron bursts

Data from the ISIS 2 satellite are used to study the spatial distribution of the suprathermal electron bursts discovered by high-altitude sounding rockets. It is found that the bursts form an integral part of the boundary plasma sheet and are linked with the formation of discrete auroral arcs. The maximum energy of the bursts is always below the peak energy in neighboring discrete forms, and it appears that if the bursts fail to reach the minimum energy found to be necessary to sustain a discrete arc they collapse within a few seconds. It is suggested that in a pass across the auroral oval the overall envelope of the maximum energy in bursts and the peak energy in discrete auroral forms reflect the variation of the available acceleration (perhaps a parallel electric potential) across the auroral oval while the energy structure in individual forms is a function of the nonlinear processes which sustain them. This explanation appears to resolve the apparent discrepancy in scale thickness between rocket and satellite inverted V's.

Johnstone, A. D.↗
NASA NTRS1987

Evolution of a localized thermal explosion in a reactive gas

Experimental observations of ignition in premixed gaseous reactants indicate that perfectly homogeneous initiation is practically unrealizable. Instead, combustion first sets in, as a rule, at small, discrete sites where inherent inhomogeneities cause chemical activity to proceed preferentially and lead to localized explosions. Combustion waves propagating away from these hot spots or reaction centers eventually envelop the remaining bulk. This study examines the spatial structure and temporal evolution of a hot spot for a model involving Arrhenius kinetics. The hot spot, characterized by peaks in pressure and temperature with little diminution in local density, is shown to have one of two possible self-similar structures. The analysis employs a combination of asymptotics and numerics, and terminates when pressure and temperature in the explosion have peaked.

Jackson, T. L.↗
NASA NTRS2019

Assessment of Buffet Forcing Function Development Process Using Unsteady Pressure Sensitive Paint

A wind tunnel test was conducted at the Ames Unitary Plan Wind Tunnel to characterize the transonic buffet environment of a generic launch vehicle forebody. The test examined a highly instrumented version of the Coe and Nute Model 11 test article first tested in the 1960s. One of the measurement techniques used during this test was unsteady pressure sensitive paint (uPSP) developed at the Arnold Engineering Development Complex. This optical measurement technique measured fluctuating pressures at over 300,000 locations on the surface of the model. The high spatial density of these measurements provided an opportunity to examine in depth the assumptions underpinning the development of buffet forcing functions (BFFs) used in the development of the Space Launch System vehicle. The comparison of discrete-measurement-based BFFs to BFFs developed by continuous surface pressure integration indicates that the current BFF development approach under predicts low frequency content of the BFFs while over predicting high frequency content. Coherence-based adjustments employed to reduce over prediction in the surface integration of discrete pressure measurements contribute to the inaccuracy of the BFFs and their implementation should be reevaluated.

Sekula, Martin K.↗
DOE OSTI2022

Improvements to the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. The code includes a variety of steady-state solvers for fixed-source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. This document summarizes the transport solver development efforts pursued during Fiscal Year 2022. We added the multiphysics transient capability for the coarse-mesh finite difference accelerated Richardson iteration for discontinuous finite element method discrete ordinates (DFEM-SN) scheme to support high-order heterogeneous transport simulations. HFEM (hybrid finite element method) - PN (spherical harmonics expansion) was completed and red-black iteration was added for solving the HFEM-PN system with both preconditioned Jacobian-free Newton Krylov and Richardson iteration solvers. The HFEM-PN scheme, as one of the low-order transport schemes, is expected for supporting routine design simulations. Pin power reconstruction capability was also designed and implemented with the Griffin ISOXML module to enhance all the low-order transport solvers for more accurate multiphysics simulations. Numerical results are presented for demonstrating the capabilities and verifying their performance, and future works are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MOOSE Optimization Module Overview: Application to Residual Stress Inversion in Nuclear Fuel Plates

Manufacturing methods, material selections and unique plate geometries used in low enriched uranium plate fuels for the High Performance Research Reactor (HPRR) have led to concerns about the residual stresses effect on bond strength of the fuel/cladding interfaces and overall fuel performance under irradiation. The current HPRR plate fuel test program performed residual stress calculations using the contour method with experimental displacement data obtained by incrementally cutting the plate fuel. In this work, we use a gradient based inverse optimization module developed within the MOOSE framework to estimate normal and shear residual stresses within the plate fuels using displacement data taken from a cut placed in the HPRR fuel plate. These results are verified against residual stress calculations obtained from the contour method. The inversion protocol used in this work utilizes two meshes to separately discretize the physics and parameter space. This allows a fine mesh to be used to fully resolve the physics and geometry of the forward problem. A coarser mesh is used to resolve the parameter space with spatially varying mesh density determined from a sensitivity analysis. The adaptive resolution of the parameter space essentially has the effect of regularization as it is designed to achieve the bias-variance tradeoff.

42 ENGINEERING↗
DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗