António Tadeu, Andreia Pereira, Luís Godinho¹

CMES-Computer Modeling in Engineering & Sciences, Vol.2, No.1, pp. 49-62, 2001, DOI:10.3970/cmes.2001.002.049

Abstract The Boundary Element Method (BEM) is used to compute the three-dimensional variation pressure field generated by a point pressure source inside a flat waveguide channel filled with a homogeneous fluid, in the presence of infinite rigid circular pipelines. The problem is solved in the frequency domain, using boundary elements to model the pipeline and an appropriate Green's function to simulate the free surface and the rigid floor of the channel. Because of the 2 ---1/2 ---D geometry of the problem, the separation of variables has been used, and the solution at each frequency is expressed in terms of waves with… More >

Paulo Santos, Julieta António, António Tadeu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.4, pp. 79-98, 2000, DOI:10.3970/cmes.2000.001.531

Abstract This paper presents the three-dimensional scattering field obtained when 2D smooth topographical deformations are subjected to a dilatational point load placed at some point in the medium. The solution is formulated using boundary elements for a wide range of frequencies and spatially harmonic line loads, which are then used to obtain time series by means of (fast) inverse Fourier transforms into space-time. The topographical surface is modeled with a number of boundary elements, defined according to the excitation frequency of the harmonic source, and in such a way that the free surface can be discretized along a sufficient distance to… More >

J.A.F. Santiago¹, L.C. Wrobel²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.3, pp. 73-80, 2000, DOI:10.3970/cmes.2000.001.375

Abstract This work presents a boundary element formulation for two-dimensional acoustic wave propagation in shallow water. It is assumed that the velocity of sound in water is constant, the free surface is horizontal, and the seabed is irregular. The boundary conditions of the problem are that the sea bottom is rigid and the free surface pressure is atmospheric.
For regions of constant depth, fundamental solutions in the form of infinite series can be employed in order to avoid the discretisation of both the free surface and bottom boundaries. When the seabed topography is irregular, it is necessary to divide the… More >

Igor Patlashenko¹, Dan Givoli²

CMES-Computer Modeling in Engineering & Sciences, Vol.1, No.2, pp. 61-70, 2000, DOI:10.3970/cmes.2000.001.221

Abstract The method of Absorbing Boundary Conditions (ABCs) is considered for the numerical solution of a class of nonlinear exterior wave scattering problems. Recently, a scheme based on the exact nonlocal Dirichlet-to-Neumann (DtN) ABC has been proposed for such problems. Although this method is very accurate, it is also highly expensive computationally. In this paper, the nonlocal ABC is replaced by a low-order local ABC, which is obtained by localizing the DtN condition in a certain "optimal'' way. The performance of the new local scheme is compared to that of the nonlocal scheme via numerical experiments in two dimensions. More >

C. S. Iorio¹, C. Perfetti¹

FDMP-Fluid Dynamics & Materials Processing, Vol.11, No.1, pp. 27-48, 2015, DOI:10.3970/fdmp.2015.011.027

Abstract In the past few years, modeling of the Acoustic StandingWaves (ASW) phenomena has become a topic of great interest due to its theoretical connections with particle/cells manipulation techniques, which represent important tools in the biotechnology field. The present paper proposes a model based on the use of moving wall boundary conditions coupled with a viscous compressible fluid in a square channel. This model successfully achieved the generation of ASWs in the square cross-section for several resonance frequencies; the corresponding acoustic potential for the fundamental resonant mode and several harmonics have also been calculated and are discussed here. More >

R. Prosser¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.4, pp. 411-424, 2009, DOI:10.3970/fdmp.2009.005.411

Abstract A novel wavelet-based method for the simulation of reacting flows on adaptive meshes is presented. The method is based on a subtraction algorithm, wherein the wavelet coefficients are calculated from the low resolution up (as opposed to the standard top-down approach). The advantage of this new method is that it allows the calculation of wavelet coefficients on sparse grids, and thus lends itself more readily to adaptive computational meshes than does the traditional wavelet algorithm. The approach is used to simulate a one-dimensional laminar pre-mixed flame with different Lewis numbers. The computational grid is adapted via the removal of grid… More >

Davood Kalantari¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.1, pp. 81-92, 2009, DOI:10.3970/fdmp.2009.005.081

Abstract In this paper a theoretical approach is elaborated for modelling the impact and ensuing spreading behaviour of a liquid droplet after its collision with a flat and rigid surface. The major outcomes of such a study can be summarized as follows: 1) The propagating-shock-wave velocity associated with the droplet is not a constant value but depends on the impact velocity and the physical and geometrical properties of the droplet. 2) The initial radial ejection velocity of the lamella is proportional to the shock-wave velocity (ua) and the impact velocity (0) according to the expression (a-u0)1/2. 3) The deceleration behaviour of… More >

Huawei Han¹, Kamiel S. Gabriel²

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.4, pp. 287-298, 2006, DOI:10.3970/fdmp.2006.002.287

Abstract Data on air-water co-current two-phase annular flow in a tube with an inner diameter of 9.525 mm (3/8 in) were previously collected at both microgravity u-g and normal gravity (1-g) conditions. The data contained measurements of pressure drop, in addition to previously published data of liquid film thickness. This paper presents the results and analysis of the influence of flow pressure gradient on interfacial wave properties of annular flow at both microgravity and normal gravity. The examined wave properties include wave base thickness, wave height (or roughness height), wave spacing, wave speed and wave frequency. It was found that, the… More >

D.E. Melnikov¹, V.M. Shevtsova²

FDMP-Fluid Dynamics & Materials Processing, Vol.3, No.4, pp. 329-338, 2007, DOI:10.3970/fdmp.2007.003.329

Abstract Thermocapillary convection in a long vertical liquid column (called liquid bridge) subjected to heating from above is considered for a three-dimensional Boussinesq fluid. The problem is solved numerically via finite-volume method. Full system of three dimensional Navier-Stokes equations coupled with the energy equation is solved for an incompressible fluid. Instability sets in through a wave propagating in axial direction with zero azimuthal wave number, which is a unique stable solution over a wide range of supercritical heating. Further increasing the applied temperature difference results in bifurcation of a second wave traveling azimuthally with a slightly higher frequency. The two waves… More >

M.K. Achour¹, S. Kaddeche², A. Gharbi², H. Ben Hadid³, D. Henry³

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.4, pp. 277-284, 2005, DOI:10.3970/fdmp.2005.001.277

Abstract The effects of an acoustic wave on the instabilities occurring in a lateral differentially heated cavity are investigated numerically. Linear stability results show that the acoustic wave affects significantly the instability characteristics of such a Hadley flow. Indeed, the sound field is found to stabilize both two dimensional transverse stationary and three dimensional longitudinal oscillatory instabilities which are the most critical modes affecting the buoyant convection in the fluid layer. Nevertheless, when stabilized by an acoustic wave, the 2D modes turn from stationary to oscillatory, with the known consequences of such a change on mass and heat transfer, especially in… More >