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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Simulated measurement of power flow in structures near to simple sources and simple boundaries

Advances in electronics technology along with the advent of low cost multichannel Fast Fourier analyzers have made it practical to use higher order central difference formulas to measure power flow in 1- and 2-D structures. The method discussed uses five point differencing for the spatial derivatives in 1-D and a thirteen point difference pattern for the spatial derivatives in 2-D plates and shells. It is assumed that the measuring transducers are accelerometers. An analytical study of the higher order differencing method and the conventional two accelerometer method was performed as a preliminary to the application of these methods to actual aircraft structures. Some classical problems were analyzed in order to simulate and compare the performance of the two methods under near field measurement conditions. These near field conditions include examples of power flows near simple sources and simple boundaries. The estimates produced by the two methods were compared to the exact solution in each example. Presented are the theory and selected results of the study, which indicate that the bias errors of the two accelerometer method under near field measurement conditions may be much larger than previous studies have suggested.

Mcgary, Michael C.↗
NASA NTRS1988

Applications of an exponential finite difference technique

An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.

Handschuh, Robert F.↗
NASA NTRS1987

On the theory of the synthesis of single and dual offset shaped reflector antennas

Two reflectors shaped so as to transforn a geometrical optics (GO) feed power pattern into a uniform amplitude and phase aperture distribution are indicated schematically. Single- and dual-shaped geometrical optics problems are solved via exact solutions obtained by progressive integration. It is concluded that an infinite set of such solutions exists.

Galindo-Israel, Victor↗
NASA NTRS1987

Nonlinear stability of a flow with bound eddies

Kovasznay (1948) obtained an exact solution of the Navier-Stokes equations that describes a flow with periodic bound eddies. The stability of this nonparallel flow with respect to three-dimensional finite amplitude disturbances is analyzed by use of the energy method. A sufficient condition for global stability is obtained. It is shown that the additional nonhomogeneity in shear stress distribution in a nonparallel flow is a destabilizing factor.

Lin, S. P.↗
NASA NTRS1988

Validation of the NESSUS probabilistic finite element analysis computer program

A computer program, NESSUS, is being developed as part of a NASA-sponsored project to develop probabilistic structural analysis methods for propulsion system components. This paper describes the process of validating the NESSUS code, as it has been developed to date, and presents numerical results comparing NESSUS and exact solutions for a set of selected problems.

Wu, Y.-T.↗
NASA NTRS1988

A composite L-band HH radar backscattering model for coniferous forest stands

The radar backscattering model developed by Richards et al. (1987), has been improved and further tested in this research. The trunk term may now be calculated from the exact solution to the electromagnetic wave equations instead of the corner reflector equation. Rough surface models have been introduced into the radar model, so that the forward reflectance and the backscattering from the ground surface are now calculated from the same model and, thus, are consistent. The number of trees in an individual pixel is assumed to be Poisson distributed, with tree height in a stand log-normally distributed. The simulated results show that the match of backscattering coefficients for eight forest stands between SIR-B image data and the simulated results are satisfying, and that the trunk term now seems to be convincingly established as the dominant term in the L-band HH radar return from coniferous forest stands.

Sun, Guoquing↗
NASA NTRS1988

Weak-wave analysis of shock interaction with a slipstream

A weak wave analysis of shock interaction with a slipstream is presented. The theory is compared to that for the acoustic case and to the exact nonlinear analysis. Sample calculations indicate that the weak wave theory yields a good approximation to the exact solution when the shock waves are sufficiently weak that the associated entropy increase is negligible. A qualitative discussion of the case of counterflowing streams is also included.

Barger, Raymond L.↗
NASA NTRS1988

A numerical study of ENO and TVD schemes for shock capturing

The numerical performance of a second-order upwind-based total variation diminishing (TVD) scheme and that of a uniform second-order essentially non-oscillatory (ENO) scheme for shock capturing are compared. The TVD scheme used is a modified version of Liou, using the flux-difference splitting (FDS) of Roe and his superbee function as the limiter. The construction of the basic ENO scheme is based on Harten, Engquist, Osher, and Chakravarthy, and the 2-D extensions are obtained by using a Strang-type of fractional-step time-splitting method. Numerical results presented include both steady and unsteady, 1-D and 2-D calculations. All the chosen test problems have exact solutions so that numerical performance can be measured by comparing the computer results to them. For 1-D calculations, the standard shock-tube problems of Sod and Lax are chosen. A very strong shock-tube problem, with the initial density ratio of 400 to 1 and pressure ratio of 500 to 1, is also used to study the behavior of the two schemes. For 2-D calculations, the shock wave reflection problems are adopted for testing. The cases presented in this report include flows with Mach numbers of 2.9, 5.0, and 10.0.

Chang, Shih-Hung↗
NASA NTRS1988

Singular perturbations and time scales in the design of digital flight control systems

The results are presented of application of the methodology of Singular Perturbations and Time Scales (SPATS) to the control of digital flight systems. A block diagonalization method is described to decouple a full order, two time (slow and fast) scale, discrete control system into reduced order slow and fast subsystems. Basic properties and numerical aspects of the method are discussed. A composite, closed-loop, suboptimal control system is constructed as the sum of the slow and fast optimal feedback controls. The application of this technique to an aircraft model shows close agreement between the exact solutions and the decoupled (or composite) solutions. The main advantage of the method is the considerable reduction in the overall computational requirements for the evaluation of optimal guidance and control laws. The significance of the results is that it can be used for real time, onboard simulation. A brief survey is also presented of digital flight systems.

Naidu, Desineni S.↗
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.↗