Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

Ergodic properties and thermodynamic behavior of elementary reversible cellular automata. I. Basic properties

This is the first part of a series devoted to the study of thermodynamic behavior of large dynamical systems with the use of a family of fully-discrete and conservative models named elementary reversible cellular automata (ERCAs). In this paper, basic properties such as conservation laws and phase space structure are investigated in preparation for the later studies. ERCAs are a family of one-dimensional reversible cellular automata having two Boolean variables on each site. Reflection and Boolean conjugation symmetries divide them into 88 equivalence classes. For each rule, additive conserved quantities written in a certain form are regarded as a kind of energy, if they exist. By the aid of the discreteness of the variables, every ERCA satisfies the Liouville theorem or the preservation of phase space volume. Thus, if an energy exists in the above sense, statistical mechanics of the model can formally be constructed. If a locally defined quantity is conserved, however, it prevents the realization of statistical mechanics. The existence of such a quantity is examined for each class and a number of rules which have at least one energy but no local conservation laws are selected as hopeful candidates for the realization of thermodynamic behavior. In addition, the phase space structure of ERCAs is analyzed by enumerating cycles exactly in the phase space for systems of comparatively small sizes. As a result, it is revealed that a finite ERCA is not ergodic, that is, a large number of orbits coexist on an energy surface. It is argued that this fact does not necessarily mean the failure of thermodynamic behavior on the basis of an analogy with the ergodic nature of infinite systems.     

The Design for Automation of Pellet-Forming Process in HL-2A Pellet Fueling System

In the ever fusion experiments in SWIP, pellet forming process was carried out through adjusting relative devices by staff member in site, which will make every pellet-forming process slight distinction and will result in pellet difference in shape, size and intensity. In the intervals of HL-2A discharges, staff member have to go site to accomplish the pellet-forming process, this wastes human power and increase the potential danger. So it is necessary to develop a remote control system to perform the pellet-forming process. The control system needs have the features of real-time, reliability and be easy to operate and maintain.     

ISLAND—Inverse Square Law Acceleration Measurement Using iNertial Drift

By using a novel free-flying, rotating, optical bench aboard the International Space Station (ISS) as the basis for a measurement of differential acceleration between two gravitating bodies, in principle the Newtonian inverse-square law and the constant of gravitation, G, can be determined at the parts in 10⁶ level.     

广义Sierpinski集

本文考察了包括平面上的各种广义 Cantor集 ,Sierpinski集和包括某些连续不可微曲线在内的广义 Sierpinski集 .由相似变换 ,导出了它们的级数表达式 ,并利用它和字符串空间的对应关系 ,计算出它们的Hausdorff维数     

A WKB analysis of the buckling condition for a cylindrical shell of arbitrary thickness subjected to an external pressure

We consider the plane-strain buckling of a cylindrical shellof arbitrary thickness which is made of a Varga material andis subjected to an external hydrostatic pressure on its outersurface. The WKB method is used to solve the eigenvalue problemthat results from the linear bifurcation analysis. We show thatthe circular cross-section buckles into a non-circular shapeat a value of µ₁ which depends on A₁/A₂ and a mode number,where A₁ and A₂ are the undeformed inner and outer radii, andµ₁ is the ratio of the deformed inner radius to A₁. Inthe large mode number limit, we find that the dependence ofµ₁ on A₁/A₂ has a boundary layer structure: it is constantover almost the entire region of 0 < A₁/A₂ < 1 and decreasessharply from this constant value to unity as A₁/A₂ tends tounity. Our asymptotic results for A₁ – 1 = O(1) and A₁– 1 = O(1/n) are shown to agree with the numerical resultsobtained by using the compound matrix method.     

Sobolev-Hardy不等式与临界双重调和问题

该文讨论一类带有奇异系数的双重调和方程{△^2u-μu/|x|^s=f(x,u),x∈Ω,u=δu/δv=0,x∈δΩ，这里Ω包含R^N是包含0的有界光滑区域，u∈H0^2(Ω），μ∈R是参数，0≤s≤2,△^2=△△表示双重拉普拉斯算子，当f(x,u)=u^p,p=2N/N-4时，上述问题就是一个临界双重调和问题，该文运用Sobolev-Hardy不等式和变分方法，得到它的解的存在性的一些结果。     

Absorbing boundary conditions for reaction-diffusion equations

We treat here of the question of absorbing boundary conditionsfor nonlinear diffusion equations. We use the conditions designedfor the linear equation, we prove them to be well posed forthe nonlinear problem, and through numerical experiments thatthey are well suited for reaction–diffusion equations.     

论恒力作用下质点的相对论运动

在相对论情况下，讨论了受恒力作用的质点的运动规律。     

Group Representation of Relativistic Quantum Crystal Optics

A group representation of radiation propagation in an anisotropic medium is developed. The system of wave equations for electromagnetic potentials, obtained from the Maxwell equations with account for the constitutive equations, has been factorized. It is shown that the linear differential operator of the factorized system is orthogonal in transparent crystals and unitary in gyrotropic ones and is represented through the momentum operator. On the basis of the commutation relations for the components of this operator, the eigenvalue problem has been solved and the expression for the change in the radiant energy in the crystal in the form of spherical waves has been obtained. The dependences of the ray and phase velocities and the polarization vectors of waves on the birefringence anisotropy and gyrotropy as well as on the angular momentum, displacement current, and bound charge determining them have been analyzed. It has been established that in the general case of gyrotropic crystals where the nonreciprocity phenomenon takes place and in magnetoelectrics Maxwell equations are represented in a form similar to the Dirac equations and the electromagnetic radiation is correctly described by means of bispinors and is quantized as fermions.