Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 451 records · Page 25

Sensitivity of actively damped structures to imperfections and modeling errors

The sensitivity of actively damped response of structures with respect to errors in the structural modeling is studied. Two ways of representing errors are considered. The first approach assumes errors in the form of spatial variations (or imperfections) in the assumed mass and stiffness properties of the structures. The second approach assumes errors due to such factors as unknown joint stiffnesses, discretization errors, and nonlinearities. These errors are represented here as discrepancies between experimental and analytical mode shapes and frequencies. The actively damped system considered here is a direct-rate feedback regulator based on a number of colocated velocity sensors and force actuators. The response of the controlled structure is characterized by the eigenvalues of the closed-loop system. The effects of the modeling errors are thus presented as the sensitivity of the eigenvalues of the closed-loop system. Results are presented for two examples: (1) a three-span simply supported beam controlled by three sensors and actuators, and (2) a laboratory structure consisting of a cruciform beam supported by cables.

Haftka, Raphael T.↗

On the monochromatic brightness variations of comets

The nonsteady development of the multispecies atmosphere of a long-period comet approaching the sun is modeled, and the associated monochromatic brightness variations are discussed in terms of the varying abundances of the atmospheric species. A simple homogeneous clathrate structure, 85% of which is H2O, provides the representation for the cometary nucleus. Monochromatic brightness variations of the emission bands due to OH, CN and C2 are found to follow the law proposed by Levin (1943) for small heliocentric distance where a quasi-steady approximation is valid. However, the large time-scales of residence for emitting species leads to a breakdown of Levin's law at large heliocentric distances.

Mendis, D. A.↗

An Inversion Analysis of Recent Variability in Natural CO2 Fluxes Using GOSAT and In Situ Observations

About one-half of the global CO2 emissions from fossil fuel combustion and deforestation accumulates in the atmosphere, where it contributes to global warming. The rest is taken up by vegetation and the ocean. The precise contribution of the two sinks, and their location and year-to-year variability are, however, not well understood. We use two different approaches, batch Bayesian synthesis inversion and variational data assimilation, to deduce the global spatiotemporal distributions of CO2 fluxes during 2009-2010. One of our objectives is to assess different sources of uncertainties in inferred fluxes, including uncertainties in prior flux estimates and observations, and differences in inversion techniques. For prior constraints, we utilize fluxes and uncertainties from the CASA-GFED model of the terrestrial biosphere and biomass burning driven by satellite observations and interannually varying meteorology. We also use measurement-based ocean flux estimates and two sets of fixed fossil CO2 emissions. Here, our inversions incorporate column CO2 measurements from the GOSAT satellite (ACOS retrieval, filtered and bias-corrected) and in situ observations (individual flask and afternoon-average continuous observations) to estimate fluxes in 108 regions over 8-day intervals for the batch inversion and at 3 x 3.75 weekly for the variational system. Relationships between fluxes and atmospheric concentrations are derived consistently for the two inversion systems using the PCTM atmospheric transport model driven by meteorology from the MERRA reanalysis. We compare the posterior fluxes and uncertainties derived using different data sets and the two inversion approaches, and evaluate the posterior atmospheric concentrations against independent data including aircraft measurements. The optimized fluxes generally resemble those from other studies. For example, the results indicate that the terrestrial biosphere is a net CO2 sink, and a GOSAT-only inversion suggests a shift in the global sink from the tropics south to the north relative to the prior and to an in-situ-only inversion. We also find a smaller terrestrial sink in higher-latitude northern regions in boreal summer of 2010 relative to 2009.

Carbon dioxide↗

Partitioned exponential methods for coupled multiphysics systems

Multiphysics problems involving two or more coupled physical phenomena are ubiquitous in science and engineering. This work develops a new partitioned exponential approach for the time integration of multiphysics problems. After a possible semi-discretization in space, the class of problems under consideration is modeled by a system of ordinary differential equations where the right-hand side is a summation of two component functions, each corresponding to a given set of physical processes. The partitioned-exponential methods proposed herein evolve each component of the system via an exponential integrator, and information between partitions is exchanged via coupling terms. Here, the traditional approach to constructing exponential methods, based on the variation-of-constants formula, is not directly applicable to partitioned systems. Rather, our approach to developing new partitioned-exponential families is based on a general-structure additive formulation of the schemes. Two method formulations are considered, one based on a linear-nonlinear splitting of the right hand component functions, and another based on approximate Jacobians. The paper develops classical (non-stiff) order conditions theory for partitioned exponential schemes based on particular families of T-trees and B-series theory. Several practical methods of third order are constructed that extend the Rosenbrock-type and EPIRK families of exponential integrators. Several implementation optimizations specific to the application of these methods to reaction-diffusion systems are also discussed. Numerical experiments reveal that the new partitioned-exponential methods can perform better than traditional unpartitioned exponential methods on some problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Hybrid Gradient Method to Designing Bayesian Experiments for Implicit Models

Bayesian experimental design (BED) aims at designing an experiment to maximize the information gathering from the collected data. The optimal design is usually achieved by maximizing the mutual information (MI) between the data and the model parameters. When the analytical expression of the MI is unavailable, e.g.,having implicit models with intractable data distributions, a neural network-based lower bound of the MI was recently proposed and a gradient ascent method was used to maximize the lower bound [1]. However, the approach in [1] requires a pathwise sampling path to compute the gradient of the MI lower bound with respect to the design variables, and such a pathwise sampling path is usually inaccessible for implicit models. In this work, we propose a hybrid gradient approach that leverages recent advances in variational MI estimator and evolution strategies (ES)combined with black-box stochastic gradient ascent (SGA) to maximize the MI lower bound. This allows the design process to be achieved through a unified scalable procedure for implicit models without sampling path gradients. Several experiments demonstrate that our approach significantly improves the scalability of BED for implicit models in high-dimensional design space.

Zhang, Jiaxin↗

A circuit-generated quantum subspace algorithm for the variational quantum eigensolver

Recent research has shown that wavefunction evolution in real and imaginary time can generate quantum subspaces with significant utility for obtaining accurate ground state energies. Inspired by these methods, we propose combining quantum subspace techniques with the variational quantum eigensolver (VQE). In our approach, the parameterized quantum circuit is divided into a series of smaller subcircuits. The sequential application of these subcircuits to an initial state generates a set of wavefunctions that we use as a quantum subspace to obtain high-accuracy groundstate energies. We call this technique the circuit subspace variational quantum eigensolver (CSVQE) algorithm. By benchmarking CSVQE on a range of quantum chemistry problems, we show that it can achieve significant error reduction in the best case compared to conventional VQE, particularly for poorly optimized circuits, greatly improving convergence rates. Furthermore, we demonstrate that when applied to circuits trapped at local minima, CSVQE can produce energies close to the global minimum of the energy landscape, making it a potentially powerful tool for diagnosing local minima.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Nonlocal Feature-Driven Exemplar-Based Approach for Image Inpainting

Here, we present a nonlocal variational image completion technique which admits simultaneous inpainting of multiple structures and textures in a unified framework. The recovery of geometric structures is achieved by using general convolution operators as a measure of behavior within an image. These are combined with a nonlocal exemplar-based approach to exploit the self-similarity of an image in the selected feature domains and to ensure the inpainting of textures. We also introduce an anisotropic patch distance metric to allow for better control of the feature selection within an image and present a nonlocal energy functional based on this metric. Finally, we derive an optimization algorithm for the proposed variational model and examine its validity experimentally with various test images.

97 MATHEMATICS AND COMPUTING↗

Benchmarking Quantum Chemistry Computations with Variational, Imaginary Time Evolution, and Krylov Space Solver Algorithms

Quantum chemistry is a key application area for noisy-intermediate scale quantum (NISQ) devices, and therefore serves as an important benchmark for current and future quantum computer performance. Previous benchmarks in this field have focused on variational methods for computing ground and excited states of various molecules, including a benchmarking suite focused on the performance of computing ground states for alkali-hydrides under an array of error mitigation methods. State-of-the-art methods to reach chemical accuracy in hybrid quantum-classical electronic structure calculations of alkali hydride molecules on NISQ devices from IBM are outlined here. Here it is demonstrated how to extend the reach of variational eigensolvers with symmetry preserving Ansätze. Next, it is outlined how to use quantum imaginary time evolution and Lanczos as a complementary method to variational techniques, highlighting the advantages of each approach. Finally, a new error mitigation method is demonstrated which uses systematic error cancellation via hidden inverse gate constructions, improving the performance of typical variational algorithms. These results show that electronic structure calculations have advanced rapidly, to routine chemical accuracy for simple molecules, from their inception on quantum computers a few short years ago, and they point to further rapid progress to larger molecules as the power of NISQ devices grows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Scattering of Low Energy Electrons and Positrons from Hydrogenic Systems and Applications

While the electron scatters from the target, the target no longer stays in its original form. One of the first methods to take into account the distortion of the target at low incident energies is the method of polarized orbitals. In this method the wave function for the process is written using the first-order perturbation theory and the equation for the scattering function is derived from the Schradinger equation. This method has been very successful in calculating the phase shifts and therefore the cross sections at various energies. The total wave function can be used to calculate photoionization cross sections. The disadvantage of this approach is that the method is not variational and therefore does not provide bounds on the phase shifts. These difficulties can be overcome by using the Feshbach projection operator formalism. This approach has been employed for the scattering of electrons and positrons from targets. Results of various calculations will be discussed.

Bhatia, Anand K.↗

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Using Particle Image Velocimetry to Determine Turbulence Model Parameters

The primary parameter of a standard k-ϵ model, C μ , was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine C μ . Considerable spatial variation was observed in C μ , with the highest values found in regions where turbulent kinetic energy was relatively low but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant C μ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of C μ based on known state variables may lead to improved simulation results without the need for calibration.

42 ENGINEERING↗

High Reynolds Number Configuration Development of a High-Lift Airfoil

An experimental program has been conducted to assess performance of a transport multielement airfoil at flight Reynolds numbers. The studies were performed at chord Reynolds numbers as high as 16 million in the NASA Langley Low Turbulence Pressure Tunnel. Sidewall boundary-layer control to enforce flow two dimensionality was provided via an endplate suction system. The basic airfoil was an 11.55 percent thick supercritical airfoil representative of the stall critical station of a new-generation transport aircraft wing. The multielement airfoil was configured as a three-element airfoil with slat and flap chord ratios of 14.48 percent and 30 percent respectively. Testing focused on the development of landing configurations with high maximum lift capability and the assessment of Reynolds and Mach number effects. Also assessed were high-lift performance effects due to devices such as drooped spoilers and trailing-edge wedges. The present experimental studies revealed significant effects on high-lift airfoil performance due to Reynolds and Mach number variations and favorable lift increments at approach angles of attack due to the use of drooped spoilers or trailing-edge wedges. However, no substantial improvements in maximum lift capability were identified. A recently developed high performance single-segment flap was also tested and results indicated considerable improvements in lift and drag performance over existing airfoils. Additionally, it was found that this new flap shape at its optimum rigging was less sensitive to Reynolds number variations than previous designs.

Walter O Valarezo↗

Modulation of cosmic-ray electrons.

An analytical approach is presented that provides qualitative insight into both observations and numerical calculations of cosmic-ray electron modulation and embodies, as a natural feature of a specific model, two cosmic-ray components whose temporal variations are uncorrelated. In the light of this approach, it is shown that the origin and time variations of the steep spectrum of electrons observed below 20 MeV may be explained by a simple method in which the spectrum of interplanetary cosmic rays is decomposed at low energies into two independently varying components.

Earl, J. A.↗

Satellite Altimetry for a Global Ocean Observing System

Space-age technologies have made satellite remote sensing a powerful new tool to study the Earth on a global scale. However, the opacity of the ocean to electromagnetic sensing has limited spaceborne measurements to the properties of the surface layer of the ocean (such as sea surface temperature and color). The radar altimetric measurement of the height of the sea surface relative to the geoid, the dynamic topography of the ocean, is a very useful quantity for studying the circulation of the ocean. The ability of measuring dynamic topography from space makes satellite altimetry a uniquely useful remote sensing technique because dynamic topography reflects oceanic processes not only at the surface but at depths as well. A simple analysis shows that a one centimeter tilt in the dynamic topography is associated with a mass transport of 1-7 Sv (1Sv= 1 million tons per second) in the open ocean depending on the vertical distribution of current velocity. Such a magnitude is an appreciable fraction of the transport of the Florida Current (circa 30 Sv), for instance. TOPEX/POSEIDON has demonstrated the capability of measuring the time variation of sea level with accuracy approaching to 2 cm when the data are averaged over boxes with several hundred kilometers on each side. The data set has been used for studying ocean circulation phenomena with a wide range of scales, ranging from fast-changing barotropic variability to seasonal and interannual variability such as El Nino and La Nina. The long record of precise measurement of global sea level has also showed great promise for monitoring the variation of mean sea level, an effective indicator of global climate change. Continuation of satellite altimetry missions with capability matching or better than that of TOPEX/POSEIDON should be included as a key component of a Global Ocean Observing System. NASA and CNES have committed to continuing the measurement of TOPEX/POSEIDON with a series of follow-on missions called Jason. The first of the series, Jason-1, is scheduled for launch in May, 2000. Such a series of missions will provide a key data stream for both research and practical applications and benefit the objectives of global programs such as CLIVAR and GODAE.

Fu, Lee-Lueng↗

Modeling the Stress Strain Behavior of Woven Ceramic Matrix Composites

Woven SiC fiber reinforced SiC matrix composites represent one of the most mature composite systems to date. Future components fabricated out of these woven ceramic matrix composites are expected to vary in shape, curvature, architecture, and thickness. The design of future components using woven ceramic matrix composites necessitates a modeling approach that can account for these variations which are physically controlled by local constituent contents and architecture. Research over the years supported primarily by NASA Glenn Research Center has led to the development of simple mechanistic-based models that can describe the entire stress-strain curve for composite systems fabricated with chemical vapor infiltrated matrices and melt-infiltrated matrices for a wide range of constituent content and architecture. Several examples will be presented that demonstrate the approach to modeling which incorporates a thorough understanding of the stress-dependent matrix cracking properties of the composite system.

Morscher, Gregory N.↗

The radiobiological implications of statistical variations in energy deposition by ionizing radiations

Traditional approaches to microdosimetry, the fundamental physics of energy deposition, the importance of statistical processes, an illustration of possible radiobiological interpretation, and modeling based on microdosimetric concepts are discussed. Emphasis is on the inadequacies in linear energy transfer (LET) theory. For many reasons, concepts based on averaging may not be applicable to ionizing radiation absorption by and damage to small biological targets.

Baily, N. A.↗

Investigation of Improved Wall Heat Flux Behavior Observed When Computing Hypersonic Flow Using a Face-averaged Nodal-gradient Approach On Tetrahedral Grids

We investigate the origin of a significant reduction in computed numerical “noise” in the spatial variation of computed wall heat flux observed when using a node-centered weighted least squares(WLSQ) gradient approach when solving the Reynolds-averaged Navier-Stokes (RANS) equations for hypersonic boundary layer flows on highly stretched tetrahedral grids. This investigation is conducted using a 2nd-order, cell-centered, finite-volume, discretization which we employ to contrast results obtained using two WLSQ gradient approaches: a conventional 3-D, cell-centered, node-neighbor, WLSQ gradient (NN-CCG) approach and a 3-D, face-averaged, node-centered gradient (F-ANG)approach. The F-ANG approach, is shown to produce significantly less “noise” in the spatial variation of wall heat flux, than the NN-CCG approach and we demonstrate the cause of this “noise” via a series of numerical experiments where we compute hypersonic laminar and turbulent boundary layer flows.We further investigate the level of grid resolution required for a tetrahedral grid to approach the fidelity of a grid resolved hexahedral grid solution obtained using an equivalent structured grid method.

Jeffery A. White↗

Engineering quantitative stomatal trait variation and local adaptation potential by cis‐regulatory editing

Summary Cis‐regulatory element editing can generate quantitative trait variation that mitigates extreme phenotypes and harmful pleiotropy associated with coding sequence mutations. Here, we applied a multiplexed CRISPR/Cas9 approach, informed by bioinformatic datasets, to generate genotypic variation in the promoter ofOsSTOMAGEN, a positive regulator of rice stomatal density. Engineered genotypic variation corresponded to broad and continuous variation in stomatal density, ranging from 70% to 120% of wild‐type stomatal density. This panel of stomatal variants was leveraged in physiological assays to establish discrete relationships between stomatal morphological variation and stomatal conductance, carbon assimilation and intrinsic water use efficiency in steady‐state and fluctuating light conditions. Additionally, promoter alleles were subjected to vegetative drought regimes to assay the effects of the edited alleles on developmental response to drought. Notably, the capacity for drought‐responsive stomatal density reprogramming instomagenand two cis‐regulatory edited alleles was reduced. Collectively our data demonstrate that cis‐regulatory element editing can generate near‐isogenic trait variation that can be leveraged for establishing relationships between anatomy and physiology, providing a basis for optimizing traits across diverse environments.

Biotechnology & Applied Microbiology↗