Source author record

William Layton

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

ResearcherUnclaimed source record

math.NA Numerical Analysis physics.flu-dyn

Catalog footprint

What is connected

5works

3topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Conditioning of linear systems arising from penalty methods

Penalizing incompressibility in the Stokes problem leads, under mild assumptions, to matrices with condition numbers $κ=\mathcal{O} (\varepsilon ^{-1}h^{-2})$, $\varepsilon =$ penalty parameter $<<1$, and $ h= $ mesh width $<1$. Although $κ=\mathcal{O}(\varepsilon ^{-1}h^{-2}) $ is large, practical tests seldom report difficulty in solving these systems. In the SPD case, using the conjugate gradient method, this is usually explained by spectral gaps occurring in the penalized coefficient matrix. Herein we point out a second contributing factor. Since the solution is approximately incompressible, solution components in the eigenspaces associated with the penalty terms can be small. As a result, the effective condition number can be much smaller than the standard condition number.

preprint2022arXiv

Stability in 3d of a sparse grad-div approximation of the Navier-Stokes equations

Inclusion of a term $-γ\nabla\nabla\cdot u$, forcing $\nabla\cdot u$ to be pointwise small, is an effective tool for improving mass conservation in discretizations of incompressible flows. However, the added grad-div term couples all velocity components, decreases sparsity and increases the condition number in the linear systems that must be solved every time step. To address these three issues various sparse grad-div regularizations and a modular grad-div method have been developed. We develop and analyze herein a synthesis of a fully decoupled, parallel sparse grad-div method of Guermond and Minev with the modular grad-div method. Let $G^{\ast}=-diag(\partial_{x}^{2},\partial_{y}^{2},\partial_{z}^{2})$ denote the diagonal of $G=-\nabla\nabla\cdot$, and $α\geq0$ an adjustable parameter. The 2-step method considered is $$\begin{eqnarray} 1 &:&\frac{\widetilde{u}^{n+1}-u^{n}}{k}+u^{n}\cdot \nabla \widetilde{u}^{n+1}+\nabla p^{n+1}-νΔ\widetilde{u}^{n+1}=f\text{ & }\nabla \cdot \widetilde{u}^{n+1}=0,\\ 2 &:&\left[ \frac{1}{k}I+(γ+α)G^{\ast }\right] u^{n+1}=\frac{1}{k }\widetilde{u}^{n+1}+\left[ (γ+α)G^{\ast }-γG\right] u^{n}. \end{eqnarray}$$ We prove its unconditional, nonlinear, long time stability in $3d$ for $α\geq0.5γ$. The analysis also establishes that the method controls the persistent size of $||\nabla\cdot u||$ in general and controls the transients in $||\nabla\cdot u||$ when $u(x,0)=0$ and $f(x,t)\neq0$ provided $α>0.5γ$. Consistent numerical tests are presented.

preprint2020arXiv

Analysis of the variable step method of Dahlquist, Liniger and Nevanlinna for fluid flow

The two-step time discretization proposed by Dahlquist, Liniger and Nevanlinna is variable step $G$-stable. (In contrast, for increasing time steps, the BDF2 method loses $A$-stability and suffers non-physical energy growth in the approximate solution.) While unexplored, it is thus ideal for time accurate approximation of the Navier-Stokes equations. This report presents an analysis, for variable time-steps, of the method's stability and convergence rates when applied to the NSE. It is proven that the method is variable step, unconditionally, long time stable and second order accurate. Variable step error estimates are also proven. The results are supported by several numerical tests.

preprint2016arXiv

Energy Dissipation in the Smagorinsky Model of Turbulence

The Smagorinsky model, unmodified, is often reported to severely overdiffuse flows. Previous estimates of the energy dissipation rate of the Smagorinsky model for shear flows reflect a blow up of model energy dissipation as Re increases. This blow up is consistent with the numerical evidence and leads to the question: Is the over dissipation due to the influence of the turbulent viscosity in boundary layers alone or is its action on small scales generated by the nonlinearity through the cascade also a contributor? This report develops model dissipation estimates for body force driven flow under periodic boundary conditions (and thus only with nonlinearity generated small scales). It is proven that the model's time averaged energy dissipation rate satisfies the same upper bound as for the NSE plus one additional term that vanishes uniformly in the Reynolds number as the Smagorinsky length scale decreases. Since this estimate is consistent with that observed for the NSE, it establishes that, without boundary layers, the Smagorinsky model does not over dissipate.

preprint2015arXiv

Algorithms and Models for Turbulence Not at Statistical Equilibrium

Standard eddy viscosity models, while robust, cannot represent backscatter and have severe difficulties with complex turbulence not at statistical equilibrium. This report gives a new derivation of eddy viscosity models from an equation for the evolution of variance in a turbulent flow. The new derivation also shows how to correct eddy viscosity models. The report proves the corrected models preserve important features of the true Reynolds stresses. It gives algorithms for their discretization including a minimally invasive modular step to adapt an eddy viscosity code to the extended models. A numerical test is given with the usual and over diffusive Smagorinsky model. The correction (scaled by $10^{-8}$ ) does successfully exhibit intermittent backscatter.