Search NASA⌕ Search

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.

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

At least 109 records · Page 6

The thermal statics of coronal loops

The thermal statics of constant pressure coronal loops is discussed, with particular emphasis on non-equilibrium and scaling relations. An analytical solution showing explicitly the occurrence of non-equilibrium in radiation dominated loops is presented. In addition, the general scaling law for hot loops is given. However, in view of the uncertainties in the coronal heating function and the observational determined loop parameters, it is suggested that scaling laws are currently of limited value.

Roberts, B.↗

Periodic solutions of the Riccati equation

It is shown that the Riccati equations of control theory are generically Morse-Smale. The existence of periodic solutions is explicitly demonstrated. The principal contribution here is thought to be the presentation of a mathematical structure in which to study the Riccati equation. Mathematical structures already in the literature are systematically exploited. It is also shown that most problems of the Riccati equation can be reduced to problems in linear algebra. A framework is developed for studying the Riccati equation, and the proof that the Morse-Smale conditions are satisfied is outlined.

Hermann, R.↗

Boundary layer effects on the flowfield about flow-through nacelles

A numerical study has been made of boundary layer effects on the flowfield about flow-through nacelles. A viscous-inviscid interacting computational model for investigating the problem was constructed by coupling a three-dimensional explicit Euler solution procedure with a compressible, 'lag-entrainment' integral boundary layer solution technique. Solutions obtained with the interacting model for long-duct turbofan engine nacelles at a free stream Mach number of 0.8 and an angle of attack of 0 deg are compared with experimental data. The numerical experiments have revealed insight into the boundary layer effects on the interactions between the internal and external flows.

Compton, W. B.↗

Finite element Euler computations in three dimensions

A two-step explicit FEM solution algorithm for the three-dimensional compressible Euler and Navier-Stokes equations based on unstructured triangular and tetrahedral grids is described and demonstrated. The method represents an extension and refinement of the algorithms presented by Loehner et al. (1984 and 1985), Peraire et al. (1987), and Morgan et al. (1987). The formulation and numerical implementation are outlined; the mesh generation, data structures, and adaptive remeshing are explained; and results for a two-dimensional airfoil, a three-dimensional engine air intake, a B747 in landing configuration, and a generic fighter aircraft are presented in extensive graphics and discussed in detail.

Peraire, Jaime↗

Calculation of the angular radiance distribution for a coupled atmosphere and canopy

The radiative transfer equations for a coupled atmosphere and canopy are solved numerically by an improved Gauss-Seidel iteration algorithm. The radiation field is decomposed into three components: unscattered sunlight, single scattering, and multiple scattering radiance for which the corresponding equations and boundary conditions are set up and their analytical or iterational solutions are explicitly derived. The classic Gauss-Seidel algorithm has been widely applied in atmospheric research. This is its first application for calculating the multiple scattering radiance of a coupled atmosphere and canopy. This algorithm enables us to obtain the internal radiation field as well as radiances at boundaries. Any form of bidirectional reflectance distribution function (BRDF) as a boundary condition can be easily incorporated into the iteration procedure. The hotspot effect of the canopy is accommodated by means of the modification of the extinction coefficients of upward single scattering radiation and unscattered sunlight using the formulation of Nilson and Kuusk. To reduce the computation for the case of large optical thickness, an improved iteration formula is derived to speed convergence. The upwelling radiances have been evaluated for different atmospheric conditions, leaf area index (LAI), leaf angle distribution (LAD), leaf size and so on. The formulation presented in this paper is also well suited to analyze the relative magnitude of multiple scattering radiance and single scattering radiance in both the visible and near infrared regions.

Liang, Shunlin↗

High performance techniques for space mission scheduling

In this paper, we summarize current research at Carnegie Mellon University aimed at development of high performance techniques and tools for space mission scheduling. Similar to prior research in opportunistic scheduling, our approach assumes the use of dynamic analysis of problem constraints as a basis for heuristic focusing of problem solving search. This methodology, however, is grounded in representational assumptions more akin to those adopted in recent temporal planning research, and in a problem solving framework which similarly emphasizes constraint posting in an explicitly maintained solution constraint network. These more general representational assumptions are necessitated by the predominance of state-dependent constraints in space mission planning domains, and the consequent need to integrate resource allocation and plan synthesis processes. First, we review the space mission problems we have considered to date and indicate the results obtained in these application domains. Next, we summarize recent work in constraint posting scheduling procedures, which offer the promise of better future solutions to this class of problems.

Smith, Stephen F.↗

A Study of Wake Development and Structure in Constant Pressure Gradients

Motivated by the application to high-lift aerodynamics for commercial transport aircraft, a systematic investigation into the response of symmetric/asymmetric planar turbulent wake development to constant adverse, zero, and favorable pressure gradients has been conducted. The experiments are performed at a Reynolds number of 2.4 million based on the chord of the wake generator. A unique feature of this wake study is that the pressure gradients imposed on the wake flow field are held constant. The experimental measurements involve both conventional LDV and hot wire flow field surveys of mean and turbulent quantities including the turbulent kinetic energy budget. In addition, similarity analysis and numerical simulation have also been conducted for this wake study. A focus of the research has been to isolate the effects of both pressure gradient and initial wake asymmetry on the wake development. Experimental results reveal that the pressure gradient has a tremendous influence on the wake development, despite the relatively modest pressure gradients imposed. For a given pressure gradient, the development of an initially asymmetric wake is different from the initially symmetric wake. An explicit similarity solution for the shape parameters of the symmetric wake is obtained and agrees with the experimental results. The turbulent kinetic energy budget measurements of the symmetric wake demonstrate that except for the convection term, the imposed pressure gradient does not change the fundamental flow physics of turbulent kinetic energy transport. Based on the turbulent kinetic energy budget measurements, an approach to correct the bias error associated with the notoriously difficult dissipation estimate is proposed and validated through the comparison of the experimental estimate with a direct numerical simulation result.

Thomas, Flint O.↗

The Dispersion Relation for the 1/sinh(exp 2) Potential in the Classical Limit

The dispersion relation for the inverse hyperbolic potential is calculated in the classical limit. This is shown for both the low amplitude phonon branch and the high amplitude soliton branch. It is shown these results qualitatively follow that previously found for the inverse squared potential where explicit analytic solutions are known.

Campbell, Joel↗

Finite element solution for energy conservation using a highly stable explicit integration algorithm

Theoretical derivation of a finite element solution algorithm for the transient energy conservation equation in multidimensional, stationary multi-media continua with irregular solution domain closure is considered. The complete finite element matrix forms for arbitrarily irregular discretizations are established, using natural coordinate function representations. The algorithm is embodied into a user-oriented computer program (COMOC) which obtains transient temperature distributions at the node points of the finite element discretization using a highly stable explicit integration procedure with automatic error control features. The finite element algorithm is shown to posses convergence with discretization for a transient sample problem. The condensed form for the specific heat element matrix is shown to be preferable to the consistent form. Computed results for diverse problems illustrate the versatility of COMOC, and easily prepared output subroutines are shown to allow quick engineering assessment of solution behavior.

Baker, A. J.↗

Explicit numerical formulas of improved stability and accuracy for the solution of parabolic equations

A class of explicit numerical formulas which involve next nearest neighbor as well as nearest neighbor points are explored in this paper. These formulas are formal approximations to the linear parabolic partial-differential equation of first order in time and second order in distance. It was found that some of these formulas can employ time steps as much as four times that for the conventional explicit technique without becoming unstable. Others showed improved accuracy for a given time step and spatial grid spacing. One formula achieved a steady-state solution of specified accuracy for an example problem in less than 4 percent of the total computational time required by the conventional explicit technique.

Olstad, W. B.↗

Time-asymptotic solutions of the Navier-Stokes equation for free shear flows using an alternating-direction implicit method

An uncoupled time asymptotic alternating direction implicit method for solving the Navier-Stokes equations was tested on two laminar parallel mixing flows. A constant total temperature was assumed in order to eliminate the need to solve the full energy equation; consequently, static temperature was evaluated by using algebraic relationship. For the mixing of two supersonic streams at a Reynolds number of 1,000, convergent solutions were obtained for a time step 5 times the maximum allowable size for an explicit method. The solution diverged for a time step 10 times the explicit limit. Improved convergence was obtained when upwind differencing was used for convective terms. Larger time steps were not possible with either upwind differencing or the diagonally dominant scheme. Artificial viscosity was added to the continuity equation in order to eliminate divergence for the mixing of a subsonic stream with a supersonic stream at a Reynolds number of 1,000.

Rudy, D. H.↗

A well posed boundary value problem in transonic gas dynamics

A boundary value problem for the Tricomi equation was studied in connection with transonic gas dynamics. The transformed equation delta u plus 1/3Y u sub Y equals 0 in canonical coordinates was considered in the complex domain of two independent complex variables. A boundary value problem was then set by prescribing the real part of the solution on the boundary of the real unit circle. The Dirichlet problem in the upper unit semicircle with vanishing values of the solution at Y = 0 was solved explicitly in terms of the hypergeometric function for the more general Euler-Poisson-Darboux equation. An explicit representation of the solution was also given for a mixed Dirichlet and Neumann problem for the same equation and domain.

Sanz, J. M.↗

Modeling the effects of latitudinal gradients in stellar winds, with application to the solar wind

A steady, axisymmetric, quasi-radial, global model previously developed for stellar winds with embedded magnetic fields has been extended to include latitudinal gradient effects on the azimuthal velocity and magnetic field. The linear results at large radii are presented for large-amplitude latitudinal variations in the radial magnetic field, mass loss rate, and radial velocity of the wind. The magnetohydrodynamic (MHD) equations predict meridional flows that develop naturally from internal magnetic stresses. The flows open flux tubes in the star's equatorial plane, redistributing mass and magnetic flux as a function of stellar latitude. The plasma spins up to conserve angular momentum in fields and plasma. The results are generally applicable to stellar winds (including radiatively driven winds), provided that the internal structure is not dominated by rotation. The asymptotic solutions do not explicitly depend on the form of the energy equation, although the assumed O(1) state which drives these solutions depends on the deposition of energy and momentum throughout the wind.

Nerney, S.↗

Introducing a Degree of Freedom in the Family of Generalized Logarithmic Spirals

The versatility of the family of generalized logarithmic spirals is improved by introducing a degree of freedom in the solution. The low-thrust acceleration profile now includes a control term that affects both the magnitude and the direction of the thrust. Exact and fully analytic solutions to the trajectory, the velocity, the time of flight, etc. are made available. Two integrals of motion are preserved. The first one is a generalization of the equation of the energy and depends on the values of the control parameter. The second one relates to the equation of the angular momentum. The problem of finding spiral transfers between two arbitrary state vectors reduces to solving one algebraic equation with one unknown. The degree of freedom allows fixing the time of flight of the transfer. If the time of flight is fixed, then there are two equations with two unknowns. No other iterative procedures are required. Coast arcs can be introduced in the solution naturally. An explicit expression for the maximum acceleration reached along the transfer is provided. Thanks to the symmetry properties of the solution a simple algorithm for generating periodic orbits is presented. An arbitrary number of intermediate nodes can be introduced to improve the flexibility of the solution when facing optimization problems. An example of a low-thrust gravity-assist Earth-Mars-Ceres trajectory shows that the solution is comparable to that obtained with other preliminary design techniques.

Roa, Javier↗

Nonlinear truncation error analysis of finite difference schemes for the Euler equations

It is pointed out that, in general, dissipative finite difference integration schemes have been found to be quite robust when applied to the Euler equations of gas dynamics. The present investigation considers a modified equation analysis of both implicit and explicit finite difference techniques as applied to the Euler equations. The analysis is used to identify those error terms which contribute most to the observed solution errors. A technique for analytically removing the dominant error terms is demonstrated, resulting in a greatly improved solution for the explicit Lax-Wendroff schemes. It is shown that the nonlinear truncation errors are quite large and distributed quite differently for each of the three conservation equations as applied to a one-dimensional shock tube problem.

Klopfer, G. H.↗

A vectorized, finite-volume, adaptive-grid algorithm for Navier-Stokes calculations

An adaptive grid, finite-volume method has been used to solve the Navier-Stokes equations for complete (forebody and afterbody) flowfields around blunt bodies. The code, which is applicable for axisymmetric or two-dimensional flows, allows the mesh to adjust during the computation to provide a closer spacing of mesh points in regions of high gradients, thus minimizing the number of required computational points. The solution technique is explicit, utilizing a maximum time-step advancement at each grid point to accelerate convergence to the steady state. The code has been fully vectorized for efficient solution on the CYBER 203 computer. A very flexible rezoning routine is used to concentrate mesh points anywhere in the field, either by a user-defined weighting function or by allowing high gradient regions to adjust the grid. The grid adjustment routine is implicit in nature and represents a very small portion of the total computational cost. Currently, the code runs in approximately 0.000016 seconds per grid point per iteration.

Gnoffo, P. A.↗