-
A Test Problem for Flow Codes
Authors:
Henry von Wahl,
L. Ridgway Scott
Abstract:
We propose a test problem for Navier-Stokes solvers based on the flow around a cylinder. We choose a range of Reynolds numbers for which the flow is time-dependent but can be characterized as essentially two-dimensional. The test problem requires accurate resolution of chaotic dynamics over a long time interval. It also requires the use of a relatively large computational domain, part of which is…
▽ More
We propose a test problem for Navier-Stokes solvers based on the flow around a cylinder. We choose a range of Reynolds numbers for which the flow is time-dependent but can be characterized as essentially two-dimensional. The test problem requires accurate resolution of chaotic dynamics over a long time interval. It also requires the use of a relatively large computational domain, part of which is curved, and it requires evaluation of derivatives of the solution and pressure on the curved boundary. We review the performance of different finite element methods for the proposed range of Reynolds numbers. These tests indicate that some of the most established methods do not capture the correct behavior.
△ Less
Submitted 26 April, 2024; v1 submitted 25 April, 2024;
originally announced April 2024.
-
Computational analysis of a contraction rheometer for the grade-two fluid model
Authors:
Sara Pollock,
L. Ridgway Scott
Abstract:
We explore the possibility of simulating the grade-two fluid model in a geometry related to a contraction rheometer, and we provide details on several key aspects of the computation. We show how the results can be used to determine the viscosity $ν$ from experimental data. We also explore the identifiability of the grade-two parameters $α_1$ and $α_2$ from experimental data. In particular, as the…
▽ More
We explore the possibility of simulating the grade-two fluid model in a geometry related to a contraction rheometer, and we provide details on several key aspects of the computation. We show how the results can be used to determine the viscosity $ν$ from experimental data. We also explore the identifiability of the grade-two parameters $α_1$ and $α_2$ from experimental data. In particular, as the flow rate varies, force data appears to be nearly the same for certain distinct pairs of values $α_1$ and $α_2$; however we determine a regime for $α_1$ and $α_2$ for which the parameters may be identifiable with a contraction rheometer.
△ Less
Submitted 4 April, 2024;
originally announced April 2024.
-
Chaotic dynamics of two-dimensional flows around a cylinder
Authors:
L. Ridgway Scott,
Rebecca Durst
Abstract:
We study flow around a cylinder from a dynamics perspective, using drag and lift as indicators. We observe that the mean drag coefficient bifurcates from the steady case when the Karman vortex street emerges. We also find a jump in the dimension of the drag/lift attractor just above Reynolds number 100. We compare the simulated drag values with experimental data obtained over the last hundred year…
▽ More
We study flow around a cylinder from a dynamics perspective, using drag and lift as indicators. We observe that the mean drag coefficient bifurcates from the steady case when the Karman vortex street emerges. We also find a jump in the dimension of the drag/lift attractor just above Reynolds number 100. We compare the simulated drag values with experimental data obtained over the last hundred years. Our simulations suggest that a vibrational resonance in the cylinder would be unlikely for Reynolds numbers greater than 1000, where the drag/lift behavior is fully chaotic.
△ Less
Submitted 13 November, 2023;
originally announced November 2023.
-
Verification and Validation of Cylinder Drag: Pressure and Stress Approximations on Curved Boundaries
Authors:
Ingeborg G. Gjerde,
L. Ridgway Scott
Abstract:
We study a technique for verification of stress and pressure computations on boundaries in flow simulations. We utilize existing experiments to provide validation of the simulations. We show that this approach can reveal critical flaws in simulation algorithms. Using the successful computational algorithms, we examine Lamb's model for cylinder drag at low Reynolds numbers. We comment on a discrepa…
▽ More
We study a technique for verification of stress and pressure computations on boundaries in flow simulations. We utilize existing experiments to provide validation of the simulations. We show that this approach can reveal critical flaws in simulation algorithms. Using the successful computational algorithms, we examine Lamb's model for cylinder drag at low Reynolds numbers. We comment on a discrepancy observed in an experimental paper, suggesting that the domain size may be a contributing factor. Our simulations on suitably large domains confirm Lamb's model. We highlight a paradox related to imposing Dirichlet (Stokes) boundary conditions on polygonal approximations of the curved surface using finite-element methods that are exactly divergence free. The finite-element simulations provide very poor representations of drag when the boundary conditions are imposed strongly. We demonstrate that relaxing the boundary conditions using Nitsche's method restores high-order approximation.
△ Less
Submitted 21 June, 2023;
originally announced June 2023.
-
Dimensions of exactly divergence-free finite element spaces in 3D
Authors:
L. Ridgway Scott,
Tabea Tscherpel
Abstract:
We examine the dimensions of various inf-sup stable mixed finite element spaces on tetrahedral meshes in 3D with exact divergence constraints. More precisely, we compare the standard Scott-Vogelius elements of higher polynomial degree and low order methods on split meshes, the Alfeld and the Worsey-Farin split. The main tool is a counting strategy to express the degrees of freedom for given polyno…
▽ More
We examine the dimensions of various inf-sup stable mixed finite element spaces on tetrahedral meshes in 3D with exact divergence constraints. More precisely, we compare the standard Scott-Vogelius elements of higher polynomial degree and low order methods on split meshes, the Alfeld and the Worsey-Farin split. The main tool is a counting strategy to express the degrees of freedom for given polynomial degree and given split in terms of few mesh quantities, for which bounds and asymptotic behavior under mesh refinement is investigated. Furthermore, this is used to obtain insights on potential precursor spaces in full de Rham complexes for finite element methods on the Worsey-Farin split.
△ Less
Submitted 18 April, 2024; v1 submitted 31 December, 2022;
originally announced January 2023.
-
Two conjectures on the Stokes complex in three dimensions on Freudenthal meshes
Authors:
Patrick E. Farrell,
Lawrence Mitchell,
L. Ridgway Scott
Abstract:
In recent years a great deal of attention has been paid to discretizations of the incompressible Stokes equations that exactly preserve the incompressibility constraint. These are of substantial interest because these discretizations are pressure-robust, i.e. the error estimates for the velocity do not depend on the error in the pressure. Similar considerations arise in nearly incompressible linea…
▽ More
In recent years a great deal of attention has been paid to discretizations of the incompressible Stokes equations that exactly preserve the incompressibility constraint. These are of substantial interest because these discretizations are pressure-robust, i.e. the error estimates for the velocity do not depend on the error in the pressure. Similar considerations arise in nearly incompressible linear elastic solids. Conforming discretizations with this property are now well understood in two dimensions, but remain poorly understood in three dimensions. In this work we state two conjectures on this subject. The first is that the Scott-Vogelius element pair is inf-sup stable on uniform meshes for velocity degree $k \ge 4$; the best result available in the literature is for $k \ge 6$. The second is that there exists a stable space decomposition of the kernel of the divergence for $k \ge 5$. We present numerical evidence supporting our conjectures.
△ Less
Submitted 11 January, 2024; v1 submitted 10 November, 2022;
originally announced November 2022.
-
Resolution of d'Alembert's paradox using slip boundary conditions: The effect of the friction parameter on the drag coefficient
Authors:
Ingeborg G. Gjerde,
L. Ridgway Scott
Abstract:
d'Alembert's paradox is the contradictory observation that for incompressible and inviscid (potential) fluid flow, there is no drag force experienced by a body moving with constant velocity relative to the fluid. This paradox can be straightforwardly resolved by considering Navier's slip boundary condition. Potential flow around a cylinder then solves the Navier--Stokes equations using friction pa…
▽ More
d'Alembert's paradox is the contradictory observation that for incompressible and inviscid (potential) fluid flow, there is no drag force experienced by a body moving with constant velocity relative to the fluid. This paradox can be straightforwardly resolved by considering Navier's slip boundary condition. Potential flow around a cylinder then solves the Navier--Stokes equations using friction parameter $β=-2ν$. This negative friction parameter can be interpreted physically as the fluid being accelerated by the cylinder wall. This explains the lack of drag.
In this paper, we introduce the Navier slip boundary condition and show that choosing the friction parameter positive resolves d'Alembert's paradox. We then further examine the effect of the friction parameter $β$ on the drag coefficient. In particular, we show that for large $β$ the drag coefficient corresponds well with experimental values. Moreover, we provide numerical evidence that the Newton continuation method (moving from small to large Reynold's numbers) requires fewer iterations to succeed. Thus the slip boundary condition is advantageous also from a computational perspective.
△ Less
Submitted 26 April, 2022;
originally announced April 2022.
-
Extrapolating the Arnoldi Algorithm To Improve Eigenvector Convergence
Authors:
Sara Pollock,
L. Ridgway Scott
Abstract:
We consider extrapolation of the Arnoldi algorithm to accelerate computation of the dominant eigenvalue/eigenvector pair. The basic algorithm uses sequences of Krylov vectors to form a small eigenproblem which is solved exactly. The two dominant eigenvectors output from consecutive Arnoldi steps are then recombined to form an extrapolated iterate, and this accelerated iterate is used to restart th…
▽ More
We consider extrapolation of the Arnoldi algorithm to accelerate computation of the dominant eigenvalue/eigenvector pair. The basic algorithm uses sequences of Krylov vectors to form a small eigenproblem which is solved exactly. The two dominant eigenvectors output from consecutive Arnoldi steps are then recombined to form an extrapolated iterate, and this accelerated iterate is used to restart the next Arnoldi process. We present numerical results testing the algorithm on a variety of cases and find on most examples it substantially improves the performance of restarted Arnoldi. The extrapolation is a simple post-processing step which has minimal computational cost.
△ Less
Submitted 15 March, 2021;
originally announced March 2021.
-
Van der Waals interactions between two hydrogen atoms: The next orders
Authors:
Éric Cancès,
Rafaël Coyaud,
L. Ridgway Scott
Abstract:
We extend a method (E. Cancès and L.R. Scott, SIAM J. Math. Anal., 50, 2018, 381--410) to compute more terms in the asymptotic expansion of the van der Waals attraction between two hydrogen atoms. These terms are obtained by solving a set of modified Slater-Kirkwood partial differential equations. The accuracy of the method is demonstrated by numerical simulations and comparison with other methods…
▽ More
We extend a method (E. Cancès and L.R. Scott, SIAM J. Math. Anal., 50, 2018, 381--410) to compute more terms in the asymptotic expansion of the van der Waals attraction between two hydrogen atoms. These terms are obtained by solving a set of modified Slater-Kirkwood partial differential equations. The accuracy of the method is demonstrated by numerical simulations and comparison with other methods from the literature. It is also shown that the scattering states of the hydrogen atom, that are the states associated with the continuous spectrum of the Hamiltonian, have a major contribution to the C$_6$ coefficient of the van der Waals expansion.
△ Less
Submitted 8 July, 2020;
originally announced July 2020.
-
A Reynolds-robust preconditioner for the Scott-Vogelius discretization of the stationary incompressible Navier-Stokes equations
Authors:
Patrick E. Farrell,
Lawrence Mitchell,
L. Ridgway Scott,
Florian Wechsung
Abstract:
Augmented Lagrangian preconditioners have successfully yielded Reynolds-robust preconditioners for the stationary incompressible Navier-Stokes equations, but only for specific discretizations. The discretizations for which these preconditioners have been designed possess error estimates which depend on the Reynolds number, with the discretization error deteriorating as the Reynolds number is incre…
▽ More
Augmented Lagrangian preconditioners have successfully yielded Reynolds-robust preconditioners for the stationary incompressible Navier-Stokes equations, but only for specific discretizations. The discretizations for which these preconditioners have been designed possess error estimates which depend on the Reynolds number, with the discretization error deteriorating as the Reynolds number is increased. In this paper we present an augmented Lagrangian preconditioner for the Scott-Vogelius discretization on barycentrically-refined meshes. This achieves both Reynolds-robust performance and Reynolds-robust error estimates. A key consideration is the design of a suitable space decomposition that captures the kernel of the grad-div term added to control the Schur complement; the same barycentric refinement that guarantees inf-sup stability also provides a local decomposition of the kernel of the divergence. The robustness of the scheme is confirmed by numerical experiments in two and three dimensions.
△ Less
Submitted 6 July, 2021; v1 submitted 20 April, 2020;
originally announced April 2020.
-
Robust multigrid methods for nearly incompressible elasticity using macro elements
Authors:
Patrick E. Farrell,
Lawrence Mitchell,
L. Ridgway Scott,
Florian Wechsung
Abstract:
We present a mesh-independent and parameter-robust multigrid solver for the Scott-Vogelius discretisation of the nearly incompressible linear elasticity equations on meshes with a macro element structure. The discretisation achieves exact representation of the limiting divergence constraint at moderate polynomial degree. Both the relaxation and multigrid transfer operators exploit the macro struct…
▽ More
We present a mesh-independent and parameter-robust multigrid solver for the Scott-Vogelius discretisation of the nearly incompressible linear elasticity equations on meshes with a macro element structure. The discretisation achieves exact representation of the limiting divergence constraint at moderate polynomial degree. Both the relaxation and multigrid transfer operators exploit the macro structure for robustness and efficiency. For the relaxation, we use the existence of local Fortin operators on each macro cell to construct a local space decomposition with parameter-robust convergence. For the transfer, we construct a robust prolongation operator by performing small local solves over each coarse macro cell. The necessity of both components of the algorithm is confirmed by numerical experiments.
△ Less
Submitted 1 November, 2021; v1 submitted 5 February, 2020;
originally announced February 2020.
-
An Adaptive Savitsky-Golay Filter for Smoothing Finite Element Computation
Authors:
Teodoro Collin,
Gordon Kindlmann,
L. Ridgway Scott
Abstract:
The smoothing technique of Savitzky and Golay is extended to data defined on multidimensional meshes. A smoothness-increasing accuracy-conserving (SIAC) filter is defined that is suitable for use with finite-element computation.
The smoothing technique of Savitzky and Golay is extended to data defined on multidimensional meshes. A smoothness-increasing accuracy-conserving (SIAC) filter is defined that is suitable for use with finite-element computation.
△ Less
Submitted 28 December, 2020; v1 submitted 2 November, 2019;
originally announced November 2019.
-
Topological Optimization of the Evaluation of Finite Element Matrices
Authors:
Robert C. Kirby,
Anders Logg,
L. Ridgway Scott,
Andy R. Terrel
Abstract:
We present a topological framework for finding low-flop algorithms for evaluating element stiffness matrices associated with multilinear forms for finite element methods posed over straight-sided affine domains. This framework relies on phrasing the computation on each element as the contraction of each collection of reference element tensors with an element-specific geometric tensor. We then pres…
▽ More
We present a topological framework for finding low-flop algorithms for evaluating element stiffness matrices associated with multilinear forms for finite element methods posed over straight-sided affine domains. This framework relies on phrasing the computation on each element as the contraction of each collection of reference element tensors with an element-specific geometric tensor. We then present a new concept of complexity-reducing relations that serve as distance relations between these reference element tensors. This notion sets up a graph-theoretic context in which we may find an optimized algorithm by computing a minimum spanning tree. We present experimental results for some common multilinear forms showing significant reductions in operation count and also discuss some efficient algorithms for building the graph we use for the optimization.
△ Less
Submitted 14 May, 2012;
originally announced May 2012.
-
Optimizing the Evaluation of Finite Element Matrices
Authors:
Robert C. Kirby,
Matthew Knepley,
Anders Logg,
L. Ridgway Scott
Abstract:
Assembling stiffness matrices represents a significant cost in many finite element computations. We address the question of optimizing the evaluation of these matrices. By finding redundant computations, we are able to significantly reduce the cost of building local stiffness matrices for the Laplace operator and for the trilinear form for Navier-Stokes. For the Laplace operator in two space dimen…
▽ More
Assembling stiffness matrices represents a significant cost in many finite element computations. We address the question of optimizing the evaluation of these matrices. By finding redundant computations, we are able to significantly reduce the cost of building local stiffness matrices for the Laplace operator and for the trilinear form for Navier-Stokes. For the Laplace operator in two space dimensions, we have developed a heuristic graph algorithm that searches for such redundancies and generates code for computing the local stiffness matrices. Up to cubics, we are able to build the stiffness matrix on any triangle in less than one multiply-add pair per entry. Up to sixth degree, we can do it in less than about two. Preliminary low-degree results for Poisson and Navier-Stokes operators in three dimensions are also promising.
△ Less
Submitted 14 May, 2012;
originally announced May 2012.
-
Unstructured Geometric Multigrid in Two and Three Dimensions on Complex and Graded Meshes
Authors:
Peter R. Brune,
Matthew G. Knepley,
L. Ridgway Scott
Abstract:
The use of multigrid and related preconditioners with the finite element method is often limited by the difficulty of applying the algorithm effectively to a problem, especially when the domain has a complex shape or adaptive refinement. We introduce a simplification of a general topologically-motivated mesh coarsening algorithm for use in creating hierarchies of meshes for geometric unstructured…
▽ More
The use of multigrid and related preconditioners with the finite element method is often limited by the difficulty of applying the algorithm effectively to a problem, especially when the domain has a complex shape or adaptive refinement. We introduce a simplification of a general topologically-motivated mesh coarsening algorithm for use in creating hierarchies of meshes for geometric unstructured multigrid methods. The connections between the guarantees of this technique and the quality criteria necessary for multigrid methods for non-quasi-uniform problems are noted. The implementation details, in particular those related to coarsening, remeshing, and interpolation, are discussed. Computational tests on pathological test cases from adaptive finite element methods show the performance of the technique.
△ Less
Submitted 5 April, 2011; v1 submitted 1 April, 2011;
originally announced April 2011.
