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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

The MHOST finite element program: 3-D inelastic analysis methods for hot section components. Volume 1: Theoretical manual

Formulations and algorithms implemented in the MHOST finite element program are discussed. The code uses a novel concept of the mixed iterative solution technique for the efficient 3-D computations of turbine engine hot section components. The general framework of variational formulation and solution algorithms are discussed which were derived from the mixed three field Hu-Washizu principle. This formulation enables the use of nodal interpolation for coordinates, displacements, strains, and stresses. Algorithmic description of the mixed iterative method includes variations for the quasi static, transient dynamic and buckling analyses. The global-local analysis procedure referred to as the subelement refinement is developed in the framework of the mixed iterative solution, of which the detail is presented. The numerically integrated isoparametric elements implemented in the framework is discussed. Methods to filter certain parts of strain and project the element discontinuous quantities to the nodes are developed for a family of linear elements. Integration algorithms are described for linear and nonlinear equations included in MHOST program.

Nakazawa, Shohei↗
NASA NTRS1984

Improved methods of vibration analysis of pretwisted, airfoil blades

Vibration analysis of pretwisted blades of asymmetric airfoil cross section is performed by using two mixed variational approaches. Numerical results obtained from these two methods are compared to those obtained from an improved finite difference method and also to those given by the ordinary finite difference method. The relative merits, convergence properties and accuracies of all four methods are studied and discussed. The effects of asymmetry and pretwist on natural frequencies and mode shapes are investigated. The improved finite difference method is shown to be far superior to the conventional finite difference method in several respects. Close lower bound solutions are provided by the improved finite difference method for untwisted blades with a relatively coarse mesh while the mixed methods have not indicated any specific bound.

Subrahmanyam, K. B.↗
DOE OSTI2017

Relation between the Variations in the Mg iiλ2798 Emission Line and 3000 Å Continuum

We investigate the relationship between the Mg II λ2798 emission line and the 3000 Å continuum variations using a sample of 68 intermediate-redshift (z∼0.65--1.50) broad-line quasars spanning a bolometric luminosity range of 44.49 erg s{sup −1}⩽logL{sub bol}⩽46.31 erg s{sup −1} (Eddington ratio from ∼0.026 to 0.862). This sample is constructed from SDSS-DR7Q and BOSS-DR12Q, each with at least two spectroscopic epochs in SDSS-I/II/III surveys. In addition, we adopt the following signal-to-noise ratio (S/N) selection criteria: (a) for Mg II and the 3000 Å continuum, S/N ⩾ 10; and (b) for narrow lines, S/N ⩾ 5. All our quasar spectra are recalibrated based on the assumption of constant narrow emission-line fluxes. In an analysis of spectrum-to-spectrum variations, we find a fairly close correlation (Spearman ρ=0.593) between the variations in broad Mg II and in the continuum. This is consistent with the idea that Mg II is varying in response to the continuum emission variations. Adopting the modified weighted least squares regression method, we statistically constrain the slopes (i.e., the responsivity α of the broad Mg II) between the variations in both components for the sources in different luminosity bins after eliminating intrinsic biases introduced by the rescaling process itself. It is shown that the responsivity is quite small (average α-bar ≈0.464) and anti-correlates with the quasar luminosity. Our results indicate that high signal-to-noise flux measurements are required to robustly detect the intrinsic variability and the time lag of the Mg II line.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Exploring precipitation pattern scaling methodologies and robustness among CMIP5 models

Pattern scaling is a well-established method for approximating modeled spatial distributions of changes in temperature by assuming a time-invariant pattern that scales with changes in global mean temperature. We compare two methods of pattern scaling for annual mean precipitation (regression and epoch difference) and evaluate which method is better in particular circumstances by quantifying their robustness to interpolation/extrapolation in time, inter-model variations, and inter-scenario variations. Both the regression and epoch-difference methods (the two most commonly used methods of pattern scaling) have good absolute performance in reconstructing the climate model output, measured as an area-weighted root mean square error. We decompose the precipitation response in the RCP8.5 scenario into a CO 2 portion and a non-CO 2 portion. Extrapolating RCP8.5 patterns to reconstruct precipitation change in the RCP2.6 scenario results in large errors due to violations of pattern scaling assumptions when this CO 2 -/non-CO 2 -forcing decomposition is applied. As a result, the methodologies discussed in this paper can help provide precipitation fields to be utilized in other models (including integrated assessment models or impacts assessment models) for a wide variety of scenarios of future climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1974

Finite element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is suited to isotropically hardening Prandtl-Reuss materials. The formulation is given in a manner which allows any conventional finite element program, for "small strain" elasticplastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. A unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures, and a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain.

Mcmeeking, R. M.↗
NASA NTRS1981

Effect of atmospheric scattering and surface reflection on upwelling solar radiation

A study is presented of the solar radiation transfer in the complete earth-atmosphere system, and numerical results are compared with satellite data obtained during the Earth Radiation Budget Experiment on Nimbus 6, in August, 1975. Emphasis is placed on the upwelling radiance distribution at the top of the atmosphere, assumed to be at 50 km. The numerical technique is based on the finite difference method, which includes azimuth and spectral variations for the entire solar wavelength range. Detailed solar properties, atmospheric physical properties, and optical properties are used. However, since the property descriptions are based on a trade-off between accuracy and computational realities, aerosol and cloud optical properties are treated with simple approximations. The radiative transfer model is in good agreement with the satellite radiance observations. The method provides a valuable tool in analyzing satellite- and ground-based radiation budget measurements and in designing instrumentation.

Suttles, J. T.↗
NASA NTRS1961

Aerodynamic Characteristics in Pitch of a Series of Cruciform-Wing Missiles with Canard Controls at a Mach Number of 2.01

An investigation has been conducted in the Langley 4- by 4-foot supersonic pressure tunnel at a Mach number of 2.01 to determine the static longitudinal stability and control characteristics of a series of missile configurations with canard controls at angles of attack up to about 28 deg. The missiles had cruciform wings and canard surfaces of delta plan form with 70 deg.swept leading edges. Five bodies having fineness ratios of 19.1, 17.7, 16.7, 15.7, and 14.8 were investigated. The results of the investigation indicated a large nonlinear variation of pitching moment with angle of attack for the body of largest fineness ratio that was progressively reduced by decreasing the fineness ratio until it was essentially eliminated for a body of fineness ratio 14.8. The increased linearity of the moment curve would make it possible to reduce the margin of stability so that, for a given canard size and deflection, a higher trim angle of attack might be obtained for the shortest missile than for the longest missile. The pitching-moment results indicated that methods of prediction which assumed linear variations with angle of attack for the wing-alone and wing-plus-interference characteristics were adequate for angles of attack up to about 12 deg. At higher angles of attack it was evident that the characteristics of these components were nonlinear and that more refined methods would be required for adequate prediction.

Spearman, M. Leroy↗
DOE OSTI2020

High-resolution analogue of time-domain phonon spectroscopy in the transmission electron microscope

Femtosecond photoexcitation of semiconducting materials leads to the generation of coherent acoustic phonons (CAPs), the behaviours of which are linked to intrinsic and engineered electronic, optical and structural properties. While often studied with pump-probe spectroscopic techniques, the influence of nanoscale structure and morphology on CAP dynamics can be challenging to resolve with these all-optical methods. Here, we used ultrafast electron microscopy (UEM) to resolve variations in CAP dynamics caused by differences in the degree of crystallinity in as-prepared and annealed GaAs lamellae. Following in situ femtosecond photoexcitation, we directly imaged the generation and propagation dynamics of hypersonic CAPs in a mostly amorphous and, following an in situ photothermal anneal, a mostly crystalline lamella. Subtle differences in both the initial hypersonic velocities and the asymptotic relaxation behaviours were resolved via construction of space-time contour plots from phonon wavefronts. Comparison to bulk sound velocities in crystalline and amorphous GaAs reveals the influence of the mixed amorphous-crystalline morphology on CAP dispersion behaviours. Further, an increase in the asymptotic velocity following annealing establishes the sensitivity of quantitative UEM imaging to both structural and compositional variations through differences in bonding and elasticity. Lastly, implications of extending the methods and results reported here to elucidating correlated electronic, optical and structural behaviours in semiconducting materials are discussed.

42 ENGINEERING↗
NASA NTRS2000

An Efficient Inverse Aerodynamic Design Method For Subsonic Flows

Computational Fluid Dynamics based design methods are maturing to the point that they are beginning to be used in the aircraft design process. Many design methods however have demonstrated deficiencies in the leading edge region of airfoil sections. The objective of the present research is to develop an efficient inverse design method which is valid in the leading edge region. The new design method is a streamline curvature method, and a new technique is presented for modeling the variation of the streamline curvature normal to the surface. The new design method allows the surface coordinates to move normal to the surface, and has been incorporated into the Constrained Direct Iterative Surface Curvature (CDISC) design method. The accuracy and efficiency of the design method is demonstrated using both two-dimensional and three-dimensional design cases.

Milholen, William E., II↗
NASA NTRS2008

Mapping Electrical Crosstalk in Pixelated Sensor Arrays

Electronic coupling effects such as Inter-Pixel Capacitance (IPC) affect the quantitative interpretation of image data from CMOS, hybrid visible and infrared imagers alike. Existing methods of characterizing IPC do not provide a map of the spatial variation of IPC over all pixels. We demonstrate a deterministic method that provides a direct quantitative map of the crosstalk across an imager. The approach requires only the ability to reset single pixels to an arbitrary voltage, different from the rest of the imager. No illumination source is required. Mapping IPC independently for each pixel is also made practical by the greater S/N ratio achievable for an electrical stimulus than for an optical stimulus, which is subject to both Poisson statistics and diffusion effects of photo-generated charge. The data we present illustrates a more complex picture of IPC in Teledyne HgCdTe and HyViSi focal plane arrays than is presently understood, including the presence of a newly discovered, long range IPC in the HyViSi FPA that extends tens of pixels in distance, likely stemming from extended field effects in the fully depleted substrate. The sensitivity of the measurement approach has been shown to be good enough to distinguish spatial structure in IPC of the order of 0.1%.

focal plane array↗
NASA NTRS1975

Finite-element formulations for problems of large elastic-plastic deformation

An Eulerian finite element formulation is presented for problems of large elastic-plastic flow. The method is based on Hill's variational principle for incremental deformations, and is ideally suited to isotropically hardening Prandtl-Reuss materials. Further, the formulation is given in a manner which allows any conventional finite element program, for 'small strain' elastic-plastic analysis, to be simply and rigorously adapted to problems involving arbitrary amounts of deformation and arbitrary levels of stress in comparison to plastic deformation moduli. The method is applied to a necking bifurcation analysis of a bar in plane-strain tension. The paper closes with a unified general formulation of finite element equations, both Lagrangian and Eulerian, for large deformations, with arbitrary choice of the conjugate stress and strain measures. Further, a discussion is given of other proposed formulations for elastic-plastic finite element analysis at large strain, and the inadequacies of some of these are commented upon.

Mcmeeking, R. M.↗
NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗
DOE OSTI2022

Variational Quantum Circuits to Prepare Low Energy Symmetry States

We explore how to build quantum circuits that compute the lowest energy state corresponding to a given Hamiltonian within a symmetry subspace by explicitly encoding it into the circuit. We create an explicit unitary and a variationally trained unitary that maps any vector output by ansatz A(α → ) from a defined subspace to a vector in the symmetry space. The parameters are trained varitionally to minimize the energy, thus keeping the output within the labelled symmetry value. The method was tested for a spin XXZ Hamiltonian using rotation and reflection symmetry and H 2 Hamiltonian within S z = 0 subspace using S 2 symmetry. We have found the variationally trained unitary gives good results with very low depth circuits and can thus be used to prepare symmetry states within near term quantum computers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

A high-order Lagrangian-decoupling method for the incompressible Navier-Stokes equations

A high-order Lagrangian-decoupling method is presented for the unsteady convection-diffusion and incompressible Navier-Stokes equations. The method is based upon: (1) Lagrangian variational forms that reduce the convection-diffusion equation to a symmetric initial value problem; (2) implicit high-order backward-differentiation finite-difference schemes for integration along characteristics; (3) finite element or spectral element spatial discretizations; and (4) mesh-invariance procedures and high-order explicit time-stepping schemes for deducing function values at convected space-time points. The method improves upon previous finite element characteristic methods through the systematic and efficient extension to high order accuracy, and the introduction of a simple structure-preserving characteristic-foot calculation procedure which is readily implemented on modern architectures. The new method is significantly more efficient than explicit-convection schemes for the Navier-Stokes equations due to the decoupling of the convection and Stokes operators and the attendant increase in temporal stability. Numerous numerical examples are given for the convection-diffusion and Navier-Stokes equations for the particular case of a spectral element spatial discretization.

Ho, Lee-Wing↗
DOE OSTI2026

Demystifying Piecewise and Localized Scatter Correction Methods

Multiplicative scatter is a common source of noise in near-infrared spectroscopy and other related instrumental techniques. A wide variety of methods are commonly used for the correction of multiplicative scatter. However, the majority of such methods assume that the parameters that describe the scatter are constant throughout the measured spectrum, which is often not the case. This work investigates a family of methods that perform scatter correction using local regions of neighboring wavelengths in order to better account for wavelength-dependent scattering. The methods in question are piecewise standard normal variate, localized standard normal variate, piecewise multiplicative scatter correction, and localized multiplicative scatter correction. This work describes the theoretical and algorithmic foundations of the family of local region-based scatter correction methods and compares their application and optimization at a qualitative and quantitative level using several datasets.

NIR spectroscopy↗
DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A matrix completion algorithm for efficient calculation of quantum and variational effects in chemical reactions

This work examines the viability of matrix completion methods as cost-effective alternatives to full nuclear Hessians for calculating quantum and variational effects in chemical reactions. The harmonic variety-based matrix completion (HVMC) algorithm, developed in a previous study (https://doi.org/10.1063/5.0018326), exploits the low-rank character of the polynomial expansion of potential energy to recover, using a small sample, vibrational frequencies (square roots of nuclear Hessian eigenvalues) constituting the reaction path. Furthermore, these frequencies are essential for calculating rate coefficients using variational transition state theory with multidimensional tunneling (VTST-MT). HVMC performance is examined for four SN2 reactions and five hydrogen transfer reactions, with each H-transfer reaction consisting of at least one vibrational mode strongly coupled to the reaction coordinate. HVMC is robust and captures zero-point energies, vibrational free energies, zero-curvature tunneling, and adiabatic ground state and free energy barriers as well as their positions on the reaction coordinate. For medium to large reactions involving H-transfer, with the exception of the most complex Ir catalysis system, less than 35% of total eigenvalue information is necessary for accurate recovery of key VTST-MT observables.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1985

Pitch Angle Distributions of Solar Energetic Particles and the Local Scattering Properties of the Interplanetary Magnetic Field

An approximate solution of the Fokker-Planck equation containing pitch angle scattering and adiabatic focusing is discussed. For modest focusing effects the omnidirectional density obeys an ordinary diffusion equation with a modified diffusion coefficient. The anisotropic part of the distribution function is properly normalized and split into an even and an odd part. The even part is determined by the ratio between the scattering mean free path and the focusing length and by the degree of polarization of the magnetic field fluctuations. The odd part is determined by the deviation of the pitch angle scattering from isotropic scattering. The method supplies a powerful tool to obtain the local characteristics of pitch angle scattering. It is insensitive to long lasting solar injections and to moderate radial variations of the mean free path. The method is applied to solar particle events observed on Helio-1 and -2.

Beeck, J.↗