A.I.Tolstykh, D.A. Shirobokov¹

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.2, pp. 207-222, 2005, DOI:10.3970/cmes.2005.007.207

Abstract A way of using RBF as the basis for PDE's solvers is presented, its essence being constructing approximate formulas for derivatives discretizations based on RBF interpolants with local supports similar to stencils in finite difference methods. Numerical results for different types of elasticity equations showing reasonable accuracy and good$h$-convergence properties of the technique are presented. Applications of the technique to problems with non-self-adjoint operators (like those for the Navier-Stokes equations) are also considered. More >

C. Shu^1,2, H. Ding², K.S. Yeo²

CMES-Computer Modeling in Engineering & Sciences, Vol.7, No.2, pp. 195-206, 2005, DOI:10.3970/cmes.2005.007.195

Abstract Local radial basis function-based differential quadrature (RBF-DQ) method was recently proposed by us. The method is a natural mesh-free approach. It can be regarded as a combination of the conventional differential quadrature (DQ) method with the radial basis functions (RBFs) by means of taking the RBFs as the trial functions in the DQ scheme. With the computed weighting coefficients, the method works in a very similar fashion as conventional finite difference schemes. In this paper, we mainly concentrate on the applications of the method to incompressible flows in the steady and unsteady regions. The multiquadric (MQ) radial basis functions are… More >

Tong-qing Guo, Zhi-liang Lu¹

CMC-Computers, Materials & Continua, Vol.2, No.2, pp. 105-112, 2005, DOI:10.3970/cmc.2005.002.105

Abstract In this paper, a rapid deforming technique is developed to generate dynamic, three-dimensional, multi-block, mesh. The second-order Runge-Kutta time-marching method is used to solve the structural equations of motion. A dual-time method and finite volume discretization are applied for the unsteady Euler/Navier-Stokes equations to calculate the aerodynamic forces, in which the physical time step is synchronous with the structural equations of motion. The Spalart-Allmaras turbulence model is adopted for a turbulent flow. Due to mass dissimilarity, exiting in flutter calculations for a compressible flow, methods of variable mass and variable stiffness are developed to calculate the dynamic pressure of flutter… More >

Zhendong Luo¹, Fei Teng²

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.1, pp. 33-58, 2014, DOI:10.3970/cmes.2014.101.033

Abstract A semi-discrete Crank-Nicolson (CN) formulation about time and a fully discrete stabilized CN finite volume element (SCNFVE) formulation based on two local Gauss integrals and parameter-free with the second-order time accuracy are established for the non-stationary Navier-Stokes equations. The error estimates of the semi-discrete and fully discrete SCNFVE solutions are derived. Some numerical experiments are presented to illustrate that the fully discrete SCNFVE formulation possesses more advantages than its stabilized finite volume element formulation with the first-order time accuracy, thus validating that the fully discrete SCNFVE formulation is feasible and efficient for finding the numerical solutions of the non-stationary Navier-Stokes… More >

Marcin Kamiński¹, Rafał Leszek Ossowski¹

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.3&4, pp. 311-334, 2011, DOI:10.3970/cmes.2011.081.311

Abstract The main aim of this article is numerical solution of the fully coupled Navier-Stokes equations with Gaussian random parameters. It is provided thanks to the specially adopted Finite Volume Method, modified using the generalized stochastic perturbation technique. This Stochastic Finite Volume Method is applied to model 3D problem with uncertainty in liquid viscosity and a coefficient of the heat conduction, separately. Probabilistic moments and characteristics of up to the fourth order are determined with the use of the Response Function Method realized numerically via the polynomial inpterpolation. Although mathematical formulation of the SFVM has been proposed in addition to the… More >

Zhiming Gao², Yichen Ma³

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.2, pp. 135-164, 2010, DOI:10.3970/cmes.2010.066.135

Abstract Shape optimization technique has an increasing role in fluid dynamics problems governed by distributed parameter systems. In this paper, we present the problem of shape optimization of two dimensional viscous flow governed by the time dependent Navier-Stokes equations. The minimization problem of the viscous dissipated energy was established in the fluid domain. We derive the structure of continuous shape gradient of the cost functional by using the differentiability of a saddle point formulation with a function space parametrization technique. Finally a gradient type algorithm with mesh adaptation and mesh movement strategies is successfully and efficiently applied. More >

Roland Martin¹, Carlos Couder-Castaneda¹

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.1, pp. 47-78, 2010, DOI:10.3970/cmes.2010.063.047

Abstract We develop an unsplit convolutional perfectly matched layer (CPML) technique to absorb efficiently compressible viscous flows and their related supersonic or subsonic regimes at the outer boundary of a distorted computational domain. More particularly subsonic outgoing flows or subsonic wall-boundary layers close to the PML are well absorbed, which is difficult to obtain without creating numerical instabilities over long time periods. This new PML (CPML) introduces the calculation of auxiliary memory variables at each time step and allows an unsplit formulation of the PML. Damping functions involving a high shift in the frequency domain allow a much better absorption of… More >

G.C. Bourantas¹, E.D. Skouras^2,3, V.C. Loukopoulos⁴, G.C. Nikiforidis¹

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.1, pp. 31-64, 2010, DOI:10.3970/cmes.2010.059.031

Abstract A meshfree point collocation method has been developed for the velocity-vorticity formulation of two-dimensional, steady state incompressible Navier-Stokes equations. Particular emphasis was placed on the application of the velocity-correc -tion method, ensuring the continuity equation. The Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are obtained for regular and irregular nodal distributions, stressing the positivity conditions that make the matrix of the system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through two representative, well-known, and… More >

V.C. Mariani¹, E.E.M. Alonso², S. Peters³

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 15-28, 2008, DOI:10.3970/cmes.2008.029.015

Abstract An application of Newton's method for linearization of advective terms given by the discretization on unstructured Voronoi meshes for the incompressible Navier-Stokes equations is proposed and evaluated in this article. One of the major advantages of the unstructured approach is its application to very complex geometrical domains and the mesh is adaptable to features of the flow. Moreover, in this work comparisons with the literature results in bi-dimensional lid-driven cavities for different Reynolds numbers allow us to assess the numerical properties of the new proposed finite-volume scheme. Results for the components of the velocity, and the pressure collocated at the… More >

E.J. Sellountos¹, A. Sequeira¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.2, pp. 127-148, 2008, DOI:10.3970/cmes.2008.023.127

Abstract In this work a hybrid velocity-vorticity scheme for the solution of the 2D Navier-Stokes equations is presented. The multi-region Local Boundary Integral Equation (LBIE) combined with Radial Basis Functions (RBF) interpolation is used for the solution of the kinematics and the multi-region BEM for the solution of the transport kinetics. The final system of equations is in band form for both methods. The issue of RBF discontinuities is resolved by constructing the RBF matrix locally in every region. The kinematics integral equation is used in three different forms, for coupling the velocity field on the boundary, on interior points and… More >