A projection method for solving incompressible viscous flows on domains with moving boundaries

In this paper, a projection method is presented for solving the flow problems in domains with moving boundaries. In order to track the movement of the domain boundaries, arbitrary‐Lagrangian–Eulerian (ALE) co‐ordinates are used. The unsteady incompressible Navier–Stokes equations on the ALE co‐ordinates are solved by using a projection method developed in this paper. This projection method is based on the Bell's Godunov‐projection method. However, substantial changes are made so that this algorithm is capable of solving the ALE form of incompressible Navier–Stokes equations. Multi‐block structured grids are used to discretize the flow domains. The grid velocity is not explicitly computed; instead the volume change is used to account for the effect of grid movement. A new method is also proposed to compute the freestream capturing metrics so that the geometric conservation law (GCL) can be satisfied exactly in this algorithm. This projection method is also parallelized so that the state of the art high performance computers can be used to match the computation cost associated with the moving grid calculations. Several test cases are solved to verify the performance of this moving‐grid projection method. Copyright © 2004 John Wiley Sons, Ltd.     

Numerical analysis of the rotating viscous flow approaching a solid sphere

A numerical simulation is performed to investigate the flow induced by a sphere moving along the axis of a rotating cylindrical container filled with the viscous fluid. Three‐dimensional incompressible Navier–Stokes equations are solved using a finite element method. The objective of this study is to examine the feature of waves generated by the Coriolis force at moderate Rossby numbers and that to what extent the Taylor–Proudman theorem is valid for the viscous rotating flow at small Rossby number and large Reynolds number. Calculations have been undertaken at the Rossby numbers (R_o) of 1 and 0.02 and the Reynolds numbers (R_e) of 200 and 500. When R_o=O(1), inertia waves are exhibited in the rotating flow past a sphere. The effects of the Reynolds number and the ratio of the radius of the sphere and that of the rotating cylinder on the flow structure are examined. When R_o ? 1, as predicted by the Taylor–Proudman theorem for inviscid flow, the so‐called ‘Taylor column’ is also generated in the viscous fluid flow after an evolutionary course of vortical flow structures. The initial evolution and final formation of the ‘Taylor column’ are exhibited. According to the present calculation, it has been verified that major theoretical statement about the rotating flow of the inviscid fluid may still approximately predict the rotating flow structure of the viscous fluid in a certain regime of the Reynolds number. Copyright © 2004 John Wiley & Sons, Ltd.     

A co‐ordinate system independent streamwise upwind method for fluid flow computation

A new co‐ordinate invariant streamwise upwind formulation for convection dominated flows is developed. The eddy diffusivity/viscosity is added directly to the equations in order to remove the oscillations in the solution. The equations then can be solved by any high‐order scheme and the solution retains the accuracy of the high‐order scheme. The accuracy and reduced lateral thickness growth rate are demonstrated with several numerical examples, including pure convective flows and lid‐driven cavity flow. The lateral spreading due to the numerical diffusion is controlled by the anisotropic tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

Passive control of the flow around a square cylinder using porous media

The passive control of bluff body flows using porous media is investigated by means of the penalization method. This method is used to create intermediate porous media between solid obstacles and the fluid in order to modify the boundary layer behaviour. The study covers a wide range of two‐dimensional flows from low transitional flow to fully established turbulence by direct numerical simulation of incompressible Navier–Stokes equations. A parametric study is performed to illustrate the effect of the porous layer permeability and thickness on the passive control. The numerical results reveal the ability of porous media to both regularize the flow and to reduce the drag forces up to 30%. Copyright © 2004 John Wiley & Sons, Ltd.     

Transient free‐surface flow of a viscoelastic fluid in a narrow channel

The interplay between inertia and elasticity is examined for transient free‐surface flow inside a narrow channel. The lubrication theory is extended for the flow of viscoelastic fluids of the Oldroyd‐B type (consisting of a Newtonian solvent and a polymeric solute). While the general formulation accounts for non‐linearities stemming from inertia effects in the momentum conservation equation, and the upper‐convected terms in the constitutive equation, only the front movement contributes to non‐linear coupling for a flow inside a straight channel. In this case, it is possible to implement a spectral representation in the depthwise direction for the velocity and stress. The evolution of the flow field is obtained locally, but the front movement is captured only in the mean sense. The influence of inertia, elasticity and viscosity ratio is examined for pressure‐induced flow. The front appears to progress monotonically with time. However, the velocity and stress exhibit typically a strong overshoot upon inception, accompanied by a plug‐flow behaviour in the channel core. The flow intensity eventually diminishes with time, tending asymptotically to Poiseuille conditions. For highly elastic liquids the front movement becomes oscillatory, experiencing strong deceleration periodically. A multiple‐scale solution is obtained for fluids with no inertia and small elasticity. Comparison with the exact (numerical) solution indicates a wide range of validity for the analytical result. Copyright © 2004 John Wiley & Sons, Ltd.     

Crystallographic report: Crystal structure of tetra(4‐methyl‐5‐imidazole‐carboxyaldehyde)zinc(II) diperchlorate

The central zinc(II) atom in the title complex is tetrahedrally coordinated by four nitrogen atoms derived from 4‐methyl‐5‐imidazolecarboxyaldehyde ligands with Zn? N in the range 2.007(3) to 2.026(4) Å. Copyright © 2003 John Wiley & Sons, Ltd.     

Crystallographic report: Crystal structure of a one‐dimensional zinc(II) coordination polymer containing 4,4′‐biphenyldicarboxylate hemihydrate

A one‐dimensional zinc(II) coordination polymer has been constructed from zinc(II), 4,4′‐biphenyldicarboxylate and pyridine in which each zinc(II) atom is coordinated by two pyridine ligands and two monodentate 4,4′‐biphenyldicarboxylate ligands that define a distorted tetrahedral geometry. Copyright © 2003 John Wiley & Sons, Ltd.     

Preparation,thermal properties,and flame retardance of epoxy–silica hybrid resins

A novel epoxy system was developed through the in situ curing of bisphenol A type epoxy and 4,4′‐diaminodiphenylmethane with the sol–gel reaction of a phosphorus‐containing trimethoxysilane (DOPO–GPTMS), which was prepared from the reaction of 9,10‐dihydro‐9‐oxa‐10‐phosphaphenanthrene‐10‐oxide (DOPO) with 3‐glycidoxypropyltrimethoxysilane (GPTMS). The preparation of DOPO–GPTMS was confirmed with Fourier transform infrared, ¹H and ³¹P NMR, and elemental analysis. The resulting organic–inorganic hybrid epoxy resins exhibited a high glass‐transition temperature (167 °C), good thermal stability over 320 °C, and a high limited oxygen index of 28.5. The synergism of phosphorus and silicon on flame retardance was observed. Moreover, the kinetics of the thermal oxidative degradation of the hybrid epoxy resins were studied. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 2354–2367, 2003     

Synthesis and characterization of wholly aromatic poly(azomethine)s containing donor–acceptor triphenylamine moieties

N‐(4‐nitrophenyl)‐4′,4″‐bisformyl‐diphenylamine was synthesized from N‐(4‐nitrophenyl)‐diphenylamine by the Vilsmeier‐Haack reaction. Soluble aromatic poly(azomethine)s (PAMs) were prepared by the solution polycondensation of N‐(4‐nitrophenyl)‐4′,4″‐bisformyl‐diphenylamine and aromatic diamine in N‐methyl‐2‐pyrrolidone (NMP) at room temperature under reduced pressure. All the PAMs are highly soluble in various organic solvents, such as N,N‐dimethylacetamide (DMAc), chloroform (CHCl₃), and tetrahydrofuran (THF). Differential scanning calorimetry (DSC) indicated that these PAMs had glass‐transition temperatures (T_gs) in the range of 170–230 °C, and a 10% weight‐loss temperatures in excess of 490 °C with char yield at 800 °C in nitrogen higher than 60%. These PAMs in NMP solution showed UV‐Vis charge‐transfer (CT) absorption at 405–421 nm and photoluminescence peaks around 462–466 nm with fluorescence quantum efficiency (Φ_F) 0.10–0.99%. The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energy levels of these PAMs can be determined from cyclic voltammograms as 4.86–5.43 and 3.31–3.34 eV, respectively. © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 4921–4932, 2007