Iskanje po katalogu digitalne knjižnice

Opis: This paper presents a mixed boundary element formulation of the boundary domain integral method (BDIM) for solving diffusion-convective transport problems. The basic idea of mixed elements is the use of a continuos interpolation polynomial for conservative field function approximation and a discontinuous interpolation polynomial for its normal derivative along the boundary element. In this way, the advantages of continuous field function approximation are retained and its conservation is preserved while the normal flux values are approximated by interpolation nodal points with a uniquely defined normal direction. Due to the use of mixed boundary elements, the final discretized matrix system is overdetermined and a special solver based on the least squares method is applied. Driven cavity, natural and forced convection in a closed cavity are studied. Driven caviaty results at Re=100, 400 and 1000 agree better with the benchmark solution than Finite Element Method of Finite Volume Method results for the same grid density with 21 x 21 degrees of freedom. The average Nusselt number values for natural convection ▫$10^3$▫▫$le$▫Ra▫$le$▫▫$10^6$▫ agree better than 0.1% with benchmark solutions for maximal calculated grid desities 61 x 61 degrees for freedom.
Ključne besede: fluid mechanics, incompressible fluid, laminar flow, velocity vorticity formulation, boundary element method, mixed boundary elements
Objavljeno v DKUM: 01.06.2012; Ogledov: 1801; Prenosov: 91
Povezava na celotno besedilo

Opis: A wavelet matrix compression technique was used to solve systems of linear equations resulting from BEM applied to fluid dynamics. The governing equations were written in velocity-vorticity formulation and solutions of the resulting systems of equations were obtained with and without wavelet matrix compression. A modification of the Haar wavelet transform, which can transformvectors of any size, is proposed. The threshold, used for making fully populated matrices sparse, was written as a product of a user defined factor and the average value of absolute matrix elements values. Numerical tests were performed to assert, that the error caused by wavelet compression depends linearly on the factor , while the dependence of the error on the share of thresholded elements in the system matrix is highly non-linear. The results also showed that the increasing non-linearity (higher Ra and Re numbervalues) limits the extent of compression. On the other hand, higher meshdensity enables higher compression ratios.
Ključne besede: fluid mechanics, computational fluid dynamics, boundary element method, wavelet transform, linear systems of equations, velocity vorticity formulation, driven cavity, natural convection, system matrix compression
Objavljeno v DKUM: 01.06.2012; Ogledov: 1806; Prenosov: 91
Povezava na celotno besedilo

Opis: The article presents new developments in boundary domain integral method (BDIM) for computation of viscous fluid flows, governed by the Navier-Stokes equations. The BDIM algorithm uses velocity-vorticity formulation and is basedon Poisson velocity equation for flow kinematics. This results in accurate determination of boundary vorticity values, a crucial step in constructing an accurate numerical algorithm for computation of flows in complex geometries, i.e. geometries with sharp corners. The domain velocity computations are done by the segmentation technique using large segments. After solving the kinematics equation the vorticity transport equation is solved using macro-element approach. This enables the use of macro-element based diffusion-convection fundamental solution, a key factor in assuring accuracy of computations for high Reynolds value laminar flows. The versatility and accuracy of the proposed numerical algorithm is shown for several test problems, including the standard driven cavity together with the driven cavity flow in an L shaped cavity and flow in a Z shaped channel. The values of Reynolds number reach 10,000 for driven cavity and 7500 for L shapeddriven cavity, whereas the Z shaped channel flow is computed up to Re = 400. The comparison of computational results shows that the developed algorithm is capable of accurate resolution of flow fields in complex geometries.
Ključne besede: fluid mechanics, numerical methods, boundary domain integral method, algorithms, incompressible fluid flow, Navier-Stokes equations, velocity vorticity formulation, segmentation technique, driven cavity flow
Objavljeno v DKUM: 01.06.2012; Ogledov: 1757; Prenosov: 83
Povezava na celotno besedilo

Opis: A natural convection phenomenon is studied in cubic and parallelepipedal inclined enclosures. The simulation of coupled laminar viscous flow and heat transfer is performed using a novel algorithm based on a combination of singledomain Boundary element method (BEM) and subdomain BEM. The algorithm solves the velocity-vorticity formulation of the incompressible Navier-Stokes equations coupled with the energy equation using the Boussinesq approximation.The subdomain BEM is used to solve the kinematics equation, the vorticity transport equation and the energy equation. The boundary vorticity values, which are needed as boundary conditions for the vorticity transport equation, are calculated by singe domain BEM solution of the kinematics equation. Simulation results are compared with benchmark results for a cubic inclined enclosure for Rayleigh number values ▫$10^3Ključne besede: podobmočna metoda robnih elementov, hitrostno-vrtinčna formulacija, laminarni tok viskozne tekočine, naravna konvekcija, nagnjena kotanja, fluid mechanics, subdomain boundary element method, velocity-vorticity formulation, laminar viscous fluid flow, natural convection, inclined enclosure
Objavljeno v DKUM: 31.05.2012; Ogledov: 2229; Prenosov: 95
Povezava na celotno besedilo

Opis: A wavelet transform based boundary element method (BEM) numerical scheme is proposed for the solution of the kinematics equation of the velocity-vorticityformulation of Navier-Stokes equations. FEM is used to solve the kinetics equations. The proposed numerical approach is used to perform two-dimensional vorticity transfer based large eddy simulation on grids with 105 nodes. Turbulent natural convection in a differentially heated enclosure of aspect ratio 4 for Rayleigh number values Ra=107-109 is simulated. Unstable boundary layer leads to the formation of eddies in the downstream parts of both vertical walls. At the lowest Rayleigh number value an oscillatory flow regime is observed, while the flow becomes increasingly irregular, non-repeating, unsymmetric and chaotic at higher Rayleigh number values. The transition to turbulence is studied with time series plots, temperature-vorticity phase diagrams and with power spectra. The enclosure is found to be only partially turbulent, what is qualitatively shown with second order statistics-Reynolds stresses, turbulent kinetic energy, turbulent heat fluxes and temperature variance. Heat transfer is studied via the average Nusselt number value, its time series and its relationship to the Rayleigh number value.
Ključne besede: numerical modelling, boundary element method, discrete wavelet transform, large eddy simulation, velocity-vertocity formulation, natural convection
Objavljeno v DKUM: 31.05.2012; Ogledov: 2162; Prenosov: 79
Povezava na celotno besedilo