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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Efficient implementation of essentially non-oscillatory shock capturing schemes

In the computation of discontinuous solutions of hyperbolic conservation laws, TVD (total-variation-diminishing), TVB (total-variation-bounded) and the recently developed ENO (essentially non-oscillatory) schemes have proven to be very useful. In this paper two improvements are discussed: a simple TVD Runge-Kutta type time discretization, and an ENO construction procedure based on fluxes rather than on cell averages. These improvements simplify the schemes considerably -- especially for multi-dimensional problems or problems with forcing terms. Preliminary numerical results are also given.

Shu, Chi-Wang↗
NASA NTRS1988

Efficient implementation of essentially non-oscillatory shock-capturing schemes

In the computation of discontinuous solutions of hyperbolic conservation laws, TVD (total-variation-diminishing), TVB (total-variation-bounded) and the recently developed ENO (essentially non-oscillatory) schemes have proven to be very useful. In this paper two improvements are discussed: a simple TVD Runge-Kutta type time discretization, and an ENO construction procedure based on fluxes rather than on cell averages. These improvements simplify the schemes considerably - especially for multi-dimensional problems or problems with forcing terms. Preliminary numerical results are also given.

Shu, Chi-Wang↗
NASA NTRS1989

Kilometric radiation power flux dependence on area of discrete aurora

Kilometer wavelength radiation, measured from distant positions over the North Pole and over the Earth's equator, was compared to the area of discrete aurora imaged by several low-altitude spacecraft. Through correlative studies of auroral kilometric radiation (AKR) with about two thousand auroral images, a stereoscopic view of the average auroral acceleration region was obtained. A major result is that the total AKR power increases as the area of the discrete auroral oval increases. The implications are that the regions of parallel potentials or the auroral plasma cavities, in which AKR is generated, must possess the following attributes: (1) they are shallow in altitude and their radial position depends on wavelength, (2) they thread flux tubes of small cross section, (3) the generation mechanism in them reaches a saturation limit rapidly, and (4) their distribution over the discrete auroral oval is nearly uniform. The above statistical results are true for large samples collected over a long period of time (about six months). In the short term, AKR frequently exhibits temporal variations with scales as short as three minutes (the resolution of the averaged data used). These fluctuations are explainable by rapid quenchings as well as fast starts of the electron cyclotron maser mechanism. There were times when AKR was present at substantial power levels while optical emissions were below instrument thresholds. A recent theoretical result may account for this set of observations by predicting that suprathermal electrons, of energies as low as several hundred eV, can generate second harmonic AKR. The indirect observations of second harmonic AKR require that these electrons have mirror points high above the atmosphere so as to minimize auroral light emissions. The results provide evidence supporting the electron cyclotron maser mechanism.

Saflekos, N. A.↗
NASA NTRS1989

Discrete cloud structure on Neptune

Recent CCD imaging data for the discrete cloud structure of Neptune shows that while cloud features at CH4-band wavelengths are manifest in the southern hemisphere, they have not been encountered in the northern hemisphere since 1986. A literature search has shown the reflected CH4-band light from the planet to have come from a single discrete feature at least twice in the last 10 years. Disk-integrated photometry derived from the imaging has demonstrated that a bright cloud feature was responsible for the observed 8900 A diurnal variation in 1986 and 1987.

Hammel, H. B.↗
NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗
NASA NTRS2019

A Study of the 20 Day Superorbital Modulation in the High-Mass X-Ray Binary IGR J16493-4348

We report on Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory(Swift) X-ray Telescope (XRT), and Swift Burst Alert Telescope (BAT) observations of IGR J16493-4348, a wind-fed supergiant X-ray binary showing significant superorbital variability. From a discrete Fourier transform of the BAT light curve,we refine its superorbital period to be 20.058±0.007 days. The BAT dynamic power spectrum and a fractional root mean square analysis both show strong variations in the amplitude of the superorbital modulation, but no observed changes in the period are found. The superorbital modulation is significantly weaker between MJD 55,700 and MJD 56,300. The joint NuSTAR and XRT observations, which were performed near the minimum and maximum of one cycle of the 20 day superorbital modulation, show that the flux increases by more than a factor of two between superorbital minimum and maximum. We find no significant changes in the 3-50 keV pulse profiles between superorbital minimum and maximum, which suggests a similar accretion regime. Modeling the pulse-phase-averaged spectra we find a possible Fe Kα emission line at 6.4 keV at superorbital maximum. This feature is not significant at superorbital minimum. While we do not observe any significant differences between the pulse-phase-averaged spectral continua apart from the overall flux change, we find that the hardness ratio near the broad main peak of the pulse profile increases from superorbital minimum to maximum. This suggests the spectral shape hardens with increasing luminosity. We discuss different mechanisms that might drive the observed superorbital modulation.

Coley, Joel B.↗
NASA NTRS1985

Cometary kilometric radio waves and plasma waves correlated with ion pick-up effect at Comet Halley

Bow-shock movements at Comet Halley are inferred from the discrete spectra of the cometary kilometric radiation (30-195 kHz); the observed emissions can be interpreted as being generated and propagating from the moving shock. The shock motion is possibly associated with the time variation of the solar wind and cometary outgassing. It is concluded that these plasma wave phenomena are manifestations of ion pick-up processes, which occur even in a remote region 7 million to 10 million km from the cometary nucleus.

Oya, H.↗
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.↗
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.↗
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.↗
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.↗