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 235 records · Page 13

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↗

Numerical Algorithms Based on Biorthogonal Wavelets

Wavelet bases are used to generate spaces of approximation for the resolution of bidimensional elliptic and parabolic problems. Under some specific hypotheses relating the properties of the wavelets to the order of the involved operators, it is shown that an approximate solution can be built. This approximation is then stable and converges towards the exact solution. It is designed such that fast algorithms involving biorthogonal multi resolution analyses can be used to resolve the corresponding numerical problems. Detailed algorithms are provided as well as the results of numerical tests on partial differential equations defined on the bidimensional torus.

Ponenti, Pj.↗

Numerical Algorithms Based on Biorthogonal Wavelets

Wavelet bases are used to generate spaces of approximation for the resolution of bidimensional elliptic and parabolic problems. Under some specific hypotheses relating the properties of the wavelets to the order of the involved operators, it is shown that an approximate solution can be built. This approximation is then stable and converges towards the exact solution. It is designed such that fast algorithms involving biorthogonal multi resolution analyses can be used to resolve the corresponding numerical problems. Detailed algorithms are provided as well as the results of numerical tests on partial differential equations defined on the bidimensional torus.

Ponenti, Pj.↗

Application of NASA General-Purpose Solver to Large-Scale Computations in Aeroacoustics

Of several iterative and direct equation solvers evaluated previously for computations in aeroacoustics, the most promising was the NASA-developed General-Purpose Solver (winner of NASA's 1999 software of the year award). This paper presents detailed, single-processor statistics of the performance of this solver, which has been tailored and optimized for large-scale aeroacoustic computations. The statistics, compiled using an SGI ORIGIN 2000 computer with 12 Gb available memory (RAM) and eight available processors, are the central processing unit time, RAM requirements, and solution error. The equation solver is capable of solving 10 thousand complex unknowns in as little as 0.01 sec using 0.02 Gb RAM, and 8.4 million complex unknowns in slightly less than 3 hours using all 12 Gb. This latter solution is the largest aeroacoustics problem solved to date with this technique. The study was unable to detect any noticeable error in the solution, since noise levels predicted from these solution vectors are in excellent agreement with the noise levels computed from the exact solution. The equation solver provides a means for obtaining numerical solutions to aeroacoustics problems in three dimensions.

Watson, Willie R.↗

Determination of stress intensity factors for interface cracks under mixed-mode loading

A simple technique was developed using conventional finite element analysis to determine stress intensity factors, K1 and K2, for interface cracks under mixed-mode loading. This technique involves the calculation of crack tip stresses using non-singular finite elements. These stresses are then combined and used in a linear regression procedure to calculate K1 and K2. The technique was demonstrated by calculating three different bimaterial combinations. For the normal loading case, the K's were within 2.6 percent of an exact solution. The normalized K's under shear loading were shown to be related to the normalized K's under normal loading. Based on these relations, a simple equation was derived for calculating K1 and K2 for mixed-mode loading from knowledge of the K's under normal loading. The equation was verified by computing the K's for a mixed-mode case with equal and normal shear loading. The correlation between exact and finite element solutions is within 3.7 percent. This study provides a simple procedure to compute K2/K1 ratio which has been used to characterize the stress state at the crack tip for various combinations of materials and loadings. Tests conducted over a range of K2/K1 ratios could be used to fully characterize interface fracture toughness.

Naik, Rajiv A.↗

On the equilibrium of a cylindrical plasma supported horizontally by magnetic fields in uniform gravity

The mechanical equilibrium of a cylinder of plasma suspended horizontally by magnetic fields in uniform gravity is examined. A set of exact solutions describing the equilibrium is derived assuming the plasma distribution to be cylindrically symmetric to obtain an exact force balance between plasma pressure, the Lorentz pressure, and gravity in space. The set of solutions considers a case of uniform temperature and cases where the temperature rises from zero at the center of the plasma cylinder to reach a constant asymptotic value outside the cylinder.

Lerche, I.↗

Finite element solutions of free surface flows

A procedure is presented for using NASTRAN to determine the flow field about arbitrarily shaped bodies in the presence of a free surface. The fundamental unknown of the problem is the velocity potential which must satisfy Laplace's equation in the fluid region. Boundary conditions on the free surface may involve second order derivatives in space and time. In cases involving infinite domains either a tractable radiation condition is applied at a truncated boundary or a series expansion is used and matched to the local finite elements. Solutions are presented for harmonic, transient, and steady state problems and compared to either exact solutions or other numerical solutions.

Zarda, P. R.↗

Effects of boundary layer refraction and fuselage scattering on fuselage surface noise from advanced turboprop propellers

An acoustic disturbance's propagation through a boundary layer is discussed with a view to the analysis of the acoustic field generated by a propfan rotor incident to the fuselage of an aircraft. Applying the parallel flow assumption, the resulting partial differential equations are reduced to an ordinary acoustic pressure differential equation by means of the Fourier transform. The methods used for the solution of this equation include those of Frobenius and of analytic continuation; both yield exact solutions in series form. Two models of the aircraft fuselage-boundary layer system are considered, in the first of which the fuselage is replaced by a flat plate and the acoustic field is assumed to be two-dimensional, while in the second the fuselage is a cylinder in a fully three-dimensional acoustic field. It is shown that the boundary layer correction improves theory-data comparisons over simple application of a pressure-doubling rule at the fuselage.

Mcaninch, G. L.↗

Solutions of the benchmark problems by the dispersion-relation-preserving scheme

The 7-point stencil Dispersion-Relation-Preserving scheme of Tam and Webb is used to solve all the six categories of the CAA benchmark problems. The purpose is to show that the scheme is capable of solving linear, as well as nonlinear aeroacoustics problems accurately. Nonlinearities, inevitably, lead to the generation of spurious short wave length numerical waves. Often, these spurious waves would overwhelm the entire numerical solution. In this work, the spurious waves are removed by the addition of artificial selective damping terms to the discretized equations. Category 3 problems are for testing radiation and outflow boundary conditions. In solving these problems, the radiation and outflow boundary conditions of Tam and Webb are used. These conditions are derived from the asymptotic solutions of the linearized Euler equations. Category 4 problems involved solid walls. Here, the wall boundary conditions for high-order schemes of Tam and Dong are employed. These conditions require the use of one ghost value per boundary point per physical boundary condition. In the second problem of this category, the governing equations, when written in cylindrical coordinates, are singular along the axis of the radial coordinate. The proper boundary conditions at the axis are derived by applying the limiting process of r approaches 0 to the governing equations. The Category 5 problem deals with the numerical noise issue. In the present approach, the time-independent mean flow solution is computed first. Once the residual drops to the machine noise level, the incident sound wave is turned on gradually. The solution is marched in time until a time-periodic state is reached. No exact solution is known for the Category 6 problem. Because of this, the problem is formulated in two totally different ways, first as a scattering problem then as a direct simulation problem. There is good agreement between the two numerical solutions. This offers confidence in the computed results. Both formulations are solved as initial value problems. As such, no Kutta condition is required at the trailing edge of the airfoil.

Tam, Christopher K. W.↗

Cross-Sectional Deformations of Monocoque Beams and Their Effects on the Natural Vibration Frequencies

The variational principle, differential equations, and boundary conditions governing the cross-sectional distortions due to inertia loading of a two-dimensional model of a thin monocoque wing are shown. A theoretical analysis of this simplified model is made in order to determine the nature of the coupling between the cross-sectional modes and the spanwise deformation modes. General solutions are obtained in finite-difference form for arbitrary cross sections and an exact solution is presented for a parabolic-arc cross section of constant cover thickness. The application of these results in evaluating the coupled frequencies of the actual structure is discussed. Frequencies evaluated for a parabolic-arc monocoque beam show good agreement with experimental values.

Thomson, Robert G.↗

Solution sensitivity and accuracy study of NASTRAN for large dynamic problems involving structural damping

Large dynamic problems involving NASTRAN SOLUTION 8 (i.e., the steady state dynamic response option wherein all response quantities vary as e sub i omega t, where omega is the driving frequency and t is time) are considered. Using a submerged steel plate with a viscoelastic layer layer as the bench mark sample, the solution sensitivity and solution accuracy are checked. The solution sensitivity is examined by running the same finite element model on different computers, different versions of NASTRAN, and different precision levels. The solution accuracy is evaluated for these same runs by comparing the NASTRAN results with the exact solution of the same problem.

Kalinowski, A. J.↗

On the dynamics of approximating schemes for dissipative nonlinear equations

Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.

Jones, Donald A.↗

Composite laminates with spatially varying fiber orientations - 'Variable stiffness panel concept'

A solution has been obtained to the plane elasticity problem for a symmetrically laminated composite panel with spatially varying fiber orientations. Since variation of the fiber angles along the length of a composite laminate results in stiffness properties that change as a function of location, the laminates are called variable stiffness panels. An analysis of the stiffness variation and its effect on the elastic response of the panel is presented here. A numerical solution has been obtained using an iterative collocation technique. Corresponding closed-form solutions are given for three different sets of boundary conditions. Two of the cases considered have exact solutions and thus serve to validate the numerical model.

Guerdal, Zafer↗

Structural response to discrete and continuous gusts of an airplane having wing bending flexibility and a correlation of calculated and flight results

An analysis is made of the structural response to gusts of an airplane having the degrees of freedom of vertical motion and wing bending flexibility and basic parameters are established. A convenient and accurate numerical solution of the response equations is developed for the case of discrete-gust encounter, an exact solution is made for the simpler case of continuous-sinusoidal-gust encounter, and the procedure is outlined for treating the more realistic condition of continuous random atmospheric turbulence, based on the methods of generalized harmonic analysis. Correlation studies between flight and calculated results are then given to evaluate the influence of wing bending flexibility on the structural response to gusts of two twin-engine transports and one four-engine bomber. It is shown that calculated results obtained by means of a discrete-gust approach reveal the general nature of the flexibility effects and lead to qualitative correlation with flight results. In contrast, calculations by means of the continuous-turbulence approach show good quantitative correlation with flight results and indicate a much greater degree of resolution of the flexibility effects.

Houbolt, John C↗

The method of Ritz applied to the equation of Hamilton

Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.

Bailey, C. D.↗

Buckling and vibration of periodic lattice structures

Lattice booms and platforms composed of flexible members or large diameter rings which may be stiffened by cables in order to support membrane-like antennas or reflector surfaces are the main components of some large space structures. The nature of these structures, repetitive geometry with few different members, makes possible relatively simple solutions for buckling and vibration of a certain class of these structures. Each member is represented by a stiffness matrix derived from the exact solution of the beam column equation. This transcendental matrix gives the current member stiffness at any end load or frequency. Using conventional finite element techniques, equilibrium equations can be written involving displacements and rotations of a typical node and its neighbors. The assumptions of a simple trigonometric mode shape is found to satisfy these equations exactly; thus the entire problem is governed by just one 6 x 6 matrix equation involving the amplitude of the displacement and rotation mode shapes. The boundary conditions implied by this solution are simple supported ends for the column type configurations.

Anderson, M. S.↗

A theoretical investigation of the sound radiation fields associated with a Bellmouth inlet

Analytical results are obtained by numerical integration of a cylindrically symmetric integral representation of the external solutions of the Helmholtz equation. The accuracy of this method is checked by comparisons of computed results with 'exact' solutions generated by the point source method. In all cases, the average error for the amplitude and phase of the points calculated in the field is found to be less than ten percent. Theoretical studies which model experiments run for the NASA Langely Bellmouth inlet configuration are presented and comparisons are made with the experimental results. In all cases very good agreement is obtained between the experimental and theoretically calculated values.

Meyer, W. L.↗

Static Load Distribution in Ball Bearings

A numerical procedure for computing the internal loading distribution in statically loaded, single-row, angular-contact ball bearings when subjected to a known combined radial and thrust load is presented. The combined radial and thrust load must be applied in order to avoid tilting between inner and outer rings. The numerical procedure requires the iterative solution of Z + 2 simultaneous nonlinear equations - where Z is the number of the balls - to yield an exact solution for axial and radial deflections, and contact angles. Numerical results for a 218 angular-contact ball bearing have been compared with those from the literature and show significant differences in the magnitudes of the ball loads, contact angles, and the extent of the loading zone.

Ricci, Mario↗