Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary value problem”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The prediction of noise radiation from supersonic elliptic jets

This paper describes the prediction of noise radiation from supersonic elliptic jets. The noise is associated with the large scale structures in the jet mixing layer. These structures are described as instability waves. The local characteristics of the instability waves are determined from a compressible, linear, analysis. The jet mean velocity and density are described in elliptic cylindrical coordinates. The local eigensolution for the instability waves is determined from a finite difference solution of the non-separable boundary value problem. This inner solution which is formulated in terms of the method of multiple scales is matched with the radiated field using the method of matched asymptotic expansions. The form of the far-field directivity is derived. Predictions are presented for the noise radiation by the several modes of instability in the elliptic jet. The radiated field is not axisymmetric and certain modes radiate strongly in the directions of the major and minor axes of the jet. The extension of the present work to other geometries and flow fields is discussed.

Morris, Philip J.↗
NASA NTRS1993

An analytical/numerical correlation study of the multiple concentric cylinder model for the thermoplastic response of metal matrix composites

The utility of a recently developed analytical micromechanics model for the response of metal matrix composites under thermal loading is illustrated by comparison with the results generated using the finite-element approach. The model is based on the concentric cylinder assemblage consisting of an arbitrary number of elastic or elastoplastic sublayers with isotropic or orthotropic, temperature-dependent properties. The elastoplastic boundary-value problem of an arbitrarily layered concentric cylinder is solved using the local/global stiffness matrix formulation (originally developed for elastic layered media) and Mendelson's iterative technique of successive elastic solutions. These features of the model facilitate efficient investigation of the effects of various microstructural details, such as functionally graded architectures of interfacial layers, on the evolution of residual stresses during cool down. The available closed-form expressions for the field variables can readily be incorporated into an optimization algorithm in order to efficiently identify optimal configurations of graded interfaces for given applications. Comparison of residual stress distributions after cool down generated using finite-element analysis and the present micromechanics model for four composite systems with substantially different temperature-dependent elastic, plastic, and thermal properties illustrates the efficacy of the developed analytical scheme.

Pindera, Marek-Jerzy↗
NASA NTRS1993

A modified Chebyshev pseudospectral method with an O(N exp -1) time step restriction

The extreme eigenvalues of the Chebyshev pseudospectral differentiation operator are O(N exp 2) where N is the number of grid points. As a result of this, the allowable time step in an explicit time marching algorithm is O(N exp -2) which, in many cases, is much below the time step dictated by the physics of the partial differential equation. A new set of interpolating points is introduced such that the eigenvalues of the differentiation operator are O(N) and the allowable time step is O(N exp -1). The properties of the new algorithm are similar to those of the Fourier method. The new algorithm also provides a highly accurate solution for non-periodic boundary value problems.

Kosloff, Dan↗
NASA NTRS1991

A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading

An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.

Norwood, D. S.↗
NASA NTRS1950

Flutter of a uniform wing with an arbitrarily placed mass according to a differential-equation analysis and a comparison with experiment

A method is presented for the calculation of the flutter speed of a uniform wing carrying an arbitrarily placed concentrated mass. The method, an extension of recently published work by Goland and Luke, involves the solution of the differential equations of motion of the wing at flutter speed and therefore does not require the assumption of specific normal modes of vibration. The order of the flutter determinant to be solved by this method depends upon the order of the system of differential equations and not upon the number of modes of vibration involved. The differential equations are solved by operational methods, and a brief discussion of operational methods as applied to boundary-value problems is included in one of two appendixes. A comparison is made with experiment for a wing with a large eccentrically mounted weight and good agreement is obtained. Sample calculations are presented to illustrate the method; and curves of amplitudes of displacement, torque, and shear for a particular case are compared with corresponding curves computed from the first uncoupled normal modes.

Runyan, Harry L↗
NASA NTRS1954

Theoretical study of the transonic lift of a double-wedge profile with detached bow wave

A theoretical study is described of the aerodynamic characteristics at small angle of attack of a thin, double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis is carried out within the framework of the transonic (nonlinear) small-disturbance theory, and the effects of angle of attack are regarded as a small perturbation on the flow previously calculated at zero angle. The mixed flow about the front half of the profile is calculated by relaxation solution of a suitably defined boundary-value problem for transonic small-disturbance equation in the hodograph plane (i.e., the Tricomi equation). The purely supersonic flow about the rear half is found by an extension of the usual numerical method of characteristics. Analytical results are also obtained, within the framework of the same theory, for the range of speed in which the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗
NASA NTRS1989

Mechanical properties of lunar regolith and lunar soil simulant

Through the Surveyor 3 and 7, and Apollo 11-17 missions a knowledge of the mechanical properties of Lunar regolith were gained. These properties, including material cohesion, friction, in-situ density, grain-size distribution and shape, and porosity, were determined by indirect means of trenching, penetration, and vane shear testing. Several of these properties were shown to be significantly different from those of terrestrial soils, such as an interlocking cohesion and tensile strength formed in the absence of moisture and particle cementation. To characterize the strength and deformation properties of Lunar regolith experiments have been conducted on a lunar soil simulant at various initial densities, fabric arrangements, and composition. These experiments included conventional triaxial compression and extension, direct tension, and combined tension-shear. Experiments have been conducted at low levels of effective confining stress. External conditions such as membrane induced confining stresses, end platten friction and material self weight have been shown to have a dramatic effect on the strength properties at low levels of confining stress. The solution has been to treat these external conditions and the specimen as a full-fledged boundary value problem rather than the idealized elemental cube of mechanics. Centrifuge modeling allows for the study of Lunar soil-structure interaction problems. In recent years centrifuge modeling has become an important tool for modeling processes that are dominated by gravity and for verifying analysis procedures and studying deformation and failure modes. Centrifuge modeling is well established for terrestrial enginering and applies equally as well to Lunar engineering. A brief review of the experiments is presented in graphic and outline form.

Steven W Perkins↗
NASA NTRS1989

Analytical modeling of structure-soil systems for lunar bases

The study of the behavior of granular materials in a reduced gravity environment and under low effective stresses became a subject of great interest in the mid 1960's when NASA's Surveyor missions to the Moon began the first extraterrestrial investigation and it was found that Lunar soils exhibited properties quite unlike those on Earth. This subject gained interest during the years of the Apollo missions and more recently due to NASA's plans for future exploration and colonization of Moon and Mars. It has since been clear that a good understanding of the mechanical properties of granular materials under reduced gravity and at low effective stress levels is of paramount importance for the design and construction of surface and buried structures on these bodies. In order to achieve such an understanding it is desirable to develop a set of constitutive equations that describes the response of such materials as they are subjected to tractions and displacements. This presentation examines issues associated with conducting experiments on highly nonlinear granular materials under high and low effective stresses. The friction and dilatancy properties which affect the behavior of granular soils with low cohesion values are assessed. In order to simulate the highly nonlinear strength and stress-strain behavior of soils at low as well as high effective stresses, a versatile isotropic, pressure sensitive, third stress invariant dependent, cone-cap elasto-plastic constitutive model was proposed. The integration of the constitutive relations is performed via a fully implicit Backward Euler technique known as the Closest Point Projection Method. The model was implemented into a finite element code in order to study nonlinear boundary value problems associated with homogeneous as well as nonhomogeneous deformations at low as well as high effective stresses. The effect of gravity (self-weight) on the stress-strain-strength response of these materials is evaluated. The calibration of the model is performed via three techniques: (1) physical identification, (2) optimized calibration at the constitutive level, and (3) optimized calibration at the finite element level (Inverse Identification). Activities are summarized in graphic and outline form.

Macari-Pasqualino, Jose Emir↗
NASA NTRS1993

A hybrid Pade-Galerkin technique for differential equations

A three-step hybrid analysis technique, which successively uses the regular perturbation expansion method, the Pade expansion method, and then a Galerkin approximation, is presented and applied to some model boundary value problems. In the first step of the method, the regular perturbation method is used to construct an approximation to the solution in the form of a finite power series in a small parameter epsilon associated with the problem. In the second step of the method, the series approximation obtained in step one is used to construct a Pade approximation in the form of a rational function in the parameter epsilon. In the third step, the various powers of epsilon which appear in the Pade approximation are replaced by new (unknown) parameters (delta(sub j)). These new parameters are determined by requiring that the residual formed by substituting the new approximation into the governing differential equation is orthogonal to each of the perturbation coordinate functions used in step one. The technique is applied to model problems involving ordinary or partial differential equations. In general, the technique appears to provide good approximations to the solution even when the perturbation and Pade approximations fail to do so. The method is discussed and topics for future investigations are indicated.

Geer, James F.↗
NASA NTRS1993

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗
NASA NTRS1993

Rational trigonometric approximations using Fourier series partial sums

A class of approximations (S(sub N,M)) to a periodic function f which uses the ideas of Pade, or rational function, approximations based on the Fourier series representation of f, rather than on the Taylor series representation of f, is introduced and studied. Each approximation S(sub N,M) is the quotient of a trigonometric polynomial of degree N and a trigonometric polynomial of degree M. The coefficients in these polynomials are determined by requiring that an appropriate number of the Fourier coefficients of S(sub N,M) agree with those of f. Explicit expressions are derived for these coefficients in terms of the Fourier coefficients of f. It is proven that these 'Fourier-Pade' approximations converge point-wise to (f(x(exp +))+f(x(exp -)))/2 more rapidly (in some cases by a factor of 1/k(exp 2M)) than the Fourier series partial sums on which they are based. The approximations are illustrated by several examples and an application to the solution of an initial, boundary value problem for the simple heat equation is presented.

Geer, James F.↗
NASA NTRS1993

A nonconforming multigrid method using conforming subspaces

For second-order elliptic boundary value problems, we develop a nonconforming multigrid method using the coarser-grid correction on the conforming finite element subspaces. The convergence proof with an arbitrary number of smoothing steps for nu-cycle is presented.

Lee, Chang Ock↗
NASA NTRS1994

Energy and maximum norm estimates for nonlinear conservation laws

We have devised a technique that makes it possible to obtain energy estimates for initial-boundary value problems for nonlinear conservation laws. The two major tools to achieve the energy estimates are a certain splitting of the flux vector derivative f(u)(sub x), and a structural hypothesis, referred to as a cone condition, on the flux vector f(u). These hypotheses are fulfilled for many equations that occur in practice, such as the Euler equations of gas dynamics. It should be noted that the energy estimates are obtained without any assumptions on the gradient of the solution u. The results extend to weak solutions that are obtained as point wise limits of vanishing viscosity solutions. As a byproduct we obtain explicit expressions for the entropy function and the entropy flux of symmetrizable systems of conservation laws. Under certain circumstances the proposed technique can be applied repeatedly so as to yield estimates in the maximum norm.

Olsson, Pelle↗
NASA NTRS1994

Fourth-order 2N-storage Runge-Kutta schemes

A family of five-stage fourth-order Runge-Kutta schemes is derived; these schemes required only two storage locations. A particular scheme is identified that has desirable efficiency characteristics for hyperbolic and parabolic initial (boundary) value problems. This scheme is competitive with the classical fourth-order method (high-storage) and is considerably more efficient and accurate than existing third-order low-storage schemes.

Carpenter, Mark H.↗
NASA NTRS1994

Parabolized stability equations

The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Herbert, Thorwald↗
NASA NTRS1994

Optimization of Residual Stresses in MMC's Using Compensating/Compliant Interfacial Layers. Part 2: OPTCOMP User's Guide

A user's guide for the computer program OPTCOMP is presented in this report. This program provides a capability to optimize the fabrication or service-induced residual stresses in uni-directional metal matrix composites subjected to combined thermo-mechanical axisymmetric loading using compensating or compliant layers at the fiber/matrix interface. The user specifies the architecture and the initial material parameters of the interfacial region, which can be either elastic or elastoplastic, and defines the design variables, together with the objective function, the associated constraints and the loading history through a user-friendly data input interface. The optimization procedure is based on an efficient solution methodology for the elastoplastic response of an arbitrarily layered multiple concentric cylinder model that is coupled to the commercial optimization package DOT. The solution methodology for the arbitrarily layered cylinder is based on the local-global stiffness matrix formulation and Mendelson's iterative technique of successive elastic solutions developed for elastoplastic boundary-value problems. The optimization algorithm employed in DOT is based on the method of feasible directions.

Pindera, Marek-Jerzy↗
NASA NTRS1994

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗
NASA NTRS1994

An Efficient Spectral Method for Ordinary Differential Equations with Rational Function Coefficients

We present some relations that allow the efficient approximate inversion of linear differential operators with rational function coefficients. We employ expansions in terms of a large class of orthogonal polynomial families, including all the classical orthogonal polynomials. These families obey a simple three-term recurrence relation for differentiation, which implies that on an appropriately restricted domain the differentiation operator has a unique banded inverse. The inverse is an integration operator for the family, and it is simply the tridiagonal coefficient matrix for the recurrence. Since in these families convolution operators (i.e. matrix representations of multiplication by a function) are banded for polynomials, we are able to obtain a banded representation for linear differential operators with rational coefficients. This leads to a method of solution of initial or boundary value problems that, besides having an operation count that scales linearly with the order of truncation N, is computationally well conditioned. Among the applications considered is the use of rational maps for the resolution of sharp interior layers.

Coutsias, Evangelos A.↗