Release oscillation in a hollow fiber – Part 1: Mathematical model and fast estimation of its frequency

Introduction

Nano scale or micro scale flows have been caught much attention due to the fast development of micro-electro-mechanical systems, ¹ micro scale devices, e.g., silicon rectangular micro-orifice, ² micropump,^3–5 microfluidics device, ⁶ three-dimensional integrated circuits,^7,8 micro scale heat exchanger, ⁹ the ion release from hollow fibers,^10,11 the drug release from nanoporous fibers,^12–15 and nanofiber membranes.^16–23 Due to capillary action, water will be sucked into the gap of a micro scale device, which will greatly affect the device’s performance. On the other hand, ions in a hollow fiber and drug in a nanofiber membrane require continuous release. But capillary action forbids a complete release, and the fluid in the tube will be oscillated due to some perturbations. The nano/micro scale fluid oscillation will greatly affect the property of such systems, especially the micro-channel cooling property, heat transfer, and drug release property.

In this paper, we will study the fluid oscillation in a nano/micro capillary tube of a hollow fiber.

Fluid oscillation in capillary tube of a hollow fiber

We consider a capillary tube as illustrated in Figure 1. The tube is immersed in water. The enclosed fluid (blue part in Figure 1) can be considered as a coolant flow, a released solution, and others.

OPEN IN VIEWER

Figure 1.

Nano/micro flow in a capillary tube, where x 0 is the equilibrium point and the mass center of the fluid in the tube.

According to Newton’s second law, the free oscillation of the fluid in nano/micro tube can be written as

m x ¨ + F = 0

(1)

The fluid mass in the tube can be calculated as

m = 2 ( x 0 + x ) S ρ

(2)

The force F can be written in the form

F = S ( P in − P out )

(3)

The isolated air in the tube follows the Boyle's law

P in V = P 0 i n V 0

(4)

Equation (4) is valid for the case when the temperature and the gas involved in the hollow fiber remain unchanged during the release process.

The air volume can be calculated as

V = S ( L − 2 x 0 − 2 x )

(5)

The air pressure P in can be expressed as

P in = P 0 i n V 0 S ( L − 2 x 0 − 2 x )

(6)

The governing equation, equation (1), becomes

2 ( x 0 + x ) S ρ x ¨ + P 0 i n V 0 L − 2 x 0 − 2 x − S P out = 0

(7)

( x 0 + x ) x ¨ + a + bx L 0 − x = 0

(8)

The initial conditions are

x ( 0 ) = A , x ˙ ( 0 ) = 0

(9)

Frequency–amplitude relation of the fluid oscillator

In this paper, we will suggest two simple but effective approaches to a fast estimation of the solution of the fluid oscillator. One is the Taylor series-based approach ²⁴ and the other is the max–min approach. ²⁵ Other simple methods are available in Tian and Liu, ²⁶ He, ²⁷ He et al., ²⁸ and Ren and Hu. ²⁹

Taylor series-based solution

Following equation (8) and the initial conditions given in equation (9), we have ²⁴

x ¨ ( 0 ) = − a + bA ( L 0 − A ) ( x 0 + A )

(10)

The Taylor series to second order is

x ( t ) = x ( 0 ) + x ˙ ( 0 ) t + 1 2 x ¨ ( 0 ) t 2 = A − a + bA 2 ( L 0 − A ) ( x 0 + A ) t 2

(11)

Considering the oscillation property, we have ²⁴

x ( T / 4 ) = A − a + bA 2 ( L 0 − A ) ( x 0 + A ) ( T 4 ) 2 = 0

(12)

where T is the period of the fluid oscillator.

Solving the period from equation (12), we have

T = 4 2 A ( L 0 − A ) ( x 0 + A ) a + bA

(13)

Its frequency reads

ω = 2 π T = π 2 a + bA 2 A ( L 0 − A ) ( x 0 + A )

(14)

or

ω = 2 π T = π 2 P 0 i n V 0 − S P out ( L − 2 x 0 ) + 2 S P out A 4 S ρ A ( L − 2 x 0 − 2 A ) ( x 0 + A )

(15)

In equation (15) accuracy is not enough, but it can depict the main factors affecting the frequency.

Max–min approach

We write the fluid oscillator in the form ²⁵

x ¨ + f ( x ) x = 0

(16)

where

f ( x ) = a + bx ( L 0 − x ) ( x 0 + x ) x

(17)

The max–min approach is to calculate the square of the frequency in a simple form ²⁵

ω 2 = f ( x ( t 0 ) ) = f ( kA )

(18)

where t 0 is a location point, k is a constant 0 < k < 1.

Choosing location point ω t = π / 3 , we obtain the square of frequency

ω 2 = a + 1 2 b A 1 2 ( L 0 − 1 2 A ) ( x 0 + 1 2 A ) A

(19)

Discussion and conclusion

The fluid oscillation will greatly affect drug release from hollow fibers or nanofiber member and micro-channel cooling property in micro devices. In practical applications, we need to improve cooling efficiency in a three-dimensional integrated circuit, to enhance heat transfer in a micro scale heat exchanger and to adjust the drug release ratio from a hollow fiber. So the exact frequency is not necessary, and the approximate one is enough for this purpose. For example, if we want to improve the drug release ratio from a hollow fiber, we should increase the frequency of the fluid oscillation, according to equation (15) or (19), the value of P out should decrease. To this end, the hollow fibers should not be placed deeply into water. Liu et al. suggested a fractional model for the ion release from the hollow fiber, ³⁰ and this theory can be easily extended to a fractional model^31–38 for ion release from a hollow fiber.

This paper elucidates the main parameters affecting the frequency in the ion release in a hollow fiber, and the theoretical results can be used to control drug release and heat transfer and other nano/micro tube flows.