Application of mesh-adaptation for pollutant transport by water flow

An adaptive finite volume method is proposed for the numerical solution of pollutant transport by water flows. The shallow water equations with eddy viscosity, bottom friction forces and wind shear stresses are used for modelling the water flow whereas, a transport-diffusion equation is used for modelling the advection and dispersion of pollutant concentration. The adaptive finite volume method uses simple centred-type discretization for the source terms, can handle complex topography using unstructured grids and satisfies the conservation property. The adaptation criteria are based on monitoring the pollutant concentration in the computational domain during its dispersion process. The emphasis in this paper is on the application of the proposed method for numerical simulation of pollution dispersion in the Strait of Gibraltar. Results are presented using different tidal conditions and wind-induced flow fields in the Strait.     

Coupling of mesoscale and microscale models—an approach to simulate scale interaction

Atmospheric flow and pollutant dispersion over built-up areas are affected by phenomena occurring at different scales. Hence, scale interactions should also be considered in the mathematical modelling of atmospheric flow and pollutant dispersion. In this paper a method is presented to couple prognostic mesoscale and microscale flow models. Results from mesoscale simulations are used to generate the initial state and boundary conditions for microscale simulations. The method comprises a three-dimensional interpolation scheme and a vertical adjustment of the interpolated quantities in the surface layer based on similarity theory. The method is applied to couple the microscale model MIMO with the mesoscale model MEMO.The coupled system MEMO–MIMO is applied to simulate the local scale flow for an industrial area in southwestern Germany. Model results are presented and compared with available measurements.     

A new finite volume method for flux-gradient and source-term balancing in shallow water equations

We present a new finite volume method for flux-gradient and source-term balancing in the numerical solution of shallow water equations on nonflat topography. The method consists of a predictor stage for discretization of gradient terms and a corrector stage for treatment of source terms. The numerical fluxes at the interfaces of each triangle are reconstructed using an upwind scheme along with slope limiters to obtain high accuracy in the method. Using a three-sided decomposition, the source terms are balanced in the corrector stage. The proposed method is implemented on unstructured meshes and verifies the exact conservation property (C-property). We also discuss a mesh adaptive procedure by monitoring the water depth in the computational domain. Several standard test examples are used to verify high accuracy, exact C-property and good resolution properties for smooth and discontinuous solutions.     

Large eddy simulation of turbulent heat transport in the Strait of Gibraltar

We develop a numerical model for large eddy simulation of turbulent heat transport in the Strait of Gibraltar. The flow equations are the incompressible Navier–Stokes equations including Coriolis forces and density variation through the Boussinesq approximation. The turbulence effects are incorporated in the system by considering the Smagorinsky model. As a numerical solver we propose a finite element semi-Lagrangian method. The solution procedure consists of combining a non-oscillatory semi-Lagrangian scheme for time discretization with the finite element method for space discretization. Numerical results illustrate a buoyancy-driven circulations along the Strait of Gibraltar and the sea-surface temperature is flushed out and move to northeast coast. The Ocean discharge and the temperature difference are shown to control the plume structure.     

An adaptive solver for the spherical shallow water equations

The shallow water equations (SWE), which describe the flow of a thin layer of fluid in two dimensions have been used by the atmospheric modelling community as a vehicle for testing promising numerical methods for solving atmospheric and oceanic problems. The SWE are important for the study of the dynamics of large-scale flows, as well for the development of new numerical schemes that are applied to more complex models. In this paper we present a finite difference p-adaptive method based on high order finite differences that is applied using an error indicator for solving the SWE on the sphere. A standard test set is used to evaluate the accuracy of the new method. The results obtained are compared with the pseudo-spectral method.     

基于物理的波浪实时模拟方法

为解决广域范围内波浪状态的实时计算和可视化问题,结合无粘滞流体力学的物理模型,提出了一种以三角形为控制单元的有限体积简化算法。该算法的优点在于：以不规则边界区域的非规则三角网格为基础,通过简化通量向量分裂方法获得三角形控制单元的边界数值通量,能快速逼近二维浅水方程的解进而模拟非规则边界浅水的实时流动。实验结果显示,所提方法能在符合现实世界物理规律的前提下较好地实现大规模波浪的实时可视化模拟。     

应用于地下水模拟的多尺度有限体积元方法

应用多尺度有限体积元方法模拟地下水流动问题,其中地下渗透场系数采用二维对数正态随机场.与传统的有限体积元法相比,多尺度有限体积元法的基函数具有能够反映单元内参数变化的优点,所以这种方法能在大尺度上捕捉解的小尺度特征获得较精确的解.文中算例分别对均匀、各向同性和各向异性对数正态随机场的二维地下水流动问题用传统数值模拟方法和多尺度有限体积元方法进行了计算.计算结果表明多尺度有限体积元方法收敛,且与传统数值模拟方法相比,多尺度有限体积元方法既节省计算量,又有较高的精度.     

Comparison of numerical techniques for 3D flutter analysis of cable-stayed bridges

This paper presents the results of a comparison study of the numerical techniques of structural and aerodynamic force models developed based on the spline finite strip method with the conventional finite element approach in three-dimensional flutter analysis of cable-stayed bridges. In the new formulation, the bridge girder is modelled by spline finite strips. The mass and stiffness properties of the torsional behaviour of complex bridge girder, which have a significant influence on the wind stability of long-span bridges, are modelled accurately in the formulation. The effects of the spatial variation of the aerodynamic forces can be taken into account in the proposed numerical model by distributing the loads to the finite strips modelling the bridge deck. The numerical example of a 423 m long-span cable-stayed bridge is presented in the comparison study. The accuracy and effectiveness of the proposed finite strip model are compared to the results obtained from the equivalent beam finite element models. The advantages and disadvantages of these different modelling schemes are discussed.     

Numerical modelling of interfacial soil erosion with viscous incompressible flows

The aim of this study is to develop a numerical model for simulating surface erosion occurring at a fluid/soil interface subject to a flow process. Balance equations with jump relations are used. A penalization procedure including a fictitious domain method is used to compute the Stokes flow around obstacles, in order to avoid body-fitted unstructured meshes and instead use fast and efficient finite volume approximations on Cartesian meshes. The evolution of the water/soil interface is described by using a level set function. The ability of the model to predict the interfacial erosion of soils is confirmed by several numerical simulations.     

Adaptive finite volume methods for displacement problems in porous media

In this paper we consider adaptive numerical simulation of miscible and immiscible displacement problems in porous media, which are modeled by single and two phase flow equations. Using the IMPES formulation of the two phase flow equation both problems can be treated in the same numerical framework. We discretise the equations by an operator splitting technique where the flow equation is approximated by Raviart-Thomas mixed finite elements and the saturation or concentration equation by vertex centered finite volume methods. Using a posteriori error estimates for both approximation schemes we deduce an adaptive solution algorithm for the system of equations and show the applicability in several examples.     

Multi-level adaptive simulation of transient two-phase flow in heterogeneous porous media

An implicit pressure and explicit saturation (IMPES) finite element method (FEM) incorporating a multi-level shock-type adaptive refinement technique is presented and applied to investigate transient two-phase flow in porous media. Local adaptive mesh refinement is implemented seamlessly with state-of-the-art artificial diffusion stabilization allowing simulations that achieve both high resolution and high accuracy. Two benchmark problems, modelling a single crack and a random porous medium, are used to demonstrate the robustness of the method and illustrate the capabilities of the adaptive refinement technique in resolving the saturation field and the complex interaction (transport phenomena) between two fluids in heterogeneous media.     

Modeling of indoor airflow and dispersion of aerosols using immersed boundary and random flow generation methods

Numerical modeling of environmental flows involves complex geometry, moving bodies, multi-phase flow, and buoyant jet effects. An in-house CFD code has been developed using finite volume and immersed boundary methods. The transport and dispersion of virus-laden aerosols in a ventilated room is investigated by this numerical code. The uniqueness of this numerical code is that it can efficiently compute small-scale turbulent flow in which most commercial CFD software will suffer from large numerical error. Random flow generation (RFG) [33] is an ideal choice for small-scale turbulence for low-Reynolds number flow in a ventilated room. In addition, Lagrangian stochastic (LS) walk model is applied to directly compute probability density function (PDF) of aerosols and estimate the risk factor of aerosol dispersion from a point source. The present study focuses on aerosols with small diameter (<10 μm) in which the effects of evaporation on the dispersion of aerosols could be neglected. Different location of aerosol sources and a typical ventilation layout are discussed in detail. The numerical results with PDF yield more useful quantitative information to assess the risk area of virus transport in a ventilated room than that shown in random trajectories of particles as widely reported in the literature. This study provides valuable information for ventilation control strategies with respiratory protection, such as enhanced air exchange, air filtration rate, and improved airflow patterns to reduce indoor infection risk via airborne virus laden droplets.     

A 3rd order upwind finite volume method for generalised Newtonian fluid flows

The non-Newtonian effects in the flow of non-Newtonian fluids that can be modelled with generalised Newtonian constitutive equations are investigated using a numerical scheme based on the finite volume formulation. Collocated arrangement of variables is used and the pressure-correction method in conjunction with the SIMPLE scheme is applied. In order to avoid the diffusive effects of low order schemes, the approximation of the convection terms is carried out using the QUICK differencing scheme, which is of third order of accuracy. For modelling viscous non-Newtonian behaviour, the Power-Law, Quemada and the modified Bingham and Casson models are employed. Validation of the code is carried out upon comparison with numerical data available in the literature. Exploiting the lid-driven cavity flow a twofold investigation is carried out regarding the non-Newtonian effects first by using different models and secondly by using the shear-thinning and shear-thickening attributes of the fluid.     

A stabilized finite element method for stochastic incompressible Navier–Stokes equations

The major emphasis of this work is the development of a stabilized finite element method for solving incompressible Navier–Stokes equations with stochastic input data. The polynomial chaos expansion is used to represent stochastic processes in the variational problem, resulting in a set of deterministic variational problems to be solved for each Wiener polynomial chaos. To obtain the chaos coefficients in the corresponding deterministic incompressible Navier–Stokes equations, we combine the modified method of characteristics with the finite element discretization. The obtained Stokes problem is solved using a robust conjugate-gradient algorithm. This algorithm avoids projection procedures and any special correction for the pressure. These numerical techniques associate the geometrical flexibility of the finite element method with the ability offered by the modified method of characteristics to solve convection-dominated problems using time steps larger than its Eulerian counterpart. Numerical results are shown for the benchmark problems of driven cavity flow and backward-facing step flow. We also present numerical results for a problem of stochastic natural convection. It is found that the proposed stabilized finite element method offers a robust and accurate approach for solving the stochastic incompressible Navier–Stokes equations, even when high Reynolds and Rayleigh numbers are used in the simulations.     

Large eddy simulation of turbulent heat transport in the Strait of Gibraltar

We develop a numerical model for large eddy simulation of turbulent heat transport in the Strait of Gibraltar. The flow equations are the incompressible Navier–Stokes equations including Coriolis forces and density variation through the Boussinesq approximation. The turbulence effects are incorporated in the system by considering the Smagorinsky model. As a numerical solver we propose a finite element semi-Lagrangian method. The solution procedure consists of combining a non-oscillatory semi-Lagrangian scheme for time discretization with the finite element method for space discretization. Numerical results illustrate a buoyancy-driven circulations along the Strait of Gibraltar and the sea-surface temperature is flushed out and move to northeast coast. The Ocean discharge and the temperature difference are shown to control the plume structure.     

Grid refinement tests of a least-squares finite element method for advective transport of reactive species

A time splitting least-squares finite element method (LSFEM) and the ‘stiff ODEs’ solver LSODE are used to simulate the advective transport of reactive species. Specifically, the rotating cone problem with chemical reactions serves as a model to test the algorithm. A non-linear filter is used to suppress spurious oscillations at each advective time step. All simulations are carried out by using linear and quadratic elements on two mesh systems. Results from the coarse mesh system suggest that the use of robust numerical methods alone will not be able to provide accurate results and point to the need of grid refinement. The grid refinement tests are evaluated by pollutant peak concentrations, pollutant concentration distributions, average relative errors, species mass conservation and distribution ratios and CPU times. Results from the fine mesh system are satisfactory and imply that accurate simulations of the transport of reactive species require adequate grid resolution and robust numerical methods for each individual advective and reactive step.     

基于自适应偏心模型的柱塞副水膜特性数值分析

柱塞副间隙形状对柱塞副的润滑特性影响很大,为提高径向柱塞泵柱塞副工作效率和可靠性,需要准确仿真柱塞副的水膜形状和压力分布。首先,建立了基于等效液阻高度和柱塞自适应偏心的柱塞副间隙泄露模型。模型考虑了柱塞在缸体内的微运动和柱塞受力平衡的因素,运用遗传算法和有限体积法对柱塞副间隙形状和压力分布进行数值仿真分析,进而讨论柱塞副间隙泄露流量。最后,从出口压力、电机转速以及初始密封间隙等方面对压力分布特性和泄露流量进行分析,深入探讨了柱塞副水膜压力分布和泄露流量的影响因素及变化规律。试验结果表明,考虑柱塞自适应偏心的模型对柱塞泄露流量的模拟准确,获取的实时柱塞偏心状态和压力分布与实际状态更为接近。     

Finite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water

We present a new hybrid numerical method for computing the transport of a passive pollutant by a flow. The flow is modeled by the Saint-Venant system of shallow water equations and the pollutant propagation is described by a transport equation. The idea behind the new finite-volume-particle (FVP) method is to use different schemes for the flow and the pollution computations: the shallow water equations are numerically integrated using a finite-volume scheme, while the transport equation is solved by a particle method. This way the specific advantages of each scheme are utilized at the right place. This results in a significantly enhanced resolution of the computed solution     

A low-dissipation DG method for the under-resolved simulation of low Mach number turbulent flows

In recent years the use of high-order Discontinuous Galerkin (DG) methods for the under-resolved direct numerical simulations (uDNS) of turbulent flows has received special attention. The suitability of the approach for this kind of applications is related to the dissipation and dispersion proprieties of the scheme: while the dispersion errors are small over a broad range of frequencies, a relevant dissipation error mainly acts at the smallest under-resolved scales, resembling a high frequency filter. Nevertheless, it was recognized (Flad and Gassner, 2017; Mengaldo et al., 2018) that the choice of the interface convective numerical flux strongly affects this dissipation behaviour and ultimately the success of the uDNS approach. In this regard, the excess of numerical dissipation caused by some upwind numerical convective fluxes must be avoided, in particular when dealing with low-speed flows, since this behaviour is exacerbated approaching the incompressibility limit. Fixes for the excess of numerical dissipation of these schemes have been proposed by several authors in the context of different numerical methods, see for example Weiss and Smith (1995). In this work a simple modification of the dissipation term of the low Mach preconditioned Roe scheme proposed by Weiss and Smith is considered. The aim is to reduce further the amount of numerical dissipation with the intent of improving the results of uDNS. A spatial DG discretization coupled with a linearly-implicit Rosenbrock-type time integrator is here considered as a numerical framework perfectly suited for the assessment and comparison of different numerical flux functions. Results on canonical turbulent flow problems as the Taylor–Green vortex and the flow in a straight sided channel are presented. The improved accuracy of the proposed flux function is demonstrated. The new low-dissipation flux can be useful also in the context of standard, lower order, finite volume methods.     

Well-balanced RKDG2 solutions to the shallow water equations over irregular domains with wetting and drying

This paper presents a new one-dimensional (1D) second-order Runge–Kutta discontinuous Galerkin (RKDG2) scheme for shallow flow simulations involving wetting and drying over complex domain topography. The shallow water equations that adopt water level (instead of water depth) as a flow variable are solved by an RKDG2 scheme to give piecewise linear approximate solutions, which are locally defined by an average coefficient and a slope coefficient. A wetting and drying technique proposed originally for a finite volume MUSCL scheme is revised and implemented in the RKDG2 solver. Extra numerical enhancements are proposed to amend the local coefficients associated with water level and bed elevation in order to maintain the well-balanced property of the RKDG2 scheme for applications with wetting and drying. Friction source terms are included and evaluated using splitting implicit discretization, implemented with a physical stopping condition to ensure stability. Several steady and unsteady benchmark tests with/without friction effects are considered to demonstrate the performance of the present model.