Numerical approximation of generalized Newtonian fluids using Powell–Sabin–Heindl elements: I. theoretical estimates

In this paper we consider the numerical approximation of steady and unsteady generalized Newtonian fluid flows using divergence free finite elements generated by the Powell–Sabin–Heindl elements. We derive a priori and a posteriori finite element error estimates and prove convergence of the method of successive approximations for the steady flow case. A priori error estimates of unsteady flows are also considered. These results provide a theoretical foundation and supporting numerical studies are to be provided in Part II. Copyright © 2003 John Wiley & Sons, Ltd.     

Förster energy transfer from poly(arylene–ethynylene)s to an erbium–porphyrin complex

We study the infrared emission at 1.54 μm of an organolanthanide complex, Er(III)-tetraphenylporphyrin [Er(TPP)acac], both as a result of direct optical excitation and via energy transfer from host π-conjugate polymers of type poly(arylene–ethynylene) [PAE]. In the first case, the emission of the neat complex is characterized in inert transparent materials and a value of the quantum yield at 1.54 μm φ_IR=4×10⁻⁴ is measured. Then, fluorescence resonance transfer is investigated in blends of Er(TPP)acac with PAEs by monitoring the quenching of the polymer fluorescence along with the enhancement of both the visible emission of the ligand and the near-infrared band of Er³⁺. These different procedures allow a detailed analysis of the transfer efficiency within a specific implementation of the Förster model for polymeric donors. The experimental values of the critical radius R₀, ranging from 1.3 to 2.5 nm for the different blends, are in good agreement with theory for a wide interval of the physical and spectroscopic parameters. This suggests that other mechanisms for excitation transfer do not play a significant role in these materials.     

Charge transport properties in carbon black/polymer composites

The charge transport properties of polymer matrix–carbon black composites are investigated in this study. Direct current conductivity is examined with varying parameters: the temperature and the conductive filler content. Conductivity data are analyzed by means of percolation theory, and both percolation threshold and critical exponent are determined at each of the examined temperatures. The temperature dependence of conductivity and the agreement of experimental results with the variable range hopping model reveal hopping conduction as the predominant transport mechanism, below and in the vicinity of the critical concentration of carbon black particles. At higher concentrations, the contribution of hopping transport to the overall conductivity is reduced and a balance between hopping and conduction via geometrical contact occurs. © 2007 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 2535–2545, 2007     

Volumetric and infrared measurements on amorphous ice structure

We have simultaneously used adsorption isotherm volumetry and Fourier transform infrared spectroscopy in order to take the investigations on amorphous ice structure a step further, especially concerning porosity and annealing-induced modifications. We have studied surface reorganization during annealing and found that the number of surface sites decreases before crystallization, their relative ratios being different for amorphous and crystalline ice. We also present results confirming that ice can have a large specific surface area and nevertheless be non-microporous.     

Combined supercritical-fluid chromatography/mass spectrometry in the analysis of diuron in plasma using on-line phase-system switching

On-line sample pretreatment by means of the phase-system switching approach is an interesting technique for the analysis of aqueous samples, e.g., plasma, by means of supercritical-fluid chromatography. In order to analyse plasma samples the following analytical procedure is used. The plasma sample is injected on to a short precolumn, which is washed with water and subsequently dried with nitrogen. Next, the solutes are desorbed with the supercritical mobile phase, analysed with packed-column supercritical-fluid chromatography and detected with either a UV detector or a mass spectrometer, equipped with a moving-belt interface. The herbicide diuron is selected as a test compound to study the feasibility of this approach. Using a selective detector the procedure is sufficiently sensitive to detect diuron in plasma, but not appropriate to detect the diuron metabolites in a post-mortem plasma sample. These have been identified with liquid chromatography/mass spectrometry. The detection limit of diuron in plasma using the procedure described is about 30 ng/mL.     

Fractional Derivatives Applied to Phase-Space Reconstructions

The concept and application of phase-space reconstructions are reviewed. Fractional derivatives are then proposed for the purpose of reconstructing dynamics from a single observed time history. A procedure is presented in which the fractional derivatives of time series data are obtained in the frequency domain. The method is applied to the Lorenz system. The ability of the method to unfold the data is assessed by the method of global false nearest neighbors. The reconstructed data is used to compute recurrences and correlation dimensions. The reconstruction is compared to the commonly used method of delays in order to assess the choice of reconstruction parameters, and also the quality of results.     

Numerical simulations of a filament in a flowing soap film

Experiments concerning the properties of soap films have recently been carried out and these systems have been proposed as experimental versions of theoretical two‐dimensional liquids. A silk filament introduced into a flowing soap film, was seen to demonstrate various stable modes, and these were, namely, a mode in which the filament oscillates and one in which the filament is stationary and aligns with the flow of the liquid. The system could be forced from the oscillatory mode into the non‐ oscillatory mode by varying the length of the filament. In this article we use numerical and computational techniques in order to simulate the strongly coupled behaviour of the filament and the fluid. Preliminary results are presented for the specific case in which the filament is seen to oscillate continuously for the duration of our simulation. We also find that the filament oscillations are strongly suppressed when we reduce the effective length of the filament. We believe that these results are reminiscent of the different oscillatory and non‐oscillatory modes observed in experiment. The numerical solutions show that, in contrast to experiment, vortices are created at the leading edge of the filament and are preferentially grown in the curvature of the filament and are eventually released from the trailing edge of the filament. In a similar manner to oscillating hydrofoils, it seems that the oscillating filaments are in a minimal energy state, extracting sufficient energy from the fluid to oscillate. In comparing numerical and experimental results it is possible that the soap film does have an effect on the fluid flow especially in the boundary layer where surface tension forces are large. Copyright © 2004 John Wiley & Sons, Ltd.     

A numerical and experimental investigation of the interactions between a non‐uniform planar array of incompressible free jets

The interaction between multiple incompressible air jets has been studied numerically and experimentally. The numerical predictions have been first validated using experimental data for a single jet configuration. The spreading features of five unequal jets in the configuration of one larger central jet surrounded by four smaller equi‐distant jets, have been studied, for different lateral spacing ratios of 1.5, 2.0 and 2.5 and a central jet Reynolds number of 1.24×105 (corresponding to a Mach number of 0.16). Flow of five equal jets has also been simulated, for the sake of comparison. The jet interactions commence at an axial distance of about 3–4 diameters and complete by an axial distance of about 10 diameters for the lowest spacing ratio of 1.5. For larger spacing ratios, the length required for the start and completion of jet interaction increase. Peripheral jets bend more towards the central jet and merge at a smaller distance, when their sizes are smaller than that of the central jet. The entrainment ratio for multiple jets is higher than that for a single jet. Excellent agreement is observed between the experimental data and theoretical predictions for both mean flow field and turbulent quantities, at regions away from the jet inlet. The potential core length and initial jet development, however, are not predicted very accurately due to differences in the assumed and actual velocity profiles at the jet inlet. Copyright © 2004 John Wiley & Sons, Ltd.     

An improved numerical simulator for multiphase flow in porous media

In basin modelling the thermodynamics of a multicomponent multiphase fluid flux are computationally too expensive when derived from an equation of state and the Gibbs equality constraints. In this article we present a novel implicit molar mass formulation technique using binary mixture thermodynamics. The two proposed solution methods, with and without cross derivative terms between components, are based on a preconditioned Newton‐GMRES scheme for each time‐step with analytical computation of the derivatives. These new algorithms reduce significantly the numerical effort for the computation of the molar masses, and we illustrate the behavior of these methods with numerical computations. Copyright © 2004 John Wiley & Sons Ltd.