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.     

HUA’S THEOREM WITH FIVE ALMOSTE QUAL PRIME VARIABLES

§1. IntroductionIn the additive theory of prime numbers, one studies the representation of positive in-tegers by powers of primes. For the quadratic case, Hua [1] proved that each large integercongruent to 5 modulo 24 can be written as the sum of ?ve squ…     

低速二维流动的粒子仿真研究

本文运用信息保存法对低速二维的流动现象进行模拟，考察了低速条件下的有限平板绕流以及微槽道气体流动问题。研究表明：在对低速流动的模拟过程中，运用IP法在能够获得较好的结果的同时，具有比DSMC方法更高的计算效率。     

Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

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     

Water wave radiation by a sphere submerged in water with an ice-cover

Using the multipoles method, we formulate the problems of radiation (both heave and sway) of water waves by a submerged sphere in deep as well as in uniform finite depth water with an ice-cover, with the ice-cover being modelled as an elastic plate of very small thickness. In each case this leads to an infinite system of linear equations which are solved numerically by standard techniques. The added-mass and damping coefficients for a heaving and swaying sphere are obtained and depicted graphically against the wave number for various values of the radius of the submerged sphere and flexural rigidity of the ice-cover to show the effect of the presence of ice-cover on these quantities. When the flexural rigidity is taken to be zero, the numerical results for the added-mass and damping coefficient for water with a free surface are recovered.     

Numerical method and models for anti-plane strain of a system with a thin elastic wedge

The paper presents a general method to find asymptotics for a (multi-)wedge system containing a thin wedge. It employs separation of the symmetric and anti-symmetric parts of the boundary displacements and tractions of the wedge. The method is applicable when the angle of the thin wedge turns to zero. A physical interpretation of the derived equations is obtained by using power expansions of non-polynomial functions, which appear after the Mellin transform. We establish that the first term in the expansion of the symmetric part corresponds to shear, while the first term of the anti-symmetric part describes deflection of the wedge axis. Numerical experiments, performed by using a code developed on the basis of the theory, show that using only the first terms of the expansions insignificantly influence accuracy: the approximate results coincide with the exact values of roots to the third significant digit even for the wedge angle of 30°.     

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.     

Diffusion via splitting and remeshing via merging in vortex methods

The technique of splitting a fat vortex element (with a core width larger than some threshold) into some thin ones in order to fix the convergence problem of the core‐spreading vortex methods is convenient and efficient. In particular, it keeps the method purely Lagrangian. In the present investigation, the splitting process is further viewed as part of the physical diffusion process. A new splitting method in which several weaker child vortices surround a thinned but still strong parent vortex is proposed. It is found that because of the survival of the parent vortex, the error arising from the splitting events can be largely reduced. The computational amount on the other hand is kept reasonably large by merging similar and close‐by vortices. The merging scheme designed herein not only involves fewer restrictions but also allows merging vortices of opposite rotations through the viewpoint of remeshing. The validity and accuracy of these techniques, proposed particularly for simulations undergoing lots of splitting and merging events, are verified by successfully simulating the interactions between two Burgers vortices under an external straining field. Copyright © 2005 John Wiley & Sons, Ltd.     

Behaviors of the positive temperature coefficient of resistance of poly(styrene‐co‐n‐butylacrylate) filled with Ni‐plated core‐shell polymeric particles

Conductive composite films of poly(styrene‐co‐n‐butylacrylate) copolymers filled with low‐density, Ni‐plated core‐shell polymeric particles were prepared and their behaviors of positive temperature coefficient of resistance (PTCR) were investigated. When the conductive fillers in the composite film were loaded beyond the critical volume, 10 up to 25 vol %, composite films exhibited a unique electrical resistant transition behavior, which the electrical resistance rapidly increased by several orders of magnitude at the critical temperature. The PTCR transition temperature, in general, occurred before the glass transition temperature of polymer matrix. Further increased the conductive filler loading to 30 vol %, the overpacked conduction paths were formed in the entire composite and the PTCR effects became blurred. While the composite film treated with thermal cycle several times from room temperature up to 120 °C, the electrical resistivity increased accompanied with the shift of the PTCR transition to lower temperature. The reason might have been caused by the formed interfacial cracks within the composite film. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 45: 322–329, 2007