SEARCH · Engineering Papers

Results for “collocated”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Error analysis for spectral approximation of the Korteweg-De Vries equation

The conservation and convergence properties of spectral Fourier methods for the numerical approximation of the Korteweg-de Vries equation are analyzed. It is proved that the (aliased) collocation pseudospectral method enjoys the same convergence properties as the spectral Galerkin method, which is less effective from the computational point of view. This result provides a precise mathematical answer to a question raised by several authors in recent years.

Maday, Y.↗

NASA NTRS1987

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

NASA NTRS1987

A new method of imposing boundary conditions for hyperbolic equations

A new method to impose boundary conditions for pseudospectral approximations to hyperbolic equations is suggested. This method involves the collocation of the equation at the boundary nodes as well as satisfying boundary conditions. Stability and convergence results are proven for the Chebyshev approximation of linear scalar hyperbolic equations. The eigenvalues of this method applied to parabolic equations are shown to be real and negative.

Funaro, D.↗

NASA NTRS1986

Plates and shells containing a surface crack under general loading conditions

The severity of the underlying assumptions of the line-spring model (LSM) are such that verification with three-dimensional solutions is necessary. Such comparisons show that the model is quite accurate, and therefore, its use in extensive parameter studies is justified. Investigations into the endpoint behavior of the line-spring model have led to important conclusions about the ability of the model to predict stresses in front of the crack tip. An important application of the LSM was to solve the contact plate bending problem. Here the flexibility of the model to allow for any crack shape is exploited. The use of displacement quantities as unknowns in the formulation of the problem leads to strongly singular integral equations, rather than singular integral equations which result from using displacement derivatives. The collocation method of solving the integral equations was found to be better and more convenient than the quadrature technique. Orthogonal polynomials should be used as fitting functions when using the LSM as opposed to simpler functions such as power series.

Joseph, Paul F.↗

NASA NTRS1987

Geoid undulation computations at laser tracking stations

Geoid undulation computations were performed at 29 laser stations distributed around the world using a combination of terrestrial gravity data within a cap of radius 2 deg and a potential coefficient set up to 180 deg. The traditional methods of Stokes' and Meissl's modification together with the Molodenskii method and the modified Sjoberg method were applied. Performing numerical tests based on global error assumptions regarding the terrestrial data and the geopotential set it was concluded that the modified Sjoberg method is the most accurate and promising technique for geoid undulation computations. The numerical computations for the geoid undulations using all the four methods resulted in agreement with the ellipsoidal minus orthometric value of the undulations on the order of 60 cm or better for most of the laser stations in the eastern United States, Australia, Japan, Bermuda, and Europe. A systematic discrepancy of about 2 meters for most of the western United States stations was detected and verified by using two relatively independent data sets. For oceanic laser stations in the western Atlantic and Pacific oceans that have no terrestrial data available, the adjusted GEOS-3 and SEASAT altimeter data were used for the computation of the geoid undulation in a collocation method.

Despotakis, Vasilios K.↗

NASA NTRS1986

Preconditioning matrices for Chebyshev derivative operators

The problem of preconditioning the matrices arising from pseudo-spectral Chebyshev approximations of first order operators is considered in both one and two dimensions. In one dimension a preconditioner represented by a full matrix which leads to preconditioned eigenvalues that are real, positive, and lie between 1 and pi/2, is already available. Since there are cases in which it is not computationally convenient to work with such a preconditioner, a large number of preconditioners were studied which were more sparse (in particular three and four diagonal matrices). The eigenvalues of such preconditioned matrices are compared. The results were applied to the problem of finding the steady state solution to an equation of the type u sub t = u sub x + f, where the Chebyshev collocation is used for the spatial variable and time discretization is performed by the Richardson method. In two dimensions different preconditioners are proposed for the matrix which arises from the pseudo-spectral discretization of the steady state problem. Results are given for the CPU time and the number of iterations using a Richardson iteration method for the unpreconditioned and preconditioned cases.

Rothman, Ernest E.↗

NASA NTRS1986

Failure-accommodating control of large flexible spacecraft

This paper considers the problem of designing failure-accommodating controllers for large flexible spacecraft when there are sector-type nonlinearities in the loops. It is proved that an LQG-type controller can be made failure tolerant by inserting appropriate gains in the actuator paths and state estimator residual paths. For the state feedback case, when the actuators are saturating type, it is proved that there exists a finite region of attraction, which is invariant in the presence of these gains. Another class of controllers, which employs collocated sensors and actuators is presented, and is shown to have excellent failure-accommodation properties in addition to its robustness properties.

Joshi, S. M.↗

NASA NTRS1986

Conforming versus non-conforming boundary elements in three-dimensional elastostatics

A critical comparison of two basic formulations in three-dimensional elastostatics, using conforming and nonconforming boundary elements, is presented. The basic structure of the boundary element method is developed. The peculiarities that both types of boundary elements present in relation to the numerical implementation are discussed. Through selected examples, key issues such as the computational advantages and disadvantages of both formulations, mesh discretization and accuracy questions, and optimal location of the collocation nodes in the case of nonconforming elements are addressed. It is shown that conforming elements are able to produce more accurate results than nonconforming ones, with substantial economy in the final size of the system equations.

Manolis, G. D.↗

NASA NTRS1987

Calculation of stresses and forces between the skin and stiffener in composite panels

A methodology for computing the stresses at the interface between the skin and stiffener in stiffened composite panels is described. The methodology is based on finite-element analyses and on elasticity solution. The finite-element analyses are standard, while the elasticity solution is based on eigenvalue expansions of the stress function. The eigenvalue expansions are assumed to be valid in the local region where the stiffener flange terminates. The local elasticity solution is coupled to the global finite-element analysis using collocation on the boundary of the local region. Accuracy and convergence of the method are discussed and several examples of its utility are presented.

Hyer, M. W.↗

NASA NTRS1987

Spectral methods in fluid dynamics

Fundamental aspects of spectral methods are introduced. Recent developments in spectral methods are reviewed with an emphasis on collocation techniques. Their applications to both compressible and incompressible flows, to viscous as well as inviscid flows, and also to chemically reacting flows are surveyed. The key role that these methods play in the simulation of stability, transition, and turbulence is brought out. A perspective is provided on some of the obstacles that prohibit a wider use of these methods, and how these obstacles are being overcome.

Hussaini, M. Y.↗

NASA NTRS1986

Narrow- and broad-band satellite measurements of shortwave radiation - Conversion simulations with a general circulation model

The Oregon State University/Lawrence Livermore National Laboratory general circulation model has been employed as a vehicle for suggesting and exploring various means of converting narrow-band measurements of reflected solar radiation from the earth-atmosphere system to broad-band quantities. For purely illustrative purposes within the model's solar radiation routine, a narrow-band filter function consisting of a square-wave window extending from 0.5 to 0.9 microns is adopted. A limitation of the model, for this sort of endeavor, is that it does not include the wavelength dependence of surface albedos. Nevertheless, the model simulations tend to mimic the calibration of a narrow-band instrument, utilizing reflected solar radiation from the earth-atmosphere system as simultaneously measured by a collocated broad-band instrument; for the model, however, this is done in terms of fluxes, in contrast to instrument-measured radiances. The model results suggest that it might be preferable to perform narrow- to broad-band conversions in terms of planetary albedo (or an equivalent quantity), rather than in terms of reflected fluxes or radiances. Further improvement is achieved if, for instruments that can differentiate between clear and overcast conditions, separate clear and overcast calibrations are performed.

Cess, Robert D.↗

NASA NTRS1987

The effects of pin elasticity, clearance, and friction on the stresses in a pin-loaded orthotropic plate

The effects of pin elasticity, clearance, and friction on the stresses in a pin loaded orthotropic plate are studied. The effects are studied by posing the problem as a planar contact elasticity problem, the pin and the plate being two elastic bodies which interact through contact. Coulomb friction is assumed, the pin loads the plate in one of its principal material directions, and the plate is infinite in extent. A collocation scheme and interaction, in conjunction with a complex variable series solution, are used to obtain numerical results. The contact region between the plate and pin is unknown and must be solved for as part of the solution. The same is true of the region of friction induced no slip. Two pin stiffnesses, two clearance levels, two friction levels and two laminates, a (0/+ or - 45/90)s and a (02/+ or - 45)s, are studied. The effects of pin elasticity, clearance, and friction on the load capacity of the plate are assessed by comparing the load capacity of the plate with the capacity when the pin is rigid, perfectly fitting, and frictionless.

Hyer, M. W.↗

NASA NTRS1987

Finite length Taylor Couette flow

Axisymmetric numerical solutions of the unsteady Navier-Stokes equations for flow between concentric rotating cylinders of finite length are obtained by a spectral collocation method. These representative results pertain to two-cell/one-cell exchange process, and are compared with recent experiments.

Streett, C. L.↗

NASA NTRS1987

Viscoelastic/damage modeling of filament-wound spherical pressure vessels

A model of the viscoelastic/damage response of a filament-wound spherical vessel used for long-term pressure containment is developed. The matrix material of the composite system is assumed to be linearly viscoelastic. Internal accumulated damage based upon a quadratic relationship between transverse modulus and maximum circumferential strain is postulated. The resulting nonlinear problem is solved by an iterative routine. The elastic-viscoelastic correspondence is employed to produce, in the Laplace domain, the associated elastic solution for the maximum circumferential strain which is inverted by the method of collocation to yield the time-dependent solution. Results obtained with the model are compared to experimental observations.

Hackett, Robert M.↗

NASA NTRS1987

Stability and transition in supersonic boundary layers

The full three-dimensional time-dependent compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation algorithm to study the stability of supersonic flows over a flat plate. Several non-linear numerical experiments suggest the existence of a secondary instability which might provide a possible route to transition. The interaction of the modes involved in this secondary instability is possibly amenable to a Floquet theory. Pertinent differences between this instability and the more common incompressible K-type instabilities are pointed out.

Erlebacher, Gordon↗

NASA NTRS1987

A numerical simulation of finite-length Taylor-Couette flow

The processes leading to laminar-turbulent transition in finite-channel-length Taylor-Couette flow are investigated analytically, solving the unsteady incompressible Navier-Stokes equations by spectral-collocation methods. A time-split algorithm, implementable in both axisymmetric and fully three-dimensional time-accurate versions, and an algorithm based on the staggered-mesh discretization of Bernardi and Maday (1986) are described in detail, and results obtained by applying the axisymmetric version of the first algorithm and a steady-state version of the second are presented graphically and compared with published experimental data. The feasibility of full three-dimensional simulations of the progression through chaotic states to turbulence under the constraints of Taylor-Couette flow is demonstrated.

Streett, C. L.↗

NASA NTRS1987

Flexible system model reduction and control system design based upon actuator and sensor influence functions

A model reduction technique based on aggregation with respect to sensor and actuator influence functions rather than modes is presented for large systems of coupled second-order differential equations. Perturbation expressions which can predict the effects of spillover on both the reduced-order plant model and the neglected plant model are derived. For the special case of collocated actuators and sensors, these expressions lead to the derivation of constraints on the controller gains that are, given the validity of the perturbation technique, sufficient to guarantee the stability of the closed-loop system. A case study demonstrates the derivation of stabilizing controllers based on the present technique. The use of control and observation synthesis in modifying the dimension of the reduced-order plant model is also discussed. A numerical example is provided for illustration.

Yam, Yeung↗

NASA NTRS1987

An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions.

Bayliss, A.↗