Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “VARIATIONAL METHODS”

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,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

System and method for focal-plane angular-spatial illuminator/detector (fasid) design for improved graded index lenses

The present disclosure relates to a method for imaging an optical signal received by a graded index (GRIN) optical element to account for known variations in a graded index distribution of the GRIN optical element. The method may involve using a plurality of optical detector elements to receive optical rays received by the GRIN optical element at a plane, where the plane forms a part of the GRIN optical element or is downstream of the GRIN optical element relative to a direction of propagation of the optical rays. The optical rays are then traced to a plurality of additional specific locations on the plane based on the known variations in the graded index distribution of the GRIN optical element. A processor may be used to determine information on both an intensity and an angle of the received optical rays at each one of the plurality of specific locations on the plane of the GRIN optical element.

Feigenbaum, Eyal↗
NASA NTRS1986

On the flight dynamics of aeroelastic vehicles

The nonlinear equations of motion for an elasic airplane are developed from first principles. Lagrange's Equation and the Principle of Virtual Work are used to generate the equations of motion and strip theory is then employed to obtain closed form integral expressions for the generalized forces. The inertial coupling is minimized by appropriate choice of the body reference axes and by making use of free vibration modes of the body. In addition, particular attention is paid to the simplifying assumptions used during the development of the equations of motion. A unique aspect of this modeling process is that since the generalized aerodynamic forces are determined from closed form, analytic expressions, this method can be used to gain insight into the effects of parameter variations not easily obtained from numerical models. A numerical example is also presented in which the modeling method is applied to a representative elastic aircraft. The model is used to address the effects of aerodynamic coupling which occurs between the rigid body degrees of freedom and the elastic degrees of freedom. Finally, model simplification is addressed and two methods are evaluated. The resulting frequency responses are compared.

Schmidt, D. K.↗
NASA NTRS1985

Variation in efficiency of parallel algorithms

The present study has the objective to investigate some iterative parallel-processor linear equation solving algorithms with respect to efficiency for analyses of typical linear engineering systems. Attention is given to a set of n linear equations, Ku = p, where K = an n x n positive definite, sparsely populated, symmetric matrix, u = an n x 1 vector of unknown responses, and p = an n x 1 vector of prescribed constants. This study is concerned with a hybrid method in which iteration is used to solve the problem, while a direct method is used on the local processor level. Variations in the efficiency of parallel algorithms are explored. Measures of the efficiency are based on computer experiments regarding the algorithms. For all the algorithms, the wall clock time is found to decrease as the number of processors increases.

Hayashi, A.↗
NASA NTRS1994

A Fully Implicit Time Accurate Method for Hypersonic Combustion: Application to Shock-induced Combustion Instability

A new fully implicit, time accurate algorithm suitable for chemically reacting, viscous flows in the transonic-to-hypersonic regime is described. The method is based on a class of Total Variation Diminishing (TVD) schemes and uses successive Gauss-Siedel relaxation sweeps. The inversion of large matrices is avoided by partitioning the system into reacting and nonreacting parts, but still maintaining a fully coupled interaction. As a result, the matrices that have to be inverted are of the same size as those obtained with the commonly used point implicit methods. In this paper we illustrate the applicability of the new algorithm to hypervelocity unsteady combustion applications. We present a series of numerical simulations of the periodic combustion instabilities observed in ballistic-range experiments of blunt projectiles flying at subdetonative speeds through hydrogen-air mixtures. The computed frequencies of oscillation are in excellent agreement with experimental data.

Yungster, Shaye↗
NASA NTRS1998

Microstructural Development of Directionally Solidified Hg(1-x)Zn(x)Se Alloys

Hg(1-x)Zn(x)Se alloys have been studied as an alternative to HgCdTe for the detection of electromagnetic radiation, because the shorter ZnSe and HgSe bonds were predicted to improve lattice stability. Previous studies showed that optical properties were stable with time and that the microhardness was higher than for HgTe based alloys. However, for this material to be commercially viable, compositional variations must be reduced, and the distribution of defects induced by the growth process must be understood. The basic structural properties of bulk materials were characterized, and dislocation etch pit densities were measured to determine whether the addition of Zn to the HgSe lattice reduced the potential for dislocation formation. Hg(0.9)Zn(0.1)Se alloys were directionally solidified using a modified Bridgman-Stockbarger method and in an applied magnetic field. Radial compositional variations were greatly reduced when solidification occurred in an applied magnetic field. The combination of a convex liquid-solid interface shape and a method for reduced wetting produced boules that were a single crystal after reaching full diameter. Observed surface features indicated ampoule wetting could be eliminated using a graphite getter. Microstructural characteristics were greatly improved over HgCdTe alloys, In six boules, a total of only one twin was observed, and no inclusions were found. A method for polishing and producing dislocation etch pits was developed for these alloys, revealing dislocation etch pit densities one to two orders of magnitude less than HgTe based alloys. A kink in the thermal profile during processing of one boule generated more dislocations than did lattice mismatch due to compositional variations. The results of this investigation indicate this alloy has improved microstructural properties and resistance to dislocation formation compared with similar II-VI alloys, and should be further investigated for the detection of electromagnetic radiation.

Cobb, S. D.↗
DOE OSTI2025

Multimodal Defect Imaging of Pure Tungsten Components Fabricated via Electron Beam Powder Bed Fusion

The utilization of additive manufacturing (AM) techniques for refractory materials in high-temperature environments has significantly expanded because of the ability to fabricate geometrically complex components. Electron beam powder bed fusion (EB-PBF), which provides lower residual stress, a cleaner vacuum environment, and better efficiency for high melting point, is one of the best-suited AM methods to produce advanced refractory components. However, the property variation attributed to the heterogeneous microstructure and process-induced defects has hindered the widespread adoption of EB-PBF-produced material like tungsten. While numerous in-situ monitoring and defect detection methods have been demonstrated for EB-PBF, a workflow that compares and evaluates process-induced abnormalities from different imaging perspectives is still limited. This study examines a feature-embedded tungsten component manufactured via the EB-PBF process to demonstrate the defect detection capabilities of a multimodal defect imaging workflow. The predefined and process-induced defects are evaluated by harnessing various imaging techniques, including in-situ electron imaging, layerwise near-infrared (NIR) imaging, post-build high-energy x-ray computed tomography (CT), and conventional destructive metallography. The results highlight the strengths and limitations of distinctive defect imaging techniques concerning specific defect types, sizes, and conditions. It was found that electron imaging can provide more abnormal detection capabilities while maintaining a higher measuring accuracy, against the conventional metallography in this case study, compared with NIR and CT imaging techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2025

3D strain field reconstruction by inversion of dynamical scattering

Strain governs not only the mechanical response of materials but also their electronic, optical, and catalytic properties. For this reason, the measurement of the 3D strain field is crucial for a detailed understanding and for further development of material properties through strain engineering. However, measuring strain variations along the electron beam direction has remained a major challenge for (scanning-) transmission electron microscopy (S/TEM). In this article, we present a method for 3D strain field determination using 4D-STEM. The method is based on the inversion of dynamical diffraction effects, which occur at strain field variations along the beam direction. We test the method against simulated data with a known ground truth and demonstrate its application to an experimental 4D-STEM dataset from an inclined pseudomorphically grown Al0.47Ga0.53N layer.

Niermann, Laura↗
NASA NTRS1987

A probabilistic Hu-Washizu variational principle

A Probabilistic Hu-Washizu Variational Principle (PHWVP) for the Probabilistic Finite Element Method (PFEM) is presented. This formulation is developed for both linear and nonlinear elasticity. The PHWVP allows incorporation of the probabilistic distributions for the constitutive law, compatibility condition, equilibrium, domain and boundary conditions into the PFEM. Thus, a complete probabilistic analysis can be performed where all aspects of the problem are treated as random variables and/or fields. The Hu-Washizu variational formulation is available in many conventional finite element codes thereby enabling the straightforward inclusion of the probabilistic features into present codes.

Liu, W. K.↗
DOE OSTI2021

Fully variational incremental CASSCF

We report the complete-active-space self-consistent field (CASSCF) method is a canonical electronic structure theory that holds a central place in conceptualizing and practicing first principles simulations. For application to realistic molecules, however, the CASSCF must be approximated to circumvent its exponentially scaling computational costs. Applying the many-body expansion - also known as the method of increments - to CASSCF (iCASSCF) has been shown to produce a polynomially scaling method that retains much of the accuracy of the parent theory and is capable of treating full valence active spaces. Due to an approximation made in the orbital gradient, the orbital parameters of the original iCASSCF formulation could not be variationally optimized, which limited the accuracy of its nuclear gradient. Herein, a variational iCASSCF is introduced and implemented, where all parameters are fully optimized during energy minimization. This method is able to recover electronic correlations from the full valence space in large systems, produce accurate gradients, and optimize stable geometries as well as transition states. Demonstrations on challenging test cases, such as the oxoMn(salen)Cl complex with 84 electrons in 84 orbitals and the automerization of cyclobutadiene, show that the fully variational iCASSCF is a powerful tool for describing challenging molecular chemistries.

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

Apparent sizes and spectral line profiles for spherical ('three-dimensional') astrophysical masers

Calculations are performed for the frequency-dependent transport of radiation in a uniform spherical maser. Coupling between the intensities at different frequencies is treated. Spherical masers are commonly utilized as idealized geometries in which to assess possible deviations (sometimes, 'three-dimensional effects') from the approximation of a linear geometry. We find that the spectral line profile rebroadens in approximately the same manner as that for a linear maser so that the relationship between the line breadth and the luminosity, beaming angle, and degree of saturation is essentially unchanged. The variation of the apparent size of a spherical maser with frequency is not found to be significant at frequencies within the line profile at which the intensity is appreciable. The 'standard approximation' is adequate for obtaining most of the basic properties of spherical masers that have been examined including the relationship between the pumping (or molecular populations) and the intensity of maser radiation. The improved methods employed here are, however, necessary to obtain variation of the apparent size with frequency in the saturated regime. These conclusions differ from those of previous investigators. The calculations here are more reliable than those of previous investigators because these are self-consistent and avoid certain inaccurate simplifications.

Emmering, Robert T.↗
NASA NTRS2001

A High-Order Finite Spectral Volume Method for Conservation Laws on Unstructured Grids

A time accurate, high-order, conservative, yet efficient method named Finite Spectral Volume (FSV) is developed for conservation laws on unstructured grids. The concept of a 'spectral volume' is introduced to achieve high-order accuracy in an efficient manner similar to spectral element and multi-domain spectral methods. In addition, each spectral volume is further sub-divided into control volumes (CVs), and cell-averaged data from these control volumes is used to reconstruct a high-order approximation in the spectral volume. Riemann solvers are used to compute the fluxes at spectral volume boundaries. Then cell-averaged state variables in the control volumes are updated independently. Furthermore, TVD (Total Variation Diminishing) and TVB (Total Variation Bounded) limiters are introduced in the FSV method to remove/reduce spurious oscillations near discontinuities. A very desirable feature of the FSV method is that the reconstruction is carried out only once, and analytically, and is the same for all cells of the same type, and that the reconstruction stencil is always non-singular, in contrast to the memory and CPU-intensive reconstruction in a high-order finite volume (FV) method. Discussions are made concerning why the FSV method is significantly more efficient than high-order finite volume and the Discontinuous Galerkin (DG) methods. Fundamental properties of the FSV method are studied and high-order accuracy is demonstrated for several model problems with and without discontinuities.

Wang, Z. J.↗
DOE OSTI2024

Error mitigation, optimization, and extrapolation on a trapped-ion testbed

Current noisy intermediate-scale quantum (NISQ) trapped-ion devices are subject to errors which can significantly impact the accuracy of calculations if left unchecked. A form of error mitigation called zero noise extrapolation (ZNE) can decrease an algorithm’s sensitivity to these errors without increasing the number of required qubits. Here we explore different methods for integrating this error mitigation technique into the Variational Quantum Eigensolver (VQE) algorithm for calculating the ground state of the HeH + molecule at 0.8 Å in the presence of experimental noise. Using the Quantum Scientific Computing Open User Testbed (QSCOUT) trapped-ion device, we test three methods of scaling noise for extrapolation: time stretching the two-qubit gates, scaling the sideband detuning parameter, and inserting two-qubit gate identity operations into the ansatz circuit. We find that time stretching and sideband detuning scaling fail to scale the noise on our particular hardware in a way that can be extrapolated to zero noise. Scaling our noise with global gate identity insertions and extrapolating after variational optimization, we achieve error suppression of 96.8%, resulting in an energy estimate within –0.004 ± 0.04 hartree of the ground state energy. This is an improvement, but still outside the chemical accuracy threshold of 0.0016 hartree. Furthermore, our results show that the efficacy of this error mitigation technique depends on choosing the correct implementation for a given device architecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2019

Evaluation of an Alternate Method for Determining Yield Strength Offset Values for Selective Laser Sintered Polymeric Materials

Due to the unique characteristics of Additively Manufactured (AM) polymeric materials, typical mechanical strength characterization methods such as those commonly used for traditionally-processed polymers or composite materials can produce results that do not accurately represent material capabilities. In order to characterize mechanical properties of these materials, new test and analysis methods are required. As part of the National Aeronautics and Space Administration (NASA) Advanced Composites Project (ACP), Boeing has evaluated true yield testing as an alternative or complimentary test to 0.2% offset yield testing for determining appropriate yield strength values of polymer materials. Previous testing has shown high strain, low modulus polymer materials such as selective laser sintered (SLS) Nylon 11 at elevated temperatures produce large variations in yield strength. The true yield test method was successful in finding the applied strain level when yield commences and appears to offer an increase in data robustness.

Henry, Chelsey↗
DOE OSTI2024

Path-BigBird: An AI-Driven Transformer Approach to Classification of Cancer Pathology Reports

PURPOSE Surgical pathology reports are critical for cancer diagnosis and management. To accurately extract information about tumor characteristics from pathology reports in near real time, we explore the impact of using domain-specific transformer models that understand cancer pathology reports. METHODS We built a pathology transformer model, Path-BigBird, by using 2.7 million pathology reports from six SEER cancer registries. We then compare different variations of Path-BigBird with two less computationally intensive methods: Hierarchical Self-Attention Network (HiSAN) classification model and an offthe-shelf clinical transformer model (Clinical BigBird). We use five pathology information extraction tasks for evaluation: site, subsite, laterality, histology, and behavior. Model performance is evaluated by using macro and micro F 1 scores. RESULTS We found that Path-BigBird and Clinical BigBird outperformed the HiSAN in all tasks. Clinical BigBird performed better on the site and laterality tasks. Versions of the Path-BigBird model performed best on the two most difficult tasks: subsite (micro F 1 score of 72.53, macro F 1 score of 35.76) and histology (micro F 1 score of 80.96, macro F 1 score of 37.94). The largest performance gains over the HiSAN model were for histology, for which a Path-BigBird model increased the micro F 1 score by 1.44 points and the macro F 1 score by 3.55 points. Overall, the results suggest that a Path-BigBird model with a vocabulary derived from wellcurated and deidentified data is the best-performing model. CONCLUSION The Path-BigBird pathology transformer model improves automated information extraction from pathology reports. Although Path-BigBird outperforms Clinical BigBird and HiSAN, these less computationally expensive models still have utility when resources are constrained.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Seismic Monitoring at the Farnsworth CO2-EOR Field Using Time-Lapse Elastic-Waveform Inversion of 3D-3C VSP Data

During the Development Phase of the U.S. Southwest Regional Partnership on Carbon Sequestration, supercritical CO2 was continuously injected into the deep oil-bearing Morrow B formation of the Farnsworth Unit in Texas for Enhanced Oil Recovery (EOR). The project injected approximately 94 kilotons of CO2 to study geologic carbon storage during CO2-EOR. A three-dimensional (3D) surface seismic dataset was acquired in 2013 to characterize the subsurface structures of the Farnsworth site. Following this data acquisition, the baseline and three time-lapse three-dimensional three-component (3D-3C) vertical seismic profiling (VSP) data were acquired at a narrower surface area surrounding the CO2 injection and oil/gas production wells between 2014 and 2017 for monitoring CO2 injection and migration. With these VSP datasets, we inverted for subsurface velocity models to quantitatively monitor the CO2 plume within the Morrow B formation. We first built 1D initial P-wave (Vp) and S-wave (Vs) velocity models by upscaling the sonic logs. We improved the deep region of the Vp and Vs models by incorporating the deep part of a migration velocity model derived from the 3D surface seismic data. We improved the shallow region of 3D Vp and Vs models using 3D traveltime tomography of first arrivals of VSP downgoing waves. We further improved the 3D baseline velocity models using elastic-waveform inversion (EWI) of the 3D baseline VSP upgoing data. Our advanced EWI method employs alternative tomographic and conventional gradients and total-variation-based regularization to ensure the high-fidelity updates of the 3D baseline Vp and Vs models. We then sequentially applied our 3D EWI method to the three time-lapse datasets to invert for spatiotemporal changes of Vp and Vs in the reservoir. Our inversion results reveal the volumetric changes of the time-lapse Vp and Vs models and show the evolution of the CO2 plume from the CO2 injection well to the oil/gas production wells.

42 ENGINEERING↗
NASA NTRSDate not supplied

Nonlinear Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites Using Mechanics of Structure Genome

The objective of this paper is to develop a micromechanics approach to homogenizing finitely deformed viscoelastic-viscoplastic composites using the mechanics of structure genome. The incremental constitutive relation for glassy polymers, formulated in the spatial configuration, is implemented in the present approach.This involves (1) pulling-back the constitutive model to the material configuration and (2)choosing the deformation gradient tensor and the first Piola–Kirchhoff stress tensor as the strain and the stress measures during homogenization, respectively. An Euler–Newton predictor–corrector method is developed for homogenization. Each step involves formulating a variational statement using the mechanics of structure genome, discretizing the statement in a finite-dimensional space, and solving the problem using an Euler/multilevel Newton method. The present approach is demonstrated by homogenizing fiber- and particle-reinforced composites undergoing uniaxial, biaxial, or shear deformation, at different stain rates.

Multi-scale modeling, High Strain Composites, Visc↗
NASA NTRS1972

Application of a one-strip integral method to the unsteady supersonic aerodynamics of an inclined flat surface

The method of integral relations is applied in a one-strip approximation to the perturbation equations governing small motions of an inclined, sharp-edged, flat surface about the mean supersonic steady flow. Algebraic expressions for low reduced-frequency aerodynamics are obtained and a set of ordinary differential equations are obtained for general oscillatory motion. Results are presented for low reduced-frequency aerodynamics and for the variation of the unsteady forces with frequency. The method gives accurate results for the aerodynamic forces at low reduced frequency which are in good agreement with available experimental data. However, for cases in which the aerodynamic forces vary rapidly with frequency, the results are qualitatively correct, but of limited accuracy. Calculations indicate that for a range of inclination angles near shock detachment such that the flow in the shock layer is low supersonic, the aerodynamic forces vary rapidly both with inclination angle and with reduced frequency.

Bennett, R. M.↗
DOE OSTI2022

Dark Energy Survey Year 3 results: marginalization over redshift distribution uncertainties using ranking of discrete realizations

ABSTRACT Cosmological information from weak lensing surveys is maximized by sorting source galaxies into tomographic redshift subsamples. Any uncertainties on these redshift distributions must be correctly propagated into the cosmological results. We present hyperrank, a new method for marginalizing over redshift distribution uncertainties, using discrete samples from the space of all possible redshift distributions, improving over simple parametrized models. In hyperrank, the set of proposed redshift distributions is ranked according to a small (between one and four) number of summary values, which are then sampled, along with other nuisance parameters and cosmological parameters in the Monte Carlo chain used for inference. This approach can be regarded as a general method for marginalizing over discrete realizations of data vector variation with nuisance parameters, which can consequently be sampled separately from the main parameters of interest, allowing for increased computational efficiency. We focus on the case of weak lensing cosmic shear analyses and demonstrate our method using simulations made for the Dark Energy Survey (DES). We show that the method can correctly and efficiently marginalize over a wide range of models for the redshift distribution uncertainty. Finally, we compare hyperrank to the common mean-shifting method of marginalizing over redshift uncertainty, validating that this simpler model is sufficient for use in the DES Year 3 cosmology results presented in companion papers.

79 ASTRONOMY AND ASTROPHYSICS↗