APCChE 2019

Hall P, evening, Day 3, Wednesday Abstracts

Session 2. Fluid and particle processing

PB301 Numerical simulation of particle suspension in a stirred tank equipped rigid-flexible impeller
Xia XIONG, Zuohua LIU, Liang WANG, Changyun TAO
School of Chemistry and Chemical Engineering, Chongqing University, Chongqing

Particle suspension characteristics are predicted computationally in a stirred tank with solid concentration of 25% and driven by rigid-flexible impeller. The software of Fluent was utilized to investigate flow field structures in stirred tank. The effects of impeller type, length of flexible material, connect way of flexible material, mixing time were investigated by using the Eulerian–Eulerian two-fluid model and the standard k–ε turbulence model. It turned out that the structure of impeller has an obvious effect on solid particle distribution. The suspension of solid particles in the rigid-flexible impeller system was more homogeneous. It was also suggested that longer flexible material and longer mixing time results in well homogenous distribution of solid particles in the tank equipped with rigid-flexible impeller. On the whole, the rigid-flexible impeller could enhance the suspension of solid particles effectively.

PB302 Numerical simulations of flow and heat transfer in a industrial hydro thermal cracking reactor of residue oil
Zhenxing ZHU, Xiaojin TANG, Zijun WANG
SINOPEC Research Institute of Petroleum Processing, Beijing, China

Hydro thermal cracking process of residue oil is an important technical route in the area of heavy oil processing. Because of a few catalyst loading, low cost, and low blocking probability, the slurry bed is better to be applied in processing residue oil, especially for heavy residue. With the increase of conversion, the stability of residue oil becomes weaker, and some unstable compounds may be crystallized easier. The precipitates will be attached at the surface of catalysts and reactor, which will lead to the deposition of catalyst and blocking. The distribution of velocity and temperature inside the reactor is essential for reducing precipitation of unstable compounds, and preventing reactor from blocking. Based on the results of a pilot plant, an industrial hydro thermal cracking reactor of residue oil was developed. Then, the flowing and distribution of fluid and deposition of catalysts were simulated by CFD. To reduce the deposition of catalysts inside the large slurry bed reactor, the position of inlet of residue oil was optimized. Furthermore, the distribution of temperature inside the reactor was predicted by CFD. It can be found that the distribution of temperature inside the reactor was very uniform, the difference in axial temperature is less than 3K, the difference in radial temperature is less than 1K.

PB303 Study on the Effects of Grinding Factor and Scaling up in Grinding Using Horizontal Dry Bead Mill
Takahiro TAMURA1, Mitsumasa KIMATA2
1 Institute of Fine Particle Technology, Ashizawa Finetech Ltd, Oyama, Japan
2 Yamagata University, Yonezawa, Japan

Though the grinding operation is a unit operation that has been known for a long time, it still often occupies an important position in industry. Thus, various grinders have been evolved and produced in a variety of shapes to meet changing needs. In principle, bead mills can provide grinding conditions under a high packing density and high energy density of the media as well as cover a broad range of output. Therefore, their future application is highly expected to respond to the requirements of various grinding operations by the current industry. However, bead mills have a short history compared to other types of grinders. Especially for dry bead mills, it is hard to say that they are well known. That is why we studied and reported the performance of horizontal dry bead mills in this presentation.
For the experiments, 2 sizes of horizontal dry bead mills with different capacities, lab and manufacturing scale, were used to study scaling up and down. The grinding experiments were performed by changing each typical grinding factor for bead mills, including the agitator peripheral speed, bead packing density, bead diameter, and material feeding speed. This set of experiments was conducted for 3 kinds of materials, silica sand of natural inorganic materials, alumina of synthetic inorganic materials, and activated carbon of organic materials.
As a result, common for all materials, similar rules were found between the particle size and specific energy, the amount of mechanical energy per unit amount of material. In addition, these rules were applicable for the different-size mills, indicating that the scale can be controlled by the rules and specific energy. Also, it is revealed that when using a single type of material, changing the bead diameter changes the energy efficiency, indicating that using small-diameter bead was effective for efficient grinding.

PB304 The facile control synthesis of hexagonal NaREF4 in a continuous microfluidic reaction system
Di LIU, Kai WANG, Yundong WANG
The State Key Laboratory of Chemical Engineering, Departmemnt of Chemical Engineering, Tsinghua University, Beijing 100084, China

In an effort to achieve the highly efficient synthesize upconversion nanoparticles (UCNPs) with controlling over properties, we report for the first time a continuous flow synthesis of hexagonal β-phase NaREF4 NPs using microflow system. The precursor solution was prepared by the reaction of rare earth oleate, sodium oleate and ammonium fluoride, which avoided the production of NaF at room temperature and speeded up the hexagonal nanoparticle growth rates at high temperature. Specifically, the separation of nucleation stage and growth stage was achieved by means of segmented temperature control in the microflow system. The nucleation and growth reaction mechanism were further studied and the experiments of using different temperature, reactant concentration and flow rate were designed to study the influence on the shape and size of β-phase NaREF4. The ultra-small (5-6nm) β-NaYF4:Yb,Er nanoparticles were prepared by Gd3+ doping in the microflow system. And the monodisperse sub-10 nm core-shell nanoparticles β-NaYF4(Gd):Yb,Er@NaYF4 were prepared successfully by mixing the core solution and the shell precursor solution via a continuous microfluidic process. The β-phase nanoparticles exhibit green and red emissions under excitation of laser at 980 nm. Our research opens a door for the synthesis of upconverting nanocrystals in the continuous flow reaction system.

PB305 Growth of a pinhole generated in a thin coated liquid film
Takumi CHIBA1, Ryota NAKAJIMA1, Atsuya TAKADA1, Shunji HOMMA1, Donghyuk KANG1, Hiroyuki HIRAHARA1, Yosuke OKUBO2, Hirokazu UMEMIYA2, Sumihisa ODA2
1 Saitama University, Saitama City, Saitama, Japan
2 Saiden Chemical Industry Co., Ltd, Saitama City, Saitama, Japan

The growth of a pinhole generated in a thin coated liquid film was simulated numerically in order to study the mechanism of dewetting on a thin coated adhesives. Navier-Stokes equations of a one-fluid formulation were solved by projection method and the motion of the liquid-air interface was traced by a VOF method. Rheological properties of the liquid were approximated by a Carreau-Yasuda model. The OpenFOAM tool box was utilized to perform the numerical simulations. Simulated results were in good agreement with experiments where a pinhole was spread on a coated film of adhesives. The growth rate was proportional to the inverse of Ohnesorge number and therefore it decreases when the viscosity of the liquid is high and the surface tension is low.

PB306 Three-dimensional simulations of laminar film condensation
Saitama University, Saitama, JAPAN

Laminar film condensation of saturated vapor on the surface of a low-fin tube was numerically simulated in order to examine the effect of the shape of the tube. Navier-Stokes equations, an energy equation, and a continuity equation including phase change were solved simultaneously. Vapor-liquid interfaces were traced by VOF method. The OpenFOAM tool box was utilized to perform the numerical simulations.
Figure 1 shows an example of the simulation results. When the fin pitch is large, the condensate drips as a tuft from the lower part of the fin. When the fin pitch is small, on the other hand, the condensate falls as a curtain after filing the condensate between the fins.
The rate of condensation was obtained from the numerical simulations. When the fin pitch is small, the rate of condensation increases with the fin pitch. However, the rate decreases with the pitch when the pitch is larger than the certain value. It has thus suggested that there is an optimum fin pitch. Numerical simulations are expected to be able to optimize the shape of low-fin tubes.

PB307 Evaluation of drying process of a droplet considering flow and temperature fields in a spray dryer
Ryuya ITANO, Tomohisa KAWAGUCHI, Koichi NAKASO, Yasushi MINO, Kuniaki GOTOH
Okayama University, Okayama, Japan

Particle generation process from droplet in a spray dryer is investigated. In the spray dryer, droplets generated from an atomizer are introduced to the dryer with hot air. While droplets are heated in the airflow, the air is cooled by exchanging heat with droplets. The atmosphere of the droplet such as temperature and humidity therefore changes with time in the dryer. Because drying behavior is influenced by the temperature and humidity around the droplet, drying of droplet should be evaluated by considering change in atmosphere of droplet. However, the drying of droplet has been studied under fixed conditions in the past researches. In this study, to observe the drying behavior, a droplet suspended on the tip of a wire is heated by the controlled air flow.
A slurry droplet of silica particles in water is studied. The effect of flow rate on the drying of droplet is investigated because flow around the droplet would affect the local heat transfer rate. Fig. 1 shows examples of the drying results with slurry droplet. At higher flow rate, mushroom-shaped particle are obtained. On the other hand, doughnut-like particle is obtained at lower flow rate. The results could be explained in terms of the particle Reynolds number.
The solution droplet of dextrin in water is also investigated. The effect of temperature on the drying is especially studied. As the drying proceeds, a solute layer is formed on the surface of the droplet. The permeability of water vapor through the solute layer is low especially for the use of dextrin. When the temperature is higher than the boiling point of the solution, boiling happens inside of the droplet then expanded particle is formed. On the other hand, the particle containing water inside of the particle is formed for the lower temperatures.

PB308 Evaluation of gas dispersion in stirred tank with fractal shaped impellers
Fukuoka University, Fukuoka, Japan

The gas dispersion in stirred tank is a widespread operation in the chemical and biochemical industry. In the process of gas dispersion, control of bubble size is important factor for the mass transfer process, and the performance relies heavily on its impeller design. Rushton turbine (RT) has been commonly used for the process because of its strong shear force. Thus, we proposed a novel kind of impellers, what we call fractal shaped impellers based on the RT. A kind of fractal structure, Koch curve was applied on blade tips of RT as shown in Figure. We expected that increase of the length of a side by iterations of fractal structure promoted shear around the wingtip and generated multiscale vortex to break up bubbles in gas dispersion.
Consequently, we investigated the effect of fractal shaped impellers on mass transfer process in gas dispersion. A stirred tank of transparent acrylic with four baffles was used for the experiment. An optical dissolved oxygen sensor was located in the front face of a baffle to evaluate the volumetric mass transfer coefficient (kLa) from variation of oxygen concentration in the tank. The bubbles discharged from a blade were visualized by high speed camera. The bubble size was measured by both photographic method and image analysis.
First, the kLa of each impellers was evaluated under the same power consumption. The kLa increased by using KF-I type impellers and KF-Ⅱtype impellers didn't make a difference. The generate bubble size also became smaller by using KF-I type impellers. Finally, we confirmed that the more small bubble generated from the blade, the more kLa increased in the tank.

PB309 Experimental investigation of gas slug shrinking behaviour inside a microchannel
Yusuke NOGUCHI1, Satoko FUJIOKA2, Koichi TERASAKA2, Simon MATTHES3, Michael SCHLÜTER3
1 School of Science for Open and Environmental Systems, Graduate School of Science and Technology, Keio University, Yokohama, Japan
2 Department of Applied Chemistry, Faculty of Science and Technology, Keio University, Yokohama, Japan
3 Institute of Multiphase Flows, Technische Universität Hamburg Harburg, Hamburg, Germany

Gas-liquid slug flow is a flow pattern in which gas and liquid flow alternately.
There are numbers of prior researches on enhancement of extraction and reaction using the mixing effect of circulating flow inside a liquid slug, however, there are few experimental examples focusing on the shape change of gas slug due to dissolution.
This investigation focuses on the shape changing of gas slug in order to control slug size by changing flow rate and channel geometry to apply it for size-controlled micro bubble production.
Liquid phase was degassed and then pressurized from a reservoir to T-junction microchannel made of PDMS with controlling flow rate. Oxygen gas was injected to the microchannel through the gas-injection nozzle by using syringe pump. Size-controlled gas-liquid slug was produced by injecting both fluid and observed with a close focus lens above the microchannel. The change of slug volume and shape effected by liquid flow rate were investigated from captured images.
Figure shows the time course of oxygen slug shrinking behaviour in condition of no flow of gas and liquid phase after slug production.
The change rate of slug length for flowing condition of liquid phase was larger than that of no flow condition, therefore, liquid phase flow promoted gas slug dissolution.
Without degassing of liquid phase, shrinking of gas slug was not observed. It is thought that dissolved gas concentration is important parameter to control gas slug size.
By measuring slug size, the possibility to control the gas slug size was indicated and precisely controlled microbubble production is expected using this method by injecting size-controlled gas slug into vessels filled with liquid.
This method for microbubble production might be the solution to overcome difficulties of small bubble detaching in conventional bubble formation from a nozzle or a slit.

PB310 Evaluation of Particle-Wall Adhesion Characteristics Based on Separation of Adhered Particle by High Speed Air Jet
Yuna ISHIDA, Maho KAGEYAMA, Yasushi MINO, Koichi NAKASO, Kuniaki GOTOH
Okayama University, Okayama, Japan

In many processes in which particles are operated, adhered particles on the wall of equipment and pipes cause problems such as contamination and blockage in feeding or transportation. Therefore, it is required to reflect the adhesion characteristics in the process design and the determination of operating conditions. However, the designing method and the prediction method of operating conditions with taking into the adhesion characteristics have not been established yet.
Because a particle behavior in these process operation is affected by fluid motion, the evaluation of adhesion characteristics based on the force induced by fluid motion is useful for establishment of the designing or the prediction method. As one of simple operation relating to the adhesion and fluid motion, we focused on the air jet removal of particles. In the air jet removal, adhered particles are separated from a wall by high-speed air flow ejected from a nozzle. The removal efficiency of the air jet removal depends on the dynamic pressure of the impinging air flow on a wall. It suggests that the separation force acting on a particle is determined by the dynamic pressure.
In this study, a slit nozzle having a width of 0.25 mm and a length of 20 mm was used as a reference nozzle, and 7-types of spherical silica particles having different particle diameter ranges were used as a reference particle. The removal efficiency of these particles from the borosilicate flat glass plate was determined experimentally. From the correlation between the diameter of reference particles and the required dynamic pressure for the removal, the separation force acting on reference particles was determined as a function of particle diameter and dynamic pressure. By means of these reference particles known to the separation force, the adhesion between the various wall surfaces and the reference particles was determined.

PB311 Droplet breakage in a pulsed disc and doughnut column
Shaowei LI, Shan JING, Hao ZHOU, Xiong YU
Institute of Nuclear and New Energy Technology, Tsinghua University, Beijing 100084, China

Lack of experimental data for droplet breakup is one of the limitations for the application of population balance model (PBM). We present here a new method to directly measure the breakage of droplet swarm in a pulsed disc and doughnut column (PDDC). A high speed camera was used to directly measure the droplet breakup frequency and daughter size distribution for a series of liquid-liquid systems. The multiple breakup was treated as an original breakup and several intermediate breakups in order to characterize the process quantitatively. The effects of pulsation intensity, dispersed phase viscosity and flow rate, interfacial tension were investigated in detail. The data covers the main factors affecting the droplet breakage in the PDDC and are sufficient to support the establishment of useful empirical correlations for breakup frequency and daughter droplet size distribution. The experimental results and correlations all indicate that the breakup frequency increases with the decreasing of the dispersed-phase viscosity or the interfacial tension, or the increasing of the pulsation intensity. The daughter droplet size distribution broadens considerably with decreasing interfacial tension. The correlation equations were then used in a CFD-PBM simulation to calculate the droplet number density, which further proved the feasibility of the correlations.

PB312 Whipping behavior of rice flour dough
Chiaki ICHIKAWA, Takahiro FUKUNAGA, Daitaro ISHIKAWA, Tomoyuki FUJII
Graduate School of Agricultural Science, Tohoku University, Sendai, Japan

The bubbles in rice flour dough were investigated under constant temperature. The change of distribution in bubble size plays an important role in controlling food texture. If the formed bubble size depended mainly on three kinds of force, such as inertial, viscous, and surface tension, then the normalized mean bubble diameter should be a function of Reynolds Number and Weber Number. We obtained our experimental data mainly using a home size whipper, and we compared the data with the results obtained from the dough with nine kinds of rice flours prepared by different milling process in this work. We also added the size effect of rice flour particle as Bond Number. Furthermore, we carried out a dynamic wettability test in order to estimate wettability of rice flour surface. The results of the dynamic wettability test were described well by Washburn equation, and dccosθ was calculated as a parameter of wettability (dc : effective diameter of capillary in powder bed, θ: contact angle). The relationship between Bo and dccosθ/dp, normalized wettability parameter by particle size (dp) of rice flour, was obeyed by a straight line. Thus, it was suggested that the powder properties of rice flour could classify by interface properties. The mean diameter of the bubbles (dbm) generated by whipping was expected to be affected by the thickness (d) of the rod of whipper, its movement speed, physical properties of material and gravitational acceleration. The dimensionless mean diameter (dbm/d) was expressed by a dimensionless equation. It was found that We was high contribution for mean diameter from the result of its power index obtained from the dimensionless equation demonstrated. The empirical equation obtained was generally applicable to the variety of materials selected. Consequently, the bubbles in rice flour dough generated by whipping process could evaluate with the interface properties of rice flour.

PB313 Rheological properties of Nitrile Rubber Latex blended with Titanium Dioxide (TiO2)
Tiam-Ting TEE1, Siew-Wai PHANG2, Tin-Sing LEE1, Soo-Tueen BEE1, Heng-Wai KANF1
1 Universiti Tunku Abdul Rahman, Kuala Lumpur, Malaysia
2 Tayloys University, Kuala Lumpur, Malaysia

Nitrile Rubber (NBR) latex is used for making various gloves. The rheological properties of Nitrile Rubber Latex blended with Titanium Dioxide (TiO2) were investigated using a rheometer. By carry out rotational and oscillatory tests, the relationships of viscoelastic properties such as viscosity, complex viscosity, storage modulus, loss modulus and phase angles were investigated.
Rotational test determine the relationship of viscosity and shear rate. Respond on the viscosity was investigated by increasing the shear rate. The samples under studied were pure NBR, compound NBR and compound NBR with different phr value of filler. In these tests, the viscosity of the NBR decreased with increasing shear rates.
Oscillatory tests were carried out to determine the relationship between complex modulus and angular frequency. The higher the angular frequency applied, the lower the complex viscosity were obtained for the samples. Viscoelastic properties such as storage modulus, loss modulus and phase angles were studied. The storage modulus of the sample decreased when angular frequency increased. Whereas the loss modulus increased with increasing angular frequency. This indicates that the fluids are more favourable on energy loosing rather than energy storing. By referring to the phase angles, the fluid's flow behaviour can be easily known. 0 °phase angle indicates the fluid exhibit perfect elastic behaviour, whereas 90 °phase angle indicates that the fluid is exhibiting Newtonian flow (water behaviour). NBR compounds shows viscoelastic behaviour

PB314 Application of DEM simulation to design of a laboratory-scale high pressure grinding roll
1 Waseda University, Shinjuku, Tokyo, Japan
2 National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, Japan
3 Furukawa Industrial Machinery Systems Co., Ltd., Tochigi, Tochigi, Japan
4 Furukawa Co., Ltd., Tsukuba, Ibaraki, Japan
5 Nittetsu Mining Co., Ltd., Nishitama, Tokyo, Japan

It is known that high pressure grinding roll (HPGR) is an effective grinding method to enhance the liberation of minerals with high energy efficiency. HPGR grinding has many kinds of operating parameters which affect the liberation state and the optimum operating parameters strongly depend on a types of ores. From the viewpoint of industrial application of HPGR in mineral processing, it is essential to decide the optimum operating conditions depending on ores. Since a batch of HPGR experiment in pilot scale requires huge amount of ore samples, it is difficult to conduct and repeat experiments in various operating conditions. Therefore, it is highly anticipated to newly develop a laboratory-scale HPGR which only requires small amounts of ores at an early stage of mine development.
The objective of this study was to develop a laboratory-scale HPGR which can reproduce grinding results similar to those obtained from a pilot-scale HPGR. To achieve this objective, the discrete element method (DEM) coupled with T10 breakage model was applied to reproduce a HPGR grinding process. First, grinding processes of ores in a pilot-scale HPGR were simulated by the DEM coupled with T10 breakage model. Comparison between simulation results and experimental results was also conducted for the validation of the simulation model. Good agreement was confirmed for throughput of ore, energy consumption and a working roll gap. After the validation of simulation model for a pilot-scale HPGR, this simulation model was applied to a laboratory-scale HPGR. An appropriate shape of a roll installed in a laboratory-scale HPGR was investigated. Simulation results indicated that the suitable roll shape for a laboratory-scale HPGR was designed by a dimensionless roll gap and dimensionless throughput. A laboratory-scale HPGR designed based on simulation results could experimentally reproduce grinding results similar to those obtained from a pilot-scale HPGR.

PB315 Effects of molecular weight on drag force of polyethyleneglycol in flows measured by a scanning probe microscope
Kobe University, Kobe, Hyogo, Japan

Very dilute high-molecular-weight polymer solutions induce complex flow behaviors. This is due to polymer deformation in a flow, which causes mechanical interaction between a deformed polymer and fluids, and between polymers in flows. We consider such polymer-fluids interaction and polymer-polymer interaction clarify the complex behaviors of polymer solution.
In this study, we measured the drag force of polyethyleneglycol (PEG) in flows by using a scanning probe microscope (SPM). Methoxy polyethyleneglycol thiol (mPEG-SH) was attached to the cantilever probe of the SPM, which was further immersed in flows of glycerol and polyethyleneglycol (PEG) solutions. Here, we have prepared several molecular weight of mPEG-SH from 10k to 40k. The mPEG-SH-bonded cantilever detects the extra force due to polymer-fluids and polymer-polymer interaction in flowing fluids. The extra force was increased with increasing the molecular weight. The conformation of the mPEG-SH polymer bonded to the probe of the cantilever was predicted, and the drag force due to the deformed mPEG-SH was calculated. The forces detected by experiments using the SPM and the forces obtained by model calculations were compared, and found to be reasonably close.

PB316 Mass transfer in alternating flow of two immiscible liquids in microchannel
Mayu NISHIYAMA, Takatomo TSUDA, Hiroyuki HIRANO, Tomohiko KIHARA, Hideaki OBATA, Kenya KUWAGI
Okayama University of Science, Okayama-shi, Okayama 700-0005, Japan

The mass transfer of the liquid-liquid extraction process in the alternating flow of two immiscible liquids in the microchannel was experimentally visualized. In this flow, each segment flows alternatively in the microchannel. In this study, water with congo red as the indicator, and cyclohexane with acetic acid were induced into the microchannel with the rugged wall of the diamond shape. The alternating flow of these two kinds of liquids was observed, and the mass transfer of acetic acid from aqueous phase into organic phase was visualized as the change of the color of aqueous phase. Figure 1 shows the experimental result. The microchannel was made of stainless-steel and placed between two quartz glass plates. The narrowest width of the microchannel is 0.3 mm, and the widest one is 1.45 mm in the diamond shape. Two kinds of liquids were fed into the microchannel with the syringe pumps at 0.07 mL/min (13 mm/s). Orange area is aqueous phase of water with congo red, white area is organic phase of cyclohexane with acetic acid. Blue area is aqueous phase of water, congo red and acetic acid which transfers from organic phase, and the color of aqueous phase changes from orange to blue as the concentration of acetic acid increases. The color of the aqueous phase completely changed from orange to blue after 9.95 mm from the confluent point. As the result, the distance from the confluent point where the color changed completely from orange to blue becomes longer as the inlet velocity increases. This result suggests that the ruggedness of the wall is not effective for this extraction process. This work was supported by the Grant for the Promotion of OUS Research Project (OUS-RP-30-4) of Okayama University of Science.

PB317 Influence of The Injection Point on The Mixing Process with A Large-type Impeller in Laminar Mixing
Keisuke SAKATA, Kazuhiro OJIMA, Kazuhiko NISHI
Chiba Institute of Technology

In laminar mixing with a large-type impeller, the initial injection point of the mixing object is known to influence mixing processes and the mixing completion time. This study investigated the relation between the injection points and the mixing performance of MAXBLEND, a large-type impeller, by CFD of tracer particles.
For this study, a 0.1 m diameter vessel was used. 1331 particles were placed on an injection point as a 53 mm3 cube. Based on the three-dimensional particle coordinates for every time, the distance from each particle to the nearest particle is calculated (Lnp). The average of these distances was calculated for every time (Lnp, ave).
In all, we investigated 48 injection points (Fig. 1). In Fig. 1, the results of the numerical simulation of mixing are shown as a relation between time and Lnp, ave . As shown there Lnp, ave gradually becomes larger with advanced mixing from 1 mm (t = 0). It becomes almost constant in a complete mixing state. For injection point No. 32, the mixing completion time was the shortest. However, the mixing completion time of No. 38 was the longest. Results showed a large difference in the complete mixing time depending on the injection point.
To investigate the scale-up effect, mixing in the 0.15 m diameter vessel was simulated. When the Reynolds number was equal using the same liquid, the relation between Lnp, ave and time in a large vessel differs from that in a small vessel.
However, these data are arranged similarly to the lines, independently of scale, using dimensionless mixing time, which is the product of the time and the impeller rotational speed, and the normalized average nearest particle distance, which is calculated by dividing Lnp, ave by the vessel diameter. Results demonstrated that mixing in the large vessel can be estimated from results of the small vessel using this relation.

PB318 Particle dispersion and aggregation under various shear fields
Katsuki MURANISHI1, Hayato MASUDA2, Takashi HASEGAWA3, Takafumi HORIE1, Naoto OHMURA1
1 Department of Chemical Science and Engineering, Kobe University, 1-1 Rokkodai, Nada, Kobe, Japan
2 School of Food and Nutritional Science, University of Shizuoka, 52-1 Yada, Suruga, Shizuoka, Japan
3 Toagosei Co., Ltd., General center of research and development, 8 Showamachi, Minato, Nagoya, Aichi, Japan

Fine particles make an important role for the production of functional materials in various industries such as foods, pharmaceuticals and paints. Because physical properties of fine particles depend on the particle size and the particle size distribution, it is important to accurately control these properties during the processing. Recent developments in functional materials have been required for more accurate control of the particles size and the particle size distribution.
Generally, fine particles are produced using a stirred tank with the aim of large-scale production. However, in the stirring tank, fine particles easily aggregate due to shear in a flow filed. Therefore, in order to control the particle size and the particle size distribution, the effect of shear on particle aggregation should be understood. However, there is the shear force (shear-rate) distribution in a stirred vessel and it differs with the impeller type. It is necessary to clarify how the difference in the impeller type affects the particle size and the particle size distribution. This study investigated the effect of impeller type on the size distribution of aggregate of fine particles.
Two types of impeller, a Rushton turbine and a 3-retreat impeller were used for particle aggregation. The size distribution of aggregates was measured by a microscope (VHX-100, Keyence Co., Ltd.) and a laser diffraction particle size analyzer (SALD-300V, Shimadzu Corp.).
It has been found that the 3-retreat impeller promotes particle aggregation as compared with the Rushton turbine impeller. Besides, energy dissipation with the 3-retreat impeller was larger than it with the Rushton turbine impeller. These results imply that an impeller with smaller energy dissipation would suppress particle aggregation. In the conference, the results using a Taylor-Couette flow apparatus having the relatively narrow shear-rate distribution will be also reported.

PB319 Delivery of hydrophilic materials into tobacco BY-2 cells using PLGA nanoparticles
Hikari AISU, Yasuhiro KONISHI, Toshiyuki NOMURA
Department of Chemical Engineering, Osaka Prefecture University, Osaka, Japan

During the past decades, nanoparticles (NPs) have been widely studied as carriers of drug delivery system. Unlike the animal cells, plant cells have a rigid cell wall, and it is physical barrier for intracellular NP uptake. However, past our study revealed that the poly lactic-co-glycolic acid (PLGA) NPs were taken into the tobacco BY-2 cells with cell wall by controlling the NP exposure condition. In this study, we attempted to deliver PLGA NPs containing hydrophilic biomaterials into the tobacco BY-2 cells. The fluorescent BSA was used as model biomaterials. BSA adsorbed PLGA NPs were prepared as follows: PLGA NP was prepared by the emulsion solvent diffusion method. In order to adsorb negatively charged BSA on the PLGA NPs, the cationic lipid DOTAP was adsorbed on the PLGA NPs and BSA was then adsorbed on them. BSA encapsulated PLGA NPs were prepared as follows: In order to encapsulate hydrophilic BSA in hydrophobic PLGA NP, BSA was hydrophobized using DOTAP and the BSA/DOTAP complex was then encapsulated in PLGA NPs. BY-2 cells and PLGA NPs were mixed in MS medium for 1 hour. Then, the localization of NP and BSA was observed using CLSM. When exposing BSA adsorbed PLGA NPs, the PLGA NPs were taken up into the cells but BSA was remained on the cell wall. This is presumably because the adsorption of BSA on the NP surface reduced the cell affinity of PLGA. In contrast, when exposing BSA encapsulated PLGA NPs, the BSA was found to be taken up into the cells when the charge of the NP surface changed from negative to positive with the increase of DOTAP. These results indicate that hydrophilic biomaterials can be delivered into plant cells by hydrophobizing BSA with DOTAP and encapsulating in PLGA NPs.

PB320 Flow polarization imaging around a rising bubble in aqueous gelatin solution under gelation process
Hiroya TOKUUMI1, Hiroki NAKANO1, Shuichi IWATA1, Ryo NAGUMO1, Hideki MORI1, Tsutomu TAKAHASHI2
1 Nagoya Institute of Technology, Nagoya, Aichi, Japan
2 Nagaoka University of Technology, Nagaoka, Japan

Flow polarization images around a rising bubble launched periodically in a cooling 3 wt% gelatin aqueous solution were captured by the two-dimensional high-speed polarization camera. The flow polarization imaging involves two dimensional retardation profile and orientation angle profile, respectively. The retardation profile correlates significantly with elasticity of the transparent solution under the validity of the optical-stress rule. The gel point of the aqueous solution of gelatin was evaluated in terms of viscoelastic properties of the gelling solution measured by the rheometer.
When the solution temperature was kept at 33 deg-C, shape of the rising bubble with diameter of 8.4 mm was sphere. Because the solution behaves Newtonian fluid of sol and the bubble Reynolds number was order of 1. There is almost no observable change for both retardation and orientation angle profiles.
After cooling down of the prepared mixture to 23 deg-C at a cooling rate of 0.254 deg-C/min for 39.3 minutes, and then the solution was kept at 23 deg-C. The growth of network structure of gelation causes the transition of the bubble shape from spherical to spherical with a cusp at about 3400 seconds. Slight change in retardation and clear orientation angle profile can be observed around the rising bubble. This suggests that at this event the elasticity of the gelling solution became dominant.
At about 3600 seconds of gelation, the bubble shape transformed to oblate shape with a long tail. The long tail was surrounded by high retardation region (i.e., high elasticity solution) in vertical orientation much higher than the others mentioned before. Therefore, the solution surrounding the long tail, which was supplied by the negative wake, was stretched vertically.

PB321 Experimental study of the dynamics of two-fluid displacement and flow through eccentric annuli
University of Adelaide, Adelaide, South Australia, Australia

We have developed a flow system for systematic study of the dynamics of two-fluid displacement and flow in annuli. This flow is relevant to many industrial applications, especially in the cementing operations in oil/gas well completions, where drilling fluids are displaced from the annulus between the casing and the well bore by a series of spacer fluids and cement slurries. The annular flow apparatus was designed such that it can be operated at various degrees of eccentricity and different angles of inclination to simulate the type of flows in oilfield drilling and cementing operations. A special feature of the flow system is that the inner cylinder can be rotated during displacement, allowing the effects of helical flow on the performance of the annular displacement process to be studied. Flow and displacement tests with various model non-Newtonian fluids representing drilling and cementing fluids have been conducted at different eccentricities, pipe inclinations, and over a range of flow rates and cylinder rotational speeds. Flow visualization results consist of video images capturing the movement and evolution of the fluid-fluid interface that clearly illustrate the effects of annular eccentricity, displacement rate, cylinder rotation and fluid rheology on the displacement process. In addition, the velocity profiles through the annulus at a fixed point downstream were determined by Particle Image Velocimetry (PIV) method to elucidate the dynamics of the displacement fluid-fluid interface. PIV results obtained have revealed interesting phenomena in the two-fluid displacement region that depend on the fluid viscosity ratio and the annular eccentricity. Pressure drop measurements produce data that can be used to develop models for predicting friction losses associated with single-fluid flow and two-fluid displacement in eccentric annuli. The results obtained will be useful for assessing the effects of important variables such as eccentricity, pipe inclination, flow rate, pipe rotation and fluid rheology, that impact on the effectiveness of the annular displacement process.

PB322 Continuous production of monodisperse water-in-oil droplets under pulsatile-free flow
Jia Yun XU, Jen Han CHIU, Jen-Huang HUANG
National Tsing Hua University, Hsinchu, Taiwan

Peristaltic pump is a most commonly used equipment to continuously introduce the fluids. However, the peristaltic pump may produce a pulsatile flow, limiting the applications for multiple purpose. To address this issue, we design a flow stabilizer connected to the outlet of the peristaltic pump. When the pulsatile flow enters the flow stabilizer, the elastic membrane is expanded, providing a buffer area for the pulsatile flow to deliver a steady flow at the outlet. The dampener was fabricated with PET and acrylic embedded with polyurethane as elastic membrane. The results show that the flow stabilizer with larger overall membrane area has better Fluctuation Reduced (FR) value. The elastic membrane arranged in parallel shows better FR than in series. Furthermore, we established a droplet generating system, which can continuously form uniform water-in-oil droplets. The system included a droplet generator, a flow stabilizer, and pumps. To continuously generate the droplets, a peristaltic pump connected with a flow stabilizer were used to feed the continuous phase. We used soybean oil with 5% w/v Polyglycerol polyricinoleate (PGPR) as the continuous phase, while pure water as the discontinuous phase. The results demonstrate that the monodisperse water-in-oil droplets with uniform size could be successfully generated. Moreover, the size of the droplets could be well controlled by altering the flow rate of the water. The droplet diameter increased from 0.311 mm to 0.519 mm (CV < 5.7 %) as the water flow rate increased from 5 μL/min to 80 μL/min. In conclusion, the integration of the flow stabilizer reduced the pulsatile flow of the peristaltic pump up to 99.25 %, allowing to continuously generate monodisperse droplets. We envision our continuous droplet generator system can apply to continuous and scalable production of chemical.

PB323 Photocatalytic degradation of volatile organic compounds (VOCs) using ultrasonically-formed droplets containing TiO2 particles
Miyu KIMURA1, Ryoichi NAKAYAMA1, Norikazu NAMIKI1, Kazuhiko SEKIGUCHI2, Naoki KAGI3
1 Kogakuin University,
2 Saitama University,
3 Tokyo Institute of Technology

It has been well-known that photocatalytic degradation of indoor air pollutants such as volatile organic compounds (VOCs) requires a high degradation efficiency for a long time without lowering the efficiency due to coking on UV-irradiated TiO2 surface.
In our previous work, we fabricated a photocatalytic reactor where photocatalytic degradation of toluene vapor took place on TiO2-containing droplets formed by ultrasonic atomization under UV irradiation, and demonstrated that a stable degradation efficiency was kept continuously under the optimal operation conditions. Therefore, the present study aimed to improve the degradation efficiency by dispensing additive particles, such as carbon black (CB) and silica (SiO2) ones, into TiO2 suspension.
Two sorts of additives were added to the suspension of TiO2 particles in the photocatalytic reactor, and then the toluene generated at a concentration of 5 ppm was fed to the reactor where TiO2-contaninung droplets irradiating with UV-A ray. The test air was sampled at the inlet and outlet of the photocatalytic reactor, and the concentration of toluene was determined by a GC-FID to determine the decomposition efficiency Ed.
As a result, Fig1 shows the change in the toluene degradation efficiency with additive particle concentrations in the TiO2 suspension. It was found that the degradation efficiency was increased by dispensing the additives, and had an optimum value at the additive concentration. The addition of CB might increases the adsorption capacity to toluene while it might be considered that CB coats the surface of TiO2 particles and UV-ray does not reach TiO2 when the amount of addition exceeds a certain level. Although the effective surface area of TiO2 particles contributing to the photocatalysis increased due to the presence of suspended TiO2 agglomerates coated with SiO2 particles, the extreme increase in SiO2 particles could enhance the ultraviolet light shielding effect of the SiO2 particles themselves.

PB324 Numerical Simulations using the MOF method of Individual Bubbles Growing and Detaching due to the Nucleate Boiling Process.
Yuya SHUDO1, Mitsuhiro OHTA1, Mark SUSSMAN2
1 Tokushima University, Tokushima 770-8506, Japan
2 Florida State University, Tallahassee, FL 32306, USA

In this study, the growth of single bubbles, due to nucleate boiling, is computationally examined using the moment-of-fluid (MOF) method. The MOF method represents sharp deforming boundaries between phases by evolving both volume fractions and centroids. The numerical simulations are carried out by solving the conservation law equations for mass, momentum, and energy, taking into account the appropriate jump conditions at the interfaces between substrate, liquid, and vapor. The substrate is modeled as a rigid body and the liquid and vapor are simulated as incompressible fluids. In order to carry out the computations as efficiently as possible, dynamic adaptive mesh refinement (AMR) is used in which the computational domain is an adaptive hierarchy of rectangular grids as opposed to a single, uniform, rectangular grid. In this study, the single bubbles periodically growing from a single nucleation site is two-dimensionally and axisymmetrically computed under the temperature difference of 7 K between the heated wall and water. The continuous generation process of nucleate boiling bubbles is successfully reproduced and thermal and flow fields in nucleate boiling are clearly shown. It is showed that the second and subsequent bubbles after the first bubble is released are stably formed. The growth rate of the second and subsequent bubble is faster than that for the first bubble and the growth rate of the second and subsequent bubble becomes constant. It is verify that our computational results for the growth rate of bubbles, the bubble departure time and the bubble volume for the second and subsequent forming bubbles agree well with the previous study. Our computational method using the MOF method can be effective for computing boiling problem.

PB325 An Experimental Study of the Morphology of a Drop Rising through a Hydrophobically Modified Alkali-Soluble Emulsion Polymer Solution
Takuya MAKITA1, Mitsuhiro OHTA1, Shuichi IWATA2
1 Tokushima University, Tokushima 770-8506, Japan,
2 Nagoya Institute of Technology, Nagoya 466-8555, Japan

The dynamic motion of a drop rising in a highly viscous hydrophobically modified alkali-soluble emulsion polymer (HASE) solution with pH = 9.0 is experimentally examined. In this study, a 1.65 wt% HASE polymer solution, which is adjusted to pH ≈ 7.0 by adding a sodium hydroxide solution, is used as the matrix fluid and salad oil is used as the drop. For sufficiently small drops, a threadlike long trailing edge at the rear of a drop is observed, similarly as for the case of a small bubble. It is newly observed that the threadlike trailing edge for the drop is much longer and thicker than that for the bubble. From experimental observations, it is found that a drop is more susceptible to the effect of viscoelastic properties of the surrounding HASE solution than the for the bubble case. For the case of bubble rise motion, the threadlike long trailing edge turned into threadlike branches for increasing bubble size; no such the tendency can be observed for the drop case. In microscopic observations for the bubble, it was reported that very unique microstructures emerged from gas-liquid interfaces in the motion of large bubbles, but on the other hand, microstructures at the liquid-liquid interface cannot be observed for drops. It is revealed that there is a pronounced difference between how the viscoelastic properties of a HASE solution effects bubble morphology versus drop morphology.

PB326 Numerical simulation of droplet breakup in droplet swarm by level set method with multi-levels
Tsinghua University, Beijing, China

In the present work, a level set method with multi-levels was developed to simulate the breakup behaviors of multi-droplets in a macro extraction device, the pulsed disc and doughnut column (PDDC). The deformation and breakage behaviors of droplets in droplet swarm were accurately captured. It was found that the shearing effect of the continuous phase flow, the collision between droplets and internals, and the Plateau–Rayleigh instability of the interface of liquid column would lead to the breakup of droplets. These breakup patterns were in good agreements with the experimental results. Moreover, the influences of pulsation intensity and droplet size in droplet breakup behaviors were investigated. It was concluded that the higher the pulsation intensity, the easier the droplets to break up. In addition, large droplets were found to be more likely to break up into multiple daughter droplets. Generally, the droplet breakup frequencies predicted by numerical simulation presented acceptable agreements with experimental results.

PB327 Determination of dynamic interfacial tension in a pulsed disc and doughnut column
Bo WANG, Shan JING, Shaowei LI*
Institute of Nuclear and New Energy Technology, Tsinghua University

Interfacial tension is an essential physical property in liquid-liquid extraction. Under mass transfer conditions, the interfacial tension changes due to mass transfer and concentration variation. The measurement of dynamic interfacial tension in such condition is of vital importance to direct the design of extraction equipment. In previous study (Zhou at al., Chemical Engineering Science, 197, 172–183), we presented the quantitative relation between the droplet breakup frequency function and interfacial tension. In this work, the droplet breakup frequency function in mass transfer process is measured at four different heights of the pulsed column using the method developed in our previous work (Zhou at al., AIChE Journal, 63(9), 4188-4200). The dynamic interfacial tension at different height of the column is then calculated by regression method. With acetic acid – water as the dispersed phase and 5%TBP – kerosene as the continuous phase, the results show that the dynamic interfacial tension increases with height decreasing. Comparing the dynamic interfacial tension with the static interfacial tension of equilibrium liquid-liquid system, we have found that the former is smaller than the latter when the concentration of acetic acid in the continuous or dispersed phase is the same. This result indicates that interphase mass transfer leads to decreasing of the interfacial tension. The decreasing extent of the dynamic interfacial tension has a positive correlation with the mass transfer flux. Furthermore, under mass transfer conditions, binary breakage occupy the dominant position at different height and up to quintuple breakage can be observed.

PB328 Volume penalization method for gas-liquid-solid flows with immersed free surface
Giang NGUYEN1, Kimiaki WASHINO1, Takuya TSUJI1, Toshitsugu TANAKA1
1 Osaka University, Suita, Osaka, Japan

Particulate flows in a mixture of gas and liquid, i.e. gas-liquid-solid three-phase flows are encountered in many powder handling processes, and it is of paramount importance to understand the complex interactions among phases. This work presents a simulation model for the study of gas-liquid-solid flows. In the given model, the solid-fluid interaction force is introduced by the Volume penalization (VP) method inherited by the study of Engles et al. (2015). Considering solid phase as a porous media, fluid is allowed to flow through the solid body with a certain value of permeability. The solution for determining the local permeability at solid surface is described in this study. The interface between gas and liquid is captured by the Volume of fluid (VOF) method coupled with the Level Set (LS) method. The interfacial interactions, such as the wetting of solid surface and surface tension force, are taken into account by using the Immersed free surface (IFS) model proposed by Fujita et al. (2015). The IFS model is originally designed for LS method, and it is adapted to the VOF framework in this work. Several simulations are performed to validate the proposed model. The resultant drag force acting on the particle as well as the computed surface tension force and contact angle are in good agreement with the theoretical values.

PB329 Molecular simulation of Amphotericin B-Ergosterol channel for the novel water channel
Hao Chen WU1, Tomohisa YOSHIOKA1, Keizo NAKAGAWA1, Takuji SHINTANI1, Toshinori TSURU2, Hideto MATSUYAMA1
1 Kobe University, Kobe, Japan
2 Hiroshima University, Hiroshima, Japan

Global water scarcity and pollution became a serious problem around the world. In order to relieve this issue, membranes processes have been widely applied to water treatment. With the increasing requirement for water, it is indispensable to continue search for membranes that can perform more efficiency on water flux and ion rejection. Although novel Aquaporin channel synchronously transports water rapidly while completely rejecting ions [1], Aquaporin is very costly and has a complex fragile molecular structure. These drawbacks make it formidable to apply or produce on an industrial scale. Some materials were searched for replacing Aquaporin channel, an Amphotericin B-Ergosterol channel (AmBEr) is a type of self-assembly organic channel [2], AmBEr is considered to be potential candidates for the biomimetic artificial channels. In this study, AmBEr was constructed and explored via molecular simulations to understand the structural characteristics and transport mechanism on a molecular scale. The chemical structures of Amphotericin B and Ergosterol monomer were shown in Fig. 1. AmBEr molecular models were constructed and investigated via molecular dynamics (MD) simulation procedures, Monte Carlo (MC) simulations and Quasi-non equilibrium MD were adopted to simulate the properties of calculation models. Structural characteristics and water sorption behaviors were introduced. All these molecular models were constructed using BIOVIA Materials Studio® commercial software.

PB330 Effects of Filling on the Segregation Patterns in a High Fill Level Rotating Drum
Chi-Yun LIN1, An-Ni HUANG1,2, Hsiu-Po KUO1,2
1 Department of Chemical and Materials Engineering, Chang Gung University, Taoyuan City 33302, Taiwan
2 Department of Otolaryngology-Head & Neck Surgery, Linkou Chang Gung Memorial Hospital, Taoyuan City 33305, Taiwan

When we mix particles of different sizes in a rotating drum, segregation may be observed. In this study, the segregation of a granular mixture consisting of 50% v/v 0.55 mm glass particles and 0.77 mm glass particles in a high fill level rotating drum is investigated. The segregation patterns appear after few revolutions, and the surface segregation band size and number become stable within 30 min. For the spatio-temporal diagram, the two smaller-particle-rich segregation bands take longer time to move toward the drum center when the fill level increases from 60%, 70% to 80%. The segregation patterns inside the bed are observed by the cross-sectional slices. Peculiar segregation patterns composed by strikes and core-in-core regions are observed for the first time.

PB331 Validation of CPFD models in for simulating hydrodynamics, heat transfer and reactions in CFB downer reactor for coal pyrolysis with binary particles
Yingya WU, Xiaogang SHI, Xingying LAN*, Jinsen GAO
State Key Laboratory of Heavy Oil Processing, China University of Petroleum, Beijing

It is well known that the low-rank coals such as lignite and sub-bituminous with higher volatile and moisture accounts for a large proportion of coal, which is the most abundant fossil energy resource in the world.Compared to the conventional fixed and Circulating Fluidized Bed (CFB) riser, the downer reactor, in which gas and solids move downwards concurrently, provides prominent advantages such as more uniform gas-solids flow structures, shorter particles residence time, less gas-solids back-mixing.These unique characteristics is particularly beneficial and suitable for coal fast pyrolysis, which need short but uniform contact time to achieve high reaction selectivity and yields of desirable products.
The CPFD approach based on the MP-PIC method was used to study the hydrodynamics, heat transfer and reactions in CFB downer reactor for coal pyrolysis with binary particles. Firstly, the comparison of solids distribution profiles of binary particles showed the CPFD approach with the cluster-based drag model is capable of accurately predicting the gas-solids flow and particle mixing characteristics. On the basis of the suitable gas-solids flow model, we further simulate the heat characteristic by coupling the heat transfer model. The simulation results show the temperature mixing index predicted by the coupled heat transfer model is in good agreement with the experimental value. By analyzing the heat transfer characteristics of particles, it can be found that the heat transfer process mainly occurs near the entrance, the interaction between gas and particles near the nozzle plays a key role. Therefore, more attention needs to be paid on the inlet of particles design of the downer. Finally, the CPFD model was further extended to reveal the coupled flow-reaction behaviors on coal pyrolysis reactor with binary particles. The results demonstrate that the CPFD model can accurately predict the hydrodynamics, heat transfer and reactions behavior of binary particles.

PB332 CFD simulations and optimization for transfer lines
Research Institute of Petroleum Processing, SINOPEC, Beijing 100083, China

Transfer lines are used for petroleum transportation in petrochemical industry, which refers to a pipeline connection between the furnace and the vacuum column in the crude oil distillation unit. A well designed transfer line can increase the evaporation rate of the feed in the vacuum distillation tower, improve oil product quality and save energy for vacuum distillation unit. To better design new pipelines or optimize the existing ones, it is necessary to describe the hydrodynamic behaviors of the multi component fluids flow in the pipelines. However, due to the complexity of the system, the model for analysis of the flow inside the transfer line has been paid less attention.
In this study, a commercial computational fluid dynamics package was used to develop a three-dimensional, fully turbulent flow model for the simulations of the multi component mixture across a pipeline. The simulated results of the pressure, velocity and Ma number distributions provided an insight into the phenomena occurring within the flow field in the transfer line. The analysis showed that the pressure drop of a transfer line is strongly affected by the structure of transition section and junction. On the basis of this, a new structure is presented to optimize the flow, and the pressure drop is lower than that of the original one. The modeling methodology and the results presented will be useful for evolving better designs of transfer line or other similar pipelines.

PB333 Influence of Hansen Solubility Parameter on Resin Flowability
Yuki KATO1, Shinichi TSUTSUMI2, Nobuyuki FUJIWARA1, Hideki YAMAMOTO1
1 Kansai University, Osaka, Japan
2 Panasonic Corporation, Osaka, Japan

The bonded magnet is a composite material formed by mixing a resin and magnetic powder and is used as a magnet for a motor or a sensor. In the case of a bonded magnet by injection molding, it can be molded into various shapes, the degree of freedom in motor design is increased. Furthermore, if the magnetic properties of the injection molded bonded magnet can be enhanced, a more sophisticated motor can be manufactured. To enhance the magnetic properties, it is necessary to increase the concentration of the magnetic powder. However, as the addition amount of the magnetic powder is increased, aggregation of the magnetic powder occurs, the flowability decreases, and there is a problem that molding becomes difficult. To address this problem, measures have been taken to treat the surface of the magnetic powder with a silane coupling agent or the like to improve the affinity to the resin and to disperse the resin highly. However, at present the selection of surface treatment agents is conducted empirically.
Therefore, I focused on Hansen Solubility Parameter (HSP). HSP is a physical property value quantitatively evaluated the affinity between substances, and it is applied in various fields.
In this study, we examined the influence of HSP of magnetic powder on the flowability after mixing. First, the magnetic powder was surface-treated with various silane coupling agents, and HSP was measured to evaluate its affinity to PA6. Thereafter, the magnetic powder was added to PA6 at the same concentration, and the viscosity was measured to measure the fluidity. As a result, as the combination of the magnetic powder and the resin with better affinity in HSP theory, the flowability was improved, and it was possible to confirm a high correlation between the HSP and the flowability.

PB334 Characterization of the Particle Cohesion by the Pressure Drop Monitoring in a Fluidized Bed
Wan-Yi HSU1, An-Ni HUANG1,2, Xuan WANG1, Hsiu-Po KUO1,2
1 Department of Chemical and Materials Engineering, Chang Gung University, Taoyuan City 33302, Taiwan
2 Department of Otolaryngology-Head & Neck Surgery, Linkou Chang Gung Memorial Hospital, Taoyuan City 33305, Taiwan

Evaluation of the flow properties of powders is essential for the designing and operations in powder-related industrial processes. In this study, the powder cohesion is evaluated by the measurement of the pressure drop across the bed during the fluidization and defluidization procedure in a FT4 Powder Rheometer (Freeman Technology). Glass beads with the diameter between 88 μm to 105 μm are fluidized in a cylindrical glass column with inner diameter of 50 mm by dried nitrogen. The pressure drop overshooting at incipient fluidization is used to characterize the powder cohesion and the bed cohesion is altered by Polyethylene glycol (PEG) addition and powder bed normal stress. When no normal stress is applied to the powder bed, the pressure drop shows a linear increasing during the overshooting period. However, when a 20 kPa normal stress is applied to the powder bed, the pressure drop shows a two-stage increasing during the overshooting period. A new cohesive powder fixed bed fluidized bed transition state is initially reported. The measured bed cohesion is ca. 45 Pa, which is comparable to the theoretical prediction of Molerus (1975) and Rumpf (1975).

PB335 Numerical study on fine particles behavior in an incinerator by coarse-graining DEM
Putri Mustika WIDARTININGSIH1, Yuki MORI1, Mikio SAKAI2
1 Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
2 Resilience Engineering Research Center, School of Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan

In industrial incinerators, hazardous substances from the combustion waste are concentrated in fine particles. Collecting and reducing fine particles is particularly challenging, and thus it becomes one of the major problems in industrial incinerators. A control plate was introduced into an industrial incinerator to reduce the emissions of hazardous fine particles. However, the efficiency of the control plate against smaller sized particles is necessary to be investigated. The coupled discrete element method and computation fluid dynamic (DEM-CFD) are employed in a solid-fluid coupling problem even where the solid particle size is fine. One of the concerns in the DEM-CFD is the requirement of excessive computer capacity when the numerous particles are applied. Therefore, a coarse grain model is adopted here to perform the calculations using a single computer. For the sake of approaching the actual condition of industrial incinerators, poly-dispersed particle system is employed in this study, where expressed by a combination of two coarse ratio. The simulation results show that the control plate works well to capture the fine particles. In addition, the velocity of particle flow is significantly decreasing when a control plate is employed.

PB336 Establishment of a novel universal index to predict particle packing ability for various slurries
Nana KITAI1, Naoya IWATA1, Takamasa MORI2,3, Hiroshi SATONE3,4
1 Graduate School of Science and Engineering, Hosei University, Koganei, Tokyo, Japan
2 Department of Chemical Science and Technology, Faculty of Bioscience and Applied Chemistry, Hosei University, Koganei, Tokyo, Japan
3 Hosei University Research Institute for Slurry Engineering, Koganei, Tokyo, Japan
4 Department of Mechanical and System Engineering, Graduate School of Engineering, University of Hyogo, Himeji, Hyogo, Japan

In wet-shaping process, it is well known that the particle dispersion state greatly affects the product qualities. Thus, it is important to evaluate and control the slurry characteristics properly. However, quantitative evaluation index that can compare the particles dispersion state under different slurry conditions such as particle concentration, primary particle size, particle density, shape, and solvent viscosity, etc. has not been established, yet. In this study, we measured the hydrostatic pressure of slurry prepared in various conditions under gravitational settling. Additionally, the time change of hydrostatic pressure for the slurry whose particles are assumed to perfectly disperse as primary particle was calculated for the slurries used in this study. In order to calculate the hydrostatic pressure curve more precisely, the particle size distribution measured by sedimentation was considered. The novel index, the ratio of measured hydrostatic pressure to calculated one, was introduced to discuss the particles dispersion state. It was shown that the index (hydrostatic pressure ratio) had a good connection to the final packing fraction of the sediment, indicating that the index should be useful to quantitatively represent the particle dispersion state regardless of the slurry conditions.

PB337 Alteration of lipid bilayer permeability and structure of liposomes using gas-liquid flow
Makoto YOSHIMOTO, Tetsuya FUJIE, Yusuke NAKANO
Department of Applied Chemistry, Yamaguchi University

Phospholipid vesicles (liposomes) are soft colloidal particles which are applicable to fabricating carriers of biopharmaceuticals, protocell models, and immobilized enzymes. The permeability of lipid bilayer membranes towards small molecular mass compounds is of significance for the above applications of liposomes. In the present work, 5(6)-carboxyfluorescein-containing large unilamellar liposomes were prepared and suspended in gas-liquid flow generated in external loop airlift bubble columns. The release rate of the dye molecules from the liposomes was examined at 40 °C under various operation conditions of the bubble columns with different configurations. The permeability coefficient of the dye molecules through the bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine lipids was clearly dependent on the superficial gas velocity of the bubble column (1-3 cm/s). The mean diameter and size distribution of liposomes were practically unaffected by being suspended in the bubble column, as evaluated with the dynamic light scattering measurements. These results indicate that the permeability of liposomes increased in the gas-liquid flow. On the other hand, the release rate of the dye molecules through the lipid bilayers incorporated with small fractions (1-5 mol%) of poly(ethylene glycol)-conjugated lipids exhibited significantly higher than that through unmodified lipid bilayers. The size distribution of the modified liposomes was altered by being suspended in the bubble column, indicating that the membrane integrity of liposomes was significantly altered. The effect of liquid shear stress on the structural stability of liposomes was separately examined by shearing the liposomes in a cone-and-plate geometry to assess the bubbles-liposomes interaction in the bubble column. All of the above results demonstrate that the structural stability of liposomes was controllable using the gas-liquid flow.

PB338 Effects of relaxation time of solutions on vortex deformation and turbulence statistics in two-dimensional turbulent flow
Kobe University, Kobe, Hyogo, Japan

A small amount of polymer in turbulent flow reduces turbulent drag, which is called drag reduction. It is said that extensional rheological properties, such as long relaxation time and extensional viscosity of the fluid are keys for the phenomena. The characteristic extensional rheological properties of the solution deform vortices in turbulent flow, thus, the turbulent flow is also affected.
In this study, in order to focus on the effects of extensional rheological properties of fluids on vortex deformation in a turbulent flow, a flowing soap film known as a two-dimensional (2D) flow with a small influence of shear stress is used. The vortex shedding in the flow were visualized using interference patterns, and the local velocity fields were measured by particle image velocimetry (PIV). The vortices shed at the grid in the 2D flow was affected by polyethylene oxide (PEO, Mw = 3,500,000) concentration. The original shape of the vortices was deformed by increasing the PEO concentration, and the vortices were disappeared close to the grid at a certain concentration. However, the vortices again appeared at even higher concentration. In this study, we categorized these vortices into three types. Indeed, the generation of these vortices were affected by the extensional relaxation time of the PEO solution. The difference of these vortices were also analyzed by the fluctuation intensity of the flow field and by the turbulent statistics. Turbulent kinetic energy transfer, production and dissipation term of the turbulent kinetic energy were calculated the local velocity. The turbulent statistics implies variation of energy transfer in 2D flows, which was caused by the viscoelasticity of PEO solution.

PB339 Study on the fixed salt removal effect by fine bubble
Yuji MIKASA1, Naoya YAMAWAKI1, Hiroaki SHIDA2, Hayato OKUMURA3, Shigenori AKAMATSU3, Yusuke NISHIUCHI3, Takashi HATA3
1 Advanced Course of Material Engineering, National Institute of Technology (Kosen), Kochi College, Nankoku, Kochi, Japan
2 Ligaric Co., Ltd, Suita, Osaka, Japan
3 Department of Social Design Engineering, National Institute of Technology (Kosen), Kochi College, Nankoku, Kochi, Japan

In recent years, it has become clear that fine bubble (micro/ultrafine bubble) have properties different from those of ordinary bubble. For example, the surface area of fine bubbles per the same volume is larger than that of normal bubbles, including chemical reactions and physical adsorption at the gas-liquid interface, a dramatic improvement in mass transport, and having a positive/negative charge (zeta potential) on the surface of fine bubble. Cleaning using the features of fine bubble in particular has been noticed, and there is also a high industrial need. For example, the use of plastic cleaning in recycling industries and anti-corrosion cleaning of structures with flying sea salt particles in coastal areas has begun.
Therefore, in this study, the effect of fixed salt removal by fine bubble was investigated. As a result, it was confirmed that the cleaning effect of fine bubble water was higher than that of control water, and the effect was dependent on the concentration of fine bubble. It was also confirmed that the cleaning effect was higher in the cleaning of ultrafine bubble water than in the control water. The removal mechanism of the fixed salt seems to be that gas in the fine bubble dissolves by the pressure change around the fixed salt, resulting in supersaturation, and the bubble nucleus is generated, and a part of the salt is lifted by the micro bubble, which promotes the fixed salt exfoliation. And, it was confirmed that the cleaning effect was higher by the constant introduction of fine bubble than the system in which fine bubble existed only in the cleaning initial stage (Fig). From these results, the usefulness of fine bubble in cleaning can be expected.

PB340 Synthesis of ε-Fe2O3 particles using Sr2+ doped magnetite-silica core-shell particles
Magnetic Powder Metallurgy Research Center, The National Institute of Advanced Industrial Science and Technology (AIST), 2266-98 Anagahora, Shimoshidami, Moriyama, Nagoya, Aichi, 463-8560, Japan

The magnetite nanoparticles including strontium were synthesized by hydrothermal method. After these particles were dispersed in the ethanol aqueous solution, and coated with the several nano meter silica shell by adding TEOS solution. The ε Fe2O3 particles were prepared by heat-treating these core-shell particles by the fixed temperature. Silica was removed from these core-shell particles by use of the NaOH solution, and the particles only of ε Fe2O3 were made. The magnetic property and the special quality of the form of particle and the particle diameter were argued to the dope amount of the strontium, the silica covering amount and the synthesis condition of the heat treatment temperature to this particle. As a result, the crystal system by heat-treatment time and the change in form of particle were checked by XRD and TEM.
Samples with 5% added Sr2+ were examined for several holding times set at 900 °C. Figure was XRD at 900 °C. The sample held for 1 hour has an aspect ratio of 1, and mostly contains ε phase and a little γ phase. Furthermore, in the sample which passed for 6 hours, although it was mainly ε phase, a little α phase was generated. The aspect ratio increased to 2 or more, and the rod width also increased to 60 nm. The crystal of silica developed with the increase of the amount of Sr2+. It is considered that this is because Sr2+ infiltrates into the silica by the doping of Sr2+, and influences the crystallization. It was found that the addition of Sr2+ causes the formation of an ε phase at a low temperature, the extension of the particles in a rod shape, and the change of the magnetic properties.

PB341 Improvement of mixing performance for highly viscous fluid using asymmetrical anchor blades
Marina YAHATA1, Takafumi HORIE1, Yoshiyuki KOMODA1, Naoto OHMURA1
1 Kobe University, Kobe, Hyogo, Japan

In chemical process industries, research and development on liquid mixing is important to save energy and resources, improve productivity and develop new materials. An appropriate impeller shape is selected to optimize mixing state in a chemical process. Very high viscosity of a fluid is an obstacle for performance of mixing processes. An anchor blade impeller is frequently used for mixing of highly viscous fluid is effective for heat transfer by scraping the tank wall but poor for the mixing. This is because the axial flow is not enough to be developed. It has been reported that mixing of highly viscous fluid can be improved by imposing unsteadiness on the flow field. This study, therefore, proposed two types of asymmetric anchor blade to improve agitation performance. Power consumption and flow patterns were experimentally and numerically observed. The asymmetric anchor blades proposed in this study were half anchor and G-shaped anchor. It is shown in Figure 1. As the results obtained by numerical simulation using a commercial CFD code (R-Flow), the power consumption using the asymmetrical anchor blades becomes smaller than that of the regular anchor blade. However, no significant difference between the half anchor and G-shaped anchor could be found. As the results of the flow pattern observation, it has been found that the mixing performance is improved in asymmetric anchor blades as compared with the regular anchor blade. It can be considered that asymmetric anchor blades make the flow field unsteady. Furthermore, it has been found that using the G-shaped anchor shows better mixing performance than the half anchor. The axial flow can be generated by the protruding portion of the G-shaped anchor. Mixing time was measured by decolorization experiments using iodine and sodium thiosulfate. The G-shaped anchor progressed the decoloring reaction most quickly and the mixing time the shortest.

PB342 Synthesis of red light emitting Si nano-particles in liquid phase
Keiichiro ISHIHARA1, Yoshiko TSUJI1,2
1 Department of Chemical System Engineering, the University of Tokyo, Tokyo, 113-0033, Japan
2 Environmental Science Center, the University of Tokyo, Tokyo, 113-0033, Japan

Semiconductor nanoparticles, known as quantum dots, are biocompatible semiconductor fluorophores and highly stable against exposure of light, which cannot be expected for organic fluorescent materials.1 Controllability of luminescence energy is also advantageous for bio-imaging since bio-permeability is good in near infrared to visible red region. However, toxic metals frequently used in quantum dots are problematic for practical use. In this context, non-toxic silicon-based nanoparticles meet the demand of bio-imaging although its stability cause the difficulty in the liquid phase synthesis, which realize large scale fabrication.
In this paper we fabricated silicon nanoparticles with two different liquid-based approaches, “Top-down” and “Bottom-up” methods. In “Top-down” method, Si powders are etched by hydrofluoric-nitric acid and decreased in its particle size to nm scale in suspension to obtain silicon nanoparticles, which show red light emitting (about 620 nm) in its PL spectra. In “Bottom-up” method, 3-Aminopropyltrimethoxysilane was used as precursors and reduced by reducing agent to obtain Si nanoparticles. The PL spectra of the nanoparticles show blue light emitting (about 420 nm), indicating that the obtained nanoparticles have smaller diameter than the nanoparticles prepared by “Top-down” method.

PB343 Optimal conditions of paddle shape and rotational speed for flow homogenization inside honeycomb
Hiroki UCHIYAMA, Takehiro ESAKI, Yousuke MATSUKUMA
Fukuoka University, Fukuoka, Japan

The dependence of the paddle rotating speed and its optimal paddle length for flow homogenization inside honeycomb has been investigated experimentally. Three types of honeycomb having various square hole size and the number of holes was used. The paddle length was varied within the range of 62, 81.6 and 120 mm and its width was 22 mm. The paddle was made of polylactic acid, which was fabricated using a 3D printer. The inlet volumetric flow rate was constant, 200 L/min. The air flowed into the header section and then into the honeycomb. The paddle attached to the motor was set to the windward side of the honeycomb. The flow rate was measured with a hot-film anemometer, which was positioned using a traverse device. The uniformity of the velocity at the honeycomb outlet was evaluated using the standard deviation of the velocity. The standard deviation is the arithmetic mean value calculated with the 169-point velocity data.
Table1 shows the experimental results. It was found that the optimum condition was almost 500 rpm of rotational speed with 120 mm paddle diameter.

PB344 Vibro-fluidized bed separation of particulate construction waste residue based on density-segregation
Tetsuaki MATSUOKA1, Takamasa YOKOUCHI1, Jun OSHITANI1, Yoshihide MAWATARI2
1 Okayama University of Science, Okayama, Japan
2 Kyushu Institute of Technology, Kitakyushu, Japan

Construction waste residue (CWR) consists of particulate organic and inorganic materials like woodchips and gravels. In Japan, the CWR should be discarded at a controlled landfill site to avoid decay of the organic materials. However, illegal dumping of the CWR is a serious trouble, because the discard at a controlled landfill site is “expensive”. If the organic materials are separated to reduce the amount to be < 5 wt%, the CWR can be discarded at an “inexpensive” least-controlled landfill site (≈ 1/5 $ of controlled landfill site). The organic materials can be easily removed from the CWR using a wet separation. Woodchips float and gravels sink in water. However, the wet separation has disadvantages; waste water treatment is necessary to avoid water pollution, a drying process is required after separation, and leakage of water may occur from the separator. Therefore, the development of dry separation is in great demand as a substitute for the wet separation. In this study, we focused on “density-segregation” in a vibro-fluidized bed. If a particulate mixture having different density is fluidized by an airflow through the bottom, lighter particles move up and heavier particles go down; this phenomenon is called “density-segregation”. If vibration is added to the fluidized bed, the density-segregation becomes clearer. We applied the vibro-fluidization to reduce the amount of organic materials in the CWR. We found that the upper layers of the fluidized bed consist of the organic materials, and the amount of organic materials in the lower layers is reduced to be < 5 wt%.

PB345 Coating of ibuprofen crystal with products from hydrothermal hydrolysis of polysaccharides
Chuo University, 1-13-27 Kasuga, Bunkyo-ku, Tokyo 112-8551, Japan

In the pharmaceutical industry, many of drugs are poorly soluble in water [1]. Ibuprofen, a nonsteroidal anti-inflammatory drug [2], is widely used to reduce fever and to treat pain or inflammation but it is poorly soluble in water. Coating with sugar is commonly used to enhance water solubility of drugs and expected to improve solubility, mask bitter taste of drugs and preserve medical effects. Due to relatively high values of the ionic product of H+ and OH- [3], mono- and oligosaccharides can be produced from the hydrothermal hydrolysis of polysaccharides [4]. The objective of this study is to demonstrate coating of less soluble drug crystal with products derived from hydrothermal hydrolysis of polysaccharides: ibuprofen was employed as a drug species, and polygalacturonic acid as a polysaccharide component. Mono- and oligo-galacturonic acids were hydrolytically obtained in a semi-batch flow reactor. The hydrolysis product solution played as a role of an anti-solvent for ibuprofen crystallization. The morphology and surface of ibuprofen crystal was observed with a scanning electron microscopy (SEM). Figure 1 (a) and (b) shows the SEM images of the neat ibuprofen crystal and coated ibuprofen crystal, respectively. The crystal structure was analyzed with the X-ray diffraction, and melting temperature of the ibuprofen crystals by differential scanning calorimetry. The effects of operating conditions on the amount of coating layer on the ibuprofen crystal surface were investigated.
[1] E. Merisko-Liversidge et al., Eur. J. Pharm. Sci., 2003, 18, 113-120
[2] A. Mauricio et al., Fluid Phase Equilib., 2013, 354, 185-190
[3] T. Funazukuri et al., J. Chem. Technol. Biotechnol., 2004, 79, 229-233
[4] T. Miyazawa et al., Ind. Eng. Chem. Res., 2004, 43, 2310-2314

PB346 Generation, identification and characterization of ultrafine bubbles in water
Takashi GOSHIMA, Susumu NII
Kagoshima University, Kagoshima City, Kagoshima, Japan

Fine bubble (FB) attracts keen attention and the appilcation is spreading in various fields. FB is classified into microbubbles (MB) and ultrafine bubbles (UFB) depending on the size. Although benefits of UFB have been recognized in many applications, there are many obstacles for deeper understanding of the effect due to the lack of methods to characterize UFB. This study aims to introduce series of UFB generators invented by the author and propose methods to identify and characterize UFB.
For identification, a method to differentiate bubbles and impurities was proposed. The dominant matter in water, either bubbles or impurities was successfully differentiated by analyzing signals obtained from a laser-based instrument, Nanosight LM-10. Comparison of the distribution patterns obtained from mixtures of UFB and various impurities (Fig. 1) suggested that UFBs are stabilized by adsorbing impurities in water. For characterization, the effect of salt addition on the number density and size distribution of UFBs was examined. Addition of NaCl solution into UFB water gave rise to UFBs' expansion and agglomeration depending on the condition. The dissolved gas concentration in NaCl solution and UFB determines the dominant change to occur. When the gas concentration in NaCl solution concentration is higher than that in UFB water, majority of the UFBs expand while the rest of UFBs agglomerate, vice versa. Furthermore, the net surface tension of single UFB in water was experimentally determined by using the invented system to estimate the volume of gas dissolved and/or dispersed in liquid. The net surface tension was lower than the surface tension of pure water.
These results provide supporting information for a suggested mechanism for stability of UFB in water by Yasui (2016) and the report of decreased surface tension of UFB water by Terasaka (2018).

PB347 Experimental study on jet asymmetry of two-component layer impinging stream mixer
Jianwei ZHANG, Yuhang YAN, Fanrong MA, Yifan ZHANG, Xinli SHA
School of Energy and Power Engineering, Shenyang University of Chemical Technology, Shenyang 110142, China

The flow field energy and mixing characteristics of the two-component layer jet asymmetric impingement mixer were studied by two-dimensional high-speed particle image testing technology (TR-PIV) and laser induced fluorescence (PLIF). In this study, the POD analysis method is used to extract the large-scale structure of the flow field and decompose the instantaneous flow field to obtain the main modal energy. The flow field energy is mainly concentrated in the first-order mode, which is easy to describe the flow characteristics of the whole flow field. Flow field reconstruction can accurately reduce the distribution of the initial flow field vortex structure. The asymmetric flow field energy formed by the up and down and left and right offsets of the stagnation point is always larger than the symmetrical flow field, and the mixing rate is lower than the symmetrical flow field. Under the same working condition, the stagnation cross flow field energy is higher than the stagnation point on the same side. The mixing efficiency is opposite; the lower tracer (stagnation point) is left (fully developed) to facilitate flow field mixing (Figure 1). Therefore, this asymmetry of the jet is beneficial to increase the flow field energy, and the initial development of the tracer promotes flow field mixing.

PB348 Modeling and Numerical Study of Ceramic Paste Extrusion
Boris GOLMAN1, Piotr Skrzypacz1, Wittaya Julklang2
1 Nazarbayev University, Nur-Sultan, Kazakhstan,
2 Suranaree University of Technology, Nakhon-Ratchasima, Thailand

The extrusion process of ceramic pastes is commonly used for the production of high-value products such as catalyst pellets for the chemical reactor, honeycomb catalyst for purifying gas exhausted from an automobile, ceramic packings for adsorption and direct heat transfer, etc. Recently, extrusion-based 3D printing was successfully applied for manufacturing of ceramic superconductors, zeolite monoliths, and porous scaffolds.
Ceramic paste extrusion is a complex process which depends on the paste rheological properties, die and extruder geometries, and operational parameters. The paste rheological properties are controlled by several factors including volume fraction of particles and their size distribution, shape, packing density and surface characteristics, as well as amount and properties of binder and other additives. Modeling and quantitative analysis of paste molding are important to design the extrusion process for the production of high-value extrudates of desired strength, shape and morphology.
The objective of this paper is to construct the mathematical model of ram extrusion of ceramic materials and use it for the optimization of the extrusion process. The paste extrusion was modeled as a flow of non-Newtonian fluid. The paste continuity and momentum equations were solved numerically for fluid based on the modified Herschel-Bulkley viscous model. Fig. 1 illustrates the velocity distribution in the extruder. The effects of materials properties, extruder and die geometries, temperature and inlet paste velocity on the distributions of paste velocity and pressure in the extruder and die were investigated numerically.
1. M. Li, L.Tang, R.G. Landers, M.C.Leu, Extrusion process modeling for aqueous-based ceramic pastes, J Manuf Sci E-T ASME, 135, 051008, 2013.
2. N.Labsi, Y.K.Benkahla, A. Boutra, Hydrodynamic and thermal characterization of the flow of a Herschel-Bulkley fluid in a pipe, ESCAPE20, 1-6, 2010.

PB349 Numerical simulation of gas-solids flow in bubbling beds with different internal structures
Jiang YUAN, SHI Xiaogang, LAN Xingying, GAO Jinsen
China University of Petroleum, Beijing, China

Bubbling beds are widely used in process industry, e.g. stripper of Fluid Catalytic Cracking (FCC) in petrochemical industry. It is important to enhance gas-solids contacting by improving the gas-solid flow, bubble distribution, bubble diameter and gas-solid backmixing in the stripper by adding internal structures. Different types of internal structures have been reported including disk and donut baffles and grid baffles. In this study, an Eulerian–Eulerian approach was used to compare the gas–solids flow behaviors in strippers with disk and donut and grid packing and in the empty stripper. Simulation results with respect to the solid volume fraction agrees well with experiments, indicating that the model is accurate in predicting gas-solids flow behaviors in bubbling beds. The CFD model also clearly predicts the local dead zones underneath the disc and donut baffles. The gas-solids flow behaviors are more uniform in the grid stripper than in other types of strippers. A machine-learning based method of capturing and characterizing the bubble was proposed. It was found that any strippers with internal structures can obtain smaller bubbles than the empty bubbling bed. The number of bubbles in the grid stripper is about 10 times more, and the bubble size is much smaller than those in the empty-type stripper. The present work indicates that the internal structures can enhance bubble breakup. The performance of the bubbling bed can further be improved by optimizing the internal structures.

PB350 Development of new synthesis method for hollow polymer particles
Yusuke KAWAI, Tetsuya YAMAMOTO
Graduate School of Materials Design Innovation Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi, Japan

Polymer hollow particles have many advantages: great light diffusivity, high heat blocking rate, low density, low cost, and good dispersion in resin. However, the conventional synthesis methods have problems with their uniformity and productivity, which prevent application for products. In this research, we developed a new method to synthesize a large amount of uniform hollow polymer particles. Polymer solid particles were synthesized by soap-free emulsion polymerization using styrene as a monomer and 2,2'-Azobis [2- (2-imidazolin-2yl) propane] dihydrochloride (VA-044) as a polymerization initiator. In the colloid of the polymer solid particles prepared above, methyl methacrylate (hydrophilic monomer) divinylbenzene (crosslinker), and ethyl acetate (good solvent for polystyrene) were added and the polymer solid particles were swollen with the monomers and the solvent. Then, VA-044 was added and the polymerization was started. Through drying the fine particles, polymer hollow particles were obtained. A TEM image of the polymer hollow particles synthesized by this method is shown in following. Monodispersed hollow particles were synthesized (Cv: 9.2%, Cv = particle diameter dispersion / particle diameter distribution × 100). In addition, it was shown that the hollow polymer particles have different optical characteristics and heat blocking properties.

PB351 Effect of geometrical configuration of baffles on mixing performance in a stirred vessel
Eri SATO1, Hiroo HORIFUCHI2, Takafumi HORIE1, Yoshiyuki KOMODA1, Naoto OHMURA1
1 Kobe University, Kobe, Hyogo, Japan
2 Sumitomo Heavy Industries Process Equipment Co., Ltd., Saijyo, Ehime, Japan

Mixing is a general process operation and plays an important role in many chemical process industries. It is used for various purposes such as gas-liquid, liquid-liquid, and dispersion of solid-liquid. In agitation in a stirred tank, a vortex is formed in a turbulent flow. When the vortex is formed, the circumferential rotational flow becomes dominant and makes the vertical mixing weaken. Consequently, performance of mixing decreases. Installing baffles can break the vortex and the flow is complicated, so that the performance of mixing is improved. Thus, the baffles play a very important role in fluid mixing. However, there are very few researches on the effect of geometrical configuration of baffles on mixing performance. Therefore, in this study, baffle length, width, clearance and the number of sheets were varied to investigate the agitating power consumption and mixing characteristics. A transparent acrylic tank with a semi-elliptical bottom was used as a stirring tank, and MAXBLEND© was used as an impeller. Five geometrical configurations of baffles with different length and width were used. Rotational torque was measured using tap water as the liquid to determine the power number, and the power diagram was obtained by plotting the power number as a function of the impeller Reynolds number. Decolorization experiments of iodine and sodium thiosulfate were conducted to determine mixing time. It is generally known that the power number in the stirred tank with baffles is independent of Reynolds number in fully developed turbulent flows. This study, therefore, tried to correlate the power number with the geometrical parameters, and a relatively good non-dimensional correlation equation with accuracy within ± 10% was successfully obtained in fully turbulent flow regions.

PB352 Formation Mechanism of Dead Water Region and Its Influence on Solute Transport in Irregular-Shaped Micro Channels
Kizuki TOYAMA, Fumina TOGI, Tetsumi KUBOTA, Shusaku HARADA
Hokkaido University, Sapporo, Japan

It is important to understand solute transport in micro-channels in various fields of engineering. For example, solute transport in natural flow channels, such as rock cracks or gaps in soils, is closely related to resource extraction or pollutant dispersion in underground. In a micro-channel, a laminar dispersion, which is known as Taylor dispersion, plays a significant role on the solute transport. However, most of the natural channels have an irregular shape and the solute dispersion is greatly influenced by the irregularity of the channel such as variance of the cross-sectional shape and the roughness of the channel wall. This is because dead water region (stagnant region) is formed in irregular-shaped channels.
The purpose of our study is to understand the formation mechanism of dead water region and its influence on solute transport and to make clear the dispersion process in various-shaped micro-channels. We examined the effect of the irregularity of wall surface on solute transport experimentally. We quantified the solute concentration field in micro-channels by image analysis based on absorption photometry and also measured the flow velocity distribution in the channels by particle image velocimetry technique. From the measurement results, we quantified the dead water region formed near the channel sidewall.
From the experimental results using rough-walled channels with various heights and wavelengths, it is found that dead water region decreases the effective aperture length and the solute transport is enhanced in flow direction. We developed a model on the formation of dead water region as a function of the height of wall roughness in micro-channels. The model describes well the formation mechanism of dead water region and the relationship between dead water region and solute dispersion in various-shaped channels.

PB353 Measurement of microbial adhesive forces using parallel plate flow chamber and atomic force microscopy
Takumi MATSUMOTO, Akinori YOSHIHARA, Yasuhiro KONISHI, Toshiyuki NOMURA
Osaka Prefecture University, Sakai, Japan

The vast majority of environmental microorganisms adhere to every interface and form highly structured microbial aggregates called biofilms. Biofilms are structured communities of microbial cells enclosed in a self-produced extracellular polymeric substance (EPS). The presence of the EPS provides up to 1,000-fold greater resistant to antimicrobial agents when compared with their planktonic counterparts. The adhered cells usually cause piping corrosion in various types of facilities. Removal of adhered cells on various surfaces is an important issue that improves productivity and ensures safety in industrial processes. In this study, the effects of the physicochemical properties of microbial cell surfaces on cell detachment from solid surfaces were investigated using a parallel plate flow chamber (PPFC) and an atomic force microscopy (AFM). Gram-positive Lactococcus lactis was used as model microbial cell, and negatively charged glass slide was used as model substrate. The cumulative number of detached cells from the model substrate as a function of flow rate in PPFC was fit with the Weibull distribution function, and the drag force distribution was calculated from the hydrodynamic force required to remove attached cells. Based on the DLVO theory, because the microbial cells and the glass substrate are negatively charged, it can be predicted that the interaction force between the microbial cells and the glass increases with an increase of the ionic strength. However, interestingly, the drag force to remove the attached cells on the surface shifted toward smaller values with increasing ionic strength. AFM imaging of microbial cells attached on the glass surface was carried out with a stepwise increase of the contact pressure. The required removal forces estimated form the lateral forces in AFM deceased with an increase of ionic strength. This was found to be agreed well with the result in PPFC.

PB354 Experiment and CFD simulation of flow characteristics in stairs model
Fukuoka University, Fukuoka City, Fukuoka, Japan

Wind flows through stairs cause various problems, such as subway trains wind, and chimney effects in stairs of high-rise buildings, but the flow structure in stairs has not been fully understood yet. In the present study, scale model of a 4-storey stairs and a 4-storey slope model having same inclination were fabricated using by a 3D printer. This model experiment was performed on the assumption that windows of the landing stairs were closed. Each model inlet was connected with a wind tunnel outlet in order to obtain uniform velocity distribution condition. The flow structures of the stairs model and the slope model were compared. Velocity was measured with the hot wire anemometer. The flow velocity was calculated with the OpenFOAM. The observed velocity was agreed with the CFD values. Therefore, the k-e model is applicable for flows of the stairs and slope model. Flow structures of the stairs were found to be similar to that of the slope.

PB355 CFD-PBM simulation of drop size distributions in a pump-mixer: Influence of the local turbulent dissipation
Han ZHOU, Shan JING, Shaowei LI
Institute of Nuclear and New Energy Technology, Tsinghua University, Beijing 100084, China

Coupling the computational fluid dynamics (CFD) and population balance model (PBM) is a promising approach to simulate the liquid-liquid two-phase flow. In this study, the influence of the local turbulent dissipation on drop size distribution (DSD) in a pump-mixer is studied using the CFD-PBM approach. The breakup kernels obtained in our previous study (Zhou et al., 2019. Measurement of droplet breakage in a pump-mixer. Chem. Eng. Sci. 195, 23–38.) were adopted and the class method is selected to solve the population balance equation (PBE). The drop size distributions (DSDs) at five different rotating speeds were simulated to investigate the influence of the local turbulent dissipation. The prediction of the DSDs was compared with the experimental results and show a good agreement, which verified the correctness of the breakup kernels.

PB356 Cluster characteristics in the jet mixing zone of fast fluidized bed
Zihan YAN, Yiping FAN, Chunxi LU, Chunming XU
China University of Petroleum-Beijing

As an efficient equipment for gas-solid reaction, mass and heat transfer, fast fluidized bed is widely used in modern chemical industry. For example, catalytic cracking, fluidized coking, coal gasification and combustion, and the gas phase polymerization. In these processes, high-speed nozzle jets are usually employed to realize a full mixing and reaction between gas and solids. The dynamic contact of feed jets with fluidized particles will directly affect the reaction results. Thus, it is important to know the cluster behaviors in this mixing process. By combining the wavelet decomposition theory with particle concentration signals, method for identifying clusters in the jet mixing zone of fast fluidized bed was proposed. Characteristics of clusters during the dynamic mixing process were obtained. Comparisons between upward and downward jets were made by analyzing the cluster time fraction. Results show that the characteristics of clusters in the upward and downward jet mixing zones are quite different. Overall, the time fractions of clusters are much higher for the upward jets, especially in the initial contact region of jets with particles. For the whole jet mixing zone, when the jets are upward, the average cluster time fraction is about 8.91%, while for the downward jets, the average cluster time fraction is only 6.19%. Besides, it is shown that the cluster time fraction is quite high near the wall of fluidized bed in the upward jets mixing zone, which results in the low efficiency of mass and heat transfer between gas and solids. In the industrial production, some undesirable phenomenon occurs. For example, the coking problem in FCC riser. On the contrary, if the jets mixing with the fluidized particles countercurrently, the time fraction of clusters decreases significantly. The gas-solid mass and heat transfer efficiency is expected to be improved.

PB357 Lagrangian analysis of solid particle concentration distribution in solid-liquid stirred tank with various particle concentration and impeller height
Masato KURODA, Shunnosuke IMAI, Ryuta MISUMI, Meguru KAMINOYAMA
Yokohama National University, Yokohama, Kanagawa, Japan

Solid-liquid mixing is widely used in crystallizer and chemical reactor with catalysis. In these processes, there is a problem that solid particles do not disperse uniformly and deposit on the bottom of the tank. So it is necessary to quantify the solid particle concentration distribution in the tank and investigate the conditions for improving particle suspension in order to solve these problems. However, it is difficult to measure the detailed particle concentration in the tank under many conditions. So, in this study, Euler-Lagrangian simulations were performed with various particle diameter, particle density, particle number, impeller speed, and impeller height to quantify the vertical distribution of solid particle concentration in the vessel and identified the factor that determines the dispersion state.
The Euler-Lagrangian approach, which tracks all particle motion, was adopted. Turbulent flow in the vessel was represented by large eddy simulation. The interactions of particle–particle and particle–solid surfaces were modeled using the distinct element method. The fluid was presumed to be water. The standard deviation of particle concentration in the vertical direction, σc, was calculated and time averaged σc, σc,av, was used as the evaluation index of dispersion state in the vessel.
Results show that particles suspension was promoted and σc,av decreased as the particle number increased because particles collisions increased around vessel bottom. And σc,av was correlated with terminal settling velocity (determined by particle size and particle density) and impeller tip velocity (determined by impeller speed and impeller diameter). Finally, when the flow pattern changed greatly by changing the impeller height in the liquid depth, H, the particle dispersion state with impeller height of H/3 and H/9 were promoted comparing to that with impeller height of H/5.

PB358 Evaluation of mixing performance in Porous Hollow Fiber Membrane device
Department of Chemical Engineering, Osaka Prefecture University

Nanoparticles are used in biomedical, optical, and electronic fields. In general, it is important for the synthesis of nanoparticles to bring about highly supersaturated reaction solution conditions. However, batch type reactions often result in wide size distributions of resultant particles and low reproducibility, because the non-uniform supersaturation condition inside the reactor cannot be avoided due to its poor mixing performance. In this study, therefore, we proposed usage of porous hollow fiber membrane (PHFM) device as a new approach of mixing process for the improvement of mixing property. PHFM devices are mainly used in the waste water treatment process. On the other hand, PHFM device have possibility to provide a more efficient mixing condition than the conventional mixing equipment, because they have large number of small pores in the wall surface of the hollow fiber. However, the mixing performance in PHFM devices have not been evaluated quantitatively. The objective of this study is characterization of mixing performance in PHFM devices by the Villermaux-Dushman method and understanding mixing mechanism in PHFM devices using the computational fluid dynamics (CFD) simulations.
Fig. 1 illustrates the champion data of the mixing time using different four mixing devices (PHFM device, batch-type mixer, T-shaped mixer, and microreactor). The PHFM devise shows much fast mixing time of approximately 0.3 ms as well as the microreactor. Moreover, PHFM devices achieved the high mixing performance despite a low Re number. This would be because the large number of small pores provide infinite water jet streams. Furthermore, the CFD simulations of the mixing process were conducted to investigate the mixing performance and the effect of pore size and feed flow rate on the mixing performance.

PB359 Experimental measurement of particle sedimentation in modified raceway ponds
Ryunosuke OSAKA1, Yoshiaki UEDA2, Yusuke SAKAI3
1 Graduate Student, Setsunan University, Neyagawa, Osaka, Japan
2 Setsunan University, Neyagawa, Osaka, Japan
3 Kumagai Gumi Co.Ltd., Shinjyuku, Tokyo, Japan

Microalgae have been paid attention for biofuels production in the bioenergy sector. Recently, an open-channel raceway pond system is commonly used from the viewpoint of the production rate and construction/operation cost. In the open-channel raceway ponds, a dead zone around the corner of a central baffle inhibits the mixing and, therefore, the dead zone volume is required to decrease in volume to increase the growth rate of microalgae. In this study, we experimentally seek to optimal geometry of a raceway pond for efficient production rate of biofuel energy. In practice, the tracer particles having the specific density of 1.02 are included in a newly proposed raceway ponds, and the sedimentation of particles are recorded every one minute with the use of a digital camera. In addition, the flow patterns in the raceway ponds are measured with the particle image velocimetry (PIV) and computational fluid dynamics (CFD) simulation. Based on the results obtained in the experiments, the physical model to design an optimal layout of the baffle in a raceway pond are proposed successfully.

PB360 Formation and Evaluation of Au/ZnO Particles by Spray Pyrolysis Method
Young-Jun LEE1, Toshiyuki FUJIMOTO2, Shinya YAMANAKA2, Yoshikazu KUGA3
1 Division of Engineering, Muroran Institute of Technology, Muroran japan
2 College of Environmental Technology, Muroran Institute of Technology, Muroran japan
3 Muroran Institute of Technology, Muroran japan

In order to expand the absorption wavelength range of zinc oxide utilized as a photocatalytic catalyst, studies of Au nanoparticles depositing on zinc oxide (Au/ZnO) have been actively carried out. Since the Au nanoparticle can absorb visible light, Au/ZnO is expected to have catalytic properties for absorbing UV-ray and visible light. The purpose of this study is to generate Au/ZnO particles by an ultrasonic spray pyrolysis method (USP).
A dilute aqueous solution of Zn(NO3)2 and HAuCl4 was used as a raw material solution. Small droplets of raw material solution generated by ultrasonic nebulizer were thermally decomposed in the reaction tube. The generated Au/ZnO particles were characterized by SEM, XRD, TEM and UV-Vis equipment. In order to prevent particle aggregation, Trion X-100, a non-ionic surfactant, was added, and its effects were also examined.
The particle size was found to be several micrometer by SEM observation. ZnO and metallic Au crystalline were found by the XRD analysis. As shown in lower part of Figure, the color of obtained Au/ZnO particles was from thick purple to pale purple. As shown in the upper part of the Figure, the TEM images of the particles showed that the Au crystals were well dispersed inside the ZnO particles. Crystallite size estimated by Scherrer's equation was agreed with that measured from TEM images: the cristallite size decreased as decreasing the additive rate of HAuCl4. The obtained particles were measured by UV-Vis equipment. The smaller the size of the Au particles, the higher the visible absorption peak near 400 nm. This is considered to be the result of resonance with light having higher energy as the atomic vibration caused by surface plasmons increases as the Au particle size is smaller. In the case of adding Trion X-100, the absorption peak shifted toward the longer wavelength side: aggregation of Au particles was enhanced by the addition of Triton X-100.

PB361 Langevin Dynamics Calculation of Brownian Coagulation Coefficient for Spherical Equal-size Aerosol Particles in Transient Regime
Toshiyuki FUJIMOTO1, Shinya YAMANAKA1, Yoshikazu KUGA2
1 College of Environmental Technology, Muroran Institute of Technology, Muroran, Japan
2 Muroran Institute of Technology, Muroran, Japan

Coagulation coefficient of aerosol particles due to Brownian motion is an important issue in aerosol science and technology. Motion of aerosol particles is diffusive in continuous region (small diffusive Knudsen number, KnD), or like free molecular motion of gaseous molecular in free molecular region (large KnD). Fuchs (1964) presented an expression of coagulation coefficient in transition regime by so-called "Flux Matching" method: physical model of this method was very clear and have been well accepted. In his method, the particles transport inside of the "limiting sphere" was assumed to be like free molecular, or diffusive outside of the sphere.
These days, some researchers presented coagulation coefficient of aerosol particles by direct calculation of motion of aerosol particles. They employed Langevin dynamics equation to represent the stochastic motion of aerosol particles. In this study, we developed new model to calculate the coagulation coefficient. As shown in Figure, our model employed spherical space in which one scavenging particles are in the center of it: the spherical space moves together with the motion of the scavenging particle. Other colliding particles move in the spherical space until the particles collide to the scavenging particle (1). The collided particle is re-generated at the surface of the spherical space (2). The particle escaped from the spherical space (3) is re-entered (4). The coagulation coefficient can be calculated from the mean time interval of collision and the concentration of colliding particle.
By using the above numerical model, we have calculated the coagulation coefficient of spherical particles of from 4 nm to 400 nm in diameter. Our results were slightly greater than that of Fuchs. The difference was significant considering the uncertainty caused by Monte Carlo Method. The results were organized by using KnD, then we presented an expression of non-dimensional diffusion coefficient.

PB362 Optimum mixing ratio of sewage sludge and food waste in cylingrical agitated vessel for composting
Young Sei LEE, Kiyoung LEE
Kyungpook National University, Sangju City, Gyeongbuk, Koera

During the past decade, the proportion of food waste in municipal waste stream is gradually increasing and hence a proper food waste management strategy needs to be devised to ensure its eco-friendly and sustainable disposal in the world. The most well-kwon approach to treat food waste has been considered compositing, which is breaking down biodegradable waste by naturally occurring micro-organisms with oxygen, in an enclosed vessel or tunnel or pit. In order to achieve optimal compositing, the mixing ratio of food waste and micro-organisms in sewage sludge and the reaction behaviors are quite important.
In this work, we studied about reaction rate, salinity, carbon/nitrogen ratio, reaction temperature, and organic substance during compositing process in the mixture of food waste and sewage sludge to find optimal composition. For our study, the sewage sludge and food wastes were collected from a rural area and garbage trucks operated by Sangju City.
Compositing experimental was carried out with various mixing ratio of the food wastes and sewage sludge such as 10:30, 30:70, 50:50, 60:40, 70:30 and 90:10 wt.%, respectively.
Finally, we reveal that the optimum ratio for compositing was 60 : 40 wt.%. Moreover, the detailed results show that the temperature was 18 ~ 22 oC at the initial stage, then it was sharply increased to 44 ~ 46 oC up to 1day, which was maintained over the fermentation period. In the case of a pH change depending on fermentation process, it was decreased to 4.5 ~ 6.2 due to production of organic acid at the initial stage of reaction. Later on, ammonia was produced that leads to increasing pH value of 7. 3~ 8.0.

PB363 Constant drying rate of wet material containing solvent with low-boiling point
Shizuoka University, Hamamatsu, Shizuoka, Japan

A constant drying rate is important to evaluate the performance of dryer. The constant drying rate is difficult to measure when a surface evaporation rate is extremely high like the material containing solvent with a low-boiling point. To measure the constant drying rate, the material surface must be wet during measuring. In this study, to keep the material surface wet, liquid (solvent) was supplied to the material (a sintered-glass ball filter was used as the material). The validity of the method was confirmed for water as the supply liquid. Ethanol and acetone were used as the solvents with the low-boiling point. The constant drying rate and the temperature of material were measured at steady state. The temperature of material obtained from the experiment was compared with the results calculated from the conventional method (Chilton-Colburn analogy).
Table 1 shows the results of constant drying rate, temperature of material and heat transfer coefficient on material surface. In the case of water, calculated temperature of material (wet-bulb temperature) is 297.8 K which is almost equal to that in experiment (297.2 K). The difference of the material temperature between the calculation results (Chilton-Colburn analogy) and experimental data was large when the material surface was assumed to be covered with the pure solvent (ethanol). It implies that the water vapor in air affects the evaporation of solvent (ethanol) on the material surface. The gas-liquid equilibrium (water-ethanol system) on the material surface is considered to predict the temperature of material, and the results are in good agreement with the experimental data (in Table 1). A heat transfer coefficient can be calculated from the experimental data (drying rate and temperature of material). The heat transfer coefficients obtained from the experiment are in orderly agreement with the values predicted (Ranz-Marshall equation).

PB364 Flow pattern analysis of self-induced sloshing due to bubble flow in a rectangular vessel
Ryohei AOKI1, Satoko FUJIOKA2, Koichi TERASAKA2
1 School of Science for Open and Environmental Systems, Graduate School of Science and Technology, Keio University
2 Department of Applied Chemistry, Faculty of Science and Technology, Keio University

Self-induced sloshing is oscillatory phenomenon of a liquid free surface caused by a flow of a fluid. This phenomenon is reported in some gas-liquid reactors, and it is important to predict and prevent its occurrence from the viewpoint of safe operation of the reactors. However, the mechanism of this phenomenon has not been clarified.
Therefore, in this study, self-induced sloshing caused by bubble flow in a rectangular vessel was experimentally examined in order to elucidate its generation mechanism. In our experiment, air bubbles were dispersed into tap water from the bottom of the rectangular vessel.
First, the effects of gas flow rate and liquid height on the oscillation of liquid surface were investigated by image photographing. Oscillation frequency and amplitude of liquid surface were analyzed. It was confirmed that self-induced sloshing occurred by increasing the gas flow rate at specific initial liquid height. The gas flow rate threshold varied depending on the initial liquid height. Pendulous oscillation of bubble flow was observed simultaneously with free surface oscillation. The frequency decreased with the increase in liquid height.
Second, PIV was carried out. Figure shows unique circulating flow observed only when self-induced sloshing occurred. Liquid circulating flow occurring on both sides of the bubble flow transitioned within one cycle of self-induced sloshing. It is suggested that the circulating flow contributes to the oscillation of bubble flow and liquid free surface.
CFD modeling was also considered for scale-up of the vessel. It is difficult to reproduce the self-induced sloshing caused by bubble flow because it is a multiscale phenomenon that includes both a small-scale bubble interface and a large-scale free surface. We tried to reproduce oscillation of bubble flow and liquid free surface simultaneously using Multi-Fluid VOF model of Eulerian multiphase flow model in ANSYS FLUENT.

PB365 Observation of pseudo plume behavior by hydrate sedimentary layer decomposition
Yusuke TAKAHASHI, Ryosuke EZURE, Shun TAKANO, Hiroyuki KOMATSU, Kazuaki YAMAGIWA, Hideo TAJIMA
Niigata University, Niigata, Japan

We are focusing on the practical use of methane hydrate, which is expected in the future as alternative energy. The surficial and shallow methane hydrate presented under the seafloor generates bubble groups (that is “plume”) by decomposition. For recovery and use as energy resource, it is necessary to consider the possibility of clogging in the recovery pipe due to rehydration of bubbles. The purpose of this research is to observe experimentally and evaluate theoretically the decomposition behavior of hydrate sedimentary layer and the rising behavior of bubbles generated by hydrate decomposition. Chlorodifluoromethane was used as a low pressure model gas of methane. Hydrate sedimentary layer was produced by cooling and pressurizing water in countercurrent contact with gas using a hydrate formation recovery device. The recovered hydrate was decomposed by the heating or depressurization method, without flowing water. Since the water was sufficiently saturated with the gas component during decomposition, the bubble radius reduction by gas dissolution was negligible. Two theoretical rising velocities were derived from the theoretical value with using the Navier-Stokes equation or the values in consideration of the bubble shape and hydrate film existence. The experimental rising velocities of small spherical bubbles radius agreed well with the theoretical value by Navier-Stokes equation. The relatively large elliptical bubbles showed a behavior close to the theoretical value of bubble with hydrate film. The rising velocity may decrease for any reason. Further investigation is needed for the influence of the bubble shape term in the theoretical formula. Under the pressure and temperature conditions closer to the hydrate equilibrium line, almost no generated bubbles could be identified. This may be because of fine bubble generation

PB366 Effect of homogenizer dispersion on particle size and size distribution in continuous production of silica microcapsules
Naoki KURODA, Takafumi HORIE, Ruri HIDEMA, Keita TANIYA, Kosuke SUZUKI, Yoshiyuki KOMODA, Yuichi ICHIHASHI, Satoru NISHIYAMA, Naoto OHMURA, Hiroshi SUZUKI
Kobe University, Kobe, Hyogo, Japan

Low temperature heat less than 423 K exhausted from factories and power plants has been wasted due to the low usability. If it is properly transported and distributed to each homes, it can be utilized for the air heating and hot water supply. Latent heat transport by using phase change slurry is expected to be used for filling this space gap. However, the latent heat transport slurry still has problems such as plugging in the pipes and increase flow resistance. The problems can be solved by encapsulating the phase change material in silica microcapsules. In this study, we developed a continuous process enabling large production of silica microcapsules for application to a latent heat transport slurry. In the first stage of this process, W/O emulsion was produced with a homogenizer, and in the second stage, W/O/W emulsion was formed with hydrolysis reaction occurring at the oil-water interface. The correlation between the water droplet diameter of the W/O emulsion and the diameter of the silica microcapsules indicated that the W/O emulsion process using the homogenizer determined the capsule diameter. It was also observed that silica microcapsules formed when the water droplets escaped from the oil droplets. Furthermore, by considering Shinner's equation, it was found that the average capsule diameter can be controlled by the rotation speed of the homogenizer. Finally, the stirred tank of the liquid-liquid dispersion process with the homogenizer was changed to the small one in order to make the residence time became uniform. Monodispersed silica microcapsules could be obtained with this continuous process.
Acknowledgements: This study was supported by JST-Mirai-project: #JPMJMI17EK.

PB367 Drag coefficients of a rising ellipsoidal droplet in a stagnant liquid
Saitama University, Saitama City, Saitama, Japan

In order to design drop dispersion systems, such as flows inside solvent extraction columns, it is necessary to know the shape and motion of the dispersed droplets. In this study, numerical analysis of the motion of ellipsoidal droplets in a stagnant liquid was performed by using a three dimensional Front-Tracking method. Drag coefficients of the ellipsoidal droplets were determined from the numerical results by two different projected areas: one calculated based on spherical droplets and the other based on numerically obtained projected area. The drag coefficients based on the numerically obtained projected area were in good agreement with the ones predicted by the semi-empirical model for spherical droplets proposed by Myint et al (J. Fluid Sci. Technol., 1, 72, 2006). The drag model was derived by correcting the projected area, and successfully predicted drag coefficients for ellipsoidal droplets.

PB368 Particle impact energy variation with the size and number of particles in a planetary ball mill
Fumie HIROSAWA, Tomohiro IWASAKI, Masashi IWATA
Osaka Prefecture University, Osaka, Japan

Dry powder grinding is an important unit operation in many industries such as mining, fine chemical, food and pharmaceutical, ranging from coarse mineral ore to submicrometer-sized fine drug powder. Recently, ultrafine dry grinding processes with high energy milling, which can produce fine particles with enhanced performance and/or improved properties, have attracted much attention. In the ultrafine dry grinding, planetary ball mills have been employed, and the milling conditions can vary the particle breakage induced by the impact energy of particles, altering the particle properties, such as size distribution, specific surface area and microscopic structure. Therefore, the impact energy must be appropriately adjusted to produce the particles with controlled properties. In order to understand the impact energy under given milling conditions such as revolution speed of the mill pot, size and density of the grinding balls and ball-to-particle filling mass ratio (BPR), the behavior of particles and balls is often simulated using the discrete element method (DEM). Many studies on the analysis of planetary mills using the DEM simulation have been published, indicating that the impact energy depends strongly on not only the milling conditions but also the particle properties. In the ultrafine dry grinding processes, both the particle size and the number of particles can drastically change as the milling time elapses, resulting in temporal variation of the impact energy. However, the impact energy in planetary mills has not been sufficiently analyzed taking into account the change in particle size and number of particles. In this paper, in order to study its effect on the impact energy, the DEM simulation of a planetary mill was performed under various conditions of particle size and number of particles. The results showed that the impact energy in the normal and tangential directions was greatly affected by these factors, suggesting that the grinding mechanism changes during milling.

PB369 Power characteristics of closed impellers with different geometries of flat blades for agitation
Masaya NAKATA, Hironori MATSUNAGA, Masanori YOSHIDA
Muroran Institute of Technology

Agitation impellers for liquid-phase mixing in chemical processes were designed based on a conventional disk turbine impeller with six flat blades which sucks liquid across the planes including the upper and lower edges of the blades. The altered impeller has the structure for the inflow entrance to be fully confined. We call this “fully closed impeller” and the prototype “open impeller”. The fully closed impeller was modified for the blade geometries, namely, the number and the axial width of the blade. An impeller with the 12 flat blades of standard width (12-STD) and an impeller with 6 flat blades of small width (6-Small) provided narrower flow paths between the neighboring blades of the rotating impeller. Additionally, an impeller with 6 flat blades of tapered width (6-Tapered) provided the flow path with a constant cross-sectional area. These flow paths were expected to produce nearly uniform liquid flows there. Liquid flows in the impeller rotational region of a baffled vessel with the closed impellers were examined for the flow path between the neighboring blades of the rotating impeller and the internal flows were compared with that with the open impeller. Energy consideration in terms of the dimensionless parameters on the basis of the flow and power measurements demonstrated that the closed impellers were successful in the head characteristic rather than the flow one. The power transmission efficiency were enhanced notably for the 6-Tapered closed impeller because of a reduced viscous loss.

PB370 Simulation analysis of the effects of adhesive force distributions in a particle bed on improving discharge particle flowability
Doshisha University, Kyotanabe, Kyoto, Japan

One techniques for improving flowability is admixing smaller particles than main particles. In this technique, an adhesive force have a different value in all contact points between main particles. This is because a combinations of smaller particles coverage states changes in all contact points. This adhesive force distribution would improve a particle flowability. In this study, we examined effects of the force distribution on improving flowability by DEM simulation.
We used DEM simulation of a two-dimensional system. The particle diameter was 60 μm. Physical properties of the particle were mainly set based on glass material including TiO2 and BaO. Two kinds of Hamaker constants, 1.0×10-20 and 1.0×10-22, were used for changing an adhesive force. Adhesive force distributions were given to particle beds by following 2 types; 1) particle surface distribution by setting different adhesive forces at a contact point on the surface, 2) spatial distribution in a particle bed by arranging positions of different adhesive forces particles. The width of particle bed container was 600 μm and the container have a discharging hole of 400 μm at the bottom. The flowability was evaluated by a particle discharging rate. As a result, when we gave a particle surface adhesive force distribution, the flowability increases with increasing surface ratios of the smaller Hamaker constant. On the other hands, in spatial adhesive force distribution results, the flowability shows the maximum value at particle number ratio 50%, which means mostly the largest adhesive force spatial distribution. This is because a particle flow structure depends on the spatial distribution. From these results, we obtained a useful knowledge, which adhesive force spatial distribution has a large effect than the surface distribution for improving discharging flowability.

PB371 Coagulation of latex particles in the aqueous dispersion of epoxy resin
Takahashi TORU, Komoda YOSHIYUKI, Ohmura NAOTO
Kobe University, Kobe, Hyogo, Japan

The mixture of functional latex particles and epoxy resin matrix is known as one of polymer composites and can improve various functions of epoxy resin. Since the remaining of particle aggregates deteriorates the appearance of epoxy composite, the production process of easily dispersed latex aggregates is now strongly required. If aggregates of weakly attached latex particles containing epoxy resin is obtained, such an aggregate may contribute to realized prompt and uniform dispersion of latex particles in epoxy resin matrix. In the present study, therefore, the dispersion state of epoxy resin in water was controlled by dispersion conditions, and the effect of the dispersed state of epoxy resin on the salt coagulation of latex particles was studied. Aqueous dispersion of epoxy resin was prepared using homogenizer or high shear device, Filmix. The droplet prepared using homogenizer had the diameter of 0.1 - 50 μm and easily settled down. In contrast, Filmix could produce highly stable dispersion of epoxy resin having the size of 0.03 μm. And then, the change of latex aggregates in the salt coagulation process was studied based on mean size and morphology. In the epoxy free salt coagulation process, latex aggregates first grew up and then gradually destroyed as increasing mixing time. Additionally, it was observed that tightly packed latex aggregates were produced. Similar trend was obtained when added unstable epoxy dispersion prepared by homogenizer. On the other hand, latex aggregates produced with stable epoxy dispersion was loosely connected and became larger as increasing mixing time. It was probably caused by the following mechanism: Fine droplets of epoxy resin covered the surface of latex particles and the surface layer was contributed to form loosely connected aggregate and increase bonding strength between particles.

PB372 Electrophoretic classification based on differences in electrophoretic mobility caused by change in the applied electric field
Kosei ONO, Tomonori FUKASAWA, Toru ISHIGAMI, Kunihiro FUKUI
Department of Chemical Engineering, Graduate School Engineering, Hiroshima University, Hiroshima, Japan

We propose a novel electrophoretic wet classification method using the applied electric field strength as an operation parameter. The influence of the applied electric field, electrolyte concentration, and particle size on the electrophoretic mobility was evaluated via direct measurements by microscopic electrophoresis. For high electrolyte concentrations or small particles, the electrophoretic mobilities decreased with decreasing applied electric field, unlike the basic theory of electrokinetic phenomena. Then, we fabricated an electrophoretic classification device and investigated the separation efficiencies for various particle sizes, electrolyte concentrations, and applied electric fields. For low electrolyte concentrations, the change in separation efficiency with respect to the applied electric field depended on the particle size. Finally, we performed a classification experiment using the applied electric field as an operation parameter. For low electrolyte concentrations, the classification performance changed according to the applied electric field. The applied electric field that provided the highest classification performance was the same as that where the difference in the separation efficiency was the highest. These results highlight the possibility of electrophoretic wet classification using the applied electric field strength as an operation parameter.
Keywords: wet classification; electrophoresis; electrophoretic mobility; particle

PB373 Direct Measurement of Fluid Pressure in Gas-Solid Fluidized Bed by Capsule-type Sensor
Yuji HONDA1, Shiori SAITO1, Tetsuya ANZAI1, Shusaku HARADA1, Shunsuke KATO2, Jun OSHITANI2, Yuya SAKAMOTO3, Takuya TSUJI3, Hirokazu KAJIWARA4, Kei MATSUOKA5
1 Hokkaido University, Sapporo, Japan
2 Okayama University of Science, Okayama, Japan
3 Osaka University, Osaka, Japan
4 Ebara Environmental Plant Co., Ltd., Tokyo, Japan
5 Ebara Corporation, Tokyo, Japan

Fundamental understanding of the floating-sinking motion of solid objects in a gas-solid fluidized bed is important for various engineering processes such as a density separation technique. Basically the vertical motion of objects in a fluidized bed is explained by “hydrostatic effect”, i.e., the objects having the density larger than the apparent density of the fluidized bed sink into the bottom of the bed. However, recent studies have suggested that the floating-sinking motion of the object is influenced by various factors such as the object shape or the fluidization state of the bed.
The purpose of this study is to examine the fluid force (pressure force), which is the most significant force in the fluidized bed, acting on various-shaped objects and to investigate the shape dependency on the floating-sinking motion. We developed capsule-type pressure sensors and put them into the fluidized beds with various fluidization states. We measured the pressure difference between the upper and the lower surface of various-shaped objects. From the measurement results, it is found that the pressure around the object greatly varies with the object shape. The pressure field around the object in the fixed bed at low air velocity almost agree with the theoretical prediction from Brinkman equation which describes the permeation flow in porous media. On the other hand, the pressure in the fluidized bed at high air velocity largely fluctuates due to the interaction with rising bubbles, however, the average pressure difference between the upper and the lower part of the object reaches almost constant and does not depend on air velocity so much.
The results of this study suggest that the fluid force acting on an object in a gas-solid fluidized bed greatly depends on its shape and the floating-sinking motion is influenced by the fluidization state of fluidized bed.

PB374 Preparation of TiO2-x/C nanocomposite by milling of titania and polypropylene
Daisuke ORITA1, Nobutaka TAKEZAWA2, Kazuhiko SAEKI2, Taiji MATSUMOTO2, Mamoru SENNA3, Yasuyoshi SEKINE4, Kyoko OKUYAMA4, Masahide SATO1, Noboru SUZUKI1
1 Utsunomiya University, Utsunomiya, Japan
2 Industrial Technology Center of Tochigi Prefecture, Utsunomiya, Japan
3 Keio University, Yokohama, Japan
4 Nippon Coke & Engineering Co., Ltd., Tochigi, Japan

Lithium-ion batteries (LIB) are used in various fields, but there is a demand for higher performance to meet the needs for further safety, cost reduction, and long life. Although titanium dioxide is characterized by its chemical stability and high theoretical anode capacity, oxygen deficiency is known to enhance performance via increased conductivity and ion diffusion rate. Thus, the titania-based materials have attracted much attention as an alternative of graphite for the anode material. Although the electrode performance is improved by introducing oxygen deficiency into titanium dioxide, there is a drawback of higher cost such as heating up to 1500 °C under pressurized hydrogen atmosphere. In this study, we propose a new simple method obtaining such a modified titania via a mechanochemical route. Mixture of anatase nanoparticles and polypropylene were milled by a high performance media agitating mill (Nippon Coke & Engineering, AL01). The mechanochemical effects were recognized as a phase transformation from anatase to rutile and color change from white to gray, attributable to the introduction of oxygen deficiencies. These effects were more intensive with smaller amount of polypropylene. Subsequent heat treatment at temperatures as low as 500 °C in the atmosphere of 5 % hydrogen/argon brought about carbon-coated and oxygen deficient titania. Their electrical properties were examined mainly by the impedance spectra measurement.

PB375 Characterization of a technique to remove particles adhered on a surface using pulse air-jet blowings
Ami ISHII1, Ryoichi NAKAYAMA1, Norikazu NAMIKI1, Koya BABA2, Tomokatsu SATO2
1 Kogakuin University,
2 Techno Ryowa Ltd.

Demands for dry precise cleaning of various product surfaces contaminated with particles has been extended at various sorts of production processes such an assembly process of lithium ion batteries.Consecutive pulse-air jet (PAJ) blowings using a nozzle with a silt-shaped opening (rectangular nozzle) are effective for removal of micrometer-sized particles from various surfaces. In the present work, we aimed to make this dry surface-cleaning technique a more practical one for removing particles on separator film surfaces.
1.5 μm-sized PSL particles generated by a Colison atomizer were deposited to a silicon wafer used as a test surface placed in an inertial impactor. The test surface was exposed to PAJ blowings formed from the rectangular nozzle with the jet interval of 20 s and its duration of 1 s. Counting of the remaining particles was performed in the unit measurement area (0.25 × 0.33 mm2) using a microscope on the test surface at each pulse of blowing to determine the removal efficiency, fr_N. The planar distribution of the removal efficiency, fr_N on the test surface was obtained using the fr_N values of 143 points of the unit area in the staggered array of the measurement area (9.50 × 9.90 mm2).
Fig. 1 compares the planar distribution of particle removal efficiency between the rectangular nozzle and the circular nozzle with the same opening area. It was found that the rectangular nozzle successfully removed particles on a wider area at a higher efficiency compared to the circular one. The area mean removal efficiencies were over 90% for the rectangular nozzle, and 47% for the circular one.

PB376 Development of non-contact evaluation method to evaluate aggregation degree of nano particles using fluorescence spectrum analysis
Yusuke ANDO, Noriaki SANO
Kyoto University, Kyotodaigaku-Katsura, Nishikyo-Ku, Kyoto 615-8510 Japan

Nanoparticles are used in various industrial fields such as electronic devices, pharmaceuticals, and catalysts. The aggregation and dispersion of particles have a great influence on product quality in their manufacturing process. For example, production of multilayer ceramic capacitors includes a process of forming a thin film of the slurry and drying it. If there is a portion with a non-uniform degree of aggregation, it will cause product failure. So, understanding and control of the degree of aggregation of nanoparticles is important. In this study, we propose a new method to inspect the degree of particle aggregation using fluorescence spectrum analysis. Because this method is rapid, non-contact, and non-destructive, it can be a promising method of evaluating the degree of agglomeration of the nanoparticles on production lines.
We examined the influence of the aggregation degree of γ-alumina nanoparticle on fluorescence. We prepared two types of dry alumina particle aggregates with different aggregation states by different drying methods from a slurry of γ-alumina particles dispersed in water. By one, the slurry was dried at room temperature and, by another, the slurry was freeze-dried. From the difference in these drying processes, the particles dried at room temperature are more firmly agglomerated than freeze dried. The results of fluorescence spectrum measurement for these particle aggregates are shown in Figure 1. It can be seen that the fluorescence spectrum of the particles dried at room temperature is shifted to the longer wavelength side than freeze dried. The shift of fluorescence indicates the difference in the degree of aggregation from drying process. Namely, as the particles are aggregated, the wavelength of the fluorescence spectrum shifts to the longer wavelength side. This phenomenon would be considered to be similar to quantum size effect.

PB377 (canceled) <101531-1>
PB378 Investigation of the necessity of bottom punch movement in a rotary die-filling system
1 Department of Nuclear Engineering and Management, School of Engineering, The University of Tokyo, Tokyo, Japan
2 Resilience Research Engineering Center, School of Engineering, The University of Tokyo, Tokyo, Japan

A rotary die-filling system is widely employed to produce pharmaceutical tablets. In order to manufacture the tablet, the powder should be fed into a cavity uniformly. For the efficient feeding of powder, the punches are equipped under the die. As the punch starts moving downward, the powder inside the shoe is fed into the cavity. The exact depth of the punches is controlled to measure the filling amount of the powder inside the cavity. Although this system is widely used, its effectiveness has not been examined so far. In the current study, the powder feeding process is conducted using the Discrete Element Method coupled with Computational Fluid Dynamics (DEM-CFD) method. For the representation of the arbitrary shaped moving wall boundaries, the Signed Distance Functions (SDF) are employed, where the scalar field of minimal distance from the nearest wall is computed at the beginning of the calculation. Besides, for the evaluation of the interaction between gas and moving wall boundaries, the Immersed Boundary Method (IBM) is coupled with the SDF, in which the tracking of the interface between gas and wall boundaries can be precisely calculated in the staggered grids. Through comparing the results between with and without the moving punches, the mechanism is revealed for the first time. Consequently, the yield rate can be drastically improved using the numerical simulation techniques.
This study was financially supported by JSPS KAKENHI (17H02825), JSPS KAKENHI (17KK0110) and by Grant-in-Aid for JSPS Research Fellow from the Japan Society for the Promotion of Science (18J12267).

PB379 Starting Torque of Two Disk Turbines
Yukiya KUDO, Takanori KANEKO, Kazuhiko NISHI
Chiba Institute of Technology

When rotation of an impeller starts, strong torque (starting torque) is generated compared to the steady state in which the flow is fully developed. The starting torque is important for the strength design of the device. This study investigated the starting torque of two disk turbines within a mixing vessel, the relation between the starting torque and the distance of the two impellers, and the impeller rotational speed. Furthermore, CFD analysis was conducted to elucidate the starting torque generation mechanism.
Three starting torques, Ts1, Ts2, and Ts3, were observed as presented in Fig. 1. Ts1 is the torque generated during acceleration of the impeller. In this study we investigated Ts2 and Ts3. In Fig. 1, the relation between Ts2, Ts3 and the distance of two impellers (c) is shown. When c is small, Ts2 is small. Results show that Ts2 increases rapidly at c = 0.02 m and becomes almost constant from 0.035 m. Actually, Ts3 shows the same tendency, but Ts3 was not detected after 0.07 m. The following were found from the results of CFD analysis. Ts2 and Ts3 were generated when a flow is strong only around impellers and the flow of the whole mixing vessel is underdeveloped. When Ts2 was generated, two impellers produce the independent flow. When Ts3 is generated, two independent flows seen in Ts2 merge to one flow. Based on these results, Ts3 is not observed when c is large in Fig. 1 because the two flows generated by two impellers in each do not merge.
The relation between impeller rotational speed (n) and each of Ts2 and Ts3 was investigated. Ts2 and Ts3 were proportional to the square of n. Results show that Ts2 and Ts3 are torques based on pressure drag attributable to the relative velocity of impellers and a fluid.

PB380 Direct measurement of interaction forces between a yeast cell and a microbubble using atomic force microscopy
Satoshi KATO, Shohei YUMIYAMA, Yasuhiro KONISHI, Toshiyuki NOMURA
Osaka Prefecture University, Sakai, Japan

The removal of adhered microorganisms on surfaces is an important issue in every industrial process. The interfacial area of the gas-liquid interface increases and the residence time of bubbles in the liquid phase becomes longer with a decrease of an air bubble size. In addition, the production cost of air bubbles is low and their environmental load is also low. These positive features make air bubbles potentially suitable for use in physical cleaning technologies, such as removal of organic materials that have adhered to surfaces and cleaning of membrane fouling. However, cleaning effects using air bubbles and their mechanisms related to interaction forces are poorly understood. Thus, it is necessary to understand the interaction forces between air bubbles and target materials in detail for the further development of industrial usage of air bubbles. Atomic force microscopy (AFM) is a powerful tool for direct measurement of the interaction forces between surfaces. However, few studies have reported on quantitative measurement of interaction forces between a single microbial cell and a natural hydrophobic bubble surface. In this study, we attempted to conduct a direct measurement of the interaction forces between a yeast cell and a microbubble (MB) using AFM. We emphasize that the adhesive force between a cell and a MB measured in this study was one order of magnitude larger than those between microorganisms and hydrophobic substrates. This means that microorganisms that have only recently adhered on surfaces can be removed by adding a hydrodynamic drag force to the adhered MBs on the microbial cells. These findings provide useful information for various applications using MBs, including cleaning technology for any industrial process.

PB381 Numerical Simulation of Binary Particles in Vibrated Fluidized bed Using DEM-CFD Model Based on Particle-Scale Similarities
Jiang ZHAOHUA, Kenta RAI, Takuya TSUJI, Kimiaki WASHINO, Toshitsugu TANAKA
Osaka University, Osaka Prefecture, Japan

Segregation and mixing of binary particles with different densities are concerned in engineering. To improve its efficiency, the usage of vibrated fluidized bed which is a combination of fluidized bed and vibrated bed has been proposed. Due to the coexistence of both aeration and vibration, the phenomena of segregation and mixing of particles in vibrated fluidized bed is complex, whose detail is more important to understand for new designs. However, it is not easy to observe internal flows experimentally, considering the opaque nature of granular systems. Coupled DEM (Discrete Element Method) - CFD (Computational Fluid Dynamics) method is an attractive numerical tool to investigate particle motions in vibrated fluidized beds in detail, while large computational time is required to consider the dynamics of large amount of particles in real engineering systems. It is necessary to find out a proper approach to enable DEM-CFD calculations with reduced computational cost. In this study, enlarged particles are used while keeping the original bed size to reduce particle number. To obtain the similar segregation and mixing behavior in bed, governing equations are non-dimensionalized by characteristic variables in particle scale. Overall flow behavior is expected to be similar when the non-dimensional variables are kept same but changing particle size. In the present study, a numerical model based on particle-scale similarities is introduced and applicability of the proposed model is verified against density segregation problem in a vibrated fluidized bed by varying model particle size.

PB382 (canceled) <101532-1>
PB383 (canceled) <100207-1>