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 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1987

A numerical solution of Duffing's equations including the prediction of jump phenomena

Numerical methodology for the solution of Duffing's differential equation is presented. Algorithms for the prediction of multiple equilibrium solutions and jump phenomena are developed. In addition, a filtering algorithm for producing steady state solutions is presented. The problem of a rigidly clamped circular plate subjected to cosinusoidal pressure loading is solved using the developed algorithms (the plate is assumed to be in the geometrically nonlinear range). The results accurately predict regions of solution multiplicity and jump phenomena.

Moyer, E. T., Jr.↗

NASA NTRS1987

On similarity solutions of a boundary layer problem with an upstream moving wall

The problem of a boundary layer on a flat plate which has a constant velocity opposite in direction to that of the uniform mainstream is examined. It was previously shown that the solution of this boundary value problem is crucially dependent on the parameter which is the ratio of the velocity of the plate to the velocity of the free stream. In particular, it was proved that a solution exists only if this parameter does not exceed a certain critical value, and numerical evidence was adduced to show that this solution is nonunique. Using Crocco formulation the present work proves this nonuniqueness. Also considered are the analyticity of solutions and the derivation of upper bounds on the critical value of wall velocity parameter.

Hussaini, M. Y.↗

NASA NTRS1988

The use of multigrid techniques in the solution of the Elrod algorithm for a dynamically loaded journal bearing

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed, utilizing a multigrid iterative technique. The code is compared with a presently existing direct solution in terms of computational time and accuracy. The model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobssen-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via liquid striations. The mixed nature of the equations (elliptic in the full film zone and nonelliptic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗

NASA NTRS1988

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗

NASA NTRS1988

Simplified solution of diffraction from a Lyot system

This paper presents a derivation of a simplified analytical solution of diffraction from a Lyot (1939) system designed for observation of the solar corona outside of eclipses. Applying the theorem of Papoulis (1986) to simplify the calculations, a simplified solution is derived which is found to be in reasonable agreement with the exact solution. The simplified solution suffices for the preliminary evaluation of the amount of diffraction reduction needed to meet certain system requirements and also serves as a guideline for further apodization.

Wang, Yaujen↗

NASA NTRS1987

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗

NASA NTRS1988

On hierarchical solutions to the BBGKY hierarchy

It is thought that the gravitational clustering of galaxies in the universe may approach a scale-invariant, hierarchical form in the small separation, large-clustering regime. Past attempts to solve the Born-Bogoliubov-Green-Kirkwood-Yvon (BBGKY) hierarchy in this regime have assumed a certain separable hierarchical form for the higher order correlation functions of galaxies in phase space. It is shown here that such separable solutions to the BBGKY equations must satisfy the condition that the clustered component of the solution has cluster-cluster correlations equal to galaxy-galaxy correlations to all orders. The solutions also admit the presence of an arbitrary unclustered component, which plays no dyamical role in the large-clustering regime. These results are a particular property of the specific separable model assumed for the correlation functions in phase space, not an intrinsic property of spatially hierarchical solutions to the BBGKY hierarchy. The observed distribution of galaxies does not satisfy the required conditions. The disagreement between theory and observation may be traced, at least in part, to initial conditions which, if Gaussian, already have cluster correlations greater than galaxy correlations.

Hamilton, A. J. S.↗

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.↗