SEARCH · Engineering Papers

Results for “analytical 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 379 records · Page 21

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Solving time-dependent two-dimensional eddy current problems

Results of transient eddy current calculations are reported. For simplicity, a two-dimensional transverse magnetic field which is incident on an infinitely long conductor is considered. The conductor is assumed to be a good but not perfect conductor. The resulting problem is an interface initial boundary value problem with the boundary of the conductor being the interface. A finite difference method is used to march the solution explicitly in time. The method is shown. Treatment of appropriate radiation conditions is given special consideration. Results are validated with approximate analytic solutions. Two stringent test cases of high and low frequency incident waves are considered to validate the results.

Lee, Min Eig↗

NASA NTRS1990

Preconditioned upwind methods to solve 3-D incompressible Navier-Stokes equations for viscous flows

A computational method for calculating low-speed viscous flowfields is developed. The method uses the implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. An extensive study of optimizing the preconditioned system is carried out for two viscous flow problems. Computed results are compared with analytical solutions and experimental data.

Hsu, C.-H.↗

NASA NTRS1991

A Markov model reduction technique for fault tolerant processor reliability analysis

A fault tolerant processor (FTP) plays a key role in many high performance, safety-critical control system applications. Realistic modeling of an FTP is crucial to gaining a high degree of confidence in the reliability and safety analysis of such a system. While fidelity is clearly a major consideration, a practical model must also be kept to a moderate size to allow its incorporation into the overall system model. This paper presents a systematic reduction technique that starts from a complex, detailed model of a triple, redundant FTP and produces a low order approximation of very high accuracy. The existence of two distinct time scales represents the key to the success of the technique. No eigenvalue solution or coordinate transformation are needed. The reduced model captures all the important features of the detailed model, is amenable to an analytical solution and provides insight into the reconfiguration behavior of an FTP.

Schor, Andrei L.↗

NASA NTRS1992

The double universal joint wrist on a manipulator: Solution of inverse position kinematics and singularity analysis

This paper presents three methods to solve the inverse position kinematics position problem of the double universal joint attached to a manipulator: (1) an analytical solution for two specific cases; (2) an approximate closed form solution based on ignoring the wrist offset; and (3) an iterative method which repeats closed form position and orientation calculations until the solution is achieved. Several manipulators are used to demonstrate the solution methods: cartesian, cylindrical, spherical, and an anthropomorphic articulated arm, based on the Flight Telerobotic Servicer (FTS) arm. A singularity analysis is presented for the double universal joint wrist attached to the above manipulator arms. While the double universal joint wrist standing alone is singularity-free in orientation, the singularity analysis indicates the presence of coupled position/orientation singularities of the spherical and articulated manipulators with the wrist. The cartesian and cylindrical manipulators with the double universal joint wrist were found to be singularity-free. The methods of this paper can be implemented in a real-time controller for manipulators with the double universal joint wrist. Such mechanically dextrous systems could be used in telerobotic and industrial applications, but further work is required to avoid the singularities.

Williams, Robert L., III↗

NASA NTRS1975

A uniformly valid solution for motion about the interior libration point of the perturbed elliptic-restricted problem

Bounded motion about Lagrange collinear libration points is considered for a perturbed elliptic-restricted problem. A practical application is the motion of a satellite near a libration point collinear with the sun and the earth-moon barycenter. Such a system is treated here as an earth-sun-satellite elliptic restricted problem with lunar perturbations. The method of dual time scales is utilized to develop a uniformly valid three-dimensional analytical solution to the satellite's equations of motion. The analytical development applies somewhat generally to that class of four-body problems where the second primary mass is much greater than the first, and the third primary mass is much greater than the second.

Richardson, D. L.↗

NASA NTRS1985

Diffusion-convection function of cosmic rays

The fundamental properties and some numerical results of the solution of the diffusion equation of an impulsive cosmic-ray point source in an uniform, unbounded and spherically symmetrical moving medium is presented. The diffusion-convection(D-C) function is an elementary composite function of the solution of the D-C equation for the particles injected impulsively from a diffusive point source into the medium. It is the analytic solution derived by the dimensional method for the propagation equation of solar cosmic rays in the heliosphere, i.e. the interplanetary space. Because of the introduction of convection effect of solar wind, a nonhomogeneous term appears in the propagation equation, it is difficult to express its solution in terms of the ordinary special functions. The research made so far has led to a solution containing only the first order approximation of the convection effect.

Zhang, G.↗

NASA NTRS1989

Efficient estimation of diffusion during dendritic solidification

A very efficient finite difference method has been developed to estimate the solute redistribution during solidification with diffusion in the solid. This method is validated by comparing the computed results with the results of an analytical solution derived by Kobayashi (1988) for the assumptions of a constant diffusion coefficient, a constant equilibrium partition ratio, and a parabolic rate of the advancement of the solid/liquid interface. The flexibility of the method is demonstrated by applying it to the dendritic solidification of a Pb-15 wt pct Sn alloy, for which the equilibrium partition ratio and diffusion coefficient vary substantially during solidification. The fraction eutectic at the end of solidification is also obtained by estimating the fraction solid, in greater resolution, where the concentration of solute in the interdendritic liquid reaches the eutectic composition of the alloy.

Yeum, K. S.↗

NASA NTRS1996

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

Kaljevic, Igor↗

NASA NTRS1983

Simplified Modeling of Tetrahedral Trusses

Tetrahedral-lattice structure exhibits great complexity in mathematical analysis, and complexity quickly increases with number of elements. Solutions for static and dynamic behavior of structure obtained with model consisting of analogous structure with fewer repeating elements. Advantage of continuum models is analytic solutions can be obtained, either in closed form or in relatively-tractable Taylor series.

Leondis, A.↗

NASA NTRS1982

Numerical computation of exponential matrices using the Cayley-Hamilton theorem

A method for computing exponential matrices, which often arise naturally in the solution of systems of linear differential equations, is developed. An exponential matrix is generated as a linear combination of a finite number (equal to the matrix order) of matrices, the coefficients of which are scalar infinite sums. The method can be generalized to apply to any formal power series of matrices. Attention is focused upon the exponential function, and the matrix exponent is assumed tri-diagonal in form. In such cases, the terms in the coefficient infinite sums can be extracted, as recursion relations, from the characteristic polynomial of the matrix exponent. Two numerical examples are presented in some detail: (1) the three dimensional infinitesimal rotation rate matrix, which is skew symmetric, and (2) an N-dimensional tri-diagonal and symmetric finite difference matrix which arises in the numerical solution of the heat conduction partial differential equation. In the second example, the known eigenvalues and eigenvectors of the finite difference matrix permit an analytical solution for the exponential matrix, through the theory of diagonalization and similarity transformations, which is used for independent verification. The convergence properties of the scalar infinite summations are investigated for finite difference matrices of various orders up to ten, and it is found that the number of terms required for convergence increases slowly with the order of the matrix.

Walden, H.↗

NASA NTRS1984

Skin-stiffener interface stresses in composite stiffened panels

A model and solution method for determining the normal and shear stresses in the interface between the skin and the stiffener attached flange were developed. An efficient, analytical solution procedure was developed and incorporated in a sizing code for stiffened panels. The analysis procedure described provides a means to study the effects of material and geometric design parameters on the interface stresses. These stresses include the normal stress, and the shear stresses in both the longitudinal and the transverse directions. The tendency toward skin/stiffener separation may therefore be minimized by choosing appropriate values for the design variables. The most important design variables include the relative bending stiffnesses of the skin and stiffener attached flange, the bending stiffness of the stiffener web, and the flange width. The longitudinal compressive loads in the flange and skin have significant effects on the interface stresses.

Wang, J. T. S.↗

NASA NTRSDate not supplied

Convolution-Based Numerical Solutions of Transient Temperature Fields during Powder Bed Fusion Additive Manufacturing: Theory, Accuracy, and Computational Cost

Powder bed fusion (PBF) additive manufacturing has found numerous applications in the aerospace domain. However, components fabricated via PBF have a complex time-temperature history that significantly impacts subsequent mechanical performance. This study examines convolution-based numerical solutions of transient temperature fields that support simulations involving arbitrary beam shapes and paths during PBF. The convolutional approach is verified through comparisons with analytical solutions of the temperature field. The computational speed and accuracy of the method are assessed through comparisons with other explicit and implicit numerical techniques. In addition, the straightforward translation of the approach from a CPU to a GPU implementation and the resultant performance improvement are presented. The role of the technique in predicting microstructure evolution during PBF (for a greater process-structure-property-performance framework) is also demonstrated. This work supports the development of computational materials methods for understanding and controlling the time-temperature history during PBF.

powder bed fusion↗

NASA NTRS1972

N-burn optimal analytic trajectories.

Derivation of N-burn analytic solutions for propellant-optimal transfer trajectories of a vehicle in a vacuum between arbitrary boundary conditions. Variational changes in the desired boundary conditions are expressed, in general, in terms of variational changes in the control vector and in the initial state vector. All coefficient matrices are computed recursively in terms of the analytic matrices established from the subarcs of the N-burn solution. The solution is applicable to shuttle ascent (exoatmospheric), rendezvous, and deorbit problems. Consideration is also given to state-variable and control-variable inequality constraints.

Jezewski, D. J.↗

NASA NTRS1974

Automatic numerical generation of body-fitted curvilinear coordinate system for field containing any number of arbitrary two-dimensional bodies

A method for automatic numerical generation of a general curvilinear coordinate system with coordinate lines coincident with all boundaries of a general multi-connected region containing any number of arbitrarily shaped bodies is presented. With this procedure the numerical solution of a partial differential system may be done on a fixed rectangular field with a square mesh with no interpolation required regardless of the shape of the physical boundaries, regardless of the spacing of the curvilinear coordinate lines in the physical field, and regardless of the movement of the coordinate system. Numerical solutions for the lifting and nonlifting potential flow about Joukowski and Karman-Trefftz airfoils using this coordinate system generation show excellent comparison with the analytic solutions. The application to fields with multiple bodies is illustrated by a potential flow solution for multiple airfoils.

Thompson, J. F.↗

NASA NTRS1985

A leaky magnetohydrodynamic waveguide model for the acceleration of high-speed solar wind streams in coronal holes

It is now reasonably well established that there is a correlation between high-speed solar wind streams and coronal holes. It has been concluded that a significant addition of momentum and/or energy in the region of supersonic flow is needed to explain the observed particle flux and flow speed observed in the high-speed streams. The most likely source of this additional momentum appears to be magnetohydrodynamic (MHD) waves propagating up from the solar surface. The present investigation is concerned with the propagation of MHD waves in a structure of finite transverse size, taking into account the consequences for the acceleration of high-speed solar wind streams. A waveguide solution for a model coronal hole is described, giving attention to a geometric or ray analysis of the slab waveguide, a wave mode analysis, an analytic solution of the dispersion relation for high-frequency waves, and the calculation of the time-averaged wave force.

Davila, J. M.↗

NASA NTRS2000

Aerodynamic Lift and Moment Calculations Using a Closed-Form Solution of the Possio Equation

In this paper, we present closed-form formulas for the lift and moment coefficients of a lifting surface in two dimensional, unsteady, compressible, subsonic flow utilizing a newly developed explicit analytical solution of the Possio equation. Numerical calculations are consistent with previous numerical tables based on series expansions or ad hoc numerical schemes. More importantly, these formulas lend themselves readily to flutter analysis, compared with the tedious table-look-up schemes currently in use.

Lin, Jensen↗

NASA NTRS1979

Computations of unsteady transonic flow governed by the conservative full potential equation using an alternating direction implicit algorithm

A development was the time linearization of the density function. This linearization reduces the solution process from solving just a single equation. Two sample cases were computed. First, a one dimensional traveling shock wave was computed and compared with the analytic solution. Second, a two dimensional case was calculated for a flow field which resulted from a thickening and subsequently, thinning airfoil. The resulting flow field, which included a traveling shock wave, was compared to the flow field obtained from the low frequency, small disturbance, transonic equation.

Goorjian, P. M.↗

NASA NTRS1990

Finite element two-dimensional panel flutter at high supersonic speeds and elevated temperature

A consistent finite element formulation for flutter of two-dimensional panels including effects of temperature variation along panel length is presented. The von Karman nonlinear strain-displacement relation is used to account for the large-deflections, and the quasi-steady first-order piston theory is employed for aerodynamic loading in the finite-element formulation. The panel flutter under combined thermal/aerodynamic loading can be separated into solving two problems in sequence. They are: (1) postbuckling deflection, stresses and buckling stability boundary, and (2) critical dynamic pressure and linear flutter boundary. The solution procedures in solving these problems are presented in detail. Finite element results of two-dimensional panels with uniform temperature change are compared with classic analytical solutions. Influence of nonuniform temperature distributions on critical buckling temperature, flutter deflection and stresses, and flutter boundary is also given.

Xue, David Y.↗