Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

ℙmax variations related to slaloms

We prove the iteration lemmata, which are the key lemmata to show that extensions by Pmax variations satisfy absoluteness for Π₂-statements in the structure 〈H (ω ₂), ∈, NS_{ω 1}, R 〉 for some set R of reals in L (ℝ), for the following statements: (1) The cofinality of the null ideal is ℵ₁. (2) There exists a good basis of the strong measure zero ideal. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Supercompactness and measurable limits of strong cardinals II: Applications to level by level equivalence

We construct models for the level by level equivalence between strong compactness and supercompactness in which for κ the least supercompact cardinal and δ ≤ κ any cardinal which is either a strong cardinal or a measurable limit of strong cardinals, 2^δ > δ ⁺ and δ is < 2^δ supercompact. In these models, the structure of the class of supercompact cardinals can be arbitrary, and the size of the power set of κ can essentially be made as large as desired. This extends and generalizes [5, Theorem 2] and [4, Theorem 4]. We also sketch how our techniques can be used to establish a weak indestructibility result. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Some geometrical aspects of semidefinite linear complementarity problems

This article studies some geometrical aspects of the semidefinite linear complementarity problem (SDLCP), which can be viewed as a generalization of the well-known linear complementarity problem (LCP). SDLCP is a special case of a complementarity problem over a closed convex cone, where the cone considered is the closed convex cone of positive semidefinite matrices. It arises naturally in the unified formulation of a pair of primal-dual semidefinite programming problems. In this article, we introduce the notion of complementary cones in the semidefinite setting using the faces of the cone of positive semidefinite matrices and show that unlike complementary cones induced by an LCP, semidefinite complementary cones need not be closed. However, under R₀-property of the linear transformation, closedness of all the semidefinite complementary cones induced by L is ensured. We also introduce the notion of a principal subtransformation with respect to a face of the cone of positive semidefinite matrices and show that for a self-adjoint linear transformation, strict copositivity is equivalent to strict semimonotonicity of each principal subtransformation. Besides the above, various other solution properties of SDLCP will be interpreted and studied geometrically.     

Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004     

Spectral Representation Theory for Dielectric Behavior of Nonspherical Cell Suspensions

Recent experiments revealed that the dielectric dispersion spectrum of fission yeast cells in a suspension was mainly composed of two sub-dispersions. The low-frequency sub-dispersion depended on the cell length, while the high-frequency one was independent of it. The cell shape effect was simulated by an ellipsoidal cell model but the comparison between theory and experiment was far from being satisfactory. Prompted by the discrepancy, we proposed the use of spectral representation to analyze more realistic cell models. We adopted a shell-spheroidal model to analyze the effects of the cell membrane. It is found that the dielectric property of the cell membrane has only a minor effect on the dispersion magnitude ratio and the characteristic frequency ratio. We further included the effect of rotation of dipole induced by an external electric field, and solved the dipole-rotation spheroidal model in the spectral representation.Good agreement between theory and experiment has been obtained.     

Numerical analysis for the periodic response of nonsmooth dynamic systems of high-dimension

This paper deals with the numerical solution of high dimensional dynamic systems with nonsmooth characteristics, such as structural elastoplasticity, mechanical clearance and dry friction. As a stable periodic response of the system attracts its transient response, it may be extracted from the transient response somehow. The paper suggests a novel iteration scheme to extrapolate the periodic response from a short time history of its transient response by curve fitting. Compared with the current schemes such as shooting and incremental harmonic balancing, the present scheme makes full use of the information within the transient response and the system characters, then its computation efficiency has increased by an order and the numerical convergence depending on the initial iteration has greatly improved.The subject supported by the Chinese Foundation of Aeronautical Science     

Different acute effects of oral and intratracheal administration of disodium arsenate and gallium arsenide on heme synthesis in rats

The acute influences of arsenic compounds on the metabolism of porphyrins and heme in various organs of rats after oral or intratracheal administration of disodium arsenate (Na₂HAsO₄) and gallium arsenide (GaAs) were examined and compared. For the oral administration experiments, 21 or 84 mg of Na₂HAsO₄, or 2 or 4 g of GaAs, per cm³ saline per kg body weight of each animal was administered to Jcl: Wistar male rats and the organs were removed after exsanguination from the vein of the right axilla under anesthesia with ether, 16 h after administration. In the case of intratracheal administration, rats given 8.2 or 16.4 mg of Na₂HAsO₄, or 0.2 or 0.4 g GaAs per cm³ saline per kg body weight were examined under the same experimental conditions as for the administration route. Increase in the body weight of rats was suppressed after intratracheal administration of the two arsenic compounds. In these rats the hematocrit value increased significantly. These changes were not shown by the orally administered rats. Elevation in δ-aminolevulinate synthase (ALA-S, EC 2.3.1.37) activity in erythroblasts by Na₂HAsO₄ was much higher after intratracheal administration than after oral administration. Suppression in the activities of δ-aminolevulinate dehydratase (ALA-D, EC 4.2.1.24) and porphobilinogen deaminase (PBG-D, EC 4.3.1.8) in peripheral erythrocytes by Na₂HAsO₄ and GaAs were stronger by intratracheal administration than by the oral route. Influences of GaAs on the activity of PBG-D in rat liver were shown to be more effective by oral administration than by the intratracheal route. Oral administration of Na₂HAsO₄ and GaAs had a stronger suppression effect on the activities of ALA-D and PBG-D in rat kidney. It seems from these results that the different extents of the influence of arsenic compounds might depend on the routes of intake.     

Physical aging of a polyetherimide: Creep and DSC measurements

Creep and differential scanning calorimetry (DSC) measurements have been used to study the physical aging behavior of a polyetherimide. Isothermal aging temperatures ranged from 160°C to T_g with aging times ranging from 10 min to 8 days. The only measurable effect of physical aging on the short-time creep curves is a shift of the creep compliance to longer times. Andrade plots of the compliance versus the cube root of time are linear at short times with the slope β decreasing with increasing aging time to a constant value once equilibrium is reached. Log β³ is related directly to the degree to which the creep curves shift to longer times with physical aging, and is used in this work as a measure of physical aging. A reduced curve of log β³ versus log aging time is obtained for the aging temperatures investigated by appropriate vertical and horizontal shifts. The enthalpy change during aging increases linearly with the logarithm of the aging time, t_a, leveling off at equilibrium at values which increase with decreasing aging temperature. Hence, both nonequilibrium and equilibrium temperature shift factors can be calculated from the DSC data. Good agreement is observed between the equilibrium temperature shift factors obtained from the creep and DSC data. The temperature dependence of the nonequilibrium temperature shift factors is found to be an order of magnitude smaller than that of the equilibrium shift factors. The time scales to reach equilibrium for enthalpy and for mechanical measurements are found to be the same within experimental error. © 1995 John Wiley & Sons, Inc.     

Synthesis of functional polymers—Vinylidene fluoride based fluorinated copolymers and terpolymers bearing bromoaromatic side‐group

The radical co‐ and terpolymerization of 4‐[(α,β,β‐trifluorovinyl)oxy]bromo benzene (TFVOBB) with 1,1‐difluoroethylene (or vinylidene fluoride, VDF, or VF₂), hexafluoropropene (HFP), perfluoromethyl vinyl ether (PMVE), and chlorotrifluroroethylene (CTFE) is presented. Although TFVOBB could be thermocyclodimerized, it could not homopolymerize under radical initiation. TFVOBB could be copolymerized in solution under a radical initiator with VDF or CTFE comonomers, while its copolymerization with HFP or PMVE were unsuccessful. The terpolymerization of TFVOBB with VDF and HFP, or VDF and PMVE, or VDF and CTFE also led to original fluorinated terpolymers bearing bromoaromatic side‐groups. The conditions of co‐ and terpolymerization were optimized in terms of the nature of the radical initiators, and of the nature of solvents (fluorinated or nonhalogenated). Various monomer concentrations in the co‐ and terpolymers were assessed by ¹⁹F and ¹H‐NMR spectroscopy. The thermal and physico chemical properties were also studied. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 5077–5097, 2004