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.     

Stress concentration and effective stiffness of aligned fiber reinforced composite with anisotropic constituents

The paper addresses the problem of calculation of the local stress field and effective elastic properties of a unidirectional fiber reinforced composite with anisotropic constituents. For this aim, the representative unit cell approach has been utilized. The micro geometry of the composite is modeled by a periodic structure with a unit cell containing multiple circular fibers. The number of fibers is sufficient to account for the micro structure statistics of composite. A new method based on the multipole expansion technique is developed to obtain the exact series solution for the micro stress field. The method combines the principle of superposition, technique of complex potentials and some new results in the theory of special functions. A proper choice of potentials and new results for their series expansions allow one to reduce the boundary-value problem for the multiple-connected domain to an ordinary, well-posed set of linear algebraic equations. This reduction provides high numerical efficiency of the developed method. Exact expressions for the components of the effective stiffness tensor have been obtained by analytical averaging of the strain and stress fields.     

Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder

This paper presents a numerical study of a two‐dimensional time‐dependent flow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time. In addition, the accuracy of these values obtained with different time stepping schemes and different finite element methods is studied. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order accurate numerical solutions of incompressible flows with the artificial compressibility method

A high‐order accurate, finite‐difference method for the numerical solution of incompressible flows is presented. This method is based on the artificial compressibility formulation of the incompressible Navier–Stokes equations. Fourth‐ or sixth‐order accurate discretizations of the metric terms and the convective fluxes are obtained using compact, centred schemes. The viscous terms are also discretized using fourth‐order accurate, centred finite differences. Implicit time marching is performed for both steady‐state and time‐accurate numerical solutions. High‐order, spectral‐type, low‐pass, compact filters are used to regularize the numerical solution and remove spurious modes arising from unresolved scales, non‐linearities, and inaccuracies in the application of boundary conditions. The accuracy and efficiency of the proposed method is demonstrated for test problems. Copyright © 2004 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

Standardization of BSE rapid test performances and experiences gathered during the implementation of large-scale testing

The accumulation of pathological prion protein is used as a diagnostic marker for transmissible spongiform encephalopathies. According to European Union (EU) regulations cattle older than 30 months of age (Germany, France, Italy, and Spain by national law >24 months) and slaughtered for human consumption must be tested by using rapid tests for bovine spongiform encephalopathy (BSE). Likewise fallen stock and clinically affected animals must be tested. This article gives a short overview of the incidence of BSE in Europe. The diagnostic hierarchy, i.e., the officially approved methodology for the confirmation of suspect rapid test cases, and the organization of the numerous laboratories involved in this large-scale testing for BSE are described. Special emphasis is given to necessary quality control measures currently in place for BSE rapid testing laboratories and to measures intended to assure a consistent performance of the commercially available rapid test kits.Presented at BERM-9—Ninth International Symposium on Biological and Environmental Reference Materials, 15–19 June 2003, Berlin, Germany.     

Standardization of Methods for Characterizing the Surface Geometry of Solids

Since a comprehensive survey published in 1999 [1] much work was done in standardizing measuring methods to characterize the surface geometry of dispersed and/or porous solids and to certify reference materials. The present paper is an extension of a short communication [2]. It gives a survey on existing standards and reports on new drafts and proposals.     

Numerical method for calculation of the incompressible flow in general curvilinear co‐ordinates with double staggered grid

A solution methodology has been developed for incompressible flow in general curvilinear co‐ordinates. Two staggered grids are used to discretize the physical domain. The first grid is a MAC quadrilateral mesh with pressure arranged at the centre and the Cartesian velocity components located at the middle of the sides of the mesh. The second grid is so displaced that its corners correspond to the centre of the first grid. In the second grid the pressure is placed at the corner of the first grid. The discretized mass and momentum conservation equations are derived on a control volume. The two pressure grid functions are coupled explicitly through the boundary conditions and implicitly through the velocity of the field. The introduction of these two grid functions avoids an averaging of pressure and velocity components when calculating terms that are generated in general curvilinear co‐ordinates. The SIMPLE calculation procedure is extended to the present curvilinear co‐ordinates with double grids. Application of the methodology is illustrated by calculation of well‐known external and internal problems: viscous flow over a circular cylinder, with Reynolds numbers ranging from 10 to 40, and lid‐driven flow in a cavity with inclined walls are examined. The numerical results are in close agreement with experimental results and other numerical data. Copyright © 2003 John Wiley & Sons, Ltd.     

Predicting the onset of DNA supercoiling using a non-linear hemitropic elastic rod

Elastic rod models provide a means to interpret single molecule DNA experiments as well as predict DNA behavior under physiological conditions. Here we use an elastic rod model to predict the stability boundary (critical torque vs. applied tension) for single molecule DNA experiments in which the molecule is subjected to applied tension and twist. We discuss the shortcomings of the usual isotropic rod model. We then derive a consistent non-linear material law from the general representation for a hemitropic (chiral) rod. Finally, we present results of a standard bifurcation analysis predicting the stability boundary. We find results from the non-linear hemitropic rod to match the data closely.