Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

NEW CHARACTERIZATIONS OF COMPLETE SPACELIKE SUBMANIFOLDS IN SEMI-RIEMANNIAN SPACE FORMS

In this paper we study n-dimensional complete spacelike submanifolds with constant normalized scalar curvature immersed in semi-Riemannian space forms. By extending Cheng-Yau`s technique to these ambients, we obtain results to such submanifolds satisfying certain conditions on both the squared norm of the second fundamental form and the mean curvature. We also characterize compact non-negatively curved submanifolds in De Sitter space of index p.

Spectral flow and iteration of closed semi-Riemannian geodesics

We introduce the notion of spectral flow along a periodic semi-Riemannian geodesic, as a suitable substitute of the Morse index in the Riemannian case. We study the growth of the spectral flow along a closed geodesic under iteration, determining its asymptotic behavior.

On the semi-Riemannian bumpy metric theorem

We prove the semi-Riemannian bumpy metric theorem using equivariant variational genericity. The theorem states that, on a given compact manifold M, the set of semi-Riemannian metrics that admit only nondegenerate closed geodesics is generic relatively to the C(k)-topology, k=2, ..., infinity, in the set of metrics of a given index on M. A higher-order genericity Riemannian result of Klingenberg and Takens is extended to semi-Riemannian geometry.

Maslov index in semi-Riemannian submersions

We study focal points and Maslov index of a horizontal geodesic gamma : I -> M in the total space of a semi-Riemannian submersion pi : M -> B by determining an explicit relation with the corresponding objects along the projected geodesic pi omicron gamma : I -> B in the base space. We use this result to calculate the focal Maslov index of a (spacelike) geodesic in a stationary spacetime which is orthogonal to a timelike Killing vector field.

Associated Family of G-Structure Preserving Minimal Immersions in Semi-Riemannian Manifolds

We prove the existence of an associated family of G-structure preserving minimal immersions into semi-Riemannian manifolds endowed with a compatible infinitesimally homogeneous G-structure. We will study in more details minimal embeddings into product of space forms.

COMPARISON RESULTS FOR CONJUGATE AND FOCAL POINTS IN SEMI-RIEMANNIAN GEOMETRY VIA MASLOV INDEX

We prove an estimate on the difference of Maslov indices relative to the choice of two distinct reference Lagrangians of a continuous path in the Lagrangian Grassmannian of a symplectic space. We discuss some applications to the study of conjugate and focal points along a geodesic in a semi-Riemannian manifold.

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.

Versatile microanalytical system with porous polypropylene capillary membrane for calibration gas generation and trace gaseous pollutants sampling applied to the analysis of formaldehyde, formic acid, acetic acid and ammonia in outdoor air

The analytical determination of atmospheric pollutants still presents challenges due to the low-level concentrations (frequently in the mu g m(-3) range) and their variations with sampling site and time In this work a capillary membrane diffusion scrubber (CMDS) was scaled down to match with capillary electrophoresis (CE) a quick separation technique that requires nothing more than some nanoliters of sample and when combined with capacitively coupled contactless conductometric detection (C(4)D) is particularly favorable for ionic species that do not absorb in the UV-vis region like the target analytes formaldehyde formic acid acetic acid and ammonium The CMDS was coaxially assembled inside a PTFE tube and fed with acceptor phase (deionized water for species with a high Henry s constant such as formaldehyde and carboxylic acids or acidic solution for ammonia sampling with equilibrium displacement to the non-volatile ammonium ion) at a low flow rate (8 3 nLs(-1)) while the sample was aspirated through the annular gap of the concentric tubes at 25 mLs(-1) A second unit in all similar to the CMDS was operated as a capillary membrane diffusion emitter (CMDE) generating a gas flow with know concentrations of ammonia for the evaluation of the CMDS The fluids of the system were driven with inexpensive aquarium air pumps and the collected samples were stored in vials cooled by a Peltier element Complete protocols were developed for the analysis in air of NH(3) CH(3)COOH HCOOH and with a derivatization setup CH(2)O by associating the CMDS collection with the determination by CE-C(4)D The ammonia concentrations obtained by electrophoresis were checked against the reference spectrophotometric method based on Berthelot s reaction Sensitivity enhancements of this reference method were achieved by using a modified Berthelot reaction solenoid micro-pumps for liquid propulsion and a long optical path cell based on a liquid core waveguide (LCW) All techniques and methods of this work are in line with the green analytical chemistry trends (C) 2010 Elsevier B V All rights reserved

Thermal analysis and compatibility studies of prednicarbate with excipients used in semi solid pharmaceutical form

Differential Scanning Calorimetry (DSC), thermogravimetry/derivative thermogravimetry (TG/DTG) and infrared spectroscopy (IR) techniques were used to investigate the compatibility between prednicarbate and several excipients commonly used in semi solid pharmaceutical form. The thermoanalytical studies of 1:1 (m/m) drug/excipient physical mixtures showed that the beginning of the first thermal decomposition stage of the prednicarbate (T (onset) value) was decreased in the presence of stearyl alcohol and glyceryl stearate compared to the drug alone. For the binary mixture of drug/sodium pirrolidone carboxilate the first thermal decomposition stage was not changed, however the DTG peak temperature (T (peak DTG)) decreased. The comparison of the IR spectra of the drug, the physical mixtures and of the thermally treated samples confirmed the thermal decomposition of prednicarbate. By the comparison of the thermal profiles of 1:1 prednicarbate:excipients mixtures (methylparaben, propylparaben, carbomer 940, acrylate crosspolymer, lactic acid, light liquid paraffin, isopropyl palmitate, myristyl lactate and cetyl alcohol) no interaction was observed.

Synthesis, characterization and adsorption properties of porous mixed valent diruthenium(II,III)-terephthalate and diruthenium(II,III)-adipate polymers

Two porous mixed valent diruthenium(II,III)-dicarboxylate compounds have been prepared and characterized by spectroscopic methods, X-ray diffraction and thermogravimetry. Crystalline solids of [Ru(2)(tere)(2)Cl] center dot 3.5H(2)O (tere=terephthalate) and [Ru(2)(adip)(2)Cl] center dot 1.5H(2)O (adip=adipate) consist of extended chains in which polymeric layers of multiply metal-metal bonded [Ru(2)](5+) cores are bridged by dicarboxylate ligands in paddlewheel type geometries. Units of [Ru(2)(dicarboxylate)(2)](n)(+) are linked by axial bridging chloride ions generating three-dimensional networks. The polymers loose non-bonded water molecules at low temperatures but do not undergo thermal decomposition below 280-300 degrees C. Both of compounds exhibit high BET surface areas, [Ru(2)(tere)(2)Cl]: 235 m(2) g(-1) and [Ru(2)(adip)(2)Cl]: 281 m(2) g(-1), and occlude similar numbers of mol of N(2) per mol of metal. The terephthalate ligand generated an organized structure with supermicropores (total pore size of 0.24 cm(3) g(-1)) while the adipate ligand led to a mesoporous structure (total pore sizes of 0.47 cm(3) g(-1)) for the corresponding diruthenium(II,III)-dicarboxylate polymers. (c) 2008 Elsevier B.V. All rights reserved.

Structure of polyaniline formed in different inorganic porous materials: A spectroscopic study

In this work, the electronic and structural characterization of polyaniline (PANI) formed in cavities of zeolites Y (ZY) and Mordenite (MOR) and montmorillonite (MMT) clay having Cu(II) as oxidant agent are presented. The formation of PANI and its structure is analyzed by Resonance Raman, UV-Vis-NIR, FT-IR and N K XANES techniques. In all cases the structure of PANT formed is different from the ""free"" polymer. The presence of azo bonds linked to phenazine-like rings are observed only for PANI-MMT composites, independent of the kind of oxidant agent employed in the synthesis. The presence of Cu(II) ions leads to the formation of Phenosafranine-like rings. The presence of these phenazine-like rings in the structure of confined PANT chains can also contribute to the enhancement of the thermal stability observed for all composites. (C) 2008 Elsevier Ltd. All rights reserved.

Environmental formaldehyde analysis by active diffusive sampling with a bundle of polypropylene porous capillaries followed by capillary zone electrophoretic separation and contactless conductivity detection

Compared to other volatile carbonylic compounds present in outdoor air, formaldehyde (CH2O) is the most toxic, deserving more attention in terms of indoor and outdoor air quality legislation and control. The analytical determination of CH2O in air still presents challenges due to the low-level concentration (in the sub-ppb range) and its variation with sampling site and time. Of the many available analytical methods for carbonylic compounds, the most widespread one is the time consuming collection in cartridges impregnated with 2,4-dinitrophenylhydrazine followed by the analysis of the formed hydrazones by HPLC. The present work proposes the use of polypropylene hollow porous capillary fibers to achieve efficient CH2O collection. The Oxyphan (R) fiber (designed for blood oxygenation) was chosen for this purpose because it presents good mechanical resistance, high density of very fine pores and high ratio of collection area to volume of the acceptor fluid in the tube, all favorable for the development of air sampling apparatus. The collector device consists of a Teflon pipe inside of which a bundle of polypropylene microporous capillary membranes was introduced. While the acceptor passes at a low flow rate through the capillaries, the sampled air circulates around the fibers, impelled by a low flow membrane pump (of the type used for aquariums ventilation). The coupling of this sampling technique with the selective and quantitative determination of CH2O, in the form of hydroxymethanesulfonate (HMS) after derivatization with HSO3-, by capillary electrophoresis with capacitively coupled contactless conductivity detection (CE-(CD)-D-4) enabled the development of a complete analytical protocol for the CH2O evaluation in air. (C) 2008 Published by Elsevier B.V.

Immobilization of Catalysts of Biological Interest on Porous Oxidized Silicon Surfaces

The present paper deals with the immobilization of redox mediators and proteins onto protected porous silicon surfaces to obtain their direct electrochemical reactions and to retain their bioactivities. This paper shows that MP-11 and viologens are able to establish chemical bonds with 3-aminopropyltriethoxylsilane-modified porous silicon surface. The functionalization of the surfaces have been fully characterized by energy dispersive X-ray analysis (EDX) and X-ray photoelectron spectroscopy (XPS) to examine the immobilization of these mediators onto the solid surface. Amperometric and open circuit potential measurements have shown the direct electron transfer between glucose oxidase and the electrode in the presence of the viologen mediator covalently linked to the 3-aminopropyltriethoxylsilane (APTES)-modified porous silicon surfaces.