Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational algorithms”

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,585 records · Page 88

Estimating Eigenenergies from Quantum Dynamics: A Unified Noise-Resilient Measurement-Driven Approach

Ground state energy estimation in physical, chemical, and materials sciences is one of the most promising applications of quantum computing. In this work, we introduce a new hybrid approach that finds the eigenenergies by collecting real-time measurements and post-processing them using the machinery of dynamic mode decomposition (DMD). From the perspective of quantum dynamics, we establish that our approach can be formally understood as a stable variational method on the function space of observables available from a quantum many-body system. We also provide strong theoretical and numerical evidence that our method converges rapidly even in the presence of a large degree of perturbative noise, and show that the method bears an isomorphism to robust matrix factorization methods developed independently across various scientific communities. Our numerical benchmarks on spin and molecular systems demonstrate an accelerated convergence and a favorable resource reduction over state-of-the-art algorithms. The DMD-centric strategy can systematically mitigate noise and stands out as a leading hybrid quantum-classical eigensolver.

Shen, Yizhi↗

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objectives of this study is to develop an analysis to predict loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordinate and the finite element method in the axial coordinate. Application of Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗

Computer modeling of the sensitivity of a laser water vapor sensor to variations in temperature and air speed

Currently, there is disagreement among existing methods of determining atmospheric water vapor concentration at dew-points below -40 C. A major source of error is wall effects which result from the necessity of bringing samples into the instruments. All of these instruments also have response times on the order of seconds. NASA Langley is developing a water vapor sensor which utilizes the absorption of the infrared radiation produced by a diode laser to estimate water vapor concentration. The laser beam is directed through an aircraft window to a retroreflector located on an engine. The reflected beam is detected by an infrared detector located near the laser. To maximize signal to noise, derivative signals are analyzed. By measuring the 2f/DC signal and correcting for ambient temperature, atmospheric pressure and air speed (which results in a Doppler shifting of the laser beam), the water vapor concentration can be retrieved. Since this is an in situ measurement there are no wall effects and measurements can be made at a rate of more than 20 per second. This allows small spatial variations of water vapor to be studied. In order to study the sensitivity of the instrument to variations in temperature and air speed, a computer program which generated the 2f, 3f, 4f, DC and 2f/DC signals of the instrument as a function of temperature, pressure and air speed was written. This model was used to determine the effect of errors in measurement of the temperature and air speed on the measured water vapor concentration. Future studies will quantify the effect of pressure measurement errors, which are expected to be very small. As a result of these studied, a retrieval algorithm has been formulated, and will be applied to data taken during the PEM-West atmospheric science field mission. Spectroscopic studies of the water vapor line used by the instrument will be used to refine this algorithm. To prepare for these studies, several lasers have been studied to determine their output frequency range and power.

Tucker, George F.↗

Buckling of imperfect, anisotropic, ring-stiffened cylinders under combined loads

The objective of this study is to develop an anlysis to predict buckling loads of ring-stiffened anisotropic cylinders subject to axial compression, torsion, and internal pressure. This structure is modeled as a branched shell. A nonlinear axisymmetric prebuckling equilibrium state is assumed which is amenable to an exact solution within each branch. Axisymmetric geometric imperfections are included. Buckling displacements are represented by a Fourier series in the circumferential coordination and the finite-element method in the axial coordinate. Application of the Trefftz criterion to the second variation of the total potential energy leads to a nonlinear eigenvalue problem for the buckling load and mode. Results are presented for both unstiffened and ring-stiffened cylinders in the form of buckling interaction diagrams. Imperfections can cause an unexpected buckling mode in the ring web which would not occur for the perfect structure, and pressurization diminishes the benefit of adding rings to the unstiffened shell to increase the buckling load. The implementation of the analysis methodology into a structural sizing algorithm is discussed.

Ley, Robert P.↗

Mechanisms for Diurnal Variability of Global Tropical Rainfall Observed from TRMM

The behavior and various controls of diurnal variability in tropical-subtropical rainfall are investigated using Tropical Rainfall Measuring Mission (TRMM) precipitation measurements retrieved from: (1) TRMM Microwave Imager (TMI), (2) Precipitation Radar (PR), and (3) TMI/PR Combined, standard level 2 algorithms for the 1998 annual cycle. Results show that the diurnal variability characteristics of precipitation are consistent for all three algorithms, providing assurance that TRMM retrievals are providing consistent estimates of rainfall variability. As anticipated, most ocean areas exhibit more rainfall at night, while over most land areas rainfall peaks during daytime ,however, various important exceptions are found. The dominant feature of the oceanic diurnal cycle is a rainfall maximum in late-evening/early-morning (LE-EM) hours, while over land the dominant maximum occurs in the mid- to late-afternoon (MLA). In conjunction with these maxima are pronounced seasonal variations of the diurnal amplitudes. Amplitude analysis shows that the diurnal pattern and its seasonal evolution are closely related to the rainfall accumulation pattern and its seasonal evolution. In addition, the horizontal distribution of diurnal variability indicates that for oceanic rainfall there is a secondary MLA maximum, co-existing with the LE-EM maximum, at latitudes dominated by large scale convergence and deep convection. Analogously, there is a preponderance for an LE-EM maximum over land, co-existing with the stronger MLA maximum, although it is not evident that this secondary continental feature is closely associated with the large scale circulation. The ocean results clearly indicate that rainfall diurnal variability associated with large scale convection is an integral part of the atmospheric general circulation.

Yang, Song↗

Derived Land Surface Emissivity From Suomi NPP CrIS

Presented here is the land surface IR spectral emissivity retrieved from the Cross-track Infrared Sounder (CrIS) measurements. The CrIS is aboard the Suomi National Polar-orbiting Partnership (NPP) satellite launched on October 28, 2011. We describe the retrieval algorithm, demonstrate the surface emissivity retrieved with CrIS measurements, and inter-comparison with the Infrared Atmospheric Sounding Interferometer (IASI) emissivity. We also demonstrate that surface emissivity from satellite measurements can be used in assistance of monitoring global surface climate change, as a long-term measurement of IASI and CrIS will be provided by the series of EUMETSAT MetOp and US Joint Polar Satellite System (JPSS) satellites. Monthly mean surface properties are produced using last 5-year IASI measurements. A temporal variation indicates seasonal diversity and El Nino/La Nina effects not only shown on the water but also on the land. Surface spectral emissivity and skin temperature from current and future operational satellites can be utilized as a means of long-term monitoring of the Earth's environment. CrIS spectral emissivity are retrieved and compared with IASI. The difference is small and could be within expected retrieval error; however it is under investigation.

Zhou, Daniel K.↗

Modern Monte Carlo methods for efficient uncertainty quantification and propagation: A survey

Uncertainty quantification (UQ) includes the characterization, integration, and propagation of uncertainties that result from stochastic variations and a lack of knowledge or data in the natural world. Monte Carlo (MC) method is a sampling-based approach that has widely used for quantification and propagation of uncertainties. However, the standard MC method is often time-consuming if the simulation-based model is computationally intensive. This article gives an overview of modern MC methods to address the existing challenges of the standard MC in the context of UQ. Specifically, multilevel Monte Carlo (MLMC) extending the concept of control variates achieves a significant reduction of the computational cost by performing most evaluations with low accuracy and corresponding low cost, and relatively few evaluations at high accuracy and corresponding high cost. Multifidelity Monte Carlo (MFMC) accelerates the convergence of standard Monte Carlo by generalizing the control variates with different models having varying fidelities and varying computational costs. Multimodel Monte Carlo method (MMMC), having a different setting of MLMC and MFMC, aims to address the issue of UQ and propagation when data for characterizing probability distributions are limited. Multimodel inference combined with importance sampling is proposed for quantifying and efficiently propagating the uncertainties resulting from small data sets. All of these three modern MC methods achieve a significant improvement of computational efficiency for probabilistic UQ, particularly uncertainty propagation. Finally, an algorithm summary and the corresponding code implementation are provided for each of the modern MC methods. The extension and application of these methods are discussed in detail.

97 MATHEMATICS AND COMPUTING↗

Toward In Situ Synchrotron Mapping of Crystal Selection Processes during Crystal Growth

In this work, we present an automated and rapid method for non-destructive mapping of crystal grains in a rod-shaped sample. The approach was designed for application to in situ float-zone crystal growth experiments at an x-ray synchrotron source, but could be useful in other applications. The methods have been tested on a TiO 2 boule grown in an optical float zone furnace. The approach applies a statistical filter to polycrystalline diffraction patterns on 2D detectors to rapidly determine the degree of powder-quality 1of the signal. When larger crystals emerge in the growth, their position, size and shape can be tracked using an automated blob-tracking algorithm that follows individual Bragg peaks as a function of position in a grid-scan, even when multiple crystals are contributing spots to diffraction images. This method is found to be robust as the same crystal shape can be independently reconstructed using different sets of Bragg reflections. Image segmentation methods are then used to map out the polycrystalline grains. We also note that other information about crystal quality, such as mosaicity or strain state, may be inferred and mapped from the intensity variation of the Bragg peaks at different locations within the sample.

36 MATERIALS SCIENCE↗

Implementation and Synthesis of Math Library Functions

Achieving speed and accuracy for math library functions like exp, sin, and log is difficult. This is because low-level implementation languages like C do not help math library developers catch mathematical errors, build implementations incrementally, or separate high-level and low-level decision making. This ultimately puts development of such functions out of reach for all but the most experienced experts. To address this, we introduce MegaLibm, a domain-specific language for implementing, testing, and tuning math library implementations. MegaLibm is safe, modular, and tunable. Implementations in MegaLibm can automatically detect mathematical mistakes like sign flips via semantic wellformedness checks, and components like range reductions can be implemented in a modular, composable way, simplifying implementations. Once the high-level algorithm is done, tuning parameters like working precisions and evaluation schemes can be adjusted through orthogonal tuning parameters to achieve the desired speed and accuracy. MegaLibm also enables math library developers to work interactively, compiling, testing, and tuning their implementations and invoking tools like Sollya and type-directed synthesis to complete components and synthesize entire implementations. MegaLibm can express 8 state-of-the-art math library implementations with comparable speed and accuracy to the original C code, and can synthesize 5 variations and 3 from-scratch implementations with minimal guidance.

97 MATHEMATICS AND COMPUTING↗

The Impact of Constrained Data Assimilation on the Forecasts of Three Convection Systems During the ARM MC3E Field Campaign

A constrained data assimilation (CDA) system based on the ensemble variational (EnVar) method and physical constraints of mass and water conservations is evaluated through three convective cases during the Midlatitude Continental Convective Clouds Experiment (MC3E) of the Atmospheric Radiation Measurement (ARM) program. Compared to the original data assimilation (ODA), the CDA is shown to perform better in the forecasted state variables and simulated precipitation. The CDA is also shown to greatly mitigate the loss of forecast skills in observation denial experiments when radar radial winds are withheld in the assimilation. Modifications to the algorithm and sensitivities of the CDA to the calculation of the time tendencies in the constraints are described.

54 ENVIRONMENTAL SCIENCES↗

Algorithm for fixed-range optimal trajectories

An algorithm for synthesizing optimal aircraft trajectories for specified range was developed and implemented in a computer program written in FORTRAN IV. The algorithm, its computer implementation, and a set of example optimum trajectories for the Boeing 727-100 aircraft are described. The algorithm optimizes trajectories with respect to a cost function that is the weighted sum of fuel cost and time cost. The optimum trajectory consists at most of a three segments: climb, cruise, and descent. The climb and descent profiles are generated by integrating a simplified set of kinematic and dynamic equations wherein the total energy of the aircraft is the independent or time like variable. At each energy level the optimum airspeeds and thrust settings are obtained as the values that minimize the variational Hamiltonian. Although the emphasis is on an off-line, open-loop computation, eventually the most important application will be in an on-board flight management system.

Lee, H. Q.↗

Modeling of composite beams and plates for static and dynamic analysis

The main purpose of this research was to develop a rigorous theory and corresponding computational algorithms for through-the-thickness analysis of composite plates. This type of analysis is needed in order to find the elastic stiffness constants for a plate and to post-process the resulting plate solution in order to find approximate three-dimensional displacement, strain, and stress distributions throughout the plate. This also requires the development of finite deformation plate equations which are compatible with the through-the-thickness analyses. After about one year's work, we settled on the variational-asymptotical method (VAM) as a suitable framework in which to solve these types of problems. VAM was applied to laminated plates with constant thickness in the work of Atilgan and Hodges. The corresponding geometrically nonlinear global deformation analysis of plates was developed by Hodges, Atilgan, and Danielson. A different application of VAM, along with numerical results, was obtained by Hodges, Lee, and Atilgan. An expanded version of this last paper was submitted for publication in the AIAA Journal.

Hodges, Dewey H.↗

An efficient design sensitivity analysis of eigenvectors

Subspace iteration has been a major advance in solving large eigen problems when only a subset of eigen-pairs is required. The essence of this method is a transformation from displacement coordinates of an n-th order eigensystem to generalized coordinates of a smaller q-th order. The eigenvalue problem is then solved in the reduced space. The method was first developed by Clint and Jennings for real symmetric systems and was then called 'simultaneous iteration'. The success of the method prompted further research along this line and there have been many improved algorithms developed. This approach has been widely used by structural engineers for extracting the most useful natural frequencies and mode shapes of large-scale dynamic systems. This paper exploits into a new direction which is in the form of iterative process for simultaneously calculating eigenvector derivatives of many eigenvectors with respect to multi-variables. The method fully uses all the available information from preceding eigenvalue solution and, thus, effectively economizes computational efforts. It iterates through two equations derived from the first variation of the two fundamental equations used in subspace iteration method. There is no expensive large matrix decomposition required and the process converges to acceptable solution in a finite number of iterations. Therefore, the procedure increases its efficiency superiority over the others as the system size or the number of interested eigenvectors become larger and larger.

Ting, T.↗

A Multi-Instrument Approach for Characterizing the Vertical Structure of Aerosol Properties: Case Studies in the Pacific Basin Troposphere

During February/March 1994, a series of aircraft-based aerosol measurements were carried out in the Pacific Basin troposphere using a differential absorption lidar system deployed by NASA Langley, and optical spectrometer probes and a wire-impactor system operated by NASA Ames. A modified Klett inversion algorithm was applied to extract altitude profiles of aerosol backscattering from the IR lidar signal. The algorithm that we have designed for this purpose utilizes the in situ aerosol measurements to normalize the lidar profile at the aircraft altitude and to supply the lidar ratio as a function of height. The lidar-derived aerosol backscattering coefficients were then compared to the backscattering coefficients calculated from the in situ measurements. During several local aircraft descents, we found good agreement between the remote lidar and in situ results for the absolute value of the aerosol backscattering coefficient and its altitude variation only when we allowed for several layers with different aerosol refractive indices. The agreement validates our lidar calibration method and provides an indication of the variation in aerosol refractive index as a function of altitude. Two of the three case studies performed in this paper reveal layers of anthropogenic aerosols transported long distances into the Pacific Basin troposphere. A third case implies the existence of a layer of dustlike aerosol particles in the lower troposphere, most likely of Asian origin.

Redemann, J.↗

Large Array Channel Capacity in the Presence of Interference

We develop a model for a large array ground receiver system for use in deep-space communications, and analyze the resulting array channel capacity. The model includes effects of array geometry, time-dependent spacecraft orbital trajectory, point and extended interference sources, and elevation-dependent noise and tropospheric channel variations. Channel capacity is expressed as the ratio of determinants of covariance matrices characterizing source, interference, and additive noise, and then reduced to a simpler quadratic form more amenable to analysis and numerical computation. This formulation facilitates inclusion of array and channel characteristics into the model, as well as comparison of optimal, suboptimal, and equivalent single antenna configurations on achievable throughput. Realistic examples of ground array channel capacity calculations are presented, demonstrating the impact of array geometry, planetary interference sources, and array combining algorithm design upon the achievable data throughput.

Vilnrotter, V.↗

A Comparison of the SOCIT and DebriSat Experiments

This paper explores the differences between, and shares the lessons learned from, two hypervelocity impact experiments critical to the update of orbital debris environment models. The procedures and processes of the fourth Satellite Orbital Debris Characterization Impact Test (SOCIT) were analyzed and related to the ongoing DebriSat experiment. SOCIT was the first hypervelocity impact test designed specifically for satellites in Low Earth Orbit (LEO). It targeted a 1960's U.S. Navy satellite, from which data was obtained to update pre-existing NASA and DOD breakup models. DebriSat is a comprehensive update to these satellite breakup models- necessary since the material composition and design of satellites have evolved from the time of SOCIT. Specifically, DebriSat utilized carbon fiber, a composite not commonly used in satellites during the construction of the US Navy Transit satellite used in SOCIT. Although DebriSat is an ongoing activity, multiple points of difference are drawn between the two projects. Significantly, the hypervelocity tests were conducted with two distinct satellite models and test configurations, including projectile and chamber layout. While both hypervelocity tests utilized soft catch systems to minimize fragment damage to its post-impact shape, SOCIT only covered 65% of the projected area surrounding the satellite, whereas, DebriSat was completely surrounded cross-range and downrange by the foam panels to more completely collect fragments. Furthermore, utilizing lessons learned from SOCIT, DebriSat's post-impact processing varies in methodology (i.e., fragment collection, measurement, and characterization). For example, fragment sizes were manually determined during the SOCIT experiment, while DebriSat utilizes automated imaging systems for measuring fragments, maximizing repeatability while minimizing the potential for human error. In addition to exploring these variations in methodologies and processes, this paper also presents the challenges DebriSat has encountered thus far and how they were addressed. Accomplishing DebriSat's goal of collecting 90% of the debris, which constitutes well over 100,000 fragments, required addressing many challenges stemming from the very large number of fragments. One of these challenges arose in identifying the foam-embedded fragments. DebriSat addressed this by X-raying all of the panels once the loose debris were removed, and applying a detection algorithm developed in-house to automate the embedded fragment identification process. It is easy to see how the amount of data being compiled would be outstanding. Creating an efficient way to catalog each fragment, as well as archiving the data for reproducibility also posed a great challenge for DebriSat. Barcodes to label each fragment were introduced with the foresight that once the characterization process began, the datasheet for each fragment would have to be accessed again quickly and efficiently. The DebriSat experiment has benefited significantly by leveraging lessons learned from the SOCIT experiment along with the technological advancements that have occurred during the time between the experiments. The two experiments represent two ages of satellite technology and, together, demonstrate the continuous efforts to improve the experimental techniques for fragmentation debris characterization.

Ausay, Erick↗

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Three-dimensional self-adaptive grid method for complex flows

A self-adaptive grid procedure for efficient computation of three-dimensional complex flow fields is described. The method is based on variational principles to minimize the energy of a spring system analogy which redistributes the grid points. Grid control parameters are determined by specifying maximum and minimum grid spacing. Multidirectional adaptation is achieved by splitting the procedure into a sequence of successive applications of a unidirectional adaptation. One-sided, two-directional constraints for orthogonality and smoothness are used to enhance the efficiency of the method. Feasibility of the scheme is demonstrated by application to a multinozzle, afterbody, plume flow field. Application of the algorithm for initial grid generation is illustrated by constructing a three-dimensional grid about a bump-like geometry.

Djomehri, M. Jahed↗