Search NASA⌕ Search

SEARCH · Search NASA

Results for “variational method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,783 records · Page 99
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Scalar conservation and boundedness in simulations of compressible flow

With the proper combination of high-order, low-dissipation numerical methods, physics-based subgrid-scale models, and boundary conditions it is becoming possible to simulate many combustion flows at relevant conditions. However, non-premixed flows are a particular challenge because the thickness of the fuel/oxidizer interface scales inversely with Reynolds number. Sharp interfaces can also be present in the initial or boundary conditions. When higher-order numerical methods are used, there are often aphysical undershoots and overshoots in the scalar variables (e.g.passive scalars, species mass fractions or progress variable). These numerical issues are especially prominent when low-dissipation methods are used, since sharp jumps in flow variables are not always coincident with regions of strong variation in the scalar fields: consequently, special detection mechanisms and dissipative fluxes are needed. Most numerical methods diffuse the interface, resulting in artificial mixing and spurious reactions. In this paper, we propose a numerical method that mitigates this issue. As a result, we present methods for passive and active scalars, and demonstrate their effectiveness with several examples.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Characterization of Deficiencies in the Frequency Domain Forced Response Analysis Technique for Supersonic Turbine Bladed Disks

Turbine blades in rocket and jet engine turbomachinery experience enormous harmonic loading conditions. These loads result from the integer number of upstream and downstream stator vanes as well as the other turbine stages. Assessing the blade structural integrity is a complex task requiring an initial characterization of whether resonance is possible and then performing a forced response analysis if that condition is met. The standard technique for forced response analysis in rocket engine turbines is to decompose a computational fluid dynamics (CFD).generated flow field into its harmonic components, and to then perform a frequency response analysis at the problematic natural frequencies using cyclically symmetric structural dynamic models. Recent CFD analysis and water-flow testing at NASA/MSFC, though, indicates that this technique may miss substantial harmonic and non ]harmonic excitation sources that become present in complex flows. This complex content can only be captured by a CFD flow field encompassing at least an entire revolution. A substantial development effort to create a series of software programs to enable application of the 360 degree forcing function in a frequency response analysis on cyclic symmetric models has been completed (to be described in a future paper), but the question still remains whether the frequency response analysis itself is capable of capturing the excitation content sufficiently. Two studies comparing frequency response analysis with transient response analysis, therefore, of bladed-disks undergoing this complex flow environment have been performed. The first is of a bladed disk with each blade modeled by simple beam elements and the disk modeled with plates (using the finite element code MSC/NASTRAN). The focus of this model is to be representative of response of realistic bladed disks, and so the dimensions are roughly equivalent to the new J2X rocket engine 1st stage fuel pump turbine. The simplicity of the model allows the CFD load to be able to be readily applied, along with analytical and experimental variations in both the temporal and spatial fourier components of the excitation. In addition, this model is a first step in identifying response differences between transient and frequency forced response analysis techniques. The second phase assesses this difference for a much more realistic solid model of a bladed-disk in order to evaluate the effect of the spatial variation in loading on blade dominated modes. Neither research on the accuracy of the frequency response method when used in this context or a comprehensive study of the effect of test-observed variation on blade forced response have been found in the literature, so this research is a new contribution to practical structural dynamic analysis of gas turbines. The primary excitation of the upstream nozzles interacts with the blades on fuel pump of the J2X causes the 5th Nodal diameter modes to be excited, as explained by Tyler and Sofrin1, so a modal analysis was first performed on the beam/plate model and the 5ND bladed-disk mode at 40167 hz was identified and chosen to be the one excited at resonance (see figure 1). The first forced response analysis with this model focuses on identifying differences between frequency and transient response analyses. A hypothesis going into the analysis was that perhaps the frequency response was enforcing a temporal periodicity that did not really exist, and so therefore it would overestimate the response. As high dynamic response was a considerable source of stress in the J2X, examining this concept could potentially be beneficial for the program.

Brown, Andrew M.↗
NASA NTRS2020

Multi Model Monte Carlo with Python (MXMCPy)

Multi Model Monte Carlo with Python (\mxmc {}) is a software package developed as a general capability for computing the statistics of outputs from an expensive, high-fidelity model by leveraging faster, low-fidelity models for speedup. Motivated by uncertainty propagation problems where classical Monte Carlo (MC) simulation is computationally intractable, various multi-model MC approaches have recently emerged that yield unbiased estimators with significantly reduced variance relative to MC for the same cost. These existing methods include multi-level Monte Carlo (MLMC), multi-fidelity Monte Carlo (MFMC), and approximate control variates (ACV). Given a fixed computational budget and a collection of models with varying cost/accuracy, each method seeks a sample allocation strategy across the models that results in an estimator with optimal variance reduction. \mxmc {} is a versatile tool that enables convenient access to many existing multi-model MC approaches within one modular and extensible package. With \mxmc {}, users can easily compare existing methods to determine the best choice for their particular problem, while developers have a basis for implementing and sharing new variance reduction approaches. This report introduces the \mxmc {} software, providing a summary of the problem-solving workflow for users as well as a brief overview of the code layout for developers.

Geoffrey F Bomarito↗
DOE OSTI2019

Comparison of biofuel life-cycle GHG emissions assessment tools: The case studies of ethanol produced from sugarcane, corn, and wheat

The use of alternative fuels, particularly bio-based fuels, has been an important strategy to achieve greenhouse gas (GHG) emission reductions compared to petroleum-based fuels. However, discrepancies between results obtained by using different attributional life-cycle assessment (LCA) tools have challenged the credibility of the individual assessments, and as result, the progress towards or compliance with GHG mitigation targets. The objective of this study was to identify the main differences and commonalities in methodological structures, calculation procedures, and assumptions for the major commercial biofuel, ethanol, across three public LCA tools, BioGrace (EU), GHGenius (Canada), and GREET (U.S.), and a research-oriented fourth, the Virtual Sugarcane Biorefinery (VSB), a Brazilian platform for sugarcane ethanol assessments. The calculated emissions across models ranged from 16 to 45 for sugarcane, 43-62 for corn, and 45-68 g CO 2 eq MJ -1 for wheat ethanol. Harmonizing the three public models with VSB assumptions for sugarcane ethanol produced in Brazil, the range was reduced to 16-17 g CO 2 eq MJ -1 for sugarcane ethanol. Agricultural production (e.g., N 2 O emissions from fertilizers; energy and fuel use; straw field-burning; and limestone application) and ethanol shipping were found to be the major causes for variations for differences calculated for sugarcane ethanol. Similarly, harmonizing BioGrace and GHGenius calculations using GREET assumptions for U.S. corn ethanol generated nearly identical results (models varied within a 3% range). The coproduct treatment method was found to be the most influential parameter in the variations calculated for both corn and wheat ethanol. The application of the tools as part of GHG emissions accounting requirements is often defined via regulations and differences and/or conflicting assumptions set-forth in these models lead to most differences observed. Our study provides recommendations for promoting transparency in LCA calculations and assumptions across the tools used in research and development or for regulatory tools regarding biofuels.

09 BIOMASS FUELS↗
NASA NTRS1982

On the detection of other planetary systems by astrometric techniques

A quantitative method for astrometrically detecting perturbations induced in a star's motion by the presence of a planetary object is described. A periodogram is defined, wherein signals observed from a star show exactly periodic variations, which can be extracted from observational data using purely statistical methods. A detection threshold is defined for the frequency of occurrence of some detectable signal, e.g., the Nyquist frequency. Possible effects of a stellar orbital eccentricity and multiple companions are discussed, noting that assumption of a circular orbit assures the spectral purity of the signal described. The periodogram technique was applied to 12 yr of astrometric data from the U.S. Naval Observatory for three stars with low mass stellar companions. Periodic perturbations were confirmed. A comparison of the accuracy of different astrometric systems shows that the detection accuracy of a system is determined by the measurement accuracy and the number of observations, although the detection efficiency can be maximized by minimizing the number of data points for the case when observational errors are proportional to the square root of the number of data points. It is suggested that a space-based astrometric telescope is best suited to take advantage of the method.

Black, D. C.↗
NASA NTRS1987

The regression intersection method of adjusting image data for band ratioing

Estimation of the combined path radiance and sensor offset terms is essential in adjusting multispectral image radiance measurements for band ratioing. Commonly applied techniques for making this estimate have required assumptions or ancillary information regarding the reflectance properties of surface materials. This paper presents a technique that is unique in that it provides absolute (not just relative) statistically-derived estimates without the use of ancillary information. It is termed the regression intersection method (RIM). RIM is based on contrasts between the spectral properties of various homogeneous areas in rugged terrain. These areas are selected by examination of the image data alone. Bispectral regression (first principal component) lines are determined for each area and are projected, in pairs, to intersection points. Ideally, the coordinates of these points must equal the measurements for zero ground radiance since that is the only condition under which spectrally different materials can have the same radiance values. The median result from several site-pair and band-pair comparisons is used in order to statistically mitigate noise and minor variations due to natural variability. Tests show that the method is successful in determining correction values for the image data that result in maximum removal of the topographic effect in ratio images.

Crippen, Robert E.↗
DOE OSTI2019

Simulation of the Formation of a Cascade of Displacements and Transient Ionization Processes in Silicon Semiconductor Structures under Neutron Exposure

The formation of a disordered defect region in bulk silicon is simulated using the molecular-dynamics method for various energies of a primary recoil atom. Variations in the volume and number of radiation-induced defects in a cluster during its formation are calculated. The generation rates of nonequilibrium carriers and amplitude-temporal dependences of pulses of ionization currents in test Schottky diodes with hyperhigh frequencies are found theoretically.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Time-lapse electrical resistivity and induced polarization monitoring of a simulated ecosystem-scale coastal flooding experiment

Accurate estimation of the spatio-temporal variation in soil saturation and salt concentration is valuable for calibrating and evaluating Earth system models as well as detecting important eco-hydrological interactions. Such representative models play a vital role in investigating the stability of ecosystems facing changing hydrologic disturbance regimes. Ecosystem-scale field experiments are a useful way to capture hydrological disturbance and gain a deeper understanding of the complex mechanisms driving environmental changes. However, previous studies have not explored the use of quantitative imaging at a large spatial scale to reproduce breakthrough curves and inform system response and recovery following hydrologic disturbance events. Therefore, in this study, non-invasive geophysical methods were used to capture the spatio-temporal variation in subsurface saturation and salt concentrations during an ecosystem-scale coastal forest flooding experiment.

coastal↗
NASA NTRS1977

A first-principles derivation of Doppler noise expected from solar wind density fluctuations

The level of Doppler noise (DN) expected from solar wind (SW) density fluctuations (DF) is derived beginning with the expression for refractive index variations. The calculation takes account of up and downlink paths and of the method actually used to produce the DN values. The usual assumptions that the DF are frozen in, that the large-scale radial variation can be separated from the DF, that the DF power spectrum is a power law with outer scale k sub 0 and that the DF are homogeneous on scales less than 2c delta t, delta t = sample time, are made. The result agrees quite well with previous observations of DN. Corrections for the finite number of points used in the actual algorithm are discussed.

Callahan, P. S.↗
NASA NTRS1994

Equilibrium, stability, and orbital evolution of close binary systems

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence star binaries. We also discuss the orbital evoltuion of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

Lai, Dong↗
NASA NTRS1991

Advanced stress analysis methods applicable to turbine engine structures

The following tasks on the study of advanced stress analysis methods applicable to turbine engine structures are described: (1) constructions of special elements which contain traction-free circular boundaries; (2) formulation of new version of mixed variational principles and new version of hybrid stress elements; (3) establishment of methods for suppression of kinematic deformation modes; (4) construction of semiLoof plate and shell elements by assumed stress hybrid method; and (5) elastic-plastic analysis by viscoplasticity theory using the mechanical subelement model.

Pian, Theodore H. H.↗
NASA NTRS1993

An analytical optimization of electric propulsion orbit transfer vehicles

Due to the electric propulsion's inherent propellant mass savings over chemical propulsion, electric propulsion orbit transfer vehicles (EPOTV's) are highly efficient mode of orbit transfer. When selecting an electric propulsion device (ion, MPD, or arcjet) and propellant for a particular mission, it is preferable to use quick, analytical system optimization methods instead of time intensive numerical integration methods. It is also of interest to determine each thruster's optimal operating characteristics for a specific mission. Analytical expressions are derived which determine the optimal specific impulse (Isp) for each type of electric thruster to maximize payload fraction for a desired thrusting time. These expressions take into account the variation of thruster efficiency with specific impulse. Verification of the method is made with representative electric propulsion values on a LEO-to-GEO mission. Application of the method to specific missions is discussed.

Oleson, Steven R.↗
DOE OSTI2025

A critical review of the history of fabricating monolithic U-Mo fuel plates

The fabrication of monolithic uranium-molybdenum alloy fuels, specifically those developed for high-performance research and test reactors, began in the early 2000s. The primary fuel form consists of uranium alloyed with 10-wt% molybdenum in a thin foil coated with zirconium and encapsulated in aluminum-6061 cladding. Over the years, the process has evolved with different types of casting, heat treatments, rolling schedules, and cladding applications. This review examines the history of these fabrication processes and their impact on microstructure and fuel-swelling performance. Even though various fabrication methods were used, we found little correlation between fabrication variation and fuel swelling. This insensitive relationship between fabrication variation and fuel swelling is primarily due to inhomogeneous microstructures that formed during casting and grain refinement that occurred during rolling. We conclude that the fabrication processes we examined produced similar microstructures, indicating that the fuel microstructure is somewhat insensitive to the fabrication parameters evaluated. However, the relatively small amount of historic data, such as those for grain sizes, limited this analysis. More recently fabricated materials, such as those from ongoing irradiation experiment, Mini-Plate-1 and Mini-Plate-2, were also excluded from this analysis and are intended to be reviewed separately. The findings, that fuel microstructure is somewhat insensitive to the fabrication parameters, do not imply that any fabrication method is acceptable, given the uncertainties in data and fuel-swelling observations. For example, only arc melting and vacuum induction melting casting processes were previously explored in the historic fabrication efforts. In conclusion, the findings should not be extrapolated to other casting processes.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

An Electrical Resistance Diagnostic for Conductivity Monitoring in Laser Powder Bed Fusion

With the growing interest in metal additive manufacturing using laser powder bed fusion (LPBF), there is a need for advanced in-situ nondestructive evaluation (NDE) methods that can dynamically monitor manufacturing process-related variations, that can be used as a feedback mechanism to further improve the manufacturing process, leading to parts with improved microstructural properties and mechanical properties. Current NDE techniques either lack sensitivity beyond build layer, are costly or time-consuming, or are not compatible for in-situ integration. In this research, we develop an electrical resistance diagnostic for in-situ monitoring of powder fused regions during laser powder bed fusion printing. The technique relies on injecting current into the build plate and detecting voltage differences from conductive variations during printing using a simple, cheap four-point electrode array directly connected to the build plate. A computational model will be utilized to determine sensitivities of the approach, and preliminary experiments will be performed during the printing process to test the overall approach.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prime factorization using quantum variational imaginary time evolution

The road to computing on quantum devices has been accelerated by the promises that come from using Shor’s algorithm to reduce the complexity of prime factorization. However, this promise hast not yet been realized due to noisy qubits and lack of robust error correction schemes. Here we explore a promising, alternative method for prime factorization that uses well-established techniques from variational imaginary time evolution. We create a Hamiltonian whose ground state encodes the solution to the problem and use variational techniques to evolve a state iteratively towards these prime factors. We show that the number of circuits evaluated in each iteration scales as \(O(n^{5}d)\) , where n is the bit-length of the number to be factorized and d is the depth of the circuit. We use a single layer of entangling gates to factorize 36 numbers represented using 7, 8, and 9-qubit Hamiltonians. We also verify the method’s performance by implementing it on the IBMQ Lima hardware to factorize 55, 65, 77 and 91 which are greater than the largest number (21) to have been factorized on IBMQ hardware.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Evaluation of the pulse-contour method of determining stroke volume in man.

The pulse-contour method for determining stroke volume has been employed as a continuous rapid method of monitoring the cardiovascular status of patients. Twenty-one patients with ischemic heart disease and 21 patients with mitral valve disease were subjected to a variety of hemodynamic interventions. The pulse-contour estimations, using three different formulas derived by Warner, Kouchoukos, and Herd, were compared with indicator-dilution outputs. A comparison of the results of the two methods for determining stroke volume yielded correlation coefficients ranging from 0.59 to 0.84. The better performing Warner formula yielded a coefficient of variation of about 20%. The type of hemodynamic interventions employed did not significantly affect the results using the pulse-contour method. Although the correlation of the pulse-contour and indicator-dilution stroke volumes is high, the coefficient of variation is such that small changes in stroke volume cannot be accurately assessed by the pulse-contour method. However, the simplicity and rapidity of this method compared to determination of cardiac output by Fick or indicator-dilution methods makes it a potentially useful adjunct for monitoring critically ill patients.

Alderman, E. L.↗
DOE OSTI2025

Exploiting 20 Ne Isotopes for Precision Characterizations of Collectivity in Small Systems

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174⁢(8) stat ⁢(31) syst for NLEFT and 1.139⁢(6) stat ⁢(39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Visualisation and outlier detection for probability density function ensembles

Abstract Exploratory data analysis (EDA) for functional data—data objects where observations are entire functions—is a difficult problem that has seen significant attention in recent literature. This surge in interest is motivated by the ubiquitous nature of functional data, which are prevalent in applications across fields such as meteorology, biology, medicine and engineering. Empirical probability density functions (PDFs) can be viewed as constrained functional data objects that must integrate to one and be nonnegative. They show up in contexts such as yearly income distributions, zooplankton size structure in oceanography and in connectivity patterns in the brain, among others. While PDF data are certainly common in modern research, little attention has been given to EDA specifically for PDFs. In this paper, we extend several methods for EDA on functional data for PDFs and compare them on simulated data that exhibit different types of variation, designed to mimic that seen in real‐world applications. We then use our new methods to perform EDA on the breakthrough curves observed in gas transport simulations for underground fracture networks.

97 MATHEMATICS AND COMPUTING↗