SEARCH · Search NASA

Results for “Local Method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Identifying and classifying localized states in gapless systems using pseudospectral methods

For decades, band theory has been used to understand the physical properties of crystals by leveraging these systems’ translational symmetries to determine the spectral range and distribution of their Bloch eigenstates. However, many important phenomena only manifest at a crystal’s defects or boundaries where this translational symmetry is broken, such as defect-localized states that appear at color centers, or boundary-localized states whose existence is guaranteed by topology. Although it is possible to extend band theoretic approaches to study these classes of states, these extensions are only applicable in gapped systems, where any localized states are spectrally isolated from the crystal’s bulk states.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

Local and overall aerodynamic coefficients for bodies in hypersonic, rarefied flow

A computational method is given for the prediction of local pressure and viscous shear stress on windward surfaces of bluff, convex, axisymmetric or quasi-axisymmetric, hypersonic bodies in the transitional, rarefied flow regime. Overall aerodynamic forces and moments are computed by integration of the local quantities. The method is based upon a correlation of local pressure and shear stress computed by the direct simulation Monte Carlo (DSMC) numerical technique for cold wall, real gas conditions and some supplemental data from low-density, hypersonic wind tunnels. The relative simplicity of the method makes it feasible to do the necessary calculations with a personal computer. Two-dimensional shapes and leeward surfaces are not included in the scope of the method as it is presented here. Results are compared with DSMC computations for both local and overall coefficients. The latter includes sphere and blunt cone drag as well as lift and pitching moment coefficients for the NASA AFE vehicle at various angles of attack. Very satisfactory agreement is shown.

Potter, J. Leith↗

DOE OSTI2021

Machine learning to identify geologic factors associated with production in geothermal fields: a case-study using 3D geologic data, Brady geothermal field, Nevada

Abstract In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fluid-flow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fluid-flow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k -means clustering (NMF k ), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes.

58 GEOSCIENCES↗

DOE OSTI2021

Machine Learning to Identify Geologic Factors Associated with Production in Geothermal Fields: A Case-Study Using 3D Geologic Data from Brady Geothermal Field and NMFk

In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fuid-fow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fuid-fow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes. This submission includes the published journal article detailing this work, the published 3D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using the GeoThermalCloud framework, which is part of SmartTensors (both are linked below).

15 GEOTHERMAL ENERGY↗

NASA NTRS2001

Signal-to-noise ratio comparison of encoding methods for hyperpolarized noble gas MRI

Some non-Fourier encoding methods such as wavelet and direct encoding use spatially localized bases. The spatial localization feature of these methods enables optimized encoding for improved spatial and temporal resolution during dynamically adaptive MR imaging. These spatially localized bases, however, have inherently reduced image signal-to-noise ratio compared with Fourier or Hadamad encoding for proton imaging. Hyperpolarized noble gases, on the other hand, have quite different MR properties compared to proton, primarily the nonrenewability of the signal. It could be expected, therefore, that the characteristics of image SNR with respect to encoding method will also be very different from hyperpolarized noble gas MRI compared to proton MRI. In this article, hyperpolarized noble gas image SNRs of different encoding methods are compared theoretically using a matrix description of the encoding process. It is shown that image SNR for hyperpolarized noble gas imaging is maximized for any orthonormal encoding method. Methods are then proposed for designing RF pulses to achieve normalized encoding profiles using Fourier, Hadamard, wavelet, and direct encoding methods for hyperpolarized noble gases. Theoretical results are confirmed with hyperpolarized noble gas MRI experiments. Copyright 2001 Academic Press.

Non-NASA Center↗

NASA NTRS2001

Meshless Petrov-Galerkin Method Applied to Axisymmetric Problems

An axisymmetric Meshless Local Petrov-Galerkin (MLPG) algorithm is presented for the potential and elasticity problems. In this algorithm the trial and test functions are chosen from different spaces. By a judicious choice of these functions, the integrals involved in the weak form can be restricted to a local neighborhood. This makes the method truly meshless. The MLPG algorithm is used to study various potential and elasticity problems for which exact solutions are available. The sensitivity and effectiveness of the MLPG algorithm to various parameters such as the weight functions, basis functions and support domain radius, etc. was studied. The MLPG algorithm yielded accurate solutions for all weight functions, basis functions and support domain radii considered for all of the problems studied.

Raju, I. S.↗

NASA NTRS2017

Insights into the Behavior of Potential Structural Failures Originating from Localized High Stress Regions in Configurations Relevant to Solid Rocket Motor Nozzles

During the structural certification effort for the Space Launch System solid rocket booster nozzle, it was identified that no consistent method for addressing local negative margins of safety in non-metallic materials had been developed. Relevant areas included bond-line terminations and geometric features in the composite nozzle liners. In order to gain understanding, analog test specimens were designed that very closely mimic the conditions in the actual full scale hardware. Different locations in the nozzle were represented by different analog specimen designs. This paper describes those tests and corresponding results. Finite element analysis results for the tests are presented. Strain gage correlation of the analysis to the test results is addressed. Furthermore, finite fracture mechanics (a coupled stress and energy failure criterion) is utilized to predict the observed crack pop-in loads for the different configurations. The finite fracture mechanics predictions are found to be within a 10% error relative to the average measured pop-in load for each of four configurations. Initiation locations, arrest behaviors, and resistances to further post-arrest crack propagation are also discussed.

McCutcheon, David Matthew↗

DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗

DOE OSTI2022

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A kinetic-based regularization method for data science applications

We propose a physics-based regularization technique for function learning, inspired by statistical mechanics. By drawing an analogy between optimizing the parameters of an interpolator and minimizing the energy of a system, we introduce corrections that impose constraints on the lower-order moments of the data distribution. This minimizes the discrepancy between the discrete and continuum representations of the data, in turn allowing to access more favorable energy landscapes, thus improving the accuracy of the interpolator. Our approach improves performance in both interpolation and regression tasks, even in high-dimensional spaces. Unlike traditional methods, it does not require empirical parameter tuning, making it particularly effective for handling noisy data. We also show that thanks to its local nature, the method offers computational and memory efficiency advantages over Radial Basis Function interpolators, especially for large datasets.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1986

A Rapid Attachment of Strain Gages

Hand-held toroidal gun concentrates heat in localized area. New method for bonding film gages eliminates time-consuming oven curing Hand-held "gun," operates on induction heating to concentrate heat in localized area. Ferritic plate added for low-reluctance or no-reluctance surfaces.

Schott, T. D.↗

NASA NTRS1994

New multigrid approach for three-dimensional unstructured, adaptive grids

A new multigrid method with adaptive unstructured grids is presented. The three-dimensional Euler equations are solved on tetrahedral grids that are adaptively refined or coarsened locally. The multigrid method is employed to propagate the fine grid corrections more rapidly by redistributing the changes-in-time of the solution from the fine grid to the coarser grids to accelerate convergence. A new approach is employed that uses the parent cells of the fine grid cells in an adapted mesh to generate successively coaser levels of multigrid. This obviates the need for the generation of a sequence of independent, nonoverlapping grids as well as the relatively complicated operations that need to be performed to interpolate the solution and the residuals between the independent grids. The solver is an explicit, vertex-based, finite volume scheme that employs edge-based data structures and operations. Spatial discretization is of central-differencing type combined with a special upwind-like smoothing operators. Application cases include adaptive solutions obtained with multigrid acceleration for supersonic and subsonic flow over a bump in a channel, as well as transonic flow around the ONERA M6 wing. Two levels of multigrid resulted in reduction in the number of iterations by a factor of 5.

Parthasarathy, Vijayan↗

DOE OSTI2021

Engineering the beat phenomenon of quasi-Rayleigh waves for regions with minimal Surface Acoustic Wave (SAW) amplitude

Surface Acoustic Waves (SAW) are prevalent in a variety of scientific fields. Typically, the ability to manipulate these surface waves has a significant impact on the efficacy of any SAW-based device. In this paper, we present a simple method for designing localized regions that have minimal SAW amplitude. These localized regions have valley-like drops in SAW amplitude and are achieved by using the spatial beat phenomenon of quasi-Rayleigh waves. Controlling the beat phenomenon allows us to manipulate where and how these localized regions appear in the waveguide. In particular, the localized regions can be designed to have variable widths and can be positioned consecutively due to minimal wave scattering. The design methodology and physics are described in detail, and several supplementary simulations are provided to demonstrate the efficacy of the proposed method. Finally, we briefly discuss how the proposed approach could potentially be exploited to implement a localized high-pass filter and enables SAW wave propagation that avoids surface obstacles.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method

We show that non-diffusive volume advection in two-dimensions is achieved with several benchmark problems using a newly developed high-order volume of fluids (VOF) interface reconstruction method. (1) A new VOF interface reconstruction method using circular/corner facets (linear facets are a degenerate case of arcs). We create a circular interface facet in each mixed zone by matching neighbor volume with a hybrid Newton’s-bisection method and the local solution is final. In the general case, the new VOF interface reconstruction has 3rd order accuracy and can be easily made seamless. The new method addresses intrinsic issues with Young’s method such as gaps between interface facets in the case of a curved interface, and inability to define curvature nor identify corners. (2) A non-diffusive volume advection scheme. In an ALE advection step, a well-defined interface can be carried over through a Lagrange step and used to compute volume distribution into a relaxed mesh. Then, an interface reconstruction step is performed to redefine the interface in the relaxed mesh. We must point out that the interface carried over is also a solution of interface reconstruction because all the volume fractions in the relaxed mesh are naturally matched. We provide an interface tracking method compatible with our reconstruction scheme, where it is granted to use the prior info as an initial guess to capture sub-mesh resolution features. As a result, we are able to treat multiple facets inside a single mixed cell and obtain highly accurate, non-diffusive solution for advection problems with rather coarse meshes. We show our solutions for two-dimensional incompressible flows with two materials with a) the X + O diagonal translation; b) the Zalesak rotational test; and c) the single vortex spiral test.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Three-dimensional radiative transfer on a massively parallel computer

We perform 3D radiative transfer calculations in non-local thermodynamic equilibrium (NLTE) in the simple two-level atom approximation on the Mas-Par MP-1, which contains 8192 processors and is a single instruction multiple data (SIMD) machine, an example of the new generation of massively parallel computers. On such a machine, all processors execute the same command at a given time, but on different data. To make radiative transfer calculations efficient, we must re-consider the numerical methods and storage of data. To solve the transfer equation, we adopt the short characteristic method and examine different acceleration methods to obtain the source function. We use the ALI method and test local and non-local operators. Furthermore, we compare the Ng and the orthomin methods of acceleration. We also investigate the use of multi-grid methods to get fast solutions for the NLTE case. In order to test these numerical methods, we apply them to two problems with and without periodic boundary conditions.

Vath, H. M.↗

NASA NTRS1992

An approximate local thermodynamic nonequilibrium radiation model for air

A radiatively coupled viscous shock layer analysis program which includes chemical and thermal nonequilibrium is used to calculate stagnation point flow profiles for typical aeroassisted orbital transfer vehicle conditions. Two methods of predicting local thermodynamic nonequilibrium radiation effects are used as a first and second order approximation to this phenomena. Tabulated results for both nitrogen and air freestreams are given with temperature, species, and radiation profiles for some air conditions. Two body solution results are shown for 45 and 60 degree hyperboloid bodies at 12 km/sec and 80 km altitude. The presented results constitute an advancement in the engineering modeling of radiating nonequilibrium reentry flows.

Gally, Thomas A.↗

NASA NTRS1992

An approximate local thermodynamic nonequilibrium radiation model for air

Gally, Thomas A.↗

DOE OSTI2022

Bandgap analysis and carrier localization in cation-disordered ZnGeN 2

The bandgap of ZnGeN 2 changes with the degree of cation site disorder and is sought in light emitting diodes for emission at green to amber wavelengths. By combining the perspectives of carrier localization and defect states, we analyze the impact of different degrees of disorder on electronic properties in ZnGeN 2 , addressing a gap in current studies, which largely focus on dilute or fully disordered systems. The present study demonstrates changes in the density of states and localization of carriers in ZnGeN 2 calculated using bandgap-corrected density functional theory and hybrid calculations on partially disordered supercells generated using the Monte Carlo method. We use localization and density of states to discuss the ill-defined nature of a bandgap in a disordered material and identify site disorder and its impact on the structure as a mechanism controlling electronic properties and potential device performance. Decreasing the order parameter results in a large reduction of the bandgap. The reduction in bandgap is due, in part, to isolated, localized states that form above the valence band continuum associated with nitrogen coordinated by more zinc than germanium. The prevalence of defect states in all but the perfectly ordered structure creates challenges for incorporating disordered ZnGeN 2 into optical devices, but the localization associated with these defects provides insight into the mechanisms of electron/hole recombination in the material.

36 MATERIALS SCIENCE↗