Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.     

Heat flow disturbed by periodic array of insulated interface cracks

The thermal effects of an interface crack between two dissimilar half-spaces is considered. The interface cracks are partially or fully insulated, and spaced in a periodic array. Using the complex variable technique, the temperature and fluxes are found in closed form, and the interactions between heat flows due to nearby cracks are determined.     

Dissipation in Turbulent Shear Flows: A Wall Function Approach for High Reynolds Numbers

Dissipation effects that are neglected in almost all numerical codes when it comes to turbulent flows are analysed systematically. It is shown how the direct and the indirect dissipation can be implemented into the codes when asymptotic considerations about the flow behaviour close to a wall are taken into account. As an example, adiabatic pipe flow is analysed with special emphasis on the distribution of energy between the mechanical and thermal part of the total energy.     

On the efficient application of weighted essentially nonoscillatory scheme

In this paper, the efficient application of high‐order weighted essentially nonoscillatory (WENO) reconstruction to the subsonic and transonic engineering problems is studied. On the basis of the physical considerations, two techniques are proposed to enhance the accuracy and efficiency of the WENO reconstruction. First, it is observed that the WENO scheme using characteristic variable has better accuracy and convergence speed than the scheme using primitive variable. For engineering problems with shock of moderate amplitude, on the basis of the Rankine–Hugoniot conditions, a simplified characteristic‐variable‐based WENO is developed. The simplified version significantly reduces the cost overhead without sacrificing the shock‐capturing capability. Second, in this work, it is found for viscous case that it is better to include the viscous effect. On the basis of a simple analysis, the viscous correction to the parameter ε in the WENO reconstruction is proposed. Numerical results indicate, with the proposed simplified characteristic‐variable‐based reconstruction and the viscous correction, that the nonlinear WENO interpolation is sharply activated in the region of shock jump, whereas in the shockless area, the WENO interpolation weights are tuned towards the designed optimal value for better accuracy. Compared with the original characteristic‐variable‐based WENO, the current implementation has similar accuracy and reduced cost. At the same time, compared with the primitive variable‐based WENO, better accuracy and convergence speed are obtained at marginal cost overhead. Several practical cases are calculated to demonstrate the accuracy and efficiency of the current methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Living fluids

One of the major emerging fields of research of the beginning of this century concerns living fluids. By “living fluids”, we mean two major categories of complex fluids: (i) fluids which are essential to life, like blood, and (ii) active fluids made of particles that are able to propel themselves in the suspending fluid by converting a form of their energy into mechanical motion. Studies on active fluids have known a considerable interest since the last decade. Blood might be viewed as an old topic, but the progresses in experimental techniques, analytical concepts and numerics, have contributed nowadays to a dramatic renewal of the interest in this field, with a great potential towards understanding physical and mechanical factors in cardiovascular diseases. These fields have considerably strengthened interdisciplinary research. The series of reviews of this dossier focus on the tremendous recent progress achieved in research on living fluids both from the experimental and theoretical points of views. These reviews present also the major open issues, making of this dossier a unique guide for future research in these fields. This project grew up thanks to the international summer school that we organized on the topic “living fluids” at the IES (Institut dʼétudes scientifiques) of Cargèse (Corsica) in 2012.     

Chemical analysis of size-tailored magnetic colloids using slurry nebulization in ICP-OES

In this work, inductively coupled plasma emission spectroscopy was utilized to directly determine the chemical composition of magnetic fluids constituted of size-sorted ferrite nanoparticles in aqueous solution. Nickel and cobalt nanoferrites were chemically synthesized following a bottom-up route and dispersed under various pH conditions. Size and structural characteristics of nanograins were investigated by X-ray diffraction using a synchrotron source. Chemical analysis was then carried out by directly introducing diluted magnetic fluid samples (slurries) into the spectrometer. To achieve reliable measurements, sample conditions and apparatus parameters were carefully investigated. Slurry stability must be optimized in order to obtain reproducible and accurate analysis. The instrument must also be calibrated to minimize the difference between the signal produced by slurries and that of aqueous ordinary solutions. Furthermore, slurry sample introduction offers many advantages over conventional sample digestion, including reduced sample pretreatment time, less possibility of contamination and the use of direct calibration with aqueous solutions.     

Numerical simulation of 2D lid‐driven cavity flow with CLEARER algorithm on extremely highly skewed grids at high Reynolds numbers

It is crucial to deal with the grid non‐orthogonality effectively in solving the flow in complex geometries, especially at high Reynolds numbers. In this study, the newly proposed Coupled and Linked Equations Algorithm Revised‐ER (CLEARER) algorithm is adopted to solve this problem successfully. In CLEARER algorithm the second relaxation factor is introduced in constructing the contravariant interface velocities, by setting it to a low value. CLEARER algorithm can overcome the severe grid non‐orthogonality and non‐linearity of equations effectively. After the numerical results with CLEARER are validated with the benchmark solutions, this algorithm is used to solve the lid‐driven flow in inclined cavity with inclination angles varying from 10 to 170^°, and Reynolds numbers varying from 5000 to 15 000. The streamlines and the centerline velocity distributions are provided in detail for all cases, which may offer some guidance for the study in this area. Copyright © 2010 John Wiley & Sons, Ltd.     

The influence of magnetic field swept rate on the ultrasonic propagation velocity of magnetorheological fluids

Structure of magnetorheological (MR) fluids depends on the strength of the magnetic field applied and on the mode of its application. The ultrasonic wave propagation velocity changes under the effect of an external magnetic field as a result of formation of clusters arranged along the direction of the field in the MR fluids. Therefore, we propose a qualitative analysis of these clustering structures by measuring properties of ultrasonic propagation. Since the MR fluids are opaque, the non-contact inspection using this ultrasonic technique can be very useful. In this study, we measured ultrasonic propagation velocity in MR fluid influenced by an external magnetic field for different swept rate precisely. With increasing magnetic field intensity, the changes of the ultrasonic wave velocity are more pronounced. Sedimentation effect takes place in certain time for different swept rate due to magnetic particle size and it follows linear relationship in log scale. Significant differences of the ultrasonic wave velocity are established between the case when the field is swept at a constant rate and the case when it is stepped up.     

Three‐dimensional computation for flow‐induced vibrations in in‐line and cross‐flow directions of a circular cylinder

We present numerical results for in‐line and cross‐flow vibrations of a circular cylinder, which is immersed in a uniform flow and is elastically supported by damper‐spring systems to compute vibrations of a rigid cylinder. In the case of a circular cylinder with a low Scruton number, it is well‐known that two types of self‐excited vibrations appear in the in‐line direction in the range of low reduced velocities. On the other hand, a cross‐flow vibration of the circular cylinder can be excited in the range of high reduced velocities. Therefore, we compute the flow‐induced vibrations of the circular cylinder in the wide range of the reduced velocities at low and high Scruton numbers and discuss about excitation mechanisms in the in‐line and cross‐flow directions. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of element order and interface reconstruction in FEM/volume‐of‐fluid incompressible flow simulation

In previous studies, the moment‐of‐fluid interface reconstruction method showed dramatic accuracy improvements in static and pure advection tests over existing methods, but this did not translate into an equivalent improvement in volume‐tracked multimaterial incompressible flow simulation using low‐order finite elements. In this work, the combined effects of the spatial discretization and interface reconstruction in flow simulation are examined. The mixed finite element pairs, Q₁Q₀ (with pressure stabilization) and Q₂P _{? 1} are compared. Material order‐dependent and material order‐independent first and second‐order accurate interface reconstruction methods are used. The Q₂P _{? 1} elements show significant improvements in computed flow solution accuracy for single material flows but show reduced convergence using element‐average piecewise constant density and viscosity in volume‐tracked simulations. In general, a refined Q₁Q₀ grid, with better material interface resolution, provided an accuracy similar to the Q₂P _{? 1} element grid with a comparable number of degrees of freedom. Moment‐of‐fluid shows more benefit from the higher‐order accurate flow simulation than the LVIRA, Youngs', and power diagram interface reconstruction methods, especially on unstructured grids, but does not recover the dramatic accuracy improvements it has shown in advection tests. Published 2012. This article is a US Government work and is in the public domain in the USA.