Lattice models of hydrogen-bonded solvents. II

We consider a lattice model of a binary mixture in which each molecule of one component can form zero, one, or two bonds to molecules of the same species on neighboring vertices of the square lattice. We allow the energy of molecules with two bonds to depend on the valence angle, thus generalizing the first paper in this series. If straight polymeric configurations are favored over all others, then a phase transition occurs for low enough temperature. On the other hand, if bent configurations are favored, there is no phase transition. Analogous results are obtained for the hexagonal lattice, where we distinguish energetically betweencis andtrans isomerism of four bonded molecules.Supported by the Danish Science Foundation under grant 511-3635.     

Inexact Newton methods for the nonlinear complementarity problem

An exact Newton method for solving a nonlinear complementarity problem consists of solving a sequence of linear complementarity subproblems. For problems of large size, solving the subproblems exactly can be very expensive. In this paper we study inexact Newton methods for solving the nonlinear, complementarity problem. In such an inexact method, the subproblems are solved only up to a certain degree of accuracy. The necessary accuracies that are needed to preserve the nice features of the exact Newton method are established and analyzed. We also discuss some extensions as well as an application. This research was based on work supported by the National Science Foundation under grant ECS-8407240.     

A note on optimality conditions for the Euclidean. Multifacility location problem

Thekey problem of the Euclidean multifacility location (EMFL) problem is to decide whether a givendead point is optimal. If it is not optimal, we wish to compute a descent direction. This paper extends the optimality conditions of Calamai and Conn and Overton to the case when the rows of the active constraints matrix are linearly dependent. We show that linear dependence occurs wheneverG, the graph of the coinciding facilities, has a cycle. In this case the key problem is formulated as a linear least squares problem with bounds on the Euclidean norms of certain subvectors.     

Application of the sequential gradient-restoration algorithm to thermal convective instability problems

The problem of the thermal stability of a horizontal incompressible fluid layer with linear and nonlinear temperature distributions is solved by using the sequential gradient-restoration algorithm developed for optimal control problems. The hydrodynamic boundary conditions for the layer include a rigid or free upper surface and a rigid lower surface. The resulting disturbing equations are solved as a Bolza problem in the calculus of variations. The results of the study are compared with the existing works in the literature.The authors acknowledge valuable discussions with Dr. A. Miele.     

一个肿瘤生长自由边界问题解的渐近性态

该文研究根据Byrne和Chaplain的思想建立的一个描述抑制物作用下无坏死核肿瘤生长的数学模型, 这个模型是一个非线性反应扩散方程组的自由边界问题. 作者运用反应扩散方程理论中的上下解方法结合自由边界问题的迭代技巧, 研究了解的渐近性态, 在营养物消耗函数f、抑制物消耗函数g和肿瘤细胞繁衍函数S的一些一般条件下,证明当常数c₁,c₂(肿瘤细胞分裂速率和营养物、抑制物扩散速率的比值)都非常小时,在一定的初边值条件下肿瘤趋于消失,在另外一些初边值条件下肿瘤半径趋于一个常数,进而时变解将趋于一个稳态解.     

A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.     

The contact Yamabe flow on K-contact manifolds

We use the contact Yamabe flow to find solutions of the contact Yamabe problem on K-contact manifolds.     

Parallel Schwarz alternating methods for anisotropic diffusion of speckled medical images

This paper deals with ultrasound medical image processing, particularly to filter the noise while preserving the edges and structures of information. The mathematical processing consists in solving by a numerical way a nonlinear evolutive boundary value problem. Several numerical semi-implicit time marching schemes are considered and analyzed. At each time step, parallel synchronous or asynchronous Schwarz alternating methods are used to solve the linear system and its convergence is studied. Lastly, the results of sequential and parallel simulations are presented.