D. Ngo-Cong^1,2, N. Mai-Duy¹, W. Karunasena², T. Tran-Cong^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.3, pp. 311-352, 2012, DOI:10.3970/cmes.2012.083.311

Abstract In this study, local moving least square - one dimensional integrated radial basis function network (LMLS-1D-IRBFN) method is presented and demonstrated with the solution of time-dependent problems such as Burgers' equation, unsteady flow past a square cylinder in a horizontal channel and unsteady flow past a circular cylinder. The present method makes use of the partition of unity concept to combine the moving least square (MLS) and one-dimensional integrated radial basis function network (1D-IRBFN) techniques in a new approach. This approach offers the same order of accuracy as its global counterpart, the 1D-IRBFN method, while the system matrix is more… More >

Peter Lucas¹, Alexander H. van Zuijlen¹, Hester Bijl¹

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.1, pp. 79-106, 2010, DOI:10.3970/cmes.2010.059.079

Abstract Despite the advances in computer power and numerical algorithms over the last decades, solutions to unsteady flow problems remain computing time intensive.
In previous work [Lucas, P.,Bijl, H., and Zuijlen, A.H. van(2010)], we have shown that a Jacobian-free Newton-Krylov (JFNK) algorithm, preconditioned with an approximate factorization of the Jacobian which approximately matches the target residual operator, enables a speed up of a factor of 10 compared to nonlinear multigrid (NMG) for two-dimensional, large Reynolds number, unsteady flow computations. Furthermore, in [Lucas, P., Zuijlen, A.H. van, and Bijl, H. (2010)] we show that this algorithm also greatly outperforms NMG for parameter… More >

R. Avila¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.54, No.1, pp. 15-64, 2009, DOI:10.3970/cmes.2009.054.015

Abstract The Meshless Local Petrov Galerkin (MLPG) method is used to solve the non-steady two dimensional Navier-Stokes equations. Transient laminar flow field calculations have been carried out in domains wherein certain surfaces have: (i) a sliding motion, (ii) a harmonic motion, (iii) an undulatory movement, and (iv) a contraction-expansion movement. The weak form of the governing equations has been formulated in a Cartesian coordinate system and taking into account the primitive variables of the flow field. A fully implicit pressure correction approach, which requires at each time step an iterative process to solve in a sequential manner the equations which govern… More >

Peter Lucas¹, Alexander H. van Zuijlen¹, Hester Bijl¹

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.2, pp. 147-176, 2009, DOI:10.3970/cmes.2009.041.147

Abstract Mesh adaptation is a fairly established tool to obtain numerically accurate solutions for flow problems. Computational efficiency is, however, not always guaranteed for the adaptation strategies found in literature. Typically excessive mesh growth diminishes the potential efficiency gain. This paper, therefore, extends the strategy proposed by [Aftosmis and Berger (2002)] to compute the refinement threshold. The extended strategy computes the refinement threshold based on a user desired number of grid cells and adaptations, thereby ensuring high computational efficiency. Because our main interest is flow around wind turbines, the adaptation strategy has been optimized for flow around wind turbine airfoils. The… More >

Marco La Mantia¹, Peter Dabnichki^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 131-154, 2008, DOI:10.3970/cmes.2008.033.131

Abstract A potential flow based boundary element method was devised to obtain the hydrodynamic forces acting on oscillating wings. A new formulation of the unsteady Kutta condition, postulating a finite pressure difference at the trailing edge of the flapping wing and proposed earlier by the authors, is implemented in the numerical procedure. A comparison with published experimental data (Read et al., 2003) is carried out and the three-dimensional computational results showed good agreement, especially if compared with a similar two-dimensional numerical approach (La Mantia and Dabnichki, 2008) and the potential analytical model of Garrick (1936). The need of considering the differences… More >

Hawa Singh¹, Paras Ram², Vikas Kumar³

FDMP-Fluid Dynamics & Materials Processing, Vol.10, No.4, pp. 521-550, 2014, DOI:10.3970/fdmp.2014.010.521

Abstract The fluctuating flow produced by magneto - hydrodynamic free convection past an impulsively started isothermal vertical plate is studied taking into account the effects of radiation and viscous dissipation. By using the similarity transformation, the governing equations are transformed into dimensionless form and then the system of nonlinear partial differential equations is solved by a perturbation technique. The considered uniform magnetic field acts perpendicular to the plate, which absorbs the fluid with a given suction velocity. A comparison is made in velocity and temperature profiles for two particular cases of real and imaginary time dependent functions. The effects of various… More >

Marcello Lappa^{1, 2}

FDMP-Fluid Dynamics & Materials Processing, Vol.1, No.1, pp. 1-20, 2005, DOI:10.3970/fdmp.2005.001.001

Abstract Possible natural transport mechanisms in cubical and shallow cavities with different heating conditions (from below or from the side) are investigated by means of numerical solution of the non-linear model equations and multiprocessor computations. Attention is focused on a variety of three-dimensional steady effects that can arise in such configurations in the case of low-Pr liquids (silicon melt) even for relatively small values of the temperature gradient due to localized boundary effects and/or true instabilities of the flow. Such aspects are still poorly known or completely ignored owing to the fact that most of the existing experiments focused on the… More >