Homogenization of two-phase flow: high contrast of phase permeabilityHomogénéisation d'écoulement diphasique : grand contraste de perméabilité d'une phase

The steady-state two-phase flow non-linear equation is considered in the case when one of phases has low effective permeability in some periodic set, while on the complementary set it is high; the second phase has no contrast of permeabilities in different zones. A homogenization procedure gives the homogenized model with macroscopic effective permeability of the second phase depending on the gradient and on the second order derivatives of the macroscopic pressure of the first phase. This effect cannot be obtained by classical (one small parameter) homogenization. To cite this article: G.P. Panasenko, G. Virnovsky, C. R. Mecanique 331 (2003).     

Bohr-Mottelson转动谱公式参数之间的新关系式

在利用Harris两参数公式研究Bohr-Mottelson转动谱公式参数之间的关系的基础上,改用Harris三参数公式,并由此提出了Bohr-Mottelson转动谱公式参数之间的新关系式,进而用I(I+1)四参数展开式计算了A～60,80,130,140,150,190区超形变偶偶核的基带和锕系和稀土区正常形变核基带,讨论了参数之间的关系,发现新关系式与实验较好地符合.     

C++QED: an object-oriented framework for wave-function simulations of cavity QED systems

We present a framework for efficiently performing Monte Carlo wave-function simulations in cavity QED with moving particles. It relies heavily on the object-oriented programming paradigm as realised in C++, and is extensible and applicable for simulating open interacting qua ntum dynamics in general. The user is provided with a number of “elements”, e.g. pumped moving particles, pumped lossy cavity modes, and various interactions to compose complex interacting systems, which contain several particles moving in electromagnetic fields of various configurations, and perform wave-function simulations on such systems. A number of tools are provided to facilitate the implementation of new elements.     

On the Constituent Structure of a Classical Mixture

Thiw work is concerned with the formulation of constituent interactions and corresponding balance relations in classical mixture theory as based on a model for the (classical) constituent structure of such a mixture.     

Protein Conformation of a Lattice Model Using Tabu Search

We apply tabu search techniques to the problem of determining the optimal configuration of a chain of protein sequences on a cubic lattice. The problem under study is difficult to solve because of the large number of possible conformations and enormous amount of computations required. Tabu search is an iterative heuristic procedure which has been shown to be a remarkably effective method for solving combinatorial optimization problems. In this paper, an algorithm is designed for the cubic lattice model using tabu search. The algorithm has been tested on a chain of 27 monomers. Computational results show that our method outperforms previously reported approaches for the same model.     

Thermal helium clusters at 3.2 Kelvin in classical and semiclassical simulations

The thermodynamic stability of⁴He_4–13 at 3.2 K is investigated with the classical Monte Carlo method, with the semiclassical path-integral Monte Carlo (PIMC) method, and with the semiclassical all-order many-body method. In the all-order many-body simulation the dipole-dipole approximation including short-range correction is used. The resulting stability plots are discussed and related to recent TOF experiments by Stephens and King. It is found that with classical Monte Carlo of course the characteristics of the measured mass spectrum cannot be resolved. With PIMC, switching on more and more quantum mechanics. by raising the number of virtual time steps results in more structure in the stability plot, but this did not lead to sufficient agreement with the TOF experiment. Only the all-order many-body method resolved the characteristic structures of the measured mass spectrum, including magic numbers. The result shows the influence of quantum statistics and quantum mechanics on the stability of small neutral helium clusters.     

Heat transfer in lattice BGK modeled fluid

The thermal lattice BGK model is a recently suggested numerical tool aiming at solving problems of thermohydrodynamics. The quality of the lattice BGK simulation is checked in this paper by calculating temperature profiles in the Couette flow under different Eckert and Mach numbers. A revised lower order model is proposed to improve the accuracy and the higher order model is proved to be advantageous in this respect, especially in the flow regime with a higher Mach number.     

Metastability of a circular o-ring due to intrinsic curvature

An o-ring takes spontaneously the shape of a chair when strong enough torsion is applied in its tangent plane. This state is metastable, since work has to be done on the o-ring to return to the circular shape. We show that this metastable state exists in a Hamiltonian where curvature and torsion are coupled via an intrinsic curvature term. If the o-ring is constrained to be planar (2d case), this metastable state displays a kink-anti-kink pair. This state is metastable if the ratio is less than , where C and A are the torsion and the bending elastic constants [#!landau!#]. In three dimensions, our variational approach shows that . This model can be generalized to the case where the bend is induced by a concentration field which follows the variations of the curvature. Received: 27 August 1997 / Revised: 23 October 1997 / Accepted: 12 November 1997     

Numerical solution of eigenvalue problems for Hamiltonian systems

This paper discusses the numerical solution of eigenvalue problems for Hamiltonian systems of ordinary differential equations. Two new codes are presented which incorporate the algorithms described here; to the best of the author’s knowledge, these are the first codes capable of solving numerically such general eigenvalue problems. One of these implements a new new method of solving a differential equation whose solution is a unitary matrix. Both codes are fully documented and are written inPfort-verifiedFortran 77, and will be available in netlib/aicm/sl11f and netlib/aicm/sl12f.     

Some properties of the bilevel programming problem

The purpose of this paper is to elaborate on the difficulties accompanying the development of efficient algorithms for solving the bilevel programming problem (BLPP). We begin with a pair of examples showing that, even under the best of circumstances, solutions may not exist. This is followed by a proof that the BLPP is NP-hard.This work was partially supported by a grant from the Advanced Research Program of the Texas Higher Education Coordinating Board.