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.     

Operator‐splitting method for the analysis of cavitation in liquid‐lubricated herringbone grooved journal bearings

This paper presents an operator‐splitting method (OSM) for the solution of the universal Reynolds equation. Jakobsson–Floberg–Olsson (JFO) pressure conditions are used to study cavitation in liquid‐lubricated journal bearings. The shear flow component of the oil film is first solved by a modified upwind finite difference method. The solution of the pressure gradient flow component is computed by the Galerkin finite element method. Present OSM solutions for slider bearings are in good agreement with available analytical and experimental results. OSM is then applied to herringbone grooved journal bearings. The film pressure, cavitation areas, load capacity and attitude angle are obtained with JFO pressure conditions. The calculated load capacities are in agreement with available experimental data. However, a detailed comparison of the present results with those predicted using Reynolds pressure conditions shows some differences. The numerical results showed that the load capacity and the critical mass of the journal (linear stability indicator) are higher and the attitude angle is lower than those predicted by Reynolds pressure conditions for cases of high eccentricities. 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.     

Grid‐free surface vorticity method applied to flow induced vibration of flexible cylinders

In order to study cross flow induced vibration of heat exchanger tube bundles, a new fluid–structure interaction model based on surface vorticity method is proposed. With this model, the vibration of a flexible cylinder is simulated at Re=2.67 × 10⁴, the computational results of the cylinder response, the fluid force, the vibration frequency, and the vorticity map are presented. The numerical results reproduce the amplitude‐limiting and non‐linear (lock‐in) characteristics of flow‐induced vibration. The maximum vibration amplitude as well as its corresponding lock‐in frequency is in good agreement with experimental results. The amplitude of vibration can be as high as 0.88D for the case investigated. As vibration amplitude increases, the amplitude of the lift force also increases. With enhancement of vibration amplitude, the vortex pattern in the near wake changes significantly. This fluid–structure interaction model is further applied to simulate flow‐induced vibration of two tandem cylinders and two side‐by‐side cylinders at similar Reynolds number. Promising and reasonable results and predictions are obtained. It is hopeful that with this relatively simple and computer time saving method, flow induced vibration of a large number of flexible tube bundles can be successfully simulated. Copyright © 2004 John Wiley & Sons, Ltd.     

Damage and fracture evaluation of granular composite materials by digital image correlation method

This paper presents the applications of digital image correlation technique to the mesoscopic damage and fracture study of some granular based composite materials including steelfiber reinforced concrete, sandstone and crystal-polymer composite. The deformation fields of the composite materials resulted from stress localization were obtained by the correlation computation of the surface images with loading steps and thus the related damage prediction and fracture parameters were evaluated. The correlation searching could be performed either directly based on the gray levels of the digital images or from the wavelet transform (WT) coefficients of the transform spectrum. The latter was developed by the authors and showed higher resolution and sensitivity to the singularity detection. Because the displacement components came from the rough surfaces of the composite materials without any coats of gratings or fringes of optical interferometry, both surface profiles and the deformation fields of the composites were visualized which was helpful to compare each other to analyze the damage of those heterogeneous materials. The project supported by the National Natural Science Foundation of China (10125211 and 10072002), the Scientific Committee of Yunnan Province for the Program of Steel Fiber Reinforced Concrete, and the Institute of Chemical Materials, CAEP at Mianyang     

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.     

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     

Rod–coil block copolymers: An iterative synthetic approach via living free-radical procedures

A modular approach has been developed for the synthesis of rod–coil block copolymers involving the initial preparation of a macroinitiator based on the rod block followed by the growth of the coil segment with living free-radical procedures. The key feature of this strategy is the utilization of an alkoxyamine group from the beginning of the synthesis, which serves as a solubilizing group and ensures that each rod block contains a single initiating fragment. Using this approach permits block copolymers based on insoluble biphenyl ester oligomers to be conveniently prepared with coil segments that range from styrenes to acrylates to 1,3-dienes. The tendency of the rod segments to crystallize is strongly dependent on the weight fraction of the rod segment and the chemical nature of the coil segment. Rod–coil molecules containing at least 25–35 wt % polystyrene or poly(n-butyl acrylate) coil segments show a two-dimensional hexagonal arrangement of rod aggregates, as characterized by transmission electron microscopy and small-angle X-ray scattering. Polyisoprene block copolymers exhibit a lamellar microstructure with short rigid domains in which the rod units lie in an interdigitated smectic C arrangement. © 2003 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 41: 3640–3656, 2003     

Isospectral theory of Euler equations

Isospectral problem of both 2D and 3D Euler equations of inviscid fluids, is investigated. Connections with the Clay problem are described. Spectral theorem of the Lax pair is studied.     

Effects of the H-bond donor structure on the properties of supramolecular side chain liquid crystalline polymers based on poly(3-carboxypropylmethylsiloxane-co-dimethylsiloxane)s and stilbazoles

Supramolecular side chain liquid crystalline polymers (SCLCPs) based on poly(3-carboxypropylmethylsiloxane-co-dimethylsiloxane) (PSIX, X=100, 76, 60, 41 or 23, denoting the mole percentage of 3-carboxypropylmethylsiloxane unit in the polymer) and stilbazole derivatives have been obtained through intermolecular hydrogen bonding (H-bonding) interactions between the carboxylic acid and the pyridyl moieties. The formation of H-bonding and self-assembly results in the formation of new mesogenic units, in which H-bonds function as molecular connectors. FTIR shows the existence of H-bonding in the complexes. The polymeric complexes behave as single component liquid crystalline polymers and exhibit stable and enantiotropic mesophases. The liquid crystalline properties of the supramolecular SCLCPs were studied using differential scanning calorimetry, polarizing optical microscopy and X-ray diffraction, and were found to exhibit smectic A phases with focal-conic textures. The thermal stability of the SCLCP increases on increasing the carboxylic acid content in the polysiloxane and the concentration of the stilbazole derivative in the complex. However, the thermal stability decreases on increasing the chain length of the stilbazole derivative. The crystal phase was not formed even on cooling to the glass transition temperature of the polymeric complex.