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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1992

BATSE - The burst and transient source experiment on the Gamma Ray Observatory

BATSE, one of the four instruments on the Gamma Ray Observatory, consists of eight detectors at the corners of the spacecraft. The large size of BATSE, its full-sky viewing and versatile data system should result in considerable advances in the study of gamma-ray bursts and other transient sources. The more significant features of BATSE include: the ability to obtain locations for weak bursts, the ability to detect spectral features and rapid spectral variations in many bursts and the ability to provide burst locations to the scientific community quicker than previously possible. The instrument also has considerable capability for the study of pulsars, solar flares, and other discrete sources at energies greater than 30 keV.

Fishman, G. J.↗

NASA NTRS2006

Accuracy Quantification of the Loci-CHEM Code for Chamber Wall Heat Fluxes in a G02/GH2 Single Element Injector Model Problem

A robust rocket engine combustor design and development process must include tools which can accurately predict the multi-dimensional thermal environments imposed on solid surfaces by the hot combustion products. Currently, empirical methods used in the design process are typically one dimensional and do not adequately account for the heat flux rise rate in the near-injector region of the chamber. Computational Fluid Dynamics holds promise to meet the design tool requirement, but requires accuracy quantification, or validation, before it can be confidently applied in the design process. This effort presents the beginning of such a validation process for the Loci- CHEM CPD code. The model problem examined here is a gaseous oxygen (GO2)/gaseous hydrogen (GH2) shear coaxial single element injector operating at a chamber pressure of 5.42 MPa. The GO2/GH2 propellant combination in this geometry represents one the simplest rocket model problems and is thus foundational to subsequent validation efforts for more complex injectors. Multiple steady state solutions have been produced with Loci-CHEM employing different hybrid grids and two-equation turbulence models. Iterative convergence for each solution is demonstrated via mass conservation, flow variable monitoring at discrete flow field locations as a function of solution iteration and overall residual performance. A baseline hybrid grid was used and then locally refined to demonstrate grid convergence. Solutions were also obtained with three variations of the k-omega turbulence model.

West, Jeff↗

NASA NTRS1985

On the existence and stability conditions for mixed-hybrid finite element solutions based on Reissner's variational principle

The extensions of Reissner's two-field (stress and displacement) principle to the cases wherein the displacement field is discontinuous and/or the stress field results in unreciprocated tractions, at a finite number of surfaces ('interelement boundaries') in a domain (as, for instance, when the domain is discretized into finite elements), is considered. The conditions for the existence, uniqueness, and stability of mixed-hybrid finite element solutions based on such discontinuous fields, are summarized. The reduction of these global conditions to local ('element') level, and the attendant conditions on the ranks of element matrices, are discussed. Two examples of stable, invariant, least-order elements - a four-node square planar element and an eight-node cubic element - are discussed in detail.

Karlovitz, L. A.↗

DOE OSTI2025

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

Discovery of cometary kilometric radiations and plasma waves at Comet Halley

The plasma-wave probe carried by the spacecraft Sakigake discovered discrete spectra of emissions from Comet Halley in the frequency range 30-195 kHz. The observed cometary kilometric radiation appears to come from moving shocks in the coma region which are possibly associated with temporal variations of the solar wind. Waves due to plasma instabilities associated with the pick-up of cometary ions by the solar wind were observed within a region almost 10 million km from the comet nucleus.

Oya, H.↗

NASA NTRS1988

Second order accurate finite difference approximations for the transonic small disturbance equation and the full potential equation

New shock-capturing finite difference approximations for solving two scalar conservation law nonlinear partial differential equations describing inviscid, isentropic, compressible flows of aerodynamics at transonic speeds are presented. A global linear stability theorem is applied to these schemes in order to derive a necessary and sufficient condition for the finite element method. A technique is proposed to render the described approximations total variation-stable by applying the flux limiters to the nonlinear terms of the difference equation dimension by dimension. An entropy theorem applying to the approximations is proved, and an implicit, forward Euler-type time discretization of the approximation is presented. Results of some numerical experiments using the approximations are reported.

Mostrel, M. M.↗

NASA NTRS1989

Significant variation from a 1/R-squared potential in the cumulative flux determined from the Giotto Comet Halley Dust Impact Experiment

The Dust Impact Detection System (DIDSY) aboard the Giotto spacecraft provided the information on the dust flux, mass spectrum, and cumulative mass distribution flux in the coma of Comet Halley. Analysis of discrete pulse height data of cometary particles for the mass range of particles between 4.0 x 10 to the -10th g and 6.0 x 10 to the -6th g registered by the Giotto DIDSY detectors 2, 3, and 4 has been completed, and a cumulative flux has been determined for this size range of particles. Inside the cometopause, anomalous peaks have been identified as deviation from a 1/R-squared curve in both pre- and postencounter measurements.

Alexander, W. M.↗

NASA NTRS1990

Some Aspects of Essentially Nonoscillatory (ENO) Formulations for the Euler Equations, Part 3

An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.

Chakravarthy, Sukumar R.↗

DOE OSTI2023

Unsupervised machine learning discovery of structural units and transformation pathways from imaging data

We show that unsupervised machine learning can be used to learn chemical transformation pathways from observational Scanning Transmission Electron Microscopy (STEM) data. To enable this analysis, we assumed the existence of atoms, a discreteness of atomic classes, and the presence of an explicit relationship between the observed STEM contrast and the presence of atomic units. With only these postulates, we developed a machine learning method leveraging a rotationally invariant variational autoencoder (VAE) that can identify the existing molecular fragments observed within a material. The approach encodes the information contained in STEM image sequences using a small number of latent variables, allowing the exploration of chemical transformation pathways by tracing the evolution of atoms in the latent space of the system. The results suggest that atomically resolved STEM data can be used to derive fundamental physical and chemical mechanisms involved, by providing encodings of the observed structures that act as bottom-up equivalents of structural order parameters. The approach also demonstrates the potential of variational (i.e., Bayesian) methods in the physical sciences and will stimulate the development of more sophisticated ways to encode physical constraints in the encoder–decoder architectures and generative physical laws and causal relationships in the latent space of VAEs.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1972

Application of transport techniques to the analysis of NERVA shadow shields

A radiation shield internal to the NERVA nuclear rocket reactor required to limit the neutron and photon radiation levels at critical components located external to the reactor was evaluated. Two significantly different shield mockups were analyzed: BATH, a composite mixture of boron carbide, aluminum and titanium hydride, and a borated steel-liquid hydrogen system. Based on the comparisons between experimental and calculated neutron and photon radiation levels, the following conclusions were noted: (1) The ability of two-dimensional discrete ordinates code to predict the radiation levels internal to and at the surface of the shield mockups was clearly demonstrated. (2) Internal to the BATH shield mockups, the one-dimensional technique predicted the axial variation of neutron fluxes and photon dose rates; however, the magnitude of the neutron fluxes was about a factor of 1.8 lower than the two-dimensional analysis and the photon dose rate was a factor of 1.3 lower.

Capo, M. A.↗

NASA NTRS1982

Hybrid state vector methods for structural dynamic and aeroelastic boundary value problems

A computational technique is developed that is suitable for performing preliminary design aeroelastic and structural dynamic analyses of large aspect ratio lifting surfaces. The method proves to be quite general and can be adapted to solving various two point boundary value problems. The solution method, which is applicable to both fixed and rotating wing configurations, is based upon a formulation of the structural equilibrium equations in terms of a hybrid state vector containing generalized force and displacement variables. A mixed variational formulation is presented that conveniently yields a useful form for these state vector differential equations. Solutions to these equations are obtained by employing an integrating matrix method. The application of an integrating matrix provides a discretization of the differential equations that only requires solutions of standard linear matrix systems. It is demonstrated that matrix partitioning can be used to reduce the order of the required solutions. Results are presented for several example problems in structural dynamics and aeroelasticity to verify the technique and to demonstrate its use. These problems examine various types of loading and boundary conditions and include aeroelastic analyses of lifting surfaces constructed from anisotropic composite materials.

Lehman, L. L.↗

NASA NTRS1990

Clutter filter design considerations for Airborne Doppler radar detection of windshear

The problem of clutter rejection when processing down-looking Doppler radar returns from a low altitude airborne platform is a paramount problem. With radar as a remote sensor for detecting and predicting windshear in the vicinity of an urban airport, dynamic range requirements can exceed 50 dB because of high clutter to signal ratios. This presentation describes signal processing considerations in the presence of distributed and/or discrete clutter interference. Previous analyses have considered conventional range cell processing of radar returns from a rigidly mounted radar platform using either the Fourier or the pulse-pair method to estimate average windspeed and windspeed variation within a cell. Clutter rejection has been based largely upon analyzing a particular environment in the vicinity of the radar and employing a variety of techniques to reduce interference effects including notch filtering, Fourier domain line editing, and use of clutter maps. For the airborne environment the clutter characteristics may be somewhat different. Conventional clutter rejection methods may have to be changed and new methods will probably be required to provide useful signal to noise ratios. Various considerations are described. A major thrust has been to evaluate the effect of clutter rejection filtering upon the ability to derive useful information from the post-filter radar data. This analysis software is briefly described. Finally, some ideas for future analysis are considered including the use of adaptive filtering for clutter rejection and the estimation of windspeed spatial gradient directly from radar returns as a means of reducing the effects of clutter on the determination of a windshear hazard.

Baxa, Ernest G., Jr.↗

DOE OSTI2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Deep groundwater↗

DOE OSTI2023

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗

DOE OSTI2023

High-Order Mesh Morphing for Boundary and Interface Fitting to Implicit Geometries

Here, we propose a method that morphs high-order meshes such that their boundaries and interfaces coincide/align with implicitly defined geometries. Our focus is particularly on the case when the target surface is prescribed as the zero isocontour of a smooth discrete function. Common examples of this scenario include using level set functions to represent material interfaces in multimaterial configurations, and evolving geometries in shape and topology optimization. The proposed method formulates the mesh optimization problem as a variational minimization of the sum of a chosen mesh-quality metric using the Target-Matrix Optimization Paradigm (TMOP) and a penalty term that weakly forces the selected faces of the mesh to align with the target surface. The distinct features of the method are use of a source mesh to represent the level set function with sufficient accuracy, and adaptive strategies for setting the penalization weight and selecting the faces of the mesh to be fit to the target isocontour of the level set field. We demonstrate that the proposed method is robust for generating boundary- and interface-fitted meshes for curvilinear domains using different element types in 2D and 3D.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Sampling inhomogeneous turbulent fields

The reconstruction of an inhomogeneous random process from a finite number of discrete samples can be performed in terms of the Karhunen-Loeve (KL) expansion for that process. The n(th) eigenfunction has n - 1 zero crossings which are the sampling points for the inhomogeneous process. The rapid variation of the KL eigenfunctions makes it unnecessary to have a high density of sampling (or grid points) near the wall. However, this result should not be construed to indicate that with spectral simulations significantly fewer grid points are required with the KL expansion as compared to other orthogonal expansions. Moin and Moser (1989) have shown that the advantage of the KL expansion over Chebychev expansion rapidly diminishes when high percentage (say 90 percent) energy recovery is demanded.

Adrian, R. J.↗

NASA NTRS1991

Analysis of the optimal laminated target made up of discrete set of materials

A new class of problems was analyzed to estimate an optimal structure of laminated targets fabricated from the specified set of homogeneous materials. An approximate description of the perforation process is based on the model of radial hole extension. The problem is solved by using the needle-type variation technique. The desired optimization conditions and quantitative/qualitative estimations of optimal targets were obtained and are discussed using specific examples.

Aptukov, Valery N.↗

DOE OSTI2023

Geothermal Play Fairway Analysis, Part 2: GIS methodology

Play Fairway Analysis (PFA) in geothermal exploration originates from a systematic methodology developed within the petroleum industry and is based on a geologic, geophysical, and hydrologic framework of identified geothermal systems. We tailored this methodology to study the geothermal resource potential of the Snake River Plain and surrounding region, but it can be adapted to other geothermal resource settings. We adapted the PFA approach to geothermal resource exploration by cataloging the critical elements controlling exploitable hydrothermal systems, establishing risk matrices that evaluate these elements in terms of both probability of success and level of knowledge, and building a code-based ‘processing model’ to process results. A geographic information system was used to compile a range of different data types, which we refer to as elements (e.g., faults, vents, heat flow, etc.), with distinct characteristics and measures of confidence. Discontinuous discrete data (points, lines, or polygons) for each element were transformed into continuous interpretive 2D grid surfaces called evidence layers. Because different data types have varying uncertainties, most evidence layers have an accompanying confidence layer which reflects spatial variations in these uncertainties. Risk layers, as defined here, are the product of evidence and confidence layers, and are the building blocks used to construct Common Risk Segment (CRS) maps for heat, permeability, and seal, using a weighted sum for permeability and heat, but a different approach with seal. CRS maps quantify the variable risk associated with each of these critical components. In a final step, the three CRS maps were combined into a Composite Common Risk Segment (CCRS) map, using a modified weighted sum, for results that reveal favorable areas for geothermal exploration. Additional maps are also presented that do not mix contributions from evidence and confidence (to allow an isolated view of evidence and confidence), as well as maps that calculate favorability using the product of components instead of a weighted sum (to highlight where all components are present). Our approach helped to identify areas of high geothermal favorability in the western and central Snake River Plain during the first phase of study and helped identify more precise local drilling targets during the second phase of work. By identifying favorable areas, this methodology can help to reduce uncertainty in geothermal energy exploration and development.

15 GEOTHERMAL ENERGY↗