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 73 records · Page 4

Efficient Calculation of NMR Shielding Constants Using Composite Method Approximations and Locally Dense Basis Sets

This paper presents a systematic study of applying composite method approximations with locally dense basis sets (LDBS) to efficiently calculate NMR shielding constants in small and medium-sized molecules. The pcSseg-n series of basis sets are shown to have similar accuracy to the pcS-n series when n ≥ 1 and can slightly reduce computational costs. We identify two different LDBS partition schemes that perform very effectively for density functional calculations. Here, we select a large subset of the recent NS372 database containing 290 H, C, N, and O shielding values evaluated by reference methods on 106 molecules to carefully assess methods of the high, medium, and low computational costs to make practical recommendations. Our assessment covers conventional electronic structure methods (density functional theory and wave function) with global basis calculations, as well as their use in one of the satisfactory LDBS approaches, and a range of composite approaches, also with and without LDBS. Altogether 99 methods are evaluated. On this basis, we recommend different methods to reach three different levels of accuracy and time requirements across the four nuclei considered.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local Mesh Refinement in the Space-Time CE/SE Method

A local mesh refinement procedure for the CE/SE method which does not use an iterative procedure in the treatments of grid-to-grid communications is described. It is shown that a refinement ratio higher than ten can be applied successfully across a single coarse grid/fine grid interface.

Chang, Sin-Chung↗

The active titration method for measuring local hydroxyl radical concentration

We are developing a method for measuring ambient OH by monitoring its rate of reaction with a chemical species. Our technique involves the local, instantaneous release of a mixture of saturated cyclic hydrocarbons (titrants) and perfluorocarbons (dispersants). These species must not normally be present in ambient air above the part per trillion concentration. We then track the mixture downwind using a real-time portable ECD tracer instrument. We collect air samples in canisters every few minutes for roughly one hour. We then return to the laboratory and analyze our air samples to determine the ratios of the titrant to dispersant concentrations. The trends in these ratios give us the ambient OH concentration from the relation: dlnR/dt = -k(OH). A successful measurement of OH requires that the trends in these ratios be measureable. We must not perturb ambient OH concentrations. The titrant to dispersant ratio must be spatially invariant. Finally, heterogeneous reactions of our titrant and dispersant species must be negligible relative to the titrant reaction with OH. We have conducted laboratory studies of our ability to measure the titrant to dispersant ratios as a function of concentration down to the few part per trillion concentration. We have subsequently used these results in a gaussian puff model to estimate our expected uncertainty in a field measurement of OH. Our results indicate that under a range of atmospheric conditions we expect to be able to measure OH with a sensitivity of 3x10(exp 5) cm(exp -3). In our most optimistic scenarios, we obtain a sensitivity of 1x10(exp 5) cm(exp -3). These sensitivity values reflect our anticipated ability to measure the ratio trends. However, because we are also using a rate constant to obtain our (OH) from this ratio trend, our accuracy cannot be better than that of the rate constant, which we expect to be about 20 percent.

Sprengnether, Michele↗

Pressure-gradient effects on hypersonic turbulent skin friction and boundary-layer profiles.

Local skin friction, total-temperature profiles, and pitot-pressure profiles were measured on the wall of a Mach-7.4 wind tunnel. The wall to adiabatic wall temperature ratio was varied from 0.3 to 0.5. Boundary-layer characteristics were compared with those predicted by a finite-difference method. Local skin friction was predicted to within 15%. Pressure-gradient effects on the temperature and Mach number distributions and the shape factor (displacement thickness/momentum thickness) were underpredicted, but the velocity distributions were closely predicted.

Hopkins, E. J.↗

The Localized Active Space Method with Unitary Selective Coupled Cluster

Here, we introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H 2 ) 2 , (H 2 ) 4 , and trans-butadiene, and the magnetic coupling constant for a bimetallic compound [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ . For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗

Application of viscous-inviscid interaction methods to transonic turbulent flows

Two different viscous-inviscid interaction schemes were developed for the analysis of steady, turbulent, transonic, separated flows over axisymmetric bodies. The viscous and inviscid solutions are coupled through the displacement concept using a transpiration velocity approach. In the semi-inverse interaction scheme, the viscous and inviscid equations are solved in an explicitly separate manner and the displacement thickness distribution is iteratively updated by a simple coupling algorithm. In the simultaneous interaction method, local solutions of viscous and inviscid equations are treated simultaneously, and the displacement thickness is treated as an unknown and is obtained as a part of the solution through a global iteration procedure. The inviscid flow region is described by a direct finite-difference solution of a velocity potential equation in conservative form. The potential equation is solved on a numerically generated mesh by an approximate factorization (AF2) scheme in the semi-inverse interaction method and by a successive line overrelaxation (SLOR) scheme in the simultaneous interaction method. The boundary-layer equations are used for the viscous flow region. The continuity and momentum equations are solved inversely in a coupled manner using a fully implicit finite-difference scheme.

Lee, D.↗

A global-local finite element method suitable for parallel computations

A global-local finite element analysis procedure is developed based on the fast convergent nature of FEM in displacement. A special scheme is used to utilize the global displacement solution as boundary conditions for local regions of interest. In the local region, a refined mesh is used for further stress analysis. This global-local procedure can be easily programmed in parallel on MIMD multi-processor computers for significant time savings. A Sequent Balance 21000 system is used for demonstrating the parallel programming.

Sun, C. T.↗

Coupling Finite Element and Meshless Local Petrov-Galerkin Methods for Two-Dimensional Potential Problems

A coupled finite element (FE) method and meshless local Petrov-Galerkin (MLPG) method for analyzing two-dimensional potential problems is presented in this paper. The analysis domain is subdivided into two regions, a finite element (FE) region and a meshless (MM) region. A single weighted residual form is written for the entire domain. Independent trial and test functions are assumed in the FE and MM regions. A transition region is created between the two regions. The transition region blends the trial and test functions of the FE and MM regions. The trial function blending is achieved using a technique similar to the 'Coons patch' method that is widely used in computer-aided geometric design. The test function blending is achieved by using either FE or MM test functions on the nodes in the transition element. The technique was evaluated by applying the coupled method to two potential problems governed by the Poisson equation. The coupled method passed all the patch test problems and gave accurate solutions for the problems studied.

Chen, T.↗

$\mathrm{BPS}$ coherent states and localization

We introduce coherent states averaged over a gauge group action to study correlators of half BPS states in $\mathcal{N}$= 4 SYM theory. The overlaps of these averaged coherent states are a generating function of correlators and can be written in terms of the Harish-Chandra-Itzykzon-Zuber (HCIZ) integral. We show that this formula immediately leads to a computation of the normalization of two point functions in terms of characters obtained originally in the work of Corley, Jevicki and Ramgoolam. We also find various generalizations for A n–1 quivers that follow directly from other solvable integrals over unitary groups. All of these can be computed using localization methods. When we promote the parameters of the generating function to collective coordinates, there is a dominant saddle that controls the effective action of these coherent states in the regime where they describe single AdS giant gravitons. We also discuss how to add open strings to this formulation. These will produce calculations that rely on correlators of matrix components of unitaries in the ensemble that is determined by the HCIZ integral to determine anomalous dimensions. We also discuss how sphere giants arise from Grassman integrals, how one gets a dominant saddle and how open strings are added in that case. The fact that there is a dominant saddle helps to understand how a 1/N expansion arises for open strings. We generalize the coherent state idea to study 1/4 and 1/8 BPS states as more general integrals over unitary groups.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electronic and thermal properties of the phase-change memory material, Ge 2 Sb 2 Te 5 , and results from spatially resolved transport calculations

Here, we report new insights into the electronic, structural, and transport (heat and charge) properties of the phase-change memory material amorphous Ge 2 Sb 2 Te 5 . Using realistic structural models of Konstantinou et al., (2019), we analyze the topology, electronic states, and lattice dynamics with density functional methods, including hybrid-functional calculations and machine-learned interatomic potentials. The Kohn–Sham orbitals near the Fermi level display a strong electron–phonon coupling, and exhibit large energy fluctuations at room temperature. The conduction tail states exhibit larger phonon-induced fluctuations than the valence tail states. To resolve transport at the atomic scale, we employ space-projected electronic conductivity and site-projected thermal conductivity methods. Local analysis of heat transport highlights the role of filamentary networks dominated by Te, with Sb and Ge making progressively smaller contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Enabling Multireference Calculations on Multimetallic Systems with Graphic Processing Units

Modeling multimetallic systems efficiently enables faster prediction of desirable chemical properties and the design of new materials. This work describes an initial implementation for performing multireference wave function method localized active-space self-consistent field (LASSCF) calculations through the use of multiple graphics processing units (GPUs) to accelerate time-to-solution. Density fitting is leveraged to reduce memory requirements, and we demonstrate the ability to fully utilize multi-GPU compute nodes. Performance improvements of 5–10x in total application runtime were observed in LASSCF calculations for multimetallic catalyst systems up to 1200 AOs and an active space of (22e,40o) using up to four NVIDIA A100 GPUs. Furthermore, written with performance portability in mind, a comparable performance is also observed in early runs on the Aurora exascale system using Intel Max Series GPUs.

Algorithms↗

Synthesis and Characterization of Uranium Complexes Supported by Substituted Aryldimethylsilylanilide Ligands

Here, we report the synthesis and characterization of substituted aryldimethylsilyldiisopropylanilide ligands and their respective bisamido complexes of U(III), (3,5-R 2 -PhMe 2 SiNDipp) 2 UI(dioxane) x (1, R = H, x= 0; 2, R = Me, x = 0; 3, R = t Bu, x = 1). We found that the steric bulk of the 3,5-R 2 -Ph ring affects the hapticity of the U–arene interaction. In the solid-state, 1 is a U–(η 6 -arene) complex, while 2 is a bis(U–(η 1 -arene)) complex. Theoretically calculated bond orders at PBE0 and PBE0-D3 levels of theory support these hapticity assignments. The 3,5- t Bu 2 -Ph rings of 3 are too bulky to interact with U and solid-state metrical parameters initially suggested a U–(η 1 -arene) interaction with one of the Dipp rings. However, bond order calculations show that this interaction is even weaker than in the previously reported ((PhMe 2 Si) 2 N) 3 U complex, leading to the conclusion that 3 is best described as a U–(η o -arene) complex. Molecular orbital analyses in conjunction with electron localization methods reveal that the U–(η 6 -arene) bonding in 1 is primarily electrostatic in nature. Some charge transfer takes place from the arene π orbitals to the U 6d/5f hybrid orbitals in addition to subtle δ-back-bonding. In 2 and 3, both π and δ interactions are substantially weaker, in agreement with the differences in the U–arene coordination modes. Surprisingly, attempts to generate less sterically bulky (3,5-R 2 -PhMe 2 SiNPh) 2 UI complexes results in disproportionation to homoleptic tetraamido (3,5-R 2 -PhMe 2 SiNPh) 4 U(IV) (4, R = H; 5, R = Me) complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Non-local detection of coherent Yu–Shiba–Rusinov quantum projections

Probing spatially confined quantum states from afar—a long-sought goal to minimize external interference—has been proposed to be feasible in condensed-matter systems through the coherent projection of the state. This can be achieved by engineering the eigenstates of the electron sea that surrounds the quantum state using cages built atom by atom, the so-called quantum corrals. However, the demonstration of the coherent nature of the projection and manipulation of its quantum composition are still important goals. Here we show this for the coherent projection of a Yu–Shiba–Rusinov quantum state that is induced by a magnetic impurity, using the eigenmodes of corrals on the surface of a superconductor. This enables us to manipulate the particle–hole composition of the projected state by tuning the corral eigenmodes through the Fermi energy. Our results demonstrate a controlled non-local method for the detection of magnet–superconductor hybrid quantum states.

quantum mechanics↗

Parallel-plate avalanche counters for heavy-ion beam tracking: History and mysteries

Despite being an old detector concept, position-sensitive parallel-plate avalanche counters (PPACs) remain widely used today for heavy-ion position and timing measurements. In modern rare isotope beam facilities and large-acceptance fragment separators, PPACs are used to characterize beam properties (diagnostics), facilitate beam delivery (tuning), and provide event-by-event beam tracking for particle identification (PID). Most popular particle localization methods in PPAC detectors are based on strip electrodes electrically connected to resistive-chain circuits or delay lines. More exotic systems include optical readouts based on recording electroluminescence light with high-granularity photodetector arrays or high-resolution resistive anode electrodes. PPACs with conventional resistive-chain or delay-line readouts achieve typical position and time resolutions of below 1 mm ( σ ) and around 150 ps ( σ ), respectively. In addition, delay-line systems are capable of counting rates above several hundred kHz for beam areas of a few millimeters square, and around 1 MHz for larger beam sizes. Resistive-chain readouts have limited rates of a few tens of kHz. A review of the operation principles and performance of PPAC detectors is presented in this paper. We will discuss decades-long experience building and operating PPACs developed at the National Superconducting Cyclotron Laboratory, and then refined at the Facility for Rare Isotope Beams (FRIB), mostly focusing on the delay-line readout technique (DLPPAC). We will also discuss problems that arise due to electric discharges at high rates and briefly describe ongoing developments at FRIB. Published by the American Physical Society 2024

Physics↗

MARBLES (Multi-scale Adaptively Refined Boltzmann LatticE Solver) [SWR-23-37]

MARBLES (Multi-scale Adaptively Refined Boltzmann LatticE Solver) is an open-source computational fluid dynamics package powered by the lattice Boltzmann equations and built on AMReX. In the lattice Boltzmann method, local collisions between meso-scale fictitious particles drive the governing equations which enables MARBLES to easily simulate flow around complex and/or moving geometry without the generation of a body-conforming mesh. Using AMReX data structures and operations ensures a high level of computational performance and parallel scaling on heterogenous architectures while also naturally supporting locally enhanced grid resolution and fidelity through automatic mesh refinement. New domains and problem definitions are easily specified through an input file with examples and guidance on all options and variables provided in the MARBLES documentation.

Henry de Frahan, Marc↗

Penalized ensemble Kalman filters for high dimensional non-linear systems

The ensemble Kalman filter (EnKF) is a data assimilation technique that uses an ensemble of models, updated with data, to track the time evolution of a usually non-linear system. It does so by using an empirical approximation to the well-known Kalman filter. However, its performance can suffer when the ensemble size is smaller than the state space, as is often necessary for computationally burdensome models. This scenario means that the empirical estimate of the state covariance is not full rank and possibly quite noisy. To solve this problem in this high dimensional regime, we propose a computationally fast and easy to implement algorithm called the penalized ensemble Kalman filter (PEnKF). Under certain conditions, it can be theoretically proven that the PEnKF will be accurate (the estimation error will converge to zero) despite having fewer ensemble members than state dimensions. Further, as contrasted to localization methods, the proposed approach learns the covariance structure associated with the dynamical system. These theoretical results are supported with simulations of several non-linear and high dimensional systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗