Source author record

L. Ridgway Scott

L. Ridgway Scott appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
14topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Resolution of d'Alembert's paradox using slip boundary conditions: The effect of the friction parameter on the drag coefficient

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.

preprint2020arXiv

Van der Waals interactions between two hydrogen atoms: The next orders

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.

preprint2016arXiv

A Nonlocal Poisson-Fermi Model for Ionic Solvent

We propose a nonlocal Poisson-Fermi model for ionic solvent that includes ion size effects and polarization correlations among water molecules in the calculation of electrostatic potential. It includes the previous Poisson-Fermi models as special cases, and its solution is the convolution of a solution of the corresponding nonlocal Poisson dielectric model with a Yukawa-type kernel function. Moreover, the Fermi distribution is shown to be a set of optimal ionic concentration functions in the sense of minimizing an electrostatic potential free energy. Finally, numerical results are reported to show the difference between a Poisson-Fermi solution and a corresponding Poisson solution.

preprint2016arXiv

A Stochastic Performance Model for Pipelined Krylov Methods

Pipelined Krylov methods seek to ameliorate the latency due to inner products necessary for projection by overlapping it with the computation associated with sparse matrix-vector multiplication. We clarify a folk theorem that this can only result in a speedup of $2\times$ over the naive implementation. Examining many repeated runs, we show that stochastic noise also contributes to the latency, and we model this using an analytical probability distribution. Our analysis shows that speedups greater than $2\times$ are possible with these algorithms.

preprint2012arXiv

Optimizing the Evaluation of Finite Element Matrices

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.

preprint2012arXiv

Topological Optimization of the Evaluation of Finite Element Matrices

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.

preprint2011arXiv

Unstructured Geometric Multigrid in Two and Three Dimensions on Complex and Graded Meshes

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.