Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Relaxation method”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Real-time optimal guidance for orbital maneuvering.

A new formulation for soft-constraint trajectory optimization is presented as a real-time optimal feedback guidance method for multiburn orbital maneuvers. Control is always chosen to minimize burn time plus a quadratic penalty for end condition errors, weighted so that early in the mission (when controllability is greatest) terminal errors are held negligible. Eventually, as controllability diminishes, the method partially relaxes but effectively still compensates perturbations in whatever subspace remains controllable. Although the soft-constraint concept is well-known in optimal control, the present formulation is novel in addressing the loss of controllability inherent in multiple burn orbital maneuvers. Moreover the necessary conditions usually obtained from a Bolza formulation are modified in this case so that the fully hard constraint formulation is a numerically well behaved subcase. As a result convergence properties have been greatly improved.

Cohen, A. O.↗
NASA NTRS1974

Input description for Jameson's three-dimensional transonic airfoil analysis program

The input parameters are presented for a computer program which performs calculations for inviscid isentropic transonic flow over three dimensional airfoils with straight leading edges. The free stream Mach number is restricted only by the isentropic assumption. Weak shock waves are automatically located where they occur in the flow. The finite difference form of the full equation for the velocity potential is solved by the method of relaxation, after the flow exterior to the airfoil is mapped to the upper half plane.

Newman, P. A.↗
NASA NTRS1974

Transonic flow about lifting wing-body combinations

Calculations for transonic flow about wing-body combinations at angles of attack of the order of the body thickness-to-length ratio are presented. These calculations show that lift can alter the flow drastically at near-sonic speeds in this angle-of-attack range so that the area rule must be modified. This effect of lift is explained from simple physical considerations. The computational procedure is based on slender-wing theory and a two-variable method of relaxation solution. A small perturbation analysis is used to show that this simple procedure is adequate for these transonic wing-body flow fields.

Barnwell, R. W.↗
NASA NTRS1976

Two inviscid computational simulations of separated flow about airfoils

Two inviscid computational simulations of separated flow about airfoils are described. The basic computational method is the line relaxation finite-difference method. Viscous separation is approximated with inviscid free-streamline separation. The point of separation is specified, and the pressure in the separation region is calculated. In the first simulation, the empiricism of constant pressure in the separation region is employed. This empiricism is easier to implement with the present method than with singularity methods. In the second simulation, acoustic theory is used to determine the pressure in the separation region. The results of both simulations are compared with experiment.

Barnwell, R. W.↗
NASA NTRS1976

A brief description of the Jameson-Caughey NYU transonic swept-wing computer program: FLO 22

A computer program for analyzing inviscid, isentropic, transonic flow past 3-D swept configurations is presented. Some basic aspects of the program are: (1) the free-stream Mach number is restricted only by the isentropic assumption; (2) weak shock waves are automatically located wherever they occur in the flow; (3) the finite-difference form of the full equation for the velocity potential is solved by the method of relaxation, after the flow exterior to the airfoil is mapped to the upper half plane; (4) the mapping procedure allows exact satisfaction of the boundary conditions and use of supersonic free stream velocities; (5) the finite difference operator is locally rotated in supersonic flow regions so as to properly account for the domain of dependence; and (6) the relaxation algorithm was stabilized using criteria from a time-like analogy.

Jameson, A.↗
NASA NTRS1981

The possible modifications of the HISSE model for pure LANDSAT agricultural data

A method for relaxing the assumption of class conditional independence of LANDSAT spectral measurements within the same patch (field) is discussed. Theoretical arguments are given which show that any significant refinement of the model beyond this proposal will not allow the reduction, essential to HISSE, of the pure data to patch summary statistics. A slight alteration of the new model is shown to be a reasonable approximation to the model which describes pure data elements from the same patch as jointly Gaussian with a covariance function which exhibits exponential decay with respect to spatial separation.

Peters, C.↗
NASA NTRS1996

Differences in Mechanical Properties of the Human and Monkey Tibia

A method which uses an instrument that detects the response of a long bone to a vibratory stimulus to quantify mechanical properties non-invasively was revised and validated for use in the tibia. Stored data from healthy men was reanalyzed and compared with values from non-human primates. The analysis uses the relationship K(sub b) = 48 EI/L(sup 3) where K(sub b) is the lateral stiffness of a beam with force applied midspan, E is the elastic modulus, I the geometric moment of inertia and L, the limb length. Values for stiffness (EI, Nm(sup2)), the Euler buckling load (P(sub cr) = EI (pi/L)(sup 2)), and bone sufficiency (S) which represents the axial load the bone can support, adjusted to BW (S=P(sub cr)/BW) were obtained. The interest precision of the method in relaxed men, 5.8%, and in sedated male monkeys, 4.3%, was based on repeated measures in the same subjects at 1 month intervals. The R tibias of 40 men, aged 38.6 +/- 7.3 yrs with BW 78.9 +/- 7.9 kg, showed average (+/- SD) L to be 35 +/- 2 cm, EI 222 +/- 71 Nm(sup 2), P(sub cr) 18.1 +/- 4.9 kN, and S 23.4 +/- 5.7 N. The R tibias of 24 Rhesus monkeys ranging in age from 2-12 years, BW 4.9 +/- 3 kg, showed L to be 14.7 +/- 1.9 cm, EI 6.0 +/- 4.8 Nm(sup 2), P(sub cr) 2.51 +/- 1.2 kN and S 57.3 N. These measurements indicate that the tibia of a terrestrial non-human primate, M. mulatta, has higher load carrying capacity for the level of body weights in the species than the human bone.

Arnaud, Sara B.↗
NASA NTRS2012

Effect of Surface Nonequilibrium Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using finite-rate gas/surface interaction model provides time-accurate solutions for multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal Response and Ablation Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas momentum conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of gas/surface interaction chemistry between air and carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was a Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kang↗
NASA NTRS2012

Effect of Non-Equilibrium Surface Thermochemistry in Simulation of Carbon Based Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver using non-equilibrium gas/surface interaction model provides time-accurate solutions for the multidimensional ablation of carbon based charring ablators. The material thermal response code used in this study is the Two-dimensional Implicit Thermal-response and AblatioN Program (TITAN), which predicts charring material thermal response and shape change on hypersonic space vehicles. Its governing equations include total energy balance, pyrolysis gas mass conservation, and a three-component decomposition model. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation (DPLR) method. Loose coupling between the material response and flow codes is performed by solving the surface mass balance in DPLR and the surface energy balance in TITAN. Thus, the material surface recession is predicted by finite-rate gas/surface interaction boundary conditions implemented in DPLR, and the surface temperature and pyrolysis gas injection rate are computed in TITAN. Two sets of nonequilibrium gas/surface interaction chemistry between air and the carbon surface developed by Park and Zhluktov, respectively, are studied. Coupled fluid-material response analyses of stagnation tests conducted in NASA Ames Research Center arc-jet facilities are considered. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator (PICA). Computational predictions of in-depth material thermal response and surface recession are compared with the experimental measurements for stagnation cold wall heat flux ranging from 107 to 1100 Watts per square centimeter.

Chen, Yih-Kanq↗
NASA NTRS1973

Relaxation factors for supercritical flows.

Relaxation procedures for solution of steady supercritical transonic flows are investigated. Von Neumann (Fourier-mode) stability analysis is used to find bounds of relaxation factors. The bounds depend on local Mach number and local mesh aspect ratio. Long wave instability of Murman-Cole implicit method is indicated. Two new relaxation procedures are introduced. Both employ central differencing exclusively. Group velocities of Fourier modes are used to study signal propagation. It was found necessary to avoid or to damp out signals propagating upstream in supersonic zones in order to obtain physically meaningful transonic solutions. It appears that requirements of high rate of convergence, of stability and accuracy, are in conflict, and that a combination of relaxation methods must be used in order not to compromise the requirements.

Kentzer, C. P.↗
NASA NTRS1973

Relaxation time measurements by an electronic method.

Description of a simple electronic system that permits the direct measurement of time constants of decaying signals. The system was used in connection with relaxation experiments on hydrogen and rubidium masers and was found to operate well. The use of a computing counter in the systems gives the possibility of making averages on several experiments and obtaining the standard deviation of the results from the mean. The program for the computing counter is given.

Brousseau, R.↗
NASA NTRS1986

Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 1: One-step method

An algorithm for solving a large class of two- and three-dimensional nonseparable elliptic partial differential equations (PDE's) is developed and tested. It uses a modified D'Yakanov-Gunn iterative procedure in which the relaxation factor is grid-point dependent. It is easy to implement and applicable to a variety of boundary conditions. It is also computationally efficient, as indicated by the results of numerical comparisons with other established methods. Furthermore, the current algorithm has the advantage of possessing two important properties which the traditional iterative methods lack; that is: (1) the convergence rate is relatively insensitive to grid-cell size and aspect ratio, and (2) the convergence rate can be easily estimated by using the coefficient of the PDE being solved.

Chang, S. C.↗
NASA NTRS1982

Modal control of relaxed static stability aircraft

A method is developed that assigns a selected portion of a closed loop system eigenstructure in accordance with certain desirable criteria. The method is applied here to a relaxed static stability aircraft, the goal being to synthesize a control law that provides the unstable aircraft with handling qualities equal to or better than those of a comparable statically stable aircraft. It is shown that by using the target system eigenstructure, good flight characteristics are achieved by the unstable aircraft. It is also shown that improved characteristics can be obtained by assigning an orthogonal eigenvector structure.

Rooney, R. H.↗
NASA NTRS1977

Implicit approximate-factorization schemes for the efficient solution of steady transonic flow problems

Implicit approximate-factorization algorithms (AF) are developed for the solution of steady-state transonic flow problems. The performance of the AF solution method is evaluated relative to that of the standard solution method for transonic flow problems, successive line over-relaxation (SLOR). Both methods are applied to the solution of the nonlinear, two-dimensional transonic small-disturbance equation. Results indicate that the AF method requires substantially less computer time than SLOR to solve the nonlinear finite-difference matrix equation for a transonic flow field. This increase in computational efficiency is achieved with no appreciable increase in computer storage or coding complexity.

Ballhaus, W. F.↗
NASA NTRS2005

Numerical Evaluation of P-Multigrid Method for the Solution of Discontinuous Galerkin Discretizations of Diffusive Equations

This paper describes numerical experiments with P-multigrid to corroborate analysis, validate the present implementation, and to examine issues that arise in the implementations of the various combinations of relaxation schemes, discretizations and P-multigrid methods. The two approaches to implement P-multigrid presented here are equivalent for most high-order discretization methods such as spectral element, SUPG, and discontinuous Galerkin applied to advection; however it is discovered that the approach that mimics the common geometric multigrid implementation is less robust, and frequently unstable when applied to discontinuous Galerkin discretizations of di usion. Gauss-Seidel relaxation converges 40% faster than block Jacobi, as predicted by analysis; however, the implementation of Gauss-Seidel is considerably more expensive that one would expect because gradients in most neighboring elements must be updated. A compromise quasi Gauss-Seidel relaxation method that evaluates the gradient in each element twice per iteration converges at rates similar to those predicted for true Gauss-Seidel.

Atkins, H. L.↗
NASA NTRS1986

Solution of elliptic partial differential equations by fast Poisson solvers using a local relaxation factor. 2: Two-step method

A two-step semidirect procedure is developed to accelerate the one-step procedure described in NASA TP-2529. For a set of constant coefficient model problems, the acceleration factor increases from 1 to 2 as the one-step procedure convergence rate decreases from + infinity to 0. It is also shown numerically that the two-step procedure can substantially accelerate the convergence of the numerical solution of many partial differential equations (PDE's) with variable coefficients.

Chang, S. C.↗