# Design and Simulation of a Fused Silica Space Cell Culture and Observation Cavity with Microfluidic and Temperature Controlling

Authors: Shangchun Fan; Jinhao Sun; Weiwei Xing; Cheng Li; Dongxue Wang
Journal: Journal of Applied Mathematics
 (2013)

---

## Abstract

We report a principle prototype of space animal cell perfusion culture and observation. Unlike previous work, our cell culture system cannot only realize microfluidic and temperature controlling, automatic observation, and recording but also meet an increasing cell culture at large scale operation and overcome shear force for animal cells. A key component in the system is ingenious structural fused silica cell culture cavity with the wedge-shaped connection. Finite volume method (FVM) is applied to calculate its multipoint flow field, pressure field, axial velocity, tangential velocity, and radial velocity. In order to provide appropriate flow rate, temperature, and shear force for space animal cell culture, a closed-loop microfluidic circuit and proportional, integrating, and differentiation (PID) algorithm are employed. This paper also illustrates system architecture and operating method of the principle prototype. The dynamic culture, autofocus observation, and recording of M763 cells are performed successfully within 72 h in the laboratory environment. This research can provide a reference for space flight mission that carries an apparatus with similar functions.


---

## Body

## 1. Introduction

The increasing demand of recombinant therapeutic proteins, monoclonal antibodies, and vaccines produced by mammalian cells has stimulated the development of space cell culture technology, for which almost all space powers have a huge interest. In order to make an in-depth research of space cell culture characteristics, some space powers have developed various space cell culture devices [1–4]. Many of these have completed space flight mission that are shown in Table 1 and Figure 1. Such space cell culture apparatus not only provides appropriate temperature, adequate nutrition and oxygen, low concentration of toxic secretions, and reasonable mechanical property but also records details of cell proliferation [5]. However, space perfusion cultures require more efficient and advanced system, the most frequently reported based on microfluidic chip with kinds of sensors [6–10] and aerospace heat transfer [11]. Rotating wall bioreactor is widely used because the internal spin-filters are conducive to perfusion rates and reproductive rate at large-scale operation [12, 13]. In this case, the defects such as high energy consumption, noise, low space utilization, and shear force have been considered as limitation of culture duration and process performance.Table 1
Space cell culture instruments development at home and abroad.
Name
Country
Flight plan
Time
Woodlawn wanderer-9
US
Salyut-3
1978
Biotherm
Switzerland
Salyut-4; Salyut-5
1982
Cytos I
France
Salyut-6
1982
Biorack
Germany
D-1
1983
Carry-on
Switzerland
STS-8
1984
Cytos II
France
Salyut-7
1986
Dynamic cell culture system
Switzerland
IML-1
1988
Static training device
China
Jian-yi-yi satellite
1989
Cell culture device
Japan
IML-2
1994
Improved CCS
China
Jian-yi-yi satellite
1994
Miniature space cell culture device (SBRI)
Switzerland
STS-65; S/MM-03
1994
Cell culture device
China
Shen Zhou 3
2002Successful applications of space cell culture instruments.
(a)
ESA/CSA cooperation of unmanned Foton M3 task tray of cell culture
(b)
NASA BTS cells culture module, for the space station
(c)
NASA BioServe culture device
(d)
NASA-CGBADespite the great importance of cell culture data after space flight, many studies have been carried out mainly to recognize the differences between original ground data before space flight experiment and new character after it. These studies apparently cannot meet the current research requirements today [14, 15]. Early efforts indicate that microfluidic technology plays an important role in large-scale cell operation, real-time observation, and recording. Dynamical culture in space biological experiments becomes increasingly needed [16, 17]. Dynamical culture can ensure that cells always get enough fresh nutrient solution and metabolic waste is moved out during culture process and construct a more suitable environment for long-term cell culture. Furthermore, compared with the conventional cell culture techniques [18], the combination of microfluidic chip and cell culture techniques will benefit long-term cell culture and improve cell biology research levels. However, only a certain aspect but not the complete system is introduced in the majority of reports in this field, not to mention considering the impact of the acceleration of gravity.The study presented in this paper involves fused silica cell culture cavity’s simulation, design, and perfusion experiment, which is carried out under realistic operating conditions of temperature and flow rate, using M763 cells. With a wedge-shaped connection, the newly designed cavity is free from the effects of gravitational acceleration before entering orbit and is better sealed. Fused silica is adopted for its outstanding optical property, insulating property, thermal conductivity, and stability, all of which enable the cavity to be isolated from harsh outside environment. In this paper, firstly the creatively designed cavity structure, simulation, and experimental verification of microflow FVM are presented. Then a description of hardware circuits and software design is given. At last, a complete space cell system is conducted and characteristics are listed. The research has some theoretical and engineering value on the design and implementation of space cell culture system.


## 2. Materials and Methods

Our automatic guided space cell culture system mainly includes cell culture unit, observation and record unit, and intelligent control unit; still all outer space conditions need to be overcome. In the outer space, the equipment is expected to executive automatically operations presetted. The experimental apparatus for cell perfusion simulations is shown in Figure2.Figure 2
Experimental set-up for cell perfusion simulations.


### 2.1. Design and Fabrication of the Cavity

Entire cell culture cavity adopts fused silica, which has the following advantages: perfect thermal stability, low deformation, resistance to high acid and contamination, high hardness, and excellent light transmittance.Cell culture cavity (Figure3) is mainly made up of two parts: the upper cover and the lower cover. The upper and the lower covers are pressed together to form a hollow cell culture cavity. Cell suspension and reagent and metabolic waste go in and out through the three wedge-shaped cylindrical pipes welded to the lower cover. Both the upper and the lower covers have 50 mm diameter circular photic zone; the rest part is made as matte at roughness 3.2 and fastened via eight through-bolts in uniform distribution. The upper cover, lower cover, micro pressure sensor cavity, and three wedge-shaped catheters are shown in the assembly drawing of cell culture cavity (Figure 4).Figure 3
A photograph of cell culture cavity.Figure 4
The assembly drawing of cell culture cavity.To make the structure better sealed, the three wedge-shaped cylindrical pipes in the lower cover and three independent catheters are connected to micro pressure sensor cavity through the wedge connection, which contributes to overcoming the overload when entering orbit. Using Polytetrafluoroethylene (PTFE) tube connects inlet and outset into a rotary value, and then the rotary valve is connected to micropump, thus making it easy for the system to control micro fluid. The pipe in between serves as a waste export as well as a biochemical detection port. The overall structural design is ingenious and highly modularized.The system allows automatically real-time observation and recording in dynamical culture process through circular photic zone, which facilitates light irradiation and benefits the microscopic observation. The practical parameters of cell culture cavity are shown in Figure5. Heating coil (Figure 6) is designed in lower cover surface to provide a suitable temperature.Figure 5
The engineering drawing for upper cover of cell culture cavity.Figure 6
The engineering drawing for lower cover of cell culture cavity.What matters most in this experiment is how to design a perfect cell culture cavity and locate it precisely. A three-dimensional servo platform with high precision is designed to cooperate with microscopic observation. The top of the servo platform is the fixed device of cell culture cavity (Figure7).Figure 7
A photograph of cell culture cavity three-dimensional servo platform.


### 2.2. The Theory of the Computational Fluid Dynamics

Continuous equation, namely, the conversation of mass equation, means that any liquidity issue must agree with the law of conversation of mass. In accordance with the law of conversation of mass, the sum of the fluid rest mass out of control interval in unit time should be equal to the reducing quality caused by density variations in control interval in the same unit time. Based on that, the differential form of fluid flow continuous equation can be derived as(1)∂ρ∂t+∂(ρux)∂x+∂(ρuy)∂y+∂(ρuz)∂z=0,
where ux, uy, and uz are, respectively, the component of velocity (m/s) for three directions, t is time(s), and ρ is density (kg/m3).The Hamiltonian differential operator is derived(2)∇=i∂∂x+j∂∂y+k∂∂z.Equation (1) can be expressed as (3) and divergence form (4):
(3)∂ρ∂t+∇·(ρu→)=0,(4)∂ρ∂t+div(ρu→)=0.The form of the continuity equation is expressed as (5) in a cylindrical coordinate system:
(5)∂ρ∂t+ρurr+∂(ρur)∂r+∂(ρuθ)r∂θ+∂(ρuz)∂z=0.Equation (5) is the normal continuity equation for general use. If for steady flow ∂ρ/∂t=0, its form should be changed to
(6)∂(ρux)∂x+∂(ρuy)∂y+∂(ρuz)∂z=0.For incompressible flow, in the flow process, the densityρ of each dot is always unchanged, that is, ∂ρ/∂t=0 and ρ≠0; the form of the continuity equation of incompressible fluid can be derived as
(7)∂ux∂x+∂uy∂y+∂uz∂z=0.Cylindrical coordinate’s form of (7) is
(8)urr+∂ur∂r+∂uθr∂θ+∂uz∂z=0.The motion equation of viscous fluid was first proposed in 1827 by Navier who only took the flow of incompressible fluid into account. Poisson proposed the equation of the compressible fluid in 1831. Saint-Venant and Stokes, respectively, in 1843 and 1845, independently proposed the form in which viscosity coefficient was a constant, now known as Navier-Stokes equation, referred to as the N-S equation.For incompressible viscous fluid, the N-S equation is(9)∂u→∂t+(u→·∇)u→=f-1ρ∇p+μρ∇2u→.In (9), u→ (vector symbol) represents velocity (m/s), p represents pressure (Pa), μ represents the viscosity of fluid (kg/m3), and the pressure gradient ∇p is expressed as and
(10)∇p=i∂p∂x+j∂p∂y+k∂p∂z,∇2 in (9) is Laplacian operator described as (11). The left part of (9) represents the inertia force; in the right part, they represent inturn mass force, pressure (or pressure gradient), and the viscous force:
(11)∇2=∂2∂x2+∂2∂y2+∂2∂z2,
N-S equation accurately describes the actual flow. Viscous fluid flow analysis usually refers to researching N-S equation. Because of its rather complex form, only a handful can find the exact solution. The study of numerical solving is usually applied. N-S equation is also the basic equation for calculating the computational fluid dynamics. It can be said that all flow problems revolve around solving N-S equation.


### 2.3. Perfusion Numerical Simulations

According to the animal cells growth requirements, we need to endow micro fluid in the cell cavity with low shear force, good transmission effect, and so on. The truth is that there are always some mutual restraint factors among these principles. For example, a fully hybrid environment is required to strengthen its transmission effect during the dynamical culture process, while the vulnerability of animal cells restrict the formation of such an environment. Therefore theoretical analysis on the nature of the fluid mechanics of the design structure is carried out to optimize the constraints among the above-mentioned principles. The authors in [19] indicate that apart from some simple cases, it is difficult to give analytical solution to the most fluid mechanics problems, so computational fluid dynamics (CFD) is adopted in their analysis.In a calculation of CFD, the cell suspension is usually regarded as incompressible and homogeneous Newtonian fluid so that Turbulent Navier-Stokes mathematical model is applied to solve [20]. Reynolds number represents the ratio of fluid inertial force of unit mass and viscous force, discriminating the laminar and turbulent flow. Since Reynolds number in microfluidic scale is relatively smaller than 0.01 (Re<0.01), the pressure-driven flow through the cell culture unit is laminar flow [21]. Reynolds number can be calculated by
(12)Re=ρVdμ,
where V stands for the average flow rate, d is equal to characteristic length, also equal to hydraulic diameter, and then equal to 4 multiplies by overcurrent cross-sectional area divided by wetted perimeter. The relationship between the flow velocity and flow rate in working fluid inlet port is described as
(13)v=wA,
where v stands for flow rate, w is the inlet flow, and A is the cross-sectional area of the micro channel.We check the validity of the number of grids. We perform numerical simulation for five different grids at 604160, 791808, 946096, 1144500, and 1529928, respectively, with same inlet total pressures. We have the curve result of the flow rate for outlet total pressure loss as shown in Figure8(b). Three curves of grid nodes at 946096, 1144500 and 1529928 were almost overlapped, as could be seen in Figure 8(b). Thus it is demonstrated that the flow field is independent of the number of grids when the grid nodes are greater than or equal to 946096. Finally, the 1144500 grid nodes are applied, and that is in line with grid independency.The grids method: (a) the side view of mesh division for cell culture internal cavity, (b) the curve result of the flow rate for outlet total pressure loss.
(a)
(b)Based on the above discussion, “Numeca” software and Fine-solver are adopted in the simulation. Finite volume method (FVM) is employed to study flow and pressure distribution of micro fluid during perfusion in the microfluidic cell culture cavity. The streamline field and pressure field distribution in the cell culture cavity under different flow velocity and inlet static pressure are calculated by applying Numeca simulation analysis software. The height of the cell culture cavity is 2 mm, diameter 50 mm, and volume 15,708 mm3. The diameter of the fluid inlet and fluid outlet pipes is 1.7 mm. During establishing a physical mesh this cavity is divided into 1,144,500 griddings as shown in Figure 8, among which 1,020,285 grid nodes are used for cavity and 41,405 grid nodes for inlet and outlet pipes. We also increase the density of the solid wall mesh, and the mesh width of first layer is 0.001 mm. For the convenience of calculation, the material parameters of the working fluid adopt the material parameters of water (density at 998.3 kg/m3; the viscosity at 1.003 × 10−3 kg/(m·s)). Apart from inlet and outlet, the rest of cavity is seen as the solid wall boundaries of adiabatic and no-slip boundary. Firstly, simulation of large flow under the temperature at 15°C is carried out at inlet total pressures of 110,000, 120,000, and 150,000 Pa and the outlet static pressure is 101,000 Pa. The result is as shown in Figure 9. At the same time the simulation of the condition in greater pressure and larger flow rate is carried out and the result is shown in Figure 10.The simulation results in large flow: (a) the middle surface streamline chart in inlet static pressure of 200,000 Pa, (b) the middle surface streamline chart in inlet static pressure of 150,000 Pa, (c) the middle surface streamline chart in inlet static pressure of 120,000 Pa, (d) the middle surface vector of cavity in inlet static pressure of 150,000 Pa.
(a)
(b)
(c)
(d)The pressure distribution graph of cell culture cavity: (a) the middle surface pressure distribution graph in inlet static pressure of 200,000 Pa, (b) the pressure field distribution when the flow rate is 2 × 10−4 kg/s.
(a)
(b)The results indicate that when the outlet static pressure and temperature keep unchanged, the flow correlates with pressure: total inlet pressure of 110,000 Pa corresponds to flow rate of 2.791 g/s, 120,000 Pa corresponds to 5.279 g/s, 150,000 Pa corresponds to 10.03 g/s, and 200,000 Pa corresponds to 15.31 g/s. Secondly, the other things being equal, a simulation analysis at the flow rate 200 mg/s is carried out again. It can be seen from Figure11 that both the pressure field and the streamline field meet the requirements. Finally, the micro flow simulation result is shown in Figure 12.The simulation analysis when the flow rate is 200 mg/s. (a) Streamline distribution, (b) velocity distribution inx direction, (c) velocity distribution in y direction, (d) velocity distribution in z direction.
(a)
(b)
(c)
(d)The simulation results of micro flow: (a) the flow field distribution when the flow rate is 3 × 10−7 kg/s, (b) the flow field distribution when the flow rate is 3 × 10−5 kg/s, (c) the flow field distribution when the flow rate is 3 × 10−4 kg/s, (d) the flow field distribution when the flow rate is 2 × 10−7 kg/s.
(a)
(b)
(c)
(d)It can be seen from the streamline distribution graph (Figure9) that ejection happens in both of the two entrances. The differential pressure (or flow rate) is so great that the fluid injects to the inner wall of the upper cover. Obviously, this is against the protection of the cell structure, because the well-distributed streamline field means establishment of the cycle process of dynamical culture. It can be seen from the middle surface pressure filed distribution (Figure 10(a)) that the max differential pressure value in the cell culture cavity is 2000 Pa, which is too great to be conducive to cell culture. When the flow rate is 2 × 10−4kg/s, the max differential pressure value is only 12 Pa (Figure 10(b)) and the pressure field in the cavity is well distributed. Sharp change on differential pressure only happens in the fluid inlet and outlet. This meets the requirements. It can be seen, from the streamline distributions (Figures 12(b), 12(c), and 12(d)), that there is no infinite loop in the whole cavity area. Infinite loop emerges in the cell cavity when the flow rate is 3 × 10−7 kg/s (Figure 12(a)) and there is hardly pressure difference, which is regarded as the minimum flow rate of the designed cavity. Of course, further calibration is supposed to be in line with the actual experimental environment and the overall equipment characteristics.The different cell culture and experimental design demand different flow rate, for they have different cell culture volumes and different days in culture. Based on the experimental experience, the simulation of the feature points sampling is carried out to work out the maximum value and the minimum value, thus estimating the flow rate control range. So we roughly select the minimum of flow rate, that is, 3 × 10−7 kg/s, and the maximum of flow rate 200 mg/s according to the animal cells growth requirements. The result of Numeca analysis software also identifies that the structure design program of cell culture cavity is reasonable and the flow field is well distributed.


### 2.4. Experimental Device for Perfusion Cell Culture

Automatic guided space cell culture system aims to realize unattended cell culture, observation, and recording with microfluidic and temperature control circuits. This system involves high-speed analog signal synchronous sampling technique, embedded technology, intelligent power technology, servo control technology, multisensors information fusion, dedicated bus technology, SCADA technology, microfluidic hardware and software technology, and intelligence analysis technique. Major circuit diagram is shown in Figure13 consisting of Synchronous differential pressure sampling circuit, controller area network (CAN) and RS485 bus circuit, embedded control board, cell culture circuit, PCI-104 data acquisition module, the serial port expansion module, and processing analysis software. Cell culture and the extended port in the back panel of observation platform are as shown in Figure 14. These ports are mainly used to connect microscope, motor, and multisensor circuit board. As is shown in Figure 15, the multisensor data is used to acquire sampling unit pressure, temperature, biochemical parameter, and lighting control; the motor with CAN bus and RS232 port is used to control valve and micropump. The microscope of observation system has three grating rulers to control displacement in three directions. When we want to move over a direction for a certain distance, the servo motor powers it. The grating rulers acquire RS232 outputs pulse to calculate the precise displacement location, so that the microscope performs autofocus observation, other direction way similar. The communication between cell culture and observation system platform and microscope observation system uses three RS232 serial ports for the displacement controlling in x, y, and z directions and one USB port for the transmission of image data. Moreover, the micro pump closed-loop controlling of microfluidic system needs pressure sensor sampling value and two motors with RS232 and CAN. According to the simulation results and discussion of cell culture cavity, the microfluidic software flow chart is designed as shown in Figure 16. Temperature controlling system consists of temperature sensor PT100, micro controller C8051F, heating coil, and so on. and adopts PID algorithm. It adjusts the cavity temperature to cell culture.Figure 13
Cell culture and observation system platform diagram.Figure 14
A photograph of cell culture and observation system platform.Figure 15
A photograph of micro pump and multisensor data acquisition unit.Figure 16
Micro fluid controlling software flow chart.The human-machine interaction interface of cell culture and observation system platform adopts a 10.5 inch touch screen and then fixes it on the front panel of cabinet. Microsoft Visual C++ is adopted for its software programming.


## 2.1. Design and Fabrication of the Cavity

Entire cell culture cavity adopts fused silica, which has the following advantages: perfect thermal stability, low deformation, resistance to high acid and contamination, high hardness, and excellent light transmittance.Cell culture cavity (Figure3) is mainly made up of two parts: the upper cover and the lower cover. The upper and the lower covers are pressed together to form a hollow cell culture cavity. Cell suspension and reagent and metabolic waste go in and out through the three wedge-shaped cylindrical pipes welded to the lower cover. Both the upper and the lower covers have 50 mm diameter circular photic zone; the rest part is made as matte at roughness 3.2 and fastened via eight through-bolts in uniform distribution. The upper cover, lower cover, micro pressure sensor cavity, and three wedge-shaped catheters are shown in the assembly drawing of cell culture cavity (Figure 4).Figure 3
A photograph of cell culture cavity.Figure 4
The assembly drawing of cell culture cavity.To make the structure better sealed, the three wedge-shaped cylindrical pipes in the lower cover and three independent catheters are connected to micro pressure sensor cavity through the wedge connection, which contributes to overcoming the overload when entering orbit. Using Polytetrafluoroethylene (PTFE) tube connects inlet and outset into a rotary value, and then the rotary valve is connected to micropump, thus making it easy for the system to control micro fluid. The pipe in between serves as a waste export as well as a biochemical detection port. The overall structural design is ingenious and highly modularized.The system allows automatically real-time observation and recording in dynamical culture process through circular photic zone, which facilitates light irradiation and benefits the microscopic observation. The practical parameters of cell culture cavity are shown in Figure5. Heating coil (Figure 6) is designed in lower cover surface to provide a suitable temperature.Figure 5
The engineering drawing for upper cover of cell culture cavity.Figure 6
The engineering drawing for lower cover of cell culture cavity.What matters most in this experiment is how to design a perfect cell culture cavity and locate it precisely. A three-dimensional servo platform with high precision is designed to cooperate with microscopic observation. The top of the servo platform is the fixed device of cell culture cavity (Figure7).Figure 7
A photograph of cell culture cavity three-dimensional servo platform.


## 2.2. The Theory of the Computational Fluid Dynamics

Continuous equation, namely, the conversation of mass equation, means that any liquidity issue must agree with the law of conversation of mass. In accordance with the law of conversation of mass, the sum of the fluid rest mass out of control interval in unit time should be equal to the reducing quality caused by density variations in control interval in the same unit time. Based on that, the differential form of fluid flow continuous equation can be derived as(1)∂ρ∂t+∂(ρux)∂x+∂(ρuy)∂y+∂(ρuz)∂z=0,
where ux, uy, and uz are, respectively, the component of velocity (m/s) for three directions, t is time(s), and ρ is density (kg/m3).The Hamiltonian differential operator is derived(2)∇=i∂∂x+j∂∂y+k∂∂z.Equation (1) can be expressed as (3) and divergence form (4):
(3)∂ρ∂t+∇·(ρu→)=0,(4)∂ρ∂t+div(ρu→)=0.The form of the continuity equation is expressed as (5) in a cylindrical coordinate system:
(5)∂ρ∂t+ρurr+∂(ρur)∂r+∂(ρuθ)r∂θ+∂(ρuz)∂z=0.Equation (5) is the normal continuity equation for general use. If for steady flow ∂ρ/∂t=0, its form should be changed to
(6)∂(ρux)∂x+∂(ρuy)∂y+∂(ρuz)∂z=0.For incompressible flow, in the flow process, the densityρ of each dot is always unchanged, that is, ∂ρ/∂t=0 and ρ≠0; the form of the continuity equation of incompressible fluid can be derived as
(7)∂ux∂x+∂uy∂y+∂uz∂z=0.Cylindrical coordinate’s form of (7) is
(8)urr+∂ur∂r+∂uθr∂θ+∂uz∂z=0.The motion equation of viscous fluid was first proposed in 1827 by Navier who only took the flow of incompressible fluid into account. Poisson proposed the equation of the compressible fluid in 1831. Saint-Venant and Stokes, respectively, in 1843 and 1845, independently proposed the form in which viscosity coefficient was a constant, now known as Navier-Stokes equation, referred to as the N-S equation.For incompressible viscous fluid, the N-S equation is(9)∂u→∂t+(u→·∇)u→=f-1ρ∇p+μρ∇2u→.In (9), u→ (vector symbol) represents velocity (m/s), p represents pressure (Pa), μ represents the viscosity of fluid (kg/m3), and the pressure gradient ∇p is expressed as and
(10)∇p=i∂p∂x+j∂p∂y+k∂p∂z,∇2 in (9) is Laplacian operator described as (11). The left part of (9) represents the inertia force; in the right part, they represent inturn mass force, pressure (or pressure gradient), and the viscous force:
(11)∇2=∂2∂x2+∂2∂y2+∂2∂z2,
N-S equation accurately describes the actual flow. Viscous fluid flow analysis usually refers to researching N-S equation. Because of its rather complex form, only a handful can find the exact solution. The study of numerical solving is usually applied. N-S equation is also the basic equation for calculating the computational fluid dynamics. It can be said that all flow problems revolve around solving N-S equation.


## 2.3. Perfusion Numerical Simulations

According to the animal cells growth requirements, we need to endow micro fluid in the cell cavity with low shear force, good transmission effect, and so on. The truth is that there are always some mutual restraint factors among these principles. For example, a fully hybrid environment is required to strengthen its transmission effect during the dynamical culture process, while the vulnerability of animal cells restrict the formation of such an environment. Therefore theoretical analysis on the nature of the fluid mechanics of the design structure is carried out to optimize the constraints among the above-mentioned principles. The authors in [19] indicate that apart from some simple cases, it is difficult to give analytical solution to the most fluid mechanics problems, so computational fluid dynamics (CFD) is adopted in their analysis.In a calculation of CFD, the cell suspension is usually regarded as incompressible and homogeneous Newtonian fluid so that Turbulent Navier-Stokes mathematical model is applied to solve [20]. Reynolds number represents the ratio of fluid inertial force of unit mass and viscous force, discriminating the laminar and turbulent flow. Since Reynolds number in microfluidic scale is relatively smaller than 0.01 (Re<0.01), the pressure-driven flow through the cell culture unit is laminar flow [21]. Reynolds number can be calculated by
(12)Re=ρVdμ,
where V stands for the average flow rate, d is equal to characteristic length, also equal to hydraulic diameter, and then equal to 4 multiplies by overcurrent cross-sectional area divided by wetted perimeter. The relationship between the flow velocity and flow rate in working fluid inlet port is described as
(13)v=wA,
where v stands for flow rate, w is the inlet flow, and A is the cross-sectional area of the micro channel.We check the validity of the number of grids. We perform numerical simulation for five different grids at 604160, 791808, 946096, 1144500, and 1529928, respectively, with same inlet total pressures. We have the curve result of the flow rate for outlet total pressure loss as shown in Figure8(b). Three curves of grid nodes at 946096, 1144500 and 1529928 were almost overlapped, as could be seen in Figure 8(b). Thus it is demonstrated that the flow field is independent of the number of grids when the grid nodes are greater than or equal to 946096. Finally, the 1144500 grid nodes are applied, and that is in line with grid independency.The grids method: (a) the side view of mesh division for cell culture internal cavity, (b) the curve result of the flow rate for outlet total pressure loss.
(a)
(b)Based on the above discussion, “Numeca” software and Fine-solver are adopted in the simulation. Finite volume method (FVM) is employed to study flow and pressure distribution of micro fluid during perfusion in the microfluidic cell culture cavity. The streamline field and pressure field distribution in the cell culture cavity under different flow velocity and inlet static pressure are calculated by applying Numeca simulation analysis software. The height of the cell culture cavity is 2 mm, diameter 50 mm, and volume 15,708 mm3. The diameter of the fluid inlet and fluid outlet pipes is 1.7 mm. During establishing a physical mesh this cavity is divided into 1,144,500 griddings as shown in Figure 8, among which 1,020,285 grid nodes are used for cavity and 41,405 grid nodes for inlet and outlet pipes. We also increase the density of the solid wall mesh, and the mesh width of first layer is 0.001 mm. For the convenience of calculation, the material parameters of the working fluid adopt the material parameters of water (density at 998.3 kg/m3; the viscosity at 1.003 × 10−3 kg/(m·s)). Apart from inlet and outlet, the rest of cavity is seen as the solid wall boundaries of adiabatic and no-slip boundary. Firstly, simulation of large flow under the temperature at 15°C is carried out at inlet total pressures of 110,000, 120,000, and 150,000 Pa and the outlet static pressure is 101,000 Pa. The result is as shown in Figure 9. At the same time the simulation of the condition in greater pressure and larger flow rate is carried out and the result is shown in Figure 10.The simulation results in large flow: (a) the middle surface streamline chart in inlet static pressure of 200,000 Pa, (b) the middle surface streamline chart in inlet static pressure of 150,000 Pa, (c) the middle surface streamline chart in inlet static pressure of 120,000 Pa, (d) the middle surface vector of cavity in inlet static pressure of 150,000 Pa.
(a)
(b)
(c)
(d)The pressure distribution graph of cell culture cavity: (a) the middle surface pressure distribution graph in inlet static pressure of 200,000 Pa, (b) the pressure field distribution when the flow rate is 2 × 10−4 kg/s.
(a)
(b)The results indicate that when the outlet static pressure and temperature keep unchanged, the flow correlates with pressure: total inlet pressure of 110,000 Pa corresponds to flow rate of 2.791 g/s, 120,000 Pa corresponds to 5.279 g/s, 150,000 Pa corresponds to 10.03 g/s, and 200,000 Pa corresponds to 15.31 g/s. Secondly, the other things being equal, a simulation analysis at the flow rate 200 mg/s is carried out again. It can be seen from Figure11 that both the pressure field and the streamline field meet the requirements. Finally, the micro flow simulation result is shown in Figure 12.The simulation analysis when the flow rate is 200 mg/s. (a) Streamline distribution, (b) velocity distribution inx direction, (c) velocity distribution in y direction, (d) velocity distribution in z direction.
(a)
(b)
(c)
(d)The simulation results of micro flow: (a) the flow field distribution when the flow rate is 3 × 10−7 kg/s, (b) the flow field distribution when the flow rate is 3 × 10−5 kg/s, (c) the flow field distribution when the flow rate is 3 × 10−4 kg/s, (d) the flow field distribution when the flow rate is 2 × 10−7 kg/s.
(a)
(b)
(c)
(d)It can be seen from the streamline distribution graph (Figure9) that ejection happens in both of the two entrances. The differential pressure (or flow rate) is so great that the fluid injects to the inner wall of the upper cover. Obviously, this is against the protection of the cell structure, because the well-distributed streamline field means establishment of the cycle process of dynamical culture. It can be seen from the middle surface pressure filed distribution (Figure 10(a)) that the max differential pressure value in the cell culture cavity is 2000 Pa, which is too great to be conducive to cell culture. When the flow rate is 2 × 10−4kg/s, the max differential pressure value is only 12 Pa (Figure 10(b)) and the pressure field in the cavity is well distributed. Sharp change on differential pressure only happens in the fluid inlet and outlet. This meets the requirements. It can be seen, from the streamline distributions (Figures 12(b), 12(c), and 12(d)), that there is no infinite loop in the whole cavity area. Infinite loop emerges in the cell cavity when the flow rate is 3 × 10−7 kg/s (Figure 12(a)) and there is hardly pressure difference, which is regarded as the minimum flow rate of the designed cavity. Of course, further calibration is supposed to be in line with the actual experimental environment and the overall equipment characteristics.The different cell culture and experimental design demand different flow rate, for they have different cell culture volumes and different days in culture. Based on the experimental experience, the simulation of the feature points sampling is carried out to work out the maximum value and the minimum value, thus estimating the flow rate control range. So we roughly select the minimum of flow rate, that is, 3 × 10−7 kg/s, and the maximum of flow rate 200 mg/s according to the animal cells growth requirements. The result of Numeca analysis software also identifies that the structure design program of cell culture cavity is reasonable and the flow field is well distributed.


## 2.4. Experimental Device for Perfusion Cell Culture

Automatic guided space cell culture system aims to realize unattended cell culture, observation, and recording with microfluidic and temperature control circuits. This system involves high-speed analog signal synchronous sampling technique, embedded technology, intelligent power technology, servo control technology, multisensors information fusion, dedicated bus technology, SCADA technology, microfluidic hardware and software technology, and intelligence analysis technique. Major circuit diagram is shown in Figure13 consisting of Synchronous differential pressure sampling circuit, controller area network (CAN) and RS485 bus circuit, embedded control board, cell culture circuit, PCI-104 data acquisition module, the serial port expansion module, and processing analysis software. Cell culture and the extended port in the back panel of observation platform are as shown in Figure 14. These ports are mainly used to connect microscope, motor, and multisensor circuit board. As is shown in Figure 15, the multisensor data is used to acquire sampling unit pressure, temperature, biochemical parameter, and lighting control; the motor with CAN bus and RS232 port is used to control valve and micropump. The microscope of observation system has three grating rulers to control displacement in three directions. When we want to move over a direction for a certain distance, the servo motor powers it. The grating rulers acquire RS232 outputs pulse to calculate the precise displacement location, so that the microscope performs autofocus observation, other direction way similar. The communication between cell culture and observation system platform and microscope observation system uses three RS232 serial ports for the displacement controlling in x, y, and z directions and one USB port for the transmission of image data. Moreover, the micro pump closed-loop controlling of microfluidic system needs pressure sensor sampling value and two motors with RS232 and CAN. According to the simulation results and discussion of cell culture cavity, the microfluidic software flow chart is designed as shown in Figure 16. Temperature controlling system consists of temperature sensor PT100, micro controller C8051F, heating coil, and so on. and adopts PID algorithm. It adjusts the cavity temperature to cell culture.Figure 13
Cell culture and observation system platform diagram.Figure 14
A photograph of cell culture and observation system platform.Figure 15
A photograph of micro pump and multisensor data acquisition unit.Figure 16
Micro fluid controlling software flow chart.The human-machine interaction interface of cell culture and observation system platform adopts a 10.5 inch touch screen and then fixes it on the front panel of cabinet. Microsoft Visual C++ is adopted for its software programming.


## 3. Results and Discussion

In the present system, the cell’s requirement for nutrients and dissolved oxygen is met by the replacement of the culture medium in the cell culture cavity. When the flow rate is too large, it will not only cause the waste of the culture medium, but also the shear force and vortex that is bad for cells. The determination of the flow rate mainly relates to nutrients and dissolved oxygen supply and shear force. The supply of nutrients and dissolved oxygen is estimated empirically according to the required culture solution replacement time, which mainly depends on the structure and capacity of the culture vessel and cell type. The size of the shear force is related to the size of the flow rate as well as the structure of the culture vessel. Therefore, the flow rate of the culture solution is first estimated empirically according to required replacement time, and then we calculate whether the shear force caused by the flow rate meets the requirements. The appropriate flow rate is ultimately determined by experiment. Since M763 cells are cultured at36.5±0.5°C, temperature controlling circuits are expected to perform with a high-speed response and in a permitted error range. Adhering to the distance metering principle, traditional focusing technology mainly employs the distance system to drive the lens to find the best observation location, while this method needs additional equipment. Considering that space flight mission asks the equipment to be as small and light as possible, the initiative technology is employed. This technology uses image information to analyze focus direction and adjust the focal length to the best position. Since the focusing quality is decided by algorithm of the initiative technology [22], we adopt the morphological wavelet decomposition pyramid for microscopy system autofocusing.


### 3.1. Measurements

Two experiments are carried out on the temperature control unit and the result shows that the temperature of the cavity can be stabled. In experiment 1, the initial temperature is set as 27.68°C, the target temperature as 36.5°C, and temperature rising curve is as shown in Figure17. It can be seen that the time from the initial temperature to the target temperature is about 38 s and the temperature error is in the range of ±0.5°C. In experiment 2, the initial temperature is set as 31.65°C, the target temperature 36.5°C, temperature rising curve is as shown in Figure 18. it can be seen that the temperature remains stable after 25 s, so it is concluded that the rising of environmental temperature leads to the reduction of rising time.Figure 17
(a) Temperature characteristic of temperature closed-loop control circuit with start temperature at 27.68°C. (b) Temperature characteristic partial enlarged detail.Figure 18
(a) Temperature characteristic of temperature closed-loop control circuit with start temperature at 31.65°C. (b) Temperature characteristic partial enlarged detail.72-hours-cell culture and observation are carried out with the aid of the culture cavity. Microscope (Nikon ECLIPSE Ti-S, Nikon Corporation, Japan), CCD (Olympus), and the cell culture and observation system platform form a cell image observation system. Figure19 is the growth records of M763 in the cavity captured by the image observation system. We can see that the cells in 24 hours, have stuck to the wall and unfolded completely. In 48 hours, apparent cell division happens. In 72 hours, the process continues. All these indicate that the system design in this paper is successful.Growth image of M763 cells for different time steps. (a), (b), and (c) are, respectively, growth images of 24 h, 48 h, and 72 h.
(a)
(b)
(c)


### 3.2. System Performance

The automatic guided space cell culture system has been calibrated. Part of its performance indexes are as shown in Table2.Table 2
The part index of automatic guided space cell culture system.
Parameter
Index
Microscopic image system resolution
0.98μm
Maximum pixel
2048
*
1536
horizontalCx = 138.7 nm/pixel;
Image pixel equivalent
VerticalCy = 139.7 nm/pixel
Frame rate
5 frame/s
Mobile platformx direction resolution
0.5μm
Mobile platformy direction resolution
0.6μm
Compression rate (for umbilical vein living cells image)
39.4 : 1
Range of recommended flow-rate
0.3 µL/s~200 µL/s
Precision of temperature control
±1°C


## 3.1. Measurements

Two experiments are carried out on the temperature control unit and the result shows that the temperature of the cavity can be stabled. In experiment 1, the initial temperature is set as 27.68°C, the target temperature as 36.5°C, and temperature rising curve is as shown in Figure17. It can be seen that the time from the initial temperature to the target temperature is about 38 s and the temperature error is in the range of ±0.5°C. In experiment 2, the initial temperature is set as 31.65°C, the target temperature 36.5°C, temperature rising curve is as shown in Figure 18. it can be seen that the temperature remains stable after 25 s, so it is concluded that the rising of environmental temperature leads to the reduction of rising time.Figure 17
(a) Temperature characteristic of temperature closed-loop control circuit with start temperature at 27.68°C. (b) Temperature characteristic partial enlarged detail.Figure 18
(a) Temperature characteristic of temperature closed-loop control circuit with start temperature at 31.65°C. (b) Temperature characteristic partial enlarged detail.72-hours-cell culture and observation are carried out with the aid of the culture cavity. Microscope (Nikon ECLIPSE Ti-S, Nikon Corporation, Japan), CCD (Olympus), and the cell culture and observation system platform form a cell image observation system. Figure19 is the growth records of M763 in the cavity captured by the image observation system. We can see that the cells in 24 hours, have stuck to the wall and unfolded completely. In 48 hours, apparent cell division happens. In 72 hours, the process continues. All these indicate that the system design in this paper is successful.Growth image of M763 cells for different time steps. (a), (b), and (c) are, respectively, growth images of 24 h, 48 h, and 72 h.
(a)
(b)
(c)


## 3.2. System Performance

The automatic guided space cell culture system has been calibrated. Part of its performance indexes are as shown in Table2.Table 2
The part index of automatic guided space cell culture system.
Parameter
Index
Microscopic image system resolution
0.98μm
Maximum pixel
2048
*
1536
horizontalCx = 138.7 nm/pixel;
Image pixel equivalent
VerticalCy = 139.7 nm/pixel
Frame rate
5 frame/s
Mobile platformx direction resolution
0.5μm
Mobile platformy direction resolution
0.6μm
Compression rate (for umbilical vein living cells image)
39.4 : 1
Range of recommended flow-rate
0.3 µL/s~200 µL/s
Precision of temperature control
±1°C


## 4. Conclusion

In this paper, a fused silica space cell culture cavity is designed, simulated, and tested. This culture cavity with circular photic zone for autoobservation and heating coil for adjusting temperature realizes in-place culturing, observing and recording, with the aid of automatic control system having microfluidic and temperature control circuits. The efficacy of the wedge-shaped connection is demonstrated; however, for future space flight mission it should be more robust to perform better leakproofness and be better resistant to the impact of the gravitational acceleration caused by the duration of carrier rocket into orbit. It has also been demonstrated that automatic guided space cell culture system with microfluidic and temperature control circuits is successful in laboratory environment.The system aims to be used in the future space experiment, and extensive ground testing has been done to verify its functions. It performs perfectly during the test. The only flaw is that the system is too bulky for a space experiment. The future work will focus on the miniaturization of the apparatus.



---

*Source: 378253-2013-12-26.xml*