APCChE 2019

Session 2 Abstracts

Session 2. Fluid and particle processing

B202 Discuss efforts against sustainable productivity by using computational fluid dynamics —A case of Japanese chemical company—
Kazutaka TAKI1, Kei SAKAKURA2, Yasushi ISHIBA3, Atsushi KAWADA4, Takeshi SUEMASU5, Shunichi KAWANAKA6, Naoki SHIMADA7
1 TOSOH Corp.
2 Idemitsu Kosan Co. Ltd.
3 Mitsubishi Chemical Corp.
4 Mitsui Chemicals Incorp.
5 Ube Industries Ltd.
6 Kuraray Co. Ltd.
7 Sumitomo Chemical Co. Ltd.

Recently, researchers and engineers in chemical companies are facing a huge variety of flows such as multiphase flows, non-isothermal flows and flows with moving boundaries. In other hand, high quality computing becomes more important. Some member in Japanese chemical companies formed a workshop for workers on the first wave of computational fluid dynamics (CFD) in 2008, which is called as chemical company workshop on CFD (CCW)(Ishiba et al., 2010, Ishiba et al., 2013 and Kawada et al. 2017).
In this presentation, we introduce our activities and current state of affairs. As regards the quality improvement of CFD, we compared numerical results in some problems:
(a) Turbulent flow in an agitation vessel
(b) Flow in a fluidization bed
(c) Vortex generation in a swirling flow
(d) A dispersed liquid-liquid two phase flow
(e) Heat transport in packed column
We have discussed that “how do CFD engineers provide better precisions for designing chemical reactors” and “what is useful for guidance of CFD” in the last decade. We need not only development of hard ware and software, but also development of qualified engineers to keep better productions. We expect that the discuss efforts of CCW is contributing sustainable productivity by using CFD.

B203 Machine learning-based framework of descriptor selection for porous material design
Tomoki YASUDA1, Shinichi OOKAWARA1,2, Shiro YOSHIKAWA1, Hideyuki MATSUMOTO1
1 Department of Chemical Science and Engineering, School of Materials and Chemical Technology, Tokyo Institute of Technology, Tokyo, Japan
2 Department of Energy Resources, Egypt-Japan University of Science and Technology, Alexandria, Egypt

Nowadays an internal 3D structure of porous material is noninvasively analyzed with sub-micron resolution using X-ray micro-CT technology for understanding of internal phenomena. It is then expected that the important properties such as reaction efficiency, effective diffusion coefficient, permeability, thermal conductivity, mechanical strength will be predicted and designed based on the relation between the microstructure and these properties. However, a geometrical descriptor, which is an index to quantify the microstructure, can be arbitrarily defined, and hence, its appropriate selection to serve the material design is challenging and needs expert knowledge.
Therefore, this study proposed a machine learning-based framework to select the descriptors of microstructure for predicting specific material properties. For its validation, the relation between the descriptors and one specific material property of more than 1000 porous material samples were analyzed. The values of defined 62 geometrical descriptors were calculated using CT scan images. The specific porous material property was obtained by numerical simulations.
The framework employed RReliefF algorithm as filter method and support vector machine as wrapper method. Although the wrapper method is computationally expensive, its numerical load was reduced using RReliefF based on a descriptor ranking of importance to the specific material property. The framework selected 8 descriptors to minimize the prediction error for permeability. Among them, the conventional descriptors such as porosity, specific surface area, tortuosity, and hydraulic radius were also included, which validated the framework. Other selected descriptors were related to chord lengths of pores in three directions, which implied that the three-dimensional shape of pores caused the difference in the permeability of porous material with similar conventional parameters' values. It should be emphasized that the proposed framework is applicable to the identification of descriptors for predicting any specific material properties, which might not be related to the microstructure of porous material in the past.

B204 Determination of mass transfer coefficient with sparged bubble analysis in stir tank bioreactors
Yusuke TOMIOKA1, Takao ITO1, Chinatsu MINAGAWA2, Koichi KAMEKURA2
1 Merck Ltd., Tokyo, Japan
2 IHI Plant Engineering Corporation, Tokyo, Japan

Stir tank bioreactors are widely used for mammalian cell culture process for biomanufacturing. From the process development to manufacturing, various size of bioreactors are used. To enable successful scale-up, Number of factors are critical to efficient cell growth and substance production. Regarding process parameters of bioreactor, liquid mixing and gas transfer efficiency are well understood for the key process parameter and evaluate those of parameter are applied as an indicator to scale up to larger process. Volumetric mass transfer coefficient (kLa) is widely used to determine gas transfer capacity of bioreactor. kLa is affected by the bioreactor design such as impeller, vessel height and diameter, gas sparger shape and capacity. The kLa have been correlated as following equation.
kLa of bioreactors that used for process development to commercial scale has been determined by the equation. Many of engineering approach are discussed to determine kLa. Bioreactors have multiple sparger to supply oxygen for various culture condition, the generated bubbles also have different diameters and size distribution. It affects the specific interfacial area, but the size of bubble and distributions has not been verified well.
In this study, we discuss about kLa determination approach with bubble diameter and size distribution analysis regarding for the series of MobiusR bioreactor. The bioreactor family is single use bioreactors by configured flexible film-based bag for all liquid contact surface as single time usage. The bioreactor series has five tank sizes, 3L vessel for early process development, 50-200L for pilot-scale manufacturing, 1000-2000L for clinical scale manufacturing. kLa performance that of the bioreactor series are already published, but the correlation of engineering parameters from modeling are not well characterized yet. We verified kLa performance of the bioreactors with general mechanical models and analyze correlation between kLa to sparged bubble specifications. also discuss the unique sparger design and characterization.

B205 [Keynote] Hard-sphere / Pseudo-particle Modeling of Reaction-Diffusion Coupling
Wei GE1,2, Chengxiang LI1
1 State Key Laboratory of Multiphase Complex Systems, Institute of Process Engineering, Chinese Academy of Sciences
2 School of Chemical Engineering, University of Chinese Academy of Sciences

While understanding the intrinsic reaction kinetics is a perpetual pursuit of chemical science, the effective reaction kinetics most frequently faced in chemical engineering is typically significantly different from it. The coupling of reaction and diffusion is a major cause of this difference and is hence a common topic for both chemistry and chemical engineering. It is, however, a challenge topic owning to its multi-scale nature with reaction at atomistic scales and diffusion ranging from micro- to macro-scales. By simplifying the particle-particle interactions in molecular dynamics (MD) with hard sphere (HS) models and improving its scalability in parallel computation with pseudo-particle modelling (PPM), the HS-PPM approach has been demonstrated effective for bridging the scales involved in reaction-diffusion coupling and advantageous over traditional MD or continuum-based simulation methods. The presentation will introduce the theories, models, algorithms, implementation and applications of this approach with example from petrochemical and coal-chemical industries.

B213 [Keynote] Energy recovery from sewage sludge by combined hydrothermal pretreatment and CO2 gasification
Ye SHEN1, Chi-Hwa WANG2*
1 NUS Environmental Research Institute, National University of Singapore, 1 Create Way, Singapore, 138602
2 Department of Chemical and Biomolecular Engineering, National University of Singapore, 4 Engineering Drive 4, Singapore 117585

Management and disposal of sewage sludge (SS) is drawing great attraction due to its large quantity and potential environmental concern. On one hand, SS contains relatively high heating value that could be recovered through further thermal treatments such as incineration and gasification. On the other hand, high nitrogen content leads to potential NOx emission during the thermal treatments. Hydrothermal carbonization (HTC) pretreatment could enhance fuel quality and remove unwanted pollutants like nitrogen and sulfur. Thus, it is popularly applied in the field of SS disposal. In this study, an integrated HTC pretreatment and CO2 gasification system was proposed and studied. Cold gas efficiency (CGE), carbon conversion factor (CCF) and ammonia emission were evaluated as key performance indicators. Nitrogen transformation behaviour from HC was studied from 400 °C to 900 °C indicating that the release of NH3 was mainly because of the decomposition of amide group. Meanwhile, gasification products from SS and hydrochar (HC) showed that NH3 (NOx precursor) emission could be reduced by 57% through HTC pretreatment although CGE and CCF were inhibited. However, this could be improved through co-gasification of HC with horticulture waste (HW). The preferred experimental condition was found to be 1:1 HC/HW mix ratio that CGE could reach 66%. The synergetic effect was believed due to the presence of calcium compounds in the HW. It was also observed that calcium element could improve both nitrogen removal efficiency during HTC and CGE during gasification.

B215 Mixing Behaviors in Mini-Fluidized Beds
Eldin Wee Chuan LIM
National University of Singapore

Fluidized beds are used extensively in the petrochemical, metallurgical, food and pharmaceutical industries. Although a rich literature focusing almost entirely on large scale fluidized bed systems has been developed, little is known about gas fluidization in mini-channels. Such mini-fluidized beds have the potential to provide large interfacial areas, high transport rates, large wall surface area for heat exchange and may find important applications for various chemical and biological processes. However, the behaviors of mini-fluidized beds differ from their large scale counterparts and correlations that have been developed for the latter have been observed to be inadequate for predicting bubbling and slugging behaviors in mini-fluidized beds. In this study, the Discrete Element Method combined with Computational Fluid Dynamics (CFD-DEM) was applied for studies of solids mixing in a mini-fluidized bed. The effects of different operating conditions, such as superficial gas velocity, continuous and pulsating flow at the inlet, interparticle van der Waals interactions, on the efficiency with which mixing takes place were investigated through the simulations.

B216 Experimental and Numerical Studies of Particle Mixing in a Rotating Drum with an Inner Particle Motion Blocking Plate
Xuan WANG1, An-Ni HUANG1,2, Wan-Yi HSU1, Hsiu-Po KUO1,2
1 Department of Chemical and Materials Engineering, Chang Gung University, Taoyuan 33302, Taiwan
2 Department of Otolaryngology-Head & Neck Surgery, Linkou Chang Gung Memorial Hospital, Taoyuan 33305, Taiwan

Small white glass particles (0.84 – 1.00 mm) and large black glass particles (1.19 – 1.41 mm) are mixed in a 40% fill, 20 rpm drum which is 100 mm in length and 98 mm in diameter. An inner particle axial motion blocking plate is placed in a rotating drum and the effect of its position on the particle mixing is studied. The plate has 78 mm inner diameter and blocks 34% particles' axial motion. The experimental results show that the degree of mixing increases by the blocking plate shearing. The shearing changes the dynamic angle of repose of the particles and the location of the plate affects the degree of mixing and segregation patterns. The particle motion in this system is realized by DEM simulation. It is found that the plate blocks particle motion deep in the bed, but enhances axial motion at the bed surface. The plate inserting increases the particle plate crossing number, the average particle axial velocity and the particle axial moving distance per revolution.

B217 Application of eulerian granular model to evaluate solid-liquid agitation
Yasuhiro SHOJI, Kentaro OTAWARA
Kureha Corporation, Ochiai 16, Nishiki-Machi, Iwaki-shi, Fukushima, Japan

The dispersion of solid particles plays an important role in solid-liquid agitation. For example, “moderate” agitation is necessary in crystallization. Insufficient agitation causes crystals to stay at the bottom of the tank; thereby particles stick together or stick to the wall of tanks. This poses handling problems such as clogging of bottom lines and insufficient quality by coarse lumps. On the other hand, overpowering agitation causes crystals to break into small particles due to a strong shearing force, causing malfunction in subsequent processes, such as, a filtration process. Thus, it is extremely important to design both a proper apparatus and proper agitating speed to obtain desirable dispersion in solid-liquid agitation.
Designing solid-liquid agitated tanks, just suspension speed NJS and segregation state are important. To estimate NJS, Zwietering's empirical formula is widely adopted. The formula requires a geometrical constant which depends on the shapes of apparatus: The geometrical constant is often unavailable so that it needs to be obtained experimentally for each apparatus which has a lot of factors such as the shape of the agitating blade, the shape of the tank, the number of blades, and the blade position, thereby taking cost and time. Thus, NJS has been estimated with Computational Fluid Dynamics (CFD) instead of experiment. In addition, the extent of segregation has been evaluated at NJS based on particle concentration distribution obtained with CFD.
NJS, obtained with CFD simulation and experiment have been found to be adequately consistent with each other. Furthermore, the shape of apparatus has been found to affect the state of segregation at NJS.

B218 Aggregation Process of Silica Micro Particles in a Two Dimensional Chaotic Mixing Field
Yusuke OCHI1, Takafumi HORIE1, Yoshiyuki KOMODA1, Kuo-Lun TUNG2, Naoto OHMURA1
1 Kobe University, Kobe, Japan
2 National Taiwan University, Taipei, Taiwan

Processes using particle coagulation are one of prominent parts of many branches of industry. With growing interest in unsteady operations, they have been applied to agglomeration processes. Coagulates produced from unsteady operations tend to have unique properties and, in terms of industrial applications, more understanding of the relationship between unsteady flow and particle agglomeration is needed. As an unsteady flow, laminar chaotic mixing induced by adding periodical perturbation has attracted a lot of attention. In fact, it has been eagerly investigated to enhance mixing efficiency and to apply to industrial processes since 1980s. There is, however, no reports on particle aggregation behavior in a chaotic mixing field. Therefore, aiming at the development of a novel coagulation process, this study investigated particle aggregation in a simplified two-dimensional chaotic flow.
A schematic of the model flow in this study is shown in Figure 1. In this model, the system consists of two rotating cylinders. The motion of cylinders can be divided into three units protocols (protocol 1a, 1b and 2). In protocol 1a and 1b, one cylinder rotates for fixed time T. On the other hand, both cylinders rotate in protocol 2. This periodic operation is supposed to perturb the fixed point which locates at the center in between the two cylinders in protocol 2 and to induce chaotic mixing.
The temporal evolution of the number of aggregates can be described using an exponential model equation. Aggregation rate constant k showed less dependency on rotation rate of cylinders ω and higher k was obtained when T was shorter. These results indicate that chaotic agitation decreases the effect of breakage and promoted coagulation. The smaller ω or T where the system has higher unsteadiness, the larger aggregates were obtained.

B221 Mixing performance of inner baffle in shear thinning fluid
Haruki FURUKAWA, Yoshito MIZUNO, Yoshihito KATO
Nagoya Institute of Technology, Nagoya, Aichi, Japan

A mixing performance of an inner baffle was investigated in shear thinning fluid. The inner baffle was placed with a clearance between baffle and vessel wall. A traditional two-bladed paddle impeller was used. Four planar baffles was vertically placed in two ways: 1) standard baffle condition 2) baffle condition with a clearance between baffle and vessel wall. A mixing pattern was visualized by decolorization method based on an oxidation-reduction reaction with sodium thiosulfate and iodine. Flow field for shear thinning fluid in mixing vessel was measured by particle image velocimetry and was shown as stream line. An isolated mixing region (IMR) like doughnut ring formed under standard baffle condition, however, IMR dissipated when baffles were placed with a clearance between baffle and vessel wall. As a result, mixing time was drastically decreased. This study suggested that inner baffle promoted mixing in shear thinning fluid.

B222 Gelation behavior of shear-thinning hydrogels fabricated via a Kenics static mixer
Athira Sreedevi MADHAVIKUTTY1, Seiichi OHTA2, Taichi ITO1,2
1 Department of Chemical System Engineering, The University of Tokyo
2 Center for Disease Biology and Integrative Medicine, The University of Tokyo

Injectable hydrogels are widely used for minimally invasive therapy. Among injectable hydrogels, shear thinning hydrogels have received much attention recently, as their viscosity decreases with increase in shear rate making it suitable for flow through long channels. Static mixers enable efficient mixing without an external energy source, and thus have been used for the fabrication of injectable hydrogels, such as dental impression materials. However, gelation behavior of shear thinning hydrogels in static mixers is still unknown. In this research, we study the shear-thinning behavior of hydrogels formed by physically cross-linking natural polysaccharides, such as xanthan gum and locust bean gums using a Kenics Static Mixer.
Mixtures of 1 wt% xanthan gum and 1 wt % locust bean gum in pure water formed hydrogel by mixing in static mixer with 16 elements at axial velocity of 100 mm/s. The mechanical strength of the formed hydrogels was measured using a rheometer. G' value for xanthan gum to locust bean gum ratio of 1:1 was found to be 18 Pa. The mechanical strength can be increased by increasing ionic strength. The pressure drop dependence on mixing conditions and hydrogel composition at different polymer ratios was characterized using the static mixer. An increase in axial flow velocity and number of mixing elements showed an increase in average pressure drop due to hydrogel formation. Furthermore, the viscosity of formed hydrogels decreased with increasing shear rates from 1.9 to 383 s-1 and was fitted to the power law model for non-Newtonian fluids. These results indicated shear-thinning property of hydrogels, which would be owed to the nature of physical crosslinkings. Pressure drop dependence of the above parameters will be discussed using non-dimensional Deborah number defined as the ratio of time taken by materials to adjust to applied stresses to characteristic time scale of an experiment.

B223 Direct numerical simulation of shear-thickening in dilute suspensions in a viscoelastic fluid
Yuki MATSUOKA1,2, Yasuya NAKAYAMA1, Toshihisa KAJIWARA1
1 Kyushu University, Fukuoka, Japan
2 Sumitomo Bakelite Co., Ltd., Shizuoka, Japan

Particle suspensions in viscoelastic fluids such as polymer solutions or polymer melts are broadly used in industrial products. Understanding the rheology of them is essential for proper manufacturing operations. Among viscoelastic fluids, so-called Boger fluids are widely utilized for separately examining the effect of elasticity of the fluids due to the constant shear viscosity. In past experimental studies, suspensions in Boger fluids are reported to show the shear-thickening in the apparent suspension viscosity even at dilute suspensions[1]. To clarify the mechanism of this thickening, numerical studies were reported[2].
We studied unclear effects of fluid elasticity on the shear-thickening in dilute suspensions using newly developed direct numerical simulation based on Smoothed-Profile Method[3]. Oldroyd-B model, which represents the shear flow rheology of Boger fluids well, is utilized in our calculations. The fluid elasticity is quantified by the Weissenberg number (Wi) and the viscosity ratio between the solvent and polymer (beta). A small beta corresponds to a high polymer concentration or a high-molecular-weight polymer (strong elasticity of fluid).
We clarified that the modulation of the velocity field of the host fluid due to the polymer stress, which is more pronounced in lower beta conditions, relax the level of the polymer stress concentration, resulting in the non-trivial weakening of Wi-dependency of the shear-thickening in the normalized suspension viscosity observed in a high Wi region in the figure. This result indicates that the two parameters measuring the fluid elasticity has counter-acting effects on the suspension rheology and are essentially important for suspension phenomena in moderate or strong viscoelastic fluids such as non-dilute polymer solutions or polymer melts.
[1] R. Scirocco et al., J. Rheol. 49, 551 (2005)
[2] M. Yang and E. S. G. Shaqfeh, J. Rheol. 62, 1363 (2018)
[3] Y. Nakayama et al., Eur. Phys. J. E 26, 361 (2008)

B224 Quantification of the flow pattern ahead of and behind stirring blades from laminar, transitional to turbulent regime.
Yokohama National University, Yokohama, Japan

In the mixing operation, the flow pattern in the turning space of the stirring blades affects the circulation flow in the tank and is directly linked to the success or failure of the mixing operation. However, the flow pattern in the turning space is not clarified because measurement in the turning space is difficult. Especially in the transitional regime, the effect of vessel size and fluid viscosity have not been understood. Although computational fluid dynamics (CFD) is very useful technique to predict fluid flow, validation with experimental data is essential.
In this study, the three-dimensional flow velocity distribution in the turning space of the six-blade paddle in a non-baffled vessel from laminar, transitional to turbulent regime was quantified by CFD using ANSYS Fluent 18.0. Laminar, transition-SST, and LES model were used for each flow regime. In addition, the velocity distribution in the turning space was quantified by two-dimensional particle image velocimetry (2D-PIV) with the transparent impeller. The flow velocity data were converted from a Cartesian coordinate system into a rotational cylindrical coordinate system synchronized with the rotation of the blades. The flow velocity distribution of CFD results were validated by PIV data.
As a result, it was found that the discharge flow of the fluid shifted from the front side to the back side of the blade as the flow regime change from laminar, transitional to turbulence. Moreover, the fluid velocity faster than the rotational speed of the blade was observed in the back space of the blade in the transitional regime, and it was found that a negative pressure area was formed. The organized vorticity motion was formed in the back space of the blade not only in the turbulence but also transitional regime.

B225 CFD analysis of solid particle concentration distribution in solid-liquid stirred tank with various particle property and stirring condition
Shunnosuke IMAI, Ryuta MISUMI, Meguru KAMINOYAMA
Yokohama National University, Yokohama, Japan

Solid-liquid mixing is widely used in crystallizer and chemical reactor with catalysis. In these processes, there is a problem that solid particles aggregate and collide with a stirring blade. The particle collision causes particle breakage and abrasion. Furthermore, when scaling up to the scale of industrial equipment, the particle dispersion state usually becomes insufficient. It is important to quantify the solid particle concentration distribution in the tank and grasp the dispersion state 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, impeller speed N, Impeller diameter D, and tank diameter 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 LES. The interactions of particle–particle and particle–solid surfaces were modeled using the distinct element method. Analysis was performed by changing the vessel diameter, impeller speed, particle diameter and particle density to 3 conditions respectively. 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.
It was found that σc,av/Cav was correlated with the value ut/ND0.8 obtained by dividing terminal velocity of single particle ut by impeller speed N and D0.8. This result showed that the uniformity of particle concentration in the vertical direction can be ordered by the ratio of ND0.8 to ut regardless of the particle property, impeller speed and vessel size.

B226 Turbulence Modeling in Side-Entry Stirred Tank Mixing Time Determination
Suci MADHANIA, Ni'am Nisbatul Fathonah, Tantular NURTONO, Sugeng WINARDI
Chemical Engineering Department, Institut Teknologi Sepuluh Nopember, Kampus ITS Sukolilo, Surabaya, 60111, Indonesia

Mixing is one of the critical processes in the industry. One of the operating units used in the mixing process is a stirred tank. The efficiency of the stirred tank is greatly influenced by several factors including stirred tank design, operating conditions, and working fluid properties. The side-entry stirred tank is widely applied in the industrial process of crude oil in the refinery industry, water-molasses mixing in the bioethanol industry, pulp stock chest in the pulp and paper industry, and anaerobic digester for biogas reactors. Mixing time is one of the key parameters used in the design of a stirred tank. This research will model mixing time in a flat bottom cylindrical side-entry stirred tank with dimensions D = 40 cm and T = 40 cm using CFD ANSYS 18.2 by applying the RKE and SKE turbulence models. The stirrer used is three-blade marine propeller d = 4 cm which is an axial type impeller. The phenomenon of mixing in a side-entry stirred tank, qualitatively described through computational prediction results in the form of flow profiles and tracer density change contours locally. Moreover, quantitatively indicated by mixing time validated using experimental data carried out by the conductometry method. The simulation results show that the mixing time modeled with the SKE turbulence model shows a similarity level of 68.16%, while the RKE turbulence model shows 31.94%.

B301 Evaluation of agglomeration and disruption for fine powder in vibrating fluidized bed
Tomonori FUKASAWA, Akihiro YAMANE, Toru ISHIGAMI, Kunihiro FUKUI
Hiroshima University, Higashi-Hiroshima, Japan

The recycling and utilization of powdery waste are becoming increasingly important for not only ensuring a sustainable environment but also decreasing the cost of waste disposal. In our previous work, focusing on the fact that the cohesion differs depending on the powder components, we report that it is possible to concentrate and separate specific components from binary mixed powders by using a vibrating fluidized bed. In order to apply this concentration and separation technology to various powders, systematically understand the aggregation and fracture characteristics of each powder in the vibrating fluidized bed are essential. In this study, the gas fluidization of three fine powder, viz. ZnO, TiO2 rutil, and TiO2 anatase were studied in a mechanically vibrating fluidized bed. The effects of the powder particle components and vibration amplitude on agglomeration and disruption in the fluidized bed were investigated. The temporal change of apparent agglomerate size was evaluated using the Ergun equation from the temporal change of bed height and pressure drop. Furthermore, the effects of powder species and vibration amplitude on aggregate size during aggregation-disruption equilibrium were evaluated based on a energy equilibrium model.

B302 Defluidization Investigation on Bubbling Fluidized Bed Particle against Alkaline Addition during Pyrolysis Process
Gunma University, Kiryu, Gunma, Japan

Bubbling fluidization process enables to enhance performance thermal conversion reaction between bed particles and raw material like biomass inside a bed column supported by gas distribution from the particle distributor at the bottom side. However, some unwanted inorganic materials like alkaline compounds contained inside biomass have the potency to create a slagging formation with oxide compounds from either inside or outside the bed column. The pyrolysis process was selected reaction to ensure the oxide compounds that were only coming from all particles placed first inside the bed column that easily associated with alkaline compounds at the high-temperature condition. The +88-125 μm silica sand and clay (bentonite) were the bed particles of Bubbling Fluidized Bed (BFB) pyrolysis process against Potassium Chloride (KCl) addition at 700oC as the use of biomass representation under superficial-N2-velocity condition. An agglomeration formation was confirmed by the decrease of average ΔP by the experiment which will theoretically should to increase by KCl additions, while the defluidization condition was initiated by the stuck up the whole particle movement that already stop moving completely even after getting more some KCl additions. Investigation on 0.05 gr and 0.1 gr encapsulated KCl addition variables within 40 samples x 15 minutes interval reaction time showed that the silica sand has a faster tendency to form initial agglomeration than clay started from 0.1136 gr KCl addition (by 0.174% of total silica sand weight). Meanwhile, the clay has more tendency to achieve defluidization formation easily than silica sand particle on variable 0.5 gr KCl addition by the rate of ΔP decrease 0.230% and 0.213% respectively. The endurance of silica sand to maintain fluidization in a longer time than clay indicates that silica sand could be considered as the better particle to be applied in this process.

B303 Dry separation of particulate mixture based on density-segregation in a vibro-fluidized bed
1 Okayama University of Science, Okayama, Japan
2 Kyushu Institute of Technology, Kitakyushu, Japan

Separation is one of the important unit operations for industry. Wet separations such as froth flotation are utilized for particulate mixture. However, the conventional wet separations, which require large amount of water, are not convenient in dry districts. In particular, the lack of water is becoming a critical issue in mining areas with drought due to global warming. Therefore, the development of dry separations to replace the commonly used wet separations is in great demand. A gas-solid fluidized bed is one of the candidates for dry separation. When a particulate mixture having density difference is fluidized, lighter particles move up and heavier particles move down in the fluidized bed. This phenomenon is called “density-segregation”. Recently, we found that the density-segregation is enhanced when vibration is added to the fluidized bed. The vibro-fluidized bed is effective to separate particulate mixtures with small density difference. Here we investigated the dry separation of particulate ore with small density difference using the vibro-fluidized bed. We found that the density-segregation does not occur using only the fluidized bed without vibration regardless of fluidizing air velocity. In contrast, the vibro-fluidized bed enables us to separate the lighter mineral and heavier mineral based on the density segregation. The fluidizing air velocity is an important key factor to enhance the density-segregation by adding the vibration to the fluidized bed.

B304 Metallic and Bimetallic Incorporated Hydroxyapatite Catalysts Prepared via Drip Thermal Treatment Using a Fluidized Bed for Oxidative Dehydrogenation of Propane
Tsutomu NAKAZATO, Yukinori KUMEDA, Takami KAI
Kagoshima University, Kagoshima, Kagoshima, Japan

Hydroxyapatite-supported metal catalysts (M/HAp) has attracted attention in recent years especially in the oxidative dehydrogenation of light alkanes because hydroxyapatite (HAp) is a promising carrier with adjustable acid-base characteristics, high stability, and various metal substitutions. In the literature, HAp itself had shown very small catalytic activity in the oxidative dehydrogenation of propane (ODHP). Recently, we found that our heat treatment method, drip thermal treatment using a fluidized bed (DTFB), could improve the catalytic activity of nickel-incorporated hydroxyapatite (Ni-HAp) catalysts in partial oxidation of methane. In this study, HAp, metallic and bimetallic incorporated HAp (M-HAp or M1-M2-HAp) were prepared via DTFB under different atmospheres, and their catalytic activity in ODHP was investigated in a fixed bed reactor at different stoichiometric ratios of gas reactants. It was found that DTFB improved the catalytic activity of HAp itself in ODHP, and that metal species such as Co, V, and Mo were useful for metallic and bi-metallic incorporated HAp catalysts.

B305 [Keynote] Fluidization-based Reaction-decoupling Fuel Conversion Fundamentals and Technologies
Guangwen XU1,2, Zhennan HAN1,2, Xi ZENG2, Junrong YUE2
1 Shenyang University of Chemical Technology, Shenyang, 110142, China
2 State Key Laboratory of Multi-Phase Complex Systems, Institute of Process Engineering, Chinese Academy of Sciences, Beijing, 100190, China

Thermochemical conversion of fuels via pyrolysis / carbonization, cracking, gasification and combustion has to involve a number of individual reactions called attribution reactions to form an intercorrelated reaction network for any conversion process. By separating one or some attribution reactions from the others to decouple their interactions existing in the reaction network, the so-called reaction decoupling enables a better understanding of the complex thermal conversion process and further the optimization of the conditions for attribution reactions as well as the entire conversion process to realize advanced performances. By far, numerous research works have been conducted about the application of reaction-decoupling method on gasification and combustion process in our research team. Some interesting findings were obtained: (1) catalytic cracking and reforming of tar on the surface of char is critical to deeply remove tar in gasification process, and there exists a loop between deactivation and regeneration of catalytic activity of char for tar removal; (2) tar generated from fuel pyrolysis has the highest activity and also capacity of NO reduction than char and pyrolysis gas do. Based on these fundamental findings, processes of fluidized bed two-stage gasification (FBTSG) and fluidized bed decoupling combustion (FBDC) have been proposed and commercially developed to break through the problems of tar generation for gasification and NOx emission reduction for combustion, respectively. After a series of tests at laboratory and pilot scales, two industrial demonstration plants with treating capacities respectively of 10000 t/a for gasifying herb residue and 60000 t/a for burning distilled spirit lees have been designed and built to verify the technical advantages of the proposed new processes. The results show that the FBTSG technology enables tar content below 50 mg/Nm3 in product gas and the FBDC technology lowers NO emission in flue gas to about 150 mg/Nm3 even burning fuel with 4 wt.% N (dye base). This report will summarize the fundamental finding, technical innvation, pilot test and industrial demonstration for both such technologies based on reaction decoupling.

B313 Reduction of particle stiffness in DEM simulation considering attraction force
Ryosuke YAMAGAMI, Kimiaki WASHINO, Ei L. CHAN, Takuya TSUJI, Toshitsugu TANAKA
Osaka University, Suita, Japan

Discrete Element Method (DEM) has been widely used to simulate powder behavior for the past several decades. DEM has many advantages compared to the continuum model; it is possible to include attraction forces such as liquid bridge, electrostatic and van der Waals forces and it accurately predicts the motion of individual particles in a system. One of the major drawbacks of DEM is the large calculation cost. Many of the particles used in industry are stiff and the contact time is short during collisions, which may require significantly small time step (in the order of micro to nanoseconds) in simulation. This makes it difficult or even impossible to complete simulations within acceptable time scale. Therefore, in order to increase the time step and reduce calculation cost, the particle stiffness is often reduced from the original material property. Although this approach is widely accepted for dry and relatively coarse particles where only contact force is dominant, it is reported that the particle behavior changes drastically when attraction forces are exerted on particles; particles with reduced stiffness become more “cohesive” than original particles. Hence, attraction forces are sometimes reduced to counter-balance this effect. In the present study, the effectiveness and limitations of this method are further discussed.

B314 Numerical simulation of particle motions in cascade impactor and human respiratory system
Osaka Prefecture University, Osaka, Japan

Dry powder inhalations (DPIs) is a dosage form for delivering powdered medicine to the lungs by an inhaled air flow of a patient, which is mainly used as a treatment for respiratory diseases such as bronchial asthma and chronic obstructive lung disease. The site of the human respiratory system where the powder inhalation formulation reaches depends on the size of the particle. Recent researchers actively conducted the design of DPIs for effective transportation in the human respiratory system, and they generally evaluated the lung reachability using a cascade impactor. Understanding particle motions and deposition behavior in not only a cascade impactor but also a respiratory system is necessary to improve the efficiency of DPI.
In the present study, we proposed the usage of a numerical calculation model combining computational fluid dynamics (CFD) and discrete element method (DEM). We constructed a cascade impactor model and a simple lung model to calculate the behavior of the fluid and particle transportation with CFD-DEM simulations. CFD-DEM simulations using a cascade impactor model demonstrated the experimental data of transportation ratios to each stage by adjusting the adhesion based on JKR theory. Furthermore, we investigated the effect of particle properties (particle sizes and particle densities) and the inhalation patterns (healthy and patient patterns) on the lung reachability of particles. As a result, we revealed that the relationship between the aerodynamic size and the lung reachability depends exponentially. The obtained relationship is essential knowledge for the design guide of DPIs because it enables ones to estimate the lung reachability of particles from the aerodynamic size.

B315 In-line particle size monitoring and control of the jet mil by laser diffraction analyzer
Hideyuki IKEDA1, Fumiaki SATO1, Takayuki HIRA2, Ryohei IKEDA2
1 Malvern Panalytical Division of Spectris Co., Ltd.
2 Seishin Enterprise Co., Ltd.

Jet mil is widely used for the process for dry powder milling. It has advantage of milling the sample without heating up the sample, because of its adiabatic expansion feature. However, in order to maximize the throughput, keeping the particle size small enough, optimizing the milling conditions is necessary. Conditions such as inlet pressure of the air, and feeding rate must be optimized by many trials before running the process for production. Those trials take lots of time and workload. Even after determining the optimized conditions, the size of the milled particles and throughput change with the variation of the lot of the materials to be milled. The flowability of the material also changes with time, which effects to the feeding rate of the feeder, and the particle size after the milling. To solve those problems. we propose real time measurement of the particle size of the milled powder, and discuss about the feed back control of the particle size.
2.Experimental method
We put a laser diffraction particle size distribution analyzer “Insitec” (Malvern Panalytical Ltd.), at the exit of a jet mill “STJ-200” (Seishin Enterprise Co., Ltd.). We used lactose for the milling sample. The inlet pressure of the jet mill was fixed to 0.6 MPa, and changed the feeding rate of the material.
We could see the increase of the particle size, when we increased the feeding rate. This is caused by the insufficient milling force due to the too many particles in the mill. Therefore, we could change the particle size and the throughput, by just changing the feeding rate. We tried to optimize the feeding rate, using the real-time measurement data of the particle size analyzer. By using technology, the optimization of the jet mill process is possible, regardless the particle size and hardness of the materials to be milled.

B316 Suppression of re-entrainment of particles captured on medium-performance filter media by loading with atomized APA solution
Norikazu NAMIKI1, Kazuma YOSHIMURA1, Ryoichi NAKAYAMA1, Takeshi IKEDA2, Hideo KUNITOMO2, Shunsuke KOBAYASHI2
1 Kogakuin University, Hachioji, Japan
2 ESPO Chemical Co., Ltd., Chiba, Japan

The entry of micrometer-sized outdoor particulate matters into indoor spaces via the re-entrainment of these particles from their sediment layer formed on medium-performance filter media used in air handling units of office buildings might have been of a certain concern. In the present study, we investigated the effect of loading of the medium-performance filter media, where micrometer-sized particles were deposited, with APA solution on the suppression of their re-entrainment by exposing to intermittent fast air-flow blowing.
We selected two types of unwoven and charged filter media as test media, the single-layer media composed of mixed fine fibers and coarse ones and the binary-layer ones composed of upper coarse fiber layer and lower fine one separately. The test media in sheeted or pleated form were loaded with Kanto loam particles at a given loading density, and then 10 ppm of APA solution at its given density. After one day for drying the test media, the media were exposed to intermittent blowings of fast air-flow at a given face velocity to count a number of particles released from the test media. For comparison the number of released particles from test media without APA solution loading or with DI water loading were counted in the same manner.
Fig. 1 compares the percentages of decrease in number of released particles from the binary-layer or single-layer filter media in sheeted or pleated form loaded with APA solution. It was found that the single-layer media had better performance on preventing re-entrainment of particles by APA-solution loading than the binary-layer ones. In addition, the pleated media were less effective in lowering the numbers of released particles from themselves compared to the sheeted ones. It was caused by the discrepancy of loading density between test particles and APA solution in the pleated media.

B321 CFD-DEM analysis of PM filtration due to adhesion force using fluidized bed
Kento YOKOO1, Tsuyoshi YAMAMOTO1, Hideki MATSUNE1, Masahiro KISHIDA1, June TATEBAYASHI2
1 Kyushu University, Fukuoka, Japan
2 Advanced Technology Institute Inc., Kobe, Japan

Particulate matter (PM) is a well-known air pollutant and most of them are soot emitted from fossil fuel combustion. PM diameter has been decreased with improvement of combustion technology. Especially, the concentration of PM2.5 (diameter less than 2.5 μm) has increased. Fine PM affects serious damage on human body. However, it is difficult for conventional method (e.g., bag filter, diesel particulate filter, and electrostatic precipitator) to remove PM from exhaust gas of combustor. In our research group, fluidized bed is used as a PM removal device. This device collects PM2.5 at high collection efficiency by adhesion force, which is dominant for small particle. On the other hand, the performance of this device strongly depends on fluidization state. Void fraction and gas-solid drag force is large at large flow rate. High void fraction makes PM difficult to approach to bed particle. PM separates from bed particle by large drag force. Thus, it is important for development of this device to understand the mechanism of PM collection.
In this paper, numerical simulation of this device has been performed based on combined computational fluid dynamics (CFD) and discrete element method (DEM) to calculate the trajectory of each individual bed particle and PM. PM is inlet with gas from the bottom side of fluidized bed and transferred upward with gas in fluidized bed while getting adhesion force such as van der Waals force, and separation force such as gas-solid drag force. When PM approaches to bed particle and adhesion force between PM and bed particle becomes larger than separation force, PM adheres to and is held on bed particle. CFD-DEM analysis also shows that adhesion force becomes dominant for PM diameter less than several tens μm at range of superficial velocity 0.3-0.5 m/s.

B322 MLA and DEM study of enhancement of copper liberation using 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

As the demand of copper has been rapidly increasing, the concentration of copper in raw ores has been decreasing, which has become a major issue for copper production. To meet a strong demand for copper production, it is highly anticipated to develop a new process which can efficiently separate copper minerals from gangue minerals. Among the available comminution techniques, high pressure grinding roll (HPGR) is expected to achieve high mineral liberation with relatively low energy consumption. It is known that HPGR grinding can selectively break grain boundary phases which have mechanical weakness. However, the mechanism of HPGR grinding is still unclear and the effect of HPGR grinding is limited.
The objective of this study was to investigate the effect of HPGR grinding on enhancement of liberation. To achieve this objective, HPGR grinding tests were conducted. The effects of operating conditions such as pressing force, feed amount and an initial roll gap, were investigated. After HPGR grinding, liberation of copper minerals was evaluated by a mineral liberation analyzer (MLA) and image analysis of pictures obtained from MLA was conducted to measure crack size and shape in grinding products. Experimental results indicated that liberation was enhanced because of the increase of cracks by HPGR grinding. These approaches also revealed how many cracks were generated and how the cracks penetrated through an ore or a boundary phase. To support experimental results, a numerical simulation was performed to evaluate forces acting on ores during HPGR grinding. The discrete element method (DEM) coupled with T10 breakage model was applied to reproduce a HPGR grinding. Simulation results suggested that the increase of force acting on ores led to the increase of crack in an ore. Therefore, experimental and simulation results revealed that the increase of the cracks in ores during HPGR grinding enhanced liberation.

B323 [Keynote] Understanding and controlling particle mixing and de-mixing
Hsiu-Po KUO1,2, An-Ni HUANG1,2
1 Department of Chemical and Materials Engineering, Chang Gung University, Taoyuan, 33302, Taiwan
2 Department of Otolaryngology—Head & Neck Surgery, Linkou Chang Gung Memorial Hospital, Taoyuan 33305, Taiwan

Particle de-mixing may occur when particles of different physical properties are mixed. When mixing particles of different sizes, segregation bands and segregation cores may occur in the axial direction and radial direction in a rotating drum, respectively. Although theories and mechanisms of particle segregation have been proposed, a clear understanding of the underlying physics is still lacking. In this talk, size-induced particle segregation phenomena in rotating drums are reported. The mechanism of the segregation structure formation, the three-dimensional segregation patterns observed by experiments and by simulation, the control of the segregation intensity and the application of the segregation patterns will be systematically discussed.

B325 Numerical investigation of combustion process in an industrial pulverized-coal boiler
An-Ni HUANG1,2, Wan-Yi HSU1, Hsiu-Po KUO1,2
1 Department of Chemical and Materials Engineering, Chang Gung University, Taoyuan 33302, Taiwan
2 Department of Otolaryngology-Head & Neck Surgery, Linkou Chang Gung Memorial Hospital, Taoyuan 33305, Taiwan

The flow field, temperature distribution and species-in-phase distributions in a 13 MW pulverized-coal boiler are numerically studied using computational fluid dynamics (CFD). The gas phase is modelled using Eulerian frame, and the trajectories of the pulverized coal particles are described by the Lagrangian scheme. As the result, the predicted exhaust temperature and components are well agreed with the experimental data. The regions of the lower oxygen fraction locate approximately at the regions of higher temperatures. The secondary air injection mostly passes to the ceiling of boiler with very short residence time, and the trajectories of coal particles are similar to the secondary air paths. CFD technique successfully predicts the pulverized-coal combustion process in an industrial boiler.

B326 Poly (sodium 4-styrenesulfonate) decorated hydrophilic carbon black particles from waste pyrolysis process
I-Feng WU1, Guan-Ting PAN2, Jeff CHEN3, Ying-Chih LIAO1
1 Department of Chemical Engineering, National Taiwan University, Taipei, Taiwan
2 Department of Chemical Engineering, National Taipei University of Technology, Taipei, Taiwan
3 SeaWay Technology Co., Ltd.

In the past few decades, plastic products from polymers have been extensively used in daily life. In company with the convenience brought by plastic materials, plastic wastes become a critical issue for human society. One practical method to deal with used plastic products is pyrolysis process, which can turn plastic materials into short chain hydrocarbon and solid residues. Carbon black, the major residue material after plastic pyrolysis, has been widely adopted in many industrial applications. However, surface modification treatments are usually required for the recycled carbon black before used in various industrial applications, such as coating or filtration. Among various surface modification techniques, ball milling process has already been proved as a straightforward and easy mechanochemical treatment to modify surface properties of carbon black.
In this study, recycled hydrophobic carbon blacks from waste tire after pyrolysis process are used, and PSS-decorated carbon blacks are prepared through planetary ball milling. The shear force from milling process not only can decorate poly(sodium 4-styrenesulfonate) (PSS) onto carbon black, but also breaks the carbon aggregates into small fragments. The PSS-decorated carbon blacks become hydrophilic, and can be well-suspended in water. Effects of process parameters, such as milling time and solvents, on the surface properties and morphologies of PSS-decorated carbon black will be carefully examined via various analytical techniques, such as FTIR, SEM and thermogravimetric analysis. Dispersion stability tests will also be performed to prove the feasibility of making black pigment ink. Moreover, the printing quality of these inks will be examined to evaluate the possibility of applying recycled materials for printing and painting technology. In summary, this research offers general guidelines for reusing waste materials and pave the way for combining recycled materials with advanced technologies.

B402 Simulation of polymer melt flow in a channel by using a hieratical approach
Takashi TANIGUCHI, Yuji HAMADA, Takeshi SATO
Department of Chemical Engineering, Kyoto University, Kyotodaigaku-Katsura, Nishikyo-Ku, Kyoto 615-8510 Japan

We have successfully applied a multiscale simulation (MSS) method [Murashima and Taniguchi, Europhys. Lett., 96, 18002 (2011)] to flows of a mono-dispersed linear entangled polymer melt in a contraction-expansion channel and in concentric cylinders. In our MSS method, a macroscopic model is coupled with a microscopic model through the velocity gradient tensor and the stress tensor. As the macroscopic model, the Smoothed Particle Hydrodynamics (SPH) method is employed, on the other hand, as the microscopic model, a slip-link model is employed. By using our MSS method, two-dimensional flows in a 4:1:4 contraction-expansion channel and axial flow in concentric cylinders are examined. From our multiscale simulations, we have evaluated detailed microscopic information that comes from the polymer chain dynamics, such as the local orientation of polymer chains, the average number of entanglements and the number density of entanglements along a polymer chain. In addition, we found that the number of entanglements generally decreases in the middle section of a polymer chain under flows. To explain this behavior, an equation for the entanglement density on a primitive path is proposed. These microscopic insights will bring us a stepping stone to design a polymer melt that has a specific property.

B403 (canceled) <100095-1>
B404 Flow characteristics of drag-reducing surfactant solution
Takashi SAEKI, Aya KAIDE
Yamaguchi University

It is well known that certain cationic surfactant solutions with a counter ion cause significant drag reduction. In comparison with the coherent structure of the normal turbulence, vortex motion of drag-reducing flow can be restrained by the rheological property of surfactant solutions; consequently flow pattern might be greatly changed. The objective of this study is to clarify the flow characteristics of drag-reducing flow and to discuss the mechanism of the phenomenon.
The surfactant and counter ion we used were Ethoquad O/12 (C18H35N+(CH3)(C2H4OH)2Cl- and sodium salicylate (HOC6H4COO- Na+), respectively. Drag reduction rate (DR%) was measured by using PVC pipeline system having the diameter of 25 mm. We also measured the flow characteristics of the surfactant solutions with different concentrations and temperatures by using particle image velocimetry (PIV). From the experimental results, DR% at average velocity of 2.0 m/s increased from zero to 73% with the increase of the temperature from 10 to 50 °C at constant concentration of the surfactant (300 mg/L). From the velocity contour plot obtained from PIV, we found that the thickness of the lower velocity region of the drag-reducing flow was thick at 10 °C, while vortex motions seems controlled at 30 °C. On the contrary, the lower velocity region thickened with the increase of the concentration of the surfactant with constant temperature. Even if the same level of drag-reducing effects are taking place, we expected that the flow patterns are quite different depending on the concentration and temperature.

B405 Investigation of the effect of nanoparticles on microscopic wetting behavior of a nanoparticle suspension droplet on a substrate using phase-shifting ellipsometer
Eita SHOJI, Tatsuya YONEMURA, Takahiro KANEKO, Masaki KUBO, Takao TSUKADA, Atsuki KOMIYA
Tohoku University, Sendai, Japan

A nanofluid is a suspension of nanoparticles in a liquid and has attracted attention in diverse engineering fields, e.g. the nanofluid is expected as nano-ink for printed electronics. In such a situation, the knowledge of the effect of nanoparticles on wetting phenomena of nanofluid on a substrate is required. Therefore, the objective of this study is to investigate the effect of nanoparticles on wetting behavior in microscopic view using a phase-shifting ellipsometry we developed.
In this study, polydimethylsiloxane (PDMS) was used as a solvent, and PDMS-modified SiO2 nanoparticles were added within 3 wt% in the solvent. The particle size distributions of the nanoparticle suspension were evaluated using dynamic light scattering method, and consequently aggregates with a micrometer size were observed although the mean diameter was approximately 600 nm. In experiment, a droplet of pure PDMS or SiO2 nanoparticle suspension was put on a silicon substrate, and the microscopic wetting behavior near the contact line was quantitatively visualized by the phase-shifting ellipsometer. The phase-shifting ellipsometer can measure the two-dimensional thickness distribution of a liquid film at the nanoscale. Relations of the contact angle and the length of precursor film with a nanometer-scale thickness ahead of the contact line with capillary number were evaluated from the obtained thickness distributions to analyze the effect of nanoparticles.
The large difference in contact angle between the pure PDMS and the suspensions at the same capillary number was observed though the viscosity change with adding the nanoparticles is small. On the other hand, there was little difference in the precursor film lengths of both fluids at the same capillary number. These results indicate that the nanoparticles did not affect the precursor film region.

B406 Evaluation method for emulsification capacity using maximum shear rate distribution of liquid given by emulsification device
Shiseido Co., Ltd.

Most milky lotions and creams as cosmetics are emulsions and non-Newtonian liquids. In the emulsifying process of these manufacturing, a rotor-stator type emulsifying device is widely used. Each emulsifying device has a different shape of the rotating blades and stators. The size of the emulsified particles of the resulting emulsions varies depending on the device and its rotational velocity of the blades. The higher the shear rate produced by the emulsifying device, the smaller the size of emulsified particles produced. That is, it is considered that the size of the emulsified particles is greatly affected by the maximum value of the shear rate that the liquid receives as it passes through the emulsifying device. The emulsifying device is designed to generate a high shear rate. But the shear rate that the liquid receives is distributed because it depends on the route through the device. So far, we have evaluated the capabilities of emulsifying devices by experiments using specific model formulations. Today, we obtain the distribution diagram of the maximum shear rate that the liquid receives by numerical analysis. This distribution diagram can be used as an indicator of the capability of the emulsifying device.
In addition, the emulsifying device works like a pump that sucks and discharges liquid by high-speed rotation of the blades. This discharge amount is the processing capacity of the apparatus itself. When it is used in a mixing kettle, it also contributes greatly to total mixing.
This paper reports the evaluation of the capacity of the emulsifying devices by calculation of their maximum shear rate distributions and discharge amounts by using numerical analysis.

B413 Effects of emulsification parameter on droplet size distribution using SPG membrane emulsification
Ayaka NAKAJIMA, Masakazu NAYA, Hiroaki MATSUKAWA, Katsuto OTAKE, Atsushi SHONO
Tokyo University of Science, Japan, Tokyo

The technology of emulsions is used in various industrial application such as cosmetics, foods and pharmaceutics. It is also well known that mono-dispersed distribution is required stability. We focused on emulsions preparation by Shirasu Porous Glass (SPG) membrane emulsification. Compared to conventional emulsification method, this emulsification method doesn't require strong shearing force to prepare emulsions. It is also noted that the emulsions prepared by SPG membrane has mono-dispersed droplet size distribution and the droplet size. SPG membrane emulsification method requires only less energy input when prepare the same droplet size.
In this study, Oil-in-water (O/W) microemulsions prepared using SPG membrane emulsification method. Soya oil and ion-exchanged water were used as dispersed and continuous phase respectively. Sodium dodecyl sulfate (SDS) was used as anionic surfactant, and it added to ion-exchanged water. Soya oil was forced into continuous phase through SPG membrane (membrane pore size is 5 or 20 μm, and its surface is hydrophilic) using syringe pump. To disperse the oil droplet, the continuous phase was stirred by magnetic bar.
The effect of volumetric flow rate, stirring speed, and concentration of surfactants were examined. In this abstract, only about changing volumetric flow rate is described. In the condition that volumetric flow rate was changed 4 steps from 0.8 ml h-1 to 6.4 ml h-1 average droplet size, average droplet size had no significant change on each condition. Interestingly, these droplet sizes were about 3-4 times as large as membrane pore size (5 μm). On the other hand, coefficient of variation (CV), which is index of mono-dispersed distribution, increased with increasing volumetric flow rate.

B414 Microdroplet generation in microfluidic flow-focusing device
Narin PAIBOON1, Suvimol SURASSMO2, Uracha Rangsardthong Ruktanonchai2, Apinan Soottitantawat1
1 Center of Excellence in Particle and Material Processing Technology, Department of Chemical Engineering, Faculty of Engineering,
Chulalongkorn University
, Bangkok 10330, Thailand
2 National Nanotechnology Center (NANOTEC), National Science and Technology Development Agency 111 Thailand Science Park, Pathumtani 12120, Thailand

Microfluidics, the potentiality of generating monodisperse droplets, is an interesting tool in the encapsulation process due to its capability of microparticles fabrication. Droplet generation which occurs by two immiscible fluid is controlled by providing of flow- or pressure-driven approach. To combining the size of droplet in both physical and operating parameters, using the dimensionless number characterize the flow motions in different systems was studied. The dependence of the droplet size/channel size on the flow rate ratio (Qc/Qd), viscosity ratio (λcd), Reynolds number (Re) and Capillary number (Ca) was investigated for a wide range of fluid properties and flow conditions. The data was collected from the two-channel size of the 100 and 190 μm etch depth of microfluidic flow-focusing devices. An increase of the flow rate ratio resulted in a decrease of the droplet size. The viscosity ratio played an important role in droplet generation. This study also reveals that the droplet size and flow regime could be correlated with the Re and Ca of both phases. A diagram of a changing of flow regime was obtained. Moreover, the results indicated that the monodisperse droplet not only occurs in a dripping regime but also could be produced in a tip-streaming regime.

B415 Preparation of magnetorheological fluid using stabilizing additives
Aya KAIDE1, Makoto KANDA1,2, Takashi SAEKI1, Hiroshi TOCHIGI2
1 Yamaguchi University, Ube, Japan
2 Cosmo Oil Lubricants Co., Ltd., Satte, Japan

Magnetorheological fluid (MRF) attracts attention as a functional fluid of which viscosity can be controlled by magnetic field. Since MRF is a mixture of a base oil and magnetic particles, we often find particle sedimentation and/or oil separation due to the density difference and familiarity of the particles and the base oil. The objective of this study is to increase the sedimentation stability of MRF by using suitable stabilizing additives. Also, we consider the mechanism of stabilization with the connection of the network structure build by the additives and interaction with the magnetic particles.
We selected three stabilizing additives; a fumed silica, a styrene-based polymer, and an organogelator (PMDA-2C8/oleyl) and prepared MRF using such additives. We measured steady state flow characteristics and time-dependency of the viscosity of MRF by using a rheometer; Rheologia A-300, Elquest. Stability test was also conducted using a graduated cylinder. Morphological observation of the stabilizing additives was conducted by using a transmission electron microscope (TEM). From the experimental results, we found thixotropic behavior both for MFR prepared with the fumed silica and PMDA-2C8/oleyl. Particularly, PMDA-2C8/oleyl was expected to build internal structure due to the network formed by self-assembly of molecules. Although the viscosity of MRF increased with the addition of the styrene-based polymer, the MRF displayed no viscosity time-dependency. The styrene-based polymer can increase the yield stress of MRF, consequently prevent the sedimentation of particles. From TEM images, we can observe that fumed silica particles connected alternately, while the self-assembly of PMDA-2C8/oleyl formed network structure in the oil.

B416 Effect of hydrophobic fumed silica on the sedimentation stability of magnetorheological fluid
Makoto KANDA1, Hiroshi TOCHIGI1, Takashi SAEKI2, Aya KAIDE2
1 Cosmo Oil Lubricants Co., Ltd., Satte-shi, Saitama, Japan
2 Yamaguchi University, Ube-shi, Yamaguchi, Japan

Magnetorheological fluid "MRF" is a slurry prepared with an oil and magnetic particles having the average particle diameter of several μm. Since the viscosity of MRF can be reversibly controlled by a magnetic field, MRF is classified as a functional fluid. The quick and large variation of the rheological property of MRF with respect to a magnetic field is convenient to the application such as dampers, brakes, and clutches of automobiles. In order to show good fluidity of MRF under no magnetic field conditions, suitable dispersant is generally added to prevent viscosity increase due to particle aggregation. Moreover, MRF has large difference of density between the oil and particles, consequently shows particle sedimentation and liquid separation of the upper layer. Therefore, the development of the technique to increase the stability of MRF is necessary.
In this study, we focused on the effect of fumed silica as a stabilizer, which is widely used for paint industries. We prepared MRF with oleic acid as a dispersant and several fumed silica additives <see table> with different surface characteristics. The rheological measurements and stability tests were conducted for thus obtained samples. From the results, we found that the degree of hydrophobization and the concentration of silanol groups of fumed silica affected to the rheological properties and stability of MRF. In addition, MRF prepared with hydrophilic fumed silica showed further reduction of viscosity by the addition of oleic acid as compared to hydrophobic fumed silica.

B417 Direct numerical simulation of flow resistivity and oil droplet coalescence on x-ray ct images of nonwoven fabrics filters
Mohammad Irwan Fatkhur ROZY, Masaki UEDA, Tomonori FUKASAWA, Toru ISHIGAMI, Kunihiro FUKUI
Department of Chemical Engineering, Graduate School of Engineering, Hiroshima University, Higashi-Hiroshima, Hiroshima, Japan

In this study, we numerically study and experimentally validate the flow resistivity of commercial nonwoven fabrics filters used for a bag filter system. A numerical method that coordinates the filter structure obtained by X-ray CT imaging with computational fluid dynamics using immersed boundary method is developed to represent a realistic flow field inside the filter during simulation. We investigated the effect of superficial velocity, porosity of the filter domain, and the type of filter then analyzed the results using Darcy's law. The results showed that the calculation results from our numerical method were quantitatively in good agreement with the experimental measurement. We also verify that the Kozeny constant can be estimated by utilizing the solid volume fraction defined by immersed boundary method. These results show that our simulation method can be used to clarify the effects of porosity, fiber arrangement, and fiber shape on the pressure drop. We will discuss the application of this method to coalescence of oil droplets in oil-water separation process.

B418 Flow characteristics of non-Newtonian fluid (non-spinnability fluid)
Katsuhide TAKENAKA1, Yoshihiro YOKOKAWA2, Aya TANAKA2, Koichiro ONISHI1
1 Sumitomo Heavy Industries Process Equipment Co., Ltd., Saijo, Ehime, Japan
2 Shiseido Co., Ltd., Osaka, Osaka, Japan

There are lots of stuff, which is subject to the rheological issue, in a market. Typical examples are grease for mechanical industries, cream and hair gel for cosmetic industries, and gelatin for food industries.
On the other hand, there is few number of the research for mixing dealing with rheology due to complex factors, in spite of facing many of mixing problems in industries. One of the important factors is fluid viscosity. Therefore, as guided by N. Harnby et al. (Mixing in the Process Industries, 2nd edition issued 1997), in case of non-Newtonian fluid, mixing Reynolds number is modified with apparent viscosity, which is converted with coefficient of shearing rate (Kagakukogaku-Binran) proposed each type of impeller.
In this study, comparing with spinnability fluid and non-spinnability fluid (no 1st normal force difference fluid), as non-Newtonian fluid has been done, in terms of flow characteristics.
The mixing power has been measured with similar viscosity, the power for spinnability fluid is proportional to cubic of impeller speed for all type of impeller we have investigated. This is typical behavior in turbulence region. On the other hand, the mixing power for non-spinnability fluid perform differently. Exponent of impeller speed change from approximately 2 to 3, depending on impeller shape. However, flow is not observed, especially at the condition obtained exponent 2 of impeller speed.

B421 Distribution of local viscosity of viscoelastic solutions measured by using optical tweezers
Ruri HIDEMA1, Hikari TAKAHASHI1, Zenji YATABE2, Hiroshi SUZUKI1
1 Kobe University, Kobe, Hyogo, Japan
2 Kumamoto University, Kumamoto, Kumamoto, Japan

Solution containing high-molecular-weight polymer or surfactant forming rod-like micelles show very complex flow behavior. Non-linear flow behavior became prominent especially under high shear rates or high extensional rates. That is due to non-uniform structure of dilute polymer or rod-like micelles solution, which is enhanced under the condition of high strain rates.
In this study, in order to detect non-uniform structure in complex fluids, optical tweezers was used to measure local viscosity of several Newtonian and viscoelastic solutions. Ethylene glycol solution, polyethyleneoxide (PEO, Mw = 3,500,000) solution and cationic surfactant solution were prepared as Newtonian and viscoelastic solutions. The cationic surfactants form rod-like micelles in the solution with the presence of counter ion supplier, which became highly viscoelastic. The average values of the local viscosity of these solutions measured by the optical tweezers were consistent with the viscosity measured by a rheometer. Then, distributions of the local viscosities were also obtained to quantify non-uniform structure of the solution. The local viscosity distributions of viscoelastic solutions were wider than that of Newtonian fluids. Especially, the distribution of rod-like micelles solution was much wider than that of the other solutions. The results suggest that the local viscosity distribution of each solution has non-uniformity information of viscoelastic complex fluids.

B422 Computations of Curtain Coating
Ryota NAKAJIMA1, Takumi CHIBA1, Shunji HOMMA1, Donghyuk KANG1, Hiroyuki HIRAHARA1, Hirokazu UMEMIYA2, Sumihisa ODA2
1 Saitama University, Saitama, Japan
2 Saiden Chemical Industry Co., Ltd., Saitama, Japan

For the production of self-adhesive labels, high speed coating of adhesives is necessary to increase the productivity of labels. It is also required that the coated adhesives is as thin as possible to reduce the use of adhesives. Curtain coating is one of the candidate to achieve high speed and thin coating. However, curtain coating can cause some troubles like air entrapment resulting in the defect of the products. In this study, three-dimensional numerical simulations of curtain coating were carried out in order to investigate the mechanism of entrapment of air to coated films. The OpenFOAM tool box with VOF method was used for the simulation.
Simulation results reproduced four different regimes in curtain coating: coatable, heel formation, air entrainment, and air entrainment with heel formation. They are mapped on the diagram of the curtain Reynolds number versus the velocity ratios between coating speed (substrate speed) and impingement one. We will discuss the mechanism of air entrapment using numerically visualized flow fields.

B423 Flow classification and its application to fluid processing
Kyushu University, Fukuoka, Japan

In fluid processes in engineering fields, flow pattern characters are essential in the process performance, and are largely controlled by the geometric structures of the channel and/or driving element. In melt mixing in polymer processing, various elements, including screw and mixing elements, defines the channel geometry, and the resulting flow patterns are responsible for different mixing characteristics. Analyzing the flow pattern should be useful in understanding the mixing characteristics in connection with the different geometric structures of the mixing elements. Mathematically, the classification of flow field can be analyzed by the eigen-analysis of the deformation rate tensor [1,2]. However, such analysis technique have not been commonly applied in fluid processing. We have developed a simplified flow classification scheme based on the invariants of the strain-rate tensor and the vorticity tensor, which is frame-invariant [3,4] . The spatial distributions of these quantities offer an essential tool in understanding the flow pattern structure, and therefore can be useful to get insights into the connection between the geometry and the process performance.
[1] Tanner, R. I., and R. R. Huilgol, Rheol. Acta 14.11 (1975) 959–962.
[2] Chong, Min S., Anthony E. Perry, and Brian J. Cantwell, Phys. Fluids A 2.5 (1990) 765–777.
[3] Nakayama, Y., et al. AIChE J. 62.7 (2016) 2563–2569.
[4] Nakayama, Y., et al., AIChE J. 64.4 (2018) 1424–1434.

B424 Effect of agitation speed on ice cream properties in batch freezer
1 School of Food and Nutritional Science, University of Shizuoka, 52-1 Yada, Suruga, Shizuoka, Japan
2 Complex Fluid and Thermal Engineering Research Center (COFTEC), Kobe University, 1-1 Rokkodai, Nada, Kobe, Japan

Ice cream is one of the most popular deserts in the world. Its physical and sensory properties depend on the microstructure of ice cream. Ice cream consists of bubble, ice, fat and the freeze-concentrated continuous phase. Because all the phases intricately interact each other, it is quite difficult to construct the precise principle for the control of ice cream properties in the production process. The microstructure of ice cream is determined in freezing process. During freezing process of the ice cream mix in an ice cream freezer, air is whipped into the mix. Air taken in the ice cream is broken down to small bubbles by the shear force. Also, the shear force induces fat globules coalescence which affects ice cream properties. Therefore, it can be considered that the shear force generated by the agitation operation dominates the microstructure of ice cream. In order to control ice cream properties from the perspective of agitation, in this study, the effect of agitation speed on the viscosity of the ice cream and the bubble size in the ice cream was investigated.
The commercial batch ice cream freezer with an impeller was used. During freezing process, the torque was monitored by a torque meter (ST3000II, SATAKE Chemical Equipment MFG. LTD.). Based on the measured torque, the apparent viscosity was estimated. The bubbles size in the ice cream was calculated from the image acquired by a digital micro scope (DMS1000B, Leica Microsystems LDT.).
The apparent viscosity decreased with the increase in the agitation speed. This means the agitation speed has an influence on the sensed creaminess during consumption. Furthermore, it was found that the increase in the agitation speed leads the increase in smaller bubbles with a narrower size distribution. These results imply that the rheological properties can be controlled by the agitation operation.

B425 [Keynote] Bubble dynamics in gas-solid two-phase fluidised beds
Siddhartha SHRETHA, Shibo KUANG, Aibing YU, Zongyan ZHOU
ARC Research Hub for Computational Particle Technology, Department of Chemical Engineering, Monash University, VIC3800, Australia

The gas-solid fluidized beds have been broadly studied, both experimentally and theoretically because they are widely applied in several industrial applications. However, most of the previous studies reported their findings based upon spherical particles while in practice, particles in gas-solid fluidized beds are generally non-spherical. In this study, the combined approach of CFD-DEM (Computational Fluid Dynamics - Discrete Element Method) is used to simulate the gas-solid fluidization process of ellipsoidal particles to analyse the effect of particle shape on bubble dynamics. Simulations for single jet fluidized bed by injecting gas jet only through central orifice showed that the bubbles for ellipsoidal particles are larger in size, have a greater bubble size distribution and higher bubble break-up frequency. The bubbles for ellipsoidal particles have lateral drift, therefore, while rising through the bed they travel away from bed centreline. Simulations for uniform gas injection demarcated that the bubble and solid flow patterns are asymmetrical for ellipsoidal particles compared to spherical particles. The bubble size is smaller and the bubble size distribution is broader for ellipsoidal particles. Both types of simulations exhibited that the difference in particle shape can result in different bubble behaviour, for example, a region of bubble formation, bubble trajectory and bubble properties. Moreover, ellipsoidal particles have preferred orientation which can depend upon the development of bubbles. The ellipsoidal particles align parallel to the fluid flow direction around the bubbles while they align perpendicular to the fluid flow direction at the top of the bubble. The findings deduced from this study can aid in the understanding of particle shape effect on gas-solid fluidized beds.

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>