Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

The Properties of Voronoi Tessellations for 2D Ising Fluid Near the Second-Order Magnetic Phase Transition

The topological and geometrical properties of Voronoi cells generated for 2D fluid of hard disks with Ising-like spins near the second-order phase transition from the paramagnetic to ferromagnetic phase are described for different disk densities. The comparison with Voronoi cells generated for the random hard disk system is given.     

The effects of flow split ratio and flow rate in manifolds

This paper reports a combined experimental and numerical investigation of three-dimensional steady turbulent flows in inlet manifolds of square cross-section. Predictions and measurements of the flows were carried out using computational fluid dynamics and laser Doppler anemometry techniques respectively. The flow structure was characterized in detail and the effects of flow split ratio and inlet flow rate were studied. These were found to cause significant variations in the size and shape of recirculation regions in the branches, and in the turbulence levels. It was then found that there is a significant difference between the flow rates through different branches. The performance of the code was assessed through a comparison between predictions and measurements. The comparison demonstrates that all important features of the flow are well represented by the predictions.     

Decay of confined,two-dimensional,spatially periodic arrays of vortices: A numerical investigation

The disarrangement of a perturbed lattice of vortices was studied numerically. The basic state is an exponentially decaying, exact solution of the Navier-Stokes equations. Square arrays of vortices with even numbers of vortex cells along each side were perturbed and their evolution was investigated. Whether the energy in the perturbation grows somewhat before it decays or decays monotonically depends on the initial strength of the vortices of the basic state, the extent of lateral confinement and the structure of the perturbation. The critical condition for temporally local instability, i.e. the critical amplitude of the basic state that must be exceeded to allow energy transfer from the basic state to the perturbation, is discussed. In the strongly confined case of a square lattice of four vortices the appearance of enchancement of global rotation is the result of energy transfer from the basic state to a temporally local unstable mode. Energy is transferred from the basic state to larger-scaled structures (inverse cascade) only if the scales of the larger structures are inherently contained in the initial structure of the perturbation. The initial structure of the double array of vortices is not maintained except for a very special form of perturbation. The facts that large scales decay more slowly than small scales and that, when non-linearities are sufficiently strong, energy is transferred from one scale to another explain the differences in the disarrangement process for different initial strengths of the vortices of the basic state. The stronger vortices, i.e. the vortices perturbed in a manner that increases their strength, tend to dominate the weaker vortices. The pairing and subsequent merging (or capture) of vortices of like sense into larger-scale vortices are described in terms of peaks in the evolution of the square root of the palinstrophy divided by the enstrophy.     

The modified characteristic finite difference fractional steps method for the coupled system of fluid dynamics in porous media and its analysis

For the coupled system of multilayer fluid dynamics in porous media, the modified characteristic finite difference fractional steps method applicable to parallel arithmetic is put forward and two‐dimensional and three‐dimensional schemes are used to form a complete set. Some techniques, such as calculus of variations, energy method, piecewise biquadratic interpolation, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates in L² norm are derived to determine the error in the approximate solution. This method has already been applied to the numerical simulation of multilayer fluid dynamics in porous media. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 665–681, 2003.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.     

Extension of multigrid methodology to supersonic/hypersonic 3-D viscous flows

A multigrid acceleration technique developed for solving the three-dimensional Navier–Stokes equations for subsonic/transonic flows has been extended to supersonic/hypersonic flows. An explicit multistage Runge–Kutta type of time-stepping scheme is used as the basic algorithm in conjunction with the multigrid scheme. Solutions have been obtained for a blunt conical frustum at Mach 6 to demonstrate the applicability of the multigrid scheme to high-speed flows. Computations have also been performed for a generic High-Speed Civil Transport configuration designed to cruise at Mach 3. These solutions demonstrate both the efficiency and accuracy of the present scheme for computing high-speed viscous flows over configurations of practical interest.     

Measurement and prediction of solubilities and diffusion coefficients of carbon dioxide in starch–water mixtures at elevated pressures

The solubility and diffusion coefficient of carbon dioxide in intermediate‐moisture starch–water mixtures were determined both experimentally and theoretically at elevated pressures up to 16 MPa at 50 °C. A high‐pressure decay sorption system was assembled to measure the equilibrium CO₂ mass uptake by the starch–water system. The experimentally measured solubilities accounted for the estimated swollen volume by Sanchez–Lacombe equation of state (S‐L EOS) were found to increase almost linearly with pressure, yielding 4.0 g CO₂/g starch–water system at 16 MPa. Moreover, CO₂ solubilities above 5 MPa displayed a solubility increase, which was not contributed by the water fraction in the starch–water mixture. The solubilities, however, showed no dependence on the degree of gelatinization (DG) of starch. The diffusion coefficient of CO₂ was found to increase with concentration of dissolved CO₂, which is pressure‐dependent, and decrease with increasing DG in the range of 50–100%. A free‐volume‐based diffusion model proposed by Areerat was employed to predict the CO₂ diffusivity in terms of pressure, temperature, and the concentration of dissolved CO₂. S‐L EOS was once more used to determine the specific free volume of the mixture system. The predicted diffusion coefficients showed to correlate well with the measured values for all starch–water mixtures. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 607–621, 2006     

Effect of ionic additives on the elution of sodium aryl sulfonates in supercritical fluid chromatography

Addition of a small amount of polar solvent (i.e., modifier) to CO2 in packed column supercritical fluid chromatography (SFC) has shown major improvements in both polar analyte solubility and interaction of the polar analyte with the stationary phase. Recently, the addition of an ionic component (i.e., additive) to the primary modifier by one of us has been shown to extend even further the application of SFC to polar analytes. In this work, the effect of various ionic additives on the elution of ionic compounds, such as sodium 4-dodecylbenzene sulfonate and sodium 4-octylbenene sulfonate, has been studied. The additives were lithium acetate, ammonium acetate, tetramethylammonium acetate, tetrabutylammonium acetate, and ammonium chloride dissolved in methanol. Three stationary phases with different degrees of deactivation were considered: conventional cyanopropyl, deltabond cyanopropyl, and bare silica. The effect of additive concentration and additive functionality on analyte retention was investigated. Sodium 4-dodecylbenzene sulfonate was successfully eluted using all the additives with good peak shape under isocratic/isobaric/isothermal conditions. Different additives, however, yielded different retention times and in some cases different peak shapes.