Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational 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 1,225 records · Page 68

A parallel, distributed memory implementation of the adaptive sampling configuration interaction method

The many-body simulation of quantum systems is an active field of research that involves several different methods targeting various computing platforms. Many methods commonly employed, particularly coupled cluster methods, have been adapted to leverage the latest advances in modern high-performance computing. Selected configuration interaction (sCI) methods have seen extensive usage and development in recent years. However, the development of sCI methods targeting massively parallel resources has been explored only in a few research works. Here, we present a parallel, distributed memory implementation of the adaptive sampling configuration interaction approach (ASCI) for sCI. In particular, we will address the key concerns pertaining to the parallelization of the determinant search and selection, Hamiltonian formation, and the variational eigenvalue calculation for the ASCI method. Load balancing in the search step is achieved through the application of memory-efficient determinant constraints originally developed for the ASCI-PT2 method. The presented benchmarks demonstrate near optimal speedup for ASCI calculations of Cr 2 (24e, 30o) with 10 6 , 10 7 , and 3 × 10 8 variational determinants on up to 16 384 CPUs. Importantly, to the best of the authors’ knowledge, this is the largest variational ASCI calculation to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Easy_PERT: a Python tool for writing PERT cards and parsing PERT card results [Slides]

This presentation begins by providing an overview of the PERT card. The PERT card uses differential operator method to compute first- and second-order tally variations due to density, composition, and reaction cross-sections. It is possible to have multiple PERT cards in one MCNP input deck to study tally variations for several sets of nuclides, reactions, and energy ranges. Furthermore, the METHOD option tells MCNP to calculate either the perturbed tally (METHOD=-1, -2, -3) or the change in the unperturbed tally (METHOD=1, 2, 3). In summation, a powerful use-case for the MCNP code PERT card is that it facilitates calculating tally sensitivities to nuclear data. Writing PERT card entries and parsing output MCTAL files is tedious and error prone. however, Easy_PERT makes use of existing tools (Faust and MCNPTools) to handle writing PERT card entries and parsing the output MCTAL files. The PERT card is early in the development process and planned upcoming capabilities include calculating sensitivities and combining MCTAL files from separate runs into one JSON file.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An efficient, direct finite difference method for computing sound propagation in arbitrarily shaped two-dimensional and axisymmetric ducts without flow

An efficient, direct finite difference method is presented for computing sound propagation in non-stepped two-dimensional and axisymmetric ducts of arbitrarily varying cross section without mean flow. The method is not restricted by axial variation of acoustic impedance of the duct wall linings. The non-uniform two-dimensional or axisymmetric duct is conformally mapped numerically into a rectangular or cylindrical computational domain using a new procedure based on a method of fast direct solution of the Cauchy-Riemann equations. The resulting Helmholtz equation in the computational domain is separable. The solution to the governing equation and boundary conditions is expressed as a linear combination of fundamental solutions. The fundamental solutions are computed only once for each duct shape by means of the fast direct cyclic reduction method for the discrete solution of separable elliptic equations. Numerical results for several examples are presented to show the applicability and efficiency of the method.

Chakravarthy, S.↗

A meshless Galerkin method for non-local diffusion using localized kernel bases

Here, we introduce a meshless method for solving both continuous and discrete variational formulations of a volume constrained, nonlocal diffusion problem. We use the discrete solution to approximate the continuous solution. Our method is nonconforming and uses a localized Lagrange basis that is constructed out of radial basis functions. By verifying that certain inf-sup conditions hold, we demonstrate that both the continuous and discrete problems are well-posed, and also present numerical and theoretical results for the convergence behavior of the method. The stiffness matrix is assembled by a special quadrature routine unique to the localized basis. Combining the quadrature method with the localized basis produces a well-conditioned, symmetric matrix. This then is used to find the discretized solution.

97 MATHEMATICS AND COMPUTING↗

Simulation of a Detector of Visible and Near-IR Electromagnetic Radiation Based on Artificial Diamond

A method for calculation of the I–V characteristic is proposed for a detector of visible and near-IR electromagnetic radiation based on artificial diamond with allowance for vertical flow of electric current. The method can be used when the free-path length of carriers in diamond is less than the scale of variations in the concentration of carriers. The method is used to determine the following detector characteristics: the dependence of current on voltage, distribution of carrier concentration, and electric potential for a particular variant of a photodetector based on polycrystalline diamond film consisting of nanosized single crystals doped with boron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

On Properties of Adjoint Systems for Evolutionary PDEs

We investigate the geometric structure of adjoint systems associated with evolutionary partial differential equations at the fully continuous, semi-discrete, and fully discrete levels and the relations between these levels. We show that the adjoint system associated with an evolutionary partial differential equation has an infinite-dimensional Hamiltonian structure, which is useful for connecting the fully continuous, semi-discrete, and fully discrete levels. We subsequently address the question of discretize-then-optimize versus optimize-then-discrete for both semi-discretization and time integration, by characterizing the commutativity of discretize-then-optimize methods versus optimize-then-discretize methods uniquely in terms of an adjoint-variational quadratic conservation law. For Galerkin semi-discretizations and one-step time integration methods in particular, we explicitly construct these commuting methods by using structure-preserving discretization techniques.

97 MATHEMATICS AND COMPUTING↗

Geostationary satellite position determination for common-view two-way time transfer measurements

In common-view two-way time transfer, each earth station receives an unwanted return signal from its own transmission as well as the desired signal from the other earth station. National Institute of Standards and Technology (NIST), the National Research Council (NRC), and the U.S. Naval Observatory (USNO) have been cooperating in a three-corner common-view two-way time transfer experiment. Some systematic effects are known to depend on the position of the satellite (Sagnac effect and the cross-correlation pulling of the pseudo-random codes). A method is presented for deriving accurate satellite ranges from each of three stations doing common-view two-way satellite time transfer measurements, when one (and only one) station also takes ranging measurements on its 'unwanted return signal' for a brief period. The method is applied to determine the variations in position of the satellite used over the course of the NIST/NRC/USNO SBS-3 experiment, with ranging data taken at NRC, where no additional hardware was required to automate the process. The fit and extrapolation which are employed in this method have an estimated precision of 2 m. If the delays of SBS-3 satellite KU band transponder and earth station equipment were measured accurately as well as the tropospheric refractions were well modeled and corrected, we would expect a ranging accuracy of 2.5 m and satellite positioning accuracy would be 200 m (latitude) 50 m (longitude) and 20 m (height above ellipsoid).

Zhuang, Qixiang↗

Numerical methods for the inverse problem of density functional theory

Here, the inverse problem of Kohn–Sham density functional theory (DFT) is often solved in an effort to benchmark and design approximate exchange-correlation potentials. The forward and inverse problems of DFT rely on the same equations but the numerical methods for solving each problem are substantially different. We examine both problems in this tutorial with a special emphasis on the algorithms and error analysis needed for solving the inverse problem. Two inversion methods based on partial differential equation constrained optimization and constrained variational ideas are introduced. We compare and contrast several different inversion methods applied to one-dimensional finite and periodic model systems.

97 MATHEMATICS AND COMPUTING↗

Vector potential methods

Vector potential and related methods, for the simulation of both inviscid and viscous flows over aerodynamic configurations, are briefly reviewed. The advantages and disadvantages of several formulations are discussed and alternate strategies are recommended. Scalar potential, modified potential, alternate formulations of Euler equations, least-squares formulation, variational principles, iterative techniques and related methods, and viscous flow simulation are discussed.

Hafez, M.↗

Combining variational autoencoders and physical bias for improved microscopy data analysis *

Electron and scanning probe microscopy produce vast amounts of data in the form of images or hyperspectral data, such as electron energy loss spectroscopy or 4D scanning transmission electron microscope, that contain information on a wide range of structural, physical, and chemical properties of materials. To extract valuable insights from these data, it is crucial to identify physically separate regions in the data, such as phases, ferroic variants, and boundaries between them. In order to derive an easily interpretable feature analysis, combining with well-defined boundaries in a principled and unsupervised manner, here we present a physics augmented machine learning method which combines the capability of variational autoencoders to disentangle factors of variability within the data and the physics driven loss function that seeks to minimize the total length of the discontinuities in images corresponding to latent representations. Our method is applied to various materials, including NiO-LSMO, BiFeO 3 , and graphene. The results demonstrate the effectiveness of our approach in extracting meaningful information from large volumes of imaging data. The customized codes of the required functions and classes to develop phyVAE is available at https://github.com/arpanbiswas52/phy-VAE.

97 MATHEMATICS AND COMPUTING↗

Results of the International Energy Agency Round Robin on Fast Pyrolysis Bio-oil Production

An international round robin study of the production of fast pyrolysis bio-oil was undertaken. Fifteen institutions in six countries contributed. Three biomass samples were distributed to the laboratories for processing in fast pyrolysis reactors. Samples of the bio-oil produced were transported to a central analytical laboratory for analysis. The round robin was focused on validating the pyrolysis community understanding of production of fast pyrolysis bio-oil by providing a common feedstock for bio-oil preparation. The round robin included: •distribution of 3 feedstock samples from a common source to each participating laboratory; •preparation of fast pyrolysis bio-oil in each laboratory with the 3 feedstocks provided; •return of the 3 bio-oil products (minimum 500 ml) with operational description to a central analytical laboratory for bio-oil property determination. The analyses of interest were: density, viscosity, dissolved water, filterable solids, CHN, S, trace element analysis, ash, total acid number, pyrolytic lignin, and accelerated aging of bio-oil. In addition, an effort was made to compare the bio-oil components to the products of analytical pyrolysis through GC/MS analysis. The results showed that clear differences can occur in fast pyrolysis bio-oil properties by applying different reactor technologies or configurations. The comparison to analytical pyrolysis method suggested that Py-GC/MS could serve as a rapid screening method for bio-oil composition when produced in fluid-bed reactors. Furthermore, hot vapor filtration generally resulted in the most favorable bio-oil product, with respect to water, solids, viscosity, and total acid number. These results can be helpful in understanding the variation in bio-oil production methods and their effects on bio-oil product composition.

09 BIOMASS FUELS↗

A theory of wave scatter from an inhomogeneous medium with a slightly rough boundary and its application to sea ice

An analytical theory of electromagnetic wave scattering from an inhomogeneous medium with a slightly rough boundary surface is formulated. The inhomogeneity in the medium is assumed to vary continuously in the vertical direction and to have a small random variation in the horizontal direction. The medium is assumed to consist of two layers. Maxwell's equations are solved by using the small perturbation method together with Fourier transform technique. The resulting differential equations are solved by using WKB and variation of parameter methods. Field amplitudes in each medium are determined by taking boundary conditions into account. The expressions for first order polarized radar backscatter cross-section are obtained. An attempt is made to apply the developed theory to compute sea ice scatter. Numerical calculations are performed for polarized radar backscatter cross-section at two frequencies, 13.3 GHz and 400 MHz. It is shown that WKB method is applicable at both of these frequencies. Theoretical results are compared with the experimental results obtained from NASA Earth Resources Program mission 126. Theoretical results and experimental results are in good agreement.

Parashar, S. K.↗

Scaling up the Single Transducer Thickness-Independent Ultrasonic Imaging Method for Accurate Characterization of Microstructural Gradients in Monolithic and Composite Tubular Structures

Ultrasonic velocity/time-of-flight imaging that uses back surface reflections to gauge volumetric material quality is highly suited for quantitative characterization of microstructural gradients including those due to pore fraction, density, fiber fraction, and chemical composition variations. However, a weakness of conventional pulse-echo ultrasonic velocity/time-of-flight imaging is that the image shows the effects of thickness as well as microstructural variations unless the part is uniformly thick. This limits this imaging method's usefulness in practical applications. Prior studies have described a pulse-echo time-of-flight-based ultrasonic imaging method that requires using a single transducer in combination with a reflector plate placed behind samples that eliminates the effect of thickness variation in the image. In those studies, this method was successful at isolating ultrasonic variations due to material microstructure in plate-like samples of silicon nitride, metal matrix composite, and polymer matrix composite. In this study, the method is engineered for inspection of more complex-shaped structures-those having (hollow) tubular/curved geometry. The experimental inspection technique and results are described as applied to (1) monolithic mullite ceramic and polymer matrix composite 'proof-of-concept' tubular structures that contain machined patches of various depths and (2) as-manufactured monolithic silicon nitride ceramic and silicon carbide/silicon carbide composite tubular structures that might be used in 'real world' applications.

Roth, Don J.↗

A Novel Ultrasonic Method for Characterizing Microstructural Gradients in Tubular Structures

Ultrasonic velocity and time-of-flight (TOF) imaging that uses back surface reflections to gage volumetric material quality is highly suited for quantitative characterization of microstructural gradients including those due to pore fraction, density, fiber fraction, and chemical composition variations. However, a weakness of conventional pulse echo ultrasonic velocity and TOF imaging is that the image shows the effects of thickness as well as microstructural variations, unless the part is uniformly thick. This limits this imaging method's usefulness in practical applications. Prior studies have described a pulse echo TOF based ultrasonic imaging method that requires using a single transducer in combination with a reflector plate placed behind samples which, eliminates the effect of thickness variation in the image. In those studies, this method was successful at isolating ultrasonic variations due to material in plate like samples of silicon nitride, metal matrix composite, and polymer matrix composite. In this study, the method is engineered for inspection of more complex shaped structures- those having (hollow) tubular or curved geometry. The experimental inspection technique and results are described as applied to a polymer matrix composite "proof of concept" tube that contains machined patches of various depths and an as manufactured monolithic silicon nitride ceramic tube that might be used in "real world" applications.

Roth, Don J.↗

In-Line Optical Transmission Imaging of Decals for Quality Control - Task 3

Quality monitoring is a critical aspect for manufacturing systems. Ideally the monitoring would be done in-line, be non-contact, non-destructive, and fast. This would enable reduced scrap and higher throughput. This poster presents an optical transmission method for evaluating and mapping coatings. With the method shown in the poster we can visualize optical variations on the macro and micro scales. This allows us to see the overall trend in loading in both the cross web and down web directions. Furthermore, we can visualize defects such dewetting spots, streaks, clumps, and pinholes where there is a lack of coating. The optical transmission signal has been found to be proportional to the IrOx loading signal using XRF measurements. Therefore, an optical transmission setup can be installed in-line and allow for a fast, non-contact method for mapping loading variations and defects.

coating uniformity↗

Surface polarimetry of the Milky Way

Two methods of determining the state of polarization of the integrated light from the Milky Way are discussed. In the symmetry method, it is assumed that the polarized foreground zodiacal light possesses north/south symmetry relative to the ecliptic. In the method of two epochs, the above assumption is not needed, but instead time variations of the zodiacal light are assumed to be insignificant; the results obtained by this method, although compatible with the values found by the symmetry method, do however show that time variations may have influenced this result. Although the evidence for a definite polarization of the Milky Way is persuasive it cannot yet be considered conclusive.

Wolstencroft, R. D.↗