Mass‐lumping or not mass‐lumping for eigenvalue problems

In this article we analyze the effect of mass‐lumping in the linear triangular finite element approximation of second‐order elliptic eigenvalue problems. We prove that the eigenvalue obtained by using mass‐lumping is always below the one obtained with exact integration. For singular eigenfunctions, as those arising in non convex polygons, we prove that the eigenvalue obtained with mass‐lumping is above the exact eigenvalue when the mesh size is small enough. So, we conclude that the use of mass‐lumping is convenient in the singular case. When the eigenfunction is smooth several numerical experiments suggest that the eigenvalue computed with mass‐lumping is below the exact one if the mesh is not too coarse. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 653–664, 2003     

The Computational Complexity of the Chow Form

We present a bounded probability algorithm for the computation of the Chowforms of the equidimensional components of an algebraic variety. In particular, this gives an alternative procedure for the effective equidimensional decomposition of the variety, since each equidimensional component is characterized by its Chow form. The expected complexity of the algorithm is polynomial in the size and the geometric degree of the input equation system defining the variety. Hence it improves (or meets in some special cases) the complexity of all previous algorithms for computing Chow forms. In addition to this, we clarify the probability and uniformity aspects, which constitutes a further contribution of the paper. The algorithm is based on elimination theory techniques, in line with the geometric resolution algorithm due to M. Giusti, J. Heintz, L. M. Pardo, and their collaborators. In fact, ours can be considered as an extension of their algorithm for zero-dimensional systems to the case of positive-dimensional varieties. The key element for dealing with positive-dimensional varieties is a new Poisson-type product formula. This formula allows us to compute the Chow form of an equidimensional variety from a suitable zero-dimensional fiber. As an application, we obtain an algorithm to compute a subclass of sparse resultants, whose complexity is polynomial in the dimension and the volume of the input set of exponents. As another application, we derive an algorithm for the computation of the (unique) solution of a generic overdetermined polynomial equation system.     

A technique to construct symmetric variable-stepsize linear multistep methods for second-order systems

Some previous works show that symmetric fixed- and variable-stepsize linear multistep methods for second-order systems which do not have any parasitic root in their first characteristic polynomial give rise to a slow error growth with time when integrating reversible systems. In this paper, we give a technique to construct variable-stepsize symmetric methods from their fixed-stepsize counterparts, in such a way that the former have the same order as the latter. The order and symmetry of the integrators obtained is proved independently of the order of the underlying fixed-stepsize integrators. As this technique looks for efficiency, we concentrate on explicit linear multistep methods, which just make one function evaluation per step, and we offer some numerical comparisons with other one-step adaptive methods which also show a good long-term behaviour.

     

Measurement of the dynamic elastic constants of short isotropic cylinders

A method based on a single test is proposed to characterize the elasticity of an isotropic homogeneous material in the shape of a cylinder of any slenderness (length-diameter) ratio. Firstly, the Rayleigh-Ritz method is used to determine the natural frequencies of the cylinders vibrating axisymmetrically. The study is focused on cylindrical samples with diameter and length of similar magnitude so that the shear modulus and the Poisson ratio can be calculated simultaneously. Subsequently, the theoretical results for cylinders of slenderness ratio between 0.1 and 3 are analyzed in order to obtain the data required to determine the elastic constants from one of the two lowest measured natural frequencies and their quotient. The analysis of the results demonstrates that any slenderness ratio is useful in the calculation of the elastic constants, although in some cases the third natural frequency should be used. Furthermore, the influence of the length-diameter quotient on the sensitivity of the method is analyzed by evaluating the systematic uncertainties for both dynamic elastic constants. Finally, the method is experimentally tested by characterizing two steel cylinders with slenderness ratios 0.1 and 1, respectively. The results demonstrate that uncertainties for both Poisson ratio and the shear modulus are smaller when the slenderness ratio is 1.     

Asymptotically balanced schemes for non-homogeneous hyperbolic systems – application to the Shallow Water equations

In this work we introduce a class of balanced numerical schemes, up to second order, for the solution of general non-homogeneous hyperbolic systems of conservation laws. We give a general technique to build such schemes. We also prove that they balance up to second order a large class of steady solutions in the whole domain but some subset whose measure tends to zero as the grid size decreases to zero. We finally present an application to Shallow Water equations that exhibit the good performances of some of the schemes introduced. To cite this article: T. Chacón Rebollo et al., C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

Lyapunov-type inequalities for nonlinear systems

In this paper, by using elementary analysis, we establish some new Lyapunov-type inequalities for nonlinear systems of differential equations, special cases of which contain the well-known equations such as Emden-Fowler-type and half-linear equations. The inequalities obtained here can be used as handy tools in the study of qualitative behaviour of solutions of the associated equations.     

Hydrophilic and hydrophobic copolymer systems based on acrylic derivatives of pyrrolidone and pyrrolidine

This article deals with the synthesis of hydrophilic methacrylic monomers derived from ethyl pyrrolidone [2‐ethyl‐(2‐pyrrolidone) methacrylate (EPM)] and ethyl pyrrolidine [2‐ethyl‐(2‐pyrrolidine) methacrylate (EPyM)] and their respective homopolymers. For the determination of their reactivity in radical copolymerization reactions, both monomers were copolymerized with methyl methacrylate (MMA), the reactivity ratios being calculated by the application of linear and nonlinear mathematical methods. EPM and MMA had ratios of r_EPM = 1.11 and r_MMA = 0.76, and this indicated that EPM with MMA had a higher reactivity in radical copolymerization processes than vinyl pyrrolidone (VP; r_VP = 0.005 and r_MMA = 4.7). EPyM and MMA had reactivity ratios of r_EPyM = 1.31 and r_MMA = 0.92, and this implied, as for the EPM–MMA copolymers, a tendency to form random or Bernoullian copolymers. The glass‐transition temperatures of the prepared copolymers were determined by differential scanning calorimetry (DSC) and were found to adjust to the Fox equation. Total‐conversion copolymers were prepared, and their behavior in aqueous media was found to be dependent on the copolymer composition. The swelling kinetics of the copolymers followed water transport mechanism case II, which is the most desirable kinetic behavior for a swelling controlled‐release material. Finally, the different states of water in the hydrogels—nonfreezing water, freezing bound water, and unbound freezing water—were determined by DSC and found to be dependent on the hydrophilic and hydrophobic units of the copolymers. © 2002 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 395–407, 2003     

Derivation of the k–ε model for locally homogeneous turbulence by homogenization techniques

We derive the incompressible and compressible k–ε model for locally homogeneous turbulence. The model is rigorously derived on formal mathematical grounds using the MPP modelling technique. This lets us calculate by either analytical or numerical means the closure constants of the model. To cite this article: T. Chacón Rebollo, D. Franco Coronil, C. R. Acad. Sci. Paris, Ser. I 337 (2003).     

Spinors, Inflation, and Non-Singular Cyclic Cosmologies

We consider toy cosmological models in which a classical, homogeneous, spinor field provides a dominant or sub-dominant contribution to the energy-momentum tensor of a flat Friedmann-Robertson-Walker universe. We find that, if such a field were to exist, appropriate choices of the spinor self-interaction would generate a rich variety of behaviors, quite different from their widely studied scalar field counterparts. We first discuss solutions that incorporate a stage of cosmic inflation and estimate the primordial spectrum of density perturbations seeded during such a stage. Inflation driven by a spinor field turns out to be unappealing as it leads to a blue spectrum of perturbations and requires considerable fine-tuning of parameters. We next find that, for simple, quartic spinor self-interactions, non-singular cyclic cosmologies exist with reasonable parameter choices. These solutions might eventually be incorporated into a successful past- and future-eternal cosmological model free of singularities. In an Appendix, we discuss the classical treatment of spinors and argue that certain quantum systems might be approximated in terms of such fields.     

Interpretation of the polyelectrolyte and antipolyelectrolyte effects of poly(N‐vinylimidazole‐co‐sodium styrenesulfonate) hydrogels

Hydrogels of N‐vinylimidazole (VI) and sodium styrenesulfonate (SSS) were synthesized in aqueous solution by radical crosslinking copolymerization with N,N′‐methylene‐bis(acrylamide) as crosslinker. Swelling in several saline solutions was measured for hydrogel samples synthesized with different comonomer concentrations (C_T = 10, 25, or 40%) and with SSS mole fractions covering a broad range (f_SSS = 0–0.7), while the crosslinker ratio was 2 wt % in all cases. The degree of swelling in aqueous solution with a specific ionic strength (μ), plotted versus the SSS composition of the feed, shows a minimum for any set of samples synthesized with a fixed C_T. The dependence of swelling on μ shows both polyelectrolyte (f_SSS beyond the minimum) and antipolyelectrolyte behaviors (in the low f_SSS limit). It was found that the nonGaussian factor of the crosslinking density and the polymer‐solvent interaction parameter increase with f_SSS for any C_T. Moreover, in the low f_SSS limit, the osmotic swelling pressure is governed not only by the ionic contribution, but also by the polymer‐solvent mixing and, the concentration of mobile counterions inside the gel is not proportional to the net fixed charge but to the addition of cationic and anionic side groups, what discards the formation of ionic pairs. The antipolyelectrolyte effect is interpreted as due to the increasing protonation of VI as μ goes up. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 1683–1693, 2007