A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

A discontinuous-Galerkin-based immersed boundary method

A numerical method to approximate partial differential equations on meshes that do not conform to the domain boundaries is introduced. The proposed method is conceptually simple and free of user-defined parameters. Starting with a conforming finite element mesh, the key ingredient is to switch those elements intersected by the Dirichlet boundary to a discontinuous-Galerkin approximation and impose the Dirichlet boundary conditions strongly. By virtue of relaxing the continuity constraint at those elements. boundary locking is avoided and optimal-order convergence is achieved. This is shown through numerical experiments in reaction-diffusion problems. Copyright (c) 2008 John Wiley & Sons, Ltd.

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 Elsevier B.V. All rights reserved.

Accurate velocity reconstruction for Discontinuous Galerkin approximations of two-phase porous media flows

We consider Discontinuous Galerkin approximations of two-phase, immiscible porous media flows in the global pressure/fractional flow formulation with capillary pressure. A sequential approach is used with a backward Euler step for the saturation equation, equal-order interpolation for the pressure and the saturation, and without any limiters. An accurate total velocity field is recovered from the global pressure equation to be used in the saturation equation. Numerical experiments show the advantages of the proposed reconstruction. To cite this article: A. Ern et al., C R. Acad. Sci. Paris, Ser. 1347 (2009). (C) 2009 Academie des sciences. Published by Elsevier Masson SAS. All rights reserved.

Discontinuous local semiflows for Kurzweil equations leading to LaSalle`s invariance principle for differential systems with impulses at variable times

In this paper, we consider an initial value problem for a class of generalized ODEs, also known as Kurzweil equations, and we prove the existence of a local semidynamical system there. Under certain perturbation conditions, we also show that this class of generalized ODEs admits a discontinuous semiflow which we shall refer to as an impulsive semidynamical system. As a consequence, we obtain LaSalle`s invariance principle for such a class of generalized ODEs. Due to the importance of LaSalle`s invariance principle in studying stability of differential systems, we include an application to autonomous ordinary differential systems with impulse action at variable times. (C) 2011 Elsevier Inc. All rights reserved.

An improved finite element space for discontinuous pressures

We consider incompressible Stokes flow with an internal interface at which the pressure is discontinuous, as happens for example in problems involving surface tension. We assume that the mesh does not follow the interface, which makes classical interpolation spaces to yield suboptimal convergence rates (typically, the interpolation error in the L(2)(Omega)-norm is of order h(1/2)). We propose a modification of the P(1)-conforming space that accommodates discontinuities at the interface without introducing additional degrees of freedom or modifying the sparsity pattern of the linear system. The unknowns are the pressure values at the vertices of the mesh and the basis functions are computed locally at each element, so that the implementation of the proposed space into existing codes is straightforward. With this modification, numerical tests show that the interpolation order improves to O(h(3/2)). The new pressure space is implemented for the stable P(1)(+)/P(1) mini-element discretization, and for the stabilized equal-order P(1)/P(1) discretization. Assessment is carried out for Poiseuille flow with a forcing surface and for a static bubble. In all cases the proposed pressure space leads to improved convergence orders and to more accurate results than the standard P(1) space. In addition, two Navier-Stokes simulations with moving interfaces (Rayleigh-Taylor instability and merging bubbles) are reported to show that the proposed space is robust enough to carry out realistic simulations. (c) 2009 Elsevier B.V. All rights reserved.

Critical discontinuous phase transition in the threshold contact process

We analyze a threshold contact process on a square lattice in which particles are created on empty sites with at least two neighboring particles and are annihilated spontaneously. We show by means of Monte Carlo simulations that the process undergoes a discontinuous phase transition at a definite value of the annihilation parameter, in accordance with the Gibbs phase rule, and that the discontinuous transition exhibits critical behavior. The simulations were performed by using boundary conditions in which the sites of the border of the lattice are permanently occupied by particles.

Computational aspects of harmonic wavelet Galerkin methods and an application to a precipitation front propagation model

This article is dedicated to harmonic wavelet Galerkin methods for the solution of partial differential equations. Several variants of the method are proposed and analyzed, using the Burgers equation as a test model. The computational complexity can be reduced when the localization properties of the wavelets and restricted interactions between different scales are exploited. The resulting variants of the method have computational complexities ranging from O(N(3)) to O(N) (N being the space dimension) per time step. A pseudo-spectral wavelet scheme is also described and compared to the methods based on connection coefficients. The harmonic wavelet Galerkin scheme is applied to a nonlinear model for the propagation of precipitation fronts, with the front locations being exposed in the sizes of the localized wavelet coefficients. (C) 2011 Elsevier Ltd. All rights reserved.

Alternative mating tactics in dimorphic males of the harvestman Longiperna concolor (Arachnida: Opiliones)

Intense male-male competition for females may drive the evolution of male morphological dimorphism, which is frequently associated with alternative mating tactics. Using modern techniques for the detection of discontinuous allometries, we describe male dimorphism in the Neotropical harvestman Longiperna concolor, the males of which use their elongated, sexually dimorphic legs IV in fights for the possession of territories where females lay eggs. We also tested three predictions related to the existence of alternative mating tactics: (1) if individuals with relatively longer legs IV (majors) are more likely to monopolize access to reproductive resources, they are expected to remain close to stable groups of females more than individuals with relatively shorter legs IV (minors) do; (2) if minors achieve fertilization by moving between territories, they are expected to be less faithful to specific sites; and (3) majors should be observed in aggressive interactions more often. We individually marked all the individuals from a population of Longiperna during the reproductive season and recorded the location of each sighting for males and females as well as the identity of males involved in fights. Majors were more likely to have harems, and large majors were even more likely to do so. Majors were more philopatric and all males involved in fights belonged to this morph. These results strongly suggest that the mating tactic of the majors is based on resource defense whereas that of the minors probably relies on sneaking into the territories of the majors and furtively copulating with females.

Separation and viability of gill and hepatopancreatic cells of a mangrove crab Ucides cordatus

The gills contain essential cells for respiration and osmoregulation, whereas the hepatopancreas is the site of digestion, absorption, and nutrients storage. The aim of this work was to separate and characterize gill and hepatopancreatic cells of the mangrove crab, Ucides cordatus. For gills, the methodology consisted of an enzymatic cellular dissociation using Trypsin at 0.5%, observation of cellular viability with Tripan Blue, and separation of cells using discontinuous sucrose gradient at concentrations of 10%, 20%, 30%, and 40%. The hepatopancreatic cells were dissociated by magnetic stirring, with posterior separation by sucrose gradient at the same concentrations above. For gills, a high cellular viability was observed (92.5 +/- 2.1%), with hemocyte cells in 10% sucrose layer (57.99 +/- 0.17%, *P < 0.05), principal cells in the 20% sucrose layer (57.33 +/- 0.18, *P < 0.05), and thick cells and pillar cells in the 30% and 40% sucrose layers, respectively (39.54 +/- 0.05%, *P < 0.05; and 41.81 +/- 0.04%, *P < 0.05). The hepatopancreatic cells also showed good viability (79.22 +/- 0.02%), with the observation of embryonic (E) cells in the 10% sucrose layer (67.87 +/- 0.06%, **P < 0.001), resorptive (R) and fibrillar (F) cells in the 20% and 30% sucrose layers (44.71 +/- 0.06%, **P < 0.001, and 43.25 +/- 0.01%, *P < 0.05; respectively), and blister (B) cells in the 40% sucrose layer (63.09 +/- 0.03%, **P < 0.001). The results are a starting point for in vitro studies of heavy metal transport in isolated cells of the mangrove crab U. cordatus, subjected to contamination by metals in the mangrove habitat where they are found.

Mixture of changing uniaxial micellar forms in lyotropic biaxial nematics

Lyotropic nematics consisting of surfactant-cosurfactant water solutions may present a biaxial phase or direct U(+) <-> U(-) transitions, in different regions of the temperature-relative concentration phase diagram, for different systems and compositions. We propose that these may be related to changes of uniaxial micellar form, which may occur either smoothly or abruptly. Smooth change of cylinder-like into disc-like shapes requires a distribution of Maier-Saupe interaction constants and we consider two limiting cases for the distribution of forms: a single Gaussian and a double Gaussian. Alternatively, an abrupt change of form is described by a discontinuous distribution of interaction constants. Our results show that the dispersive distributions yield a biaxial phase, while an abrupt change of shape leads to coexistence of uniaxial phases. Fitting the theory to the experiment for the ternary system KL/decanol/D2O leads to transition lines in very good agreement with experimental results. In order to rationalise the results of the comparison, we analyse temperature and concentration form dependence, which connects micellar and experimental macroscopic parameters. Physically consistent variations of micellar asymmetry, amphiphile partitioning and volume are obtained. To the best of the authors` knowledge, this is the first truly statistical microscopic approach that is able to model experimentally observed lyotropic biaxial nematic phases.

Compressible Sherrington-Kirkpatrick spin-glass model

We introduce a Sherrington-Kirkpatrick spin-glass model with the addition of elastic degrees of freedom. The problem is formulated in terms of an effective four-spin Hamiltonian in the pressure ensemble, which can be treated by the replica method. In the replica-symmetric approximation, we analyze the pressure-temperature phase diagram, and obtain expressions for the critical boundaries between the disordered and the ordered (spin-glass and ferromagnetic) phases. The second-order para-ferromagnetic border ends at a tricritical point, beyond which the transition becomes discontinuous. We use these results to make contact with the temperature-concentration phase diagrams of mixtures of hydrogen-bonded crystals.

Antiferromagnetic spherical spin-glass model

We study the thermodynamic properties and the phase diagrams of a multi-spin antiferromagnetic spherical spin-glass model using the replica method. It is a two-sublattice version of the ferromagnetic spherical p-spin glass model. We consider both the replica-symmetric and the one-step replica-symmetry-breaking solutions, the latter being the most general solution for this model. We find paramagnetic, spin-glass, antiferromagnetic and mixed or glassy antiferromagnetic phases. The phase transitions are always of second order in the thermodynamic sense, but the spin-glass order parameter may undergo a discontinuous change.

Conservative ensembles for nonequilibrium lattice-gas systems

We show how to set up a constant particle ensemble for the steady state of nonequilibrium lattice-gas systems which originally are defined on a constant rate ensemble. We focus on nonequilibrium systems in which particles are created and annihilated on the sites of a lattice and described by a master equation. We consider also the case in which a quantity other than the number of particle is conserved. The conservative ensembles can be useful in the study of phase transitions and critical phenomena particularly discontinuous phase transitions.

Culture-area relation in Axelrod`s model for culture dissemination

Axelrod`s model for culture dissemination offers a nontrivial answer to the question of why there is cultural diversity given that people`s beliefs have a tendency to become more similar to each other`s as they interact repeatedly. The answer depends on the two control parameters of the model, namely, the number F of cultural features that characterize each agent, and the number q of traits that each feature can take on, as well as on the size A of the territory or, equivalently, on the number of interacting agents. Here, we investigate the dependence of the number C of distinct coexisting cultures on the area A in Axelrod`s model, the culture-area relationship, through extensive Monte Carlo simulations. We find a non-monotonous culture-area relation, for which the number of cultures decreases when the area grows beyond a certain size, provided that q is smaller than a threshold value q (c) = q (c) (F) and F a parts per thousand yen 3. In the limit of infinite area, this threshold value signals the onset of a discontinuous transition between a globalized regime marked by a uniform culture (C = 1), and a completely polarized regime where all C = q (F) possible cultures coexist. Otherwise, the culture-area relation exhibits the typical behavior of the species-area relation, i.e., a monotonically increasing curve the slope of which is steep at first and steadily levels off at some maximum diversity value.