Source author record

Shuxian Xu

Shuxian Xu 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

Catalog footprint

What is connected

3works

2topics

3close 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

A Doubly Adaptive Penalty Method for the Navier Stokes Equations

We develop, analyze and test adaptive penalty parameter methods. We prove unconditional stability for velocity when adapting the penalty parameter, $ε,$ and stability of the velocity time derivative under a condition on the change of the penalty parameter, $ε(t_{n+1})-ε(t_n)$. The analysis and tests show that adapting $ε(t_{n+1})$ in response to $\nabla\cdot u(t_n)$ removes the problem of picking $ε$ and yields good approximations for the velocity. We provide error analysis and numerical tests to support these results. We supplement the adaptive-$ε$ method by also adapting the time-step. The penalty parameter $ε$ and time-step are adapted independently. We further compare first, second and variable order time-step algorithms. Accurate recovery of pressure remains an open problem.

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.

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint