Search NASA⌕ Search

SEARCH · Search NASA

Results for “local solution bounds”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Acoustic wave propagation in the solar atmosphere 1. Rediscussion of the linearized theory including nonstationary solutions

The normal dispersion analysis for linear adiabatic wave propagation in stratified atmospheres adopts a real frequency and solves for the complex vertical wavenumber. We show that an exponentially stratified atmosphere does not have any spatially bounded normal modes for real frequencies. The usual treatment involves a representation where the imaginary part of the vertical wavenumber yields a rho(sup -1/2) dependence of the velocity amplitude which diverges as the absolute value of z approaches infinity. This solution includes a cutoff frequency below which acoustic modes cannot propagate. The standard dispersion analysis is a local representation of the wave behavior in both space and time but which is assumed to represent the motion throughout - infinity is less than t is less than infinity and 0 is less than infinity. However, any solution which has a purely sinusoidal time dependence extends through this full domain and is divergent due to the rho(sup -1/2) dependence. We show that a proper description is in terms of a near field of a boundary piston which is driven arbitrarily as a function of space and time. The atmosphere which responds to this piston is a semi-infinite layer which has an initially constant sound speed but which has the usual gravitational stratification. In a restricted domain of space and time above this boundary, the wavelike behavior of the medium may be described by frequencies and vertical wavenumbers which are both complex. When both parameters are allowed to have imaginary components, a new range of solutions is found for which there is virtually no cutoff frequency. We show that vertical energy propagation can take place through the solar atmosphere as a result of oscillations below the nominal cutoff frequency. Previously, the largest amplitude oscillations which generally have low frequencies were dropped from the calculation of energy flux becuase their frequencies are below the cutoff frequency. This new family of near-field waves permits these modes to carry energy vertically outward and raises the possibility that the largest amplitude 5 minute oscillations play a substantial role in the transport of acoustic energy to the chromosphere.

Wang, Zhengzhi↗
NASA NTRS1991

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations was solved for in the steady state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time step, while the turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positively. Low Reynolds number modifications to the original two equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations.

Mavriplis, D. J.↗
NASA NTRS1991

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been solved for in the steady-state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time-step, while the turbulence equations are advanced in a point-implicit scheme with a time-step which guarantees stability and positivity. Low Reynolds number modifications to the original two-equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations.

Mavriplis, D. J.↗
NASA NTRS1972

Inviscid analysis of jet injection between two moving streams

An analytical method is developed for determining the flow interaction when a two-dimensional jet is injected between two moving streams. The jet is flowing out of channel and is turned as it enters between the external streams. The local velocity variation resulting from the flow interaction provides a static pressure variation along the jet bounding streamlines that is a priori unknown. Hense, the flow must be obtained by coupling the three flow regions (the jet and the free stream on either side) along the jet boundaries. Both external streams have the same total pressure, which is different from that in the jet. The solution is for the condition that the total pressure in the jet does not differ from the free-stream value by a large amount compared with the free-stream dynamic head. Results are given for the shape of the jet boundaries for various injection configurations.

Goldstein, M. E.↗
NASA NTRS1990

Emittance bounds for transient radiative cooling of a scattering rectangular region

Transient cooling was analyzed for a two-dimensional gray rectangular region that emits, absorbs, and isotropically scatters radiation. The region, initially at uniform temperature, is placed in surroundings at a much lower temperature. The cooling analysis yields two simultaneous equations for the transient temperature and scattering source function distributions. The region starts to cool with the medium at uniform temperature. For this condition, the local emittance around the perimeter was obtained by a numerical solution of the radiative equations, and the rate of heat loss from the entire medium was evaluated. As cooling continues, the emissive ability decreases because of the lower temperatures in the outer regions of the rectangle. Based on some previous work, a lower limit for the transient overall emittance was found by obtaining a similarity solution. For some ranges of optical thickness and scattering albedo, the lower limit is only a small amount below the initial emittance. For these conditions, the initial emittance can be used to compute the entire transient resulting in a considerable simplification.

Siegel, Robert↗
NASA NTRS1995

Gamma-ray transfer and energy deposition in supernovae

Solutions to the energy-independent (gray) radiative transfer equations are compared to results of Monte Carlo simulations of the Ni-56 and Co-56 decay gamma-ray energy deposition in supernovae. The comparison shows that an effective, purely absorptive, gray opacity, kappa(sub gamma) approximately (0. 06 +/- 0.01)Y(sub e) sq cm/g, where Y is the total number of electrons per baryon, accurately describes the interaction of gamma-rays with the cool supernova gas and the local gamma-ray energy deposition within the gas. The nature of the gamma-ray interaction process (dominated by Compton scattering in the relativistic regime) creates a weak dependence of kappa(sub gamma) on the optical thickness of the (spherically symmetric) supernova atmosphere: The maximum value of kappa(sub gamma) applies during optically thick conditions when individual gamma-rays undergo multiple scattering encounters and the lower bound is reached at the phase characterized by a total Thomson optical depth to the center of the atmosphere tau(sub e) approximately less than 1. Gamma-ray deposition for Type Ia supernova models to within 10% for the epoch from maximum light to t = 1200 days. Our results quantitatively confirm that the quick and efficient solution to the gray transfer problem provides an accurate representation of gamma-ray energy deposition for a broad range of supernova conditions.

Swartz, Douglas A.↗
NASA NTRS1994

Effect of magnetospheric convection on thermal plasma in the inner magnetosphere

The effects of E x B convection on the distribution of plasma parameters in the inner magnetosphere have been examined. Analytical solutions describing density distributions along convective trajectories in the equatorial plane have been found. These solutions suggest the following dependence of plasma concentration n on the magnetic field B along convective trajectories: n proportional to B(exp alpha), where the parameter alpha varies between alpha is less than or = between 4/3 and 2. The alpha = 2 case corresponds to disregarding transport parallel to B. The lower bound alpha = 4/3 describes the situation when parallel transport dominates over convective motions perpendicular to B. A solution has also been obtained describing the effect of convection on diffusive equilibrium profiles in the dayside plasmasphere. In addition, ion temperature variations due to adiabatic effects associated with plasma convection have been analyzed in detail. Convective drifts lead to ion temperature anisotropies with the value and sign of the anisotropy contingent on density and temperature variations, local time, and the location of the convective trajectory with respect to the location of the plasmapause. We have also found that convective motions lead to a substantial exchange of energy between the ionosphere and the plasmasphere through electron heat fluxes.

Khazanov, G. V.↗
NASA NTRS1989

The formation of sharp edges in planetary rings by nearby satellites

Equations are derived which govern the shapes of the perturbed streamlines near the 'sharp edge' boundaries between regions of high and low planetary ring optical depth; these are maintained by the shepherd satellites, which transfer angular momentum to and from ring particles. The results obtained by these equations' solution with a simple numerical model, whose parameters resemble those of the Encke division, are found to faithfully reproduce the sharp edges bounding the division; they imply that the ring thickness in the unperturbed regions far from the edges is of the order of 10 m, and that the angle-averaged surface density varies on a much shorter radial length scale than that over which the satellite torque is applied. This feature's relationship to the local reversal of angular momentum viscous transport, in the most strongly perturbed regions, is demonstrated.

Borderies, Nicole↗
NASA NTRS1988

Definitive IGRF models for 1945, 1950, 1955, and 1960

Definitive internal geomagnetic field models are obtained for the epochs 1945, 1950, 1955, and 1960, with the solution including local biases for the fixed observatories. In order to derive a priori models for the epochs, the International Geomagnetic Reference Field (12/83) model (based on Magsat data) was projected backward in time using a spline fit to annual secular-variation models. An uncertainty model was developed which included estimates of the effects of crustal and core fields not represented by the model. Differences between biases from year to year are found to be within the predicted bounds.

Langel, R. A.↗
NASA NTRS2014

Bifurcation Analysis Using Rigorous Branch and Bound Methods

For the study of nonlinear dynamic systems, it is important to locate the equilibria and bifurcations occurring within a specified computational domain. This paper proposes a new approach for solving these problems and compares it to the numerical continuation method. The new approach is based upon branch and bound and utilizes rigorous enclosure techniques to yield outer bounding sets of both the equilibrium and local bifurcation manifolds. These sets, which comprise the union of hyper-rectangles, can be made to be as tight as desired. Sufficient conditions for the existence of equilibrium and bifurcation points taking the form of algebraic inequality constraints in the state-parameter space are used to calculate their enclosures directly. The enclosures for the bifurcation sets can be computed independently of the equilibrium manifold, and are guaranteed to contain all solutions within the computational domain. A further advantage of this method is the ability to compute a near-maximally sized hyper-rectangle of high dimension centered at a fixed parameter-state point whose elements are guaranteed to exclude all bifurcation points. This hyper-rectangle, which requires a global description of the bifurcation manifold within the computational domain, cannot be obtained otherwise. A test case, based on the dynamics of a UAV subject to uncertain center of gravity location, is used to illustrate the efficacy of the method by comparing it with numerical continuation and to evaluate its computational complexity.

Smith, Andrew P.↗
NASA NTRS1975

Upwind second-order difference schemes and applications in unsteady aerodynamic flows

Explicit second-order upwind difference schemes in combination with spatially symmetric schemes can produce larger stability bounds and better numerical resolution than symmetric schemes alone. However, if conservation form is essential, a special operator is required for transition between schemes. An operational approach has been devised for deriving transition operators so that strict conservation and local consistency are maintained. Various aspects of hybrid schemes are studied numerically for model linear and nonlinear equations. To demonstrate the utility of combining two different algorithms, MacCormack's explicit, noncentered, second-order method is combined with a completely upwind version, and numerical solutions of the Euler equations are obtained for two-dimensional, transonic flows with embedded supersonic regions and shock waves. The general utility of the operational approach for combining schemes is emphasized by deriving a second-order conservative scheme for the steady transonic small-disturbance potential equation.

Warming, R. F.↗
NASA NTRS2004

Discovering Planetary Nebula Geometries: Explorations with a Hierarchy of Models

Astronomical objects known as planetary nebulae (PNe) consist of a shell of gas expelled by an aging medium-sized star as it makes its transition from a red giant to a white dwarf. In many cases this gas shell can be approximately described as a prolate ellipsoid. Knowledge of the physics of ionization processes in this gaseous shell enables us to construct a model in three dimensions (3D) called the Ionization-Bounded Prolate Ellipsoidal Shell model (IBPES model). Using this model we can generate synthetic nebular images, which can be used in conjunction with Hubble Space Telescope (HST) images of actual PNe to perform Bayesian model estimation. Since the IBPES model is characterized by thirteen parameters, model estimation requires the search of a 13-dimensional parameter space. The 'curse of dimensionality,' compounded by a computationally intense forward problem, makes forward searches extremely time-consuming and frequently causes them to become trapped in local solutions. We find that both the speed and of the search can be improved by judiciously reducing the dimensionality of the search space. Our basic approach employs a hierarchy of models of increasing complexity that converges to the IBPES model. Earlier studies establish that a hierarchical sequence converges more quickly, and to a better solution, than a search relying only on the most complex model. Here we report results for a hierarchy of five models. The first three models treat the nebula as a 2D image, while the last two models explore its characteristics as a 3D object and enable us to characterize the physics of the nebula. This five-model hierarchy is applied to HST images of ellipsoidal PNe to estimate their geometric properties and gas density profiles.

Huyser, Karen A.↗
NASA NTRS2017

GTOC9: Methods and Results from the Jet Propulsion Laboratory Team

The removal of 123 pieces of debris from the Sunsynchronous LEO environment is accomplished by a 10-spacecraft campaign wherein the spacecraft, flying in succession over an 8-yr period, rendezvous with a series of the debris objects, delivering a de-orbit package at each one before moving on to the next object by means of impulsive manoeuvres. This was the GTOC9 problem, as posed by the European Space Agency. The methods used by the Jet Propulsion Laboratory team are described, along with the winning solution found by the team. Methods include branch-and-bound searches that exploit the natural nodal drift to compute long chains of rendezvous with debris objects, beam searches for synthesising campaigns, ant colony optimisation, and a genetic algorithm. Databases of transfers between all bodies on a fine time grid are made, containing an easyto- compute yet accurate estimate of the transfer V . Lastly, a final non-linear programming optimisation is performed to ensure the trajectories meet all the constraints and are locally optimal in initial mass.

Sims, Jon↗
NASA NTRS1993

Fault geometries and extension in the Valles Marineris, Mars

The central troughs of the Valles Marineris system on Mars are generally regarded to be of tectonic origin because they are linear and bounded by faults. Their radial orientation to the Tharsis rise suggests that at least their incipient formation was influenced by the formation of the rise, but the exact role of Tharsis is not yet well understood. Collapse and erosion may have significantly modified the original tectonic troughs. To better understand the magnitude and mechanism of crustal extension across the Valles Marineris, the attitudes of trough-wall faults that parallel the long axes of the troughs were studied. Our results will help to determine if the extension is related to the Tharsis rise, to a local arch along which the troughs lie, or to other mechanisms. Fault-plane attitudes were obtained by performing a photogrammetric study of fault traces in the walls of Candor and Melas Chasmata. If the exposed wall faults are not major trough-bounding faults, it is assumed that these major faults are buried at the base of the walls and have similar geometries to the exposed faults. Three dimensional coordinates of several points along each fault trace were obtained and three-point solutions to their geometries were constructed. As expected, the results show that the faults dip toward the interior of the troughs, suggesting that they are indeed normal faults associated with trough formation.

Chadwick, D. J.↗
NASA NTRS1995

An approach to the development of numerical algorithms for first order linear hyperbolic systems in multiple space dimensions: The constant coefficient case

Two methods for developing high order single step explicit algorithms on symmetric stencils with data on only one time level are presented. Examples are given for the convection and linearized Euler equations with up to the eighth order accuracy in both space and time in one space dimension, and up to the sixth in two space dimensions. The method of characteristics is generalized to nondiagonalizable hyperbolic systems by using exact local polynominal solutions of the system, and the resulting exact propagator methods automatically incorporate the correct multidimensional wave propagation dynamics. Multivariate Taylor or Cauchy-Kowaleskaya expansions are also used to develop algorithms. Both of these methods can be applied to obtain algorithms of arbitrarily high order for hyperbolic systems in multiple space dimensions. Cross derivatives are included in the local approximations used to develop the algorithms in this paper in order to obtain high order accuracy, and improved isotropy and stability. Efficiency in meeting global error bounds is an important criterion for evaluating algorithms, and the higher order algorithms are shown to be up to several orders of magnitude more efficient even though they are more complex. Stable high order boundary conditions for the linearized Euler equations are developed in one space dimension, and demonstrated in two space dimensions.

Goodrich, John W.↗
NASA NTRS2013

Smart-Divert Powered Descent Guidance to Avoid the Backshell Landing Dispersion Ellipse

A smart-divert capability has been added into the Powered Descent Guidance (PDG) software originally developed for Mars pinpoint and precision landing. The smart-divert algorithm accounts for the landing dispersions of the entry backshell, which separates from the lander vehicle at the end of the parachute descent phase and prior to powered descent. The smart-divert PDG algorithm utilizes the onboard fuel and vehicle thrust vectoring to mitigate landing error in an intelligent way: ensuring that the lander touches down with minimum- fuel usage at the minimum distance from the desired landing location that also avoids impact by the descending backshell. The smart-divert PDG software implements a computationally efficient, convex formulation of the powered-descent guidance problem to provide pinpoint or precision-landing guidance solutions that are fuel-optimal and satisfy physical thrust bound and pointing constraints, as well as position and speed constraints. The initial smart-divert implementation enforced a lateral-divert corridor parallel to the ground velocity vector; this was based on guidance requirements for MSL (Mars Science Laboratory) landings. This initial method was overly conservative since the divert corridor was infinite in the down-range direction despite the backshell landing inside a calculable dispersion ellipse. Basing the divert constraint instead on a local tangent to the backshell dispersion ellipse in the direction of the desired landing site provides a far less conservative constraint. The resulting enhanced smart-divert PDG algorithm avoids impact with the descending backshell and has reduced conservatism.

Carson, John M.↗
NASA NTRS1993

On the accurate long-time solution of the wave equation in exterior domains: Asymptotic expansions and corrected boundary conditions

We consider the solution of scattering problems for the wave equation using approximate boundary conditions at artificial boundaries. These conditions are explicitly viewed as approximations to an exact boundary condition satisfied by the solution on the unbounded domain. We study the short and long term behavior of the error. It is provided that, in two space dimensions, no local in time, constant coefficient boundary operator can lead to accurate results uniformly in time for the class of problems we consider. A variable coefficient operator is developed which attains better accuracy (uniformly in time) than is possible with constant coefficient approximations. The theory is illustrated by numerical examples. We also analyze the proposed boundary conditions using energy methods, leading to asymptotically correct error bounds.

Hagstrom, Thomas↗
NASA NTRS2018

NICER Discovers the Ultracompact Orbit of the Accreting Millisecond Pulsar IGR J17062-6143

We present results of recent Neutron Star Interior Composition Explorer (NICER) observations of the accreting millisecond X-ray pulsar (AMXP) IGR J17062-6143 that show that it resides in a circular, ultracompact binary with a 38-minute orbital period. NICER observed the source for 26 kiloseconds over a 5.3-day span in 2017 August, and again for 14 and 11 kiloseconds in 2017 October and November, respectively. A power spectral analysis of the August exposure confirms the previous detection of pulsations at 163.656 Hertz in Rossi X-ray Timing Explorer (RXTE) data, and reveals phase modulation due to orbital motion of the neutron star. A coherent search for the orbital solution using the Z squared method finds a best-fitting circular orbit with a period of 2278.21 seconds (37.97 minutes), a projected semimajor axis of 0.00390 lt-s (Localization Test Statistic), and a barycentric pulsar frequency of 163.6561105 Hertz. This is currently the shortest known orbital period for an AMXP. The mass function is 9.12 times 10 (sup minus 8) solar mass, presently the smallest known for a stellar binary. The minimum donor mass ranges from approximately 0.005 to 0.007 times the solar mass for a neutron star mass from 1.2 to 2 times the solar mass. Assuming mass transfer is driven by gravitational radiation, we find donor mass and binary inclination bounds of 0.0175-0.0155 times the solar mass and 19 degrees less than i less than 27.5 degrees, where the lower and upper bounds correspond to 1.4 and 2 times the solar mass neutron stars, respectively. Folding the data accounting for the orbital modulation reveals a sinusoidal profile with fractional amplitude 2.04 plus or minus 0.11 percent (0.3-3.2 kiloelectronvolts).

Strohmayer, T. E.↗