An efficient numerical method for compressible flows of a real gas using arithmetic averaging

An efficient numerical method is presented for the solution of the Euler equations governing the compressible flow of a real gas. The scheme is based on the approximate solution of a specially constructed set of linearised Riemann problems. An average of the flow variables across the interface between cells is required, and this is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual square root averaging. The scheme is applied to a test problem for five different equations of state.

An efficient shock capturing scheme for the pseudo-unsteady equations representing steady inviscid flow

An efficient algorithm is presented for the solution of the steady Euler equations of gas dynamics. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The scheme is applied to a standard test problem of flow down a channel containing a circular arc bump for three different mesh sizes.

An efficient numerical scheme for the Euler equations

A finite difference scheme based on flux difference splitting is presented for the solution of the Euler equations for the compressible flow of an ideal gas. A linearised Riemann problem is defined, and a scheme based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency, leading to arithmetic averaging. This is in contrast to the usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. The scheme is applied to a shock tube problem and a blast wave problem. Each approximate solution compares well with those given by other schemes, and for the shock tube problem is in agreement with the exact solution.

An efficient numerical scheme for the two-dimensional shallow water equations using arithmetic averaging

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations of ideal fluid flow. A linearised problem, analogous to that of Riemann for gas dynamics is defined, and a scheme, based on numerical characteristic decomposition is presented for obtaining approximate solutions to the linearised problem, and incorporates the technique of operator splitting. An average of the flow variables across the interface between cells is required, and this average is chosen to be the arithmetic mean for computational efficiency leading to arithmetic averaging. This is in contrast to usual ‘square root’ averages found in this type of Riemann solver, where the computational expense can be prohibitive. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids nonphysical, spurious oscillations. An extension to the two-dimensional equations with source terms is included. The scheme is applied to the one-dimensional problems of a breaking dam and reflection of a bore, and in each case the approximate solution is compared to the exact solution of ideal fluid flow. The scheme is also applied to a problem of stationary bore generation in a channel of variable cross-section. Finally, the scheme is applied to two other dam-break problems, this time in two dimensions with one having cylindrical symmetry. Each approximate solution compares well with those given by other authors.

An algorithm for the shallow water equations with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations in a generalised coordinate system. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The scheme has good jump capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for flow past a circular obstruction.

An isentropic flow algorithm with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the three-dimensional equations of isentropic flow in a generalised coordinate system, and with a general convex gas law. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The algorithm requires only one function evaluation of the gas law in each computational cell. The scheme has good shock capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for Mach 3 flow of air past a circular cylinder. Furthermore, the algorithm also applies to shallow water flows by employing the familiar gas dynamics analogy.

Prediction of supercritical flow in open channels

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow water equations in open channels. A linearised problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearised problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Cylindrically symmetric flows in a Lagrangian frame

An approximate Riemann solver, in a Lagrangian frame of reference, is presented for the compressible flow equations with cylindrical and spherical symmetry, including flow in a duct of variable cross section. The scheme is applied to a cylindrically symmetric problem involving the interaction of shocks.

Flux difference splitting for open channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.

Strong-segregation limit of the self-consistent field theory for diblock copolymer melts

The self-consistent field theory (SCFT) introduced by Helfand for diblock copolymer melts is expected to converge to the strong-segregation theory (SST) of Semenov in the asymptotic limit, $\chi N \rightarrow \infty$. However, past extrapolations of the lamellar/cylinder and cylinder/sphere phase boundaries, within the standard unit-cell approximation, have cast some doubts on whether or not this is actually true. Here we push the comparison further by extending the SCFT calculations to $\chi N = 512,000$, by accounting for exclusion zones in the coronae of the cylindrical and spherical unit cells, and by examining finite-segregation corrections to SST. In doing so, we provide the first compelling evidence that SCFT does indeed reduce to SST.

Condition number estimates for combined potential integral operators in acoustics and their boundary element discretisation

We consider the classical coupled, combined-field integral equation formulations for time-harmonic acoustic scattering by a sound soft bounded obstacle. In recent work, we have proved lower and upper bounds on the $L^2$ condition numbers for these formulations, and also on the norms of the classical acoustic single- and double-layer potential operators. These bounds to some extent make explicit the dependence of condition numbers on the wave number $k$, the geometry of the scatterer, and the coupling parameter. For example, with the usual choice of coupling parameter they show that, while the condition number grows like $k^{1/3}$ as $k\to\infty$, when the scatterer is a circle or sphere, it can grow as fast as $k^{7/5}$ for a class of `trapping' obstacles. In this paper we prove further bounds, sharpening and extending our previous results. In particular we show that there exist trapping obstacles for which the condition numbers grow as fast as $\exp(\gamma k)$, for some $\gamma>0$, as $k\to\infty$ through some sequence. This result depends on exponential localisation bounds on Laplace eigenfunctions in an ellipse that we prove in the appendix. We also clarify the correct choice of coupling parameter in 2D for low $k$. In the second part of the paper we focus on the boundary element discretisation of these operators. We discuss the extent to which the bounds on the continuous operators are also satisfied by their discrete counterparts and, via numerical experiments, we provide supporting evidence for some of the theoretical results, both quantitative and asymptotic, indicating further which of the upper and lower bounds may be sharper.

Synthesis, structure and magnetic properties of mono- and di-nuclear nickel(II) thiocyanate complexes with tridentate N3 donor Schiff bases

The 1:1 condensation of N-methyl-1,3-diaminopropane and N,N-diethyl-1,2-diminoethane with 2-acetylpyridine, respectively at high dilution gives the tridentate mono-condensed Schiff bases N-methyl-N'-(1-pyridin-2-yl-ethylidene)-propane-1,3-diamine (L-1) and N,N-diethyl-N'-(1-pyridin-2-yl-ethylidene)-ethane-1,2-diamine (L-2). The tridentate ligands were allowed to react with methanol solutions of nickel(II) thiocyanate to prepare the complexes [Ni(L-1)(SCN)(2)(OH2) (1) and [{Ni(L-2)(SCN)}(2)] (2). Single crystal X-ray diffraction was used to confirm the structures of the complexes. The nickel(II) in complex 1 is bonded to three nitrogen donor atoms of the ligand L-1 in a mer orientation, together with two thiocyanates bonded through nitrogen and a water molecule, and it is the first Schiff base complex of nickel(II) containing both thiocyanate and coordinated water. The coordinated water initiates a hydrogen bonded 2D network. In complex 2, the nickel ion occupies a slightly distorted octahedral coordination sphere, being bonded to three nitrogen atoms from the ligand L-2, also in a mer orientation, and two thiocyanate anions through nitrogen. In contrast to 1, the sixth coordination site is occupied by a sulfur atom from a thiocyanate anion in an adjacent molecule, thus creating a centrosymmetric dimer. A variable temperature magnetic study of complex 2 indicates the simultaneous presence of zero-field splitting, weak intramolecular ferromagnetic coupling and intermolecular antiferromagnetic interactions between the nickel(II) centers.

Rheology of magmas with bimodal crystal size and shape distributions: insights from analog experiments

Magmas in volcanic conduits commonly contain microlites in association with preexisting phenocrysts, as often indicated by volcanic rock textures. In this study, we present two different experiments that inves- tigate the flow behavior of these bidisperse systems. In the first experiments, rotational rheometric methods are used to determine the rheology of monodisperse and polydisperse suspensions consisting of smaller, prolate particles (microlites) and larger, equant particles (phenocrysts) in a bubble‐free Newtonian liquid (silicate melt). Our data show that increasing the relative proportion of prolate microlites to equant pheno- crysts in a magma at constant total particle content can increase the relative viscosity by up to three orders of magnitude. Consequently, the rheological effect of particles in magmas cannot be modeled by assuming a monodisperse population of particles. We propose a new model that uses interpolated parameters based on the relative proportions of small and large particles and produces a considerably improved fit to the data than earlier models. In a second series of experiments we investigate the textures produced by shearing bimodal suspensions in gradually solidifying epoxy resin in a concentric cylinder setup. The resulting textures show the prolate particles are aligned with the flow lines and spherical particles are found in well‐organized strings, with sphere‐depleted shear bands in high‐shear regions. These observations may explain the measured variation in the shear thinning and yield stress behavior with increasing solid fraction and particle aspect ratio. The implications for magma flow are discussed, and rheological results and tex- tural observations are compared with observations on natural samples.

Synthesis, X-ray crystal structure and catalytic activity of the manganese(II) complex [Mn(bdoa)(H2O)3](bdoaH2 =benzene-1,2-dioxyacetic acid)

Benzene-1,2-dioxyacetic acid (bdoaH2) reacts with Mn(CH3CO2)2·4H2O in an ethanol-water mixture to give the manganese(II) complex [Mn(bdoa)(H2O)3]. The X-ray crystal structure of the complex shows the metal to be pseudo seven-coordinate. The quadridentate bdoa2− dicar☐ylate ligand forms an essentially planar girdle around the metal, being strongly bondedtransoid by a car☐ylate oxygen atom from each of the two car☐ylate moieties (mean MnO 2.199A˚) and also weakly chelated by the two internal ether oxygen atoms (mean MnO 2.413A˚). The coordination sphere about the manganese is completed by three water molecules (mean MnO 2.146A˚) lying in a meridional plane orthogonal to that of the bdoa2− ligand. Magnetic, conductivity and voltammetry data for the complex are given, and its use as a catalyst for the disproportionisation of H2O2 is described.

Effect of water-soluble polymers, polyethylene glycol and poly(vinylpyrrolidone),on the gelation of aqueous micellar solutions of Pluronic copolymer F127

The micellization of F127 (E98P67E98) in dilute aqueous solutions of polyethylene glycol (PEG6000 and PEG35000) and poly(vinylpyrrolidone) (PVP K30 and PVP K90) is studied. The average hydrodynamic radius (rh,app) obtained from the dynamic light scattering technique increased with increase in PEG concentration but decreased on addition of PVP, results which are consistent with interaction of the micelles with PEG and the formation of micelles clusters, but no such interaction occurs with PVP. Tube inversion was used to determine the onset of gelation. The critical concentration of F127 for gelation increased on addition of PEG and of PVP K30 but decreased on addition of PVP K90. Small-angle X-ray scattering (SAXS) was used to show that the 30 wt% F127 gel structure (fcc) was independent of polymer type and concentration, as was the d-spacing and so the micelle hard-sphere radius. The maximum elastic modulus (G0 max) of 30 wt% F127 decreased from its value for water alone as PEG was added, but was little changed by adding PVP. These results are consistent with the packed-micelles in the 30 wt% F127 gel being effectively isolated from the polymer solution on the microscale while, especially for the PEG, being mixed on the macroscale.