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 451 records · Page 25

All sources NASA NTRS DOE OSTI

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

NASA NTRS2003

Further Developments in the MLPG Method for Beam Problems

An accurate and yet simple Meshless Local Petrov-Galerkin (MLPG) formulation for analyzing beam problems is presented. In the formulation, simple weight functions are chosen as test functions as in the conventional MLPG method. Linear test functions are also chosen, leading to a variation of the MLPG method that is computationally efficient compared to the conventional implementation. The MLPG method is evaluated by applying the formulation to a variety of patch tests, thin beam problems, and problems with load discontinuities. The formulation successfully reproduces exact solutions to machine accuracy when higher order power and spline functions are chosen as test functions or when the linear test function is used, and when constructing the trial functions, the order of the basis function is properly balanced by the order of the weight function. For mixed boundary value problems, deflections, slopes, moments, and shear forces are calculated to the same accuracy by the MLPG method without the use of elaborate post-processing techniques. Problems with load discontinuities require special care - when a reasonable number of nodes are used, the method yields very accurate results.

Raju, I. S.↗

NASA NTRS1992

A two-dimensional MHD global coronal model - Steady-state streamers

A 2D, time-dependent, numerical, MHD model for the simulation of coronal streamers from the solar surface to 15 solar is presented. Three examples are given; for dipole, quadrupole and hexapole (Legendre polynomials P1, P2, and P3) initial field topologies. The computed properties are density, temperature, velocity, and magnetic field. The calculation is set up as an initial-boundary value problem wherein a relaxation in time produces the steady state solution. In addition to the properties of the solutions, their accuracy is discussed. Besides solutions for dipole, quadrupole, and hexapole geometries, the model use of realistic values for the density and Alfven speed while still meeting the requirement that the flow speed be super-Alfvenic at the outer boundary by extending the outer boundary to 15 solar radii.

Wang, A.-H.↗

NASA NTRS2008

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗

NASA NTRS1988

Singular perturbation analysis of the atmospheric orbital plane change problem

A three-state model is presented for the aeroassisted orbital plane change problem. A further model order reduction to a single state model is examined using singular perturbation theory. The optimal solution for this single state model compares favorably with the exact numerical solution using a four-state model; however, a separate boundary layer solution is required to satisfy the terminal constraint on altitude. This, in general, involves the solution of a two-point boundary value problem, but for a two-state model. An approximation is introduced to obtain an analytical control solution for lift and bank angle. Included are numerical simulation results of a guidance law derived from this analysis, along with comparison to earlier work by other researchers.

Calise, A. J.↗

NASA NTRS1996

First-Order System Least-Squares for Second-Order Elliptic Problems with Discontinuous Coefficients

The first-order system least-squares methodology represents an alternative to standard mixed finite element methods. Among its advantages is the fact that the finite element spaces approximating the pressure and flux variables are not restricted by the inf-sup condition and that the least-squares functional itself serves as an appropriate error measure. This paper studies the first-order system least-squares approach for scalar second-order elliptic boundary value problems with discontinuous coefficients. Ellipticity of an appropriately scaled least-squares bilinear form of the size of the jumps in the coefficients leading to adequate finite element approximation results. The occurrence of singularities at interface corners and cross-points is discussed. and a weighted least-squares functional is introduced to handle such cases. Numerical experiments are presented for two test problems to illustrate the performance of this approach.

Manteuffel, Thomas A.↗

NASA NTRS1980

A mixed-mode crack analysis of rectilinear anisotropic solids using conservation laws of elasticity

A very simple and convenient method of analysis for studying two-dimensional mixed-mode crack problems in rectilinear anisotropic solids is presented. The analysis is formulated on the basis of conservation laws of anisotropic elasticity and of fundamental relationships in anisotropic fracture mechanics. The problem is reduced to a system of linear algebraic equations in mixed-mode stress intensity factors. One of the salient features of the present approach is that it can determine directly the mixed-mode stress intensity solutions from the conservation integrals evaluated along a path removed from the crack-tip region without the need of solving the corresponding complex near-field boundary value problem. Several examples with solutions available in the literature are solved to ensure the accuracy of the current analysis. This method is further demonstrated to be superior to other approaches in its numerical simplicity and computational efficiency. Solutions of more complicated and practical engineering problems dealing with the crack emanating from a circular hole in composites are presented also to illustrate the capacity of this method.

Wang, S. S.↗

NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗

NASA NTRS1984

Stability of finite difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaires or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗

NASA NTRS1987

The dynamics and control of large flexible space structures X, part 1

The effect of delay in the control system input on the stability of a continuously acting controller which is designed without considering the delay is studied. The stability analysis of a second order plant is studied analytically and verified numerically. For this example it is found that the system becomes unstable for a delay which is equivalent to only 16 percent of its natural period of motion. It is also observed that even a small amount of natural damping in the system can increase the amount of delay that can be tolerated before the onset of instability. The delay problem is formulated in the discrete time domain and an analysis procedure suggested. The maximum principle from optimal control theory is applied to minimize the time required for the slewing of a general rigid spacecraft. The slewing motion need not be restricted to a single axis maneuver. The minimum slewing time is calculated based on a quasi-linearization algorithm for the resulting two point boundary value problem. Numerical examples based on the rigidized in-orbit model of the SCOLE also include the more general reflector line-of-sight slewing maneuvers.

Bainum, Peter M.↗

NASA NTRS1975

Constrained optimal controller for linear systems with state and control dependent disturbance

The problem is posed with the additional constraints that the dynamic controller uses only noise-corrupted outputs, and that its dimension is significantly lower than that of a Kalman filter. The unknown disturbance is viewed as an adversary which tries to maximize a performance criterion: a criterion that the controller gains attempt to minimize. The optimal controller gains are determined by solving a nonlinear matrix two-point boundary value problem.

Basuthakur, S.↗

NASA NTRS1979

Multi-level adaptive finite element methods. 1: Variation problems

A general numerical strategy for solving partial differential equations and other functional problems by cycling between coarser and finer levels of discretization is described. Optimal discretization schemes are provided together with very fast general solvers. It is described in terms of finite element discretizations of general nonlinear minimization problems. The basic processes (relaxation sweeps, fine-grid-to-coarse-grid transfers of residuals, coarse-to-fine interpolations of corrections) are directly and naturally determined by the objective functional and the sequence of approximation spaces. The natural processes, however, are not always optimal. Concrete examples are given and some new techniques are reviewed. Including the local truncation extrapolation and a multilevel procedure for inexpensively solving chains of many boundary value problems, such as those arising in the solution of time-dependent problems.

Brandt, A.↗