1. BEM simulation of compressible fluid flow in an enclosure induced by thermoacoustic wavesLeopold Škerget, Jure Ravnik, 2009, original scientific article Abstract: The problem of unsteady compressible fluid flow in an enclosure induced by thermoacoustic waves is studied numerically. Full compressible set of Navier-Stokes equations are considered and numerically solved by boundary-domain integral equations approach coupled with wavelet compression and domain decomposition to achieve numerical efficiency. The thermal energy equation is written in its most general form including the Rayleigh and reversible expansion rate terms. Both, the classical Fourier heat flux model and wave heat conduction model are investigated. The velocity-vorticity formulation of the governing Navier-Stokes equations is employed, while the pressure field is evaluated from the corresponding pressure Poisson equation. Material properties are taken to be for the perfect gas, and assumed to be pressure and temperature dependent. Keywords: compressible fluid flow, boundary element method, thermoacoustic waves, velocity-vorticity fomulation Published: 31.05.2012; Views: 1296; Downloads: 59 Link to full text |
2. Fast single domain-subdomain BEM algorithm for 3D incompressible fluid flow and heat transferJure Ravnik, Leopold Škerget, Zoran Žunič, 2009, original scientific article Abstract: In this paper acceleration and computer memory reduction of an algorithm for the simulation of laminar viscous flows and heat transfer is presented. The algorithm solves the velocity-vorticity formulation of the incompressible Navier-Stokes equations in 3D. It is based on a combination of a subdomain boundary element method (BEM) and single domain BEM. The CPU time and storage requirements of the single domain BEM are reduced by implementing a fast multipole expansion method. The Laplace fundamental solution, which is used as a special weighting function in BEM, is expanded in terms of spherical harmonics. The computational domain and its boundary are recursively cut up forming a tree of clusters of boundary elements and domain cells. Data sparse representation is used in parts of the matrix, which correspond to boundary-domain clusters pairs that are admissible for expansion. Significant reduction of the complexity is achieved. The paper presents results of testing of the multipole expansion algorithm by exploring its effect on the accuracy of the solution and its influence on the non-linear convergence properties of the solver. Two 3D benchmark numerical examples are used: the lid-driven cavity and the onset of natural convection in a differentially heated enclosure. Keywords: boundary element method, fast multipole method, fluid flow, heat transfer, velocity-vorticity fomulation Published: 31.05.2012; Views: 1307; Downloads: 47 Link to full text |
3. SPLETNA APLIKACIJA ZA GENERIRANJE POLJUBNO STRUKTURIRANIH DATOTEKČrtomir Drofenik, 2011, bachelor thesis/paper Abstract: V diplomski nalogi obravnavamo rešitve na področju pretvarjanja oz. prenašanja podatkov iz različnih podatkovnih virov v različne tipe datotek. V osrednjem delu se osredotočimo na izdelavo spletne aplikacije, s katero uporabniku omogočamo določanje strukture podatkov, ki jih bodo shranili v določen tip končne datoteke. Podrobneje opišemo tudi podatkovno strukturo, ki smo jo ustvarili, kakor tudi uporaba posebenega prejavalnika za generiranje datotek. Keywords: spletna aplikacija, generiranje datotek, razpoznavanje in generiranje podatkovnih struktur, Apache Velocity Published: 25.09.2011; Views: 1223; Downloads: 75 Full text (1,62 MB) |
4. Two-dimensional velocity-vorticity based LES for the solution of natural convection in a differentially heated enclosure by wavelet transform based BEM and FEMJure Ravnik, Leopold Škerget, Matjaž Hriberšek, 2006, original scientific article Abstract: 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. Keywords: numerical modelling, boundary element method, discrete wavelet transform, large eddy simulation, velocity-vertocity formulation, natural convection Published: 31.05.2012; Views: 1499; Downloads: 49 Link to full text |
5. Improved method for shot particle velocity measurement within a shotblasting chamberAleš Hribernik, Gorazd Bombek, 2006, original scientific article Abstract: This paper presents the application of two axially-spaced impact detecting sensors for correlation-based shot velocity measurement. The sensors were placed directly into the flow of shot particles for sensing the particles' impact. The implementation of cross-correlation for phase shift determination, however, was not an easy task. Several extensive improvements were necessary, before the original measurement system, already presented in this journal in December 2003, was capable of reliable velocity measurements. This paper presents the development of the system and optimisation of its components. The construction of a sensor built from a small microphone encapsulated in a metal housing and placed under a protective shield was optimised for a high ratio of frontal to side impact signal alteration. High vibration damping was also necessary. It was achieved by the combination of a well-weighted design for supporting elements and by a careful selection of used materials. These measures improved the similarity of the signals generated by the particlesć impact, significantly, and allowed the use of cross-correlation. The accuracy of the obtained results increased drastically by the application of cross-correlation, signal filtering was unnecessary and the signal processing was much faster. Keywords: velocity measurements, shotblasting Published: 30.05.2012; Views: 1231; Downloads: 66 Link to full text |
6. Boundary domain integral method for the study of double diffusive natural convection in porous mediaJanja Kramer Stajnko, Renata Jecl, Leopold Škerget, 2007, original scientific article Abstract: The main purpose of this paper is to present a boundary domain integral method(BDIM) for the solution of natural convection in porous media driven by combining thermal and solutal buoyancy forces. The Brinkman extension of the classical Darcy equation is used for the momentum conservation equation. The numerical scheme was tested on a natural convection problem within a square porous cavity, where different temperature and concentration values are applied on the vertical walls, while the horizontal walls are adiabatic and impermeable. The results for different governing parameters (Rayleigh number, Darcy number, buoyancy ratio and Lewis number) are presented and compared withpublished work. There is a good agreement between those results obtained using the presented numerical scheme and reported studies using other numerical methods. Keywords: double diffusive natural convection, porous medium, velocity-vorticity formulation, Brinkman extended Darcy formulation, boundary domain integrated method Published: 31.05.2012; Views: 1097; Downloads: 54 Link to full text |
7. Velocity vorticity-based large eddy simulation with the bounadr element methodJure Ravnik, Leopold Škerget, Matjaž Hriberšek, 2006, published scientific conference contribution (invited lecture) Abstract: A large eddy simulation using the velocity-vorticity formulation of the incompressible Navier-Stokes equations in combination with the turbulent heat transfer equation is proposed for the solution of the turbulent natural convection drive flow in a 1:4 enclosure. The system of equations is closed by an enthropy-based subgrid scale model.The Prandtl turbulent number is used to estimate turbulent diffusion in the heat transfer equation. The boundary element method is used to solve the kinematics equation and estimate the boundary vorticity values. The vorticity transport equation is solved by FEM. The numerical example studied in this paper is the onset of a turbulent flow regime occuring at high Rayleigh number values ▫$(Ra=10^7-10^10)$▫. The formation of vortices in the boundary layer is observed, along with buoyancy driven diffusive convective transport. Quantitative comparison with the laminar flow model and the worh of other authors is also presented in terms of Nusselt number value oscillations. Keywords: fluid mechanics, incompressible viscous fluid, turbulent flow, velocity vorticity formulation, finite element method, large eddy simulation Published: 31.05.2012; Views: 1227; Downloads: 18 Link to full text |
8. Velocity-vorticity formulation for 3D natural convection in an inclined enclosure by BEMJure Ravnik, Leopold Škerget, Zoran Žunič, 2008, original scientific article Abstract: 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^3Keywords: 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 Published: 31.05.2012; Views: 1519; Downloads: 61 Link to full text |
9. |
10. Boundary domain integral method for high Reynolds viscous fluid flows in complex planar geometriesMatjaž Hriberšek, Leopold Škerget, 2005, original scientific article Abstract: 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. Keywords: fluid mechanics, numerical methods, boundary domain integral method, algorithms, incompressible fluid flow, Navier-Stokes equations, velocity vorticity formulation, segmentation technique, driven cavity flow Published: 01.06.2012; Views: 1118; Downloads: 46 Link to full text |