Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

On c‐Bhaskar Rao designs with block size 4

Several new families of c‐Bhaskar Rao designs with block size 4 are constructed. The necessary conditions for the existence of a c‐BRD (υ,4,λ) are that: (1)λ_min=?λ/3 ≤ c ≤ λ and (2a) c≡λ (mod 2), if υ > 4 or (2b) c≡ λ (mod 4), if υ = 4 or (2c) c≠ λ ? 2, if υ = 5. It is proved that these conditions are necessary, and are sufficient for most pairs of c and λ; in particular, they are sufficient whenever λ?c ≠ 2 for c > 0 and whenever c ? λ_min≠ 2 for c < 0. For c < 0, the necessary conditions are sufficient for υ> 101; for the classic Bhaskar Rao designs, i.e., c = 0, we show the necessary conditions are sufficient with the possible exception of 0‐BRD (υ,4,2)'s for υ≡ 4 (mod 6). © 2002 Wiley Periodicals, Inc. J Combin Designs 10: 361–386, 2002; Published online in Wiley InterScience ( www.interscience.wiley.com ). DOI 10.1002/jcd.10009     

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.     

A high‐resolution scheme for low Mach number flows

A method for computing low Mach number flows using high‐resolution interpolation and difference formulas, within the framework of the Marker and Cell (MAC) scheme, is presented. This increases the range of wavenumbers that are properly resolved on a given grid so that a sufficiently accurate solution can be obtained without extensive grid refinement. Results using this scheme are presented for three problems. The first is the two‐dimensional Taylor–Green flow which has a closed form solution. The second is the evolution of perturbations to constant‐density, plane channel flow for which linear stability solutions are known. The third is the oscillatory instability of a variable density plane jet. In this case, unless the sharp density gradients are resolved, the calculations would breakdown. Under‐resolved calculations gave solutions containing vortices which grew in place rather than being convected out. With the present scheme, regular oscillations of this instability were obtained and vortices were convected out regularly. Stable computations were possible over a wider range of sensitive parameters such as density ratio and co‐flow velocity ratio. Copyright © 2004 John Wiley Sons, Ltd.     

LINEAR ACTIVE STRUCTURES AND MODES (Ⅱ)--DISCRETE SYSTEMS AND BEAMS

The basic concepts about the active structures and some attributes of the modes were presented in paper "Liner Active Structures and Modes ( Ⅰ ) ". The characteristics of the active discrete systems and active beams were discussed, especially, the stability of the active structures and the orthogonality of the eigenvectors. The notes about modes were portrayed by a model of a seven-storeyed building with sensors and actuators. The concept of the adjoint active structure was extended from the discrete systems to the beams that were the representations of the continuous structures. Two types of beams with different placements of the measuring and actuating systems were discussed in detail. One is the beam with the discrete sensors and actuators, and the other is the beam with distributed sensor and actuator function. The orthogonality conditions were derived with the modal shapes of the active beam and its adjoint active beam. An example shows that the variation of eigenvalues with feedback amplitude for the homo-configuration and non-homo-configuration active structures.     

A Note on Geometric Structures of Linear Ordered Orthogonal Arrays and (T,M,S)-Nets of Low Strength

The equivalent geometrical configurations of linear ordered orthogonal arrays are determined when their strengths are 3 and 4. Existence of such geometrical configurations is investigated. They are also useful in the study of (T, M, S)-nets.     

Adaptive remeshing for ductile fracture prediction in metal formingRemaillage adaptatif pour la prévision de la rupture ductile en mise en forme

The analysis of mechanical structures using the Finite Element Method in the framework of large elastoplastic strain, needs frequent remeshing of the deformed domain during computation. Indeed, the remeshing is due to the large geometrical distortion of finite elements and the adaptation to the physical behavior of the solution. This paper gives the necessary steps to remesh a mechanical structure during large elastoplastic deformations with damage. An important part of this process is constituted by geometrical and physical error estimates. The proposed method is integrated in a computational environment using the ABAQUS/Explicit solver and the BL2D-V2 adaptive mesher. To cite this article: H. Borouchaki et al., C. R. Mecanique 330 (2002) 709–716.     

Orders of Indescribable Sets

We extract some properties of Mahlo’s operation and show that some other very natural operations share these properties. The weakly compact sets form a similar hierarchy as the stationary sets. The height of this hierarchy is a large cardinal property connected to saturation properties of the weakly compact ideal.     

Local homology and dimensional full-valuedness

A natural criterion for dimensional full-valuedness of locally compact spaces with finitely generated local homology is given.