Search NASA⌕ Search

SEARCH · Search NASA

Results for “collocation”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

An analytical and experimental stress analysis of a practical mode II fracture-test specimen

A boundary collocation method has been employed to determine the Mode II stress intensity factors for a pair of through-the-thickness edge cracks in a finite isotropic plate. An elastostatic analysis has been carried out in terms of the complete Williams stress function employing both even and odd components. The results of the numerical analysis were verified by a two-step procedure whereby the symmetric and antisymmetric portions of the solution were independently compared with existing solutions. The complete solution was verified by comparison with a photoelastic analysis. A compact shear specimen (CSS) of Hysol epoxy resin was loaded in a photoelastic experiment designed to study the isochromatic fringe patterns resulting from the Mode II crack tip stress distribution. The experiment verified that a pure Mode II stress distribution existed in the neighborhood of the crack tips and confirmed the accuracy of the boundary collocation solution for the Mode II stress intensity factors.

Chisholm, D. B.↗
NASA NTRS1977

Polynomial interpolation methods for viscous flow calculations

Higher-order collocation procedures which result in block-tridiagonal matrix systems are derived from (1) Taylor series expansions and from (2) polynomial interpolation, and the relationships between the two formulations, called respectively Hermite and spline collocation, are investigated. A Hermite block-tridiagonal system for a nonuniform mesh is derived, and the Hermite approach is extended in order to develop a variable-mesh sixth-order block-tridiagonal procedure. It is shown that all results obtained by Hermite development can be recovered by appropriate spline polynomial interpolation. The additional boundary conditions required for these higher-order procedures are also given. Comparative solutions using second-order accurate finite difference and spline and Hermite formulations are presented for the boundary layer on a flat plate, boundary layers with uniform and variable mass transfer, and the viscous incompressible Navier-Stokes equations describing flow in a driven cavity.

Rubin, S. G.↗
NASA NTRS1977

Geodetic linear estimation technique and the norm choice problem

In this work the mathematical and probabilistic background of standard linear estimation techniques used in geodesy is clarified, and their interrelationship is revealed with the help of best approximation theory and the normal equations. Emphasis is given to the separation of the deterministic solution to the approximation problem from the probabilistic justification of the metric of the approximation. Least squares prediction has been related to deterministic (exact) collocation, and minimum error bound has been identified as a prediction optimality criterion in the latter. Criteria for the optimal choice of norm in Hilbert space collocation are proposed for gravimetric geodesy problems.

Dermanis, A.↗
NASA NTRS1977

Compliance calibration of specimens used in the R-curve practice

The compliance calibrations for the compact (CS) and crack-line-wedge-loaded (CLWL) specimens have been determined by experimental measurements and by boundary-collocation analysis. The CS and CLWL specimen configurations were modeled more accurately than those used in previous analytical investigations. Polynomial expressions for the compliance at various stations along the crack line for CS and CLWL specimens are presented. The compliance calibrations for the center-crack tension (CCT) specimen have been determined theoretically by boundary-collocation and finite-element analysis. The calculated compliance values for the CCT specimen are compared with values obtained from the Irwin-Westergaard expression and from a modification to the Irwin-Westergaard expression proposed by Eftis and Liebowitz. The Eftis-Liebowitz expression was found to be in good agreement (plus or minus 2 percent) with both analyses for crack aspect ratios up to 0.8 and for gage half-span to specimen width ratios up to 0.5.

Mccabe, D. E.↗
NASA NTRS1979

Computation of aerodynamic interference effects on oscillating airfoils with controls in ventilated subsonic wind tunnels

Lift interference effects are discussed based on Bland's (1968) integral equation. A mathematical existence theory is utilized for which convergence of the numerical method has been proved for general (square-integrable) downwashes. Airloads are computed using orthogonal airfoil polynomial pairs in conjunction with a collocation method which is numerically equivalent to Galerkin's method and complex least squares. Convergence exhibits exponentially decreasing error with the number n of collocation points for smooth downwashes, whereas errors are proportional to 1/n for discontinuous downwashes. The latter can be reduced to 1/n to the m+1 power with mth-order Richardson extrapolation (by using m = 2, hundredfold error reductions were obtained with only a 13% increase of computer time). Numerical results are presented showing acoustic resonance, as well as the effect of Mach number, ventilation, height-to-chord ratio, and mode shape on wind-tunnel interference. Excellent agreement with experiment is obtained in steady flow, and good agreement is obtained for unsteady flow.

Fromme, J. A.↗
NASA NTRS1979

Intercomparison of GEOS 3 tracking systems

Measurements and time biases of the GEOS 3 tracking systems in the Atlantic calibration area were determined in relation to the Goddard lasers. In collocation tests the lasers agreed with each other to about 10 cm. Early RAMLAS laser data contained an apparent time bias of 647 microsec. Later, in a laser collocation test, RAMLAS differed from MOBLAS 2 by +11.3 cm and -28 microsec. An overlooked correction reduced the -28 microsec to -8 microsec. The geoceivers initially had range rate biases of 3-6 m/s and constant time biases of several milliseconds. During the first month of GEOS 3 operations the range rate biases decreases monotonically by 2.5 + or - 1.0 m/s. The C band radar biases relative to the lasers were usually within several meters and several tenths of milliseconds. The unified S band system relative biases were usually within + or - 4 mm/s and + or - 0.1 ms. The average and scatter of the altimeter biases recovered on 18 laser reference arcs was - 5.3 + or - 0.6 m.

Berbert, J. H.↗
NASA NTRS1979

Analysis of Combustion Instability in Liquid Fuel Rocket Motors

The development of a technique to be used in the solution of nonlinear velocity-sensitive combustion instability problems is described. The orthogonal collocation method was investigated. It found that the results are heavily dependent on the location of the collocation points and characteristics of the equations, so the method was rejected as unreliabile. The Galerkin method, which has proved to be very successful in analysis of the pressure sensitive combustion instability was found to work very well. It was found that the pressure wave forms exhibit a strong second harmonic distortion and a variety of behaviors are possible depending on the nature of the combustion process and the parametric values involved. A one-dimensional model provides further insight into the problem by allowing a comparison of Galerkin solutions with more exact finite-difference computations.

Wong, K. W.↗
NASA NTRS1983

Spectral methods for the Euler equations

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collocation methods are put in historical context. The basic concepts of both Fourier and Chebyshev spectral collocation methods are provided. Filtering strategies for both shock-fitting and shock-capturing approaches are also presented. Fourier shock capturing techniques are evaluated using a one-dimensional, periodic astrophysical 'nozzle' problem. Examples of shock-fitting approaches include a shock/acoustic wave interaction, shock/vortex interaction, and the classical blunt body problem. While the shock capturing spectral method does not yet show a clear advantage over second-order finite differences, equivalent accuracy can be obtained using shock fitting with far fewer grid points.

Hussaini, M. Y.↗
NASA NTRS1984

Spectral methods for the Euler equations: Fourier methods and shock-capturing

Spectral methods for compressible flows are introduced in relation to finite difference and finite element techniques within the framework of the method of weighted residuals. Current spectral collocation methods are put in historical context. The basic concepts of Fourier spectral collocation methods are provided. Filtering strategies for shock-capturing approaches are also presented. Fourier shock capturing techniques are evaluated using a one dimensional, periodic astrophysical ""nozzle'' problem.

Hussaini, M. Y.↗
NASA NTRS1984

Precise control of flexible manipulators

Experimental apparatus were developed for physically testing control systems for pointing flexible structures, such as limber spacecraft, for the case that control actuators cannot be collocated with sensors. Structural damping ratios are less than 0.003, each basic configuration of sensor/actuator noncollocation is available, and inertias can be halved or doubled abruptly during control maneuvers, thereby imposing, in particular, a sudden reversal in the plant's pole-zero sequence. First experimental results are presented, including stable control with both collocation and noncollocation.

Cannon, R. H., Jr.↗
NASA NTRS1984

Contact stresses in pin-loaded orthotropic plates

The effects of pin elasticity, friction, and clearance on the stresses near the hole in a pin-loaded orthotropic plate are described. The problem is modeled as a contact elasticity problem using complex variable theory, the pin and the plate being two elastic bodies interacting through contact. This modeling is in contrast to previous works which assumed that the pin is rigid or that it exerts a known cosinusoidal radial traction on the hole boundary. Neither of these approaches explicitly involves a pin. A collocation procedure and iteration were used to obtain numerical results for a variety of plate and pin elastic properties and various levels of friction and clearance. Collocation was used to enforce the boundary and iteration was used to find the contact and no-slip regions on the boundary. Details of the numerical scheme are discussed.

Hyer, M. W.↗
NASA NTRS1979

Solution of 3-dimensional time-dependent viscous flows. Part 1: Investigation of candidate algoriths

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow to aid in the design of helicopter rotors. Numerical algorithms are examined to determine their overall suitability for the efficient and routine solution of an appropriate system of partial differential equations. It is concluded that a consistently split time linearized block implicit scheme using either quintic B spline collocation or the generalized operator compact implicit approach to generate a fourth order accurate algorithm is particularly well suited for use on the present problem. High cell Reynolds number behavior leads to favoring the generalized operator compact implicit approach over the quintic B spline collocation method.

Weinberg, B. C.↗
NASA NTRS1984

Dugdale plastic zone size and CTOD equations for the compact specimen

With the aim of applying the Dugdale model to the compact specimen, the equations for stress intensity factor and crack surface displacement are obtained for a pair of pin loads and a uniform stress acting on a segment of the crack surface. The plastic zone size is calculated from the 'finiteness' condition of Dugdale, and the results agree well with collocation results from Terada (1983). The load that causes incipient yielding at the compression point of the compact specimen is calculated from a finite element strip yield analysis for an elastic-perfectly plastic material. The crack tip opening displacement for the compact specimen is calculated by adding the displacement at the tip of the physical crack length due to the pin load and due to the uniform stress. The results are within 1.5 percent of current collocation results.

Newman, J. C., Jr.↗
NASA NTRS1984

A spectral element method for fluid dynamics - Laminar flow in a channel expansion

A spectral element method that combines the generality of the finite element method with the accuracy of spectral techniques is proposed for the numerical solution of the incompressible Navier-Stokes equations. In the spectral element discretization, the computational domain is broken into a series of elements, and the velocity in each element is represented as a high-order Lagrangian interpolant through Chebyshev collocation points. The hyperbolic piece of the governing equations is then treated with an explicit collocation scheme, while the pressure and viscous contributions are treated implicitly with a projection operator derived from a variational principle. The implementation of the technique is demonstrated on a one-dimensional inflow-outflow advection-diffusion equation, and the method is then applied to laminar two-dimensional (separated) flow in a channel expansion. Comparisons are made with experiment and previous numerical work.

Patera, A. T.↗
NASA NTRS1986

A spectral multi-domain technique with application to generalized curvilinear coordinates

Spectral collocation methods have proven to be efficient discretization schemes for many aerodynamic and fluid mechanic problems. The high order accuracy and resolution shown by these methods allows one to obtain engineering accuracy solutions on coarse meshes, or alternatively, to obtain solutions with very small error. One drawback to these techniques was the requirement that a complicated physical domain must map into a simple computational domain for discretization. This mapping must be smooth if the high order accuracy and expontential convergence rates associated with spectral methods are to be preserved. Additionally even smooth stretching transformations can decrease the accuracy of a spectral method, if the stretching is severe. A further difficulty with spectral methods was in their implementation on parallel processing computers, where efficient spectral algorithms were lacking. The above restrictions are overcome by splitting the domain into regions, each of which preserve the advantages of spectral collocation, and allow the ratio of the mesh spacing between regions to be several orders of magnitude higher than allowable in a single domain. Such stretchings would be required to resolve the thin viscous region in an external aerodynamic problem. Adjoining regions are interfaced by enforcing a global flux balance which preserves high-order continuity of the solution, regardless of the type of the equations being solved.

Macaraeg, M. G.↗
NASA NTRS1986

Preconditioning for first-order spectral discretization

Efficient solution of the equations from spectral discretizations is essential if the high-order accuracy of these methods is to be realized. Direct solution of these equations is rarely feasible, thus iterative techniques are required. A preconditioning scheme for first-order Chebyshev collocation operators is proposed herein, in which the central finite difference mesh is finer than the collocation mesh. Details of the proper techniques for transferring information between the meshes are given here, and the scheme is analyzed by examination of the eigenvalue spectra of the preconditioned operators. The effect of artificial viscosity required in the inversion of the finite difference operator is examined. A second preconditioning scheme, involving a high-order upwind finite difference operator of the van Leer type is also analyzed to provide a comparison with the present scheme. Finally, the performance of the present scheme is verified by application to several test problems.

Streett, C. L.↗
NASA NTRS1985

Contact stresses in pin-loaded orthotropic plates

The effects of pin elasticity, friction, and clearance on the stresses near the hole in a pin-loaded orthotropic plate are described. The problem is modeled as a contact elasticity problem using complex variable theory, the pin and the plate being two elastic bodies interacting through contact. This modeling is in contrast to previous works which assumed that the pin is rigid or that it exerts a known cosinusoidal radial traction on the hole boundary. Neither of these approaches explicitly involves a pin. A collocation procedure and iteration were used to obtain numerical results for a variety of plate and pin elastic properties and various levels of friction and clearance. Collocation was used to enforce the boundary and iteration was used to find the contact and no-slip regions on the boundary. Details of the numerical scheme are discussed.

Hyer, M. W.↗