Researcher profile

Shuxian Xu

Shuxian Xu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

math.NA Numerical Analysis

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author

3works

0followers

2topics

3close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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