Reconstruction of an unsteady flow from incomplete boundary data:Part I. Theory

Purpose – To propose and investigate a stable numerical procedure for the reconstruction of the velocity of a viscous incompressible fluid flow in linear hydrodynamics from knowledge of the velocity and fluid stress force given on a part of the boundary of a bounded domain. Design/methodology/approach – Earlier works have involved the similar problem but for stationary case (time-independent fluid flow). Extending these ideas a procedure is proposed and investigated also for the time-dependent case. Findings – The paper finds a novel variation method for the Cauchy problem. It proves convergence and also proposes a new boundary element method. Research limitations/implications – The fluid flow domain is limited to annular domains; this restriction can be removed undertaking analyses in appropriate weighted spaces to incorporate singularities that can occur on general bounded domains. Future work involves numerical investigations and also to consider Oseen type flow. A challenging problem is to consider non-linear Navier-Stokes equation. Practical implications – Fluid flow problems where data are known only on a part of the boundary occur in a range of engineering situations such as colloidal suspension and swimming of microorganisms. For example, the solution domain can be the region between to spheres where only the outer sphere is accessible for measurements. Originality/value – A novel variational method for the Cauchy problem is proposed which preserves the unsteady Stokes operator, convergence is proved and using recent for the fundamental solution for unsteady Stokes system, a new boundary element method for this system is also proposed.

The structure and dynamics of 2-dimensional fluids in swelling clays

The interlayer pores of swelling 2:1 clays provide an ideal 2-dimensional environment in which to study confined fluids. In this paper we discuss our understanding of the structure and dynamics of interlayer fluid species in expanded clays, based primarily on the outcome of recent molecular modelling and neutron scattering studies. Counterion solvation is compared with that measured in bulk solutions, and at a local level the cation-oxygen coordination is found to be remarkably similar in these two environments. However, for the monovalent ions the contribution to the first coordination shell from the clay surfaces increases with counterion radius. This gives rise to inner-sphere (surface) complexes in the case of potassium and caesium. In this context, the location of the negative clay surface charge (i.e. arising from octahedral or tetrahedral substitution) is also found to be of major importance. Divalent cations, such as calcium, eagerly solvate to form outer-sphere complexes. These complexes are able to pin adjacent clay layers together, and thereby prevent colloidal swelling. Confined water molecules form hydrogen bonds to each other and to the clays' surfaces. In this way their local environment relaxes to close to the bulk water structure within two molecular layers of the clay surface. Finally, we discuss the way in which the simple organic molecules methane, methanol and ethylene glycol behave in the interlayer region of hydrated clays. Quasi-elastic neutron scattering of isotopically labelled interlayer CH 3OD and (CH2OD)2 in deuterated clay allows us to measure the diffusion of the CH3- and CH2-groups in both clay and liquid environments. We find that in both the one-layer methanol solvates and the two-layer glycol solvates the diffusion of the most mobile organic molecules is close to that in the bulk solution.

A new fast generalized sphere decoding algorithm in MIMO systems

A new generalized sphere decoding algorithm is proposed for underdetermined MIMO systems with fewer receive antennas N than transmit antennas M. The proposed algorithm is significantly faster than the existing generalized sphere decoding algorithms. The basic idea is to partition the transmitted signal vector into two subvectors x and x with N - 1 and M - N + 1 elements respectively. After some simple transformations, an outer layer Sphere Decoder (SD) can be used to choose proper x and then use an inner layer SD to decide x, thus the whole transmitted signal vector is obtained. Simulation results show that Double Layer Sphere Decoding (DLSD) has far less complexity than the existing Generalized Sphere Decoding (GSDs).