Existence and boundary behavior for singular nonlinear differential equations with arbitrary boundary conditions

Existence theory is developed for the equation ?(u)=F(u), where ? is a formally self-adjoint singular second-order differential expression and F is nonlinear. The problem is treated in a Hilbert space and we do not require the operators induced by ? to have completely continuous resolvents. Nonlinear boundary conditions are allowed. Also, F is assumed to be weakly continuous and monotone at one point. Boundary behavior of functions associated with the domains of definitions of the operators associated with ? in the singular case is investigated. A special class of self-adjoint operators associated with ? is obtained.     

Lower-dimensional linear complementarity problem approaches to the solution of a bi-obstacle problem

A globally convergent Broyden-like method for solving a bi-obstacle problem is proposed based on its equivalent lower-dimensional linear complementarity problem. A suitable line search technique is introduced here. The global and superlinear convergence of the method is verified under appropriate assumptions.     

A second order numerical scheme for the solution of the one-dimensional Boussinesq equation

A predictor–corrector (P-C) scheme is applied successfully to a nonlinear method arising from the use of rational approximants to the matrix-exponential term in a three-time level recurrence relation. The resulting nonlinear finite-difference scheme, which is analyzed for local truncation error and stability, is solved using a P-C scheme, in which the predictor and the corrector are explicit schemes of order 2. This scheme is accelerated by using a modification (MPC) in which the already evaluated values are used for the corrector. The behaviour of the P-C/MPC schemes is tested numerically on the Boussinesq equation already known from the bibliography free of boundary conditions. The numerical results are derived for both the bad and the good Boussinesq equation and conclusions from the relevant known results are derived.      

An alternative approach for non-linear optimal control problems based on the method of moments

We propose an alternative method for computing effectively the solution of non-linear, fixed-terminal-time, optimal control problems when they are given in Lagrange, Bolza or Mayer forms. This method works well when the nonlinearities in the control variable can be expressed as polynomials. The essential of this proposal is the transformation of a non-linear, non-convex optimal control problem into an equivalent optimal control problem with linear and convex structure. The method is based on global optimization of polynomials by the method of moments. With this method we can determine either the existence or lacking of minimizers. In addition, we can calculate generalized solutions when the original problem lacks of minimizers. We also present the numerical schemes to solve several examples arising in science and technology.     

A new alternating‐direction finite element method for hyperbolic equation

A modified backward difference time discretization is presented for Galerkin approximations for nonlinear hyperbolic equation in two space variables. This procedure uses a local approximation of the coefficients based on patches of finite elements with these procedures, a multidimensional problem can be solved as a series of one‐dimensional problems. Optimal order H₀¹ and L² error estimates are derived. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Finite element computations of two-dimensional arterial flow in the presence of a transverse magnetic field

A finite element solution of the Navier-Stokes equations for steady flow under the magnetic effect through a double-branched two-dimensional section of a three-dimensional model of the canine aorta is discussed. The numerical scheme involves transforming the physical co-ordinates to a curvilinear boundary-fitted co-ordinate system. The shear stress at the wall is calculated for a Reynolds number of 1000 with the branch-to-main aortic flow rate ratio as a parameter. The results are compared with earlier works involving experimental data and found to be in reasonable qualitative agreement. The steady flow, shear stress and branch flow under the effect of a magnetic field have been discussed in detail.     

Modified Householder iterative method free from second derivatives for nonlinear equations

In this paper, we suggest and analyze a new two-step predictor–corrector type iterative method free from second derivatives for solving nonlinear equations of the type f(x)=0. This new method includes the two-step Newton method as a special case. We prove that the new iterative method is of fourth-order. Several examples are given to illustrate the efficiency of this new method and its comparison with other iterative methods. This method can be considered as a significant improvement of the Newton method and its variant forms.     

Modified Kirchhoff flow with a partially penetrable obstacle and its application to the efficiency of free flow turbines

An explicitly solvable analog of the Kirchhoff flow for the case of a semipenetrable obstacle is considered. Its application to estimating the efficiency of free flow turbines is discussed.     

A combined programming and iteration algorithm for finite element analysis of three-dimensional contact problems

Comparing with two-dimensional contact problems, three-dimensional frictional contact problems are more difficult to deal with, because of the unknown slip direction of the tangential force and enormous computing time. In order to overcome these difficulties, a combined PQP (Parametric Quadratic Programming) and iteration method is derived in this paper. The iteration algorithm, which alleviates the difficulty of unknown slip direction, is used along with the PQP method to cut down computing costs. Numerical example is given to demonstrate the validity of the present algorithm. The project supported by the Machinary and Electronics Ministry of China     

ROMBERG方法的改进

本文研究了全部根构成一个等比数列的多项式的性质，并将所得的结果用到数值积分的Ｒｏｍｂｅｒｇ方法＾「１」，得到了相应遥改进。