SEARCH · Search NASA

Results for “explicit solution”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Multidimensional difference schemes with fourth-order accuracy

An explicit finite-difference algorithm is presented for the solution of quasilinear divergence free multidimensional hyperbolic systems. The method consists of four steps per time level. The resulting scheme is fourth-order accurate in both space and time, though the intermediate steps are only first-order accurate. The family of schemes introduced is dissipative, and hence, suitable for both smooth flows and flows containing shocks. This method is compared, in several numerical examples, with both second-order schemes and others that are fourth order in space, but second order in time.

Turkel, E.↗

NASA NTRS1981

Remarks on nonstandard least-squares problems

It is shown that after treating a least-squares problem in the standard fashion, the imposition of stochastic restrictions on certain parameters can be accomplished essentially by inverting a matrix of an order no larger than the total number of restricted parameters, and that generally there is an explicit term for correcting the unrestricted to the restricted solution. In addition, any known type of nonstandard least-squares adjustment problem can be reduced to an equivalent standard linear problem, that is, one in which there is exactly one observation per condition equation and in which the observations are not correlated.

Eichhorn, H.↗

NASA NTRS1984

Multiple steady states for characteristic initial value problems

The time dependent, isentropic, quasi-one-dimensional equations of gas dynamics and other model equations are considered under the constraint of characteristic boundary conditions. Analysis of the time evolution shows how different initial data may lead to different steady states and how seemingly anamolous behavior of the solution may be resolved. Numerical experimentation using time consistent explicit algorithms verifies the conclusions of the analysis. The use of implicit schemes with very large time steps leads to erroneous results.

Salas, M. D.↗

NASA NTRS1989

Second-order accurate nonoscillatory schemes for scalar conservation laws

Explicit finite difference schemes for the computation of weak solutions of nonlinear scalar conservation laws is presented and analyzed. These schemes are uniformly second-order accurate and nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time.

Huynh, Hung T.↗

NASA NTRS1990

Characteristic-based algorithms for flows in thermo-chemical nonequilibrium

A generalized finite-rate chemistry algorithm with Steger-Warming, Van Leer, and Roe characteristic-based flux splittings is presented in three-dimensional generalized coordinates for the Navier-Stokes equations. Attention is placed on convergence to steady-state solutions with fully coupled chemistry. Time integration schemes including explicit m-stage Runge-Kutta, implicit approximate-factorization, relaxation and LU decomposition are investigated and compared in terms of residual reduction per unit of CPU time. Practical issues such as code vectorization and memory usage on modern supercomputers are discussed.

Walters, Robert W.↗

NASA NTRS1989

Second-order accurate nonoscillatory schemes for scalar conservation laws

Huynh, Hung T.↗

NASA NTRS1991

Numerical identification of boundary conditions on nonlinearly radiating inverse heat conduction problems

An explicit and unconditionally stable finite difference method for the solution of the transient inverse heat conduction problem in a semi-infinite or finite slab mediums subject to nonlinear radiation boundary conditions is presented. After measuring two interior temperature histories, the mollification method is used to determine the surface transient heat source if the energy radiation law is known. Alternatively, if the active surface is heated by a source at a rate proportional to a given function, the nonlinear surface radiation law is then recovered as a function of the interface temperature when the problem is feasible. Two typical examples corresponding to Newton cooling law and Stefan-Boltzmann radiation law respectively are illustrated. In all cases, the method predicts the surface conditions with an accuracy suitable for many practical purposes.

Murio, Diego A.↗

NASA NTRS1992

Klein-Gordon equation and reflection of Alfven waves in nonuniform media

A new analytical approach is presented for assessing the reflection of linear Alfven waves in smoothly nonuniform media. The general one-dimensional case in Cartesian coordinates is treated. It is shown that the wave equations, upon transformation into the form of the Klein-Gordon equation, display a local critical frequency for reflection. At any location in the medium, reflection becomes strong as the wave frequency descends past this characteristic frequency set by the local nonuniformity of the medium. This critical frequecy is given by the transformation as an explicit function of the Alfven velocity and its first and second derivatives, and hence as an explicit spatial function. The transformation thus directly yields, without solution of the wave equations, the location in the medium at which an Alfven wave of any given frequency becomes strongly reflected and has its propagation practically cut off.

Musielak, Z. E.↗

NASA NTRS1990

Diffuse reflection from a stochastically bounded, semi-infinite medium

In order to determine the diffuse reflection from a medium bounded by a rough surface, the problem of radiative transfer in a boundary layer characterized by a statistical distribution of heights is considered. For the case that the surface is defined by a multivariate normal probability density, the propagation probability for rays traversing the boundary layer is derived and, from that probability, a corresponding radiative transfer equation. A solution of the Eddington (two stream) type is found explicitly, and examples are given. The results should be applicable to reflection from the regoliths of solar system bodies, as well as from a rough ocean surface.

Lumme, K.↗

NASA NTRS1991

The placement of equipment in the Space Station Freedom using constraint based reasoning

This paper describes the Rack Equipment Placement and Optimization System. The primary objective of this system is to assist engineers with the placement of equipment into the racks of the modules of Space Station Freedom. It accomplishes this by showing a user where equipment placement is possible and by generating potential layouts. The system uses an explicit representation of integration constraints to search for potential solutions for individual rack equipment items. A simulated annealing process is being evaluated for total solution generation as well. Versions of this system are in use now and are assisting with the development of the Space Station Freedom at the Marshall Space Flight Center in Huntsville, Alabama.

Tanner, Steve↗

NASA NTRS1990

Improved boundary layer heat transfer calculations near a stagnation point

A thermal design of a solar receiver has been developed for the solutions of problems involving phase-change thermal energy storage and natural convection loss. Two dimensional axisymmetrical solidification and melting of materials contained between two concentric cylinders of finite length has been studied for thermal energy storage analysis. For calculation of free convection loss inside receiver cavity, two dimensional axisymmetrical, laminar, transient free convection including radiation effects has been studied using integral/finite difference method. Finite difference equations are derived for the above analysis subject to constant or variable material properties, initial conditions, and boundary conditions. The validity of the analyses has been substantiated by comparing results of the present general method with available analytic solutions or numerical results reported in the literature. Both explicit and implicit schemes are tested in phase change analysis with different number of nodes ranging from 4 to 18. The above numerical methods have been applied to the existing solar receiver analyzing computer code as additional subroutines. The results were computed for one of the proposed Brayton cycle solar receiver models running under the actual environmental conditions. Effect of thermal energy storage on the thermal behavior of the receiver has been estimated. Due to the thermal energy storage, about 65% reduction on working gas outlet temperature fluctuation has been obtained; however, maximum temperature of thermal energy storage containment has been increased about 18%. Also, effect of natural convection inside a receiver cavity on the receiver heat transfer has been analyzed. The finding indicated that thermal stratification occurs during the sun time resulting in higher receiver temperatures at the outlet section of the gas tube, and lower temperatures at the inlet section of the gas tube when compared with the results with no natural convection. Due to heat supply from the air during the shade time, minimum temperature has been increased, while maximum temperature has been reduced due to convection loss to air. Consequently, cyclic temperature fluctuation has been reduced 29% for working gas and 16% for thermal energy storage containment. On the other hand, despite the presence of the natural convection the time-averaged temperatures for receiver components were found to be similar for two cases with/without natural convection (maximum difference was 1.8%).

Ahn, Kyung Hwan↗

NASA NTRS1996

Utilization of STRAT Data

One of the goals of the STRAT mission was to improve our understanding of transport in the stratosphere. We have used our coupled 2-D model to analyse transport between the tropics and midlatitudes using aircraft measurements of (CO2) and (N2O) made during the SPADE and ASHOE/MAESA campaigns (Boering et al. 1994) for model validation. We have further used the model to estimate errors in the experimental determination of the age of stratospheric air that arise from nonlinear increases of real tracers such as (CO2) and (SF6). The interpretation of measurements made during STRAT is still in progress. Most of the work performed under this contract has been described in the second year progress report, submitted in April of 1996. In the following, the research is summarized and some additional work that has become relevant in the meantime is described. We have developed an interactive model of the dynamics, chemistry and radiation of the stratosphere with partial support from this contract. The dynamics module integrates the primitive equations on a sphere. Lower boundary zonal winds and eddy geopotential heights are specified from observations. Chemistry and tracer transport are done in two dimensions. The chemistry model is a family model. Long lived species are transported and short lived species are determined in intervals of 5 to 10 days by assuming photochemical steady state. Full diurnal integrations with the fast species are performed periodically and diurnal coefficients, derived from the explicit integration are used to determine the steady state solution. Heating rates in the stratosphere are calculated from model ozone and the model temperatures and circulation fields are used in the chemistry transport module.

Source record↗

NASA NTRS2004

OTIS 3.2 Software Released

Trajectory, mission, and vehicle engineers concern themselves with finding the best way for an object to get from one place to another. These engineers rely upon special software to assist them in this. For a number of years, many engineers have used the OTIS program for this assistance. With OTIS, an engineer can fully optimize trajectories for airplanes, launch vehicles like the space shuttle, interplanetary spacecraft, and orbital transfer vehicles. OTIS provides four modes of operation, with each mode providing successively stronger optimization capability. The most powerful mode uses a mathematical method called implicit integration to solve what engineers and mathematicians call the optimal control problem. OTIS 3.2, which was developed at the NASA Glenn Research Center, is the latest release of this industry workhorse and features new capabilities for parameter optimization and mission design. OTIS stands for Optimal Control by Implicit Simulation, and it is implicit integration that makes OTIS so powerful at solving trajectory optimization problems. Why is this so important? The optimization process not only determines how to get from point A to point B, but it can also determine how to do this with the least amount of propellant, with the lightest starting weight, or in the fastest time possible while avoiding certain obstacles along the way. There are numerous conditions that engineers can use to define optimal, or best. OTIS provides a framework for defining the starting and ending points of the trajectory (point A and point B), the constraints on the trajectory (requirements like "avoid these regions where obstacles occur"), and what is being optimized (e.g., minimize propellant). The implicit integration method can find solutions to very complicated problems when there is not a lot of information available about what the optimal trajectory might be. The method was first developed for solving two-point boundary value problems and was adapted for use in OTIS. Implicit integration usually allows OTIS to find solutions to problems much faster than programs that use explicit integration and parametric methods. Consequently, OTIS is best suited to solving very complicated and highly constrained problems.

Riehl, John P.↗

NASA NTRS1994

New multigrid approach for three-dimensional unstructured, adaptive grids

A new multigrid method with adaptive unstructured grids is presented. The three-dimensional Euler equations are solved on tetrahedral grids that are adaptively refined or coarsened locally. The multigrid method is employed to propagate the fine grid corrections more rapidly by redistributing the changes-in-time of the solution from the fine grid to the coarser grids to accelerate convergence. A new approach is employed that uses the parent cells of the fine grid cells in an adapted mesh to generate successively coaser levels of multigrid. This obviates the need for the generation of a sequence of independent, nonoverlapping grids as well as the relatively complicated operations that need to be performed to interpolate the solution and the residuals between the independent grids. The solver is an explicit, vertex-based, finite volume scheme that employs edge-based data structures and operations. Spatial discretization is of central-differencing type combined with a special upwind-like smoothing operators. Application cases include adaptive solutions obtained with multigrid acceleration for supersonic and subsonic flow over a bump in a channel, as well as transonic flow around the ONERA M6 wing. Two levels of multigrid resulted in reduction in the number of iterations by a factor of 5.

Parthasarathy, Vijayan↗

NASA NTRS1986

Computation of three-dimensional, rotational flow through turbomachinery blade rows for improved aerodynamic design studies

The performance of a three dimensional computer code developed for predicting the flowfield in stationary and rotating turbomachinery blade rows is described in this study. The four stage Runge-Kutta numerical integration scheme is used for solving the governing flow equations and yields solution to the full, three dimensional, unsteady Euler equations in cylindrical coordinates. This method is fully explicit and uses the finite volume, time marching procedure. In order to demonstrate the accuracy and efficiency of the code, steady solutions were obtained for several cascade geometries under widely varying flow conditions. Computed flowfield results are presented for a fully subsonic turbine stator and a low aspect ratio, transonic compressor rotor blade under maximum flow and peak efficiency design conditions. Comparisons with Laser Anemometer measurements and other numerical predictions are also provided to illustrate that the present method predicts important flow features with good accuracy and can be used for cost effective aerodynamic design studies.

Subramanian, S. V.↗

NASA NTRS2003

Route Monopolie and Optimal Nonlinear Pricing

To cope with air traffic growth and congested airports, two solutions are apparent on the supply side: 1) use larger aircraft in the hub and spoke system; or 2) develop new routes through secondary airports. An enlarged route system through secondary airports may increase the proportion of route monopolies in the air transport market.The monopoly optimal non linear pricing policy is well known in the case of one dimension (one instrument, one characteristic) but not in the case of several dimensions. This paper explores the robustness of the one dimensional screening model with respect to increasing the number of instruments and the number of characteristics. The objective of this paper is then to link and fill the gap in both literatures. One of the merits of the screening model has been to show that a great varieD" of economic questions (non linear pricing, product line choice, auction design, income taxation, regulation...) could be handled within the same framework.VCe study a case of non linear pricing (2 instruments (2 routes on which the airline pro_ddes customers with services), 2 characteristics (demand of services on these routes) and two values per characteristic (low and high demand of services on these routes)) and we show that none of the conclusions of the one dimensional analysis remain valid. In particular, upward incentive compatibility constraint may be binding at the optimum. As a consequence, they may be distortion at the top of the distribution. In addition to this, we show that the optimal solution often requires a kind of form of bundling, we explain explicitly distortions and show that it is sometimes optimal for the monopolist to only produce one good (instead of two) or to exclude some buyers from the market. Actually, this means that the monopolist cannot fully apply his monopoly power and is better off selling both goods independently.We then define all the possible solutions in the case of a quadratic cost function for a uniform distribution of agent types and explain the implications for airlines in terms of service differentiation.

Tournut, Jacques↗

NASA NTRS1986

Interface procedures for overlapping grids

The solution of partial differential equations by means of either implicit or conservative explicit finite-difference methods, on the one hand, or multicomponent overlapping grid schemes on the other, is presently subjected to interpolation at grid boundaries in an attempt to eliminate the time lag. The accuracy of the transient solution of a hyperbolic or parabolic partial differential equation is dependent on the procedures used to transfer information between grids in a composite grid system. The present techniques reduce numerical solution errors.

Mastin, C. Wayne↗

NASA NTRS1973

Martian thermal boundary layers: Subhourly variations induced by radiative-conductive heat transfer within the dust-laden atmosphere-ground system

Radiative-conductive heat transfer has been investigated for the ground-atmosphere system of the planet Mars. The basic goal was the quantitative determination of time dependent vertical distributions of temperature and static stability for Southern-Hemispheric summer season and middle and polar latitudes, for both dust-free and dust-laden atmospheric conditions. The numerical algorithm which models at high spatial and temporal resolution the thermal energy transports in the dual ground-atmosphere system, is based on solution of the applicable heating rate equation, including radiative and molecular-conductive heat transport terms. The two subsystems are coupled by an internal thermal boundary condition applied at the ground-atmosphere interface level. Initial data and input parameters are based on Mariner 4, 6, 7, and 9 measurements and the JPL Mars Scientific Model. Numerical experiments were run for dust-free and dust-laden conditions in the midlatitudes, as well as ice-free and ice-covered polar regions. Representative results and their interpretation are presented. Finally, the theoretical framework of the generalized problem with nonconservative Mie scattering and explicit thermal-convective heat transfer is formulated, and applicable solution algorithms are outlined.

Pallmann, A. J.↗