Numerical solution of the Monge–Ampère equation by a Newton's algorithm

We solve numerically the Monge–Ampère equation with periodic boundary condition using a Newton's algorithm. We prove convergence of the algorithm, and present some numerical examples, for which a good approximation is obtained in 10 iterations. To cite this article: G. Loeper, F. Rapetti, C. R. Acad. Sci. Paris, Ser. I 340 (2005).     

Convergence of multistep discretizations of DAEs

Standard ODE methods such as linear multistep methods encounter difficulties when applied to differential-algebraic equations (DAEs) of index greater than 1. In particular, previous results for index 2 DAEs have practically ruled out the use of all explicit methods and of implicit multistep methods other than backward difference formulas (BDFs) because of stability considerations. In this paper we embed known results for semi-explicit index 1 and 2 DAEs in a more comprehensive theory based on compound multistep and one-leg discretizations. This explains and characterizes the necessary requirements that a method must fulfill in order to be applicable to semi-explicit DAEs. Thus we conclude that the most useful discretizations are those that avoid discretization of the constraint. A freer use of e.g. explicit methods for the non-stiff differential part of the DAE is then possible.Dedicated to Germund Dahlquist on the occasion of his 70th birthdayThis author thanks the Centro de Estadística y Software Matemático de la Universidad Simón Bolivar (CESMa) for permitting her free use of its research facilities.Partial support by the Swedish Research Council for Engineering Sciences TFR under contract no. 222/91-405.     

Mathematical and numerical analysis of a simplified time-dependent viscoelastic flow

A time-dependent model corresponding to an Oldroyd-B viscoelastic fluid is considered, the convective terms being disregarded. Global existence in time is proved in Banach spaces provided the data are small enough, using the implicit function theorem and a maximum regularity property for a three fields Stokes problem. A finite element discretization in space is then proposed. Existence of the numerical solution is proved for small data, so as a priori error estimates, using again an implicit function theorem. Supported by the Swiss National Science Foundation. Fellowship PBEL2–114311.     

An alternative approach for non-linear optimal control problems based on the method of moments

We propose an alternative method for computing effectively the solution of non-linear, fixed-terminal-time, optimal control problems when they are given in Lagrange, Bolza or Mayer forms. This method works well when the nonlinearities in the control variable can be expressed as polynomials. The essential of this proposal is the transformation of a non-linear, non-convex optimal control problem into an equivalent optimal control problem with linear and convex structure. The method is based on global optimization of polynomials by the method of moments. With this method we can determine either the existence or lacking of minimizers. In addition, we can calculate generalized solutions when the original problem lacks of minimizers. We also present the numerical schemes to solve several examples arising in science and technology.     

Anodic stripping coulometry with collection for absolute trace analysis of lead, using flow-through electrochemical cells with porous working electrodes

Summary Flow-through electrochemical cells with porous working electrodes made of crushed reticulated vitreous carbon and plated with mercury were used for absolute analysis of trace amounts of lead by anodic stripping coulometry with collection (ASCWC) in a flow system. The role of mercury coating, flow rate and pH were investigated. The coulombic content of the collection peak corresponded to the theoretical values calculated by Faraday's law in a concentration range from about 10^–9 to 10^–6 mol/l. The relative error and the relative standard deviation was +0.15% and 0.8%, respectively for 2×10^–6 mol/l analyte concentration. The absolute detection limit (3 s) was 0.1 ng of Pb, the linear response range 7×10⁴. One leave from: Department of Analytical Chemistry, Slovak Technical University, CS-812 37 Bratislava, Czechoslovakia     

Left, right, and inner socles of associative systems

We investigate the basic properties of the different socles that can be considered in not necessarily semiprime associative systems. Among other things, we show that the socle defined as the sum of minimal (or minimal and trivial) inner ideals is always an ideal. When trivial inner ideals are included, this inner socle contains the socles defined in terms of minimal left or right ideals. This revised version was published online in June 2006 with corrections to the Cover Date.     

Speed of adaptation in structured populations

Adaptation of populations takes place with the occurrence and subsequent fixation of mutations that confer some selective advantage to the individuals which acquire it. For this reason, the study of the process of fixation of advantageous mutations has a long history in the population genetics literature. Particularly, the previous investigations aimed to find out the main evolutionary forces affecting the strength of natural selection in the populations. In the current work, we investigate the dynamics of fixation of beneficial mutations in a subdivided population. The subpopulations (demes) can exchange migrants among their neighbors, in a migration network which is assumed to have either a random graph or a scale-free topology. We have observed that the migration rate drastically affects the dynamics of mutation fixation, despite of the fact that the probability of fixation is invariant on the migration rate, accordingly to Maruyama's conjecture. In addition, we have noticed a topological dependence of the adaptive evolution of the population when clonal interference becomes effective.     

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     

Diagramme polythermique du systeme ternaire H2O—Al(NO3)3—Mg(NO3)2

The solid—liquid equilibria of the ternary system H₂O—Al(NO₃)₃—Mg(NO₃)₂ were studied at –30, –20, –10 and 0°C by using a synthetic method which allows to detemine all the characteristic points of isothermal sections. The stable solid phases which appear are respectively: ice, Al(NO₃)₃·9H₂O, Mg(NO₃)₂·9H₂O and Mg(NO₃)₂·6H₂O. Neither double salts nor mixed crystals are observed in the temperature and composition field studied. Polytherm diagram layout show two invariant transformations correspond with an eutectic point and a peritectic point.This revised version was published online in November 2005 with corrections to the Cover Date.     

Multipole-Bound Molecular Negative Ions

 Within the framework of a simple electrostatic model, as compared to recent experimental results, we here discuss the stability of very weakly bound molecular negative ions. In contrast with the case of conventional valence anions, the excess electron is then located in a very diffuse orbital and is mainly bound by electrostatic dipolar, quadrupolar, and polarization forces, at large distances from the neutral molecular core. By fitting a single repulsion parameter of the model to the available experimental data, it is possible to make quantitative predictions of the excess-electron binding energies in these species. Critical values of the dipole moment, quadrupole moment or polarizability required for the observation of stable multipole-bound negative ions are predicted and compared to available experimental data and ab initio calculations. Received October 24, 2001; accepted for publication November 16, 2001