• 한국자동차공학회논문집
• /
• 제11권3호
• /
• pp.115-123
• /
• 2003

With the advent of high performance computers and more efficient numerical algorithms, computational fluid dynamics(CFD) has come out as a modem alternative for reducing the use of wind tunnels test in automotive engineering. However, in spite of the fact that many competent researchers have made all their talents in developing turbulence model over since the past dozen or more years, it has been an important impediment in using the CFD effectively to design machinery and to diagnose or to improve engineering problems in the industry since the turbulence model has been acting as the Achilles' tendon in aspect of the reliability even to this time. In this study, Reynolds-averaged Wavier-Stokes equations were solved to simulate an incompressible turbulent flow around a bus-like bluff body near ground plane. In order to investigate the effect of the discretisation schemes and turbulence model on the aerodynamic forces several turbulence models with five convective difference schemes are adopted. From the results of this study, it is clear that choice of turbulence model and discretisation scheme profoundly affects the computational outcome. The results also show that the adoption of RNG $k-\varepsilon$ turbulence model and nonlinear quadratic turbulence model with the second order accurate discretisation scheme predicts fairly well the aerodynamic coefficients.

• Advances in aircraft and spacecraft science
• /
• 제7권5호
• /
• pp.387-404
• /
• 2020

The numerical simulation of shock waves in supersonic flows is challenging because of several instabilities which can affect the solution. Among them, the carbuncle phenomenon can introduce nonphysical perturbations in captured shock waves. In the present work, a hybrid numerical flux is proposed for the evaluation of the convective fluxes that avoids carbuncle and keeps high-accuracy on shocks and boundary layers. In particular, the proposed flux is a combination between an upwind approximate Riemann problem solver and the Local Lax-Friedrichs scheme. A simple strategy to mix the two fluxes is proposed and tested in the framework of a discontinuous Galerkin discretisation. The approach is investigated on the subsonic flow in a channel, on the supersonic flow around a cylinder, on the supersonic flow on a flat plate and on the flow in a overexpanded rocket nozzle.

• Korea-Australia Rheology Journal
• /
• 제21권1호
• /
• pp.1-12
• /
• 2009

This paper reports the suitability of a domain decomposition technique for the hybrid simulation of dilute polymer solution flows using Eulerian Brownian dynamics and Radial Basis Function Networks (RBFN) based methods. The Brownian Configuration Fields (BCF) and RBFN method incorporates the features of the BCF scheme (which render both closed form constitutive equations and a particle tracking process unnecessary) and a mesh-less method (which eliminates element-based discretisation of domains). However, when dealing with large scale problems, there appear several difficulties: the high computational time associated with the Stochastic Simulation Technique (SST), and the ill-condition of the system matrix associated with the RBFN. One way to overcome these disadvantages is to use parallel domain decomposition (DD) techniques. This approach makes the BCF-RBFN method more suitable for large scale problems.

• 한국전산유체공학회지
• /
• 제3권2호
• /
• pp.17-26
• /
• 1998

The three-dimensional incompressible Navier-Stokes equations have been solved by a node-centered finite volume method with unstructured hybrid grids. The pressure-velocity coupling is handled by the artificial compressibility algorithm and convective fluxes are obtained by Roe's flux difference splitting scheme with linear reconstruction of the solutions. Euler implicit method with Jacobi matrix solver is used for the time-integration. The viscous terms are discretised in a manner to handle any kind of grids such as tetragedra, prisms, pyramids, hexahedra, or mixed-element grid. Inviscid bump flow is solved to check the accuracy of high order convective flux discretisation. And viscous flows around a circular cylinder and a sphere are studied to show the efficiency and accuracy of the solver.

• Structural Engineering and Mechanics
• /
• 제72권1호
• /
• pp.71-81
• /
• 2019

The objective of this paper is to develop a size-dependent nonlinear model of beams for fluid-conveying nanotubes with an initial deflection. The nonlinear frequency response of the nanotube is analysed via an Euler-Bernoulli model. Size influences on the behaviour of the nanosystem are described utilising the nonlocal strain gradient theory (NSGT). Relative motions at the inner wall of the nanotube is taken into consideration via Beskok-Karniadakis model. Formulating kinetic and elastic energies and then employing Hamilton's approach, the nonlinear motion equations are derived. Furthermore, Galerkin's approach is employed for discretisation, and then a continuation scheme is developed for obtaining numerical results. It is observed that an initial deflection significantly alters the frequency response of NSGT nanotubes conveying fluid. For small initial deflections, a hardening nonlinearity is found whereas a softening-hardening nonlinearity is observed for large initial deflections.

• Steel and Composite Structures
• /
• 제49권5호
• /
• pp.587-602
• /
• 2023

This article presents the numerical modelling of transient heat transfer in highly heterogeneous composite materials where the thermal conductivity, specific heat and density are assumed to be directional-dependent. This article uses a coupled finite element-finite difference scheme to perform the transient heat transfer analysis of unidirectional (1D) and multidirectional (2D/3D) functionally graded composite panels. Here, 1D/2D/3D functionally graded structures are subjected to nonuniform heat source and inhomogeneous boundary conditions. Here, the multidirectional functionally graded materials are modelled by varying material properties in individual or in-combination of spatial directions. Here, fully spatial-dependent material properties are evaluated using Voigt's micromechanics scheme via multivariable power-law functions. The weak form is obtained through the Galerkin method and solved further via the element-space and time-step discretisation through the 2D-isoparametric finite element and the implicit backward finite difference schemes, respectively. The present model is verified by comparing it with the previously reported results and the commercially available finite element tool. The numerous illustrations confirm the significance of boundary conditions and material heterogeneity on the transient temperature responses of 1D/2D/3D functionally graded panels.

• Steel and Composite Structures
• /
• 제9권2호
• /
• pp.131-158
• /
• 2009

This paper presents a modelling technique for the nonlinear analysis of composite steel-concrete beams with partial shear interaction. It extends the applicability of two stiffness elements previously derived by the authors using the direct stiffness method, i.e. the 6DOF and the 8DOF elements, to account for material nonlinearities. The freedoms are the vertical displacement, the rotation and the slip at both ends for the 6DOF stiffness element, as well as the axial displacement at the level of the reference axis for the 8DOF stiffness element. The solution iterative scheme is based on the secant method, with the convergence criteria relying on the ratios of the Euclidean norms of both forces and displacements. The advantage of the approach is that the displacement and force fields of the stiffness elements are extremely rich as they correspond to those required by the analytical solution of the elastic partial interaction problem, thereby producing a robust numerical technique. Experimental results available in the literature are used to validate the finite element proposed in the paper. For this purpose, those reported by Chapman and Balakrishnan (1964), Fabbrocino et al. (1998, 1999) and Ansourian (1981) are utilised; these consist of six simply supported beams with a point load applied at mid-span inducing positive bending moment in the beams, three simply supported beams with a point load applied at mid-span inducing negative bending moment in the beams, and six two-span continuous composite beams respectively. Based on these comparisons, a preferred degree of discretisation suitable for the proposed modelling technique expressed as a function of the ratio between the element length and depth is proposed, as is the number of Gauss stations needed. This allows for accurate prediction of the nonlinear response of composite beams.