Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary-value problems”

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 109 records · Page 6

An MHD model of the earth's magnetosphere

It is pointed out that the earth's magnetosphere arises from the interaction of the solar wind with the earth's geomagnetic field. A global magnetohydrodynamics (MHD) model of the earth's magnetosphere has drawn much attention in recent years. In this model, MHD equations are used to describe the solar wind interaction with the magnetosphere. In the present paper, some numerical aspects of the model are considered. Attention is given to the ideal MHD equations, an equation of state for the plasma, the model as an initial- and boundary-value problem, the shock capturing technique, computational requirements and techniques for global MHD modeling, a three-dimensional mesh system employed in the global MHD model, and some computational results.

Wu, C. C.↗

Collocation for an integral equation arising in duct acoustics

A mathematical model is developed to describe the effect of aircraft-engine inlet geometry on the reflected and radiated acoustic field without flow, as studied experimentally using a spinning-mode synthesizer by Silcox (1983). The acoustic pressure in the inlet interior and exterior is modeled by a pure cylindrical azimuthal mode for the Helmholtz equation with hardwall boundary and by the Helmholtz equation and the radiation condition at infinity, respectively. The analytical approach to the solution of the resulting boundary-value problem and the program implementation are explained; numerical results are presented in tables and graphs; and the uniqueness of the problem is demonstrated.

Moss, W. F.↗

Reducing the number of variational equations in the implementation of multiple shooting

The standard method of multiple shooting for a system of n first-order differential equations with k unknown initial conditions requires the integration of k sets of variational equations on the first shot and n sets of variational equations on every shot thereafter. This paper describes a variant of multiple shooting that requires the solution of k sets of variational equations on every shot. The technique applies to both linear and nonlinear boundary-value problems. Techniques to deal with difficulties unique to the solution of nonlinear problems are suggested.

Krogh, F. T.↗

Theoretical Efficiencies of Microwave Diode Triplers

Report discusses computer simulations of 300- to 900-GHz triplers using nonideal GaAs Schottky diodes operating in varistor mode. Nonlinear boundary-value problem posed by diode state equations not solved in closed form. Consequently, such computer simulations needed to optimize tripler configuration and operating conditions.

Frerking, M. A.↗

Optimal torque control for SCOLE slewing maneuvers

The Spacecraft Control Laboratory Experiment (SCOLE) was slewed from one attitude to the required attitude and an integral performance index which involves the control torques was minimized. Kinematic and dynamical equations, optimal control, two-point boundary-value problems, and estimation of unknown boundary conditions are presented.

Bainum, P. M.↗

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↗

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.↗

Minimum-time pointing control of a two-link manipulator

Minimum-time pointing control for the end-effector of a planar, two-link manipulator is developed. Minimum-time pointing control is a new area of research for multilink manipulators, which can be applied to rapid retargeting control of a multibody spacecraft. The minimum-time control problem of aligning the second link of the two-link manipulator with a given target point is considered. A numerical method called the minimizing-boundary-condition method is used to determine optimal solutions for the two-point boundary-value problem associated with first-order necessary conditions. Minimum-time solutions for different models of pointing systems are compared. The results of the comparison show that a two-link manipulator with two degree-of-freedom performs better than a conventional one-link system for minimum-time pointing.

Wie, Bong↗

Boundary force method analyses of notched plates with cracks

The boundary-force method developed by Tan (1985) and Tan et al. (1986) is extended to two-dimensional mixed boundary-value problems in fracture mechanics. The formulation is explained in detail, and results are presented in tables and graphs for sample problems involving an edge crack emanating from a semielliptical and a semicircular notch and an edge crack emanating from a V notch. The accuracy of the approach is demonstrated by comparisons with collocation solutions for an edge-clamped center-cracked tension specimen.

Newman, J. C., Jr.↗

Methods for analysis of cracks in three-dimensional solids

Various methods used for determining stress-intensity factors for cracked three-dimensional bodies are reviewed. The review is limited to the determination of mode-I stress-intensity factors. Some exact solutions for cracks in infinite solids are presented, and various approximate methods that have been used for the solution of the boundary-value problem of finite solids with cracks are addressed. The techniques used to extract the stress-intensity factors from these solutions are considered, and various methods used to calculate the stress-intensity factors for through-the-thickness cracks, semielliptical surface cracks, and quarter-elliptical corner cracks at holes are compared for the case of remote tensile loading.

Raju, I. S.↗

Use of the quasilinearization algorithm for the simulation of LSS slewing

The use of the Maximum Principle for the large angle slewing of large space structures (LSS) usually results in the so-called two-point boundary-value problem, in which many requirements (e.g., minimum time, small amplitude, and limited control power, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical algorithm. There are many candidate algorithms available for this problem (e.g., quasilinearization, gradient, etc.). Here researchers discuss only the quasilinearization method which has been used for several cases of large angle slewing of LSS. The basic idea of this algorithm is to make a series of successive approximations of the solution from a particular solvable case (linear or nonlinear) to a more general practical case. For the rigid spacecraft slewing problem with no constraints on the controls, the solution procedure can be found in the literature. This procedure needs to be modified if a minimum time for the slewing problem is desired with control limits given. Recently, an indirect method for finding the minimum time was developed to meet all these requirements. For the general mixed (including both rigid and flexible parts) problem, an additional constraint of small vibrational amplitude on the flexible parts is imposed.

Li, Feiyue↗

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↗

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↗

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗

Solving time-dependent two-dimensional eddy current problems

Transient eddy current calculations are presented for an EM wave-scattering and field-penetrating case in which a two-dimensional transverse magnetic field is incident on a good (i.e., not perfect) and infinitely long conductor. The problem thus posed is of initial boundary-value interface type, where the boundary of the conductor constitutes the interface. A potential function is used for time-domain modeling of the situation, and finite difference-time domain techniques are used to march the potential function explicitly in time. Attention is given to the case of LF radiation conditions.

Lee, Min Eig↗

Singular trajectories for time-optimal half-loop maneuvers of a high alpha fighter aircraft

Consideration is given to the problem of deriving a time-optimal open-loop control for the half-loop maneuver of a high-alpha aircraft, with initial conditions Mach 0.6 and 15,000 feet. Pontriagin's maximum principle is used to derive candidate optimal solutions. Using the two-point boundary-value algorithm, the flight path angle is maximized for various increasing specified final times until a final time of 13.6 sec yields a 180-deg flight-path angle. As the final time increased from 0.0 to 13.6 sec, the optimization process revealed 13 distinct switching structures of the control law, of which 11 contained singular arcs, and two had double singular arcs.

Hoffman, Eric↗

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗