Search NASA⌕ Search

SEARCH · Search NASA

Results for “Variational discretization”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Measurement of auroral rays from the space shuttle

Videotapes of auroras from two space shuttle missions were analyzed to determine the locations of auroral rays. One hundred and forty auroral ray locations were determined by triangulating from pairs of images taken several seconds (20-100 km) apart. The rays, observed over a range of 18 hours in local magnetic time between 62 deg and 85 deg magnetic latitude, were mostly confined to the Feldstein oval in the evening and midnight sectors but were uniformly scattered well into the polar cap in the morning sector. Plots of ray border altitude versus magnetic time show a steady increase in lower border altitude, suggesting a decrease in energy from near 10 keV in the early evening to near 2 keV by midmorning. This variation for the energy in discrete rayed arcs is opposite to previous reports for the average energy of auroral electrons and suggests that the acceleration mechanism within discrete arcs differs from that pertaining to the auroral oval in general. Comparison of upper and lower border altitudes suggests a nearly monoenergetic electron flux shortly before midnight and a more extended distribution at both earlier and later times. The apparently independent variation in maximum energy and energy spread is difficult to reconcile with models involving acceleration of auroral electrons by parallel electric fields. There was only a slight dependence of border altitudes on magnetic latitude within the auroral oval, while rays in the polar cap (northward B(sub z)) were significantly higher.

Jack, Todd M.↗
NASA NTRS2000

A national survey of policies on disclosure of conflicts of interest in biomedical research

BACKGROUND: Conflicts of interest pose a threat to the integrity of scientific research. The current regulations of the U.S. Public Health Service and the National Science Foundation require that medical schools and other research institutions report the existence of conflicts of interest to the funding agency but allow the institutions to manage conflicts internally. The regulations do not specify how to do so. METHODS: We surveyed all medical schools (127) and other research institutions (170) that received more than $5 million in total grants annually from the National Institutes of Health or the National Science Foundation; 48 journals in basic science and clinical medicine; and 17 federal agencies in order to analyze their policies on conflicts of interest. RESULTS: Of the 297 institutions, 250 (84 percent) responded by March 2000, as did 47 of the 48 journals and 16 of the 17 federal agencies. Fifteen of the 250 institutions (6 percent)--5 medical schools and 10 other research institutions--reported that they had no policy on conflicts of interest. Among the institutions that had policies, there was marked variation in the definition and management of conflicts. Ninety-one percent had policies that adhered to the federal threshold for disclosure ($10,000 in annual income or equity in a relevant company or 5 percent ownership), and 9 percent had policies that exceeded the federal guidelines. Only 8 percent had policies requiring disclosure to funding agencies, only 7 percent had such policies regarding journals, and only 1 percent had policies requiring the disclosure of information to the relevant institutional review boards or to research subjects. Twenty journals (43 percent) reported that they had policies requiring disclosure of conflicts of interest. Only four federal agencies had policies that explicitly addressed conflicts of interest in extramural research, and all but one of the agencies relied primarily on institutional discretion. CONCLUSIONS: There is substantial variation among policies on conflicts of interest at medical schools and other research institutions. This variation, combined with the fact that many scientific journals and funding agencies do not require disclosure of conflicts of interest, suggests that the current standards may not be adequate to maintain a high level of scientific integrity.

Human↗
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↗
DOE OSTI2024

Role of momentum in the generator-coordinate method applied to barrier penetration

Nuclear fission at barrier-top energies is conventionally modeled by a one-dimensional Schrödinger equation applied to internal fission channels, but that treatment is hard to justify in the configuration-interaction approach to nuclear Hamiltonians. Here we show that inclusion of states of finite momentum by the generator coordinate method (GCM) considerably extends the range of energies at which GCM-based Hamiltonians could reproduce the Schrödinger treatment. Furthermore, the transmission probabilities for crossing the barrier are calculated by a discrete version of Kohn's variational method, which may also be useful for other systems of interacting fermions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1982

Satellite observations of suprathermal electron bursts

Data from the ISIS 2 satellite are used to study the spatial distribution of the suprathermal electron bursts discovered by high-altitude sounding rockets. It is found that the bursts form an integral part of the boundary plasma sheet and are linked with the formation of discrete auroral arcs. The maximum energy of the bursts is always below the peak energy in neighboring discrete forms, and it appears that if the bursts fail to reach the minimum energy found to be necessary to sustain a discrete arc they collapse within a few seconds. It is suggested that in a pass across the auroral oval the overall envelope of the maximum energy in bursts and the peak energy in discrete auroral forms reflect the variation of the available acceleration (perhaps a parallel electric potential) across the auroral oval while the energy structure in individual forms is a function of the nonlinear processes which sustain them. This explanation appears to resolve the apparent discrepancy in scale thickness between rocket and satellite inverted V's.

Johnstone, A. D.↗
NASA NTRS1985

Finite difference operators on unstructured triangular meshes

A new form of the Laplace operator is derived which is shown to be second-order accurate when calculated at mesh triangle nodes if the cells corresponding to the nodes have a high degree of symmetry. It is demonstrated that the Laplace operator directly derived from Stokes' integral theorem is only zeroth-order accurate on these same cells. The construction of gradient and Laplace operators on unstructured triangular grids are discussed from the variational point of view. Two discrete representations of the Laplacian are compared to each other, and conclusions are drawn regarding their accuracy on a pointwise basis. Finally, geometric relations valid on arbitrary polygons are given for completeness in an appendix.

Erlebacher, G.↗
NASA NTRS1992

Conceptual analysis of a lunar base transportation system

Important to the planning for a lunar base is the development of transportation requirements for the establishment and maintenance of that base. This was accomplished as part of a lunar base systems assessment study conducted by the NASA Langley Research Center in conjunction with the NASA Johnson Space Center. Lunar base parameters are presented using a baseline lunar facility concept and timeline of developmental phases. Masses for habitation and scientific modules, power systems, life support systems, and thermal control systems were generated, assuming space station technology as a starting point. The masses were manifested by grouping various systems into cargo missions and interspersing manned flights consistent with construction and base maintenance timelines. A computer program that sizes the orbital transfer vehicles (OTV's), lunar landers, lunar ascenders, and the manned capsules was developed. This program consists of an interative technique to solve the rocket equation successively for each velocity correction (delta V) in a mission. The delta V values reflect integrated trajectory values and include gravity losses. As the program computed fuel masses, it matched structural masses from General Dynamics' modular space-based OTV design. Variables in the study included the operation mode (i.e., expendable vs. reusable and single-stage vs. two-stage OTV's), cryogenic specific impulse, reflecting different levels of engine technology, and aerobraking vs. all-propulsive return to Earth orbit. The use of lunar-derived oxygen was also examined for its general impact. For each combination of factors, the low-Earth orbit (LEO) stack masses and Earth-to-orbit (ETO) lift requirements are summarized by individual mission and totaled for the developmental phase. In addition to these discrete data, trends in the variation of study parameters are presented.

Hoy, Trevor D.↗
NASA NTRS2006

Isotope Variations in Terrestrial Carbonates and Thermal Springs as Biomarkers: Analogs for Martian Processes

Stable isotope measurements of carbonate minerals contained within ALH84001 [1] suggest that fluids were present at 3.9 Gy on Mars [2, 3, 4, 5]. Both oxygen and carbon isotopes provide independent means of deciphering paleoenvironmental conditions at the time of carbonate mineral precipitation. In terrestrial carbonate rocks oxygen isotopes not only indicate the paleotemperature of the precipitating fluid, but also provide clues to environmental conditions that affected the fluid chemistry. Carbon isotopes, on the other hand, can indicate the presence or absence of organic compounds during precipitation (i.e. biogenically vs. thermogenically-generated methane), thus serving as a potential biomarker. We have undertaken a study of micro scale stable isotope variations measured in some terrestrial carbonates and the influence of organic compounds associated with the formation of these carbonates. Preliminary results indicate that isotope variations occur within narrow and discrete intervals, providing clues to paleoenvironmental conditions that include both biological and non-biological activity. These results carry implications for deciphering Martian isotope data and therefore potential biological prospecting on the planet Mars. Recently, Fourier Transform Spectrometer observations have detected methane occurring in the Martian atmosphere [6] that could be attributed to a possible biogenic source. Indeed, Mars Express has detected the presence of methane in the Martian atmosphere [7], with evidence indicating that methane abundances are greatest above those basins with high water concentrations.

Socki, Richard A.↗
NASA NTRS2011

Covariant C and O Isotope Trends in Some Terrestrial Carbonates and ALH 84001: Possible Linkage Through Similar Formation Processes

Carbonate minerals found on the surface of Mars and in martian meteorites indicate that liquid water has played a significant role in the planet's history. These findings have raised questions regarding the history of the martian hydrosphere and atmosphere as well as the possibility of life. Sunset Crater, Arizona is a dry environment with relatively high evaporation and brief periods of precipitation. This environment resembles Mars and may make Sunset Crater a good analog to martian carbonates. In this study we sought to identify discrete micro-scale isotopic variation within the carbonate crusts in Sunset Crater to see if they resembled the micro-scale isotope variation found in ALH 84001 carbonates. Sunset Crater carbonate formation may be used as a martian analog and ultimately provide insight into carbonate formation in ALH 84001.

Volk, Kathryn E.↗
DOE OSTI2023

Linearization errors in discrete goal-oriented error estimation

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.

42 ENGINEERING↗
NASA NTRS1991

Cometary activity, discrete outgassing areas, and dust-jet formation

Conceptual models for various types of features observed in cometary comae (jets, spirals, halos, fans, etc.), their computer simulation, and the hydrodynamic models for jet formation are critically reviewed, and evidence for anisotropic, strongly collimated flows of ejecta emanating from discrete active regions (vents) on the rotating cometary nuclei is presented. Techniques employed to generate synthetic comet images that simulate the features observed are described, and their relevance to the primary objects of coma-morphology studies is discussed. Modeling of temporal variations in the water emission from discrete active regions suggests that production curves asymmetric with respect to perihelion should be commonplace. Critical comparisons with the activity profiles of Enke's comet and with light curves of disappearing comets and comets that undergo outbursts are presented. Recent developments in the understanding of the processes that cause the nongravitational perturbations of cometary motions are reviewed, and the observed discontinuities are identified with the birth of new sources and/or deactivation of old vents.

Sekanina, Z.↗
NASA NTRS1993

Variational approaches to CFD: Applications to potential Euler and Navier-Stokes equations

Flow simulations based on variational techniques are discussed. Variational techniques are applied for the problem formulation, for the construction of discrete schemes, and for the solution procedures. Inviscid transonic flows, as well as viscous incompressible and compressible flows, are analyzed by several methods. Variational techniques to solve the discrete systems of equations are well known. The most popular one is the conjugate gradient method. Extensions to nonsymmetric and nonpositive definite systems were under investigation for the last two decades. Variational techniques were used for convergence acceleration of iterative procedures. Some of these ideas based on recent work are examined. The lectures are organized in three parts including numerical simulation of inviscid as well as viscous incompressible and compressible flows.

Hafez, M.↗
NASA NTRS2001

Discrete-State Simulated Annealing For Traveling-Wave Tube Slow-Wave Circuit Optimization

Algorithms based on the global optimization technique of simulated annealing (SA) have proven useful in designing traveling-wave tube (TWT) slow-wave circuits for high RF power efficiency. The characteristic of SA that enables it to determine a globally optimized solution is its ability to accept non-improving moves in a controlled manner. In the initial stages of the optimization, the algorithm moves freely through configuration space, accepting most of the proposed designs. This freedom of movement allows non-intuitive designs to be explored rather than restricting the optimization to local improvement upon the initial configuration. As the optimization proceeds, the rate of acceptance of non-improving moves is gradually reduced until the algorithm converges to the optimized solution. The rate at which the freedom of movement is decreased is known as the annealing or cooling schedule of the SA algorithm. The main disadvantage of SA is that there is not a rigorous theoretical foundation for determining the parameters of the cooling schedule. The choice of these parameters is highly problem dependent and the designer needs to experiment in order to determine values that will provide a good optimization in a reasonable amount of computational time. This experimentation can absorb a large amount of time especially when the algorithm is being applied to a new type of design. In order to eliminate this disadvantage, a variation of SA known as discrete-state simulated annealing (DSSA), was recently developed. DSSA provides the theoretical foundation for a generic cooling schedule which is problem independent, Results of similar quality to SA can be obtained, but without the extra computational time required to tune the cooling parameters. Two algorithm variations based on DSSA were developed and programmed into a Microsoft Excel spreadsheet graphical user interface (GUI) to the two-dimensional nonlinear multisignal helix traveling-wave amplifier analysis program TWA3. The algorithms were used to optimize the computed RF efficiency of a TWT by determining the phase velocity profile of the slow-wave circuit. The mathematical theory and computational details of the DSSA algorithms will be presented and results will be compared to those obtained with a SA algorithm.

Wilson, Jeffrey D.↗
NASA NTRS1994

A Remarkable Auroral Event on Jupiter Observed in the Ultraviolet with the Hubble Space Telescope

Two sets of ultraviolet images of the Jovian north aurora were obtained with the Faint Object Camera on board the Hubble Space Telescope. The first series shows an intense discrete arc in near corotation with the planet. The maximum apparent molecular hydrogen emission rate corresponds to an electron precipitation of approximately 1 watt per square meter, which is about 30,000 times larger than the solar heating by extreme ultraviolet radiation. Such a particle heating rate of the auroral upper atmosphere of Jupiter should cause a large transient temperature increase and generate strong thermospheric winds. Twenty hours after initial observations, the discrete arc had decreased in brightness by more than one order of magnitude. The time scale and magnitude of the change in the ultraviolet aurora leads us to suggest that the discrete Jovian auroral precipitation is related to large-scale variations in the current system as is the case of for Earth's discrete aurorae.

Gerard, J. C.↗
DOE OSTI2023

Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces

Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used to generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, including many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process using exact (as opposed to variational) analysis. We relate the theory to the existing continuous-state Gaussian diffusion as well as other approaches to discrete diffusion, and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce “Blackout Diffusion”, which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10, Binarized MNIST, and CelebA datasets confirm the feasibility of our approach. Generalizing from specific (Gaussian) forward processes to discrete-state processes without a variational approximation sheds light on how to interpret diffusion models, which we discuss.

Santos, Javier E.↗
DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗
DOE OSTI2022

A Differentiable Hydrology Approach for Modeling With Time-Varying Parameters

In the analysis of observational data with mathematical models, it is often desirable to let one or more model parameters vary with time to account for changing environmental conditions or to offer more flexibility. A standard approach for dynamic models is to allow for the number of parameters to grow as a function of time, leading to a high-dimensional inference problem for long data records. We propose using Hamiltonian Monte Carlo, a gradient-based Markov chain Monte Carlo (MCMC) method, to perform Bayesian inference for time-varying parameters (TVP) in hydrology models. As derivatives of model error functions with regard to parameters are not available in closed form, we implemented the GR4J rainfall-runoff models in Theano/PyMC3, and Jax/NumPyro, allowing for the application of automatic differentiation techniques. A simulation experiment assessing the viability of these methods for recovering underlying temporal variation in parameters indicates that a discrete Gaussian random walk-based prior appears to be best suited across a range of scenarios while the continuous-time Gaussian process compares relatively poorly. In an analysis of 20 years of daily real-world streamflow records from the Model Parameter Estimation Experiment (MOPEX) hydrology data set with increasing data sparsity and at multiple temporal resolutions, we find that TVP inference with gradient-based MCMC is a flexible and reliable approach for analyzing hydrology models with dynamic parameter sets.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.

Gnoffo, Peter A.↗