Abstract: In Computational Fluid Dynamics (CFD), there are a variety of numerical methods, of which some depend on macroscopic model representatives. These models can be solved by finite-volume, finite-element or finite-difference methods on a microscopic description. However, the lattice Boltzmann method (LBM) is considered to be a mesoscopic particle method, with its scale lying between the macroscopic and microscopic scales. The LBM works well for solving incompressible flow problems, but certain limitations arise from solving compressible flows, particularly at high Mach numbers. An improved lattice Boltzmann model for compressible flow problems is presented in this research study. A higher-order Taylor series expansion of the Maxwell equilibrium distribution function is used to overcome limitations in LBM when solving high-Mach-number flows. Large eddy simulation (LES) is implemented in LBM to simulate turbulent jet flows. The results have been validated with available experimental data for turbulent compressible free jet flow at subsonic speeds.
Abstract: In this study, we have analyzed the transport of analytes
under a two dimensional steady incompressible flow of power-law
fluids through rectangular nanochannel. A mathematical model
based on the Cauchy momentum-Nernst-Planck-Poisson equations is
considered to study the combined effect of mixed electroosmotic
(EO) and pressure driven (PD) flow. The coupled governing
equations are solved numerically by finite volume method. We
have studied extensively the effect of key parameters, e.g., flow
behavior index, concentration of the electrolyte, surface potential,
imposed pressure gradient and imposed electric field strength on
the net average flow across the channel. In addition to study
the effect of mixed EOF and PD on the analyte distribution
across the channel, we consider a nonlinear model based on
general convective-diffusion-electromigration equation. We have also
presented the retention factor for various values of electrolyte
concentration and flow behavior index.
Abstract: In designing a low-energy-consuming buildings, the heat transfer through a large glass or wall becomes critical. Multiple layers of the window glasses and walls are employed for the high insulation. The gravity driven air flow between window glasses or wall layers is a natural heat convection phenomenon being a key of the heat transfer. For the first step of the natural heat transfer analysis, in this study the development and application of a finite volume method for the numerical computation of viscous incompressible flows is presented. It will become a part of the natural convection analysis with high-order scheme, multi-grid method, and dual-time step in the future. A finite volume method based on a fully-implicit second-order is used to discretize and solve the fluid flow on unstructured grids composed of arbitrary-shaped cells. The integrations of the governing equation are discretised in the finite volume manner using a collocated arrangement of variables. The convergence of the SIMPLE segregated algorithm for the solution of the coupled nonlinear algebraic equations is accelerated by using a sparse matrix solver such as BiCGSTAB. The method used in the present study is verified by applying it to some flows for which either the numerical solution is known or the solution can be obtained using another numerical technique available in the other researches. The accuracy of the method is assessed through the grid refinement.
Abstract: Thermal enhancement of a single mini channel in
Proton Exchange Membrane Fuel Cell (PEMFC) cooling plate is
numerically investigated. In this study, low concentration of Al2O3 in
Water - Ethylene Glycol mixtures is used as coolant in single channel
of carbon graphite plate to mimic the mini channels in PEMFC
cooling plate. A steady and incompressible flow with constant heat
flux is assumed in the channel of 1mm x 5mm x 100mm. Nano
particle of Al2O3 used ranges from 0.1, 0.3 and 0.5 vol %
concentration and then dispersed in 60:40 (water: Ethylene Glycol)
mixture. The effect of different flow rates to fluid flow and heat
transfer enhancement in Re number range of 20 to 140 was observed.
The result showed that heat transfer coefficient was improved by
18.11%, 9.86% and 5.37% for 0.5, 0.3 and 0.1 vol. % Al2O3 in 60:40
(water: EG) as compared to base fluid of 60:40 (water: EG). It is also
showed that the higher vol. % concentration of Al2O3 performed
better in term of thermal enhancement but at the expense of higher
pumping power required due to increase in pressure drop
experienced. Maximum additional pumping power of 0.0012W was
required for 0.5 vol % Al2O3 in 60:40 (water: EG) at Re number 140.
Abstract: This paper presents a fully Lagrangian coupled
Fluid-Structure Interaction (FSI) solver for simulations of
fluid-structure interactions, which is based on the Moving Particle
Semi-implicit (MPS) method to solve the governing equations
corresponding to incompressible flows as well as elastic structures.
The developed solver is verified by reproducing the high velocity
impact loads of deformable thin wedges with three different materials
such as mild steel, aluminium and tin during water entry. The present
simulation results for aluminium are compared with analytical solution
derived from the hydrodynamic Wagner model and linear Wan’s
theory. And also, the impact pressure and strain on the water entry
wedge with three different materials, such as mild steel, aluminium
and tin, are simulated and the effects of hydro-elasticity are discussed.
Abstract: This paper is devoted to the study of a viscous
incompressible flow around a circular cylinder performing harmonic
oscillations, especially the steady streaming phenomenon. The
research methodology is based on the asymptotic explanation method
combined with the computational bifurcation analysis. The research
approach develops Schlichting and Wang decomposition method.
Present studies allow to identify several regimes of the secondary
streaming with different flow structures. The results of the research
are in good agreement with experimental and numerical simulation
data.
Abstract: The practical application of the Computational Fluid Dynamics (CFD), for predicting the flow pattern around Multipurpose Amphibious Vehicle (MAV) hull has made much progress over the last decade. Today, several of the CFD tools play an important role in the land and water going vehicle hull form design. CFD has been used for analysis of MAV hull resistance, sea-keeping, maneuvering and investigating its variation when changing the hull form due to varying its parameters, which represents a very important task in the principal and final design stages. Resistance analysis based on CFD (Computational Fluid Dynamics) simulation has become a decisive factor in the development of new, economically efficient and environmentally friendly hull forms. Three-dimensional finite volume method (FVM) based on Reynolds Averaged Navier-Stokes equations (RANS) has been used to simulate incompressible flow around three types of MAV hull bow models in steady-state condition. Finally, the flow structure and streamlines, friction and pressure resistance and velocity contours of each type of hull bow will be compared and discussed.
Abstract: A method for simulating flow around the solid bodies has been presented using hybrid meshfree and mesh-based schemes. The presented scheme optimizes the computational efficiency by combining the advantages of both meshfree and mesh-based methods. In this approach, a cloud of meshfree nodes has been used in the domain around the solid body. These meshfree nodes have the ability to efficiently adapt to complex geometrical shapes. In the rest of the domain, conventional Cartesian grid has been used beyond the meshfree cloud. Complex geometrical shapes can therefore be dealt efficiently by using meshfree nodal cloud and computational efficiency is maintained through the use of conventional mesh-based scheme on Cartesian grid in the larger part of the domain. Spatial discretization of meshfree nodes has been achieved through local radial basis functions in finite difference mode (RBF-FD). Conventional finite difference scheme has been used in the Cartesian ‘meshed’ domain. Accuracy tests of the hybrid scheme have been conducted to establish the order of accuracy. Numerical tests have been performed by simulating two dimensional steady and unsteady incompressible flows around cylindrical object. Steady flow cases have been run at Reynolds numbers of 10, 20 and 40 and unsteady flow problems have been studied at Reynolds numbers of 100 and 200. Flow Parameters including lift, drag, vortex shedding, and vorticity contours are calculated. Numerical results have been found to be in good agreement with computational and experimental results available in the literature.
Abstract: In this paper a numerical algorithm is described for solving the boundary value problem associated with axisymmetric, inviscid, incompressible, rotational (and irrotational) flow in order to obtain duct wall shapes from prescribed wall velocity distributions. The governing equations are formulated in terms of the stream function ψ (x,y)and the function φ (x,y)as independent variables where for irrotational flow φ (x,y)can be recognized as the velocity potential function, for rotational flow φ (x,y)ceases being the velocity potential function but does remain orthogonal to the stream lines. A numerical method based on the finite difference scheme on a uniform mesh is employed. The technique described is capable of tackling the so-called inverse problem where the velocity wall distributions are prescribed from which the duct wall shape is calculated, as well as the direct problem where the velocity distribution on the duct walls are calculated from prescribed duct geometries. The two different cases as outlined in this paper are in fact boundary value problems with Neumann and Dirichlet boundary conditions respectively. Even though both approaches are discussed, only numerical results for the case of the Dirichlet boundary conditions are given. A downstream condition is prescribed such that cylindrical flow, that is flow which is independent of the axial coordinate, exists.
Abstract: In this article, various models of surface tension force (CSF, CSS and PCIL) for interfacial flows have been applied to dynamic case and the results were compared. We studied the Kelvin- Helmholtz instabilities, which are produced by shear at the interface between two fluids with different physical properties. The velocity inlet is defined as a sinusoidal perturbation. When gravity and surface tension are taking into account, we observe the development of the Instability for a critic value of the difference of velocity of the both fluids. The VOF Model enables to simulate Kelvin-Helmholtz Instability as dynamic case.
Abstract: Along with increasing development of generation of supersonic planes especially fighters and request for increasing the performance and maneuverability scientists and engineers suggested the delta and double delta wing design. One of the areas which was necessary to be researched, was the Aerodynamic review of this type of wings in high angles of attack at low speeds that was very important in landing and takeoff the planes and maneuvers. Leading Edges of the wings,cause the separation flow from wing surface and then formation of powerful vortex with high rotational speed which studing the mechanism and location of formation and also the position of the vortex breakdown in high angles of attack is very important. In this research, a double delta wing with 76o/45o sweep angles at high angle of attack in steady state and incompressible flow were numerically analyzed with Fluent software. With analaysis of the numerical results, we arrived the most important characteristic of the double delta wings which is keeping of lift at high angles of attacks.
Abstract: The present work is concerned with the free
convective two dimensional flow and heat transfer, in isotropic fluid
filled porous rectangular enclosure with differentially heated walls for
steady state incompressible flow have been investigated for non-
Darcy flow model. Effects of Darcy number (0.0001 £Da£ 10),
Rayleigh number (10 £Ra£ 5000), and aspect ratio (0.25 £AR£ 4), for
a range of porosity (0.4 £e£ 0.9) with and without moving lower wall
have been studied. The cavity was insulated at the lower and upper
surfaces. The right and left heated surfaces allows convective
transport through the porous medium, generating a thermal
stratification and flow circulations. It was found that the Darcy
number, Rayleigh number, aspect ratio, and porosity considerably
influenced characteristics of flow and heat transfer mechanisms. The
results obtained are discussed in terms of the Nusselt number,
vectors, contours, and isotherms.
Abstract: A new conserving approach in the context of Immersed Boundary Method (IBM) is presented to simulate one dimensional, incompressible flow in a moving boundary problem. The method employs control volume scheme to simulate the flow field. The concept of ghost node is used at the boundaries to conserve the mass and momentum equations. The Present method implements the conservation laws in all cells including boundary control volumes. Application of the method is studied in a test case with moving boundary. Comparison between the results of this new method and a sharp interface (Image Point Method) IBM algorithm shows a well distinguished improvement in both pressure and velocity fields of the present method. Fluctuations in pressure field are fully resolved in this proposed method. This approach expands the IBM capability to simulate flow field for variety of problems by implementing conservation laws in a fully Cartesian grid compared to other conserving methods.
Abstract: The three-dimensional incompressible flow past a
rectangular open cavity is investigated, where the aspect ratio of the
cavity is considered as 4. The principle objective is to use large-eddy
simulation to resolve and control the large-scale structures, which are
largely responsible for flow oscillations in a cavity. The flow past an
open cavity is very common in aerospace applications and can be a
cause of acoustic source due to hydrodynamic instability of the shear
layer and its interactions with the downstream edge. The unsteady
Navier-stokes equations have been solved on a staggered mesh using
a symmetry-preserving central difference scheme. Synthetic jet has
been used as an active control to suppress the cavity oscillations in
wake mode for a Reynolds number of ReD = 3360. The effect of
synthetic jet has been studied by varying the jet amplitude and
frequency, which is placed at the upstream wall of the cavity. The
study indicates that there exits a frequency band, which is larger than
a critical value, is effective in attenuating cavity oscillations when
blowing ratio is more than 1.0.
Abstract: The steady incompressible flow has been solved in cylindrical coordinates in both vapour region and wick structure. The governing equations in vapour region are continuity, Navier-Stokes and energy equations. These equations have been solved using SIMPLE algorithm. For study of parameters variation on heat pipe operation, a benchmark has been chosen and the effect of changing one parameter has been analyzed when the others have been fixed.
Abstract: In this paper, the local grid refinement is focused by
using a nested grid technique. The Cartesian grid numerical method is
developed for simulating unsteady, viscous, incompressible flows
with complex immersed boundaries. A finite volume method is used in
conjunction with a two-step fractional-step procedure. The key aspects
that need to be considered in developing such a nested grid solver are
imposition of interface conditions on the inter-block and accurate
discretization of the governing equation in cells that are with the
inter-block as a control surface. A new interpolation procedure is
presented which allows systematic development of a spatial
discretization scheme that preserves the spatial accuracy of the
underlying solver. The present nested grid method has been tested by
two numerical examples to examine its performance in the two
dimensional problems. The numerical examples include flow past a
circular cylinder symmetrically installed in a Channel and flow past
two circular cylinders with different diameters. From the numerical
experiments, the ability of the solver to simulate flows with
complicated immersed boundaries is demonstrated and the nested grid
approach can efficiently speed up the numerical solutions.
Abstract: In this paper, fluid flow patterns of steady incompressible flow inside shear driven cavity are studied. The numerical simulations are conducted by using lattice Boltzmann method (LBM) for different Reynolds numbers. In order to simulate the flow, derivation of macroscopic hydrodynamics equations from the continuous Boltzmann equation need to be performed. Then, the numerical results of shear-driven flow inside square and triangular cavity are compared with results found in literature review. Present study found that flow patterns are affected by the geometry of the cavity and the Reynolds numbers used.