Incremental Unknowns Method for Solving Three-Dimensional Convection-Diffusion Equations

We use the incremental unknowns method in conjunction with the iterative methods to approximate the solution of the nonsymmetric and positive-definite linear systems generated from a multilevel discretization of three-dimensional convection-diffusion equations. The condition numbers of incremental unknowns matrices associated with the convection-diffusion equations and the number of iterations needed to attain an acceptable accuracy are estimated. Numerical results are presented with two-level approximations, which demonstrate that the incremental unknowns method when combined with some iter- ative methods is very effcient.     

Preconditioning analysis of the one dimensional incremental unknowns method on nonuniform meshes

The condition number of the incremental unknowns matrix on nonuniform meshes associated to the elliptic problem is analyzed. Comparing to the usual nodal unknowns matrix, the condition number of the incremental unknowns matrix is reduced significantly even if the meshes are nonuniform. Furthermore, if a diagonal scaling is used, the condition number of the preconditioned incremental unknowns matrix comes out to be O(1). Numerical experiments are performed respectively on Shishkin mesh and Chebyshev mesh. Computational results with respect to the two particular nonuniform meshes confirm our theoretical analysis.     

Connection between the spectral problem for linear matrix pencils and some problems of algebra

We suggest a method by which the solution of systems of nonlinear algebraic equations in one and two variables can be reduced to the spectral problem for linear pencils of two matrices and for a system of two matrix pencils of two matrices, respectively. This method is substantially different from the traditional methods of solution; at the same time it is useful for the study of the solvability and the determination of the number of solutions of such systems. We propose a modification of the well-known elimination method for the solution of nonlinear algebraic systems of two equations in two unknowns which provides a new approach to studying and solving the problem.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 80, pp. 98–116, 1978.     

A simple algebraic approximation to the Erlang loss system

We use Palm calculus to derive a simple, intuitive system of two linear-quadratic equations and two unknowns, whose algebraic solution yields Harel’s (1988) upper bound to the Erlang loss probability. We then derive a sequence of progressively stronger systems of equations, which eventually become exact. We provide two example applications.     

An existence result for nonisothermal immiscible incompressible 2‐phase flow in heterogeneous porous media

In the present article, we study the temperature effects on two‐phase immiscible incompressible flow through a porous medium. The mathematical model is given by a coupled system of 2‐phase flow equations and an energy balance equation. The model consists of the usual equations derived from the mass conservation of both fluids along with the Darcy‐Muskat and the capillary pressure laws. The problem is written in terms of the phase formulation; ie, the saturation of one phase, the pressure of the second phase, and the temperature are primary unknowns. The major difficulties related to this model are in the nonlinear degenerate structure of the equations, as well as in the coupling in the system. Under some realistic assumptions on the data, we show the existence of weak solutions with the help of an appropriate regularization and a time discretization. We use suitable test functions to obtain a priori estimates. We prove a new compactness result to pass to the limit in nonlinear terms.     

Global existence for systems of wave equations with nonresonant nonlinearities and null forms

We consider the Cauchy problem for systems of nonlinear wave equations with different propagation speeds in three space dimensions. We prove global existence of small amplitude solutions for systems with some nonresonant nonlinearities which may depend on both of the unknowns and their derivatives. Our method here can be also adopted to treat the null forms.     

Asymmetric invariants for a class of strictly hyperbolic systems including the Timoshenko beam

We introduce a set of conserved quantities of energy‐type for a strictly hyperbolic system of two coupled wave equations in one space dimension. The system is subject to mechanical boundary conditions. Some of these invariants are asymmetric in the sense that their defining quadratic form contains second order derivatives in only one of the unknowns. We study their independence with respect to the usual energies and characterize their sign. In many cases, our results provide sharp well‐posedness and stability results. Finally, we apply some of our conservation laws to the study of a singular perturbation problem previously considered by J. Lagnese and J. L. Lions. Copyright © 2007 John Wiley & Sons, Ltd.     

The convergence of V-cycle multigrid algorithms for axisymmetric Laplace and Maxwell equations

We investigate some simple finite element discretizations for the axisymmetric Laplace equation and the azimuthal component of the axisymmetric Maxwell equations as well as multigrid algorithms for these discretizations. Our analysis is targeted at simple model problems and our main result is that the standard V-cycle with point smoothing converges at a rate independent of the number of unknowns. This is contrary to suggestions in the existing literature that line relaxations and semicoarsening are needed in multigrid algorithms to overcome difficulties caused by the singularities in the axisymmetric Maxwell problems. Our multigrid analysis proceeds by applying the well known regularity based multigrid theory. In order to apply this theory, we prove regularity results for the axisymmetric Laplace and Maxwell equations in certain weighted Sobolev spaces. These, together with some new finite element error estimates in certain weighted Sobolev norms, are the main ingredients of our analysis.

     

A new approach to improve ill-conditioned parabolic optimal control problem via time domain decomposition

In this paper we present a new steepest-descent type algorithm for convex optimization problems. Our algorithm pieces the unknown into sub-blocs of unknowns and considers a partial optimization over each sub-bloc. In quadratic optimization, our method involves Newton technique to compute the step-lengths for the sub-blocs resulting descent directions. Our optimization method is fully parallel and easily implementable, we first presents it in a general linear algebra setting, then we highlight its applicability to a parabolic optimal control problem, where we consider the blocs of unknowns with respect to the time dependency of the control variable. The parallel tasks, in the last problem, turn “on” the control during a specific time-window and turn it “off” elsewhere. We show that our algorithm significantly improves the computational time compared with recognized methods. Convergence analysis of the new optimal control algorithm is provided for an arbitrary choice of partition. Numerical experiments are presented to illustrate the efficiency and the rapid convergence of the method.     

La relation de possion pour l'equation des ondes dans un ouvert non borne application a la theorie de la diffusion

We study in this note the solutions of two types of hyperbolic systems of conservation laws with oscillating data. The first one (a2?2 system) has only one linearly degenerate eigenvalue. Using the results of R.J.Di Perna related to genuinely nonlinear fields, one can describe the propagatrion of oscillations (which appear in only one direction) with an integro differential system for which one of the two unknowns is a field depending, of y ? ]0,1[, in addition of x and t. The second system is a linearly degenerate 3?3 system. We apply theory of compensated compactness, due to L. Tartar and F. Murat and, in the same way as above, we show that the initial oscillations can propagate; this propagation is then described witha a relaxed system of 3 unknowns     

Analysis of the coupling of BEM, FEM and mixed-FEM for a two-dimensional fluid–solid interaction problem

In this paper we consider a fluid–solid interaction problem posed in the plane. We employ a mixed variational formulation in the obstacle, in which the Cauchy stress tensor and the rotation are the only unknowns. This new mixed formulation is coupled, through suitable transmission conditions on the wet interface, with a Helmholtz equation satisfied by the pressure of the fluid in the unbounded domain. We use a traditional primal variational formulation in this part of the domain and incorporate the far field information through boundary integral equations. We approximate the resulting weak formulation by a Galerkin scheme based on PEERS in the solid and on a FEM-BEM approach in the fluid part. We show that our scheme is uniquely solvable and convergent, and then provide optimal error estimates. Finally, we illustrate our analysis with some computational experiments.     

Problems and methods of averaged optimization

Typical extremal problems containing the mean values of unknowns or of some functions of unknowns are considered. Relationships between these problems and cyclic modes of dynamical systems are revealed, and optimality conditions for such modes are found.     

Study of a numerical scheme for miscible two‐phase flow in porous media

We study the convergence of a finite volume scheme for a model of miscible two‐phase flow in porous media. In this model, one phase can dissolve into the other one. The convergence of the scheme is proved thanks to an estimate on the two pressures, which allows to prove some estimates on the discrete time derivative of some nonlinear functions of the unknowns. Monotony arguments allow to show some properties on the limits of these functions. A key point in the scheme is to use particular averaging formula for the dissolution function arising in the space term. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 723–748, 2014     

Polygonal heat conductors with a stationary hot spot

We consider a convex polygonal heat conductor whose inscribed circle touches every side of the conductor. Initially, the conductor has constant temperature and, at every time, the temperature of its boundary is kept at zero. The hot spot is the point at which temperature attains its maximum at each given time. It is proved that, if the hot spot is stationary, then the conductor must satisfy two geometric conditions. In particular, we prove that these geometric conditions yield some symmetries provided the conductor is either pentagonal or hexagonal. This research was partially supported by Grants-in-Aid for Scientific Research (B) (# 12440042) and (B) (# 15340047) of Japan Society for the Promotion of Science, and by a Grant of the Italian MURST.     

Optimal second order rectangular elasticity elements with weakly symmetric stress

We present new second order rectangular mixed finite elements for linear elasticity where the symmetry condition on the stress is imposed weakly with a Lagrange multiplier. The key idea in constructing the new finite elements is enhancing the stress space of the Awanou’s rectangular elements (rectangular Arnold–Falk–Winther elements) using bubble functions. The proposed elements have only 18 and 63 degrees of freedom for the stress in two and three dimensions, respectively, and they achieve the optimal second order convergence of errors for all the unknowns. We also present a new simple a priori error analysis and provide numerical results illustrating our analysis.     

Analysis of lumped models with contact and friction

We consider two mathematical models that describe the vibrations of spring-mass-damper systems with contact and friction. In the first model, both the contact and frictional boundary conditions are described with subdifferentials of nonconvex functions. In the second model, the contact is modeled with a Lipschitz continuous function, and the restitution force is described by a differential equation involving a Volterra integral term. The two models lead to second-order differential inclusions with and without an integral term, in which the unknowns are the positions of the masses. For each model, we prove the existence of a solution by using an abstract result for first-order differential inclusions in finite dimensional spaces. For the second model, in addition, we prove the uniqueness of the solution by using a fixed point argument. Finally, we provide examples of systems with contact and friction conditions for which our results are valid.     

Interpreting linear systems of equalities and inequalities. Application to the water supply problem

This paper shows how the mathematical and the engineering points of view are complementary and help to model real problems that can be stated as systems of linear equations and inequalities. The paper is devoted to point out these relations and making them explicit for the readers to realize about the new possibilities that arise when contemplating the compatibility conditions or the set of general solutions from the dual perspective. After reviewing an orthogonally based powerful algorithm to analyse the compatibility of linear systems of equations and solving them, a water supply problem is used to illustrate its mathematical and engineering multiple aspects, including the optimal statement of the problem in terms of an adequate selection and numbering of equations and unknowns, an analysis of the compatibility conditions and a physical interpretation of the general solution, together with that of each individual generators of the affine space. The possibilities of removing unknowns without altering the compatibility of the problem is also analysed. Next, the Γ‐algorithm to analyse the compatibility of linear systems of inequalities and solving them is described and then, the water supply problem is revisited adding some constraints, such as capacity limits for the pipes and retention valves, and discussed as to how they affect the resulting general solution and other aspects. Finally, some conclusions are derived. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of lumped models with contact and friction

We consider two mathematical models that describe the vibrations of spring-mass-damper systems with contact and friction. In the first model, both the contact and frictional boundary conditions are described with subdifferentials of nonconvex functions. In the second model, the contact is modeled with a Lipschitz continuous function, and the restitution force is described by a differential equation involving a Volterra integral term. The two models lead to second-order differential inclusions with and without an integral term, in which the unknowns are the positions of the masses. For each model, we prove the existence of a solution by using an abstract result for first-order differential inclusions in finite dimensional spaces. For the second model, in addition, we prove the uniqueness of the solution by using a fixed point argument. Finally, we provide examples of systems with contact and friction conditions for which our results are valid.     

Non-isothermal lubrication problem with Tresca fluid–solid interface law. Part II Asymptotic behavior of weak solutions

We study in this paper the asymptotic analysis of an incompressible Newtonian and non-isothermal problem, when one dimension of the fluid domain tends to zero. We prove the strong convergence of the unknowns which are the temperature, the velocity and the pressure of the fluid, we obtain the limit problem with the specific Reynolds equation, and we also prove the uniqueness of the limit temperature velocity and pressure distributions.