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.     

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.     

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     

Apparatus for High-Resolution Solid Hydrogen Matrix Isolation Spectroscopy and Its Application

Apparatus integrating a Fourier transform-infrared (FT-IR) spectrometer and a mid-infrared difference frequency generation (DFG) laser spectrometer was built for the study of the matrix isolation spectrum in solid molecular hydrogen. A 3-cm-long molecular hydrogen crystal was grown in a liquid-helium Dewar, and its infrared absorption spectrum in the 1-5 μm region was recorded to test the system. The W₀(0) (ν=0←0, J=6←0) line around 2410 cm⁻¹ of solid hydrogen was investigated with the DFG laser spectroscopy. High-resolution matrix isolation spectrum of CO² co-deposited with hydrogen on a BaF₂ cold plate at liquid-helium temperature was studied.     

A retractable structure based on Bricard linkages and rotating rings of tetrahedra

Retractable structures which alter their geometries according to practical requirements are widely used in roof structures of stadiums for their versatility. A retractable structure based on threefold-symmetric Bricard linkages and rotating rings of tetrahedra is proposed and developed in this paper. By replacing each link of a threefold-symmetric Bricard linkage with a tetrahedron, a retractable structure is obtained which can repeatedly open and close by rotating the tetrahedra. The geometric relationship between the retractable structure and the threefold-symmetric Bricard linkage is derived and ranges of main geometric parameters are determined to ensure a continuous and smooth movement in the deployment of the structure. Then the relative displacement of the support and the open rate of the structure, two main parameters concerned in practical applications, are investigated. Preferred ranges of main geometric parameters are proposed and a physical model is manufactured to verify them. A discussion on some possible improvements and modifications of the structure is also presented.     

Axisymmetric frictionless contact of functionally graded materials

The main interest of this study is a new method to solve the axisymmetric frictionless contact problem of functionally graded materials (FGMs). Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into a series of sub-layers with shear modulus varying linearly in each sub-layer and continuous at the sub-interfaces. With this model, the axisymmetric frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Hankel integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically. An erratum to this article can be found at     

G(F4′)条件与加强形式的停止定理

首次提出了Ｇ（Ｆ４′）条件，给出了Ｇ（Ｆ４′）条件的两个等价关系，得到了加强形式的停止定理和强鞅的一个等价关系．     

Improvement of least-squares finite difference method on simulating slider air bearings in head disk systems

Least-squares finite difference (LSFD) method, one of mesh-free methods, is used to solve slider air bearings problem through discritizing the generalized Reynolds equation into nonlinear systems of algebraic equations. Two approximation schemes for the linearization of these equations are presented and compared. And, some new techniques to search supporting points for the reference node in the mesh-free method were proposed and explored. Therefore, these improvements eliminate some potential limitation of the LSFD method previously published and further facilitate its employment in complex slider models. Advanced step slider as an example of negative pressure sliders is simulated and verified using the improved LSFD mesh-free method in head disk systems.     

Fiber formation via solution spinning of the cellulose/ammonia/ammonium thiocyanate system

Fiber formation via the cellulose/ammonia/ammonium thiocyanate system by wet spinning has been investigated. This report presents a characterization of the structure and tensile properties of fibers spun under various coagulation conditions. Microscopic observations showed that the molecular size of coagulant was the dominant factor governing the crosssectional shape of the fibers. Density, birefringence, and crystallinity data indicated that a higher cellulose concentration and lower coagulation temperature favored development of a fiber with a denser and more oriented structure. Under optimum conditions, a welldefined fibrillar structure was obtained. Fiber tensile property measurements suggested the existence of a linear relationship between the fiber breaking tenacity and the product of the square of the Hermans' orientation factor and the infrared crystallinity index.