Handling solid-fluid interfaces for viscous flows: Explicit jump approximation vs. ghost cell approaches

NCUIR > College of Earth Sciences > Graduate Institute of Hydrological and Oceanic Sciences > journal & Dissertation > Item 987654321/50482

Please use this identifier to cite or link to this item: http://ir.lib.ncu.edu.tw/handle/987654321/50482

Title:	Handling solid-fluid interfaces for viscous flows: Explicit jump approximation vs. ghost cell approaches
Authors:	Zhang,QH;Liu,PLF
Contributors:	水文與海洋科學研究所
Keywords:	EMBEDDED BOUNDARY METHOD;TIME-DEPENDENT BOUNDARY;NAVIER-STOKES EQUATIONS;ELLIPTIC-EQUATIONS;DISCONTINUOUS COEFFICIENTS;IRREGULAR DOMAINS;SINGULAR SOURCES;DIFFERENTIAL-EQUATIONS;POISSONS-EQUATION;HEAT-EQUATION
Date:	2010
Issue Date:	2012-03-27 17:33:08 (UTC+8)
Publisher:	國立中央大學
Abstract:	The ghost cell approaches (GCA) for handling stationary solid boundaries, regular or irregular, are first investigated theoretically and numerically for the diffusion equation with Dirichlet boundary conditions. The main conclusion of this part of investigation is that the approximation for the diffusion term has to be second-order accurate everywhere in order for the numerical solution to be rigorously second-order accurate. Violating this principle, the linear and quadratic GCAs have the following shortcomings: (1) restrictive constraints on grid size when the viscosity is small; (2) susceptibleness to instability of a time-explicit formulation for strongly transient flows; (3) convergence deterioration to zeroth- or first-order for solutions with high-frequency modes. Therefore, the widely-used linear extrapolation for enforcing no-slip boundary conditions should be avoided, even for regular solid boundaries. As a remedy, a simple method based on explicit jump approximation (EJA) is proposed. EJA hinges on the idea that a velocity-derivative jump at the boundary reduces to the value of the velocity-derivative at the fluid side because the velocity of the stationary boundary is zero. Although the time-marching linear system of EJA is not symmetric, it is strictly diagonal dominant with positive diagonal entries. Numerical results show that, over a large range of viscosity and grid sizes, EJA performs much better than GCAs in terms of stability and accuracy. Furthermore, the second-order convergence of EJA does not depend on viscosity and the spectrum of the solution, as those of GCAs do. This paper is written with enough details so that one can reproduce the numerical results. (C) 2010 Elsevier Inc. All rights reserved.
Relation:	JOURNAL OF COMPUTATIONAL PHYSICS
Appears in Collections:	[Graduate Institute of Hydrological and Oceanic Sciences] journal & Dissertation

Files in This Item:

File	Description	Size	Format
index.html		0Kb	HTML	545	View/Open

社群 sharing

Loading...