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.

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

Vectorizable algorithms for adaptive schemes for rapid analysis of SSME flows

An initial study into vectorizable algorithms for use in adaptive schemes for various types of boundary value problems is described. The focus is on two key aspects of adaptive computational methods which are crucial in the use of such methods (for complex flow simulations such as those in the Space Shuttle Main Engine): the adaptive scheme itself and the applicability of element-by-element matrix computations in a vectorizable format for rapid calculations in adaptive mesh procedures.

Oden, J. Tinsley↗

NASA NTRS1987

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗

NASA NTRS1986

An Investigation of Energy Transmission Due to Flexural Wave Propagation in Lightweight, Built-Up Structures

A technique to measure flexural structure-borne noise intensity is investigated. Two accelerometers serve as transducers in this cross-spectral technique. The structure-borne sound power is obtained by two different techniques and compared. In the first method, a contour integral of intensity is performed from the values provided by the two-accelerometer intensity technique. In the second method, input power is calculated directly from the output of force and acceleration transducers. A plate and two beams were the subjects of the sound power comparisons. Excitation for the structures was either band-limited white noise or a deterministic signal similar to a swept sine. The two-accelerometer method was found to be sharply limited by near field and transducer spacing limitations. In addition, for the lightweight structures investigated, it was found that the probe inertia can have a significant influence on the power input to the structure. In addition to the experimental investigation of structure-borne sound energy, an extensive study of the point harmonically forced, point-damped beam boundary value problem was performed to gain insight into measurements of this nature. The intensity formulations were also incorporated into the finite element method. Intensity mappings were obtained analytically via finite element modeling of simple structures.

Mickol, John Douglas↗

NASA NTRS1987

Investigation of plasma contactors for use with orbiting wires

The proposed Shuttle-based short tether experiments with hollow cathodes have the potential for providing important data that will not be obtained in long tether experiments. A critical property for hollow cathode effectiveness as a plasma contactor is the cross magnetic field conductivity of the emitted plasma. The different effects of hollow cathode cloud overlap in the cases of motion-driven and battery-driven operation are emphasized. The calculations presented on the size and shape of the hollow cathode cloud improve the qualitative picture of hollow cathodes in low Earth orbit and provide estimates of time constants for establishing the fully-expanded cloud. The magnetic boundary value problem calculations indicate the way in which the magnetic field will effect the shape of the cloud by resisting expansion in the direction perpendicular to the field. The large-scale interactions of the system were also considered. It was concluded that recent plasma chamber experiments by Stenzel and Urrutia do not model an electrodynamic tether well enough to apply the results to tethered system behavior. Orbiting short tether experiments on hollow cathodes will provide critical information on hollow cathode performance and the underlying physics that cannot be obtained any other way. Experiments should be conducted as soon as funding and a suitable space vehicle are available.

Estes, Robert D.↗

NASA NTRS1987

Galerkin/Runge-Kutta discretizations for parabolic equations with time dependent coefficients

A new class of fully discrete Galerkin/Runge-Kutta methods is constructed and analyzed for linear parabolic initial boundary value problems with time dependent coefficients. Unlike any classical counterpart, this class offers arbitrarily high order convergence while significantly avoiding what has been called order reduction. In support of this claim, error estimates are proved, and computational results are presented. Additionally, since the time stepping equations involve coefficient matrices changing at each time step, a preconditioned iterative technique is used to solve the linear systems only approximately. Nevertheless, the resulting algorithm is shown to preserve the original convergence rate while using only the order of work required by the base scheme applied to a linear parabolic problem with time independent coefficients. Furthermore, it is noted that special Runge-Kutta methods allow computations to be performed in parallel so that the final execution time can be reduced to that of a low order method.

Keeling, Stephen L.↗

NASA NTRS1985

Notes and comments on computational elastoplasticity - Some new models and their numerical implementation

The following topics are discussed in this paper: (1) the basic interactive nature of classical elasto-plasticity and a redefinition of elastic and plastic processes that facilitates numerical calculations, (2) generalized mid-point or end-point algorithms to determine the stress increment in an elastic-plastic solid from a given strain increment, (3) an endochronic (internal time) rate theory of time-independent plasticity which encompasses various multiple-yield-surface theories and nonlinear kinematic hardening theories as its specializations, and (4) comments on finite element and boundary element methods for solving boundary value problems in elasto-plasticity.

Atluri, S. N.↗

NASA NTRS1985

Analytical results for post-buckling behaviour of plates in compression and in shear

The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.

Stein, M.↗

NASA NTRS1986

A pseudospectral Legendre method for hyperbolic equations with an improved stability condition

A new pseudospectral method is introduced for solving hyperbolic partial differential equations. This method uses different grid points than previously used pseudospectral methods: in fact the grid points are related to the zeroes of the Legendre polynomials. The main advantage of this method is that the allowable time step is proportional to the inverse of the number of grid points 1/N rather than to 1/n(2) (as in the case of other pseudospectral methods applied to mixed initial boundary value problems). A highly accurate time discretization suitable for these spectral methods is discussed.

Tal-Ezer, Hillel↗

NASA NTRS1986

Stress patterns in an interplate shear zone - An effective anisotropic model and implications for the Transverse Ranges, California

The pervasive Quaternary strike-slip (or wrench) structure that occurs in southern California is analyzed by solving a series of planar, linear elastic boundary value problems of tractional type that assume pure shear forces at infinity. It is concluded that the cross-cutting tectonic fabric of the Transverse Ranges in southern California causes the Ranges to act as a stress concentrator. In addition, it is found that rigid rotation of the Transverse Ranges is mechanically supported by an asymmetric distribution of stress and elastic strain encircling it. It is noted that for Quaternary tectonic orientation, a clockwise sense to this rotation is expected.

Ivins, E. R.↗

NASA NTRS1986

MHD simulation of the 'geoeffectiveness' of interplanetary disturbances

A series of examples is used to evaluate the 'geoeffectiveness' prediction technique that is based on the classical, initial boundary value problem of MHD. For the nonplanar simulation of a corotating stream, it is shown that the technique can compute relevant solar wind parameters at the earth's location. For the nonplanar simulation of a complex series of events, qualitative agreement is found for large-scale structures, but the phasing and amplitudes are not satisfactory. Preliminary work with a fully three-dimensional, time-dependent simulation of a flare-generated interplanetary shock wave shows that the effect of the latitudinal variation produces attenuation of the heliolongitudinal component of the IMF and of the radial velocity.

Dryer, M.↗

NASA NTRS1986

Vortical Flows Research Program of the Fluid Dynamics Research Branch

The research interests of the staff of the Fluid Dynamics Research Branch in the general area of vortex flows are summarized. A major factor in the development of enchanced maneuverability and reduced drag by aerodynamic means is the use of effective vortex control devices. The key to control is the use of emerging computational tools for predicting viscous fluid flow in close coordination with fundamental experiments. In fact, the extremely complex flow fields resulting from numerical solutions to boundary value problems based on the Navier-Stokes equations requires an intimate relationship between computation and experiment. The field of vortex flows is important in so many practical areas that a concerted effort in this area is justified. A brief background of the research activity undertaken is presented, including a proposed classification of the research areas. The classification makes a distinction between issues related to vortex formation and structure, and work on vortex interactions and evolution. Examples of current research results are provided, along with references where available. Based upon the current status of research and planning, speculation on future research directions of the group is also given.

Source record↗

NASA NTRS1988

Modeling and estimation of a low degree geopotential model from terrestrial gravity data

The development of appropriate modeling and adjustment procedures for the estimation of harmonic coefficients of the geopotential, from surface gravity data was studied, in order to provide an optimum way of utilizing the terrestrial gravity information in combination solutions currently developed at NASA/Goddard Space Flight Center, for use in the TOPEX/POSEIDON mission. The mathematical modeling was based on the fundamental boundary condition of the linearized Molodensky boundary value problem. Atmospheric and ellipsoidal corrections were applied to the surface anomalies. Terrestrial gravity solutions were found to be in good agreement with the satellite ones over areas which are well surveyed (gravimetrically), such as North America or Australia. However, systematic differences between the terrestrial only models and GEMT1, over extended regions in Africa, the Soviet Union, and China were found. In Africa, gravity anomaly differences on the order of 20 mgals and undulation differences on the order of 15 meters, over regions extending 2000 km in diameter, occur. Comparisons of the GEMT1 implied undulations with 32 well distributed Doppler derived undulations gave an RMS difference of 2.6 m, while corresponding comparison with undulations implied by the terrestrial solution gave RMS difference on the order of 15 m, which implies that the terrestrial data in that region are substantially in error.

Pavlis, Nikolaos K.↗

NASA NTRS1988

Numerical methods for incompressible viscous flows with engineering applications

A numerical scheme has been developed to solve the incompressible, 3-D Navier-Stokes equations using velocity-vorticity variables. This report summarizes the development of the numerical approximation schemes for the divergence and curl of the velocity vector fields and the development of compact schemes for handling boundary and initial boundary value problems.

Rose, M. E.↗

NASA NTRS1987

On the solution of integral equations with strongly singular kernels

Some useful formulas are developed to evaluate integrals having a singularity of the form (t-x) sup-m, m greater than or equal 1. Interpreting the integrals with strong singularities in Hadamard sense, the results are used to obtain approximate solutions of singular integral equations. A mixed boundary value problem from the theory of elasticity is considered as an example. Particularly for integral equations where the kernel contains, in addition to the dominant term (t-x) sup-m, terms which become unbounded at the end points, the present technique appears to be extremely effective to obtain rapidly converging numerical results.

Kaya, A. C.↗

NASA NTRS1988

Minimum time attitude slewing maneuvers of a rigid spacecraft

The minimum time attitude slewing motion of a rigid spacecraft with its controls provided by torques and forces, which have their upper and lower limits prescribed, is considered. The two-point boundary-value problem is derived by applying the Pontriagin's Maximum Principle to the system and solved by using a quasi-linearization algorithm. The nominal solutions to the problem as well as the starting values of the total slewing time and the unknown initial costates for this algorithm are generated by using Euler's eigenaxis rotation theorem. It is pointed out that one of the four initial costates associated with the quaternions can be arbitrarily selected without affecting the optimal controls and, thus, simplifying the computation. The minimum slewing time is determined by shortening the total slewing time until at least one of the controls becomes a bang-bang type. Several numerical tests for the rigidized SCOLE model are presented to show the applications of the methods.

Li, Feiyue↗

NASA NTRS1987

The role of damage-softened material behavior in the fracture of composites and adhesives

The failure mechanism of polymer composites and adhesives with high strain in the zone near and ahead of the crack tip is investigated analytically, summarizing the findings of Ungsuwarungsri (1986). A double-cantilever-beam specimen with the nonlinear material confined to a thin strip between two wedge-loaded elastic beams is modeled as a beam on a nonlinear foundation; the two-point boundary-value problem for stationary and propagating cracks is solved numerically; and an FEM approach is applied to study the behavior of the nonlinear strip in detail. Typical results are presented in extensive graphs, and the accuracy and efficiency of the present method are shown to be superior to those of the Berry (1963) procedure.

Ungsuwarungsri, T.↗

NASA NTRS1987

A parallelized elliptic solver for reacting flows

A modified Newton algorithm for the solution of nonlinear elliptic boundary value problems via finite discretization methods is presented. A serial implementation of this algorithm which has recently been applied successfully to the computation of an axisymmetric over-ventilated subsonic laminar methane-air jet diffusion flame is described. Parallel implementation issues and a complexity theory are presented. Included as well are actual performance data for model systems obtained on the Intel Hypercube and a discussion of its implications for modeling realistic systems.

Keyes, David E.↗

NASA NTRS1988

The dynamics and control of the in-orbit SCOLE configuration

The study of the dynamics of the Spacecraft Control Laboratory Experiment (SCOLE) is extended to emphasize the synthesis of control laws for both the linearized system as well as the large amplitude slewing maneuvers required to rapidly reorient the antenna line of sight. For control of the system through small amplitude displacements from the nominal equilibrium position LQR techniques are used to develop the control laws. Pontryagin's maximum principle is applied to minimize the time required for the slewing of a general rigid spacecraft system. The minimum slewing time is calculated based on a quasi-linearization algorithm for the resulting two point boundary value problem. The effect of delay in the control input on the stability of a continuously acting controller (designed without considering the delay) is studied analytically for a second order plant. System instability can result even for delays which are only a small fraction of the natural period of motion.

Bainum, Peter M.↗