Stokes flows are a type of fluid flow where convective forces are small in comparison with viscous forces, and momentum transport is entirely due to viscous diffusion. Besides being routinely used as benchmark test cases in numerical fluid dynamics, Stokes flows are relevant in several applications in science and engineering including porous media flow, biological flows, microfluidics, microrobotics, and hydrodynamic lubrication. The present study concerns the discretization of the equations of motion of Stokes flows in three dimensions utilizing the MINI mixed finite element, focusing on the superconvergence of the method which was investigated with numerical experiments using five purpose-made benchmark test cases with analytical solution. Despite the fact that the MINI element is only linearly convergent according to standard mixed finite element theory, a recent theoretical development proves that, for structured meshes in two dimensions, the pressure superconverges with order 1.5, as well as the linear part of the computed velocity with respect to the piecewise-linear nodal interpolation of the exact velocity. The numerical experiments documented herein suggest a more general validity of the superconvergence in pressure, possibly to unstructured tetrahedral meshes and even up to quadratic convergence which was observed with one test problem, thereby indicating that there is scope to further extend the available theoretical results on convergence.
相關內容
Several widely-used first-order saddle point optimization methods yield an identical continuous-time ordinary differential equation (ODE) to that of the Gradient Descent Ascent (GDA) method when derived naively. However, their convergence properties are very different even on simple bilinear games. We use a technique from fluid dynamics called High-Resolution Differential Equations (HRDEs) to design ODEs of several saddle point optimization methods. On bilinear games, the convergence properties of the derived HRDEs correspond to that of the starting discrete methods. Using these techniques, we show that the HRDE of Optimistic Gradient Descent Ascent (OGDA) has last-iterate convergence for general monotone variational inequalities. To our knowledge, this is the first continuous-time dynamics shown to converge for such a general setting. Moreover, we provide the rates for the best-iterate convergence of the OGDA method, relying solely on the first-order smoothness of the monotone operator.
In this paper we prove upper and lower bounds on the minimal spherical dispersion. In particular, we see that the inverse $N(\varepsilon,d)$ of the minimal spherical dispersion is, for fixed $\varepsilon>0$, up to logarithmic terms linear in the dimension $d$. We also derive upper and lower bounds on the expected dispersion for points chosen independently and uniformly at random from the Euclidean unit sphere.
The generalized Biot-Brinkman equations describe the displacement, pressures and fluxes in an elastic medium permeated by multiple viscous fluid networks and can be used to study complex poromechanical interactions in geophysics, biophysics and other engineering sciences. These equations extend on the Biot and multiple-network poroelasticity equations on the one hand and Brinkman flow models on the other hand, and as such embody a range of singular perturbation problems in realistic parameter regimes. In this paper, we introduce, theoretically analyze and numerically investigate a class of three-field finite element formulations of the generalized Biot-Brinkman equations. By introducing appropriate norms, we demonstrate that the proposed finite element discretization, as well as an associated preconditioning strategy, is robust with respect to the relevant parameter regimes. The theoretical analysis is complemented by numerical examples.
We develop a lowest-order nonconforming virtual element method for planar linear elasticity in the pure traction formulation, which can be viewed as an extension of the idea in Falk (1991) to the virtual element method, with the family of polygonal meshes satisfying a very general geometric assumption. A discrete Korn's inequality is established. And the method is shown to be uniformly convergent for the nearly incompressible case with an optimal convergence rate.
A new local discontinuous Galerkin (LDG) method for convection-diffusion equations on overlapping meshes with periodic boundary conditions was introduced in \cite{Overlap1}. With the new method, the primary variable $u$ and the auxiliary variable $p=u_x$ are solved on different meshes. In this paper, we will extend the idea to convection-diffusion equations with non-periodic boundary conditions, i.e. Neumann and Dirichlet boundary conditions. The main difference is to adjust the boundary cells. Moreover, we study the stability and suboptimal error estimates. Finally, numerical experiments are given to verify the theoretical findings.
The Virtual Element Method (VEM) is a Galerkin approximation method that extends the Finite Element Method (FEM) to polytopal meshes. In this paper, we present a conforming formulation that generalizes the Scott-Vogelius finite element method for the numerical approximation of the Stokes problem to polygonal meshes in the framework of the virtual element method. In particular, we consider a straightforward application of the virtual element approximation space for scalar elliptic problems to the vector case and approximate the pressure variable through discontinuous polynomials. We assess the effectiveness of the numerical approximation by investigating the convergence on a manufactured solution problem and a set of representative polygonal meshes. We numerically show that this formulation is convergent with optimal convergence rates except for the lowest-order case on triangular meshes, where the method coincides with the $P_1-P_0$ Scott-Vogelius scheme, and on square meshes, which are situations that are well-known to be unstable.
In communication complexity the Arthur-Merlin (AM) model is the most natural one that allows both randomness and non-determinism. Presently we do not have any super-logarithmic lower bound for the AM-complexity of an explicit function. Obtaining such a bound is a fundamental challenge to our understanding of communication phenomena. In this article we explore the gap between the known techniques and the complexity class AM. In the first part we define a new natural class, Small-advantage Layered Arthur-Merlin (SLAM), that has the following properties: - SLAM is (strictly) included in AM and includes all previously known subclasses of AM with non-trivial lower bounds. - SLAM is qualitatively stronger than the union of those classes. - SLAM is a subject to the discrepancy bound: in particular, the inner product function does not have an efficient SLAM-protocol. Structurally this can be summarised as SBP $\cup$ UAM $\subset$ SLAM $\subseteq$ AM $\cap$ PP. In the second part we ask why proving a lower bound of $\omega(\sqrt n)$ on the MA-complexity of an explicit function seems to be difficult. Both of these results are related to the notion of layer complexity, which is, informally, the number of "layers of non-determinism" used by a protocol.
This paper focuses on the strong convergence of the truncated $\theta$-Milstein method for a class of nonautonomous stochastic differential delay equations whose drift and diffusion coefficients can grow polynomially. The convergence rate, which is close to one, is given under the weaker assumption than the monotone condition. To verify our theoretical findings, we present a numerical example.
In a companion study \cite{patterson2020computing2D}, we present a numerical method for simulating 2D viscous flow through an open compliant closed channel, drive by pressure gradient. We consider the highly viscous regime, where fluid dynamics is described by the Stokes equations, and the less viscous regime described by the Navier-Stokes equations. In this study, we extend the method to 3D tubular flow. The problem is formulated in axisymmetric cylindrical coordinates, an approach that is natural for tubular flow simulations and that substantially reduces computational cost. When the elastic tubular walls are stretched or compressed, they exert forces on the fluid. These singular forces introduce unsmoothness into the fluid solution. As in the companion 2D study \cite{patterson2020computing2D}, we extend the immersed interface method to an open tube, and we compute solution to the model equations using the resulting method. Numerical results indicate that this new method preserves sharp jumps in the solution and its derivatives, and converges with second-order accuracy in both space and time.
In this work, we consider the distributed optimization of non-smooth convex functions using a network of computing units. We investigate this problem under two regularity assumptions: (1) the Lipschitz continuity of the global objective function, and (2) the Lipschitz continuity of local individual functions. Under the local regularity assumption, we provide the first optimal first-order decentralized algorithm called multi-step primal-dual (MSPD) and its corresponding optimal convergence rate. A notable aspect of this result is that, for non-smooth functions, while the dominant term of the error is in $O(1/\sqrt{t})$, the structure of the communication network only impacts a second-order term in $O(1/t)$, where $t$ is time. In other words, the error due to limits in communication resources decreases at a fast rate even in the case of non-strongly-convex objective functions. Under the global regularity assumption, we provide a simple yet efficient algorithm called distributed randomized smoothing (DRS) based on a local smoothing of the objective function, and show that DRS is within a $d^{1/4}$ multiplicative factor of the optimal convergence rate, where $d$ is the underlying dimension.
注冊地址： 北京市海淀區羊坊店路18號2幢3層301-191