Abstract
We analyze the stability and convergence of first-order accurate and second-order accurate timestepping schemes for the Navier--Stokes equations with variable viscosity. These schemes are characterized by a mixed implicit/explicit treatment of the viscous term, in which a numerical parameter, $\lambda$, determines the degree of splitting between the implicit and explicit contributions. The reason for this splitting is that it avoids the need to solve computationally expensive linear systems that may change at each timestep. Provided the parameter $\lambda$ is within a permissible range, we prove that the first-order accurate and second-order accurate schemes are convergent. We show further that the efficiency of the second-order accurate scheme depends on how $\lambda$ is chosen within the permissible range, and we discuss choices that work well in practice. We use parameters motivated by this analysis to simulate internal gravity waves, which arise in stratified fluids with variable density. We examine h...
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call