Search NASA⌕ Search

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.

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

At least 595 records · Page 33

Beyond the two-point correlation: Constraining primordial non-Gaussianity with density-perturbation moments

Constraining primordial non-Gaussianity (PNG) on the large-scale cosmic structure (LSS) is an important step in understanding properties of the early Universe, specifically in distinguishing between different inflationary models. Measuring PNG relies on evaluating the scale-dependent correlations in the density field. New summary statistics beyond the two- and three-point correlation functions in configuration space and their Fourier-space counterparts, the power- and bispectrum may provide increased sensitivity. We introduce a new method for extracting the PNG signal imprinted on the LSS by using the first three Gaussian moments of the normalized correlation in density perturbations, evaluated on varying distance scales. We aim to assess this method’s sensitivity to local PNG, parameterized by f NL . We performed spherical convolutions on a range of scales on dark-matter-halo simulations to measure the scale-dependent correlations in the density field. From these, we computed the first three moments and compared them to a model expectation vector, parameterized to the second power in f NL . Our method provides about 21% improvement in sensitivity to f NL with respect to using the two-point correlation function alone. Notably, we find that the second moment alone carries nearly as much constraining power as the mean, highlighting the potential of higher order statistics. Given its simplicity and efficiency, this framework is well suited for application to current and upcoming large-scale surveys such as the Dark Energy Spectroscopic Instrument (DESI).

early universe↗

Experimental Evaluation of Stagnation Point Collection Efficiency of the NACA 0012 Swept Wing Tip

This paper presents the experimental work of a number of icing tests conducted in the Icing Research Tunnel at NASA Glenn Research Center to develop a test method for measuring the local collection efficiency of an impinging cloud at the leading edge of a NACA 0012 swept wing and with the data obtained to further calibrate a proposed correlation for such impingement efficiency calculation as a function of the modified inertia parameter and the sweep angle. The preliminary results showed that there could be some limitation of the test method due to the ice erosion problem when encountered, and also found that, for conditions free of such problem, the stagnation point collection efficiency measurement for sweep angles up to 45 could be well approximated by the proposed correlation. Further evaluation of this correlation is recommended in order to assess its applicability for swept-wing icing scaling analysis.

Tsao, Jen-Ching↗

Dual Channel Dual Staging: Hierarchical and Portable Staging for GPU-Based In-Situ Workflow

In-situ workflows have emerged as an attractive approach for addressing data movement challenges at very large scales. Since GPU-based architectures dominate the HPC landscapes, porting these in-situ workflows, and, specifically, the inter-application data exchange, to GPU-based systems can be challenging. Technologies such as GPUDirect RDMA (GDR), which is typically used for I/O in GPU applications as an optimization that circumvents the CPU overhead, can be leveraged to support bulk data exchanges between GPU applications. However, current GDR design often lacks performance portability across HPC clusters built with different hardware configurations. Furthermore, the local CPU may also be effectively used as an auxiliary communication mechanism to offload data exchanges. In this paper, we present a dual channel dual staging approach for efficient, scalable, and performance-portable inter-application data exchange for in-situ workflows. This approach exploits the data access pattern within in-situ workflows along with the inherent execution asynchrony to accelerate data exchanges and, at the same time, improve performance portability. Specifically, the dual channel dual staging method leverages both the local CPU and the remote data staging server to build a hierarchical joint staging area and uses this staging area to transform blocking inter-application bulk data exchanges into best-effort local data movements between GPU and CPU. The dual channel dual staging is implemented as a portability extension of the Dataspaces-GPU staging framework. We present an experimental evaluation of its performance, portability, and scalability using this implementation on three leadership GPU clusters. The evaluation results demonstrate that the dual channel dual staging method saves up to 75% in data-exchange time compared to host-based, GDR, and alternate portable designs, while maintaining scalability (up to 512 GPUs) and performance portability across the three platforms.

Zhang, Bo [University of Utah]↗

Quantifying VOC-Reaction Tracers, Ozone Production, and Continuing Aerosol Production Rates in Urban and Far-Downwind Atmospheres

We have found a surprisingly informative decomposition of the complex question of smoggy ozone production (basically, [HO2] in a more locally determined field of [NO]) in the process of linked investigations of modestly smoggy Eastern North America (by NASA aircraft, July 2004) and rather polluted Flushing, NYC (Queens College, July, 2001). In both rural and very polluted situations, we find that a simple contour graph parameterization of the local principal ozone production rate can be estimated using only the variables [NO] and j(sub rads) [HCHO]: Po(O3) = c (j(sub rads) [HCHO])(sup a) [HCHO](sup b). Here j(sub rads) is the photolysis of HCHO to radicals, presumably capturing many harder-UV photolytic processes and the principle ozone production is that due to HO2; mechanisms suggest that ozone production due to RO2 is closely correlated, often suggesting a limited range of different proportionality factors. The method immediately suggests a local interpretation for concepts of VOC limitation and NOx limitation. We believe that the product j(sub rads) [HCHO] guages the oxidation rate of observed VOC mixtures in a way that also provides [HO2] useful for the principle ozone production rate k [HO2] [NO], and indeed, all ozone chemical production. The success of the method suggests that dominant urban primary-HCHO sources may transition to secondary plume-HCHO sources in a convenient way. Are there other, simple, near-terminal oxidized VOC's which help guage ozone production and aerosol particle formation? Regarding particles, we report on, to the extent NASA Research resources allow, on appealing relationships between far-downwind (Atlantic PBL) HCHO and very fine aerosol (including sulfate. Since j(sub rads) [HCHO] provides a time-scale, we may understand distant-plume particle production in a more quantitative manner. Additionally we report on a statistical search in the nearer field for relationships between glyoxals (important near-terminal aromatic and isoprene reaction products) and aerosol production, looking for VOC's that might be most implicated. All three variables j(sub rads), [HCHO], and [NO] are relatively easily measured in widespread air pollution monitoring networks, and all are deducible form space-borne observations, though estimation of [NO] from [NO2] (the species observable from space) may require care. We report also on airborne and surface observations of HCHO, suggesting that concentrated (urban) and more diffuse (forest) sources may be distinguishable from space. The use of the 3.58 micron microwindow for HCHO remote sensing should allow much sharper resolution of HCHO than the UV. UV sensing requires large and expensive instruments, but even these seem justified since formaldehyde is so informative.

Chatfield, Robert↗

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

Kaljevic, Igor↗

Development and application of a local linearization algorithm for the integration of quaternion rate equations in real-time flight simulation problems

High angular rates encountered in real-time flight simulation problems may require a more stable and accurate integration method than the classical methods normally used. A study was made to develop a general local linearization procedure of integrating dynamic system equations when using a digital computer in real-time. The procedure is specifically applied to the integration of the quaternion rate equations. For this application, results are compared to a classical second-order method. The local linearization approach is shown to have desirable stability characteristics and gives significant improvement in accuracy over the classical second-order integration methods.

Barker, L. E., Jr.↗

Privacy-preserving federated learning: Application to behind-the-meter solar photovoltaic generation forecasting

Here, the growing usage of decentralized renewable energy sources has made accurate estimation of their aggregated generation crucial for maintaining grid flexibility and reliability. However, the majority of distributed photovoltaic (PV) systems are behind-the-meter (BTM) and invisible to utilities, leading to three challenges in obtaining an accurate forecast of their aggregated output. Firstly, traditional centralized prediction algorithms used in previous studies may not be appropriate due to privacy concerns. There is therefore a need for decentralized forecasting methods, such as federated learning (FL), to protect privacy. Secondly, there has been no comparison between localized, centralized, and decentralized forecasting methods for BTM PV production, and the trade-off between prediction accuracy and privacy has not been explored. Lastly, the computational time of data-driven prediction algorithms has not been examined. This article presents a FL power forecasting method for PVs, which uses federated learning as a decentralized collaborative modeling approach to train a single model on data from multiple BTM sites. The machine learning network used to design this FL-based BTM PV forecasting model is a multi-layered perceptron, which ensures privacy and security of the data. Comparing the suggested FL forecasting model to non-private centralized and entirely private localized models revealed that it has a high level of accuracy, with an RMSE that is 18.17% lower than localized models and 9.9% higher than centralized models.

14 SOLAR ENERGY↗

Analysis of Fiber Clustering in Composite Materials Using High-Fidelity Multiscale Micromechanics

A new multiscale micromechanical approach is developed for the prediction of the behavior of fiber reinforced composites in presence of fiber clustering. The developed method is based on a coupled two-scale implementation of the High-Fidelity Generalized Method of Cells theory, wherein both the local and global scales are represented using this micromechanical method. Concentration tensors and effective constitutive equations are established on both scales and linked to establish the required coupling, thus providing the local fields throughout the composite as well as the global properties and effective nonlinear response. Two nondimensional parameters, in conjunction with actual composite micrographs, are used to characterize the clustering of fibers in the composite. Based on the predicted local fields, initial yield and damage envelopes are generated for various clustering parameters for a polymer matrix composite with both carbon and glass fibers. Nonlinear epoxy matrix behavior is also considered, with results in the form of effective nonlinear response curves, with varying fiber clustering and for two sets of nonlinear matrix parameters.

composite materials↗

A generalized scheme for three-dimensional hyperbolic grid generation

A hyperbolic grid generation scheme formulated from grid orthogonality and cell volume specification is generalized such that high quality three-dimensional grids can be obtained for a wide variety of geometries. The speed of the scheme is one to two orders of magnitude faster than typical elliptic grid generation methods. The robustness of the scheme is significantly enhanced by several new techniques. By using a spatially variable smoothing coefficient, orthogonality and smoothness of the grid are maintained around complex concave and convex surface topologies. A metric correction procedure is employed to guarantee that the grid marches out of a corner by bisecting the angles subtended by the neighboring points. Extra robustness at severe convex corners is achieved by special local treatments. Different extrapolation methods are used to provide smoothness at floating edges and axis regions. The versatility of the new hyperbolic grid generation scheme is demonstrated by three-dimensional grids generated for external components of the intergrated Space Shuttle vehicle and the SOFIA telescope.

Chan, William M.↗

Type Ia supernova magnitude step from the local dark matter environment

ABSTRACT Residuals in the Hubble diagram at optical wavelengths and host galaxy stellar mass are observed to correlate in type Ia supernovae (SNe Ia) (‘magnitude step’). Among possible progenitor channels for the associated explosions, those based on dark matter (DM) have attracted significant attention, including our recent proposal that ‘normal’ SNe Ia from bare detonations in sub-Chandrasekhar white dwarf stars are triggered by the passage of asteroid-mass primordial black holes (PBHs): the magnitude step could then originate from a brightness dependence on stellar properties, on DM properties, or both. Here, we present a method to estimate the local DM density and velocity dispersion of the environment of SN Ia progenitors. We find a luminosity step of $0.52\pm 0.11\,$ mag corresponding to bins of high versus low DM density in a sample of 222 low-redshift events from The Open Supernova Catalog. We investigate whether the magnitude step can be attributed to local DM properties alone, assuming asteroid-mass PBHs. Given the inverse correlation between SN Ia brightness and PBH mass, an intriguing explanation is a spatially inhomogeneous PBH mass function. If so, a strong mass segregation in the DM density-dependent PBH mass scale is needed to explain the magnitude step. While mass segregation is observed in dense clusters, it is unlikely to be realized on galactic scales. Therefore, if DM consists of asteroid-mass PBHs, the magnitude step is more likely to exist, and dominantly to be attributed to local stellar properties.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating Quantum Materials with Half-Polarized Diffraction and magnetic PDF analysis at the HB-2A Neutron Powder Diffractometer

Local magnetic ordering and anisotropy is often central to the emergent behavior and subsequent functional properties in quantum materials and beyond. Neutron powder diffraction provides a straightforward yet extremely powerful technique for quantitative measurements of microscopic magnetic properties. The HB-2A powder diffractometer located at the High Flux Isotope Reactor in ORNL is traditionally utilized for long-range magnetic structure determination. Recently these capabilities have been extended to include methods aimed at accessing local magnetism: Half- polarized neutron powder diffraction (pNPD) and magnetic pair distribution function (mPDF) analysis. These two distinct techniques are possible on HB-2A due to the versatility of the instrument’s reciprocal space coverage, resolution and novel ultra-low temperature multi-sample changers that operate down to dilution refrigerator temperatures. This provides unique capabilities not found on any powder diffraction instrument and is particularly well suited to investigations of magnetic quantum materials. The development and implementation of these techniques will be discussed with a series of science case examples ranging from geometric frustrated magnets to magnetic metal-organic frameworks. Data reduction and analysis tools will be presented that enable the extraction of the local site susceptibility tensor and local spin-spin correlations in real space. Finally, potential combinations of these techniques in the form of half-polarized magnetic pair distribution function (pmPDF) analysis will be considered. Looking forward, HB-2A is undergoing a detector upgrade that will be in the user program by 2026. This will offer an order of magnitude increase in count rates to further aid the development of these often low signal measurements and provide new scientific capabilities.

Neutron Scattering↗

Assessing sensitivities of climate model weighting to multiple methods, variables, and domains in the south-central United States

Given the increasing use of climate projections and multi-model ensemble weighting for a diverse array of applications, this project assesses the sensitivities of climate model weighting strategies and their resulting ensemble means to multiple components, such as the weighting schemes, climate variables, or spatial domains of interest. The purpose of this study is to assess the sensitivities associated with multi-model weighting strategies. The analysis makes use of global climate models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) and their statistically downscaled counterparts created with the localized constructed analogs (LOCA) method. This work focuses on historical and projected future mean precipitation and daily high temperatures of the south-central United States. Results suggest that the model weights and the corresponding weighted model means can be sensitive to the weighting strategy that is applied. For instance, when estimating model weights based on Louisiana precipitation, the weighted projections show a wetter and cooler south-central domain in the future compared to other weighting strategies. Alternatively, for example, when estimating model weights based on New Mexico temperature, the weighted projections show a drier and warmer south-central domain in the future. However, when considering the entire south-central domain in estimating the model weights, the weighted future projections show a compromise in the precipitation and temperature estimates. As for uncertainty, our matrix of results provided a more certain picture of future climate compared to the spread in the original model ensemble. If future impact assessments utilize weighting strategies, then our findings suggest that how the specific weighting strategy is used with climate projections may depend on the needs of an impact assessment or adaptation plan.

54 ENVIRONMENTAL SCIENCES↗

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Intersection-distribution-based remapping between arbitrary meshes for staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics

In this paper, we present a new intersection-distribution-based remapping method between arbitrary polygonal meshes for indirect staggered multi-material arbitrary Lagrangian-Eulerian hydrodynamics. All cell-centered material quantities are conservatively remapped using intersections between the Lagrangian (old, source) mesh and the rezoned (new, target) mesh. The new nodal masses are obtained by conservative distribution of all material masses in each new cell to the cell's corners and then collecting those corner masses at new nodes. This distribution is done using a local constrained optimization approach for each cell in the new mesh. In order to remap nodal momentum we first define cell-centered momentum for each cell in the old mesh, conservatively remap this to the new mesh and then conservatively distribute the new zonal momentum to each cell's bounding nodes, again using local constrained optimization. Our method also conserves total energy by applying a new nodal kinetic energy correction that relies on a process similar to that used for remapping nodal mass and momentum. Cell-centered kinetic energy is computed, conservatively remapped and then distributed to nodes. The discrepancy between this conservatively remapped and actual nodal kinetic energy is then conservatively distributed to the internal energies of the materials in the cells surrounding each node. Unlike conventional cell-based corrections of this type, this new nodal kinetic energy correction has not been observed to drive material internal energy negative in any of our testing. Unlike flux based remapping, our new intersection-distribution method can be applied to remapping between source and target meshes that are arbitrarily different, which provides superior flexibility in the rezoning strategy. Our method is accurate, essentially conservative and essentially bounds preserving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Local Structure Analysis and Modelling of Lignin‐Based Carbon Composites through the Hierarchical Decomposition of the Radial Distribution Function

Abstract Carbonized lignin has been proposed as a sustainable and domestic source of activated, amorphous, graphitic, and nanostructured carbon for many industrial applications as the structure can be tuned through processing conditions. However, the inherent variability of lignin and its complex physicochemical structure resulting from feedstock and pulping selection make the Process‐Structure‐Property‐Performance (PSPP) relationships hard to define. In this work, radial distribution functions (RDFs) from synchrotron X‐ray and neutron scattering of lignin‐based carbon composites (LBCCs) are investigated using the Hierarchical Decomposition of the Radial Distribution Function (HDRDF) modelling method to characterize the local atomic environment and develop quantitative PSPP relationships. PSPP relationships for LBCCs defined by this work include crystallite size dependence on lignin feedstock as well as increasing crystalline volume fraction, nanoscale composite density, and crystallite size with increasing reduction temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Renormalization and conformal invariance of non-local quantum electrodynamics

We study renormalization group flow in a non-local version of quantum electrodynamics (QED). We determine the regime in which the theory flows to a local theory in the infrared and study a possible UV completion of four-dimensional QED. In addition, we find that there exist non-local conformal theories with a one-dimensional conformal manifold and non-local deformations of QED in three dimensions that are exactly marginal. Along the way we develop methods for coupling non-local derivatives to external sources and discuss unitarity and conformal vs. scale invariance of these theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗