Accuracy of least-squares methods for the Navier-Stokes equations

Recently there has been substantial interest in least-squares finite element methods for velocity-vorticity-pressure formulations of the incompressible Navier-Stokes equations. The main cause for this interest is the fact that algorithms for the resulting discrete equations can be devised which require the solution of only symmetric, positive definite systems of algebraic equations. On the other hand, it is well-documented that methods using the vorticity as a primary variable often yield very poor approximations. Thus, here we study the accuracy of these methods through a series of computational experiments, and also comment on theoretical error estimates. It is found, despite the failure of standard methods for deriving error estimates, that computational evidence suggests that these methods are, at the least, nearly optimally accurate. Thus, in addition to the desirable matrix properties yielded by least-squares methods, one also obtains accurate approximations

Prediction of impact force and duration during low velocity impact on circular composite laminates

Two simple and improved models--energy-balance and spring-mass--were developed to calculate impact force and duration during low velocity impact of circular composite plates. Both models include the contact deformation of the plate and the impactor as well as bending, transverse shear, and membrane deformations of the plate. The plate was transversely isotropic graphite/epoxy composite laminate and the impactor was a steel sphere. Calculated impact forces from the two analyses agreed with each other. The analyses were verified by comparing the results with reported test data

A dynamic inverse for nonlinear maps

We consider the problem of estimating the time-varying root of a time-dependent nonlinear map. We introduce a "dynamic inverse" of a map, another generally time-dependent map which one composes with the original map to form a nonlinear vector-field. The flow of this vector field decays exponentially to the root. We then show how a dynamic inverse may be determined dynamically while being used simultaneously to find a root. We construct a continuous-time analog computational paradigm around the dynamic inverse

A solution procedure for behavior of thick plates on a nonlinear foundation and postbuckling behavior of long plates

Approximate solutions for three nonlinear orthotropic plate problems are presented: (1) a thick plate attached to a pad having nonlinear material properties which, in turn, is attached to a substructure which is then deformed; (2) a long plate loaded in inplane longitudinal compression beyond its buckling load; and (3) a long plate loaded in inplane shear beyond its buckling load. For all three problems, the two dimensional plate equations are reduced to one dimensional equations in the y-direction by using a one dimensional trigonometric approximation in the x-direction. Each problem uses different trigonometric terms. Solutions are obtained using an existing algorithm for simultaneous, first order, nonlinear, ordinary differential equations subject to two point boundary conditions. Ordinary differential equations are derived to determine the variable coefficients of the trigonometric terms

Difference Methods and Deferred Corrections for Ordinary Boundary Value Problems

Compact as possible difference schemes for systems of nth order equations are developed. Generalizations of the Mehrstellenverfahren and simple theoretically sound implementations of deferred corrections are given. It is shown that higher order systems are more efficiently solved as given rather than as reduced to larger lower order systems. Tables of coefficients to implement these methods are included and have been derived using symbolic computations

Numerical optimization in Hilbert space using inexact function and gradient evaluations

Trust region algorithms provide a robust iterative technique for solving non-convex unstrained optimization problems, but in many instances it is prohibitively expensive to compute high accuracy function and gradient values for the method. Of particular interest are inverse and parameter estimation problems, since function and gradient evaluations involve numerically solving large systems of differential equations. A global convergence theory is presented for trust region algorithms in which neither function nor gradient values are known exactly. The theory is formulated in a Hilbert space setting so that it can be applied to variational problems as well as the finite dimensional problems normally seen in trust region literature. The conditions concerning allowable error are remarkably relaxed: relative errors in the gradient error condition is automatically satisfied if the error is orthogonal to the gradient approximation. A technique for estimating gradient error and improving the approximation is also presented

Modeling of Overland Flow by the Diffusion Wave Approach

One of the major issues of present times, i.e. water quality degradation and a need for precise answers to transport of pollutants by overland flow, is addressed with special reference to the evaporator pits located adjacent to streams in the oil-producing regions of Eastern Kentucky. The practical shortcomings of the state-of-the-art kinematic wave are discussed and a new mathematical modeling-approach for overland flows using the more comprehensive diffusion wave is attempted as the first step in solving this problem. A Fourier series representation of the solution to the diffusion wave is adopted and found to perform well. The physically justified boundary conditions for steep slopes is considered and both numerical and analytical schemes are developed. The zero-depth-gradient lower condition is used and found to be adequate. The steady state analysis for mild slopes is found to be informative and both analytical and numerical solutions are found. The effect of imposing transients on the steady state solution are considered. Finally the cases for which these techniques can be used are presented

Solving two-point seismic-ray tracing problems in a heterogeneous medium. Part 1. A general adaptive finite difference method

A study of two-point seismic-ray tracing problems in a heterogeneous isotropic medium and how to solve them numerically will be presented in a series of papers. In this Part 1, it is shown how a variety of two-point seismic-ray tracing problems can be formulated mathematically as systems of first-order nonlinear ordinary differential equations subject to nonlinear boundary conditions. A general numerical method to solve such systems in general is presented and a computer program based upon it is described. High accuracy and efficiency are achieved by using variable order finite difference methods on nonuniform meshes which are selected automatically by the program as the computation proceeds. The variable mesh technique adapts itself to the particular problem at hand, producing more detailed computations where they are needed, as in tracing highly curved seismic rays. A complete package of programs has been produced which use this method to solve two- and three-dimensional ray-tracing problems for continuous or piecewise continuous media, with the velocity of propagation given either analytically or only at a finite number of points. These programs are all based on the same core program, PASVA3, and therefore provide a compact and flexible tool for attacking ray-tracing problems in seismology. In Part 2 of this work, the numerical method is applied to two- and three-dimensional velocity models, including models with jump discontinuities across interfaces

Weak Sharp Minima on Riemannian Manifolds

This is the first paper dealing with the study of weak sharp minima for constrained optimization problems on Riemannian manifolds, which are important in many applications. We consider the notions of local weak sharp minima, boundedly weak sharp minima, and global weak sharp minima for such problems and obtain their complete characterizations in the case of convex problems on finite-dimensional Riemannian manifolds and their Hadamard counterparts. A number of the results obtained in this paper are also new for the case of conventional problems in linear spaces. Our methods involve appropriate tools of variational analysis and generalized differentiation on Riemannian and Hadamard manifolds developed and efficiently implemented in this paper