On conforming finite element methods for the inhomogeneous stationary Navier-Stokes equations

Summary We consider the stationary Navier-Stokes equations, written in terms of the primitive variables, in the case where both the partial differential equations and boundary conditions are inhomogeneous. Under certain conditions on the data, the existence and uniqueness of the solution of a weak formulation of the equations can be guaranteed. A conforming finite element method is presented and optimal estimates for the error of the approximate solution are proved. In addition, the convergence properties of iterative methods for the solution of the discrete nonlinear algebraic systems resulting from the finite element algorithm are given. Numerical examples, using an efficient choice of finite element spaces, are also provided.Supported, in part, by the U.S. Air Force Office of Scientific Research under Grant No. AF-AFOSR-80-0083Supported, in part, by the same agency under Grant No. AF-AFOSR-80-0176-A. Both authors were also partially supported by NASA Contract No. NAS1-15810 while they were in residence at the Institute for Computer Applications in Science and Engineering, NASA Langley Research Center, Hampton, VA 23665, USA     

A Shadowing Theorem for ordinary differential equations

A new notion of shadowing of a pseudo orbit, an approximate solution, of an autonomous system of ordinary differential equations by an associated nearby true orbit is introduced. Then a general theorem which guarantees the existence of shadowing of pseudo orbits in compact hyperbolic sets is proved.Supported in part by the Air Force.Supported in part by NSF grant DMS 9201951.     

The empirical performance of a polynomial algorithm for constrained nonlinear optimization

In [6], a polynomial algorithm based on successive piecewise linear approximation was described. The algorithm is polynomial for constrained nonlinear (convex or concave) optimization, when the constraint matrix has a polynomial size subdeterminant. We propose here a practical adaptation of that algorithm with the idea of successive piecewise linear approximation of the objective on refined grids, and the testing of the gap between lower and upper bounds. The implementation uses the primal affine interior point method at each approximation step. We develop special features to speed up each step and to evaluate the gap. Empirical study of problems of size up to 198 variables and 99 constraints indicates that the procedure is very efficient and all problems tested were terminated after 171 interior point iterations. The procedure used in the implementation is proved to converge when the objective is strongly convex.Supported in part by the Office of Naval Research under Grant No. N00014-88-K-0377 and Grant No. ONR N00014-91-J-1241.     

Weak and Strong Solutions for the Stokes Approximation of Non-homogeneous Incompressible Navier-Stokes Equations

In this paper,the Dirichlet problem of Stokes approximate of non-homogeneous incompressibleNavier-Stokes equations is studied.It is shown that there exist global weak solutions as well as global andunique strong solution for this problem,under the assumption that initial density ρ_0(x)is bounded away from0 and other appropriate assumptions(see Theorem 1 and Theorem 2).The semi-Galerkin method is applied toconstruct the approximate solutions and a prior estimates are made to elaborate upon the compactness of theapproximate solutions.     

Critical point theorems for indefinite functionals

A variational principle of a minimax nature is developed and used to prove the existence of critical points for certain variational problems which are indefinite. The proofs are carried out directly in an infinite dimensional Hilbert space. Special cases of these problems previously had been tractable only by an elaborate finite dimensional approximation procedure. The main applications given here are to Hamiltonian systems of ordinary differential equations where the existence of time periodic solutions is established for several classes of Hamiltonians.Supported in part by the U.S. Army under Contract No. DAAG-29-75-C-0024 and by the Conseglio Nazionale delle Ricerche-Gruppo Nazionale Analisi Funzionale e ApplicazioneSupported in part by the J.S. Guggenheim Memorial Foundation, and by the Office of Naval Research under Contract No. N00014-76-C-0300. Reproduction in whole or in part is permitted for any purpose of the U.S. Government     

A relaxation method for nonconvex quadratically constrained quadratic programs

We present an algorithm for finding approximate global solutions to quadratically constrained quadratic programming problems. The method is based on outer approximation (linearization) and branch and bound with linear programming subproblems. When the feasible set is non-convex, the infinite process can be terminated with an approximate (possibly infeasible) optimal solution. We provide error bounds that can be used to ensure stopping within a prespecified feasibility tolerance. A numerical example illustrates the procedure. Computational experiments with an implementation of the procedure are reported on bilinearly constrained test problems with up to sixteen decision variables and eight constraints.This research was supported in part by National Science Foundation Grant DDM-91-14489.     

A general parametric analysis approach and its implication to sensitivity analysis in interior point methods

Adler and Monteiro (1992) developed a parametric analysis approach that is naturally related to the geometry of the linear program. This approach is based on the availability of primal and dual optimal solutions satisfying strong complementarity. In this paper, we develop an alternative geometric approach for parametric analysis which does not require the strong complementarity condition. This parametric analysis approach is used to develop range and marginal analysis techniques which are suitable for interior point methods. Two approaches are developed, namely the LU factorization approach and the affine scaling approach. Presented at the ORSA/TIMS, Nashville, TN, USA, May 1991. Supported by the National Science Foundation (NSF) under Grant No. DDM-9109404 and Grant No. DMI-9496178. This work was done while the author was a faculty member of the Systems and Industrial Engineering Department at The University of Arizona. Supported in part by the GTE Laboratories and the National Science Foundation (NSF) under Grant No. CCR-9019469.     

Approximate analysis of arbitrary configurations of open queueing networks with blocking

An algorithm for analyzing approximately open exponential queueing networks with blocking is presented. The algorithm decomposes a queueing network with blocking into individual queues with revised capacity, and revised arrival and service processes. These individual queues are then analyzed in isolation. Numerical experience with this algorithm is reported for three-node and four-node queueing networks. The approximate results obtained were compared against exact numerical data, and they seem to have an acceptable error level.Supported in part by a grant from CAIP Center, Rutgers University.Supported in part by the National Science Foundation under Grant DCR-85-02540.     

Linear stochastic differential equations with boundary conditions

Summary We study linear stochastic differential equations with affine boundary conditions. The equation is linear in the sense that both the drift and the diffusion coefficient are affine functions of the solution. The solution is not adapted to the driving Brownian motion, and we use the extended stochastic calculus of Nualart and Pardoux [16] to analyse them. We give analytical necessary and sufficient conditions for existence and uniqueness of a solution, we establish sufficient conditions for the existence of probability densities using both the Malliavin calculus and the co-aera formula, and give sufficient conditions that the solution be either a Markov process or a Markov field.Supported in part by NSF Grant No. MCS-8301880The research was carried out while this author was visiting the Institute for Advanced Study, Princeton NJ, and was supported by a grant from the RCA corporation     

A Uniformly-Valid Asymptotic Solution to a Matrix System of Ordinary Differential Equations and a Proof of its Validity

A new multiple-scale perturbation technique is employed to find the approximate solution to a fairly general matrix system of ordinary differential equations. This system includes a linear part given by a slowly-varying matrix and a small nonlinear part. The general proof of the method given in previous work is used to show rigorously that the present approximate solution is indeed asymptotic to the solution of the differential system. Some typical special cases of the general solution are also given.     

Initial-value methods for second-order singularly perturbed boundary-value problems

In this paper, we present a numerical method for solving linear and nonlinear second-order singularly perturbed boundary-value-problems. For linear problems, the method comes from the well-known WKB method. The required approximate solution is obtained by solving the reduced problem and one or two suitable initial-value problems, directly deduced from the given problem. For nonlinear problems, the quasilinearization method is applied. Numerical results are given showing the accuracy and feasibility of the proposed method.This work was supported in part by the Consiglio Nazionale delle Ricerche (Contract No. 86.02108.01 and Progetto Finalizzatto Sistemi Informatia e Calcolo Paralello, Sottoprogetto 1), and in part by the Ministero della Pubblica Istruzione, Rome, Italy.     

A finite element method for the Tricomi problem

Summary We consider the numerical solution of the Tricomi problem. Using a weak formulation based on different spaces of test and trial functions, we construct a new Galerkin procedure for the Tricomi problem. Existence, uniqueness, and uniform stability of the approximate solution is proven, and a priori error bounds are given.Research supported in part by the Department of Energy under contract DOE E(40-1)3443     

A class of singularly perturbed,nonlinear, fixed-endpoint control problems

Singular perturbation techniques are applied to a class of nonlinear, fixed-endpoint control problems to decompose the full-order problem into three lower-order problems, namely, the reduced problem and the left and right boundary-layer problems. The boundary-layer problems are linear-quadratic and, contrary to previous singular perturbation works, the reduced problem has a simple formulation. The solutions of these lower-order problems are combined to yield an approximate solution to the full nonlinear problem. Based on the properties of the lower-order problems, the full problem is shown to possess an asymptotic series solution.This work was supported in part by the National Science Foundation under Grant No. ENG-47-20091 and in part by the US Air Force under Grant No. AFOSR-73-2570.The author acknowledges the helpful suggestions of Professor P. V. Kokotovic, University of Illinois, Urbana, Illinois.     

Series Nash solution of two-person,nonzero-sum,linear-quadratic differential games

It is well-known that the Nash equilibrium solution of a two-person, nonzero-sum, linear differential game with a quadratic cost function can be expressed in terms of the solution of coupled generalized Riccati-type matrix differential equations. For high-order games, the numerical determination of the solution of the nonlinear coupled equations may be difficult or even impossible when the application dictates the use of small-memory computers. In this paper, a series solution is suggested by means of a parameter imbedding method. Instead of solving a high-order matrix-Riccati equation, a lower-order matrix-Riccati equation corresponding to a zero-sum game is solved. In addition, lower-order linear equations have to be solved. These solutions to lower-order equations are the coefficients of the series solution for the nonzero-sum game. Cost functions corresponding to truncated solutions are compared with those for exact Nash equilibrium solutions.This research was supported in part by the National Science Foundation under Grant No. GK-3893, in part by the Air Force under Grant No. AFOSR-68-1579B, and in part by the Joint Services Electronics Program under Contract No. DAAB-07-67-C-0199 with the Coordinated Science Laboratory, University of Illinois, Urbana, Illinois.     

Successive approximation method in optimum distributed-parameter systems

The problem of controlling a linear distributed-parameter system with a nonquadratic error measure is discussed. The calculus of variations approach is used to derive an algorithm based on the first variation for theN-dimensional linear diffusion process. The procedure for determining whether the resulting solution is optimum is discussed. Extension of the algorithm for other linear distributed-parameter systems is indicated.This research was supported in part by the National Science Foundation, Grant No. GK-304.     

A quadratically convergent O(\sqrt n L)-iteration algorithm for linear programming

Recently, Ye, Tapia and Zhang (1991) demonstrated that Mizuno—Todd—Ye's predictor—corrector interior-point algorithm for linear programming maintains the O( L)-iteration complexity while exhibiting superlinear convergence of the duality gap to zero under the assumption that the iteration sequence converges, and quadratic convergence of the duality gap to zero under the assumption of nondegeneracy. In this paper we establish the quadratic convergence result without any assumption concerning the convergence of the iteration sequence or nondegeneracy. This surprising result, to our knowledge, is the first instance of a demonstration of polynomiality and superlinear (or quadratic) convergence for an interior-point algorithm which does not assume the convergence of the iteration sequence or nondegeneracy.Supported in part by NSF Grant DDM-8922636 and NSF Coop. Agr. No. CCR-8809615, the Iowa Business School Summer Grant, and the Interdisciplinary Research Grant of the University of Iowa Center for Advanced Studies.Supported in part by NSF Coop. Agr. No. CCR-8809615, AFOSR 89-0363, DOE DEFG05-86ER25017 and ARO 9DAAL03-90-G-0093.Supported in part by NSF Grant DMS-9102761 and DOE Grant DE-FG05-91ER25100.     

Laguerre polynomial approach for solving linear delay difference equations

This paper presents a numerical method for the approximate solution of mth-order linear delay difference equations with variable coefficients under the mixed conditions in terms of Laguerre polynomials. The aim of this article is to present an efficient numerical procedure for solving mth-order linear delay difference equations with variable coefficients. Our method depends mainly on a Laguerre series expansion approach. This method transforms linear delay difference equations and the given conditions into matrix equation which corresponds to a system of linear algebraic equation. The reliability and efficiency of the proposed scheme are demonstrated by some numerical experiments and performed on the computer algebraic system Maple.     

On the existence of Nash strategies and solutions to coupled riccati equations in linear-quadratic games

The existence of linear Nash strategies for the linear-quadratic game is considered. The solvability of the coupled Riccati matrix equations and the stability of the closed-loop matrix are investigated by using Brower's fixed-point theorem. The conditions derived state that the linear closed-loop Nash strategies exist, if the open loop matrixA has a sufficient degree of stability which is determined in terms of the norms of the weighting matrices. WhenA is not necessarily stable, sufficient conditions for existence are given in terms of the solutions of auxiliary problems using the same procedure.This work was supported in part by the Joint Services Electronics Program (US Army, US Navy, and US Air Force) under Contract No. DAAG-29-78-C-0016, in part by the National Science Foundation under Grant No. ENG-74-20091, and in part by the Department of Energy, Electric Energy Systems Division, under Contract No. US-ERDA-EX-76-C-01-2088.