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 397 records · Page 22

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Computation of multi-dimensional viscous supersonic jet flow

A new method has been developed for two- and three-dimensional computations of viscous supersonic flows with embedded subsonic regions adjacent to solid boundaries. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases relevant to internal supersonic flow is presented and compared with data. Comparison between data and computation are in general excellent thus indicating that the computational technique has great promise as a tool for calculating supersonic flow with embedded subsonic regions. Finally, a User's Manual is presented for the computer code used to perform the calculations.

Kim, Y. N.↗

NASA NTRS1986

Computation of multi-dimensional viscous supersonic flow

A method has been developed for two- and three-dimensional computations of viscous supersonic jet flows interacting with an external flow. The approach employs a reduced form of the Navier-Stokes equations which allows solution as an initial-boundary value problem in space, using an efficient noniterative forward marching algorithm. Numerical instability associated with forward marching algorithms for flows with embedded subsonic regions is avoided by approximation of the reduced form of the Navier-Stokes equations in the subsonic regions of the boundary layers. Supersonic and subsonic portions of the flow field are simultaneously calculated by a consistently split linearized block implicit computational algorithm. The results of computations for a series of test cases associated with supersonic jet flow is presented and compared with other calculations for axisymmetric cases. Demonstration calculations indicate that the computational technique has great promise as a tool for calculating a wide range of supersonic flow problems including jet flow. Finally, a User's Manual is presented for the computer code used to perform the calculations.

Buggeln, R. C.↗

NASA NTRS1996

New Nonlinear Multigrid Analysis

The nonlinear multigrid is an efficient algorithm for solving the system of nonlinear equations arising from the numerical discretization of nonlinear elliptic boundary problems. In this paper, we present a new nonlinear multigrid analysis as an extension of the linear multigrid theory presented by Bramble. In particular, we prove the convergence of the nonlinear V-cycle method for a class of mildly nonlinear second order elliptic boundary value problems which do not have full elliptic regularity.

Xie, Dexuan↗

NASA NTRS1987

Numerical solution of the radiation transport equation in disk geometry

An efficient numerical method for solving the problem of radiation transport in a dusty medium with two dimensional (2-D) disk geometry is described. It is a generalization of the one-dimensional quasi-diffusion method in which the transport equation is cast in diffusion form and then solved as a boundary value problem. The method should be applicable to a variety of astronomical sources, the dynamics of which are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, interstellar molecular clouds, accretion disks, and disk galaxies. The computational procedure and practical considerations for implementing the method are described in detail. To illustrate the effects of 2-D radiation transport, some model results (dust temperature distributions and IR flux spectra) for externally heated, interstellar dust clouds with spherically symmetric and disk geometry are compared.

Spagna, George F., Jr.↗

NASA NTRS1989

Thermo-viscoelastic analysis of composite materials

The thermo-viscoelastic boundary value problem for anisotropic materials is formulated and a numerical procedure is developed for the efficient analysis of stress and deformation histories in composites. The procedure is based on the finite element method and therefore it is applicable to composite laminates containing geometric discontinuities and complicated boundary conditions. Using the present formulation, the time-dependent stress and strain distributions in both notched and unnotched graphite/epoxy composites have been obtained. The effect of temperature and ply orientation on the creep and relaxation response is also studied.

Lin, Kuen Y.↗

NASA NTRS1985

Optimal landing of a helicopter in autorotation

Gliding descent in autorotation is a maneuver used by helicopter pilots in case of engine failure. The landing of a helicopter in autorotation is formulated as a nonlinear optimal control problem. The OH-58A helicopter was used. Helicopter vertical and horizontal velocities, vertical and horizontal displacement, and the rotor angle speed were modeled. An empirical approximation for the induced veloctiy in the vortex-ring state were provided. The cost function of the optimal control problem is a weighted sum of the squared horizontal and vertical components of the helicopter velocity at touchdown. Optimal trajectories are calculated for entry conditions well within the horizontal-vertical restriction curve, with the helicopter initially in hover or forwared flight. The resultant two-point boundary value problem with path equality constraints was successfully solved using the Sequential Gradient Restoration Technique.

Lee, A. Y. N.↗

NASA NTRS1991

Air-breathing hypersonic vehicle guidance and control studies; An integrated trajectory/control analysis methodology: Phase 1

A tool which generates optimal trajectory/control histories in an integrated manner is generically adapted to the treatment of single-stage-to-orbit air-breathing hypersonic vehicles. The methodology is implemented as a two point boundary value problem solution technique. Its use permits an assessment of an entire near-minimum-fuel trajectory and desired control strategy from takeoff to orbit while satisfying physically derived inequality constraints and while achieving efficient propulsive mode phasing. A simpler analysis strategy that partitions the trajectory into several boundary condition matched segments is also included to construct preliminary trajectory and control history representations with less computational burden than is required for the overall flight profile assessment. A demonstration was accomplished using a tabulated example (winged-cone accelerator) vehicle model that is combined with a newly developed multidimensional cubic spline data smoothing routine. A constrained near-fuel-optimal trajectory, imposing a dynamic pressure limit of 1000 psf, was developed from horizontal takeoff to 20,000 ft/sec relative air speed while aiming for a polar orbit. Previously unspecified propulsive discontinuities were located. Flight regimes demanding rapid attitude changes were identified, dictating control effector and closed-loop controller authority was ascertained after evaluating effector use for vehicle trim. Also, inadequacies in vehicle model representations and specific subsystem models with insufficient fidelity were determined based on unusual control characteristics and/or excessive sensitivity to uncertainty.

Hattis, Philip D.↗

NASA NTRS1974

Low thrust space vehicle trajectory optimization using regularized variables

Optimizing the trajectory of a low thrust space vehicle usually means solving a nonlinear two point boundary value problem. In general, accuracy requirements necessitate extensive computation times. In celestial mechanics, regularizing transformations of the equations of motion are used to eliminate computational and analytical problems that occur during close approaches to gravitational force centers. It was shown in previous investigations that regularization in the formulation of the trajectory optimization problem may reduce the computation time. In this study, a set of regularized equations describing the optimal trajectory of a continuously thrusting space vehicle is derived. The computational characteristics of the set are investigated and compared to the classical Newtonian unregularized set of equations. The comparison is made for low thrust, minimum time, escape trajectories and numerical calculations of Keplerian orbits. The comparison indicates that in the cases investigated for bad initial guesses of the known boundary values a remarkable reduction in the computation time was achieved. Furthermore, the investigated set of regularized equations shows high numerical stability even for long duration flights and is less sensitive to errors in the guesses of the unknown boundary values.

Schwenzfeger, K. J.↗

NASA NTRS1982

Solution to the backward-Kolmogorov equation for a nonstationary oscillation problem

The transition probability density function of a Markovian approximation of the response amplitude of an oscillator under nonstationary excitation is determined in an analytical form. A solution is presented for the associated backward-Kolmogorov equation by transforming the equation into a form amenable to solution by the method of the separation of variables. This procedure results in a boundary value problem which is then solved by using an infinite series of Laguerre polynomials. It is found that the infinite series solution is equivalent to a closed-form solution involving a Bessel function.

Solomos, G. P.↗

NASA NTRS1987

Acoustic gravity waves: A computational approach

This paper discusses numerical solutions of a hyperbolic initial boundary value problem that arises from acoustic wave propagation in the atmosphere. Field equations are derived from the atmospheric fluid flow governed by the Euler equations. The resulting original problem is nonlinear. A first order linearized version of the problem is used for computational purposes. The main difficulty in the problem as with any open boundary problem is in obtaining stable boundary conditions. Approximate boundary conditions are derived and shown to be stable. Numerical results are presented to verify the effectiveness of these boundary conditions.

Hariharan, S. I.↗

NASA NTRS1988

Acoustic gravity waves - A computational approach

Hariharan, S. I.↗

NASA NTRS1984

On some properties of force-free magnetic fields in infinite regions of space

Techniques for solving boundary value problems (BVP) for a force free magnetic field (FFF) in infinite space are presented. A priori inequalities are defined which must be satisfied by the force-free equations. It is shown that upper bounds may be calculated for the magnetic energy of the region provided the value of the magnetic normal component at the boundary of the region can be shown to decay sufficiently fast at infinity. The results are employed to prove a nonexistence theorem for the BVP for the FFF in the spatial region. The implications of the theory for modeling the origins of solar flares are discussed.

Aly, J. J.↗

NASA NTRS1980

The control and estimation of large space structures

A shape control or estimation problem for a large space structure can be modeled by a partial differential equation which represents changes in shape with respect to space and time, together with spatially discrete forcing functions or observations which represent the placement of actuators or sensors at discrete points along the structure. The use of Green's functions to convert boundary value problems into integral equations provides a convenient treatment of this mixture of continuous and discrete mathematics. Control and estimation algorithms are developed for the one-dimensional static beam to illustrate this technique.

Weeks, C.↗

NASA NTRS1985

Stability of interfaces with mesh refinement

A study is conducted of the stability of mesh refinement in space and time for several different interface equations and finite-difference approximations. First, a root condition which implies stability for the initial-boundary value problem for this type of interface is derived. From the root condition, the stability of several interface equations is proved, using the maximum principle. In some cases, the final verification steps can be done analytically; in other cases, a simple computer program has been written to check the condition for values of a parameter along the boundary of the unit circle. Using this method, stability for Lax-Wendroff with all the interface conditions considered, and for Leapfrog with interpolation interface conditions when the fine and coarse grids overlap is proved.

Berger, M. J.↗

NASA NTRS1981

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

NASA NTRS1982

Solidification of a binary mixture

The time dependent concentration and temperature profiles of a finite layer of a binary mixture are investigated during solidification. The coupled time dependent Stefan problem is solved numerically using an implicit finite differencing algorithm with the method of lines. Specifically, the temporal operator is approximated via an implicit finite difference operator resulting in a coupled set of ordinary differential equations for the spatial distribution of the temperature and concentration for each time. Since the resulting differential equations set form a boundary value problem with matching conditions at an unknown spatial point, the method of invariant imbedding is used for its solution.

Antar, B. N.↗

NASA NTRS1984

Feedback laws for fuel minimization for transport aircraft

The Theoretical Mechanics Branch has as one of its long-range goals to work toward solving real-time trajectory optimization problems on board an aircraft. This is a generic problem that has application to all aspects of aviation from general aviation through commercial to military. Overall interest is in the generic problem, but specific problems to achieve concrete results are examined. The problem is to develop control laws that generate approximately optimal trajectories with respect to some criteria such as minimum time, minimum fuel, or some combination of the two. These laws must be simple enough to be implemented on a computer that is flown on board an aircraft, which implies a major simplification from the two point boundary value problem generated by a standard trajectory optimization problem. In addition, the control laws allow for changes in end conditions during the flight, and changes in weather along a planned flight path. Therefore, a feedback control law that generates commands based on the current state rather than a precomputed open-loop control law is desired. This requirement, along with the need for order reduction, argues for the application of singular perturbation techniques.

Price, D. B.↗

NASA NTRS1982

Integrating matrix solution of the hybrid state vector equations for beam vibration

A simple, versatile, and efficient computational technique has been developed for dynamic analysis of linear elastic beam and rod type of structures. Moreover, the method provides a rather general solution approach for two-point boundary value problems that are described by a single independent spatial variable. For structural problems, the method is implemented by a mixed state vector formulation of the differential equations, combined with an integrating matrix solution procedure. Highly accurate solutions are easily achieved with this approach. Example solutions are given for beam vibration problems including discontinuous stiffness and mass parameters, elastic restraint boundary conditions, concentrated inertia loading, and rigid body modes

Lehman, L. L.↗