Search NASA⌕ Search

SEARCH · Search NASA

Results for “consistent 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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Improving Kinetic Consistency Across Elastic Model Transitions With Quadratic Inequality Constrained Weighted Least Squares (WLSQI)

Flexible body modeling presents a large challenge in the development of simulations to aid in design of flight control systems for launch and landing vehicles. Typically, the flexible body model is not a single continuous model but are rather discrete sets of Linear Time Invariant (LTI) Finite Element Models (FEM) incremented by propellant levels. This introduces the problem of smoothly transitioning modal and physical states of the vehicle when switching from one FEM to the next. This paper introduces a new approach to optimally transition flexible body states with weighted least squares, building off previous methods.

GN&C↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗
NASA NTRS2001

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗
NASA NTRS1984

Discrete orthogonal function expansions for non-uniform grids using the fast Fourier transform

A technique for applying discrete Fourier series to infinite domains is presented. The technique uses mappings designed to minimize truncation error and can be applied to solve mixed initial boundary value problems among others. The method is alias-free and yields consistent differentiation and integration operators. The mapping-induced truncation error is explicitly expressible and small in nearly all cases of interest. The method is illustrated for three problems involving convection, diffusion, and vortex interaction.

Cain, A. B.↗
NASA NTRS1994

The stability of perfect elliptic disks. 1: The maximum streaming case

Self-consistent distribution functions are constructed for two-dimensional perfect elliptic disks (for which the potential is exactly integrable) in the limit of maximum streaming; these are tested for stability by N-body integration. To obtain a discrete representation for each model, simulated annealing is used to choose a set of orbits which sample the distribution function and reproduce the required density profile while carrying the greatest possible amount of angular momentum. A quiet start technique is developed to place particles on these orbits uniformly in action-angle space, making the initial conditions as smooth as possible. The roundest models exhibit spiral instabilities similar to those of cold axisymmetric disks; the most elongated models show bending instabilities like those seen in prolate systems. Between these extremes, there is a range of axial ratios 0.25 approximately less than b/a approximately less than 0.6 within which these models appear to be stable. All the methods developed in this investigation can easily be extended to integrable potentials in three dimensions.

Levine, Stephen E.↗
NASA NTRS2007

The Prelude to and Aftermath of the Giant Flare of 2004 December 27: Persistent and Pulsed X-ray Properties of SGR 1806-20 from 1993 to 2005

We report on the evolution of key spectral and temporal parameters of SGR 1806-20 prior to and following the highly energetic giant flare of 2004 December 27. Using RXTE, we track the pulse frequency of the SGR and find that the spin-down rate varied erratically in the months before and after the flare. Contrary to the giant flare in SGR 1900+14, we find no evidence for a discrete jump in spin frequency at the time of the December 27th flare (Absolute value of (Delta * upsilon/upsilon) < 5 10 (exp -6)). In the months surrounding the flare, we find a strong correlation between pulsed flux and torque consistent with the model for magnetar magnetosphere electrodynamics proposed. As with the flare in SGR 1900+14, the pulse morphology of SGR 1806-20 changes drastically following the flare. Using Chandra and other publicly available imaging X-ray detector observations, we construct a spectral history of SGR 1806-20 from 1993 to 2005. The usual magnetar persistent emission spectral model of a power law plus a blackbody provides an excellent fit to the data. We confirm the earlier finding of increasing spectral hardness of SGR 1806-20 between 1993 and 2004. However, our results indicate significant differences in the temporal evolution of the spectral hardening. Rather than a direct correlation between torque and spectral hardness, we find evidence for a sudden torque change that preceded a gradual hardening of the energy spectrum on a timescale of years. Interestingly, the spectral hardness, spin-down rate, phase-averaged flux, and pulsed flux of SGR 1806-20 all peak months before the flare epoch.

Woods, Peter M.↗
NASA NTRS2009

Electron Backscatter Diffraction (EBSD) Analysis and U-Pb Geochronology of the Oldest Lunar Zircon: Constraining Early Lunar Differentiation and Dating Impact-Related Deformation

The evolution of the early moon was dominated by two processes (i) crystallization of the Lunar Magma Ocean (LMO) and differentiation of potassium-rare earth element-phosphorous-rich residual magma reservoir referred to as KREEP, and (ii) an intense meteorite bombardment referred to as lunar cataclysm . The exact timing of these processes is disputed, and resolution relies on collection and interpretation of precise age data. This study examines the microstructure and geochronology of zircon from lunar impact breccias collected during the Apollo 17 mission. A large zircon clast within lunar breccia 72215,195 shows sector zoning in optical microscopy, cathodoluminescence (CL) imaging and Raman mapping, and indicates that it was a relict fragment of a much larger magmatic grain. Sensitive high resolution ion microprobe (SHRIMP) U-Pb analysis of the zircon shows that U and Th concentration correlate with sector zoning, with darkest CL domains corresponding with high-U and Th (approx.150 and approx.100 ppm respectively), and the brightest-CL sectors containing approx.30-50 ppm U and approx.10-20 ppm Th. This indicates that variations in optical CL and Raman properties correspond to differential accumulation of alpha-radiation damage in each sector. Electron backscatter diffraction (EBSD) mapping shows that the quality of electron backscatter patterns (band contrast) varies with sector zoning, with the poorest quality patterns obtained from high-U and Th, dark-CL zones. EBSD mapping also reveals a deformation microstructure that is cryptic in optical, CL and Raman imaging. Two orthogonal sets of straight discrete and gradational low-angle boundaries accommodate approx.12 misorientation across the grain. The deformation bands are parallel to the crystallographic {a}-planes of the zircon, have misorientation axes parallel to the c-axis, and are geometrically consistent with formation by dislocation creep associated with <100>{010} slip. The deformation bands are unlike curved morphology of crystal-plastic microstructures in tectonically deformed terrestrial zircon, and geometrically similar to dislocation microstructures reported in experimentally shocked zircon. We interpret these crystal-plastic deformation microstructures to have resulted from a significant impact, either directly from impact shock, or during ductile flow directly following the impact. The deformation bands appear to continue undeflected through the non-indexed, radiation-damaged areas of the grain, which suggests that the orientation variation predates any significant mechanical weakening from radiation damage in the grain, and therefore occurred early in its history.

Timms, Nick↗
DOE OSTI2024

Mechanical and biochemical recovery of landfill waste in an underserved community

Historically in the United States, waste collected for recycling has been sold and shipped to processors in China. In 2013 and 2018, China introduced the Green Fence and National Sword policies which restricts the import of contaminated materials and banned the import of many recyclables. The cost of recycling in the United States has increased following these policy changes, which has led to many communities reducing their recycling programs or halting them altogether. Rural and underserved communities that don’t have resources to afford sophisticated recycling programs have been heavily impacted. Previous work at INL demonstrated that MSW is a potentially viable feedstock for both biochemical and thermochemical conversion. The goal of this project is to assess preprocessing tools that can produce consistent feedstocks that meet conversion specifications, remove problematic contaminants, and reduce the amount of waste that is landfilled. Municipal solid waste was collected from an underserved community in southeast Idaho, contaminants were characterized, and mechanically separated into two discrete fractions. The unit operations identified during mechanical separation trials will be mobilized to on-site with a goal of 50% recovery of paper and plastic waste.

09 - BIOMASS FUELS↗
DOE OSTI2024

Utah FORGE: 2024 Discrete Fracture Network Model Data

The Utah FORGE 2024 Discrete Fracture Network (DFN) Model dataset provides a set of files representing discrete fracture network modeling for the FORGE site near Milford, Utah. The dataset includes four distinct DFN model file sets, each corresponding to different time frames and modeling approaches in 2024. These models characterize both natural and induced fractures in the geothermal reservoir, which consists of crystalline granitic and metamorphic rock approximately 8,000 feet below the ground surface. The dataset includes a reference DFN model from February 2024 that incorporates planar fractures and well trajectories, as well as upscaled permeability, porosity, compressibility, and storage values on specified grids. Additionally, there are models based on new microseismic (MEQ) data from May and July 2024, including fracture planes fitted to the latest MEQ catalog datasets, tensile fractures from hydraulic stimulation, and an alternative connected DFN for modeling purposes. Coordinate data is provided in both global and local frames, with detailed instructions on the transformations used to align with principal stress orientations. The dataset also includes notes and calculation files for estimating fracture sizes and differences between various fracture sets. There are subfolders for Global Coordinates and Local Coordinates. To move from the global to the local coordinate frame, fractures and wells were a) rotated 20 degrees counterclockwise looking down about the global point (335376.400482041, 4263189.99998761, 250.093546450195) to better align with the principal stresses; and b) translated by (-335408.68, -4263010.9, 1150). Upscaled permeability values using the _XYZ suffix show directions with respect to the global XYZ coordinate frame, while those using the _IJK suffix are aligned with local coordinate frame.

15 GEOTHERMAL ENERGY↗
NASA NTRS1975

Fully unsteady subsonic and supersonic potential aerodynamics for complex aircraft configurations for flutter applications

A general theory for study, oscillatory or fully unsteady potential compressible aerodynamics around complex configurations is presented. Using the finite-element method to discretize the space problem, one obtains a set of differential-delay equations in time relating the potential to its normal derivative which is expressed in terms of the generalized coordinates of the structure. For oscillatory flow, the motion consists of sinusoidal oscillations around a steady, subsonic or supersonic flow. For fully unsteady flow, the motion is assumed to consist of constant subsonic or supersonic speed for time t or = 0 and of small perturbations around the steady state for time t 0.

Tseng, K.↗
NASA NTRS1977

Differential altimetry for satellite orbit determination

Differential altimetry is concerned with the employment of differenced satellite altimeter measurements at orbit ground trace intersections. The employment of this procedure makes it possible to eliminate two of the major error sources found in direct altimetry. Previous applications have not included the appropriate dynamic constraints required to account for correlations due to satellite orbit motion. A description is given of an investigation in which these correlations are included. The methodology produced is consistent with the dynamic environment. The regional or local limitations of previous approaches are overcome by extending the technique to the global scale. Attention is given to the description of the data type, the geometric topography height, altimeter errors, discretization errors, an approximate orbit determination problem, and a comparison of differenced altimeter measurements for retrograde and prograde orbits.

Hagar, H., Jr.↗
NASA NTRS1981

Atmospheric emissions photometric imaging experiment /AEPT/ for Spacelab 1

The atmospheric emissions photometric imaging experiment (AEPI) to be flown on Spacelab 1 is designed to study faint natural and artificial atmospheric emission phenomena. Optical emissions are imaged in the region 2150 A to 7320 A using a television system consisting of two optical channels, one wide-angle and one telephoto. The detection system is an image-enhanced SEC vidicon. A third optical channel images onto the photocathode of a microchannel plate photomultiplier tube that has 100 discrete anodes. Photons are counted for each discrete anode, providing a direct measure of the luminosity of an object viewed by the TV telephoto lens, albeit with low spatial resolution. The AEPI detector is mounted on a two-axis gimbal comprised of a Modified Apollo Telescope Mount Star Tracker (MAST), which provides experiment pointing over a 40-deg x 80-deg range, exclusive of restrictions due to the proximity of other experiments. The pointing stability is 1 arcmin with respect to the spacecraft coordinate system for an exposure of 1 second. The tracking capability is 3.5 deg/s with a stability of 1 arcmin. The detector and pointing system are located on the Spacelab pallet. The experiment is controlled by stored programs resident in the Dedicated Experiment Processor located in the Spacelab module.

Sandie, W. G.↗
NASA NTRS1991

Discrete Kalman filtering equations of second-order form for control-structure interaction simulations

A second-order form of discrete Kalman filtering equations is proposed as a candidate state estimator for efficient simulations of control-structure interactions in coupled physical coordinate configurations as opposed to decoupled modal coordinates. The resulting matrix equation of the present state estimator consists of the same symmetric, sparse N x N coupled matrices of the governing structural dynamics equations as opposed to unsymmetric 2N x 2N state space-based estimators. Thus, in addition to substantial computational efficiency improvement, the present estimator can be applied to control-structure design optimization for which the physical coordinates associated with the mass, damping and stiffness matrices of the structure are needed instead of modal coordinates.

Park, K. C.↗
NASA NTRS1995

Evaluation of Aircraft Platforms for SOFIA by Computational Fluid Dynamics

The selection of an airborne platform for the Stratospheric Observatory for Infrared Astronomy (SOFIA) is based not only on economic cost, but technical criteria, as well. Technical issues include aircraft fatigue, resonant characteristics of the cavity-port shear layer, aircraft stability, the drag penalty of the open telescope bay, and telescope performance. Recently, two versions of the Boeing 747 aircraft, viz., the -SP and -200 configurations, were evaluated by computational fluid dynamics (CFD) for their suitability as SOFIA platforms. In each configuration the telescope was mounted behind the wings in an open bay with nearly circular aperture. The geometry of the cavity, cavity aperture, and telescope was identical in both platforms. The aperture was located on the port side of the aircraft and the elevation angle of the telescope, measured with respect to the vertical axis, was 500. The unsteady, viscous, three-dimensional, aerodynamic and acoustic flow fields in the vicinity of SOFIA were simulated by an implicit, finite-difference Navier-Stokes flow solver (OVERFLOW) on a Chimera, overset grid system. The computational domain was discretized by structured grids. Computations were performed at wind-tunnel and flight Reynolds numbers corresponding to one free-stream flow condition (M = 0.85, angle of attack alpha = 2.50, and sideslip angle beta = 0 degrees). The computational domains consisted of twenty-nine(29) overset grids in the wind-tunnel simulations and forty-five(45) grids in the simulations run at cruise flight conditions. The maximum number of grid points in the simulations was approximately 4 x 10(exp 6). Issues considered in the evaluation study included analysis of the unsteady flow field in the cavity, the influence of the cavity on the flow across empennage surfaces, the drag penalty caused by the open telescope bay, and the noise radiating from cavity surfaces and the cavity-port shear layer. Wind-tunnel data were also available to compare to the CFD results; the data permitted an assessment of CFD as a design tool for the SOFIA program.

Klotz, S. P.↗
NASA NTRS1998

A Novel Four-Node Quadrilateral Smoothing Element for Stress Enhancement and Error Estimation

A four-node, quadrilateral smoothing element is developed based upon a penalized-discrete-least-squares variational formulation. The smoothing methodology recovers C1-continuous stresses, thus enabling effective a posteriori error estimation and automatic adaptive mesh refinement. The element formulation is originated with a five-node macro-element configuration consisting of four triangular anisoparametric smoothing elements in a cross-diagonal pattern. This element pattern enables a convenient closed-form solution for the degrees of freedom of the interior node, resulting from enforcing explicitly a set of natural edge-wise penalty constraints. The degree-of-freedom reduction scheme leads to a very efficient formulation of a four-node quadrilateral smoothing element without any compromise in robustness and accuracy of the smoothing analysis. The application examples include stress recovery and error estimation in adaptive mesh refinement solutions for an elasticity problem and an aerospace structural component.

Tessler, A.↗
NASA NTRS2005

Airborne Sunphotometer and Solar Spectral Flux Radiometer Measurements During INTEX/ITCT 2004

During the period 12 July - 8 August 2004, the NASA Ames 14-channel Airborne Tracking Sunphotometer (MTS-14) and Solar Spectral Flux Radiometer (SSFR) were operated aboard a Jetstream 31 (J31) aircraft and acquired measurements during nineteen science flights (approx. 53 flight hours) over the Gulf of Maine in support of the INTEX-NA (INtercontinental chemical Transport Experiment-North America) and ITCT (Intercontinental Transport and Chemical Transformation of anthropogenic pollution) field studies. In this paper, we will present results from analyses of those data sets. AATS-14 measures the direct solar beam transmission at fourteen discrete wavelengths (354-2138 nm), and provides instantaneous measurements of aerosol optical depth (AOD) spectra and water vapor column content, in addition to vertical profiles of aerosol extinction and water vapor density during suitable aircraft ascents and descents. SSFR consists of separate nadir and zenith viewing hemispheric FOV sensors that yield measurements of up- and downwelling solar irradiance at a spectral resolution of approx. 8-12 nm over the wavelength range 300-1700 nm. The objectives of the J31-based measurements during INTEX/ITCT were to provide AOD data for the evaluation of MODIS (MODerate-resolution Imaging Spectroradiometer) and MISR (Multi-angle Imaging Spectro-Radiometer) AOD retrievals, quantify sea surface spectral albedo (which can contribute the largest uncertainty to satellite aerosol retrievals for low aerosol loading), test closure (consistency) among suborbital results, test chemical-transport models using AOD profiles, and assess regional radiative forcing by combining satellite and suborbital results. Specific J31 flight patterns were designed to achieve these objectives, and they included a mixture of vertical profiles (spiral and ramped ascents and descents) and constant altitude horizontal transects at a variety of altitudes. Additional information is included in the original extended abstract.

Russell, P. B.↗
DOE OSTI2025

Ising on $\mathbb{S}^2$-- The Affine Conjecture

We review the recent construction \cite{brower2024isingmodelmathbbs2} of the 2d Ising model on a triangulated sphere $\mathbb{S}^2$. Surprisingly, this led to a precise map of the lattice couplings to the target geometry in order to reach the conform field theory (CFT) in the continuum limit. For the integrable 2d Ising CFT, the map was found analytically \cite{Brower_2023}. Here we conjecture how this might be generalized. The discrete geometry is implemented by the piecewise flat triangulation introduced by Regge in 1960 for the Einstein Hilbert action \cite{Regge1961GeneralRW}. Then following our Ising example, we posit the existence of a smooth map of lattice couplings in affine parameters consistent with quantum correlators. A sequence of theoretical investigations and numerical simulations are recommended to test this conjecture. They begin with non-integrable CFT's -- the 2d $\phi^4$ theory on $\mathbb{S}^2$; the 3d Ising model on $\mathbb{S}^3$ and $\mathbb{R} \times \mathbb{S}^2$; QED3 on $\mathbb{R} \times \mathbb{S}^{2} $ as an intermediate step to 4d non-Abelian lattice gauge theory on $\mathbb{R} \times \mathbb{S}^3$.

Brower, Richard C. [Boston U.] (ORCID:000000028185↗
NASA NTRS1986

Analysis of mixed-mode crack propagation using the boundary integral method

Crack propagation in a rotating inner raceway of a high speed roller bearing is analyzed using the boundary integral equation method. The method consists of an edge crack in a plate under tension, upon which varying Hertzian stress fields are superimposed. A computer program for the boundary integral equation method was written using quadratic elements to determine the stress and displacement fields for discrete roller positions. Mode I and Mode II stress intensity factors and crack extension forces G sub 00 (energy release rate due to tensile opening mode) and G sub r0 (energy release rate due to shear displacement mode) were computed. These calculations permit determination of that crack growth angle for which the change in the crack extension forces is maximum. The crack driving force was found to be the alternating mixed-mode loading that occurs with each passage of the most heavily loaded roller. The crack is predicted to propagate in a step-like fashion alternating between radial and inclined segments, and this pattern was observed experimentally. The maximum changes DeltaG sub 00 and DeltaG sub r0 of the crack extension forces are found to be good measures of the crack propagation rate and direction.

Mendelson, A.↗