Source author record

Tsorng-Whay Pan

Tsorng-Whay Pan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

preprint2016arXiv

A numerical study of the transition to oscillatory flow in 3D lid-driven cubic cavity flows

In this article, three dimensional (3D) lid-driven cubic cavity flows have been studied numerically for various values of Reynolds number ($Re$). The numerical solution of the Navier-Stokes equations modeling incompressible viscous fluid flow in a cubic cavity is obtained via a methodology combining a first order accurate operator-splitting, $L^2$-projection Stokes solver, a wave-like equation treatment of the advection and finite element methods. The numerical results obtained for Re$=$400, 1000, and 3200 show a good agreement with available numerical and experimental results in literature. Simulation results predict that the critical Re$_{cr}$ for the transition from steady flow to oscillatory (a Hopf bifurcation) is somewhere in [1870, 1875] for the mesh size $h=1/96$. Via studying the flow field distortion of fluid flow at Re before and after Re$_{cr}$, the occurrence of the first pair of Taylor-Görtler-like vortices is connected to the flow field distortion at the transition from steady flow to oscillatory flow in 3D lid-driven cubic cavity flows for Re $< 2000$.

preprint2016arXiv

Dynamics of particle sedimentation in viscoelastic fluids: A numerical study on particle chain in two-dimensional narrow channel

In this article we present a numerical method for simulating the sedimentation of circular particles in two-dimensional channel filled with a viscoelastic fluid of FENE-CR type, which is generalized from a domain/distributed Lagrange multiplier method with a factorization approach for Oldroyd-B fluids developed in [J. Non-Newtonian Fluid Mech. 156 (2009) 95]. Numerical results suggest that the polymer extension limit L for the FENE-CR fluid has no effect on the final formation of vertical chain for the cases of two disks and three disks in two-dimensional narrow channel, at least for the values of L considered in this article; but the intermediate dynamics of particle interaction before having a vertical chain can be different for the smaller values of L when increasing the relaxation time. For the cases of six particles sedimenting in FENE-CR type viscoelastic fluid, the formation of chain of 4 to 6 disks does depend on the polymer extension limit L. For the smaller values of L, FENE-CR type viscoelastic fluid can not bring them together like the case of these particles settling in a vertical chain formation in Oldroyd-B fluid; but two separated chains of three disks are formed. Similar results for the case of ten disks are also obtained. The numerical results of several more particle cases suggest that for smaller values of L, the length of the vertical chain is shorter and the size of cluster is smaller.

preprint2014arXiv

The dynamics of inextensible capsules in shear flow under the effect of the nature state

The effect of the nature state on the motion of an inextensible capsule in simple shear flow has been studied in this paper. Besides the viscosity ratio of the internal fluid and external fluid of the capsule, the nature state effect also plays a role for having the transition between two well known motions, tumbling and tank-treading (TT) with the long axis oscillating about a fixed inclination angle (a swinging mode), when varying the shear rate. The intermittent region between tumbling and TT with a swinging mode of the capsule with a biconcave rest shape has been obtained in a narrow range of the capillary number. In such region, the dynamics of the capsule is a mixture of tumbling and TT with a swinging mode; when having the tumbling motion, the membrane tank-tread backward and forward within a small range. As the capillary number is very close to and below the threshold for the pure TT with a swinging mode, the capsule tumbles once after several TT periods in each cycle. The number of TT periods in one cycle decreases with respect to the decreasing of the capillary number, until the capsule has one tumble and one TT period alternatively and such alternating motion exists over a range of the capillary number; and then the capsule performs more tumbling between two consecutive TT periods when reducing the capillary number further, and finally shows pure tumbling. The critical value of the swelling ratio for having the intermittent region has been estimated.

preprint2014arXiv

The motion of a neutrally buoyant particle of an elliptic shape in two dimensional shear flow: a numerical study

In this paper, we investigate the motion of a neutrally buoyant cylinder of an elliptic shape freely moving in two dimensional shear flow by direct numerical simulation. An elliptic shape cylinder in shear flow, when initially being placed at the middle between two walls, either keeps rotating or has a stationary inclination angle depending on the particle Reynolds number $Re=G_r r_a^2/ν$, where $G_r$ is the shear rate, $r_a$ is the semi-long axis of the elliptic cylinder and $ν$ is the kinetic viscosity of the fluid. The critical particle Reynolds number $Re_{cr}$ for the transition from a rotating motion to a stationary orientation depends on the aspect ratio $AR=r_b/r_a$ and the confined ratio $K=2r_a/H$ where $r_b$ is the semi-short axis of the elliptic cylinder and $H$ is the distance between two walls. Although the increasing of either parameters makes an increase in $Re_{cr}$, the dynamic mechanism is distinct. The $AR$ variation causes the change of geometry shape; however, the $K$ variation influences the wall effect. The stationary inclination angle of non-rotating slender elliptic cylinder with smaller confined ratio seems to depend only on the value of $Re-Re_{cr}$. An expected equilibrium position of the cylinder mass center in shear flow is the centerline between two walls, but when placing the particle away from the centerline initially, it migrates either toward an equilibrium height away from the middle between two walls or back to the middle depending on the confined ratio and particle Reynolds number.

preprint2013arXiv

An oscillating motion of a red blood cell and a neutrally buoyant particle in Poiseuille flow in a narrow channel

Two motions of oscillation and vacillating breathing (swing) of a red blood cell have been observed in bounded Poiseuille flows (Phys. Rev. E 85, 16307 (2012)). To understand such motions, we have studied the oscillating motion of a neutrally buoyant rigid particle of the same shape in Poiseuille flow in a narrow channel and obtained that the crucial point is to have the particle interacting with Poiseuille flow with its mass center moving up and down in the channel central region. Since the mass center of the cell migrates toward the channel central region, its oscillating motion of the inclination angle is similar to the aforementioned motion as long as the cell keeps the shape of long body. But as the up-and-down oscillation of the cell mass center damps out, the oscillating motion of the inclination angle also damps out and the cell inclination angle approaches to a fixed angle.

preprint2013arXiv

Circular band formation for incompressible viscous fluid--rigid particle mixtures in a rotating cylinder

In this paper we have investigated a circular band formation of fluid-rigid particle mixtures in a fully filled cylinder horizontally rotating about its cylinder axis by direct numerical simulation. These phenomena are modeled by the Navier-Stokes equations coupled to the Euler-Newton equations describing the rigid solid motion of the non-neutrally particles. The formation of circular bands studied in this paper is not resulted by mutual interaction between the particles and the periodic inertial waves in the cylinder axis direction (as suggested in Phys. Rev. E, 72, 021407 (2005)), but due to the interaction of particles. When a circular band is forming, the part of the band formed by the particles moving downward becomes more compact due to the particle interaction strengthened by the downward acceleration from the gravity. The part of a band formed by the particles moving upward is always loosening up due to the slow down of the particle motion by the counter effect of the gravity. To form a compact circular band (not a loosely one), enough particles are needed to interact among themselves continuously through the entire circular band at a rotating rate so that the upward diffusion of particles can be balanced by the compactness process when these particles moving downward.

preprint2013arXiv

Numerical simulation of red blood cell suspensions behind a moving interface in a capillary

Computational modeling and simulation are presented on the motion of red blood cells behind a moving interface in a capillary. The methodology is based on an immersed boundary method and the skeleton structure of the red blood cell (RBC) membrane is modeled as a spring network. The computational domain is moving with either a designated RBC or an interface in an infinitely long two-dimensional channel with an undisturbed flow field in front of the domain. The tanking-treading and the inclination angle of a cell in a simple shear flow are briefly discussed for the validation purpose. We then present the results of the motion of red blood cells behind a moving interface in a capillary, which show that the RBCs with higher velocity than the interface speed form a concentrated slug behind the interface.

preprint2012arXiv

A numerical study of the motion of a neutrally buoyant cylinder in two dimensional shear flow

In this paper, we investigate the motion of a neutrally buoyant cylinder of circular or elliptic shape in two dimensional shear flow of a Newtonian fluid by direct numerical simulation. The numerical results are validated by comparisons with existing theoretical, experimental and numerical results, including a power law of the normalized angular speed versus the particle Reynolds number. The centerline between two walls is an expected equilibrium position of the cylinder mass center in shear flow. When placing the particle away from the centerline initially, it migrates toward another equilibrium position for higher Reynolds numbers due to the interplay between the slip velocity, the Magnus force, and the wall repulsion force.