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.

Generalized solutions of a nonlinear parabolic equation with generalized functions as initial data

In [H. Brezis, A. Friedman, Nonlinear parabolic equations involving measures as initial conditions, J. Math. Pure Appl. (9) (1983) 73-97.] Brezis and Friedman prove that certain nonlinear parabolic equations, with the delta-measure as initial data, have no solution. However in [J.F. Colombeau, M. Langlais, Generalized solutions of nonlinear parabolic equations with distributions as initial conditions, J. Math. Anal. Appl (1990) 186-196.] Colombeau and Langlais prove that these equations have a unique solution even if the delta-measure is substituted by any Colombeau generalized function of compact support. Here we generalize Colombeau and Langlais` result proving that we may take any generalized function as the initial data. Our approach relies on recent algebraic and topological developments of the theory of Colombeau generalized functions and results from [J. Aragona, Colombeau generalized functions on quasi-regular sets, Publ. Math. Debrecen (2006) 371-399.]. (C) 2009 Elsevier Ltd. All rights reserved.

Surface Activity of the Triflate Ion at the Air/Water Interface and Properties of N,N,N-Trimethyl-N-Dodecylammonium Triflate Aqueous Solutions

The surface activity of salts added to water is Air orders of magnitude lower than that of surfactants. Sodium trifluoromethanesulfonate (NaTf) produced a change in surface tension. with concentration, Delta gamma/Delta c, of -13.2 mN.L/m.mol. This value is ca. 4-fold larger than those of simple salts and that of methanesulfonate. This unexpected surface effect suggested that positively charged micelles containing Tf could exhibit interesting properties. Dodecyltrimethylammonium triflate (DTATf) had a higher Kraft temperature (37 degrees C) and a lower cmc (5 x 10(-3)M) and degree of dissociation (0.11) than the chloride and bromide salts of DTA. Above the Kraft temperature, at a characteristic temperature t(1), the addition of NaTf above 0.05 M. to a DTATf solution induced phase separation. By increasing the temperature of the two-phase system to above t(1), a homogeneous, transparent solution was obtained at a characteristic temperature t(2). These results, together with well-known triflate properties, led us to suggest that the Tf ion pairs With DTA and that the -CF(3) group may be dehydrated in the interfacial region, resulting in new and interesting self-aggregated structures.

Effect of SO(2) on the Transport Properties of an Imidazolium Ionic Liquid and Its Lithium Solution

Transport coefficients have been measured as a function of the concentration of sulfur dioxide, SO(2), dissolved in 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)-imide, [BMMI][Tf(2)N], as well as in its lithium salt solution, Li[Tf(2)N]. The SO(2) reduces viscosity and density and increases conductivity and diffusion coefficients in both the neat [BMMI] [Tf(2)N] and the [BMMI][Tf(2)N]-Li[Tf(2)N] solution. The conductivity enhancement is not assigned to a simple viscosity effect; the weakening of ionic interactions upon SO(2) addition also plays a role. Microscopic details of the SO(2) effect were unraveled using Raman spectroscopy and molecular dynamics (MD) simulations. The Raman spectra suggest that the Li(+)-[Tf(2)N] interaction is barely affected by SO(2), and the SO(2)-[Tf(2)N] interaction is weaker than previously observed in an investigation of an ionic liquid containing the bromide anion. Transport coefficients calculated by MD simulations show the same trend as the experimental data with respect to SO(2) content. The MD simulations provide structural information on SO(2) molecules around [Tf(2)N], in particular the interaction of the sulfur atom of SO(2) with oxygen and fluorine atoms of the anion. The SO(2)-[BMMI] interaction is also important because the [BMMI] cations with above-average mobility have a larger number of nearest-neighbor SO(2) molecules.

Room Temperature Ionic Liquid-Lithium Salt Mixtures: Optical Kerr Effect Dynamical Measurements

The addition of lithium salts to ionic liquids causes an increase in viscosity and a decrease in ionic mobility that hinders their possible application as an alternative solvent in lithium ion batteries. Optically heterodyne-detected optical Kerr effect spectroscopy was used to study the change in dynamics, principally orientational relaxation, caused by the addition of lithium bis(trifluoromethylsulfonyl)imide to the ionic liquid 1-buty1-3-methylimidazolium bis(trifluoromethylsulfonyl)imide. Over the time scales studied (1 ps-16 ns) for the pure ionic liquid, two temperature-independent power laws were observed: the intermediate power law (1 ps to similar to 1 ns), followed by the von Schweidler power law. The von Schweidler power law is followed by the final complete exponential relaxation, which is highly sensitive to temperature. The lithium salt concentration, however, was found to affect both power laws, and a discontinuity could be found in the trend observed for the intermediate power law when the concentration (mole fraction) of lithium salt is close to chi(LiTf(2)N) = 0.2. A mode coupling theory (MCT) schematic model was also used to fit the data for both the pure ionic liquid and the different salt concentration mixtures. It was found that dynamics in both types of liquids are described very well by MCT.

Effect of solvent on the adsorption behavior and on the surface properties of Sterculia striata polysaccharide

Characterization of Sterculia striate polysaccharide (SSP) films adsorbed onto Si wafers from solutions prepared in ethyl methyl imidazolium acetate (EmimAc), water or NaOH 0.01 mol/L was systematically studied by means of ellipsometry, atomic force microscopy and contact angle measurements. SSP adsorbed from EmimAc onto Si wafer as homogeneous monolayers (similar to 0.5 nm thick), while from water or NaOH 0.01 mol/L SSP formed layers of similar to 4.0 nm and similar to 1.5 nm thick, respectively. Surface energy values found for SSP adsorbed from EmimAc or water were 68 +/- 2 mJ/m(2) and 65 +/- 2 mJ/m(2), respectively, whereas from NaOH it amounted to 57 +/- 3 mJ/m(2). The immobilization of lysozyme (LYS) onto SSP films was also investigated. The mean thickness of LYS (d(LYS)) immobilized onto SSP films adsorbed from each solvent tended to increase with the decrease of gamma(P)(S) and gamma(total)(S). However, the enzymatic activity of LYS molecules was higher when they were immobilized onto SSP films with higher gamma(P)(S) and gamma(total)(S) values. (C) 2010 Elsevier Ltd. All rights reserved.

Unraveling Dynamical Heterogeneity in the Ionic Liquid 1-Butyl-3-methylimidazolium Chloride

Heterogeneous dynamics within a time range of nanoseconds was investigated by molecular dynamics (MD) simulations of 1 -butyl-3-methylimidazolium chloride ([bmim]Cl). After identifying groups of fast and slow ions, it was shown that the separation between the location of the center of mass and the center of charge of cations, d(CMCC), is a signature of such difference in ionic mobility. The distance d(CMCC) can be used as a signature because it relaxes in the time window of the dynamical heterogeneity. The relationship between the parameter dcmcc and conformations of the side alkyl chain in [bmim] is discussed. Since the relatively slow relaxation of dcmcc is a consequence of coexisting polar and nonpolar domains in the bulk, the MD simulations reveal a subtle interplay between structural and dynamical heterogeneity in ionic liquids.

Platinum nanoparticle-modified electrodes, morphologic, and electrochemical studies concerning electroactive materials deposition

The present work describes the synthesis of platinum nanoparticles followed by their electrophoretic deposition onto transparent fluorine-doped tin oxide electrodes. The nano-Pt-modified electrodes were characterized by voltammetric studies in acidic solutions showing a great electrocatalytic behavior towards H(+) reduction being very interesting for fuel cell applications. Morphological characterization was performed by atomic force microscopy on different modified electrodes showing a very rough surface which can be tuned by means of time of deposition. Also, nickel hydroxide thin films were galvanostatically grown onto these electrodes showing an interesting electrochemical behavior as sharper peaks, indicating a faster ionic exchange from the electrolyte to the film.

Application of Microelectrode Voltammetry to Study the Properties of Surfactant Solutions: Alkyltrimethylammonium Bromides

Microelectrode cyclic voltammetry (MV) has been employed to investigate the micellar properties of solutions of homologous alkyltrimethylammonium bromides, RMe(3)ABr, R = C(10), C(12), and C(14), in water and in the presence of added NaBr. The micellar self-diffusion coefficient was calculated from the limiting current for the reversible electron transfer of micelle-bound ferrocene. From the values of this property, other parameters were calculated, including the micellar hydrodynamic radius, RH, and aggregation number, N(agg); the latter was also theoretically calculated. We determined the values of the diffusion coefficient as a function of various experimental variables and observed the following trends: The diffusion coefficient decreases as a function of increasing surfactant concentration (no additional electrolyte added); it decreases as a function of increasing surfactant concentration at fixed NaBr concentration; and it shows a complex dependence (increase then decrease) on the NaBr concentration at a fixed RMe(3)ABr concentration. The value of the intermicellar interaction parameter decreases and then increases as a function of increasing NaBr concentration. These results are discussed in terms of intermicellar,interactions and the effect of NaBr on the micellar surface charge density and sphere-to-rod geometry change. The NaBr concentration required to induce the latter change increases rapidly as a function of decreasing the length of R: no geometry change was detected for C(10)Me(3)ABr. Values of N(agg) increase as I function of increasing the length of R and are in good agreement with both literature values and values that were calculated theoretically. Thus, MV is a convenient and simple technique for obtaining fundamental properties of surfactant solutions, including additive-induced changes of micellar parameters (N(agg)) and morphology changes.

Nickel hydroxide electrodes as amperometric detectors for carbohydrates in flow injection analysis and liquid chromatography

The present paper describes the utilization of nickel hydroxide modified electrodes toward the catalytic oxidation of carbohydrates (glucose, fructose, lactose and sucrose) and their utilization as electrochemical sensor. The modified electrodes were employed as a detector in flow injection analysis for individual carbohydrate detection, and to an ionic column chromatography system for multi-analyte samples aiming a prior separation step. Kinetic studies were performed on a rotating disk electrode (RDE) in order to determine both the heterogeneous rate constant and number of electrons transferred for each carbohydrate. Many advantages were found for the proposed system including fast and easy handling of the electrode modification, low cost procedure, a wide range of linearity (0.5-50 ppm), low detection limits (ppb level) and high sensitivities. (C) 2009 Elsevier B.V. All rights reserved.

V(2)O(5) nanoparticles obtained from a synthetic bariandite-like vanadium oxide: Synthesis, characterization and electrochemical behavior in an ionic liquid

Highly stable and crystalline V(2)O(5) nanoparticles with an average diameter of 15 nm have been easily prepared by thermal treatment of a bariandite-like vanadium oxide, V(10)O(24)center dot 9H(2)O. Their characterization was carried out by powder X-ray diffractometry (XRD). Fourier transform infrared (FT-IR) and Raman spectroscopies, and transmission electron microscopy (TEM). The fibrous and nanostructured film obtained by electrophoretic deposition of the V(2)O(5) nanoparticles showed good electroactivity when submitted to cyclic voltammetry in an ionic liquid-based electrolyte. The use of this film for the preparation of a nanostructured electrode led to an improvement of about 50% in discharge capacity values when compared with similar electrodes obtained by casting of a V(2)O(5) xerogel. (C) 2009 Elsevier Inc. All rights reserved.

PVP superabsorbent nanogels

Monomer free hydrogel nanoparticles (nanogels) were prepared by crosslinking preformed poly(N-vinyl-2-pyrrolidone) (PVP) entrapped in the aqueous pool of hexadecyltrimethylammonium bromide reverse micelles using the Fenton reaction. The PVP nanoparticles were spherical with a dry diameter of 27 nm. The diameter of the swollen particles was ten times higher, i.e., a swelling ratio, Q, above 900, characterizing this preparation as superabsorbent. PVP nanogel swelling was dependent on bound Fe(3+) and varied with pH and ionic strength. Nanogel deswelling by salt followed the anions lyotropic series, i.e., SCN(-) < HSO(3)(-) < NO(3)(-) < I(-) < Cl(-) < CH(3)COO(-) < CF(3)SO(3)(-). The value of Q reached 6,000 in iron-free PVP nanoparticles at low pH, making this nanogel one of the most efficient swelling systems so far described.

Prigogine-Defay Ratio for an Ionic Glass-Former: Molecular Dynamics Simulations

The pressure dependence of the glass-transition temperature, T(g)(P), of the ionic glass-former 2Ca(NO(3))(2) center dot 3KNO(3), CKN, has been obtained by molecular dynamics (MD) simulations The liquid-glass difference of thermal expansivity, Delta alpha, heat capacity, Delta C(p), and isothermal compressibility, Delta kappa, have been calculated as a function of pressure. It has been found that the Ehrenfest relation dT(g)/dP = TV Delta alpha/Delta C(p) predicts the pressure dependence of T, but the other Ehrenfest relation, dT(g)/dP = Delta kappa/Delta alpha, does not. Consequently, the Prigogine-Defay ratio, Pi = Delta C(p)Delta kappa/TV Delta alpha(2), is Pi similar to 1.2 at low pressures, but increases 1 order of magnitude at high pressures. The pressure dependence of the Prigogine-Defay ratio is interpreted in light of recent explanations for the finding Pi > 1.

Alkoxy Chain Effect on the Viscosity of a Quaternary Ammonium Ionic Liquid: Molecular Dynamics Simulations

The viscosity of ionic liquids based on quaternary ammonium cations is reduced when one of the alkyl chains is replaced by an alkoxy chain (Zhou et al. Chem. Eur. J. 2005, 11, 752.). A microscopic picture of the role played by the ether function in decreasing the viscosity of quaternary ammonium ionic liquids is provided here by molecular dynamics (MD) simulations. A model for the ionic liquid N-ethyl-N,N-dimethyl-N-(2-methoxyethyl)ammonium bis(trifluoromethanesulfonyl)imide, MOENM(2)E TFSI, is compared to the tetraalky-lammonium counterpart. The alkoxy derivative has lower viscosity, higher ionic diffusion coefficients, and higher conductivity than the tetraalkyl system at the same density and temperature. A clear signature of the ether function on the liquid structure is observed in cation-cation correlations, but not in anion-anion or anion-cation correlations. In both the alkyl and the alkoxy ionic liquids, there is aggregation of long chains of neighboring cations within micelle-like structures. The MD simulations indicate that the less effective assembly between the more flexible alkoxy chains, in comparison to alkyl chains, is the structural reason for higher ionic mobility in MOENM(2)E TFSI.

Characterization of anodic silicon oxide films grown in room temperature ionic liquids

The electroformation of silicon oxide was performed in two room temperature ionic liquids (RTIL), 1-butyl-3-methyl-imidazolium bis(trifluoromethane sulfonyl) imide (BMITFSI) and N-n-butyl-N-methylpiperidinium bis(trifluoromethane sulfonyl) imide (BMPTFSI). This phenomenon was studied by electrochemical techniques and it was observed that the oxide growth follows a high-field mechanism. X-ray Photoelectron Spectroscopy experiments have shown that a non-stoichiometric oxide film was formed, related to the low water content present in both RTILs (< 30 ppm). The roughness values obtained by using AFM technique of the silicon surface after etching with HF was 1.5 nm (RMS). The electrochemical impedance spectroscopy at low frequencies range was interpreted as a resistance in parallel with a CPE element, the capacitance obtained was associated with the dielectric nature of the oxide formed and the resistance was interpreted considering the chemical dissolution of the oxide by the presence of the TFSI anion. The CPE element was associated with the surface roughness and the very thin oxide film obtained. (C) 2007 Elsevier Ltd. All rights reserved.