Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXACT 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Scattering by a two dimensional groove in a ground plane

Higher order boundary conditions involve derivatives of the fields beyond the first and were recently shown to be more effective than traditional first order conditions in modeling dielectric coatings and layers. In this report an application of a third order generalized boundary condition to scattering by a filled rectangular groove is presented. Deficiencies of such higher order boundary conditions are addressed and a correction is proposed for the present case. As part of the process of examining and improving the accuracy of the proposed generalized boundary conditions, an exact solution is developed and a comparison is provided with a solution based on the standard impedance boundary condition.

Barkeshli, K.↗
NASA NTRS1989

Crack-face displacements for embedded elliptic and semi-elliptical surface cracks

Analytical expressions for the crack-face displacements of an embedded elliptic crack in infinite solid subjected to arbitrary tractions are obtained. The tractions on the crack faces are assumed to be expressed in a polynomial form. These displacements expressions complete the exact solution of Vijayakumar and Atluri, and Nishioki and Atluri. For the special case of an embedded crack in an infinite solid subjected to uniform pressure loading, the present displacements agree with those by Green and Sneddon. The displacement equations derived were used with the finite-element alternating method (FEAM) for the analysis of a semi-elliptic surface crack in a finite solid subjected to remote tensile loading. The maximum opening displacements obtained with FEAM are compared to those with the finite-element method with singularity elements. The maximum crack opening displacements by the two methods showed good agreement.

Raju, I. S.↗
NASA NTRS1989

Optimal placement of excitations and sensors by simulated annealing

The optimal placement of discrete actuators and sensors is posed as a combinatorial optimization problem. Two examples for truss structures were used for illustration; the first dealt with the optimal placement of passive dampers along existing truss members, and the second dealt with the optimal placement of a combination of a set of actuators and a set of sensors. Except for the simplest problems, an exact solution by enumeration involves a very large number of function evaluations, and is therefore computationally intractable. By contrast, the simulated annealing heuristic involves far fewer evaluations and is best suited for the class of problems considered. As an optimization tool, the effectiveness of the algorithm is enhanced by introducing a number of rules that incorporate knowledge about the physical behavior of the problem. Some of the suggested rules are necessarily problem dependent.

Salama, Moktar↗
NASA NTRS1989

Control of the errors of discretization and idealization in finite element analysis

Understanding of the basic principles which control errors of discretization in finite element analysis has increased very substantially since 1980. The main milestones were: (1) development of the theoretical basis of p-extensions (1981); (2) understanding of the proper interplay between mesh design and assignment of polynomial degree to elements. Practical realization of exponential convergence rates, independently of the smoothness of the exact solution (1984); and (3) industrial experience with the new finite element technology known as the p- or hp-version of the finite element method: General Dynamics reported thirty- to forty-fold savings in terms of human time and large savings in computer time (1986). Lockheed reported favorably on their evaluation of error estimation and quality control capabilities of the p-version in industrial settings (1987). The gains in our understanding of how to control the errors of discretization represent only half of the control necessary to ensure that a numerical model is in fact an accurate representation of the corresponding physical system. Control of the errors of idealization is equally important. A brief overview of the main ideas of how to ensure the quality and reliability of mathematical models of structural systems is presented.

Szabo, Barna A.↗
NASA NTRS1988

Methods for analysis of cracks in three-dimensional solids

Various methods used for determining stress-intensity factors for cracked three-dimensional bodies are reviewed. The review is limited to the determination of mode-I stress-intensity factors. Some exact solutions for cracks in infinite solids are presented, and various approximate methods that have been used for the solution of the boundary-value problem of finite solids with cracks are addressed. The techniques used to extract the stress-intensity factors from these solutions are considered, and various methods used to calculate the stress-intensity factors for through-the-thickness cracks, semielliptical surface cracks, and quarter-elliptical corner cracks at holes are compared for the case of remote tensile loading.

Raju, I. S.↗
NASA NTRS1989

Assessment of numerical techniques for unsteady flow calculations

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

Hsieh, Kwang-Chung↗
NASA NTRS1989

A weak Hamiltonian finite element method for optimal guidance of an advanced launch vehicle

A temporal finite-element method based on a mixed form of the Hamiltonian weak principle is presented for optimal control problems. The mixed form of this principle contains both states and costates as primary variables, which are expanded in terms of nodal values and simple shape functions. Time derivatives of the states and costates do not appear in the governing variational equation; the only quantities whose time derivatives appear therein are virtual states and virtual costates. Numerical results are presented for an elementary trajectory optimization problem; they show very good agreement with the exact solution along with excellent computational efficiency and self-starting capability. The feasibility of this approach for real-time guidance applications is evaluated. A simplified model for an advanced launch vehicle application that is suitable for finite-element solution is presented.

Hodges, Dewey H.↗
NASA NTRS1990

Dynamics of column stability with partial end restraints

The dynamic behavior of columns with partial end restraints and loads consisting of a dead load and a pulsating load are investigated. The differential equation is solved using a lumped impulse recurrence formula relative to time coupled with a finite difference discretization along the member length. A computer program is written from which the first critical frequencies are found as a function of end stiffness. The case of a pinned ended column compares very well with the exact solution. Also, the natural frequency and buckling load formulas are derived for equal and unequal end restraints.

Gregory, Peyton B.↗
NASA NTRS1990

Differential equation based method for accurate approximations in optimization

A method to efficiently and accurately approximate the effect of design changes on structural response is described. The key to this method is to interpret sensitivity equations as differential equations that may be solved explicitly for closed form approximations, hence, the method is denoted the Differential Equation Based (DEB) method. Approximations were developed for vibration frequencies, mode shapes and static displacements. The DEB approximation method was applied to a cantilever beam and results compared with the commonly-used linear Taylor series approximations and exact solutions. The test calculations involved perturbing the height, width, cross-sectional area, tip mass, and bending inertia of the beam. The DEB method proved to be very accurate, and in most cases, was more accurate than the linear Taylor series approximation. The method is applicable to simultaneous perturbation of several design variables. Also, the approximations may be used to calculate other system response quantities. For example, the approximations for displacements are used to approximate bending stresses.

Pritchard, Jocelyn I.↗
NASA NTRS1990

Concerning the interaction of non-stationary cross-flow vortices in a three-dimensional boundary layer

Recently there has been much work devoted to considering some of the many and varied interaction mechanisms which may be operative in three-dimensional boundary layer flows. This paper is concerned with resonant triads of crossflow vortices. The effects of interactions upon resonant triads is examined where each member of the triad has the property of being linearly neutrally stable so that the importance of the interplay between modes can be relatively easily assessed. Modes within the boundary layer flow above a rotating disc are investigated because of the similarity between this disc flow and many important practical flows and, secondly, because the selected flow is an exact solution of the Navier-Stokes equations which makes its theoretical analysis especially attractive. It is demonstrated that the desired triads of linearly neutrally stable modes can exist within the chosen boundary layer flow. Evolution equations are obtained to describe the development of the amplitudes of these modes once the interaction mechanism is accounted for. It is found that the coefficients of the interaction terms within the evolution equations are, in general, given by quite intricate expressions although some elementary numerical work shows that the evaluation of these coefficients is practicable. The basis of the work lends itself to generalization to more complicated boundary layers, and effects of detuning or non-parallelism could be provided for within the asymptotic framework.

Bassom, Andrew P.↗
NASA NTRS1990

Radiative configuration factors from cylinders to coaxial axisymmetric bodies

Exact solutions are obtained for the radiative configuration factor between differential elements of arbitrary orientation and cylinders. A general method is then proposed for calculating the view factor from a cylinder to a coaxial axisymmetric body using only a single numerical integration. The method is illustrated for axisymmetric bodies with function generators described by a power law equation. The method may be useful in calculating radiative heat transfer between cylindrical bodies and high-density exhaust gases and between annular radiative fins and their bases.

Saltiel, C.↗
NASA NTRS1990

Crack-path effect on material toughness

The main features of a toughening mechanism associated with a curvilinear crack path are examined using a model consisting of a macrocrack in a brittle solid with a curvilinear segment at the crack tip. A numerical procedure for finite and semiinfinite cracks is formulated and evaluated using an example which has an exact solution (a finite crack in the form of a circular arc in a uniform stress field). It is shown that, for a relatively small amplitude of crack path oscillations, the toughening ratio can be taken equal to the ratio of the corresponding crack path lengths.

Rubinstein, Asher A.↗
NASA NTRS1990

Differential equation based method for accurate approximations in optimization

This paper describes a method to efficiently and accurately approximate the effect of design changes on structural response. The key to this new method is to interpret sensitivity equations as differential equations that may be solved explicitly for closed form approximations, hence, the method is denoted the Differential Equation Based (DEB) method. Approximations were developed for vibration frequencies, mode shapes and static displacements. The DEB approximation method was applied to a cantilever beam and results compared with the commonly-used linear Taylor series approximations and exact solutions. The test calculations involved perturbing the height, width, cross-sectional area, tip mass, and bending inertia of the beam. The DEB method proved to be very accurate, and in msot cases, was more accurate than the linear Taylor series approximation. The method is applicable to simultaneous perturbation of several design variables. Also, the approximations may be used to calculate other system response quantities. For example, the approximations for displacement are used to approximate bending stresses.

Pritchard, Jocelyn I.↗
NASA NTRS1990

Ballooning mode instability at the plasmapause

The ballooning mode instability, which can excite hydromagnetic waves at the plasmapause, is studied in the presence of azimuthal plasma flows induced during geomagnetically disturbed periods. A general sufficient criterion for the ballooning mode stability is derived, for a change in the potential energy greater than or equal to 0, which involves the integration over an entire field line. A local stability analysis at the equatorial plasmapause region shows that the ballooning modes could be spontaneously generated via instability under at least two conditions: one is similar to the usual interchange condition, and the second to the quasi-interchange modes. Both of these local instability conditions can be derived from the general stability criterion. Finally an exact solution for the equilibrium state with flow is derived analytically, and the change in the potential energy is computed numerically. It is found that, in the cases studied, the flow does not spontaneously excite the ballooning modes; it only further stabilizes (or destabilizes) the ballooning spectrum if originally the system is stable (or unstable). The analysis would be useful for the interpretation of some of the low-frequency modes observed at the ground and near the equatorial plasmapause.

Lakhina, G. S.↗
NASA NTRS1990

Crack-face displacements for an embedded elliptic crack in an infinite solid subjected to arbitrary tractions

Analytical expression for the crack-face displacements of an embedded elliptic crack in infinite solid subjected to arbitrary traction are obtained. The tractions on the crack faces are assumed to be expressed in a polynomial form. These displacements expressions complete the exact solution of Vijayakumar and Atluri, and Nishioki and Atluri. For the special case of an embedded crack in an infinite solid subjected to uniform pressure loading, the present displacements agree with those by Green and Sneddon. The displacement equations derived were used with the finite-element alternating method (FEAM) for the analysis of a semi-elliptic surface crack in a finite solid subjected to remote tensile loading. The maximum opening displacements obtained with FEAM are compared to those with the finite-element method with singularity elements. The maximum crack opening displacements by the two methods showed good agreement.

Raju, I. S.↗
NASA NTRS1990

The free-wake computation of rotor-body flows

In this paper a method is described for predicting the compressible, free-wake, flow about a lifting rotor-body configuration. The method is an extension of a unique vorticity embedded full-potential method used to calculate free-wake rotor hover performance. An unusual feature of this method is that it obviates the requirement for multiple grids to treat the rotor-body problem. The approach used to treat the body is similar to that used to include the rotor wake in the full potential calculation. The body is modeled as a structured circulation sheet and the strength of this sheet is determined in an iterative manner. Initially the method is tested to compute the flow past simple isolated bodies like cylinders and spheres. After a comparison of these simple computations with exact solutions this procedure is included into the HELIX-I, free-wake rotor code, to compute the flow around a rotor mounted on a large whirl tower. The effects of the tower on rotor wake geometry and load distribution are presented.

Ramachandran, K.↗
NASA NTRS1990

TE scattering by a two-dimensional groove in a ground plane using higher order boundary conditions

Higher order boundary conditions involve derivatives of the fields beyond the first and were recently shown to be more effective than traditional first order conditions in modeling dielectric coatings and layers. In this report an application of a third order generalized boundary condition to scattering by a filled rectangular groove is presented. Deficiencies of such higher-order boundary conditions are addressed and a correction is proposed for the present case. As part of the process of examining and correcting the accuracy of the proposed generalized boundary conditions, an exact solution is developed and a comparison is provided with a solution based on the standard impedance boundary condition.

Barkeshli, Kasra↗
NASA NTRS1991

Programs For Modeling Fault-Tolerant Computing Systems

Pade Approximation with Scaling, (PAWS) and Scaling Taylor Exponential Matrix (STEM) computer programs are software tools for design and validation. Provide flexible, user-friendly, language-based interface for input of Markov mathematical methods describing behaviors of fault-tolerant computer systems. Markov models include both recovery from faults via reconfiguration and behaviors of such systems when faults occur. PAWS and STEM produce exact solutions of probability of system failure and provide conservative estimate of number of significant digits in solution. Written in PASCAL and FORTRAN.

Butler, Ricky W.↗