Search results

The purpose of this study is twofold: first, to derive a consistent model free-surface runup flow problems over deformable beds. The authors couple the nonlinear one-dimensional shallow water equations, including friction terms for the water free-surface and the two-dimensional second-order solid elastostatic equations for the bed deformation. Second, to develop a robust hybrid finite element/finite volume method for solving free-surface runup flow problems over deformable beds. The authors combine the finite volume for free-surface flows and the finite element method for bed elasticity.

The authors propose a new model for wave runup by static deformation on seabeds. The model consists of the depth-averaged shallow water system for the water free-surface coupled to the second-order elastostatic formulation for the bed deformation. At the interface between the water flow and the seabed, transfer conditions are implemented. Here, hydrostatic pressure and friction forces are considered for the elastostatic equations, whereas bathymetric forces are accounted for in the shallow water equations. As numerical solvers, the authors propose a well-balanced finite volume method for the flow system and a stabilized finite element method for elastostatics.

The developed coupled depth-averaged shallow water system and second-order solid elastostatic system is well suited for modeling wave runup by deformation on seabeds. The derived coupling conditions at the interface between the water flow and the bed topography resolve well the condition transfer between the two systems. The proposed hybrid finite volume element method is accurate and efficient for this class of models. The novel technique used for wet/dry treatment accurately captures the moving fronts in the computational domain without generating nonphysical oscillations. The presented numerical results demonstrate the high performance of the proposed methods.

Enhancing modeling and computations for wave runup problems is at an early stage in the literature, and it is a new and exciting area of research. To the best of our knowledge, solving wave runup problems by static deformation on seabeds using a hybrid finite volume element method is presented for the first time. The results of this research study, and the research methodologies, will have an important influence on a range of other scientists carrying out research in related fields.

The extrudate swell phenomenon is analysed by solving, simultaneously, the Navier‐Stokes equations along with the continuity equation by means of a finite volume method. In this work, the planar jet flows of incompressible viscous Newtonian and power‐law fluids for Reynolds numbers as high as 75 are simulated. The method uses the velocity components and pressure as the primitive variables and employs an unstructured triangular grid and triangular or polygonal control volume for each separate variable. The numerical results show good agreement with previously reported experimental and numerical results. Shear thickening results in an increase in swelling ratio, while the introduction of surface tension results in a describes in swelling ratio.

This paper aims to focus on the temporal and spatial fourth‐order finite volume discretization of the incompressible form of the Navier‐Stokes equations on structured uniform grids. The main purpose of the paper is to assess the accuracy enhancement with the inclusion of a high‐order reconstruction of the point‐wise velocity field on a fourth‐order accurate numerical scheme for the solution of the unsteady incompressible Navier‐Stokes equations.

The present finite volume method uses a fractional time‐step for decoupling velocity and pressure. A Runge‐Kutta integration scheme is implemented for integrating the momentum equation along with a polynomial interpolation and Simpson formula for space‐integration. The formulation is based on step‐by‐step de‐averaging process applied to the velocity field.

The reconstruction of the point‐wise velocity field on a higher‐order basis is essential to obtain solutions that effectively stand for a fourth‐order approximation of the point‐wise one. Results are provided for the Taylor vortex decay problem and for co‐ and counter‐rotating vortices to assess the increase in accuracy promoted by the inclusion of the high‐order de‐averaging procedure.

High‐order reconstruction of the point‐wise velocity field should be considered in high‐order finite volume methods for the solution of the unsteady incompressible form of the Navier‐Stokes equations on structured grids.

The inclusion of a high‐order reconstruction of the point‐wise velocity field is a simple and effective method of enhancing the accuracy of a finite volume code for the computational fluid dynamics analysis.

The paper develops an improved version of a fourth‐order accurate finite volume projection method with the inclusion of a high‐order reconstruction step.

Topology optimization is a method used for developing optimized geometric designs by distributing material pixels in a given design space that maximizes a chosen quantity of interest (QoI) subject to constraints. The purpose of this study is to develop a problem-agnostic automatic differentiation (AD) framework to compute sensitivities of the QoI required for density distribution-based topology optimization in an unstructured co-located cell-centered finite volume framework. Using this AD framework, the authors develop and demonstrate the topology optimization procedure for multi-dimensional steady-state heat conduction problems.

Topology optimization is performed using the well-established solid isotropic material with penalization approach. The method of moving asymptotes, a gradient-based optimization algorithm, is used to perform the optimization. The sensitivities of the QoI with respect to design variables, required for optimization algorithm, are computed using a discrete adjoint method with a novel AD library named residual automatic partial differentiator (Rapid).

Topologies that maximize or minimize relevant quantities of interest in heat conduction applications are presented. The efficacy of the technique is demonstrated using a variety of realistic heat transfer applications in both two and three dimensions, in conjugate heat transfer problems with finite conductivity ratios and in non-rectangular/non-cuboidal domains.

In contrast to most published work which has either used finite element methods or Cartesian finite volume methods for transport applications, the topology optimization procedure is developed in a general unstructured finite volume framework. This permits topology optimization for flow and heat transfer applications in complex design domains such as those encountered in industry. In addition, the Rapid library is designed to provide a problem-agnostic pathway to automatically compute all required derivatives to machine accuracy. This obviates the necessity to write new code for finding sensitivities when new physics are added or new cost functions are considered and permits general-purpose implementations of topology optimization for complex industrial applications.

This paper sets out to implement bounded high‐order (HO) resolution schemes in a hybrid finite volume/finite element method for the solution of the radiative transfer equation.

The hybrid finite volume/finite element method had formerly been developed using the step scheme, which is only first‐order accurate, for the spatial discretization. Here, several bounded HO resolution schemes, namely the MINMOD, CLAM, MUSCL and SMART schemes, formulated using the normalized variable diagram, were implemented using the deferred correction procedure.

The results obtained reveal an interaction between spatial and angular discretization errors, and show that the HO resolution schemes yield improved accuracy over the step scheme if the angular discretization error is small.

Although the HO resolution schemes reduce the spatial discretization error, they do not influence the angular discretization error. Therefore, the global error is only reduced if the angular discretization error is also small.

The use of HO resolution schemes is only effective if the angular refinement yields low‐angular discretization errors. Moreover, spatial and angular refinement should be carried out simultaneously.

The paper extends a methodology formerly developed in computational fluid dynamics, and aimed at the improvement of the solution accuracy, to the hybrid finite volume/finite element method for the solution of the radiative transfer equation.

The purpose of this paper is to develop Triple Finite Volume Method (tFVM), the author discretizes incompressible Navier-Stokes equation by tFVM, which leads to a special linear system of saddle point problem, and most computational efforts for solving the linear system are invested on the linear solver GMRES.

In this paper, by recently developed preconditioner Hermitian/Skew-Hermitian Separation (HSS) and the parallel implementation of GMRES, the author develops a quick solver, HSS-pGMRES-tFVM, for fast solving incompressible Navier-Stokes equation.

Computational results show that, the quick solver HSS-pGMRES-tFVM significantly increases the solution speed for saddle point problem from incompressible Navier-Stokes equation than the conventional solvers.

Altogether, the contribution of this paper is that the author developed the quick solver, HSS-pGMRES-tFVM, for fast solving incompressible Navier-Stokes equation.

A new finite volume scheme to solve Maxwell’s equations is presented. The approach is based on a leapfrog time scheme and a centered flux formula. This method is well suited for handling complex geometries, and therefore we can use unstructured grids. It is also able to capture the discontinuities of the electromagnetic fields through different media, without producing spurious oscillations. Owing to these properties, we can treat difficult problems, such a computing a scattered wave across complex objects. An analysis of the scheme is presented and numerical experiments are performed.

The finite volume method for radiative heat transfer calculations has been parallelized using two strategies, the angular domain decomposition and the spatial domain decomposition. In the first case each processor performs the calculations for the whole domain and for a subset of control angles, while in the second case each processor deals with all the control angles but only treats a spatial subdomain. The method is applied to three‐dimensional rectangular enclosures containing a grey emitting‐absorbing medium. The results obtained show that the number of iterations required to achieve convergence is independent of the number of processors in the angular decomposition strategy, but increases with the number of processors in the domain decomposition method. As a consequence, higher parallel efficiencies are obtained in the first case. The influence of the angular discretization, grid size and absorption coefficient of the medium on the parallel performance is also investigated.

An existing two‐dimensional finite volume technique is modified by introducing a correction term to increase the accuracy of the method to second order. It is well known that the accuracy of the finite volume method strongly depends on the order of the approximation of the flux term at the control volume (CV) faces. For highly orthotropic and anisotropic media, first order approximations produce inaccurate simulation results, which motivates the need for better estimates of the flux expression. In this article, a new approach to approximate the flux term at the CV face is presented. The discretisation involves a decomposition of the flux and an improved least squares approximation technique to calculate the derivatives of the dependent function on the CV faces for estimating both the cross diffusion term and a correction for the primary flux term. The advantage of this method is that any arbitrary unstructured mesh can be used to implement the technique without considering the shapes of the mesh elements. It was found that the numerical results well matched the available exact solution for a representative transport equation in highly orthotropic media and the benchmark solutions obtained on a fine mesh for anisotropic media. Previously proposed CV techniques are compared with the new method to highlight its accuracy for different unstructured meshes.

The purpose of this paper is to investigate the effect of numerical boundary condition implementation on local error and convergence in L2-norm of a finite volume discretization of the transient heat conduction equation subject to several boundary conditions, and for cases with volumetric heat generation, using both fully implicit and Crank-Nicolson time discretizations. The goal is to determine which combination of numerical boundary condition implementation and time discretization produces the most accurate solutions with the least computational effort.

The paper studies several benchmark cases including constant temperature, convective heating, constant heat flux, time-varying heat flux, and volumetric heating, and compares the convergence rates and local to analytical or semi-analytical solutions.

The Crank-Nicolson method coupled with second-order expression for the boundary derivatives produces the most accurate solutions on the coarsest meshes with the least computation times. The Crank-Nicolson method allows up to 16X larger time step for similar accuracy, with nearly negligible additional computational effort compared with the implicit method.

The findings can be used by researchers writing similar codes for quantitative guidance concerning the effect of various numerical boundary condition approximations for a large class of boundary condition types for two common time discretization methods.

The paper provides a comprehensive study of accuracy and convergence of the finite volume discretization for a wide range of benchmark cases and common time discretization methods.