Numerical three-dimensional free surface circulation model for the South Biscayne Bay,Florida

A three– dimensional, time dependent free surface model has been developed for predicting circulation and surface height variations in a tidal bay. An explicit finite difference numerical solution is obtained by transforming the vertical coordinate in the governing model equations. The transformed geometry consists of a fixed, flat free surface and a constant depth basin for easy computation. The ocean–bay interface open boundary condition is incorporated into this hydrodynamic model without approximation, and yields rather accurate results for the bay circulation and tide level variations. The numerical method employs a staggered grid Richardson lattice, which has the inherent property of not requiring calculation of the tangential velocity components on solid surfaces. The momentum equations ignore horizontal diffusion which is small for the South Biscayne Bay, for which vertical diffusion and advection dominate.     

Continuous source in tidal flow: A numerical study of the transport equation

A detailed comparative numerical study of a continuous source discharging in a tidal flow has been performed. The advective term in the one-dimensional transport equation consists of a constant freshwater velocity and an isotropic oscillating component. Various finite element solutions are investigated. Compared to the analytical solution of the dynamic steady state concentration, the numerical results for typical estuarine conditions clearly indicate the superiority of the collocation method with Hermite basis functions. An extended Fourier series analysis that accounts for the source condition is developed to explain the numerical behaviour of the different schemes.     

Stabilized Crouzeix–Raviart element for Darcy–Forchheimer model

In this article, a stabilized mixed finite element method for steady Darcy–Forchheimer flow is introduced, in which the velocity and pressure are approximated by nonconforming Crouzeix–Raviart element and piecewise constant, respectively. A discrete inf‐sup condition and a priori error estimates are derived. An iterative scheme is given for practical computation. Finally, some numerical examples are carried out to verify the theoretical analysis and a comparison between two discretizations is given to demonstrate that one of the discretizations has better properties. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1568–1588, 2015     

The finite volume method based on the Crouzeix–Raviart element for a fracture model

In this paper, we present and analyze a finite volume method based on the Crouzeix–Raviart element for the coupled fracture model, where the fluid flow is governed by Darcy's law in the one‐dimensional fracture and two‐dimensional surrounding matrix. In the numerical scheme, the pressure in the matrix and fracture is respectively approximated by the Crouzeix–Raviart elements and piecewise constant functions, and then the velocity is calculated by piecewise constant functions element by element. The existence and uniqueness of the numerical solution are discussed, and optimal order error estimates for both the pressure p and the velocity u are proved on general triangulations. We finally carry out numerical experiments, and results confirm our theoretical analysis.     

A convergent FEM-DG method for the compressible Navier–Stokes equations

This paper presents a new numerical method for the compressible Navier–Stokes equations governing the flow of an ideal isentropic gas. To approximate the continuity equation, the method utilizes a discontinuous Galerkin discretization on piecewise constants and a basic upwind flux. For the momentum equation, the method is a new combined discontinuous Galerkin and finite element method approximating the velocity in the Crouzeix–Raviart finite element space. While the diffusion operator is discretized in a standard fashion, the convection and time-derivative are discretized using discontinuous Galerkin on the element average velocity and a Lax–Friedrich type flux. Our main result is convergence of the method to a global weak solution as discretization parameters go to zero. The convergence analysis constitutes a numerical version of the existence analysis of Lions and Feireisl.     

A simple finite element method for the Stokes equations

The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.     

曲率流的参数化有限元逼近

许多物理现象可以在数学上描述为受曲率驱动的自由界面运动,例如薄膜和泡沫的演变、晶体生长,等等.这些薄膜和界面的运动常依赖于其表面曲率,从而可以用相应的曲率流来描述,其相关自由界面问题的数值计算和误差分析一直是计算数学领域中的难点.参数化有限元法是曲率流的一类有效计算方法,已经能够成功模拟一些曲面在几类基本的曲率流下的演化过程.本文重点讨论曲率流的参数化有限元逼近,它的产生、发展和当前的一些挑战.     

粘弹性方程全离散化有限体积元格式及数值模拟

本文利用有限体积元方法研究二维粘弹性方程, 给出一种时间二阶精度的全离散化有限体积元格式, 并给出这种全离散化有限体积元解的误差估计, 最后用数值例子验证数值结果与理论结果是相吻合的. 通过与有限元方法和有限差分方法相比较, 进一步说明了全离散化有限体积元格式是求解二维粘弹性方程数值解的最有效方法之一.     

A Mixed Analytical/Numerical Method for Velocity and Heat Transfer of Laminar Power-Law Fluids

This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.     

解非线性有限元方程的逐层校正迭代法

本文提出一种求解非线性有限元方程的逐层校正迭代法.有关数值分析表明,当网格分划较细,网格分划参数h_j较小时,仅需一次简单的迭代和校正步骤就可满足数值计算的要求,使用该方法的计算复杂性是最佳阶的,即为O(N_j),其中N_j为最细网格层上离散结点变量的数目.     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A finite volume element formulation and error analysis for the non-stationary conduction–convection problem

The non-stationary conduction–convection problem including the velocity vector field and the pressure field as well as the temperature field is studied with a finite volume element (FVE) method. A fully discrete FVE formulation and the error estimates between the fully discrete FVE solutions and the accuracy solution are provided. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the FVE method is feasible and efficient for finding the numerical solutions of the non-stationary conduction–convection problem and is one of the most effective numerical methods by comparing the results of the numerical simulations of the FVE formulation with those of the numerical simulations of the finite element method and the finite difference scheme for the non-stationary conduction–convection problem.     

A New Reduced-Order fve Algorithm Based on POD Method for Viscoelastic Equations

A proper orthogonal decomposition (POD) technique is used to reduce the finite volume element (FVE) method for two-dimensional (2D) viscoelastic equations. A reduced-order fully discrete FVE algorithm with fewer degrees of freedom and sufficiently high accuracy based on POD method is established. The error estimates of the reduced-order fully discrete FVE solutions and the implementation for solving the reduced-order fully discrete FVE algorithm are provided. Some numerical examples are used to illustrate that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced-order fully discrete FVE algorithm is one of the most effective numerical methods by comparing with corresponding numerical results of finite element formulation and finite difference scheme and that the reduced-order fully discrete FVE algorithm based on POD method is feasible and efficient for solving 2D viscoelastic equations.     

非定常不可压Navier-Stokes方程两重网格方法的稳定性和L2误差估计

讨论了二维非定常不可压Navier-Stokes方程的两重网格方法．此方法包括在粗网格上求解一个非线性问题,在细网格上求解一个Stokes问题．采用一种新的全离散(时间离散用Crank-Nicolson格式,空间离散用混合有限元方法)格式数值求解N-S方程．证明了该全离散格式的稳定性．给出了L2误差估计．对比标准有限元方法,在保持同样精度的前提下,TGM能节省大量的计算量．     

Mixed convection in the stagnation-point flow of a Maxwell fluid towards a vertical stretching surface

In the present analysis, we study the steady mixed convection boundary layer flow of an incompressible Maxwell fluid near the two-dimensional stagnation-point flow over a vertical stretching surface. It is assumed that the stretching velocity and the surface temperature vary linearly with the distance from the stagnation-point. The governing nonlinear partial differential equations have been reduced to the coupled nonlinear ordinary differential equations by the similarity transformations. Analytical and numerical solutions of the derived system of equations are developed. The homotopy analysis method (HAM) and finite difference scheme are employed in constructing the analytical and numerical solutions, respectively. Comparison between the analytical and numerical solutions is given and found to be in excellent agreement. Both cases of assisting and opposing flows are considered. The influence of the various interesting parameters on the flow and heat transfer is analyzed and discussed through graphs in detail. The values of the local Nusselt number for different physical parameters are also tabulated. Comparison of the present results with known numerical results of viscous fluid is shown and a good agreement is observed.     

Assessment of subgrid-scale models for the incompressible Navier-Stokes equations

In this paper, we assess two kinds of subgrid finite element methods for the two-dimensional (2D) incompressible Naver-Stokes equations (NSEs). These methods introduce subgrid-scale (SGS) eddy viscosity terms which do not act on the large flow structures. The eddy viscous terms consist of the fluid flow fluctuation strain rate stress tensors. The fluctuation tensor can be calculated by a elliptic projection or a simple L² projection (projective filter) in finite element spaces. The finite element pair P₂/P₁ is adopted to numerically implement analysis and computation. We give a complete error analysis based on the assumptions of some regularity conditions. On the part of numerical tests, the numerical computations for the stationary flows show that the numerical results agree with some benchmark solutions and theoretical analysis very well. Furthermore, the given SGS models are applied to the non-stationary fluid flows.     

Finite element analysis of flood discharge atomization based on water–air two-phase flow

Flood discharge atomization is a phenomenon of water fog diffusion caused by the discharge of water from a spillway structure, which brings strong wind and heavy rainfall. These unnatural winds and rainfall are harmful for the safe operation of hydropower stations with high water heads. Compared to the method of prototype observations, physical models and mathematical models, which are semi-theoretical and semi-empirical, numerical simulation methods have the advantage of being not limited by a similar scale and are more economical. A finite element model is presented to simulate flood discharge atomization based on water–air two-phase flow in this paper. Equations governing flood discharge atomization are composed of partial differential equations of mass and momentum conservation laws with unknowns for pressure, velocity and the water concentration. The finite element method is used to solve the governing equations by adopting appropriate solution strategies to increase the convergence and numerical stability. Then, the finite element model is applied to a practical project, the Shuibuya hydropower station, which experienced a flood discharge in 2016. Simulation results show that the proposed model can simulate flood discharge atomization with efficient convergence and numerical stability in three dimensions, and good agreement was observed between numerical simulations and prototype observational data. Based on the simulation results, the mechanism of flood discharge atomization was analyzed.     

A Least‐Squares Formulation for the Shallow Water Equations with Small Viscosity

A least‐squares finite element method for the shallow water equations with viscosity parameter μ > 0 is proposed and studied. The shallow water equations are reformulated as a first order system by adding a new variable for the velocity flux. The reformulated first order system is combined with a characteristic‐based time discretization and a least squares approach. For the correct boundary treatment in the limit case μ → 0, a trace theorem is presented. For the numerical computation of the velocity, the finite element spaces introduced recently by Mardal, Tai and Winther (SIAM Journal on Numerical Analysis 40, pp. 1605–1631) are used. The degrees of freedom in these spaces can be identified with the normal and tangential components, respectively. Numerical results for some test examples are shown. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

轴对称非线性磁场分布的有限元分析

一、引言 轴对称电磁结构在轴对称电流作用下的磁场分布问题本是三维问题,利用轴对称关系可将三维麦克斯韦方程化为形如下面(1)式的二维非线性椭圆型方程.它与标准的平面磁场方程有所不同.而且v(B)不是B的单增函数.实际工程中的大型圆柱变压器、平     

油气资源数值模拟的变网格交替方向特征有限元格式和分析

石油科学在有机地球化学、石油的生成、运移、聚集研究取得重大进展,在评价油气资源时,对盆地发育史尤其对流体流动规律和受热变化历史的计算是非常重要的.其数学模型是三维空间非线性偶合偏微分方程组的初边值问题.从实际出发,考虑了流体的压缩性和三维问题大规模科学与工程计算的特征, 提出了一类变网格交替方向特征有限元格式,应用变分形式、算子分裂、广义$L^2$投影、能量方法、负模估计、微分方程先验估计的理论和技巧,得到最佳阶$L^2$误差估计. 此方法已成功应用到油气资源评估数值模拟的生产实践中,成功解决了这一重要问题.