Search NASA⌕ Search

SEARCH · Search NASA

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Multiple Coarse Grid Multigrid Methods for Solving Elliptic Problems

In this paper we describe some classes of multigrid methods for solving large linear systems arising in the solution by finite difference methods of certain boundary value problems involving Poisson's equation on rectangular regions. If parallel computing systems are used, then with standard multigrid methods many of the processors will be idle when one is working at the coarsest grid levels. We describe the use of Multiple Coarse Grid MultiGrid (MCGMG) methods. Here one first constructs a periodic set of equations corresponding to the given system. One then constructs a set of coarse grids such that for each grid corresponding to the grid size h there are four grids corresponding to the grid size 2*h. Multigrid operations such as restriction of residuals and interpolation of corrections are done in parallel at each grid level. For suitable choices of the multigrid operators the MCGMG method is equivalent to the Parallel Superconvergent MultiGrid (PSMG) method of Frederickson and McBryan. The convergence properties of MCGMG methods can be accurately analyzed using spectral methods.

Xiao, Shengyou↗
NASA NTRS1974

Viscous-fluid reaction to a torsionally oscillating spheroid Linearized steady-state weak-separation solution

For both the linearized prolate and oblate problems, perturbation series for all orders of eccentricity are obtained. The 'weak separation' method of solving partial-differential-equation boundary-value problems assumes that each linearly independent solution is a sum of products of single-variable functions. The method is applied to a hierarchy of partial differential equations. The weak-separation solutions are constructed by generalizing only the Bessel functions of the homogeneous solution.

Tompkins, D. R., Jr.↗
NASA NTRS1990

The dynamics and control of large flexible space structures - 13

The optimal control of three-dimensional large angle maneuvers and vibrations of a Shuttle-mast-reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam subject to three-dimensional deformations. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem is then solved by using the quasilinearization algorithm and the method of particular solutions. The study of the large angle maneuvering of the Shuttle-beam-reflector spacecraft in the plane of a circular earth orbit is extended to consider the effects of the structural offset connection, the axial shortening, and the gravitational torque on the slewing motion. Finally the effect of additional design parameters (such as related to additional payload requirement) on the linear quadratic regulator based design of an orbiting control/structural system is examined.

Bainum, Peter M.↗
NASA NTRS1991

Maneuver simulations of flexible spacecraft by solving TPBVP

The optimal control of large angle rapid maneuvers and vibrations of a Shuttle mast reflector system is considered. The nonlinear equations of motion are formulated by using Lagrange's formula, with the mast modeled as a continuous beam. The nonlinear terms in the equations come from the coupling between the angular velocities, the modal coordinates, and the modal rates. Pontryagin's Maximum Principle is applied to the slewing problem, to derive the necessary conditions for the optimal controls, which are bounded by given saturation levels. The resulting two point boundary value problem (TPBVP) is then solved by using the quasilinearization algorithm and the method of particular solutions. In the numerical simulations, the structural parameters and the control limits from the Spacecraft Control Lab Experiment (SCOLE) are used. In the 2-D case, only the motion in the plane of an Earth orbit or the single axis slewing motion is discussed. In the 3-D slewing, the mast is modeled as a continuous beam subjected to 3-D deformations. The numerical results for both the linearized system and the nonlinear system are presented to compare the differences in their time response.

Bainum, Peter M.↗
NASA NTRS1991

Optimal Control Problems with Switching Points

The main idea of this report is to give an overview of the problems and difficulties that arise in solving optimal control problems with switching points. A brief discussion of existing optimality conditions is given and a numerical approach for solving the multipoint boundary value problems associated with the first-order necessary conditions of optimal control is presented. Two real-life aerospace optimization problems are treated explicitly. These are altitude maximization for a sounding rocket (Goddard Problem) in the presence of a dynamic pressure limit, and range maximization for a supersonic aircraft flying in the vertical, also in the presence of a dynamic pressure limit. In the second problem singular control appears along arcs with active dynamic pressure limit, which in the context of optimal control, represents a first-order state inequality constraint. An extension of the Generalized Legendre-Clebsch Condition to the case of singular control along state/control constrained arcs is presented and is applied to the aircraft range maximization problem stated above. A contribution to the field of Jacobi Necessary Conditions is made by giving a new proof for the non-optimality of conjugate paths in the Accessory Minimum Problem. Because of its simple and explicit character, the new proof may provide the basis for an extension of Jacobi's Necessary Condition to the case of the trajectories with interior point constraints. Finally, the result that touch points cannot occur for first-order state inequality constraints is extended to the case of vector valued control functions.

Seywald, Hans↗
NASA NTRS1966

Coupled thermoelasticity of bars.

Boundary value problems for semiinfinite bars made of coupled thermoelastic material solved by new functions that are corrections to classic tabulated functions

ERROR FUNCTION↗
NASA NTRS1986

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

Failure mechanisms of materials under very high strains experienced at and ahead of the crack tip such as formation, growth, and interaction of microvoids in ductile materials, microcracks in brittle solids or crazes in polymers and adhesives are represented by one-dimensional, nonlinear stress-strain relations possessing different ways by which the material loses capacity to carry load up to fracture or total separation. A double cantilever beam (DCB) type specimen is considered. The nonlinear material is confined to a thin strip between the two elastic beams loaded by a wedge. The problem is first modeled as a beam on a nonlinear foundation. The pertinent equation is solved numerically as a two-point boundary value problem for both the stationary and the quasi-stationay propagating crack. A finite element model is then used to model the problem in more detail in order to assess the adequacy of the beam model for the reduction of experimental data to determine in-situ properties of the thin interlayer.

Ungsuwarungsri, T.↗
NASA NTRS1979

Application of the pseudostate theory to the three-body Lambert problem

The pseudostate theory, which approximates three-body trajectories by overlapping the conic effects of both massive bodies on the third body, has been used to solve boundary value problems. Frequently, the approach to the secondary is quite close, as in interplanetary gravity assist trajectories or satellite tour trajectories. In this case the orbit with respect to the primary is radically changed so that perturbation techniques are time consuming, yet higher accuracy than point-to-point conics (V-infinity matching) is necessary. This method reduces the solution of the three-body Lambert problem to solving two conic Lambert problems and inverting a 7 x 7 matrix, the components of which are all found analytically. Typically 90-95% of the point-to-point conic error with respect to an integrated trajectory is eliminated.

Byrnes, D. V.↗
NASA NTRS1977

Numerical determination of the fundamental eigenvalue for the Laplace operator on a spherical domain

Methods for obtaining approximate solutions for the fundamental eigenvalue of the Laplace-Beltrami operator (i.e., the membrane eignevalue problem for the vibration equation) on the unit spherical surface are developed. Two types of spherical surface domains are considered: the interior of a spherical triangle, and the exterior of a great circle arc extending for less than pi radians (a spherical surface with a slit). In both cases, zero boundary conditions are imposed. In order to solve the resulting second-order elliptic partial differential equations in two independent variables, a finite difference approximation is employed. The fundamental eigenvalue is approximated by iteration utilizing the power method and point successive overrelaxation. Some numerical results are given and compared, in certain special cases, with analytical solutions to the eigenvalue problem. The significance of the numerical eigenvalue results is discussed in terms of the singularities in the solution of three-dimensional boundary-value problems near a polyhedral corner of the domain.

Walden, H.↗
NASA NTRS1986

Exact solution of three-dimensional transport problems using one-dimensional models

Several parameters of certain three-dimensional semiconductor devices including diodes, transistors, and solar cells can be determined without solving the actual boundary-value problem. The recombination current, transit time, and open-circuit voltage of planar diodes are emphasized here. The resulting analytical expressions enable determination of the surface recombination velocity of shallow planar diodes. The method involves introducing corresponding one-dimensional models having the same values of these parameters.

Misiakos, K.↗
NASA NTRS1989

Application of the pseudostate theory to the three-body Lambert problem

The pseudostate theory, which approximates three-body trajectories by overlapping the conic effects of both massive bodies on the third body, has been used to solve boundary-value problems. Frequently, the approach to the secondary is quite close, as in interplanetary gravity-assist or satellite-tour trajectories. In this case, the orbit with respect to the primary is radically changed so that perturbation techniques are time consuming, yet higher accuracy than point-to-point conics is necessary. This method reduces the solution of the three-body Lambert problem to solving two conic Lambert problems and inverting a 7 x 7 matrix, the components of which are all found analytically. Typically 90-95 percent of the point-to-point conic error, with respect to an integrated trajectory, is eliminated.

Byrnes, Dennis V.↗
NASA NTRS1996

Numerical Conformal Mapping Using Cross-Ratios and Delaunay Triangulation

We propose a new algorithm for computing the Riemann mapping of the unit disk to a polygon, also known as the Schwarz-Christoffel transformation. The new algorithm, CRDT, is based on cross-ratios of the prevertices, and also on cross-ratios of quadrilaterals in a Delaunay triangulation of the polygon. The CRDT algorithm produces an accurate representation of the Riemann mapping even in the presence of arbitrary long, thin regions in the polygon, unlike any previous conformal mapping algorithm. We believe that CRDT can never fail to converge to the correct Riemann mapping, but the correctness and convergence proof depend on conjectures that we have so far not been able to prove. We demonstrate convergence with computational experiments. The Riemann mapping has applications to problems in two-dimensional potential theory and to finite-difference mesh generation. We use CRDT to produce a mapping and solve a boundary value problem on long, thin regions for which no other algorithm can solve these problems.

Driscoll, Tobin A.↗
NASA NTRS1995

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.↗
NASA NTRS1984

Stress analysis of multilayered plates around circular holes

The elasticity problem for a thick plate consisting of two bonded dissimilar homogeneous layers is considered. It is assumed that the plate is infinite, contains a circular hole, and is subjected to axisymmetric external loads. A technique suitable to the formulation of three-dimensional boundary value problems is introduced and applied to the three-dimensional elasticity solution of an infinite plate with a hole. The effect of Poisson's ratio is studied in some detail, and the axisymmetric problem of a laminated plate is solved. The stresses at the interface are calculated and are observed to become unbounded at the hole boundary. An approximate treatment of the singularity is presented and the stress intensity factors are computed. Some of the results are compared with those found from the finite element method. The agreement between the two sets with regard to trend and magnitude is found to be very satisfactory.

Delale, F.↗
NASA NTRS1979

Multi-level adaptive computations in fluid dynamics

The multi-level adaptive technique (MLAT) is a general strategy of solving continuous problems by cycling between coarser and finer levels of discretization. It provides very fast solvers together with adaptive, nearly optimal discretization schemes to general boundary-value problems in general domains. Here the state of the art is surveyed, emphasizing steady-state fluid dynamics applications, from slow viscous flows to transonic ones. Various new techniques are briefly discussed, including distributive relaxation schemes, the treatment of evolution problems, the combined use of upstream and central differencing, local truncation extrapolations, and other 'super-solver' techniques.

Brandt, A.↗
NASA NTRS1981

Accuracy and convergence of a finite element algorithm for laminar boundary layer flow

The Galerkin-weighted residuals formulation is employed to derive an implicit finite element solution algorithm for a generally non-linear initial-boundary value problem. Solution accuracy and convergence with discretization refinement are quantized in several error norms, for the non-linear parabolic partial differential equation system governing laminar boundary layer flow, using linear, quadratic and cubic functions. Richardson extrapolation is used to isolate integration truncation error in all norms, and Newton iteration is employed for all equation solutions performed in double-precision. The mathematical theory supporting accuracy and convergence concepts for linear elliptic equations appears extensible to the non-linear equations characteristic of laminar boundary layer flow.

Soliman, M. O.↗
NASA NTRS1985

Convergence to steady state of solutions of Burgers' equation

Consider the initial boundary value problem for Burgers' equation. It is shown that its solutions converge, in time, to a unique steady state. The speed of the convergence depends on the boundary conditions and can be exponentially slow. Methods to speed up the rate of convergence are also discussed.

Kreiss, G.↗
NASA NTRS1987

Current distribution in a three-dimensional IC analyzed by a perturbation method. Part 1: A simple steady state theory

The steady state current distribution in a three dimensional integrated circuit is presented. A device physics approach, based on a perturbation method rather than an equivalent lumped circuit approach, is used. The perturbation method allows the various currents to be expressed in terms of elementary solutions which are solutions to very simple boundary value problems. A Simple Steady State Theory is the subtitle because the most obvious limitation of the present version of the analysis is that all depletion region boundary surfaces are treated as equipotential surfaces. This may be an adequate approximation in some applications but it is an obvious weakness in the theory when applied to latched states. Examples that illustrate the use of these analytical methods are not given because they will be presented in detail in the future.

Edmonds, Larry D.↗