76 research outputs found

    The Investigation of Efficiency of Physical Phenomena Modelling Using Differential Equations on Distributed Systems

    Get PDF
    This work is dedicated to development of mathematical modelling software. In this dissertation numerical methods and algorithms are investigated in software making context. While applying a numerical method it is important to take into account the limited computer resources, the architecture of these resources and how do methods affect software robustness. Three main aspects of this investigation are that software implementation must be efficient, robust and be able to utilize specific hardware resources. The hardware specificity in this work is related to distributed computations of different types: single CPU with multiple cores, multiple CPUs with multiple cores and highly parallel multithreaded GPU device. The investigation is done in three directions: GPU usage for 3D FDTD calculations, FVM method usage to implement efficient calculations of a very specific heat transferring problem, and development of special techniques for software for specific bacteria self organization problem when the results are sensitive to numerical methods, initial data and even computer round-off errors. All these directions are dedicated to create correct technological components that make a software implementation robust and efficient. The time prediction model for 3D FDTD calculations is proposed, which lets to evaluate the efficiency of different GPUs. A reasonable speedup with GPU comparing to CPU is obtained. For FVM implementation the OpenFOAM open source software is selected as a basis for implementation of calculations and a few algorithms and their modifications to solve efficiency issues are proposed. The FVM parallel solver is implemented and analyzed, it is adapted to heterogeneous cluster Vilkas. To create robust software for simulation of bacteria self organization mathematically robust methods are applied and results are analyzed, the algorithm is modified for parallel computations

    The numerical solution of banded linear systems by generallized factorization procedures

    Get PDF
    The numerical solution of banded linear systems by generallized factorization procedure

    Lattice Boltzmann simulation of flow and heat transfer in random porous media constructed by simulated annealing algorithm

    Get PDF
    In this article, the lattice Boltzmann (LB) method for transport phenomena is combined with the simulated annealing (SA) algorithm for digitized porous-medium construction to study flow and heat transfer in random porous media. Importantly, in contrast to previous studies which simplify porous media as arrays of regularly shaped objects or effective pore networks, the LB + SA method in this article can model statistically meaningful random porous structures in irregular morphology, and simulate pore-scale transport processes inside them. Pore-scale isothermal flow and heat conduction in a set of constructed random porous media characterized by statistical descriptors were then simulated through use of the LB + SA method. The corresponding averages over the computational volumes and the related effective transport properties were also computed based on these pore scale numerical results. Good agreement between the numerical results and theoretical predictions or experimental data on the representative elementary volume scale was found. The numerical simulations in this article demonstrate combination of the LB method with the SA algorithm is a viable and powerful numerical strategy for simulating transport phenomena in random porous media in complex geometries

    Seventh Copper Mountain Conference on Multigrid Methods

    Get PDF
    The Seventh Copper Mountain Conference on Multigrid Methods was held on April 2-7, 1995 at Copper Mountain, Colorado. This book is a collection of many of the papers presented at the conference and so represents the conference proceedings. NASA Langley graciously provided printing of this document so that all of the papers could be presented in a single forum. Each paper was reviewed by a member of the conference organizing committee under the coordination of the editors. The vibrancy and diversity in this field are amply expressed in these important papers, and the collection clearly shows the continuing rapid growth of the use of multigrid acceleration techniques

    Numerical Boundary Condition Procedures

    Get PDF
    Topics include numerical procedures for treating inflow and outflow boundaries, steady and unsteady discontinuous surfaces, far field boundaries, and multiblock grids. In addition, the effects of numerical boundary approximations on stability, accuracy, and convergence rate of the numerical solution are discussed

    Simulation of the Steady-State Transport of Radon from Oil intoHouses with Basements under Constant Negative Pressure

    Full text link

    Adaptive mesh refinement for computational aeroacoustics

    Get PDF
    UNIVERSITY OF SOUTHAMPTON ABSTRACT FACULTY OF ENGINEERING, SCIENCE & MATHEMATICS SCHOOL OF ENGINEERING SCIENCES Doctor of Philosophy ADAPTIVE MESH REFINEMENT FOR COMPUTATIONAL AEROACOUSTICS by Xun HuangThis thesis describes a parallel block-structured adaptive mesh refinement (AMR) method that is employed to solve some computational aeroacoustic problems with the aim of improving the computational efficiency. AMR adaptively refines and coarsens a computational mesh along with sound propagation to increase grid resolution only in the area of interest. While sharing many of the same features, there is a marked difference between the current and the established AMR approaches. Rather than low-order schemes generally used in the previous approaches, a high-order spatial difference scheme is employed to improve numerical dispersion and dissipation qualities. To use a high-order scheme with AMR, a number of numerical issues associated with fine-coarse block interfaces on an adaptively refined mesh, such as interpolations, filter and artificial selective damping techniques and accuracy are addressed. In addition, the asymptotic stability and the transient behaviour of a high-order spatial scheme on an adaptively refined mesh are also studied with eigenvalue analysis and pseudospectra analysis respectively. In addition, the fundamental AMR algorithm is simplified in order to make the work of implementation more manageable. Particular emphasis has been placed on solving sound radiation from generic aero-engine bypass geometry with mean flow. The approach of AMR is extended to support a body-fitted multi-block mesh. The radiation from an intake duct is modelled by the linearised Euler equations, while the radiation from an exhaust duct is modelled by the extended acoustic perturbation equations to suppress hydrodynamic instabilities generated in a sheared mean flow. After solving the near-field sound solution, the associated far-field sound directivity is estimated by solving the Ffowcs Williams-Hawkings equation. The overall results demonstrate the accuracy and the efficiency of the presented AMR method, but also reveal some limitations. The possible methods to avoid these limitations are given at the end of this thesis

    Properties of approximate inverses and adaptive control concepts for preconditioning [online]

    Get PDF

    Discrétisation Espace-Temps d'Équations d'Ondes Élasto-Acoustiques dans des Bases Trefftz-DG Polynomiales

    Get PDF
    Discontinuous Finite Element Methods (DG FEM) have proven flexibility and accuracy for solving wave problems in complex media. However, they require a large number of degrees of freedom, which increases the corresponding computational cost compared with that of continuous finite element methods. Among the different variational approaches to solve boundary value problems, there exists a particular family of methods, based on the use of trial functions in the form of exact local solutions of the governing equations. The idea was first proposed by Trefftz in 1926, and since then it has been further developed and generalized. A Trefftz-DG variational formulation applied to wave problems reduces to surface integrals that should contribute to decreasing the computational costs.Trefftz-type approaches have been widely used for time-harmonic problems, while their implementation for time-dependent simulations is still limited. The feature of Trefftz-DG methods applied to time-dependent problems is in the use of space-time meshes. Indeed, standard DG methods lead to the construction of a semi-discrete system of ordinary differential equations in time which are integrated by using an appropriate scheme. But Trefftz-DG methods applied to wave problems lead to a global matrix including time and space discretizations which is huge and sparse. This significantly hampers the deployment of this technology for solving industrial problems.In this work, we develop a Trefftz-DG framework for solving mechanical wave problems including elasto-acoustic equations. We prove that the corresponding formulations are well-posed and we address the issue of solving the global matrix by constructing an approximate inverse obtained from the decomposition of the global matrix into a block-diagonal one. The inversion is then justified under a CFL-type condition. This idea allows for reducing the computational costs but its accuracy is limited to small computational domains. According to the limitations of the method, we have investigated the potential of Tent Pitcher algorithms following the recent works of Gopalakrishnan et al. It consists in constructing a space-time mesh made of patches that can be solved independently under a causality constraint. We have obtained very promising numerical results illustrating the potential of Tent Pitcher in particular when coupled with a Trefftz-DG method involving only surface terms. In this way, the space-time mesh is composed of elements which are 3D objects at most. It is also worth noting that this framework naturally allows for local time-stepping which is a plus to increase the accuracy while decreasing the computational burden.Les méthodes d'éléments finis de type Galerkine discontinu (DG FEM) ont démontré précision et efficacité pour résoudre des problèmes d'ondes dans des milieux complexes. Cependant, elles nécessitent un très grand nombre de degrés de liberté, ce qui augmente leur coût de calcul en comparaison du coût des méthodes d'éléments finis continus. Parmi les différentes approches variationnelles pour résoudre les problèmes aux limites, se distingue une famille particulière, basée sur l'utilisation de fonctions tests qui sont des solutions locales exactes des équations à résoudre. L'idée vient de E.Trefftz en 1926 et a depuis été largement développée et généralisée. Les méthodes variationnelles de type Trefftz-DG appliquées aux problèmes d'ondes se réduisent à des intégrales de surface, ce qui devrait contribuer à réduire les coûts de calcul.Les approches de type Trefftz ont été largement développées pour les problèmes harmoniques, mais leur utilisation pour des simulations en domaine transitoire est encore limitée. Quand elles sont appliquées dans le domaine temporel, les méthodes de Trefftz utilisent des maillages qui recouvrent le domaine espace-temps. C'est une des paraticularités de ces méthodes. En effet, les méthodes DG standards conduisent à la construction d'un système semi-discret d'équations différentielles ordinaires en temps qu'on intègre avec un schéma en temps explicite. Mais les méthodes de Trefftz-DG appliquées aux problèmes d'ondes conduisent à résoudre une matrice globale, contenant la discrétisation en espace et en temps, qui est de grande taille et creuse. Cette particularité gêne considérablement le déploiement de cette technologie pour résoudre des problèmes industriels.Dans ce travail, nous développons un environnement Trefftz-DG pour résoudre des problèmes d'ondes mécaniques, y compris les équations couplées de l'élasto-acoustique. Nous prouvons que les formulations obtenues sont bien posées et nous considérons la difficulté d'inverser la matrice globale en construisant un inverse approché obtenu à partir de la décomposition de la matrice globale en une matrice diagonale par blocs. Cette idée permet de réduire les coûts de calcul mais sa précision est limitée à de petits domaines de calcul. Etant données les limitations de la méthode, nous nous sommes intéressés au potentiel du "Tent Pitcher", en suivant les travaux récents de Gopalakrishnan et al. Il s'agit de construire un maillage espace-temps composé de macro-éléments qui peuvent être traités indépendamment en faisant une hypothèse de causalité. Nous avons obtenu des résultats préliminaires très encourageants qui illustrent bien l'intérêt du Tent Pitcher, en particulier quand il est couplé à une méthode de Trefftz-DG formulée à partir d'intégrales de surface seulement. Dans ce cas, le maillage espace-temps est composé d'éléments qui sont au plus de dimension 3. Il est aussi important de noter que ce cadre se prête à l'utilisation de pas de temps locaux ce qui est un plus pour gagner en précision avec des coûts de calcul réduits

    A Numerical Performance Analysis of Heat Recovery Ventilators with Staggered-Baffle Channels

    Get PDF
    This thesis presents a 2D numerical analysis of a cross flow compact energy recovery ventilator (HRV) with staggered baffle channels, using finite difference iterative method. The model was developed by discretization of the momentum and continuity equation into convective-diffusive terms and applying the power law scheme. Solving the system of equations required the combination of the Gauss-Seidel iterative method and the tridiagonal-matrix direct method applied to a staggered velocity and pressure grid. The Semi-Implicit method for the pressure-linked equations algorithm (SIMPLE) was chosen for coding of the program. The simulation was carried out to determine the effects of baffle height h/Dh, baffle spacing S/Dh and Reynolds number on thermal and flow performance, therefore it was necessary to vary one of the parameters within a range while all others were kept constant. The results showed that the baffle height has the greatest effect on the overall performance. At greater baffle heights the pressure drop, the Nusselt number, and heat transfer effectiveness was increased. The baffle pitch has the least effect on the overall performance; with increasing baffle pitch the pressure drop decreased while the Nusselt number had a slight increase prior to stabilization. The change of Reynolds increased the pressure drop and Nusselt number but reduced the residence time in the channel, diminishing the heat transfer effectiveness
    corecore