SEARCH · Engineering Papers

Results for “Solutions”

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 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Benchmark solutions for the galactic ion transport equations with spatial and energy coupling

In order to anticipate future space shielding requirements, NASA has initiated an effort to formulate computational methods to simulate radiation effects in space. As part of the program, numerical transport algorithms have been developed for the deterministic Boltzman equation describing galactic cosmic ray (GCR) interactions with matter. It thus becomes necessary to assess the accuracy of proposed deterministic algorithms. For this reason, analytical benchmark solutions to mathematically tractable galactic cosmic ray equations have recently been obtained. Even though these problems involve simplifying assumptions of the associated physics, they still contain the essential features of the basic transport processes. The solutions obtained are features of the basic transport processes. The solutions obtained are compared to results from numerical algorithms in order to ensure proper coding and to provide a measure of the accuracy of the numerical methods used in the algorithm. For the first time, mathematical methods have been applied to the galactic ion transport (GIT) equations in the straight ahead approximation with constant nuclear properties. The approach utilizes a Laplace transforms inversion yielding a closed form benchmark solution which is also computationally efficient.

Ganapol, Barry D.↗

NASA NTRS1989

Numerical solution of unsteady incompressible viscous flows in generalized moving coordinate systems

A solution method of the time-accurate, incompressible Navier-Stokes equations in generalized curvilinear moving coordinate systems is presented in this paper. Accuracy is achieved by a conservative finite-volume discretization which satisfies the geometric conservation laws in generalized moving coordinate systems. The solution method is second-order accurate in space and first-order accurate in time. A fractional step solution method is used to efficiently solve the discrete equations. The unknowns, namely the pressure and the volume-fluxes, are chosen to facilitate the formulation of a consistent Poisson equation and to obtain a robust Poisson solver with favorable convergence properties. The method is validated by comparisons to other numerical and experimental solutions. The comparisons show good agreement.

Rosenfeld, Moshe↗

NASA NTRS1989

Error analysis of finite element solutions for postbuckled cylinders

A general method of error analysis and correction is investigated for the discrete finite-element results for cylindrical shell structures. The method for error analysis is an adaptation of the method of successive approximation. When applied to the equilibrium equations of shell theory, successive approximations derive an approximate continuous solution from the discrete finite-element results. The advantage of this continuous solution is that it contains continuous partial derivatives of an order higher than the basis functions of the finite-element solution. Preliminary numerical results are presented in this paper for the error analysis of finite-element results for a postbuckled stiffened cylindrical panel modeled by a general purpose shell code. Numerical results from the method have previously been reported for postbuckled stiffened plates. A procedure for correcting the continuous approximate solution by Newton's method is outlined.

Sistla, Rajaram↗

NASA NTRS1989

Adaptive numerical solutions of the Euler equations in 3D using finite elements

The development of an adaptive mesh solution for a flow involving shock interaction on a swept cylinder and an initial solution for a flow past a complex fighter configuration is reported. The finite element solution algorithm, the mesh generation, and the adaptivity of the solution are described. Sample results for the flow past an F-18 configuration at Mach 0.9 and alpha of 3 deg and for shock interaction on a swept cylinder at Mach 8.04 are summarized.

Peraire, J.↗

NASA NTRS1989

A block-corrected subdomain solution procedure for recirculating flow calculations

This paper describes a robust and efficient subdomain solution procedure for two-dimensional recirculating flows. The solution domain is divided into a number of overlapping subdomains, and a direct fully coupled solution is obtained for each subdomain using a sparse matrix form of LU decomposition. An effective parabolic block correction procedure, which calculates global corrections to the tentative solution by a marching technique similar to that used for boundary layer flows, is used to accelerate the convergence of the basic procedure. The use of effective block correction is found to be essential for the success of the subdomain approach on strongly recirculating flows. In a number of laminar two-dimensional flows, the new block-corrected method performed extremely well, rivaling the best direct methods in execution time, while requiring substantially less computer storage. The new method proved to be from two to ten times faster than conventional iterative methods, while requiring only a moderate increase in storage.

Braaten, M. E.↗

NASA NTRS1989

Closed-form solution for loop transfer recovery via reduced-order observers

A well-known property of the reduced-order observer is exploited to obtain the controller solution of the loop transfer recovery problem. In that problem, the controller is sought that generates some desired loop shape at the plant's input or output channels. Past approaches to this problem have typically yielded controllers generating loop shapes that only converge pointwise to the desired loop shape. In the proposed approach, however, the solution (at the input) is obtained directly when the plant's first Markov parameter is full rank. In the more general case when the plant's first Markov parameter is not full rank, the solution is obtained in an analogous manner by appending a special set of input and output signals to the original set. A dual form of the reduced-order observer is shown to yield the LTR solution at the output channel.

Bacon, Barton J.↗

NASA NTRS1989

Maximum principle solutions for time-optimal half-loop maneuvers of a high alpha fighter aircraft

An investigation was conducted of maximum principle solutions for an initial 0.6 Mach number and 15,000-ft altitude. The authors generate these solutions for a family of prescribed final times tf, starting with tf = 0.5 s. Using a nonlinear wind-tunnel model they construct maximum principle solutions. Above tf = 1.2 s some small nonlinear variations in the aerodynamic pitching moment coefficient presented difficulty with respect to numerical convergence. This was circumvented by fitting analytical models to the aerodynamic coefficients of the wind-tunnel model at Mach 0.4. Maximum principle solutions of the analytical model are shown to compare well with those obtained for tf of less than 1.2 s. Using the analytical model the authors extended the prescribed final time to a value of 13.65 s at which time the aircraft completes the half-loop maneuver. This is 0.53 s longer than that obtained using the singular perturbation feedback control law.

Stalford, Harold↗

NASA NTRS1990

High accuracy solutions of incompressible Navier-Stokes equations

In recent years, high accuracy finite difference approximations were developed for partial differential equations of elliptic type, with particular emphasis on the convection-diffusion equation. These approximations are of compact type, have a local truncation error of fourth order, and allow the use of standard iterative schemes to solve the resulting systems of algebraic equations. These high accuracy approximations are extended to the solution of Navier-Stokes equations. Solutions are obtained for the model problem of driven cavity and are compared with solutions obtained using other approximations and those obtained by other authors. It is discovered that the high order approximations do indeed produce high accuracy solutions and have a potential for use in solving important problems of viscous fluid flows.

Gupta, Murli M.↗

NASA NTRS1990

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

NASA NTRS1990

Application of finite-element-based solution technologies for viscoplastic structural analyses

Finite-element solution technology developed for use in conjunction with advanced viscoplastic models is described. The development of such solution technology is necessary for performing stress/life analyses of engineering structural problems where the complex geometries and loadings make the conventional analytical solutions difficult. The versatility of the solution technology is demonstrated by applying it to viscoplastic models possessing different mathematical structures and encompassing isotropic and anisotropic material. The computational results qualitatively replicate deformation behavior observed in experiments on prototypical structural components.

Arya, V. K.↗

NASA NTRS1989

Solution of steady-state one-dimensional conservation laws by mathematical programming

Solution techniques for a class of steady-state scalar conservation laws are developed analytically. Discretization by finite-volume formulas is employed to obtain an overdetermined system of algebraic equations, which are then perturbed nonsingularly (with perturbation coefficient = epsilon) and solved using the l(1) mathematical-programming algorithm of Seneta and Steiger (1984); this approach limits the matrix bandwidth to two, so that an explicit solution can be found efficiently. It is shown that, for small values of epsilon, the l(1) solutions exhibit sharp correctly located shocks and are nonoscillatory O(epsilon) approximations of the physically relevant solutions.

Lavery, J. E.↗

NASA NTRS1989

Comparison of Burnett, super-Burnett and Monte Carlo solutions for hypersonic shock structure

The continuum Navier-Stokes, Burnett, and super-Burnett equations are solved for one-dimensional shock structures in various monoatomic gases. Solutions for a hard sphere gas, argon, and a Maxwellian gas from Mach 1.3 to Mach 50 are obtained. A numerical method utilizing the complete time-dependent continuum equations and obtaining the steady-state shock structure by allowing the system to relax from arbitrary initial conditions is employed. Shock density, velocity, temperature, and entropy profiles are also obtained using the direct simulation Monte Carlo method, and these results are used as bases for comparison for continuum solution profiles. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Solutions to the super-Burnett equations with coefficients as presently derived, however, are considered to be inferior to those of the Burnett equations.

Fiscko, Kurt A.↗

NASA NTRS1990

A block-based algorithm for the solution of compressible flows in rotor-stator combinations

A block-based solution algorithm is developed for the solution of compressible flows in rotor-stator combinations. The method allows concurrent solution of multiple solution blocks in parallel machines. It also allows a time averaged interaction at the stator-rotor interfaces. Numerical results are presented to illustrate the performance of the algorithm. The effect of the interaction between the stator and rotor is evaluated.

Akay, H. U.↗

NASA NTRS1990

Nitric acid uptake by sulfuric acid solutions under stratospheric conditions - Determination of Henry's Law solubility

The uptake of nitric acid by sulfuric acid solutions representative of stratospheric particulate at low temperatures was measured to determine the solubility of nitric acid in sulfuric acid solutions as a function of H2SO4 concentration and solution temperature. Solubilities are reported for sulfuric acid solutions ranging from 58 to 87 wt pct H2SO4 over a temperature range from 188 to 240 K, showing that, in general, the solubility of nitric acid increases with decreasing sulfuric acid concentration and with decreasing temperature. The measured solubilities indicate that nitric acid in the global stratosphere will be found predominantly in the gas phase.

Reihs, Christa M.↗

NASA NTRS1991

Closed-form solutions for linear regulator design of mechanical systems including optimal weighting matrix selection

Vibration in modern structural and mechanical systems can be reduced in amplitude by increasing stiffness, redistributing stiffness and mass, and/or adding damping if design techniques are available to do so. Linear Quadratic Regulator (LQR) theory in modern multivariable control design, attacks the general dissipative elastic system design problem in a global formulation. The optimal design, however, allows electronic connections and phase relations which are not physically practical or possible in passive structural-mechanical devices. The restriction of LQR solutions (to the Algebraic Riccati Equation) to design spaces which can be implemented as passive structural members and/or dampers is addressed. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical system. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist.

Hanks, Brantley R.↗

NASA NTRS1991

Development of iterative techniques for the solution of unsteady compressible viscous flows

Efficient iterative solution methods are being developed for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. Thus, the extra work required by iterative schemes can also be designed to perform efficiently on current and future generation scalable, missively parallel machines. An obvious candidate for iteratively solving the system of coupled nonlinear algebraic equations arising in CFD applications is the Newton method. Newton's method was implemented in existing finite difference and finite volume methods. Depending on the complexity of the problem, the number of Newton iterations needed per step to solve the discretized system of equations can, however, vary dramatically from a few to several hundred. Another popular approach based on the classical conjugate gradient method, known as the GMRES (Generalized Minimum Residual) algorithm is investigated. The GMRES algorithm was used in the past by a number of researchers for solving steady viscous and inviscid flow problems with considerable success. Here, the suitability of this algorithm is investigated for solving the system of nonlinear equations that arise in unsteady Navier-Stokes solvers at each time step. Unlike the Newton method which attempts to drive the error in the solution at each and every node down to zero, the GMRES algorithm only seeks to minimize the L2 norm of the error. In the GMRES algorithm the changes in the flow properties from one time step to the next are assumed to be the sum of a set of orthogonal vectors. By choosing the number of vectors to a reasonably small value N (between 5 and 20) the work required for advancing the solution from one time step to the next may be kept to (N+1) times that of a noniterative scheme. Many of the operations required by the GMRES algorithm such as matrix-vector multiplies, matrix additions and subtractions can all be vectorized and parallelized efficiently.

Sankar, Lakshmi N.↗

NASA NTRS1991

Exact and approximate solutions to the oblique shock equations for real-time applications

The derivation of exact solutions for determining the characteristics of an oblique shock wave in a supersonic flow is investigated. Specifically, an explicit expression for the oblique shock angle in terms of the free stream Mach number, the centerbody deflection angle, and the ratio of the specific heats, is derived. A simpler approximate solution is obtained and compared to the exact solution. The primary objectives of obtaining these solutions is to provide a fast algorithm that can run in a real time environment.

Hartley, T. T.↗

NASA NTRS1991

Analytical solution for boundary heat fluxes from a radiating rectangular medium

Reference is made to the work of Shah (1979) which demonstrated the possibility of partially integrating the radiative equations analytically to obtain an 'exact' solution. Shah's solution was given as a double integration of the modified Bessel function of order zero. Here, it is shown that the 'exact' solution for a rectangular region radiating to cold black walls can be conveniently derived, and expressed in simple form, by using an integral function, Sn, analogous to the exponential integral function appearing in plane-layer solutions.

Siegel, R.↗