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 847 records · Page 47

Model Error Estimation for the CPTEC Eta Model

Statistical data assimilation systems require the specification of forecast and observation error statistics. Forecast error is due to model imperfections and differences between the initial condition and the actual state of the atmosphere. Practical four-dimensional variational (4D-Var) methods try to fit the forecast state to the observations and assume that the model error is negligible. Here with a number of simplifying assumption, a framework is developed for isolating the model error given the forecast error at two lead-times. Two definitions are proposed for the Talagrand ratio tau, the fraction of the forecast error due to model error rather than initial condition error. Data from the CPTEC Eta Model running operationally over South America are used to calculate forecast error statistics and lower bounds for tau.

Tippett, Michael K.↗

Complex Correlation Calculation of e(-) - H Total Cross Sections

Calculation of e(-) - H total and elastic partial wave cross sections is being carried out using the complex correlation variational T-matrix method. In this preliminary study, elastic partial wave phase shifts are calculated with the correlation functions which are confined to be real. In that case the method reduces to the conventional optical potential approach with 2 projection operators. The number of terms in the Hylleraas-type wave function for the S-1 phase shifts is 95 while for the S-3 it is 56, except for k = 0.8 where it is 84. Our results, which are rigorous lower bounds, are seen to be in general agreement with those of Schwartz, but they are of greater accuracy and outside of his error limits for k = 0.3 and 0.4 for S-1. The main aim of this approach is the application to higher energy scattering. By virtue of the complex correlation functions, the T-matrix is not unitary so that elastic and total scattering cross sections are independent of each other. Our results will be compared specifically with those of Bray and Stelbovics.

Bhatia, A. K.↗

Complex Correlation Calculation of e-H Total Cross Sections

Calculation of e - H total and elastic partial wave cross sections is being carried out using the complex correlation variational T-matrix method. In this preliminary study, elastic partial wave phase shifts are calculated with the correlation functions which are confined to be real. In that case the method reduces to the conventional optical potential approach with projection operators.

Bhatia, A. K.↗

Complex Correlation Calculation of e-H Total Cross Sections

Calculation of e-H total and elastic partial wave cross sections is being carried out using the complex correlation variational T-matrix method. In this preliminary study, elastic partial wave phase shifts are calculated with the correlation functions which are confined to be real. In that case the method reduces to the conventional optical potential approach with projection operators. The number of terms in the Hylleraas-type wave function for the S phase shifts is 95 while for the S it is 56, except for k=0.8 where it is 84. Our results, which are rigorous lower bounds, are given. They are seen to be in general agreement with those of Schwartz, but they are of 0 greater accuracy and outside of his error limits for k=0.3 and 0.4 for S. The main aim of this approach' is the application to higher energy scattering. By virtue of the complex correlation functions, the T matrix is not unitary so that elastic and total scattering cross sections are independent of each other. Our results will be compared specifically with those of Bray and Stelbovics.

Bhatia, A. K.↗

Trajectories for High Specific Impulse High Specific Power Deep Space Exploration

Flight times and deliverable masses for electric and fusion propulsion systems are difficult to approximate. Numerical integration is required for these continuous thrust systems. Many scientists are not equipped with the tools and expertise to conduct interplanetary and interstellar trajectory analysis for their concepts. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. An analytical method derived in the companion paper was also evaluated. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

Simulation of Trajectories for High Specific Impulse Deep Space Exploration

Difficulties in approximating flight times and deliverable masses for continuous thrust propulsion systems have complicated comparison and evaluation of proposed propulsion concepts. These continuous thrust propulsion systems are of interest to many groups, not the least of which are the electric propulsion and fusion communities. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. The analytical method derived in the companion paper was also used to simulate the trajectory. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

Acridine orange staining reaction as an index of physiological activity in Escherichia coli

The assumption that the acridine orange (AO) color reaction may be used as an index of physiological activity was investigated in laboratory grown Escherichia coli. Spectrofluorometric observations of purified nucleic acids, ribosomes and the microscopic color of bacteriophage-infected cells stained with AO confirmed the theory that single-stranded nucleic acids emit orange to red fluorescence while those that are double-stranded fluoresce green in vivo. Bacteria growing actively in a rich medium could be distinguished from cells in stationary phase by the AO reaction. Cells from log phase appeared red, whereas those in stationary phase were green. However, this differentiation was not seen when the bacteria were grown in a minimal medium or when a variation of the staining method was used. Also, shifting bacteria in stationary phase to starvation conditions rapidly changed their AO staining reaction. Boiling and exposure to lethal concentrations of azide and formalin resulted in stationary-phase cells that appeared red after staining but bacteria killed with chlorine remained green. These findings indicate that the AO staining reaction may be suggestive of physiological activity under defined conditions. However, variables in staining and fixation procedures as well as uncertainties associated with mixed bacterial populations in environmental samples may produce results that are not consistent with the classical interpretation of this reaction. The importance of validating the putative physiological implications of this staining reaction is stressed.

NASA Discipline Number 04-10↗

Computation of Variability in the Average Thermal and Mechanical Properties of a Melt-infiltrated Sic/sic Composite

Thermal conductivity and tensile properties (elastic modulus, proportional limit strength, inplane tensile strength, and strain to failure) of a SiC/SiC composite were experimentally determined at 816 and 1204 C. Tests were performed at room temperature on the same material to obtain interlaminar shear and tensile strength. For each thermal and mechanical property 24 tests were conducted to capture the variation. A random sampling method was used to quantify the variability exhibited by the mean values of the thermal and mechanical properties. The minimum number of tests required to characterize the mean value of each thermal or mechanical property for the composite was determined by varying the group size of the random sample.

Kalluri, Sreeramesh↗

Asymptotic Energies and QED Shifts for Rydberg States of Helium

This paper reviews progress that has been made in obtaining essentially exact solutions to the nonrelativistic three-body problem for helium by a combination of variational and asymptotic expansion methods. The calculation of relativistic and quantum electrodynamic corrections by perturbation theory is discussed, and in particular, methods for the accurate calculation of the Bethe logarithm part of the electron self energy are presented. As an example, the results are applied to the calculation of isotope shifts for the short-lived 'halo' nucleus He-6 relative to He-4 in order to determine the nuclear charge radius of He-6 from high precision spectroscopic measurements carried out at the Argonne National Laboratory. The results demonstrate that the high precision that is now available from atomic theory is creating new opportunities to create novel measurement tools, and helium, along with hydrogen, can be regarded as a fundamental atomic system whose spectrum is well understood for all practical purposes.

Drake, G.W.F.↗

EMU Suit Performance Simulation

Introduction: Designing a planetary suit is very complex and often requires difficult trade‐offs between performance, cost, mass, and system complexity. To verify that new suit designs meet requirements, full prototypes must be built and tested with human subjects. However, numerous design iterations will occur before the hardware meets those requirements. Traditional draw‐prototype‐test paradigms for research and development are prohibitively expensive with today's shrinking Government budgets. Personnel at NASA are developing modern simulation techniques that focus on a human‐centric design paradigm. These new techniques make use of virtual prototype simulations and fully adjustable physical prototypes of suit hardware. This is extremely advantageous and enables comprehensive design down‐selections to be made early in the design process. Objectives: The primary objective was to test modern simulation techniques for evaluating the human performance component of two EMU suit concepts, pivoted and planar style hard upper torso (HUT). Methods: This project simulated variations in EVA suit shoulder joint design and subject anthropometry and then measured the differences in shoulder mobility caused by the modifications. These estimations were compared to human‐in‐the‐loop test data gathered during past suited testing using four subjects (two large males, two small females). Results: Results demonstrated that EVA suit modeling and simulation are feasible design tools for evaluating and optimizing suit design based on simulated performance. The suit simulation model was found to be advantageous in its ability to visually represent complex motions and volumetric reach zones in three dimensions, giving designers a faster and deeper comprehension of suit component performance vs. human performance. Suit models were able to discern differing movement capabilities between EMU HUT configurations, generic suit fit concerns, and specific suit fit concerns for crewmembers based on individual anthropometry

Cowley, Matthew S.↗

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

Crossing The Gap Using Variational Quantum Eigensolver: A Comparative Study

Within the evolving domain of quantum computational chemistry, the Variational Quantum Eigensolver (VQE) has been developed to explore not only the ground state but also the excited states of molecules. In this study, we compare the performance of Variational Quantum Deflation (VQD) and Subspace-Search Variational Quantum Eigensolver (SSVQE) methods in determining the low-lying excited states of $LiH$. Our investigation reveals that while VQD exhibits a slight advantage in accuracy, SSVQE stands out for its efficiency, allowing the determination of all low-lying excited states through a single parameter optimization procedure. We further evaluate the effectiveness of optimizers, including Gradient Descent (GD), Quantum Natural Gradient (QNG), and Adam optimizer, in obtaining $LiH$'s first excited state, with the Adam optimizer demonstrating superior efficiency in requiring the fewest iterations. Moreover, we propose a novel approach combining Folded Spectrum VQE (FS-VQE) with either VQD or SSVQE, enabling the exploration of highly excited states. We test the new approaches for finding all three $H_4$'s excited states. Folded Spectrum SSVQE (FS-SSVQE) can find all three highly excited states near $-1.0$ Ha with only one optimizing procedure, but the procedure converges slowly. In contrast, although Folded spectrum VQD (FS-VQD) gets highly excited states with individual optimizing procedures, the optimizing procedure converges faster.

Chen, I-Chi↗

The cosmic ray differential diurnal variation dependences on the zenith angle and the geomagnetic disturbance

Simultaneous and continuous muon measurements in two opposite azimuthal directions under equal zenith angles demonstrated the importance of this method for cosmic ray diurnal variation investigations. Lately these measurements were extended by means of improved telescopes. The obtained cosmic ray diurnal variations were presented as intensity differential curves. Theoretical investigations connected the properties of these curves with some interplanetary spece parameters. The harmonics of these curves were interpreted physically. Some order difference curves were introduced. In earlier works some dependences between the parameters characterizing the first and the second harmonics of the differential intensity curves and the geomagnetic activity were found. Then all measurements were carried out under only one zenith angle. The results of investigations of similar dependences using data of simultaneous measurements under three different zenith angles are presented.

Kavlakov, S.↗

Evaluation of a strained-coordinate perturbation procedure - Nonlinear subsonic and transonic flows

An evaluation is made of a perturbation method devised to obtain highly accurate approximations to families of strongly nonlinear solutions which are either continuous or discontinuous, and which represent variations in some arbitrary parameter. The method first defines a unit perturbation by using two nonlinear solutions which differ from one another by a nominal change in some geometric or flow parameter, then employs this unit perturbation to predict a family of related nonlinear solutions over a range of parameter variation. Coordinate straining is incorporated into this perturbation method for determining the unit perturbation to account for the movement of discontinuities and maxima of high-gradient regions due to the perturbation. Attention is given to transonic and subsonic flows. Comparisons of the perturbation results with the corresponding 'exact' nonlinear solutions show a remarkable accuracy and range of validity of the perturbation method across the spectrum of examples considered.

Stahara, S. S.↗

Spatial decomposition of magnetic anisotropy in magnets: Application to doped Fe 16 N 2

We propose a scheme of decomposition of the total relativistic energy in solids to intra- and interatomic contributions. The method is based on a site variation of such fundamental constant as the speed of light. As a practical illustration of the method, we tested such decomposition in the case of a spin-orbit interaction variation for the decomposition of the magnetic anisotropy energy (MAE) in CoPt. We further studied the α" - Fe 16 N 2 magnet doped by Bi, Sb, Co, and Pt atoms. It was found that the addition of Pt atoms can enhance the MAE by as much as five times while Bi and Sb substitutions double the total MAE. Using the proposed technique, we demonstrate the spatial distribution of these enhancements. Our studies also suggest that Sb, Pt, and Co substitutions could be synthesized by experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Riley County Water Resources Project Summary - Comparing Runoff Curve Calculation Methods to Inform Local Resiliency Initiatives in Riley County, Kansas

Riley County, Kansas, has observed increased levels of flooding, potentially due to changes in land use/land cover (LULC) and seasonal vegetation variation. This study contrasts two methods of generating runoff curve numbers (CN) from 2006-2020. (1) The traditional Soil Conservation Service CN calculation method uses a look-up table and tracked LULC to determine runoff changes. These tables allow for land cover-specific CN and account for various farming techniques but lack flexibility in calculations for various seasons or plant health. (2) A dynamic method employs normalized difference vegetation index (NDVI) compiled over the rainy season each year to calculate CN using seasonal vegetation. This method allows for a more precise analysis of runoff variability within and between rainy seasons because it can be updated with greater temporal detail and captures higher spatial resolutions by using NDVI as a proxy for LULC. This study further uses inputs from the United States Geological Survey (USGS) National Land Cover Database (NLCD), the United States Department of Agriculture (USDA) Cropland Data Layer, and Landsat imagery to create more precise LULC raster datasets including both urban cover and crop-specific land use and curve number maps of the area. Results can guide decision makers in the City of Manhattan, Riley County Department of Planning and Development, Riley County Conservation District, the Kansas Forest Service, and the Kansas Department of Health and Environment toward informed decisions on resiliency strategies to address future flooding.

DEVELOP Project Summary↗

Diurnal variations of cosmic ray geomagnetic cut-off threshold rigidities

The spectrographic global survey method was used to investigate the rigidity variations Rc of geomagnetic cut-off as a function of local time and the level of geomagnetic disturbance for a number of stations of the world wide network. It is shown that geomagnetic cut-off threshold rigidities undergo diurnal variations. The diurnal wave amplitude decreases with increasing threshold rigidity Rc, and the wave maximum occurs at 2 to 4 hr LT. The amplitude of diurnal variations increases with increasing geomagnetic activity. The results agree with those from trajectory calculations made for an asymmetric model of the magnetosphere during different geomagnetic disturbance conditions.

Dvornikov, V. M.↗

Improved heavy-ion PID using scintillation light detector with neural network analysis: a Monte Carlo simulation study

The photon collection efficiency of gaseous scintillator detectors varies according to the position of the impinging charged particles in the medium that generates scintillation light. Thus, when impinging particles are distributed over a large area, the intrinsic photon-number resolution of the system is affected by a large variation. This work presents and discusses a method for adjusting the total number of detected photons to account for variation in the photon collection efficiency as a function of the position of the light source within the scintillating medium. The method was developed and validated by processing data from systematic simulation studies based on GEANT4 that model the response of the Energy Loss Optical Scintillation System (ELOSS) detector. The position of the charged particle is calculated using a deep neural network algorithm. This is accomplished by analyzing the distribution of scintillation light recorded by the array of photosensors. The estimated particle position is then used to calculate the correction factor and adjust the amount of captured light to account for variations in the photon collection efficiency. The neural network algorithm provides excellent tracking capabilities, achieving sub-millimeter position resolution and an angular resolution of 12 mrad, approaching the performance of traditional tracking detectors (e.g., drift chambers). The present method can be generalized to any optical scintillation system where the photon collection efficiency depends on the position of the impinging particle.

Heavy-ion detectors↗