Energy consistent discontinuous Galerkin methods for a quasi-incompressible diffuse two phase flow model

We design consistent discontinuous Galerkin finite element schemes for the approximation of a quasi-incompressible two phase flow model of Allen–Cahn/Cahn–Hilliard/Navier–Stokes–Korteweg type which allows for phase transitions. We show that the scheme is mass conservative and monotonically energy dissipative. In this case the dissipation is isolated to discrete equivalents of those effects already causing dissipation on the continuous level, that is, there is no artificial numerical dissipation added into the scheme. In this sense the methods are consistent with the energy dissipation of the continuous PDE system.

Non-ergodicity of inviscid two-dimensional flow on a beta-plane and on the surface of a rotating sphere

It is shown that, for a sufficiently large value of β, two-dimensional flow on a doubly-periodic beta-plane cannot be ergodic (phase-space filling) on the phase-space surface of constant energy and enstrophy. A corresponding result holds for flow on the surface of a rotating sphere, for a sufficiently rapid rotation rate Ω. This implies that the higher-order, non-quadratic invariants are exerting a significant influence on the statistical evolution of the flow. The proof relies on the existence of a finite-amplitude Liapunov stability theorem for zonally symmetric basic states with a non-vanishing absolute-vorticity gradient. When the domain size is much larger than the size of a typical eddy, then a sufficient condition for non-ergodicity is that the wave steepness ε < 1, where ε = 2[surd radical]2Z/βU in the planar case and $\epsilon = 2^{\frac{1}{4}} a^{\frac{5}{2}}Z^{\frac{7}{4}}/\Omega U^{\frac{5}{2}}$ in the spherical case, and where Z is the enstrophy, U the r.m.s. velocity, and a the radius of the sphere. This result may help to explain why numerical simulations of unforced beta-plane turbulence (in which ε decreases in time) seem to evolve into a non-ergodic regime at large scales.

Liquid crystal phase formation by biopolymers

This review discusses liquid crystal phase formation by biopolymers in solution. Lyotropic mesophases have been observed for several classes of biopolymer including DNA, peptides, polymer/peptide conjugates, glycopolymers and proteoglycans. Nematic or chiral nematic (cholesteric) phases are the most commonly observed mesophases, in which the rod-like fibrils have only orientational order. Hexagonal columnar phases are observed for several systems (DNA, PBLG, polymer/peptide hybrids) at higher concentration. Lamellar (smectic) phases are reported less often, although there are examples such as the layer arrangement of amylopectin side chains in starch. Possible explanations for the observed structures are discussed. The biological role of liquid crystal phases for several of these systems is outlined. Commonly, they may serve as a template to align fibrils for defined structural roles when the biopolymer is extruded and dried, for instance in the production of silk by spiders or silkworms, or of chitin in arthropod shells. In other cases, liquid crystal phase formation may occur in vivo simply as a consequence of high concentration, for instance the high packing density of DNA within cell nuclei.

Sensitivity analysis of a cyanobacterial growth and movement model under two different flow regimes

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling cyanobacterial behaviour in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes, reservoirs and rivers. A new deterministic–mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including light, nutrients and temperature. A parameter sensitivity analysis using a one-at-a-time approach was carried out. There were two objectives of the sensitivity analysis presented in this paper: to identify the key parameters controlling the growth and movement patterns of cyanobacteria and to provide a means for model validation. The result of the analysis suggested that maximum growth rate and day length period were the most significant parameters in determining the population growth and colony depth, respectively.

The role of phase separation and feed cycle length in leach beds coupled to methanogenic reactors for digestion of a solid substrate (Part 2): Hydrolysis, acidification and methanogenesis in a two-phase system

The effect of phase separation and batch duration on the trophic stages of anaerobic digestion was assessed for the first time in leach beds coupled to methanogenic reactors digesting maize (Zea mays). The system was operated for consecutive batches of 7, 14 and 28 days for ~120 days. Hydrolysis rate was higher the shorter the batch, reaching 8.5 gTSdestroyed d-1 in the 7-day system. Phase separation did not affect acidification but methanogenesis was enhanced in the short feed cycle leach beds. Phase separation was inefficient on the 7-day system, where ~89% of methane was produced in the leach bed. Methane production rate increased with shortening the feed cycle, reaching 3.523 l d-1 average in the 7-day system. Low strength leachate from the leach beds decreased methanogenic activity of methanogenic reactors’ sludges. Enumeration of cellulolytic and methanogenic microorganisms indicated a constant inoculation of leach beds and methanogenic reactors through leachate recirculation.

Recovery of lactoferrin and lactoperoxidase from sweet whey using colloidal gas aphrons (CGAs) generated from an anionic surfactant, AOT

The recovery of lactoferrin and lactoperoxidase from sweet whey was studied using colloidal gas aphrons (CGAs), which are surfactant-stabilized microbubbles (10-100 mum). CGAs are generated by intense stirring (8000 rpm for 10 min) of the anionic surfactant AOT (sodium bis-2-ethylhexyl sulfosuccinate). A volume of CGAs (10-30 mL) is mixed with a given volume of whey (1 - 10 mL), and the mixture is allowed to separate into two phases: the aphron (top) phase and the liquid (bottom) phase. Each of the phases is analyzed by SDS-PAGE and surfactant colorimetric assay. A statistical experimental design has been developed to assess the effect of different process parameters including pH, ionic strength, the concentration of surfactant in the CGAs generating solution, the volume of CGAs and the volume of whey on separation efficiency. As expected pH, ionic strength and the volume of whey (i.e. the amount of total protein in the starting material) are the main factors influencing the partitioning of the Lf(.)Lp fraction into the aphron phase. Moreover, it has been demonstrated that best separation performance was achieved at pH = 4 and ionic strength = 0.1 mol/L i.e., with conditions favoring electrostatic interactions between target proteins and CGAs (recovery was 90% and the concentration of lactoferrin and lactoperoxidase in the aphron phase was 25 times higher than that in the liquid phase), whereas conditions favoring hydrophobic interactions (pH close to pI and high ionic strength) led to lower performance. However, under these conditions, as confirmed by zeta potential measurements, the adsorption of both target proteins and contaminant proteins is favored. Thus, low selectivity is achieved at all of the studied conditions. These results confirm the initial hypothesis that CGAs act as ion exchangers and that the selectivity of the process can be manipulated by changing main operating parameters such as type of surfactant, pH and ionic strength.

Preferred structures in large-scale circulation and the effect of doubling greenhouse gas concentration in HadCM3

Preferred structures in the surface pressure variability are investigated in and compared between two 100-year simulations of the Hadley Centre climate model HadCM3. In the first (control) simulation, the model is forced with pre-industrial carbon dioxide concentration (1×CO2) and in the second simulation the model is forced with doubled CO2 concentration (2×CO2). Daily winter (December-January-February) surface pressures over the Northern Hemisphere are analysed. The identification of preferred patterns is addressed using multivariate mixture models. For the control simulation, two significant flow regimes are obtained at 5% and 2.5% significance levels within the state space spanned by the leading two principal components. They show a high pressure centre over the North Pacific/Aleutian Islands associated with a low pressure centre over the North Atlantic, and its reverse. For the 2×CO2 simulation, no such behaviour is obtained. At higher-dimensional state space, flow patterns are obtained from both simulations. They are found to be significant at the 1% level for the control simulation and at the 2.5% level for the 2×CO2 simulation. Hence under CO2 doubling, regime behaviour in the large-scale wave dynamics weakens. Doubling greenhouse gas concentration affects both the frequency of occurrence of regimes and also the pattern structures. The less frequent regime becomes amplified and the more frequent regime weakens. The largest change is observed over the Pacific where a significant deepening of the Aleutian low is obtained under CO2 doubling.

Self-assembly of two-component gels: stoichiometric control and component selection

Two-component systems capable of self-assembling into soft gel-phase materials are of considerable interest due to their tunability and versatility. This paper investigates two-component gels based on a combination of a L-lysine-based dendron and a rigid diamine spacer (1,4-diaminobenzene or 1,4-diaminocyclohexane). The networked gelator was investigated using thermal measurements, circular dichroism, NMR spectroscopy and small angle neutron scattering (SANS) giving insight into the macroscopic properties, nanostructure and molecular-scale organisation. Surprisingly, all of these techniques confirmed that irrespective of the molar ratio of the components employed, the "solid-like" gel network always consisted of a 1:1 mixture of dendron/diamine. Additionally, the gel network was able to tolerate a significant excess of diamine in the "liquid-like" phase before being disrupted. In the light of this observation, we investigated the ability of the gel network structure to evolve from mixtures of different aromatic diamines present in excess. We found that these two-component gels assembled in a component-selective manner, with the dendron preferentially recognising 1,4-diaminobenzene (>70%). when similar competitor diamines (1,2- and 1,3-diaminobenzene) are present. Furthermore, NMR relaxation measurements demonstrated that the gel based oil 1,4-diaminobenzene was better able to form a selective ternary complex with pyrene than the gel based oil 1,4-diaminocyclohexane, indicative of controlled and selective pi-pi interactions within a three-component assembly. As such, the results ill this paper demonstrate how component selection processes in two-component gel systems call control hierarchical self-assembly.

Chapter 17: X-Ray rheology of liquid crystal polymers

X-ray Rheology is an experimental technique which uses time-ressolved x-ray scattering as probe of the molecular level structural reorganisation which accompanies flow. It provides quantitative information on the direction alignment and on the level of global orientation. This information is very helpful in interpreting the classic rheological data on liquid crystal polymers. In this research we use data obtained from a cellulose derivate which exhibits a thermotropic liquid crystal phase. We show how increased shear rates lead to a rapid rise in the global orientation and we related this to therories of flow in liquid crystal polymers from the literature. We show that the relaxation time is independent of the prior shear rate.

An HMV mapping study of a pulsating jet system

A new system for the generation of hydrodynamic modulated voltammetry (HMV) is presented. This system consists of an oscillating jet produced through the mechanical vibration of a large diaphragm. The structure of the cell is such that a relatively small vibration is transferred to a large fluid flow at the jet outlet. Positioning of an electrode (Pt, 0.5 mm or 25 mu m diameter) over the exit of this jet enables the detection of the modulated flow of liquid. While this flow creates modest mass transfer rates (time averaged similar to 0.015 cm s(-1)) it can also be used to create a HMV system where a 'lock-in' approach is adopted to investigate the redox chemistry in question. This is demonstrated for the Fe(CN)(6)(3-/4-) redox system. Here 'lock-in' to the modulated hydrodynamic signal is achieved through the deployment of bespoke software. The apparatus and procedure is shown to produce a simple and efficient way to obtain the desired signal. In addition the spatial variation of the HMV signal, phase correction and time averaged current with respect to the jet orifice is presented. (C) 2008 Elsevier B.V. All rights reserved.

A moving-mesh finite element method and its application to the numerical solution of phase-change problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase. The algorithm uses a distributed conservation principle to determine nodal mesh velocities, which are then used to move the nodes. The nodal values are obtained from an ALE (Arbitrary Lagrangian-Eulerian) equation, which represents a generalization of the original algorithm presented in Applied Numerical Mathematics, 54:450--469 (2005). Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and, for the first time, two-phase Stefan problems in one and two space dimensions, paying particular attention to the implementation of the interface boundary conditions. Results are presented to demonstrate the accuracy and the effectiveness of the method, including comparisons against analytical solutions where available.

Particle image velocimetry measurements of turbulent flow within outdoor and indoor urban scale models and flushing motions in urban canopy layers

We investigate the spatial characteristics of urban-like canopy flow by applying particle image velocimetry (PIV) to atmospheric turbulence. The study site was a Comprehensive Outdoor Scale MOdel (COSMO) experiment for urban climate in Japan. The PIV system captured the two-dimensional flow field within the canopy layer continuously for an hour with a sampling frequency of 30 Hz, thereby providing reliable outdoor turbulence statistics. PIV measurements in a wind-tunnel facility using similar roughness geometry, but with a lower sampling frequency of 4 Hz, were also done for comparison. The turbulent momentum flux from COSMO, and the wind tunnel showed similar values and distributions when scaled using friction velocity. Some different characteristics between outdoor and indoor flow fields were mainly caused by the larger fluctuations in wind direction for the atmospheric turbulence. The focus of the analysis is on a variety of instantaneous turbulent flow structures. One remarkable flow structure is termed 'flushing', that is, a large-scale upward motion prevailing across the whole vertical cross-section of a building gap. This is observed intermittently, whereby tracer particles are flushed vertically out from the canopy layer. Flushing phenomena are also observed in the wind tunnel where there is neither thermal stratification nor outer-layer turbulence. It is suggested that flushing phenomena are correlated with the passing of large-scale low-momentum regions above the canopy.

Electro-optical properties of liquid crystal copolymers and their relationship to structural order

A range of side chain liquid crystal copolymers have been prepared using mesogenic and non-mesogenic units. It is found that high levels of the non-mesogenic moieties may be introduced without completely disrupting the organization of the liquid crystal phase. Incorporation of this comonomer causes a marked reduction in the glass transition temperature (Tg), presumably as a result of enhanced backbone mobility and a corresponding lowering of the nematic transition temperature, thereby restricting the temperature range for stability of the liquid crystal phase. The effect of the interactions between the various components of these side-chain polymers on their electro-optic responses is described. Infrared (i.r.) dichroism measurements have been made to determine the order parameters of the liquid crystalline side-chain polymers. By identifying a certain band (CN stretching) in the i.r. absorption spectrum, the order parameter of the mesogenic groups can be obtained. The temperature and composition dependence of the observed order parameter are related to the liquid crystal phase transitions and to the electro-optic response. It is found that the introduction of the non-mesogenic units into the polymer chain lowers the threshold voltage of the electro-optic response over and above that due to the reduction in the order parameter. The dynamic electro-optic responses are dominated by the temperature-dependent viscosity and evidence is presented for relaxation processes involving the polymer backbone which are on a time scale greater than that for the mesogenic side-chain units.