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 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Correcting GOES-R Magnetometer Data for Stray Fields

Time-varying spacecraft magnetic fields, i.e. stray fields, are a problem for magnetometer systems. While constant fields can be removed by calibration, stray fields are difficult to distinguish from ambient field variations. Putting two magnetometers on a long boom and solving for both the ambient and stray fields can help, but this gradiometer solution is more sensitive to noise than a single magnetometer. As shown here for the R-series Geostationary Operational Environmental Satellites (GOES-R), unless the stray fields are larger than the noise, simply averaging the two magnetometer readings gives a more accurate solution. If averaging is used, it may be worthwhile to estimate and remove stray fields explicitly. Models and estimation algorithms to do so are provided for solar array, arcjet and reaction wheel fields.

Magnetometer↗

NASA NTRS1976

Solutions to divergence form equations using the method of partial implicitization

Partial implicitization is well suited for use on vector-processing computers since the method is an explicit unconditionally stable numerical technique. So far equations cast in nondivergence form have been used, though equations cast in the divergence form would have inherently greater accuracy. Since the method of partial implicitization is applicable only to relaxing problems to their steady state and is not applicable to true transient problems, a simple mathematical manipulation can be carried out which will produce the desired partial implicitization solution to divergence form equations.

Graves, R. A., Jr.↗

NASA NTRS1987

Supercomputing of supersonic flows using upwind relaxation and MacCormack schemes

The performance of two numerical solution schemes, (1) an implicit upwind relaxation with a finite-volume discretization (Thomas and Walters, 1985) and (2) an explicit-implicit MacCormack (1981) scheme with a finite-difference discretization, is compared in two-dimensional simulations of supersonic flow past a flat plate with leading edge, a rearward-facing step, a 10-deg compression corner, a NACA 0012 airfoil at high angle of attack, and a cavity. The algebraic turbulence model, the solution methods, and the boundary conditions and SIMD coding are explained, and the results are presented in tables and graphs and characterized with reference to published experimental data. Scheme (1) is found to converge more rapidly and to give more accurate results than (2) in a wide range of problem types.

Baysal, O.↗

NASA NTRS1993

A detailed study of mean-flow solutions for stability analysis of transitional flows

A finite-volume upwind-difference parabolized Navier-Stokes code is utilized to obtain laminar mean-flow solutions at Mach 3.5 on a half-angle cone of 5 deg at an angle-of-attack of 2 deg. A detailed study is conducted on this configuration; the main focus is the velocity profiles in the leeward and windward symmetry planes at various axial locations. Comparisons of the solution profiles are made with both a central-difference code that incorporates scalar and matrix dissipation models and another state-of-the-art upwind-difference finitevolume code. The results obtained emphasize the importance of using matrix dissipation models for schemes that require explicit artificial dissipation. These results also illustrate the accuracy and efficiency of the planeby-plane marching procedure for computing mean-flow solutions for predicting the onset of transition with linear instability.

Ramakrishnan, R.↗

NASA NTRS1985

Accelerated solution of the steady Euler equations

The present paper is concerned with two methods for the accelerated solution of the steady Euler equations. One method makes use of a second-order embedding to facilitate the derivation of the relaxation solution of the steady equations of motion, while the other method employs a multile-gridding concept to accelerate the convergence of a simple, explicit, time-marching scheme applied to the unsteady equations. It is pointed out that the surrogate equation technique provides a means for formulating problems involving the full steady Euler equations in such a way as to allow the use of relaxation solution procedures. It is, therefore, possible to solve either irrotational or rotational flow problems spanning the entire spectrum of subsonic, transonic, and supersonic conditions. The solutions can be obtained without an employement of either derived dependent variables, semidirect methods, or an unsteady formulation.

Johnson, G. M.↗

NASA NTRS1985

Numerical solution of the two-dimensional Euler equations by second-order upwind difference schemes

Two time-level, five-point explicit and implicit upwind difference schemes based on the characteristic flux difference splitting concept have been developed for the two-dimensional Euler equations. The method is conservative, second-order accurate in time and space, and general coordinate systems are used to treat complex geometries. Nonlinear flux limiters are employed to yield oscillation free sharp shock profiles. Upstream interpolation is used to yield a class of higher-order upwind schemes which closely mimic the locally one-dimensional method of characteristics (with fixed time intervals) through operator splitting. Numerical results have been obtained for a plane shock reflection and for flow over a circular arc in a channel. Characteristics of upwind TVD schemes, as applied to two-dimensional flows with embedded shocks are discussed.

Yang, J. Y.↗

NASA NTRS1990

A package for 3-D unstructured grid generation, finite-element flow solution and flow field visualization

A set of computer programs for 3-D unstructured grid generation, fluid flow calculations, and flow field visualization was developed. The grid generation program, called VGRID3D, generates grids over complex configurations using the advancing front method. In this method, the point and element generation is accomplished simultaneously, VPLOT3D is an interactive, menudriven pre- and post-processor graphics program for interpolation and display of unstructured grid data. The flow solver, VFLOW3D, is an Euler equation solver based on an explicit, two-step, Taylor-Galerkin algorithm which uses the Flux Corrected Transport (FCT) concept for a wriggle-free solution. Using these programs, increasingly complex 3-D configurations of interest to aerospace community were gridded including a complete Space Transportation System comprised of the space-shuttle orbitor, the solid-rocket boosters, and the external tank. Flow solutions were obtained on various configurations in subsonic, transonic, and supersonic flow regimes.

Parikh, Paresh↗

NASA NTRS2018

Analytic Expressions for Derivatives from Series Solutions to the Three Body Problem

This paper presents a notation system to facilitate to solution of diﬀerential equations via Taylor series expansions and applies it to solve the circular restricted three body problem. Unlike previous Taylor series methods in the astrodynamics literature, computer algebra solvers are not used. Instead the notation system allows one to solve a system of diﬀerential equations analytically “by hand” without resorting to computer algebra software. This method produces recurrence relations explicitly in terms of a sequence of derivatives of the state with respect to time for the coeﬃcients of Taylor Series solutions that can be evaluated numerically or manipulated further to investigate properties of the solution. For example, additional derivatives with respect to other parameters may also be found, including those that describe the dependence of the solution on initial conditions.

Strange, Nathan↗

NASA NTRS1992

Short-arc orbit determination using coherent X-band ranging data

The use of X-band frequencies in ground-spacecraft and spacecraft-ground telecommunication links for current and future robotic interplanetary missions makes it possible to perform ranging measurements of greater accuracy than previously obtained. It is shown that ranging data of sufficient accuracy, when acquired from multiple stations, can sense the geocentric angular position of a distant spacecraft. The application of high-accuracy S/X-band and X-band ranging to orbit determination with relatively short data arcs is investigated in planetary approach and encounter scenarios. Actual trajectory solutions for the Ulysses spacecraft constructed from S/X-band ranging and Doppler data are presented; error covariance calculations are used to predict the performance of X-band ranging and Doppler data. The Ulysses trajectory solutions indicate that the aim point for the spacecraft's February 1992 Jupiter encounter was predicted to a geocentric accuracy of 0.20 to 0.23/microrad. Explicit modeling of range bias parameters for each station pass is shown to largely remove systematic ground system calibration errors and transmission media effects from the Ulysses range measurements, which would otherwise corrupt the angle finding capabilities of the data. The Ulysses solutions were found to be reasonably consistent with the theoretical results, which suggest that angular accuracies of 0.08 to 0.1/microrad are achievable with X-band ranging.

Thurman, S. W.↗

NASA NTRS2011

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗

NASA NTRS1973

The electric field in the vicinity of a photo-emitting plate in a plasma

Knowledge of the electric field in the vicinity of a spacecraft provides necessary corrections to measurements of ambient fields and charged particle fluxes and a required parameter to understand the motion of dust particles and grains around the spacecraft. The complex problem of the entire electric field, needed for grain transport calculations, is treated for the one-dimensional case. By approximations to the analysis of a photo-emitting plate immersed in a dilute plasma, explicit expressions for the potential and the electric field were derived. The two classes of implicit solutions for the potential reported by them are shown to yield essentially identical results for the conditions in our solar system. The values were obtained following assumptions of a later paper by Guernsey and Fu, in which the emitted photoelectron distribution is taken as Maxwellian. These results are compared with a simple model due to Grand and Tunaley, wherein the effect of the dilute plasma is neglected entirely.

Barengoltz, J.↗

NASA NTRS1980

Relating viscosity to polymer concentration

Equation developed by VanDijk and first applied to viscosity of Newtonian suspension of rigid particles by Eilers is rearranged to yield intrinsic viscosity as explicit function of polymer concentration in polymer solvent system. Experiments have shown relationship valid for polymer solutions having relative viscosities ranging from 1 to 100.

Fedors, R. F.↗

NASA NTRS1984

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradients through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallelepiped computational domain, and an explicit MacCormack (1972) explicit time-split predictor-corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗

NASA NTRS1987

Calculation of viscous effects on transonic flow for oscillating airfoils and comparisons with experiment

A method is described for calculating unsteady transonic flow with viscous interaction by coupling a steady integral boundary-layer code with an unsteady, transonic, inviscid small-disturbance computer code in a quasi-steady fashion. Explicit coupling of the equations together with viscous -inviscid iterations at each time step yield converged solutions with computer times about double those required to obtain inviscid solutions. The accuracy and range of applicability of the method are investigated by applying it to four AGARD standard airfoils. The first-harmonic components of both the unsteady pressure distributions and the lift and moment coefficients have been calculated. Comparisons with inviscid calcualtions and experimental data are presented. The results demonstrate that accurate solutions for transonic flows with viscous effects can be obtained for flows involving moderate-strength shock waves.

Howlett, James T.↗

NASA NTRS1990

A finite element solver for 3-D compressible viscous flows

Computation of the flow field inside a space shuttle main engine (SSME) requires the application of state of the art computational fluid dynamic (CFD) technology. Several computer codes are under development to solve 3-D flow through the hot gas manifold. Some algorithms were designed to solve the unsteady compressible Navier-Stokes equations, either by implicit or explicit factorization methods, using several hundred or thousands of time steps to reach a steady state solution. A new iterative algorithm is being developed for the solution of the implicit finite element equations without assembling global matrices. It is an efficient iteration scheme based on a modified nonlinear Gauss-Seidel iteration with symmetric sweeps. The algorithm is analyzed for a model equation and is shown to be unconditionally stable. Results from a series of test problems are presented. The finite element code was tested for couette flow, which is flow under a pressure gradient between two parallel plates in relative motion. Another problem that was solved is viscous laminar flow over a flat plate. The general 3-D finite element code was used to compute the flow in an axisymmetric turnaround duct at low Mach numbers.

Reddy, K. C.↗

NASA NTRS1992

On the stability and instantaneous velocity of grasped frictionless objects

A quantitative test for form closure valid for any number of contact points is formulated as a linear program, the optimal objective value of which provides a measure of how far a grasp is from losing form closure. Another contribution of the study is the formulation of a linear program whose solution yields the same information as the classical approach. The benefit of the formulation is that explicit testing of all possible combinations of contact interactions can be avoided by the algorithm used to solve the linear program.

Trinkle, Jeffrey C.↗

NASA NTRS1993

Two-integral distribution functions for axisymmetric galaxies

The new method presented for finding distribution functions, which depend only on the classical integrals of energy and angular momentum for stellar systems with known axisymmetric densities, is the analog for the axisymmetric case of Eddington's classical solution for the isotropic distribution function, depending only on energy, of a known spherical density. It is required that density be expressed as a function of the potential and of a radial coordinate. Our solution is also an integral which is derived directly from the density, and hence can be used with complicated densities. A numerical quadrature is generally required to evaluate this solution, but contour integrals can be computed accurately by numerical quadrature; this is preferable to an explicit evaluation if the latter is an infinite series, such as is obtained using Fricke's method. We give several examples, including some for which our distribution functions are new. Our method can be extended simply to the related problems of finding anisotropic distribution functions for spherical or disk systems.

Hunter, C.↗

NASA NTRS2000

Parallel CE/SE Computations via Domain Decomposition

This paper describes the parallelization strategy and achieved parallel efficiency of an explicit time-marching algorithm for solving conservation laws. The Space-Time Conservation Element and Solution Element (CE/SE) algorithm for solving the 2D and 3D Euler equations is parallelized with the aid of domain decomposition. The parallel efficiency of the resultant algorithm on a Silicon Graphics Origin 2000 parallel computer is checked.

Himansu, Ananda↗