SEARCH
検索詳細
林 公祐大学院工学研究科 機械工学専攻教授
研究活動情報
■ 受賞- 2023年08月 日本混相流学会, 論文賞, 鉛直管におけるフラッディング状態での管内流動特性
- 2022年08月 日本混相流学会, 論文賞, 変形気泡に働くせん断誘起揚力に関する研究
- 2021年10月 神戸大学工学部, 工学部優秀教育賞
- 2021年03月 IOP publishing, Outstanding Reviewer Award 2020 (Fluid Dynamics Research)
- 2020年11月 マツダ財団, マツダ研究助成奨励賞
- 2017年08月 日本混相流学会, 日本混相流学会「論文賞」, 鉛直管における気液対向流制限日本国その他の賞
- 2013年05月 15th International Topical Meeting on Nuclear Reactor Thermalhydraulics (NURETH-15), Best Poster Award, Void Distribution and Bubble Motion in Bubbly Flows in a 4 X 4 Rod BundleThe “best” poster based on the structure of the presentation (visual quality), イタリア共和国国際学会・会議・シンポジウム等の賞
- 2013年03月 日本学術会議IUTAM分科会, 第62回理論応用力学講演会 優秀講演賞, 粘性力及び表面張力の評価法が界面運動の数値予測に及ぼす影響36歳未満の若手研究者による優秀な講演に対して授与される国内学会・会議・シンポジウム等の賞
- 2010年07月 日本混相流学会, 日本混相流学会 奨励賞, 円管内気泡及び液滴の終端上昇速度に関する研究「奨励賞」は、本会会誌または他の学術誌、あるいは本会主催・共催の国際会議論文集に掲載された論文または技術報告（原則として３年以内）を対象とし、際立った独創性、新規性、将来性などが期待され40歳以下の本会会員に贈呈する。国内外の国際的学術賞
- 2006年02月 日本機械学会, 第18回計算力学講演会 学生優秀講演賞, 表面張力評価へのカーン‐ヒリヤード方程式の応用国内学会・会議・シンポジウム等の賞
- 2004年09月 日本混相流学会, 日本混相流学会年会講演会2004 学生優秀講演賞, 改良SCAによる流体粒子の界面追跡計算国内学会・会議・シンポジウム等の賞
- Abstract Experiments on contaminated Taylor flows in a square microchannel were carried out to investigate the effects of surfactant on the bubble shape in the nose and tail regions for different surfactant properties. The nose curvature was found to be proportional to the bubble length at low surfactant concentrations, while it was independent of the concentration at high concentrations. The rate of increase in the nose curvature at the former concentrations can be expressed in terms of the surface coverage ratio. The bubble velocity decreased with increasing the nose curvature, whereas the surface tension reduced by surfactant adsorption worked better to correlate the velocity data. The curvature of the bubble tail increased steeply at low concentrations as a consequence of the early coverage due to interfacial advection. The tail curvature also had a strong correlation with the surface coverage ratio.Springer Science and Business Media LLC, 2025年01月, Microfluidics and Nanofluidics, 29, #8, 英語研究論文（学術雑誌）
- 2024年12月, 混相流, 38(4) (4), 404 - 414, 日本語鉛直円管内旋回環状流の液膜厚さおよび圧力勾配の一次元三流体モデルによる予測[査読有り]研究論文（学術雑誌）
- 2024年12月, 混相流, 38(4) (4), 423 - 434, 英語Effects of Diameters on Flow Characteristics in Vertical Pipes under Flooding Conditions[査読有り]研究論文（学術雑誌）
- 2024年11月, 20th Brazilian Congress of Thermal Sciences and Engineering (ENCIT2024), 1 - 8, 英語Velocity Field External to a Tube Bundle Applying Spatial Filter Velocimetry and Machine Learning[査読有り]研究論文（国際会議プロシーディングス）
- 2024年11月, 20th Brazilian Congress of Thermal Sciences and Engineering (ENCIT2024), 1 - 8, 英語Experimental Analysis of Velocity Fields within the Bubble Wakes[査読有り]研究論文（国際会議プロシーディングス）
- 2024年11月, The 34th International Symposium on Transport Phenomena, 1 - 5, 英語Experimental Analysis of Bubble Terminal Velocity in the Presence of Fibers in Aqueous Solutions[査読有り]研究論文（国際会議プロシーディングス）
- 2024年11月, The 34th International Symposium on Transport Phenomena, 1 - 4, 英語Fully-Contaminated State of Single Bubbles in Surfactant-laden Water at Concentrations beyond CMC[査読有り]研究論文（国際会議プロシーディングス）
- 2024年11月, The 34th International Symposium on Transport Phenomena, 1 - 4, 英語Effects of Surfactant on Shape of Taylor Bubbles in a Square Microchannel[査読有り]研究論文（国際会議プロシーディングス）
- Abstract Interface tracking simulations of gas–liquid Taylor flow in horizontal square microchannels were carried out to understand the relation between the pressure drop in the bubble part and the curvatures at the nose and tail of a bubble. Numerical conditions ranged for 0.00159 ≤ CaT ≤ 0.0989, 0.0817 ≤ WeT ≤ 25.4, and 8.33 ≤ ReT ≤ 791, where CaT, WeT, and ReT are the capillary, Weber, and Reynolds numbers based on the total volumetric flux. The dimensionless pressure drop in the bubble part increased with increasing the capillary number and the Weber number. The curvature at the nose of a bubble increased and that at the tail of a bubble decreased as the capillary number increased. The variation of the curvature at the tail of a bubble was more remarkable than that at the nose of a bubble due to the increase in the Weber number, which was the main cause of large pressure drop in the bubble part at the same capillary number. The relation between the bubble velocity and the total volumetric flux was also discussed. The distribution parameter of the drift-flux model without inertial effects showed a simple relation with the capillary number. A correlation of the distribution parameter, which is expressed in terms of the capillary number and the Weber number, was developed and was confirmed to give good predictions of the bubble velocity.Springer Science and Business Media LLC, 2024年07月, Microfluidics and Nanofluidics, 28(8) (8)[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2024年07月, Chemical Engineering Science, 120557 - 120557[査読有り]研究論文（学術雑誌）
- Begell House, 2024年06月, Multiphase Science and Technology, 36(2) (2), 41 - 54, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2024年06月, Multiphase Science and Technology, 36(2) (2), 55 - 73, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2024年04月, Chemical Engineering Research and Design, 204, 343 - 353, 英語[査読有り]研究論文（学術雑誌）
- BACKGROUND The treatment of vertebrobasilar junction (VBJ) aneurysms is challenging. Although flow diverters (FDs) are a possible treatment option, geometrical conditions hinder intervention. VBJ aneurysms possess dual inflow vessels from the bilateral vertebral arteries (VAs), one of which is ideally occluded prior to FD treatment. However, it remains unclear which VA should be occluded. OBSERVATIONS A 75-year-old male with a growing VBJ complex aneurysm exhibiting invagination toward the brainstem and causing perifocal edema required intervention. Preoperative computational fluid dynamics (CFD) analysis demonstrated that left VA occlusion would result in more stagnant flow and less impingement of flow than right VA occlusion. According to the simulated strategy, surgical clipping of the left VA just proximal to the aneurysm was performed, followed by FD placement from the basilar artery trunk to the right VA. The patient demonstrated tolerance of the VA occlusion, and follow-up computed tomography angiography at 18 months after FD treatment confirmed the disappearance of the aneurysm. LESSONS Preoperative flow dynamics simulations using CFD analysis can reveal an optimal treatment strategy involving a hybrid surgery that combines FD placement and direct surgical occlusion for a VBJ complex aneurysm.Journal of Neurosurgery Publishing Group (JNSPG), 2024年03月, Journal of Neurosurgery: Case Lessons, 7(10) (10), CASE23736, 英語[査読有り]研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2023年12月, 混相流, 37(4) (4), 401 - 411, 日本語[査読有り]研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2023年12月, 混相流, 37(4) (4), 392 - 400, 日本語[査読有り]研究論文（学術雑誌）
- 2023年11月, 12th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF 2023), 1 - 3, 英語Effects of Surfactant on Bubbly Flow in Narrow Rectangular Column[査読有り]研究論文（国際会議プロシーディングス）
- 2023年11月, 12th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF 2023), 1 - 4, 英語Measurements of Surfactant Amount at the Gas-Liquid Interface in Contaminated Taylor Flows in a Microchannel[査読有り]研究論文（国際会議プロシーディングス）
- 2023年11月, 12th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF 2023), 1 - 3, 英語Drag Coefficients of Particle-Covered Bubbles Produced Using a Microchannel[査読有り]研究論文（国際会議プロシーディングス）
- Elsevier BV, 2023年11月, Chemical Engineering Research and Design[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2023年10月, Nuclear Engineering and Design, 412, 112476 - 112476[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2023年09月, Experimental Thermal and Fluid Science, 147, 110947 - 110947, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2023年05月, Chemical Engineering Research and Design, 193, 75 - 84, 英語[査読有り]研究論文（学術雑誌）
- Informa UK Limited, 2023年03月, Heat Transfer Engineering, 1 - 14, 英語[査読有り][招待有り]研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2023年01月, JAPANESE JOURNAL OF MULTIPHASE FLOW, 37(1) (1), 128 - 137, 英語研究論文（学術雑誌）
- Begell House, 2023年, Multiphase Science and Technology, 35(1) (1), 55 - 68, 英語
Effects of surfactants on lift coefficients, <i>C<sub>L</sub></i>, of single ellipsoidal bubbles rising through linear shear flows were investigated. Two types of surface-active agents, Triton X-100 and 1-octanol, were used. The liquid properties except for the surface tension were identical to those in a clean system of logM = -5.5, where M is the Morton number. The range of the bubble Reynolds number was 0.1 < Re < 70. Bubble shapes were either spherical or ellipsoidal. Comparing with clean bubbles, less deformation of contaminated bubbles was confirmed. A shape correlation without taking the dimensionless shear rate into account gave good evaluations of the bubble aspect ratio, which means that the shear rate is not a dominant factor causing the change of shape deformation. However, drag coefficients were affected by the shear rate. A new correlation of drag coefficients was deduced, which agreed well with the experimental data. Both clean and contaminated <i>C<sub>L</sub></i> data showed a similar tendency: after a drastic decrease in <i>C<sub>L</sub></i>, <i>C<sub>L</sub></i> slightly increases and then decreases to a negative value while increasing the bubble Reynolds number, Re. A difference in concentration of Trion X-100 resulted in only a slight change in <i>C<sub>L</sub></i> at high Re regime. Different types of surfactant resulted in noticeably different values of <i>C<sub>L</sub></i> especially at low Re. The <i>C<sub>L</sub></i> of small bubbles in contaminated systems could be reproduced by a correlation for solid particles, supporting that fully contaminated spherical bubbles behave like solid spheres. For deformed bubbles, <i>C<sub>L</sub></i> can be expressed by relating the negative lift force due to shape deformation with the drag force.
[査読有り]研究論文（学術雑誌） - Begell House, 2023年, Multiphase Science and Technology, 35(2) (2), 29 - 41, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2022年12月, Nuclear Engineering and Design, 399, 112001 - 112001, 英語[査読有り]研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2022年11月, 混相流, 36(4) (4), 440 - 451, 英語[査読有り]研究論文（学術雑誌）
- 2022年11月, 2022 ANS Winter Meeting and Technology Expo, 1 - 2, 英語Effects of Upper and Lower Pipe Ends on Interfacial and Wall Friction Factors under Counter-Current Flow Limitation in a Vertical Pipe[査読有り]
- 2022年11月, 2022 ANS Winter Meeting and Technology Expo, 1 - 2, 英語Effects of Liquid Viscosity on Friction Factors of Gas-Liquid Swirling Annular Flows in a Vertical Pipe[査読有り]研究論文（国際会議プロシーディングス）
- Begell House, 2022年09月, Multiphase Science and Technology, 34(4) (4), 41 - 55, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2022年08月, Multiphase Science and Technology, 34(4) (4), 1 - 13, 英語[査読有り]研究論文（学術雑誌）
- Abstract : The authors present the development of the lattice Boltzmann method (LBM) for application to gas-liquid flows, gas mixture, and gas-solid flows. The LBM scheme is constructed using a simple procedure that takes into account collision and propagation processes. The fluid density is calculated using the sum of the distribution function with respect to number of the discrete velocity. Since all distribution functions exist in the computational region with adequate boundary conditions, the LBM is effective for mass conservation and can be adjusted to the multiphase flow simulation. For the multiphase flow simulation, the color-gradient model is first proposed, and then the Shan-Chen and free energy models were implemented using the LBM. Then, the phase-field model, which is commonly used in LBMs, was investigated. The phase-field method easily captures the complicated interface and calculates multiphase flows with high density and viscosity ratios. For multicomponent flows, the LBM can calculate the Stefan--Maxwell equation. The present model calculates multicomponent phenomena, including uphill diffusion, without interpolating the distribution function. For particulate flows, we list the momentum exchange method and the immersed boundary -lattice Boltzmann method (IB-LBM) and explain the occurrence of boundary slip. Finally, we introduce models for wettability, including three phases: gas, liquid, and solid. Keywords : Lattice Boltzmann method, Multiphase flow, Multicomponent flow, Particulate flow, Wettability Cite this article : Seta, T.; Yamamoto, K.; Yoshino, M.; Takada, N.; Matsukuma, Y.; Yamamoto, K.; Hayashi, K.; Kimura, H.; Tomiyama, A. Lattice Boltzmann method for multiphase and multicomponent flows: A review. Muitiph. Sci. Technol. 2022, Vol. 34, in press. https://doi.org/10.1615/MultScienTechn.2022044039Begell House, 2022年07月, Multiphase Science and Technology, 34, in press., 英語, 国際誌, 国際共著していない[査読有り]研究論文（学術雑誌）
- 2022年07月, International Journal of Multiphase Flow, 153, 104134, 英語Fragmentation of Drops Falling through a Miscible Liquid with and without Drop-Drop Interactions[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2022年06月, World Neurosurgery, 162, 42 - 42, 英語, 国際誌[査読有り]研究論文（学術雑誌）
- 2022年05月, 2022 Japan-U.S. Seminar on Two-Phase Flow Dynamics, 1 - 6, 英語EFFECTS OF SWIRL INTENSITY ON INTERFACIAL AND WALL FRICTION FACTORS OF ANNULAR FLOWS IN A VERTICAL PIPE[査読有り]研究論文（国際会議プロシーディングス）
- Elsevier BV, 2022年05月, Computers & Mathematics with Applications, 114, 180 - 187, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2022年03月, Multiphase Science and Technology, 英語[査読有り]研究論文（学術雑誌）
- Numerical simulations of contaminated spherical drops falling through a stagnant liquid at low Reynolds numbers are carried out using the finite difference method. The numerical results are used to describe the behavior of the surfactant concentrations and to understand the surfactant effects on the fluid motions in detail. The predicted interfacial surfactant concentration, Γ, is almost zero for angles, θ, below a certain value (the stagnant-cap angle, θcap), whereas it steeply increases and reaches a large value for θ>θcap (the stagnant-cap region). The increase in the initial surfactant concentration, C0, in the drop enhances the adsorption from the drop to the interface, which results in the increase in Γ and the decrease in θcap. Peaks appear in the predicted Marangoni stresses around θcap, which causes similar peaks in the pressure distribution. The high-pressure spots prevent the fluid motion along the interface, which results in the formation of the stagnant-cap region and the attenuation of the tangential velocity in the continuous phase. The surfactant flux from the bulk to the interface decreases C in the vicinity of the interface for θ<θcap and the weak diffusion cannot compensate for the reduction in C by adsorption, which results in C at the interface smaller than C0. The pattern of the low C region is determined by the advection and does not smear out because of a small diffusive flux.MDPI AG, 2022年01月, Fluids, 7(2) (2), 55 - 55, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2022年01月, Multiphase Science and Technology, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2021年12月, Chemical Engineering Science, 245, 116852 - 116852, 英語[査読有り]研究論文（学術雑誌）
- 2021年10月, Nuclear Technology, 英語Void Fraction and Interfacial Friction in Vertical Circular Pipes with the Square Top End under Flooding Conditions[査読有り]
- The Japanese Society for Multiphase Flow, 2021年09月, 混相流, 35(4) (4), 533 - 542, 日本語[査読有り]研究論文（学術雑誌）
- Lattice Boltzmann simulations and a velocity measurement of flows in a cerebral aneurysm reconstructed from MRA (magnetic resonance angiography) images of an actual aneurysm were carried out and the numerical results obtained using the bounce-back schemes were compared with the experimental data to discuss the effects of the numerical treatment of the no-slip boundary condition of the complex boundary shape of the aneurysm on the predictions. The conclusions obtained are as follows: (1) measured data of the velocity in the aneurysm model useful for validation of numerical methods were obtained, (2) the numerical stability of the quadratic interpolated bounce-back scheme (QBB) in the flow simulation of the cerebral aneurysm is lower than those of the half-way bounce-back (HBB) and the linearly interpolated bounce-back (LBB) schemes, (3) the flow structures predicted using HBB and LBB are comparable and agree well with the experimental data, and (4) the fluctuations of the wall shear stress (WSS), i.e., the oscillatory shear index (OSI), can be well predicted even with the jaggy wall representation of HBB, whereas the magnitude of WSS predicted with HBB tends to be smaller than that with LBB.MDPI AG, 2021年09月, Fluids, 6(10) (10), 338 - 338, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2021年09月, International Journal of Multiphase Flow, 142, 103653 - 103653, 英語[査読有り]研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2021年08月, 混相流, 35(3) (3), 463 - 472, 日本語[査読有り]研究論文（学術雑誌）
- BACKGROUND: The thin-walled regions (TIWRs) of intracranial aneurysms have a high risk of rupture during surgical manipulation. They have been reported to be predicted by wall shear stress and pressure (PS) based on computational fluid dynamics analysis, although this remains controversial. In this study, we investigated whether the oscillatory shear index (OSI) can predict TIWRs. METHODS: Twenty-five unruptured aneurysms were retrospectively analyzed; the position and orientation of the computational fluid dynamics color maps were adjusted to match the intraoperative micrographs. The red area on the aneurysm wall was defined as TIWR, and if most of the regions on the color map corresponding to TIWR were OSI low (lower quartile range), time-averaged wall shear stress (TAWSS) high, or PS high (upper quartile range), each region was defined as a matched region and divided by the total number of TIWRs to calculate the match rate. In addition, the mean values of OSI, TAWSS, and PS corresponding to TIWRs were quantitatively compared with those in adjacent thick-walled regions. RESULTS: Among 27 TIWRs of 25 aneurysms, 23, 10, and 14 regions had low OSI, high TAWSS, and high PS regions (match rate: 85.2%, 37.0%, and 51.9%), respectively. Receiver operating characteristic curve analysis demonstrated that OSI was the most effective hemodynamic parameter (area under the curve, 0.881), followed by TAWSS (0.798). Multivariate analysis showed that OSI was a significant independent predictor of TIWRs (odds ratio, 18.30 [95% CI, 3.2800-102.00], P < 0.001). CONCLUSIONS: OSI may be a unique predictor for TIWRs. Low OSI strongly corresponds with TIWRs of intracranial aneurysms.2021年08月, World neurosurgery, 152, e377-e386, 英語, 国際誌研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2021年03月, 混相流, 35(1) (1), 85 - 92, 日本語[査読有り][招待有り]研究論文（学術雑誌）
- Elsevier BV, 2021年03月, Nuclear Engineering and Design, 373, 111020 - 111020, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2021年, Multiphase Science and Technology, 33(3) (3), 17 - 31, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2021年, Multiphase Science and Technology, 33(3) (3), 53 - 66, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2021年, Multiphase Science and Technology, 33(4) (4), 1 - 13, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2021年01月, Nuclear Engineering and Design, 371, 110951 - 110951, 英語[査読有り]研究論文（学術雑誌）
- Wiley, 2021年01月, Chemie Ingenieur Technik, 93(1-2) (1-2), 247 - 259, 英語[査読有り][招待有り]研究論文（学術雑誌）
- Begell House, 2021年, Multiphase Science and Technology, 33(2) (2), 87 - 101, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2021年, Multiphase Science and Technology, 33(2) (2), 69 - 85, 英語[査読有り]研究論文（学術雑誌）
- The Japanese Society for Multiphase Flow, 2020年11月, 混相流, 35(1) (1), 159 - 167, 日本語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2020年10月, International Journal of Multiphase Flow, 131, 103403 - 103403, 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2020年08月, International Journal of Multiphase Flow, 129, 103350 - 103350, 英語[査読有り]研究論文（学術雑誌）
- Abstract This paper presents an analysis of the performance of the Spatial Filter Velocimetry (SFV) technique applied to external flows across tube bundles to evaluate the benefits of using low particle concentration in a complex geometry. The experimental vector velocity fields were obtained in a tube bundle test section made of acrylic with 20 rows of 4 tubes of 20mm O.D. mounted in a normal triangular configuration and transversal pitch per diameter ratio of 1.25. The SFV results reliability was evaluated by comparing the volumetric flow rates estimated by the experimental vector velocity fields with the volumetric flow rates measured by the flowmeter, showing good agreement. Additionally, the experimental results are used to check the validity of numerical results obtained through CFD simulation. The experiments were conducted for Reynolds numbers of 909, 1842 and 3902. The velocity profiles obtained through the SFV technique agree with the CFD simulations results, which indicate consistent results for both analyzes. Therefore, the numerical and experimental results indicate that the SFV technique provides accurate velocity data and is suitable to be applied to complex geometries, especially in regions away of the test section wall due to the good visibility provided by low particle concentration.American Society of Mechanical Engineers, 2020年07月, ASME 2020 Fluids Engineering Division Summer Meeting (FEDSM2020), 1 - 12, 英語[査読有り]研究論文（国際会議プロシーディングス）
- 2020年07月, ASME 2020 Fluids Engineering Division Summer Meeting, 英語Numerical Simulations of Flows in Cerebral Aneurysms Using the Lattice Boltzmann Method[査読有り]研究論文（国際会議プロシーディングス）
- Elsevier BV, 2020年06月, International Journal of Multiphase Flow, 127, 103270 - 103270, 英語[査読有り]研究論文（学術雑誌）
- Springer Science and Business Media LLC, 2020年04月, Flow, Turbulence and Combustion, 105, 1277 - 1301, 英語[査読有り]研究論文（学術雑誌）
- Springer Science and Business Media LLC, 2020年04月, Journal of the Brazilian Society of Mechanical Sciences and Engineering, 42(4) (4), 英語[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2020年03月, International Journal of Multiphase Flow, 124, 103207 - 103207, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2020年, Multiphase Science and Technology, 32(4) (4), 325 - 340, 英語[査読有り]研究論文（学術雑誌）
- Begell House, 2020年, Multiphase Science and Technology, 32(3) (3), 221 - 236, 英語[査読有り]研究論文（学術雑誌）
- 2019年11月, 11th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF)Application of Spatiotemporal Filter Velocimetry to Evaluation of Drag Force Acting on Single Drop[査読有り]研究論文（国際会議プロシーディングス）
- 2019年11月, Nuclear Engineering and Design, 353, 110223, 英語[査読有り]研究論文（学術雑誌）
- 2019年09月, Multiphase Science and Technology, 31(3) (3), 215 - 234, 英語[査読有り]研究論文（学術雑誌）
- 2019年07月, AJK-Fluids, 英語Evaluation of Drag Force Acting on Single Contaminated Drop based on Velocity Distribution[査読有り]研究論文（国際会議プロシーディングス）
- Elsevier BV, 2019年06月, International Journal of Heat and Mass Transfer, 136, 521 - 530, 英語[査読有り]研究論文（学術雑誌）
- 2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 2 pages, 英語Viscous Stress and Pressure Acting on Interface of Contaminated Drop[査読有り]研究論文（国際会議プロシーディングス）
- 2019年05月, Computers & Mathematics with Applications, 78, 2746 - 2760, 英語[査読有り]研究論文（学術雑誌）
- 2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 2 pages, 英語Effects of Surfactant on Bubble Breakup at T-Junction and Characteristics of Taylor Flows in Mini and Micro Channels[査読有り]研究論文（国際会議プロシーディングス）
- Begell House, 2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 33(2) (2), 6 pages - 101, 英語[査読有り]研究論文（国際会議プロシーディングス）
- 2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 9 pages, 英語Combined Effects of Surfactant and Electrolyte on Mass Transfer from Single Carbon-Dioxide Bubbles in Vertical Pipes[査読有り]研究論文（国際会議プロシーディングス）
- 2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 9 pages, 英語Characterization of Velocity Field External to Tube Bundle Using Spatial Filter Velocimetry based on Variable Meshing Scheme[査読有り]研究論文（国際会議プロシーディングス）
- 2019年03月, Multiphase Science and Technology, 31(1) (1), 61 - 71, 英語[査読有り]研究論文（学術雑誌）
- 2019年02月, Internation Journal of Multiphase Flow, 111, 285 - 293, 英語[査読有り]研究論文（学術雑誌）
- 2019年02月, ISIJ International, 59(2) (2), 209 - 215, 英語[査読有り]研究論文（学術雑誌）
- OBJECTIVE: Cerebral aneurysm growth often precedes rupture. Definite contributors to aneurysm growth have not been determined even by means of recently developed commercially available computational fluid dynamics (CFD) software. We developed an original CFD tool that can analyze data from time-of-flight magnetic resonance angiography (TOF-MRA) before growth in the growing aneurysms and investigate possible factors for aneurysm growth in the near future. METHODS: We retrospectively reviewed unruptured aneurysms that were treated at our institute because of aneurysm growth (growing group) between April 2013 and March 2017. Stable aneurysms that had demonstrated no growth for more than 5 years were selected (stable group). TOF-MRA data of these aneurysms were retrospectively converted to 3-dimensional vessel geometric data; 3 hemodynamic indices including streamline, wall shear stress (WSS), and oscillatory shear index were calculated by our original CFD tool using the lattice Boltzmann method to quantitatively compare the 2 groups. RESULTS: Six growing aneurysms and 6 stable aneurysms were analyzed. Of the 6 growing aneurysms, WSS on the focal aneurysmal sac increased temporally in the vicinity of the constant low WSS area at the peak systolic phase. By contrast, WSS did not increase during any part of the cardiac cycle in 3 of the 6 stable aneurysms. The peak values of WSS were significantly different between the 2 groups. CONCLUSIONS: A focal increase in WSS in the peak systolic phase may be a risk factor for aneurysm enlargement in the near future.2019年02月, World neurosurgery, 122, e1439-e1448 - e1448, 英語, 国際誌[査読有り]研究論文（学術雑誌）
- OBJECTIVE: Thin-walled regions (TIWRs) within cerebral aneurysms have a high risk of rupture during surgical manipulation. Previous reports have demonstrated specific changes in the parameters of computational fluid dynamics in TIWRs; however, they have not been fully evaluated. We identified and investigated a novel parameter, wall shear stress vector cycle variation (WSSVV), with user-friendly software that could predict TIWRs. METHODS: Twelve unruptured cerebral aneurysms were analyzed. TIWRs were defined as reddish areas compared with the normal-colored parent artery on intraoperative views. The position and orientation of these clinical images were adjusted to match the WSSVV color maps. TIWRs and thick-walled regions (TKWRs) were marked and compared with the corresponding regions on WSSVV maps. The default images obtained from WSSVV imaging required appropriate maximum color bar value (MCBV) adjustment for predicting TIWRs. Sensitivity and specificity analyses were performed by changing the MCBV from 300 to 700 at intervals of 100. With the optimal MCBV, the WSSVV values were quantitatively compared. RESULTS: All of the selected 18 TIWRs and 16 TKWRs corresponded to low- and high-value regions of the WSSVV color maps at the adjusted MCBV, respectively. The mean optimal MCBV was 483.3 ± 167.50 (range, 300-700). According to receiver operating characteristic analysis, the best MCBV for predicting TIWRs was 500 (highest sensitivity, 0.89; specificity, 0.94). Under this condition, the quantitative values of the computational fluid dynamics color maps for TIWRs and TKWRs were significantly different (P < 0.01). CONCLUSIONS: Low WSSVV values may indicate TIWRs within cerebral aneurysms.2019年01月, World neurosurgery, 121, e287-e295 - e295, 英語, 国際誌[査読有り]研究論文（学術雑誌）
- Elsevier BV, 2018年12月, International Journal of Heat and Fluid Flow, 74, 28 - 35, 英語[査読有り]研究論文（学術雑誌）
- Aspect ratios of ellipsoidal bubbles in clean and contaminated systems were measured to investigate the effects of surfactant on the bubble aspect ratio. The bubble diameter ranged from 0.80 to 5.8 mm. Air was used for the gas phase and the glycerol-water solutions were used for the liquid phase. The experimental ranges of the Morton number, M, the bubble Reynolds number, Re, the Eötvös number, Eo, the Weber number, We, and the Tadaki number, Ta, were as follows: 10−7 < M < 10−3, 3.6 × 10−1 < Re < 1.7 × 102, 1.1 × 10−1 < Eo < 5.9, 5.0 × 10−3 < We < 3.7 and 5.3 × 10−2 < Ta < 3.6. Triton X-100 and 1-octanol were used for surfactant. The terminal velocities of bubbles were confirmed to be independent of the surfactant concentration, and therefore, the bubbles were fully-contaminated from the point of view of the terminal velocity. By making use of the aspect ratio data, the applicability of available aspect ratio correlations was examined and an empirical correlation was proposed. The applicability of the proposed correlation was also examined for sodium dodecyl sulphate (SDS) and 1-decanol. As a result, the following conclusions were obtained: (1) the aspect ratio database, which is of use for validating and developing shape correlations, was obtained, (2) the available aspect ratio correlations are not applicable to the present data for fully-contaminated ellipsoidal bubbles, (3) the aspect ratios are well correlated in terms of a combination of the Eötvös and bubble Reynolds numbers, which proves that all the relevant forces, i.e. the viscous, inertial, buoyant and surface tension forces, should be taken into account to correlate bubble shapes, and (4) the empirical correlation of the aspect ratio developed by using the Triton X-100 and 1-octanol data also gives good evaluations for bubbles contaminated with SDS and 1-decanol in spite of largely different adsorption-desorption characteristics.Elsevier Inc., 2018年09月, Experimental Thermal and Fluid Science, 96, 460 - 469, 英語[査読有り]研究論文（学術雑誌）
- Elsevier, 2018年09月, Experimental Thermal and Fluid Science, 96, 397 - 405, 英語Evaluation of Adsorption of Surfactant at a Moving Interface of a Single Spherical Drop[査読有り]研究論文（学術雑誌）
- 2018年08月, Nuclear Science and Engineering, 193, 147 - 159, 英語[査読有り][招待有り]研究論文（学術雑誌）
- 2018年07月, 19th Symposium on Applications of Laser and Imaging Techniques to Fluid Mechanics, Lisbon 2016, 英語Diagnostics of flow about a contaminated single drop using spatiotemporal filter velocimetry[査読有り]研究論文（国際会議プロシーディングス）
- Experiments on air–water two-phase swirling annular flows in a vertical pipe of 40 mm diameter were carried out at atmospheric pressure and room temperature to investigate interfacial and wall friction factors, fi and fw. The friction factors were evaluated using measured pressure drops and void fractions. Measurements of liquid film thickness and flow observation were also conducted. The ranges of the gas and liquid volume fluxes, JG and JL, were 12.5 ≤ JG ≤ 20.0 m/s and 0.03 ≤ JL ≤ 0.11 m/s, respectively. The main conclusions obtained are as follows: (1) fi and fw in swirling annular flows are several times larger than those in non-swirling flows, (2) fi is well correlated in terms of the liquid volume fraction and the gas Reynolds number, ReG, (3) ReG and the liquid Reynolds number, ReL, are required for correlating fw, and (4) the liquid film thicknesses in two-phase swirling flows in a one-fifth model of a BWR separator are well predicted using the two-fluid model and the correlations of fi and fw developed based on the experimental data.Elsevier Ltd, 2018年04月, Nuclear Engineering and Design, 330, 97 - 105, 英語[査読有り]研究論文（学術雑誌）
- 2018年04月, 8th European-Japanese Twp-Phase Flow Group Meeting, EJTPFGM-8, 1, 英語Wall and Interfacial Shear Stress of Swirling Annular Flow[査読有り]研究論文（国際会議プロシーディングス）
- 2018年04月, 8th European-Japanese Twp-Phase Flow Group Meeting, EJTPFGM-8, 1, 英語Numerical Simulation of a Drop in Contaminated Systems from low concentration to CMC[査読有り]研究論文（国際会議プロシーディングス）
- We analytically and numerically investigate the boundary slip, including the velocity slip and the temperature jump, in immersed boundary-thermal lattice Boltzmann methods (IB-TLBMs) with the two-relaxation-time collision operator. We derive the theoretical equation for the relaxation parameters considering the effect of the advection velocity on the temperature jump of the IB-TLBMs. The analytical and numerical solutions demonstrate that the proposed iterative correction methods without the computational cost of the sparse matrix solver reduce the boundary slip and boundary-value deviation as effectively as the implicit correction method for any relaxation time. Because the commonly used multi-direct forcing method does not consider the contributions of the body force to the momentum flux, it cannot completely eliminate the boundary slip because of the numerical instability for a long relaxation time. Both types of proposed iterative correction methods are more numerically stable than the implicit correction method. In simulations of flow past a circular cylinder and of natural convection, the present iterative correction methods yield adequate results without the errors of the velocity slip, the temperature jump, and the boundary-value deviation for any relaxation time parameters and for any number of Lagrangian points per length. The combination of the present methods and the two-relaxation-time collision operator is suitable for simulating fluid flow with thermal convection in the multiblock method in which the relaxation time increases in inverse proportion to the grid size.John Wiley and Sons Ltd, 2018年03月, International Journal for Numerical Methods in Fluids, 86(7) (7), 454 - 490, 英語[査読有り]研究論文（学術雑誌）
- A semi-empirical correlation for CCFL (counter-current flow limitation) in vertical pipes was derived from one-dimensional momentum equations. Available correlations such as the empirical correlations proposed by Wallis (1969) and Zapke and Kröger (1996) can be deduced from the derived fundamental functional form of the semi-empirical correlation in limiting flow conditions. Comparisons between the semi-empirical correlation with available experimental data of CCFL taking place at the sharp-edged lower end of a vertical pipe showed that the correlation is applicable for various fluid properties and pipe diameters. The fundamental functional form of the correlation was also transformed so as to express the characteristics of CCFL occurring at the sharp-edged upper end of a vertical pipe. The present fundamental functional form of CCFL correlation is useful not only to express CCFL data but also to understand how relevant parameters play their roles in the CCFL characteristics.Elsevier Ltd, 2018年03月, Nuclear Engineering and Design, 328, 182 - 187, 英語[査読有り]研究論文（学術雑誌）
- 2018年02月, International Journal of Multiphase Flow, 99, 86 - 93, 英語[査読有り]研究論文（学術雑誌）
- Experiments on aeration from a bubble diffuser pipe having five aeration holes were carried out to in vestigate the effects of the azimuthal angle of the holes and liquid viscosity on flows inside and outside the pipe. The azimuthal angle was varied by rotating the pipe. When the azimuthal angle of the hole exceeded a certain angle, the liquid height inside the pipe was fixed just below the holes at any liquid viscosity, which resulted in preventing of the form ation of liquid slugs inside the pipe. Hence, uniform aeration was realized just by rotating the pipe. The downstream gas flow rates from holes were slightly higher than those upstream, which is due to the pressure recovery in the direction of the stream.Society of Chemical Engineers, Japan, 2018年, 化学工学論文集, 44(1) (1), 59 - 66, 日本語[査読有り]研究論文（学術雑誌）
- Mass transfer rates, k(L), of single carbon-dioxide (CO2) bubbles rising through vertical pipes filled with electrolyte aqueous solutions were measured to investigate the effects of the presence of electrolyte on k(L) Sodium chloride (NaCI) was used for electrolyte and its concentration was varied from 3.5 to 14 wt.%. The pipe diameters, D, were 12.5, 18.2 and 25.0 mm and the bubble diameter, d, ranged from 5 to 25 mm. The diameter ratio, lambda(= d/D), was varied from 0.2 to 1.7, to cover various bubble shapes, i.e. ellipsoidal, cap, semi-Taylor bubbles and Taylor bubbles. The conclusions obtained are as follows: (1) k(L) in the NaCI aqueous solutions decreases with increasing the NaCI concentration mainly due to the reduction of the diffusion coefficient of CO2 in the liquid phase, (2) the Sherwood numbers, Sh, of Taylor bubbles of L/D > 1 in the NaCI aqueous solutions can be well evaluated using the available Sherwood number correlation for clean Taylor bubbles, where L is the length of a Taylor bubble, (3) Sh of ellipsoidal, cap, semi-Taylor bubbles and Taylor bubbles of LID <= 1 in the clean water and NaCI aqueous solutions can be well correlated in terms of the Peclet number and the dimensionless bubble diameter d* = d/d(T), where d(T) is the bubble diameter at the transition from the ellipsoidal-cap bubble regime to the semi-Taylor bubble regime, and (4) the Sh correlations give good predictions for long-term bubble dissolution processes in NaCI aqueous solutions. (c) 2017 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2017年12月, INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER, 115, 663 - 671, 英語[査読有り]研究論文（学術雑誌）
- Spatiotemporal filter velocimetry (SFV) was extended to Lagrangian measurements with boundary-fitted measurement areas, and was applied to flows about single spherical drops of glycerol-water solution falling in stagnant silicon oil under clean and contaminated conditions to examine its applicability to the estimation of the Marangoni stress and surfactant concentration at a moving interface. Effects of bulk concentration of surfactant on the velocity field, the Marangoni stress and the surface concentration of surfactant were discussed from the measured data. As a result, we confirmed that accurate velocity distribution in the vicinity of the interface measured by SFV enables us to evaluate interfacial velocity and interfacial shear stresses and to estimate the Marangoni stress, interfacial tension and surfactant concentration at the interface with the assumption of negligible surface viscosity. The flow inside the drop and the interfacial velocity become weak due to the Marangoni stress caused by the gradient of surfactant concentration at the interface as the bulk concentration of surfactant increases. These results demonstrate that SFV is of great use in experimental analysis of adsorption and desorption kinetics at a moving interface.Elsevier Ltd, 2017年12月, International Journal of Multiphase Flow, 97, 157 - 167, 英語[査読有り]研究論文（学術雑誌）
- 2017年12月, 10th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF2017), 英語Evaluation of Adsorption of Surfactant at a Drop Interface using Spatiotemporal Filter Velocimetry[査読有り]研究論文（国際会議プロシーディングス）
- Experiments on aeration from a bubble diffuser pipe having five aeration holes are carried out to investigate effects of the azimuthal angle of the holes on flows inside and outside the pipe. The azimuthal angle is varied by rotating the pipe. When the hole angle becomes larger than a certain angle, the liquid height inside the pipe is fixed just below the holes, which results in the prevention of slugging inside the pipe. The effects of hole angle on bubble generation mode and bubble size are small except for the large hole angles, at which generated bubbles are to break up due to an interaction between the bubbles and the pipe wall. Hence uniform aeration is easily realized just by rotating the pipe to some extent. Uniform aeration is also realized with a longer diffuser pipe having ten aeration holes by rotating the pipe. The power in aeration is evaluated from the pressure difference and the total gas inflow. Although the decrease in the hole diameter is also effective to realize the uniform aeration, the increase in the azimuthal angle is superior to the former from the point of view of energy saving.ELSEVIER SCIENCE INC, 2017年12月, EXPERIMENTAL THERMAL AND FLUID SCIENCE, 89, 90 - 97, 英語[査読有り]研究論文（学術雑誌）
- 2017年11月, ANS Winter Meeting and Nuclear Technology Expo, 2 pages, 英語Spatiotemporal Flow Structure of Counter-Current Two-Phase Annular Flows[査読有り]研究論文（国際会議プロシーディングス）
- Lift coefficients, C-L, of single bubbles in linear shear flows are measured to investigate effects of the bubble shape, the liquid velocity gradient and the fluid property on C-L. The range of the Morton number, M, tested is from logM = -6.6 to -3.2. The shapes of bubbles are spherical and ellipsoidal. A correlation of bubble aspect ratio for single bubbles in infinite stagnant liquids proposed in our previous study can give good evaluations for bubbles in the linear shear flows. The C-L of spherical bubbles at low bubble Reynolds numbers, Re, depend on the dimensionless shear rate Sr and Re and decrease with increasing Re. These characteristics agree with the Legendre-Magnaudet correlation. The use of a single dimensionless group such as Re, the Eotvos number, the Weber number and the Capillary number cannot correlate C-L of non-spherical bubbles. The trend of the critical Re for the reversal of the sign of C-L is the same as that for the onset of oscillation of bubble motion, which supports the mechanism proposed by Adoua et al., at least within the range of -6.6 <= logM <= -3.2. An experimental database of C-L is provided for validation of available C-L models and CFD. (C) 2017 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2017年11月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 96, 113 - 122, 英語[査読有り]研究論文（学術雑誌）
- 微小気泡の流動を利用した球底へのキプリス幼生の着生阻害海水を水槽に連続供給しながら，水槽内に設置した丸底フラスコの下方へ微小気泡群を30日間継続供給した結果，着生し成長したフジツボ固体数は，35～45個体であった．微小気泡群を供給しなかった場合，個体の総数は664個体であった．丸底フラスコ下半球付近では，気泡径51～100 μmの微小気泡群の頻度が最も高く，32.6～40.2%であった．気泡径151～300 μmの微小気泡群の頻度は，フラスコ底部から95 mm下方位置で3.4%であったが，フラスコ下半球近傍では11.3～17.0%へ増加した．これらの結果は，気泡径51～300 μmの微小気泡群に，キプリス幼生の着生を阻害する機能があることを示唆している．NP2モデルを使用して微小気泡群の流動計算を行った結果，気泡径75～175 μmの微小気泡群が丸底フラスコ近傍に供給されやすいことが確認できた．キプリス日本海水学会, 2017年10月, 日本海水学会誌, 71(5) (5), 315 - 320, 日本語[査読有り]研究論文（学術雑誌）
- The swirler and the 1st pick-off ring (POR) of a gas-liquid separator are improved to reduce the pressure drop, DPT, in the separator while keeping the flow rate of separated liquid, WST, high. Two swirler designs are presented, i.e. V8 (swirler with eight vanes) and V6 (swirler with six vanes), whose diameters of the hubs supporting the vanes are smaller than that of the normal swirler, N8, and the nose shape of vanes are coned. The swirlers are installed into a one-fifth scale separator with three PORs. The gap width of the 1st POR is reduced to the value of the liquid film thickness in the barrel of the separator and the rear shape of the 1st POR is smoothed. The performance of the downscaled separator with the proposed swirlers and the 1st POR are evaluated under the nominal operating condition of the Hyper BWR and a higher quality condition. As a result, the following conclusions are obtained in the present experimental range: (1) DPT in the separator is reduced by 23% using V8, and V8 gives the same WST as N8, confirming that the decrease in the hub diameter and the coned-nose shape of vanes are effective, (2) V6 can realize further reduction of DPT (38%) without deteriorating the liquid-separation rate, (3) the reduced gap width of the 1st POR and a tapered shape of the downstream side of the 1st POR are also effective to decrease the pressure drop, i.e. the improved POR decreases DPT by 42% and 54% by combining with V8 and V6, respectively, and (4) the improved separators also show good performances in DPT and WST at the higher quality condition. (C) 2017 Elsevier B.V. All rights reserved.ELSEVIER SCIENCE SA, 2017年10月, NUCLEAR ENGINEERING AND DESIGN, 322, 360 - 367, 英語[査読有り]研究論文（学術雑誌）
- 2017年09月, 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics(NURETH-17), 14 pages, 英語On CCFL in Collider Test Facility: CCFL Results at PWR Hot-Leg Pipe Geometry with 190 mm Diameter[査読有り]研究論文（国際会議プロシーディングス）
- 日本混相流学会, 2017年06月, 混相流, 31(1) (1), 37 - 46, 日本語
One-region computations with the annular flow model were done for counter-current flow limitation (CCFL) at the sharp-edged lower end of vertical pipes to evaluate effects of diameters and fluid properties on CCFL characteristics. CCFL characteristics computed with several correlations for interfacial friction factors were compared with CCFL data and the correlation proposed by Bharathan et al. (which is a function of void fraction) was selected. The adjustment factors were obtained to give good agreement between CCFL characteristics computed with the correlation by Bharathan et al. and CCFL data and were correlated with the viscosity ratio of gas and liquid phases and the dimensionless diameter. By using the correlation for interfacial friction factors by Bharathan et al. modified with the viscosity ratio of gas and liquid phases and the dimensionless diameter, effects of diameters and fluid properties on CCFL characteristics were computed.
[査読有り]研究論文（学術雑誌） - 日本混相流学会, 2017年06月, 混相流, 31(2) (2), 152 - 161, 日本語
To evaluate effects of the diameter on countercurrent flow limitation (CCFL) in vertical pipes, in our previous study, we classified CCFL into CCFL-L at the sharp-edged lower end, CCFL-U at the sharp-edged upper end, and CCFL-P inside the vertical pipe with round-edged upper and lower ends, and we found that the characteristic length in the Wallis parameter, w = D(1-β)Lβ (where D and L respectively denote the diameter and the Laplace capillary length), is respectively β = 0, 1 and β ≈ 0.5 for CCFL-L, CCFL-U and CCFL-P. In this study, we evaluated effects of fluid properties on CCFL in vertical pipes by using existing CCFL data. The most effective fluid properties is the liquid viscosity and we selected the viscosity ratio of gas and liquid (μG/μL) as a dimensionless parameter. From effects of μG/μL on the slope m and constant C in the Wallis correlation, we classified three regions of large positive interrelations of m and C with μG/μL for small μL (room-temperature to high-temperature water), small interrelations of m and C with μG/μL for medium μL (low-concentration glycerol water solution), and the large negative interrelation of m with μG/μL for large μL (high-concentration glycerol water solution). We obtained exponent functions of μG/μL for m and C in the case of the small μL (room-temperature to high-temperature water) for CCFL-L and CCFL-U except CCFL-P without high-temperature CCFL data.
[査読有り]研究論文（学術雑誌） - 2017年06月, 9th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamic (ExHFT-9), 9 pages, 英語Terminal Velocities of Single Bubble Rising through Stagnant Liquids and Slurry in a Hele-Shaw Cell[査読有り]研究論文（国際会議プロシーディングス）
- 2017年06月, 9th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamic (ExHFT-9), 8 pages, 英語Pressure Drop in Air-Water Two-Phase Flows in Horizontal U-Bends[査読有り]研究論文（国際会議プロシーディングス）
- 2017年06月, 9th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamic (ExHFT-9), 10 pages, 英語Evaluation of Adsorption of Surfactant at a Moving Interface of a Single Spherical Drop[査読有り]研究論文（国際会議プロシーディングス）
- Mass transfer rates, k(L), of single carbon dioxide bubbles contaminated with alcohols in vertical pipes of 12.5, 18.2 and 25.0 mm diameters were measured to investigate effects of the carbon chain length of the alcohols on kL. Four straight-chain alcohols, i.e. 1-pentanol, 1-heptanol, 1-octanol and 1-decanol, were used. A wide range of bubble diameter, d, covered ellipsoidal and Taylor bubbles. Effects of the alcohol concentration, C-sol, on k(L) were also investigated to make clear the values of C-sol corresponding to fully contaminated conditions. As a result, the following conclusions were obtained. (1) The kL of bubbles contaminated with alcohols decrease with increasing C-sol, and they become constant at high C-sol, i.e. bubbles in the latter C-sol ranges are fully-contaminated from the point of view of mass transfer. On the other hand, the rising velocities decrease with increasing C-sol due to the reduction of surface tension even when k(L) are independent of Cob (2) The k(L). of ellipsoidal bubbles fully-contaminated with alcohols decrease with increasing the carbon chain length, whereas the reduction rate of kL becomes smaller and smaller as the carbon chain length increases. The k(L). of fully-contaminated ellipsoidal bubbles can be expressed by accounting for the adsorption constant of alcohol in an available kL correlation. (3) The increase in the carbon chain length also decreases k(L) of fully-contaminated Taylor bubbles, whereas the effect of carbon chain length on k(L), becomes smaller with increasing d. (C) 2017 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2017年05月, INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER, 108(Part B) (Part B), 1991 - 2001, 英語[査読有り]研究論文（学術雑誌）
- 2017年05月, Third International Symposium on Multiscale Multiphase Process Engineering (MMPE2017), 6 pages, 英語Effects of electrolyte on mass transfer from single carbon-dioxide bubbles in vertical pipes[査読有り]研究論文（国際会議プロシーディングス）
- Experimental Investigation of Temperature Distribution in Condensing Steam-Air Flow in a Circular TubeExperiments on a steam-air mixture flowing in a vertical tube were carried out to understand a temperature distribution in the tube and to accumulate a database of condensation heat transfer of the steam-air mixture. The condensation tube used was 49.5 mm in diameter and 510 mm long. The ratio of inlet mass fluxes of air and steam was 3:2. They were mixed at a mixing section and flowed into the condensation tube at fully-developed conditions through a straight tube. The water flowing in the outside of the condensation tube cooled the steam-air mixture. The pressure was 0.125 MPa at the inlet of the tube. The experiments were carried out at two different inlet air temperatures and at three different Reynolds numbers. The temperature was measured at 11 points in the radial direction and 7 levels in the vertical direction in the condensation tube. The local heat fluxes were evaluated from temperature distributions measured in the tube wall. Experimental data showed that the heat transfer coefficient h decreased in the flow direction except when flow at the inlet was in the saturated condition. Vierow & Schrock's correlation well evaluated the decrease trend of hc in the flow direction, however, it failed in evaluating dependence of h on the Reynolds number Rem. Siddique's correlation takes the radial distribution of mass fraction of non-condensable gas into account, which well evaluated the dependence of h on Re . This result implies that the hc correlations should take into account axial and radial distributions of the mass fraction of non-condensable gas as well as liquid film thickness. c c c m2017年04月, 2017 International Congress on Advances in Nuclear Power Plants (ICAPP 2017), 英語[査読有り]研究論文（国際会議プロシーディングス）
- Experiments on the total gas holdup, alpha(G), in air-water cylindrical bubble columns were carried out to investigate effects of the column diameter, D-H, and the initial liquid height, H-0, on alpha(G). Ranges of D-H and H-0 were 160 <= D-H <= 2000 mm and 400 <= H-0 <= 4000 mm, respectively. The superficial gas velocity, J(G), was varied from 0.025 to 0.35 m/s. The characteristics of gas holdup showed that all the flows in the present experiments were pure heterogeneous. The following conclusions were obtained for alpha(G) in air-water bubble columns: (1) the effects of D-H and H-0 on alpha(G) are negligible when scaling up from small to large bubble columns, provided that alpha(G) in the small columns are obtained for D-H >= 200 mm and Ho 2200 mm. The height-to-diameter ratio is useless in evaluation of the critical height, above which aG does not depend on H-0, (2) for the above ranges of D-H and H-0, Akita-Yoshida's and Koide's correlations can give good evaluations of alpha(G) for a wide range of J(G) by tuning the model constants, (3) for D-H < 200 mm, the decrease in D-H increases the population of large bubbles, which results in the decrease in alpha(G), and (4) for H-0 less than or similar to 2200 mm and D-H >= 200 mm, alpha(G) at a constant J(G) decreases with increasing H-0 and approaches an asymptotic value, and the Froude number using J(G) and H-0 as the characteristic scales well correlates alpha(G) in this regime. (C) 2016 Elsevier Inc. All rights reserved.ELSEVIER SCIENCE INC, 2017年04月, EXPERIMENTAL THERMAL AND FLUID SCIENCE, 82, 359 - 366, 英語[査読有り]研究論文（学術雑誌）
- 鉛直管内での気液対向流に対する界面摩擦係数We have been experimentally and analytically evaluated countercurrent flow limitation (CCFL), which is defined by the relationship between volumetric fluxes of updraft gas and falling liquid (J and J ) under flooding conditions, in the hot legs, the pressurizer surge line and the heat transfer tubes of a steam generator, and derived Wallis-type CCFL correlations for them. To generalize a prediction method for CCFL in vertical pipes, in this study, we did one-region computations for CCFL inside vertical pipes (CCFL-P) based on the annular flow model with several correlations for the interfacial friction factors and compared the computed results with air-water CCFL data. From the comparison, we selected the correlation proposed by Bharathan et al. (1978), which is a function of the void fraction. Then we obtained adjustment factors for the diameters of D = 18, 25.4, 31.8 and 50.8 mm by fitting the results computed with the correlation by Bharathan et al. (1978) to the air-water CCFL data and derived an adjustment function with the dimensionless diameter (D ). The results showed that the interfacial friction factors should be large for a large diameter pipe. It is difficult, however, to evaluate effects of fluid properties on the interfacial friction factors due to lack of steam-water CCFL data. G L *日本混相流学会, 2017年03月, 混相流, 31(1) (1), 37 - 46, 日本語[査読有り]研究論文（学術雑誌）
- Lift coefficients, CL, of single bubbles in linear shear flows are measured to investigate effects of the bubble shape, the liquid velocity gradient and the fluid property on CL. The range of the Morton number, M, tested is from logM = − 6.6 to − 3.2. The shapes of bubbles are spherical and ellipsoidal. A correlation of bubble aspect ratio for single bubbles in infinite stagnant liquids proposed in our previous study can give good evaluations for bubbles in the linear shear flows. The CL of spherical bubbles at low bubble Reynolds numbers, Re, depend on the dimensionless shear rate Sr and Re and decrease with increasing Re. These characteristics agree with the Legendre-Magnaudet correlation. The use of a single dimensionless group such as Re, the Eötvös number, the Weber number and the Capillary number cannot correlate CL of non-spherical bubbles. The trend of the critical Re for the reversal of the sign of CL is the same as that for the onset of oscillation of bubble motion, which supports the mechanism proposed by Adoua et al., at least within the range of −6.6 ≤ logM ≤ −3.2. An experimental database of CL is provided for validation of available CL models and CFD.Elsevier Ltd, 2017年, International Journal of Multiphase Flow, 96, 113 - 122, 英語[査読有り]研究論文（学術雑誌）
- 日本混相流学会, 2016年12月, 混相流, 30(4) (4), 392 - 401, 日本語
Characteristics of counter-current flow limitation (CCFL) inside vertical pipes with sharp-edged upper end (CCFL-P/SE) and with round-edged upper and lower ends (CCFL-P) were evaluated to clarify the dependence of the limiting locations. The CCFL is defined by the relationship between the gas and liquid volumetric fluxes (JG and JL) under a quasi-steady state after onset of flooding. CCFL correlations were constructed for CCFL-P/SE and CCFL-P based on the previous experimental data. The correlations were compared to the characteristics of CCFL at the vertical pipes of the sharp-edged upper end (CCFL-U) and the sharp-edged lower end (CCFL-L). CCFL-P has the highest liquid volumetric flux when CCFL occurs, and the lowest is CCFL-L. For the Wallis correlation, selection of the characteristic length of w = D(1-β)Lβ (where D and L respectively denote the diameter and the Laplace capillary length) is quite important. Our results have indicated that CCFL-P/SE and CCFL-P are also expressed by β ≈ 0.5, CCFL-U is expressed by β = 1, and CCFL-L is expressed by β = 0.
[査読有り]研究論文（学術雑誌） - Experiments on air-water two-phase swirling flows in a one-fifth scale model of a steam separator with three pick-off rings (PORs) are carried out to investigate characteristics of two-phase swirling flows and separator performance. In addition to flow observation, liquid film thickness, flow rates of separated liquid and pressure drops are measured. The ranges of the gas and liquid volume fluxes, J(G) and J(L), tested are 12.0 <= J(G) <= 17.8 m/s and 0.05 <= J(L) <= 0.11 m/s, respectively. The main conclusions obtained are as follows: (1) the liquid film thicknesses at the PORs are smaller than the gap widths of the PORs under the nominal operating condition, which results in a high liquid-separation rate, (2) the increase in J(G) decreases the flow rate of unseparated liquid, whereas the flow rate of unseparated liquid is independent of J(L) except for low J(G) conditions, under which the liquid film thicknesses at the 2nd and 3rd PORs tend to be larger than the gap widths of PORs, (3) the presence of the 2nd and 3rd PORs realizes high liquid-separation rates by capturing the liquid unseparated at the 1st POR and the pressure drops at the 2nd and 3rd PORs are much smaller than that at the 1st POR, and (4) the pressure drops at the swirler and the 1st POR, which are the main source of the pressure drop of the separator under the nominal operating condition, increase with J(G), whereas J(L) increases only the latter. (C) 2016 Elsevier B.V. All rights reserved.ELSEVIER SCIENCE SA, 2016年11月, NUCLEAR ENGINEERING AND DESIGN, 308, 205 - 213, 英語[査読有り]研究論文（学術雑誌）
- 2016年11月, ANS Conference 2016 Winter Meeting and Nuclear Technology Expo, 3 pages, 英語A Semi-Empirical Correlation for Counter-Current Flow Limitation in a Vertical Pipe[査読有り]研究論文（国際会議プロシーディングス）
- Total gas holdups alpha(G) in a cylindrical slurry bubble column were measured at various values of the superficial gas velocity J(G), the mean particle volumetric concentration C-S and the initial slurry height H-0 to investigate their effects on alpha(G). The column diameter D-H and height were 200 mm and 2000 mm, respectively. The gas, liquid and solid phases were air, water and hydrophilic silica particles of 100 mu m in mean diameter, respectively. Experimental conditions were 0.025 <= J(G) <= 0.40 m/s, 0 <= C-S <= 0.50 and 1.5 <= H-0* <= 5.0, where H-0* = H-0/D-H. The conclusions obtained are as follows: (1) alpha(G) decreases with increasing H0 and becomes independent of H-0 for H-0* > 4 at low J(G), whereas it depends on H-0 at higher J(G) even for H-0* > 4, (2) the increase in C-S decreases alpha(G) up to C-S similar to 0.40, whereas alpha(G) becomes independent of C-S at larger C-S, and (3) alpha(G) in the slurry bubble column at various C-S, H-0 and J(G) are well correlated in terms of the Froude number Fr-H using H-0 as a characteristic length.SOC CHEMICAL ENG JAPAN, 2016年09月, JOURNAL OF CHEMICAL ENGINEERING OF JAPAN, 49(9) (9), 824 - 830, 英語[査読有り]研究論文（学術雑誌）
- 2016年05月, 9th International Conference on Multiphase Flow, No. 262, 6 pages, 英語Numerical Simulation of Surfactant Concentration at an Interface of Single Spherical Drop[査読有り]研究論文（国際会議プロシーディングス）
- 2016年05月, 9th International Conference on Multiphase Flow, No. 245, 6 pages, 英語Mass Transfer from Single Carbon Dioxide Bubbles in Water-Alchol Solutions in a Verical Pipe[査読有り]研究論文（国際会議プロシーディングス）
- 2016年05月, 9th International Conference on Multiphase Flow, No. 262, 6 pages, 英語Experimental Evaluation of Surfacant Concentration at an Interface of Single Spherical Drop[査読有り]研究論文（国際会議プロシーディングス）
- 2016年05月, 9th International Conference on Multiphase Flow, No. 168, 6 pages, 英語Effects of Azimuthal Angle of Aeration Holes on Flows inside and outside a Bubble Diffuser Pipe[査読有り]研究論文（国際会議プロシーディングス）
- In this study, we measured counter-current flow limitation (CCFL) characteristics in an inverted U-tube (18.4 mm diameter and 1.0 m straight-part length) simulating steam generator (SG) U-tubes under conditions of steam condensation at pressures of 0.1-0.14 MPa. Differential pressure P between the top of the inverted U-tube and the lower tank was measured, and the flow patterns wave estimated by comparing the waveforms of P with those in air-water experiments. As a result, we classified the flow patterns under CCFL conditions into CCFL-P, CCFL-L and CCFL-T. The falling water flow rate under CCFL conditions slightly increased as the pressure increased and the cooling water temperature decreased (subcooling of cooling water increased). In the case of CCFL-L, CCFL characteristics in the inverted U-tube were between those in air-water and saturated steam-water experiments at 0.1 MPa. Furthermore, we derived a Wallis type CCFL correlation and its uncertainty from CCFL data, including previously measured data, i.e., J*(1/2)(G) + 0.88J(L)*(1/2) = 0.76 +/- 0.05.TAYLOR & FRANCIS LTD, 2016年04月, JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY, 53(4) (4), 486 - 495, 英語[査読有り]研究論文（学術雑誌）
- Effects of the initial liquid height, H-0, of air-water bubble columns on the total gas holdup, alpha(G), were investigated in this study. Systematic databases of alpha(G) in a rectangular and a cylindrical column of the hydraulic diameter DH of 200 mm were obtained. An image processing method was applied to highspeed video images of the liquid height to obtain accurate gas holdup data. Ranges of the superficial gas velocity, J(G), and the dimensionless liquid height, H-0* (= H-0/D-H), were 0.025 <= J(G) <= 0.40 m/s and 1.5 <= H-0(*) <= 5.0, respectively. The bubbly flows observed in these ranges could be classified into either the heterogeneous bubbly flow consisting of bubbles much smaller than D-H or that with huge bubbles of the column-width scale. The main parameter governing the flow regime transition was J(G). The gradient, d(alpha G)/d(JG), of alpha(G) with respect to J(G) was of use in flow regime identification. The increase in H-0* decreased alpha(G) because a long bubble residence time at a high lit) leads to the increase in the mean bubble diameter due to bubble coalescence. The Froude number using H-0. as a characteristic length well correlates alpha(G) at various H-0 An empirical correlation of alpha(G) in terms of the Froude number was then proposed. Comparing the correlation with the alpha(G) data showed that the correlation can give good evaluations of alpha(G) in the rectangular and cylindrical columns by tuning model parameters. (C) 2015 Elsevier Inc. All rights reserved.ELSEVIER SCIENCE INC, 2016年04月, Experimental Thermal and Fluid Science, 72, 67 - 74, 英語[査読有り]研究論文（学術雑誌）
- Aspect ratios E of ellipsoidal bubbles in infinite stagnant clean liquids are measured for log M = -6.6, -5.5, -4.9 and -3.9, where M is the Morton number. An empirical correlation of E applicable to a wide range of the Morton number is proposed by making use of the present data and Sugihara's data at log M = -11 (2007). The aspect ratio in this correlation is expressed in terms of the combination of the Eotvos number and the bubble Reynolds number to account for the effects of the inertial, viscous, buoyant and surface tension forces on E. Terminal velocities of ellipsoidal bubbles are accurately predicted by using the proposed correlation and a drag correlation proposed by Rastello et al. (2011). (C) 2015 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2016年03月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 79, 23 - 30, 英語[査読有り]研究論文（学術雑誌）
- 2015年12月, 混相流, 29(4) (4), 343 - 350, 日本語気水分離器における液相分離部が気液分離性能に及ぼす影響[査読有り]研究論文（学術雑誌）
- 2015年11月, ANS Conference 2015 Winter Meeting and Nuclear Technology Expo, 2 pages, 英語Two-Phase Swirling Flow in a Gas-Liquid Separator[査読有り]研究論文（国際会議プロシーディングス）
- Mass transfer from single fully contaminated carbon dioxide bubbles rising through vertical pipes was measured to investigate the effects of surfactants. The diameter ratio of the bubble diameter to the pipe diameter was varied to cover various bubble shapes from ellipsoidal to Taylor bubbles. Triton X-100 and 1-octanol served as surfactants. The mass transfer rates, k(L), of contaminated bubbles in the ellipsoidal regime were smaller than those of clean bubbles due to the interface immobilization. The k(L) of Taylor bubbles contaminated with Triton X-100 and 1-octanol showed varying trends because of the difference in the surfactant distributions at the bubble interfaces which strongly depend on the Hatta number. The Sherwood numbers of contaminated bubbles were well correlated by introducing the bubble diameter, at which the bubble shape transits from ellipsoidal to Taylor bubbles.WILEY-V C H VERLAG GMBH, 2015年11月, CHEMICAL ENGINEERING & TECHNOLOGY, 38(11) (11), 1955 - 1964, 英語[査読有り][招待有り]研究論文（学術雑誌）
- 2015年10月, 7th European-Japanese Two-Phase Flow Group Meeting, 3 pages, 英語Evaluation of Interfacial Shear Stress of Single Drop Using Spatiotemporal Filter Velocimetry[査読有り]研究論文（国際会議プロシーディングス）
- Effects of liquid viscosity on flows inside and outside a bubble diffuser pipe were investigated in this study. Glycerol-Water solutions and air were used for the liquid and gas phases, respectively. A circular pipe having five aeration holes on the topside and larger opening on the bottom side was placed in a rectangular water tank. The pipe diameter was 20 mm and the hole diameter was 5 mm. The liquid viscosity mu(L) ranged from 1 to 100 mPa s and the total gas flow rate Q(IN) ranged from 3 x 10(-4) to 7 x 10(-4) m(3)/s. The flows inside and outside the pipe were observed using a high-speed video camera and the gas flow rates from each aeration hole were measured by capturing bubbles generated from each hole. The liquid entered into the pipe through the aeration hole and the opening, resulting in wavy flow pattern with intermittent slugging. The conclusiohs obtained are as follows: (1) the increases in Q(IN) and mu(L) decrease the frequency of slugging disturbing aeration, so that the gas flow rate from each hole becomes uniform at high Q(IN) and mu(L), (2) the main triggers of slugging are interfacial waves at low mu(L) and interface deformation caused by liquid drops falling from the aeration holes into the pipe at high mu(L), (3) a slugging criterion for inviscid fluids proposed by Mishima and Ishii (1980) is applicable not only to low mu(L) liquids but also to high mu(L) liquids because the wave growth rate is much larger than the viscous damping rate even at high viscosities tested in the present experiments, and (4) the Davidson-Schuler correlation (Davidson and Schuler, 1960) gives reasonable estimations of the bubble diameter not only for low viscosity liquids but also for high viscosity liquids, provided that the influence of slugging is small. (C) 2015 Elsevier Inc. All rights reserved.ELSEVIER SCIENCE INC, 2015年09月, EXPERIMENTAL THERMAL AND FLUID SCIENCE, 66, 197 - 205, 英語[査読有り]研究論文（学術雑誌）
- 2015年09月, 9th International Symposium on Measurement Techniques for Multiphase Flow, 2 pages, 英語Effects of Liquid Height on Gas Holdup in a Rectangular Bubble Column[査読有り]研究論文（国際会議プロシーディングス）
- 2015年07月, ASME-JSME-KSME Joint Fluids Engineering Conference 2015 (AJK2015-FED), AJK2015 - 21048, 2 pages, 英語Motion of Contaminated Drops Falling onto Flat Plate[査読有り]研究論文（国際会議プロシーディングス）
- The purpose of this study is to derive a counter-current flow limitation (CCFL) correlation and evaluate its uncertainty for steam generator (SG) U-tubes in a pressurized water reactor (PWR). Experiments were conducted to evaluate effects of the liquid viscosity on CCFL characteristics using air-40wt% or air-60wt% glycerol water solution and saturated steam-water at atmospheric pressure with vertical pipes simulating the lower part of the SG U-tubes. The steam-water experiments confirmed that CCFL characteristics could be expressed in terms of the Wallis parameters (J(G)* and J(L)*) for the pipe diameters of D = 14, 20, and 27mm. A CCFL correlation was derived using the ratio mu(G)/mu(L) of the viscosities of the gas and liquid phases, mu(G) and mu(L), as a correction term representing effects of fluid properties, where J(G)(*1/2)(mu(G)/mu(L))(-0.07) was expressed by a cubic function of J(L)(*1/2)(mu(G)/mu(L))(0.1). In the correlation, the constant C indicating the value of J(G)(*1/2)(mu(G)/mu(L))(-0.07) at J(L)* = 0 was (1.04 +/- 0.05), and this uncertainty of +/- 0.05 would cover most of the previous experimental data including the ROSA-IV/LSTF data at 1, 3, and 7 MPa.TAYLOR & FRANCIS LTD, 2015年06月, JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY, 52(6) (6), 887 - 896, 英語[査読有り]研究論文（学術雑誌）
- Implicit Temperature-Correction-based Immersed-Boundary Thermal Lattice-Boltzmann Method with Two Relaxation TimesIn the present paper, we verify the effectiveness of the two-relaxation-time (TRT) collision operator in reducing boundary slip of temperature computed by the immersed boundary-thermal lattice Boltzmann method (IB-TLBM). In the linear collision operator of the TRT, we decompose the distribution function into symmetric and antisymmetric components and define the relaxation parameters for each part. The Chapman-Enskog expansion indicates that one relaxation time for the antisymmetric component is related to the thermal conductivity. We derive the theoretical relation between a temperature slip at the boundary and reveal that the relaxation time for the symmetric part controls the temperature at the boundary and boundary slip of temperature computed by the IB-TLBM. We apply the IB-TLBM based on the implicit correction method with two relaxation times for the natural convection in a square enclosure containing a circular cylinder. The streamline, isotherms, and average Nusselt number calculated by the proposed method agree well with those of previous numerical studies.一般社団法人日本機械学会, 2015年05月, 23th International Conference on Nuclear Engineering (ICONE23), 2015(23) (23), 10 pages, ICONE23 - 1039, 英語[査読有り]研究論文（国際会議プロシーディングス）
- Mass transfer from single carbon dioxide bubbles rising through contaminated water in a vertical pipe of 12.5 mm diameter was measured to investigate effects of surfactant. The bubble diameter was widely varied to cover various bubble shapes such as spheroidal, wobbling, cap and Taylor bubbles. The gas and liquid phases were 99.9% purity carbon dioxide and a surfactant solution made of purified water and Triton X-100. Comparison of mass transfer rates between contaminated and clean bubbles made clear that the surfactant decreases the mass transfer rates of small bubbles. The Sherwood numbers of small bubbles in the extreme cases, i.e. zero and the highest surfactant concentrations, are well correlated in terms of the bubble Reynolds number, Schmidt number and the ratio, 2, of the bubble diameter to pipe diameter. The Sherwood numbers at intermediate surfactant concentrations, however, are not well correlated using available correlations. The mass transfer rates of Taylor bubbles also decrease with increasing the surfactant concentration. They however increase with the diameter ratio and approach those of clean Taylor bubbles as 2 increases. The main cause of this tendency was revealed by interface tracking simulations, i.e. the surfactant adsorbs only in the bubble tail region and the nose-to-side region is almost clean at high A. (C) 2014 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2015年04月, INTERNATIONAL JOURNAL OF HEAT AND MASS TRANSFER, 83, 652 - 658, 英語[査読有り]研究論文（学術雑誌）
- The time, t(C), elapsed from bubble contact to the rupture of the liquid film between two bubbles in a quasi two-dimensional column was measured to investigate the effects of the particle diameter on bubble coalescence. The particle diameter ranged from 60 to 150 mu m and the particle volumetric concentration ranged from 0 to 0.50 (50%). The effects of particle diameter on bubble coalescence were evaluated through the experiments and were implemented into a bubble lcoalescence model for a multi-fluid model. Distributions of the gas holdup in a three-dimensional slurry bubble column were also measured using an electrical conductivity probe to obtain experimental data for validation of the coalescence model. The conclusions obtained are as follows: (1) t(C) decreases, in other words, bubble coalescence is enhanced with decreasing the particle diameter, which results in the reduction of the gas holdup of the slurry bubble column, (2) the particle-effect multiplier to t(C) is of great use to take into account the effects of the particle concentration and diameter in multi-fluid simulations of slurry bubble columns, and (3) the interaction between the bubble interfaces and particles in the bubble coalescence process after bubble contact is a local phenomenon and rarely depends on the macroscopic behavior of the interface, which allows us to develop a correlation for particle-induced enhancement of bubble coalescence for three-dimensional bubble columns even with a small experimental setup like a Hele-Shaw cell.SOC CHEMICAL ENG JAPAN, 2015年03月, JOURNAL OF CHEMICAL ENGINEERING OF JAPAN, 48(3) (3), 181 - 189, 英語[査読有り]研究論文（学術雑誌）
- Effects of numerical treatments for the surface tension evaluation on predictions of the motions of droplets ranging from micron to sub-micron meters were investigated. Various combinations of schemes for evaluating the normal to the interface and interface curvature were examined, i.e. the ALE (arbitrary Lagrangian-Eulerian) like scheme and BFA (balanced-force algorithm) for the normal vector and CSF (continuum surface force) and HF (height function) for the interface curvature. The interface motion was predicted using THAINC (tangent of hyperbola with adaptive slope for interface capturing) proposed in our previous study. Numerical errors in pressure and velocity were examined for neutrally buoyant drops of 1 mm in radius to validate the code, which confirmed that the results were similar to those reported in literature: the combination of BFA and HF gave the lowest errors. The droplet size was reduced to 0.1 μm to investigate the accuracy of the schemes for droplet sizes found in industrial coating processes. The static contact angle was then taken into account in the code. The effect of implementation on the errors was examined. The reduction of droplet sizes and implementation of contact angle had no substantial effect on the order of errors. A model for the dynamic contact angle was also implemented and the wetting behaviour of a drop of 1.14 mm in radius was well predicted. Finally a simulation of the wetting behaviour of a sub-micron meter droplet demonstrated that the present code combining BFA, HF and the dynamic contact angle model is accurate in predicting the motion of sub-micron meter droplets.Multi-Science Publishing Co. Ltd, 2015年03月, Journal of Computational Multiphase Flows, 7(1) (1), 15 - 32, 英語[査読有り]研究論文（学術雑誌）
- Sound generated by vapor bubbles condensing in subcooled water was measured to investigate effects of the degree, ΔT, of subcooling. Saturated steam from a boiler was injected into a horizontal water pipe flow through an orifice of 1 mm in diameter. The duct diameter was 16 mm and the Reynolds number of water was about 5 × 104. The subcooling was varied from 10 to 40 K. The sound level increased with increasing ΔT for ΔT ≤ 20 K, whereas it was independent of ΔT at higher ΔT. Video images of condensing bubbles revealed that the interfacial velocity of a condensing bubble had the same trend as that of the sound level, and therefore, the interfacial velocity in condensation is the key in the sound level. The rapid condensation at high ΔT induced emission of microbubbles and secondary pressure pulse generation. The latter also contributes to the increase in the sound level.Begell House Inc., 2015年, Multiphase Science and Technology, 27(1) (1), 67 - 75, 英語[査読有り]研究論文（学術雑誌）
- Multi-fluid simulations of heterogeneous bubbly flows in an air-water bubble column were carried out to verify the speculation that a heterogeneous bubbly flow is predictable without turbulence models such as k-epsilon and LES models, provided that the velocity fluctuation caused by large-scale vortical flow structures prevails over the bubble-induced and shear-induced turbulences. Experiments on the heterogeneous air-water bubbly flows in a rectangular bubble column were also carried out to obtain experimental data of the mean velocity, fluctuation velocity and void fraction for validation of the numerical method. A small LDV probe developed in our previous study was utilized to measure the liquid velocity in the column at high spatial and high temporal resolutions. The distribution of void fraction was measured using an electrical conductivity probe. The conclusions obtained under the present experimental conditions are as follows: (1) the small LDV probe is of great use in measuring the distributions of mean and fluctuation velocities of the liquid phase at high spatial and high temporal resolutions, (2) the velocity fluctuation in the heterogeneous regime in the bubble column is mainly due to large-scale vortical structures, and (3) the multi-fluid model can give good predictions of a heterogeneous bubbly flow without using the turbulence models, provided that large-scale vortical structures in the flow prevail over the bubble-induced and shear-induced turbulences. (C) 2014 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2014年12月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 67, 111 - 121, 英語[査読有り]研究論文（学術雑誌）
- 2014年11月, ANS Conference 2014 Winter Meeting and Nuclear Technology Expo, CD-ROM, 2 pages, 英語CCFL of Steam-Condensate Two-Phase Flow in a U-Tube[査読有り]研究論文（国際会議プロシーディングス）
- 2014年09月, Proc. 2nd International Symposium on Multiscale Multiphase Process Engineering, pp. 94 - 99, 英語Effects of Surfactant on Mass Transfer from Single Carbon Dioxide Bubbles in Vertical Pipes[査読有り]研究論文（国際会議プロシーディングス）
- 2014年09月, Proc. 2nd International Symposium on Multiscale Multiphase Process Engineering, pp. 230 - 235, 英語Effects of Height-To-Diameter Ratio on Gas Holdup in a Slurry Bubble Column[査読有り]研究論文（国際会議プロシーディングス）
- 2014年09月, Proc. 2nd International Symposium on Multiscale Multiphase Process Engineering, pp. 106 - 115, 英語An Interface Capturing Scheme for Free-Surface Flows in a Flow Channel Consisting of Solid Obstacles[査読有り]研究論文（国際会議プロシーディングス）
- 2014年08月, Proc. 15th International Heat Transfer Conference, IHTC-15, IHTC15 - 08960, 11 pages, 英語Bubble Growth in Microgravity under the Action of Electric Forces: Experiments and Numerical Simulation[査読有り]研究論文（国際会議プロシーディングス）
- 2014年07月, International Conference on Numerical Methods for Multiphase Flow 2014, ICNMMF2014, 2 pages, 英語Numerical Simulations of Clean and Contaminated Bubbles in Linear Shear Flows[査読有り]研究論文（国際会議プロシーディングス）
- 2014年07月, 17th International Symposium on Applications of Laser Techniques to Fluid Mechamics (LISBON 2014), 英語Evaluation of Interfacial Shear Stress of Single Drop using Spatiotemporal Filter Velocimetry[査読有り]研究論文（国際会議プロシーディングス）
- Effects of numerical treatment of viscous and surface tension forces on predicted motion of an interface are investigated. The viscous force and tangential surface tension force (Marangoni force) are treated in two different ways: one is a smeared-out interface method and the other is a ghost fluid method. In the smeared-out interface method, the arithmetic mean and harmonic mean are tested for evaluation of the viscous stress. Linear shear flows, single oscillating drops and surface tension waves are simulated using these methods. These benchmark tests show that the ghost fluid method gives the most accurate evaluations of the viscous and the Marangoni forces. The harmonic mean can give good evaluations of the tangential viscous stress at clean interface, whereas large errors in the viscous stress are caused for contaminated interface. Although the arithmetic mean causes some errors for clean and contaminated interfaces compared to the ghost fluid method, the error can be reduced by increasing spatial resolution and this method can give more reasonable predictions for contaminated interface than the harmonic mean. The arithmetic mean is therefore recommended rather than the harmonic mean when using the smeared-out interface method.Multi-Science Publishing Co. Ltd, 2014年06月, Journal of Computational Multiphase Flows, 6(2) (2), 111 - 126, 英語[査読有り]研究論文（学術雑誌）
- 2014年02月, Journal of Nuclear Science and Technology, 51(5) (5), 580 - 589, 英語Void distribution and bubble motion in bubbly flows in a 4×4 rod bundle. Part II: numerical simulation[査読有り]研究論文（学術雑誌）
- Lack of local void fraction data in a rod bundle makes it difficult to validate a numerical method for predicting gas-liquid two-phase flow in the bundle. Distributions of local void fraction and bubble velocity in each subchannel in a 4x4 rod bundle were, therefore, measured using a double-sensor conductivity probe. Liquid velocity in the subchannel was also measured using laser Doppler velocimetry (LDV) to obtain relative velocity between bubbles and the liquid phase. The size and pitch of rods were 10 and 12.5mm, respectively. Air and water at atmospheric pressure and room temperature were used for the gas and liquid phases, respectively. The volume fluxes of gas and liquid phases ranged from 0.06 to 0.15m/s and from 0.9 to 1.5m/s, respectively. Experimental results showed that the distributions of void fraction in inner and side subchannels depend not only on lift force acting on bubbles but also on geometrical constraints on bubble dynamics, i.e. the effects of rod walls on bubble shape and rise velocity. The relative velocity between bubbles and the liquid phase in the subchannel forms a non-uniform distribution over the cross-section, and the relative velocity becomes smaller as bubbles approach the wall due to the wall effects.TAYLOR & FRANCIS LTD, 2014年02月, JOURNAL OF NUCLEAR SCIENCE AND TECHNOLOGY, 51(2) (2), 220 - 230, 英語[査読有り]研究論文（学術雑誌）
- In the present paper, we verify the effectiveness of the two-relaxation-time (TRT) collision operator in reducing boundary slip computed by the immersed boundary-lattice Boltzmann method (IB-LBM). In the linear collision operator of the TRT, we decompose the distribution function into symmetric and antisymmetric components and define the relaxation parameters for each part. The Chapman-Enskog expansion indicates that one relaxation time for the symmetric component is related to the kinematic viscosity. Rigorous analysis of the symmetric shear flows reveals that the relaxation time for the antisymmetric part controls the velocity gradient, the boundary velocity, and the boundary slip velocity computed by the IB-LBM. Simulation of the symmetric shear flows, the symmetric Poiseuille flows, and the cylindrical Couette flows indicates that the profiles of the numerical velocity calculated by the TRT collision operator under the IB-LBM framework exactly agree with those of the multirelaxation time (MRT). The TRT is as effective in removing the boundary slip as the MRT. We demonstrate analytically and numerically that the error of the boundary velocity is caused by the smoothing technique using the delta function used in the interpolation method. In the simulation of the flow past a circular cylinder, the IB-LBM based on the implicit correction method with the TRT succeeds in preventing the flow penetration through the solid surface as well as unphysical velocity distortion. The drag coefficient, the wake length, and the separation points calculated by the present IB-LBM agree well with previous studies at Re = 10, 20, and 40.AMER PHYSICAL SOC, 2014年02月, Physical Review E, 89(2) (2), 023307, 英語[査読有り]研究論文（学術雑誌）
- Interface tracking simulations of mass transfer from Taylor bubbles were carried out to investigate effects of shape oscillation on the mass transfer. Mass transfer from carbon dioxide Taylor bubbles in a glycerol-water solution was also measured to obtain experimental data for the validation of the numerical method. A high spatial resolution was used to resolve thin concentration boundary layers on the bubble interface, which enabled us to capture the agitation of boundary layer due to interfacial waves. The predicted mass transfer coefficients were in good agreements with the experimental data, provided that the spatial resolution was high enough to capture thin concentration boundary layers of high Schmidt number Taylor bubbles. The simulations have made it clear that the effect of agitation of concentration boundary layer due to interfacial waves on the total mass transfer rate is small and a fluctuation of the Sherwood number is caused by a fluctuation of the bubble surface area. This result implies that interface tracking methods using boundary layer approximations can give accurate predictions for mass transfer from bubbles even when the flow field is not fully resolved when the wave effect on the net mass transfer is negligible. (C) 2013 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2014年01月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 58, 236 - 245, 英語[査読有り]研究論文（学術雑誌）
- To investigate the effects of hydrophilic particles on slurry bubble flows in a bubble column, distributions of the local gas holdup and the bubble frequency are measured using an electric conductivity probe. Particles are made of silica and their diameter is 100 mu m. The particle volumetric concentration C-S is varied from 0 to 0.40. The measured data imply that the presence of particles promotes bubble coalescence. The film drainage time for two coalescing bubbles in a quasi two-dimensional bubble flow in a small vessel is also measured to quantitatively evaluate the particle effect on coalescence. A particle-effect multiplier is introduced into a coalescence efficiency model by taking into account the data of film drainage time and is implemented into a multi-fluid model. The main conclusions obtained are as follows: (1) the local gas holdup and bubble frequency in slurry bubble flows decrease with increasing the particle concentration, (2) the hydrophilic particles enhance bubble coalescence and the enhancement saturates at C-S similar or equal to 0.45, (3) the particle effect on coalescence is well accounted for by introducing the particle-effect multiplier to the film drainage time, and (4) the multi-fluid model can give good predictions for the distribution of the local gas holdup in the slurry bubble column. (C) 2013 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2014年01月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 58, 154 - 167, 英語[査読有り]研究論文（学術雑誌）
- An interface-capturing method for simulating two-phase flows in complex geometries is proposed. A structured orthogonal grid is used, and the presence of stationary solid boundaries is taken into account by using the volume fraction of the solid phase in a computational cell. The transport equation of the fluid volume fraction in the presence of solid boundaries is derived. The fluid phases are transported by using the THAINC (tangent of hyperbola with adaptive slope for interface capturing) method. An immersed boundary method is implemented to accurately calculate the volume flux of each fluid phase through a computational cell face, a part of which is blocked by the solid phase. Transportation of a fluid square in the presence of solid boundaries is carried out, to con firm that the errors in shape and volume conservation are low. The applicability of the proposed method to two-phase flows in complex geometries is examined through simulations of a flow about cylindrical tubes in a staggered arrangement and a dam break problem with an obstacle. In the former the pressure drop in the single-phase condition agrees well with the available correlation and the gas-liquid behavior in the tube geometry is qualitatively well predicted. The predicted liquid flow in the dam break problem also agrees well with available experimental data.SOC CHEMICAL ENG JAPAN, 2014年, JOURNAL OF CHEMICAL ENGINEERING OF JAPAN, 47(3) (3), 230 - 240, 英語[査読有り]研究論文（学術雑誌）
- MASS TRANSFER FROM SINGLE CARBON DIOXIDE BUBBLES IN CONTAMINATED WATERMass transfer from single carbon dioxide bubbles rising through contaminated water in a vertical pipe of 12.5 mm diameter was measured to investigate effects of surfactant. The bubble diameter was widely varied to cover various bubble shapes such as spheroidal, wobbling, cap and Taylor bubbles. The gas and liquid phases were 99.9 % purity carbon dioxide and a surfactant solution made of purified water and Triton X-100. Comparison of mass transfer rates between contaminated and clean bubbles made clear that the surfactant decreases the mass transfer rates of small bubbles. The Sherwood number of small bubbles in the extreme cases, i.e. zero and the highest surfactant concentrations, is well correlated in terms of the bubble Reynolds number, Schmidt number and the ratio, lambda, of the bubble diameter to pipe diameter. The Sherwood numbers at intermediate surfactant concentration, however, are not well correlated using available correlations. The mass transfer rates of Taylor bubbles also decrease with increasing the surfactant concentration. They however increase with the diameter ratio and approaches that of clean Taylor bubbles as 2, increases. The main cause of this tendency was revealed by interface tracking simulations, i.e. the surfactant adsorbs only in the bubble tail region and the nose-to-side region is almost clean at high lambda.AMER SOC MECHANICAL ENGINEERS, 2014年, ASME FLUIDS ENGINEERING DIVISION SUMMER MEETING - 2014, VOL 1D: SYMPOSIA, FEDSM2014 - 21103, 8 pages, 英語[査読有り]研究論文（国際会議プロシーディングス）
- In this paper, we investigated the demulsification behavior of oil-in-water (O/W) emulsions during membrane permeation in the oil-water separation process using a numerical simulation approach. To accurately deal with the large deformation of the oil-water interface by coalescence and wetting, and to estimate the volume of the coalesced oil droplet, the coupled level set and volume-of-fluid method was used as the interface capturing method. We applied the simulation model to the permeation of O/W emulsions through a membrane pore, and then investigated the effects of the wettability of the membrane surface, filtration flux, and pore size on the demulsification efficiency. The results showed that oil droplets were likely to coalesce on the outlet membrane surface. High wettability on the membrane surface and low fluid velocity inside the pore increased the demulsification efficiency. This is the first work to numerically simulate the demulsification behavior of emulsions through membranes in the oil-water separation process.ROYAL SOC CHEMISTRY, 2014年, SOFT MATTER, 10(40) (40), 7985 - 7992, 英語[査読有り]研究論文（学術雑誌）
- Experiments on bubbly flows in a rectangular column are carried out to investigate effects of inlet condition on How structure. An air diffuser having 35 nozzles is used. The inlet gas How rates from 34 nozzles are uniform and the lift coefficients of bubbles are small at this flow rate, whereas the gas inflow from the remaining one nozzle is varied to change the sign and magnitude of lift coefficients. The main conclusions obtained are as follows: (1) bubbly flow in a bubble column is apt to be heterogeneous even with small non-uniformity in the distribution of gas inflow; (2) when the nozzle causing non-uniform gas inflow is located in the center of the diffuser plate and the inflow distribution is symmetric, heterogeneous structure is formed mainly due to the lift-induced flow instability pointed out by Lucas et al. (2005); (3) when the nozzle causing non-uniform gas inflow is located near the side wall of the column and the inflow distribution is asymmetric, heterogeneous structure is induced by the Coanda effect; and (4) a multi-fluid model has a potential of predicting the effects of inlet condition on flow structure in a bubble column. (C) 2013 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2013年12月, CHEMICAL ENGINEERING SCIENCE, 104, 166 - 176, 英語[査読有り]研究論文（学術雑誌）
- Immersed boundary-lattice Boltzmann methods (IB-LBM) with the single-relaxation time (SRT) cause non-physical distortion in fluid velocity when the Reynolds number is low, i.e. the relaxation time τ is high, and IB-LBM requires high spatial resolution to stably simulate high Reynolds number flows. An immersed boundary-finite difference lattice Boltzmann method (IB-FDLBM) using two-relaxation times (TRT) is therefore proposed in this study to simulate low and high Reynolds number flows stably and accurately. Benchmark problems such as circular Couette flows, flows past a circular cylinder and a sphere at various Reynolds numbers are carried out for validation. The main conclusions obtained are as follows: (1) TRT reduces numerical errors causing non-physical distortion in the fluid velocity at low Reynolds numbers, and accurate predictions are obtained when the parameter Λ, which is a function of the two relaxation times, is low, (2) for stable simulation the parameter Λ should be decreased as the Reynolds number increases, (3) implementation of TRT and the implicit direct forcing method into IB-FDLBM can solve two problems in simulation of low Reynolds number flows, i.e. non-physical velocity distortion and non-physical penetration of flow into the solid body, and (4) IB-FDLBM with TRT gives good predictions of the drag coefficients of a circular cylinder and a sphere in uniform flows for a wide range of the Reynolds number, Re, i.e., 0.1 ≤ Re ≤ 1x104.JSME, 2013年09月, Journal of Fluid Science and Technology, 8(3) (3), 262 - 276, 英語[査読有り]研究論文（学術雑誌）
- 2013年09月, 9th UK - Japan Seminar on Multi-Phase Flow, 日本語Gas-Liquid Two-Phase Flows Inside and Outside an Air Diffuser Pipe[査読有り]研究論文（国際会議プロシーディングス）
- 2013年09月, 9th UK-Japan Seminar on Multi-Phase Flow, 英語Distribution of Void Fraction and Liquid Velocity in Air-Water Bubble Column[査読有り]研究論文（国際会議プロシーディングス）
- In the present paper, a study of bubble dynamics in adiabatic conditions (i.e. with no heat and mass transfer) is carried out in order to understand the fundamental physics ruling the evolution of their interface, and open the way to the comprehension of more complex heat and mass transfer issues involving them. An external electric field is applied to the system to investigate the potentiality of this technique to enhance phase separation and heat transfer. A new numerical method for EHD action on a growing bubble has been developed combining VOF and level-set methods. To experimentally validate it, well-defined and simple boundary conditions have been selected: all the forms of heat and mass transfer have been excluded from the system, and an axisymmetric electric field configuration has been chosen in order to adopt a 2-D simulation. In the experimental apparatus a steady bubble stemming from an orifice in a still fluid (FC-72) has been produced, and a dc potential up to 20. kV was applied to a ring-shaped electrode coaxial with the orifice. In this way, in a first instance, all the dynamical effects have been eliminated, and a clear location of the three phase contact line, which is pinned to the orifice, has been defined. The bubble profile obtained from the experiment has been digitized and compared successfully with the outcomes of the numerical method. Furthermore, the study of the local curvature of the interface allowed the development of a new theory of the local electrical stress to be included in the capillary equation, which compares well with experimental data. Finally, the validated numerical method has been used to predict the bubble shape in the absence of gravity, in preparation of microgravity experiments. © 2013 Elsevier Inc.2013年09月, Experimental Thermal and Fluid Science, 49, 160 - 168, 英語[査読有り]研究論文（学術雑誌）
- 2013年06月, 8th World Conference on Experimental Heat Transfer, Fluid Mechanics, and Thermodynamics (ExHFT8), USB memory, 8 pages, 英語Effects of Inlet Condition on Flow Structure of Bubbly Flow in Rectangular Column[査読有り]研究論文（国際会議プロシーディングス）
- 2013年05月, 15th International Topical Meeting on Nuclear Reactor Thermal - Hydraulics, NURETH-15, NURETH15 - 244, USB memory, 12 pages, 英語Void Distribution and Bubble Motion in Bubbly Flows in a 4 X 4 Rod Bundle[査読有り]研究論文（国際会議プロシーディングス）
- 2013年05月, 15th International Topical Meeting on Nuclear Reactor Thermal - Hydraulics, NURETH-15, NURETH15 - 245, USB memory, 12 pages, 英語NUMERICAL SIMULATION OF BUBBLY FLOWS IN A 4 X 4 ROD BUNDLE USING A MULTI-FLUID MODEL[査読有り]研究論文（国際会議プロシーディングス）
- 2013年05月, 8th International Conference on Multiphase Flow, ICMF2013, ICMF2013 - 229, USB memory, 9 pages, 英語Mass Transfer from Single Carbon Dioxide Bubbles in Glycerol-Water Solution[査読有り]研究論文（国際会議プロシーディングス）
- 2013年05月, 8th International Conference on Multiphase Flow, ICMF2013, ICMF2013 - 286, USB memory, 8 pages, 英語Effects of Shape Oscillation on Mass Transfer from a Taylor Bubble[査読有り]研究論文（国際会議プロシーディングス）
- 2013年05月, 8th International Conference on Multiphase Flow, ICMF2013, ICMF2013 - 157, USB memory, 13 pages, 英語Effects of Hydrophilic Particles on Bubble Coalescence in Slurry Bubble Column[査読有り]研究論文（国際会議プロシーディングス）
- Terminal velocities and shapes of drops rising through vertical pipes in clean and fully-contaminated systems are measured by using a high-speed video camera and an image processing method. Silicon oils and glycerol water solutions are used for the dispersed and continuous phases, respectively. Triton X-100 is used for surfactant. Clean and contaminated drops take either spherical, spheroidal or deformed spheroidal shapes when the diameter ratio λ is less than a critical value, λC, whereas they take bullet shapes for λ> λC (Taylor drops). The applicability of available drag and Froude number correlations is examined through comparisons with the measured data. Effects of surfactant on the shape and terminal velocity of a Taylor drop are also discussed based on the experimental data and interface tracking simulations. The conclusions obtained are as follows: (1) drag and Froude number correlations proposed so far give reasonable estimations of the terminal velocities of clean drops at any λ, (2) the terminal velocities of contaminated drops are well evaluated by making the viscosity ratio μ* infinity in the drag correlation for clean drops in the viscous force dominant regime, (3) the effects of surfactant on the shape and terminal velocity of a Taylor drop become significant as the Eötvös number, EoD, decreases and μ* increases, and (4) the reduction in surface tension due to the addition of surfactant would be the cause of the increase in the terminal velocity and elongation of a contaminated Taylor drop. © 2012 Elsevier Ltd.2013年03月, International Journal of Multiphase Flow, 49, 8 - 23, 英語[査読有り]研究論文（学術雑誌）
- The applicability of the immersed boundary-finite difference lattice Boltzmann method (IB-FDLBM) to high Reynolds number flows about a circular cylinder is examined. Two-dimensional simulations of flows past a stationary circular cylinder are carried out for a wide range of the Reynolds number, Re, i.e., 1 Re 1 105. An immersed boundary-lattice Boltzmann method (IB-LBM) is also used for comparison. Then free-falling circular cylinders are simulated to demonstrate the feasibility of predicting moving particles at high Reynolds numbers. The main conclusions obtained are as follows: (1) steady and unsteady flows about a stationary cylinder are well predicted with IB-LBM and IB-FDLBM, provided that the spatial resolution is high enough to satisfy the conditions of numerical stability, (2) high spatial resolution is required for stable IB-LBM simulation of high Reynolds number flows, (3) IB-FDLBM can stably simulate flows at very high Reynolds numbers without increasing the spatial resolution, (4) IB-FDLBM gives reasonable predictions of the drag coefficient for 1 Re1105, and (5) IB-FDLBM gives accurate predictions for the motion of free-falling cylinders at intermediate Reynolds numbers.Multi-Science Publishing, 2013年03月, Journal of Computational Multiphase Flows, 5(1) (1), 27 - 41, 英語[査読有り]研究論文（学術雑誌）
- 2012年12月, 計算数理工学論文集, Vol. 12, pp. 37-42, 日本語Two-Relaxation-Time を用いたIB-LBMによる熱流動解析[査読有り]研究論文（学術雑誌）
- 日本計算数理工学会, 2012年12月, ICOME2012/JASCOME2012, 12, 37 - 42, 日本語Numerical analysis of thermal-hydraulics by immersed boundary-lattice Boltzmann method using two-relaxation-time研究論文（国際会議プロシーディングス）
- Immersed boundary-lattice Boltzmann methods (IB-LBM) with the single-relaxation time (SRT) cause non-physical distortion in fluid velocity when the Reynolds number is low, i.e. the relaxation time τ is high, and IB-LBM requires high spatial resolution to stably simulate high Reynolds number flows. An immersed boundary-finite difference lattice Boltzmann method (IB-FDLBM) using two-relaxation times (TRT) is therefore proposed in this study to simulate low and high Reynolds number flows stably and accurately. Benchmark problems such as circular Couette flows, flows past a circular cylinder and a sphere at various Reynolds numbers are carried out for validation. The main conclusions obtained are as follows: (1) TRT reduces numerical errors causing non-physical distortion in the fluid velocity at low Reynolds numbers, and accurate predictions are obtained when the parameter Λ, which is a function of the two relaxation times, is low, (2) for stable simulation the parameter Λ should be decreased as the Reynolds number increases, (3) implementation of TRT and the implicit direct forcing method into IB-FDLBM can solve two problems in simulation of low Reynolds number flows, i.e. non-physical velocity distortion and non-physical penetration of flow into the solid body, and (4) IB-FDLBM with TRT gives good predictions of the drag coefficients of a circular cylinder and a sphere in uniform flows for a wide range of the Reynolds number, Re, i.e., 0.1 ≤ Re ≤ 1x104.The Japan Society of Mechanical Engineers, 2012年12月, 9th International Conference on CFD in the Minerals and Process Industries,, CD-ROM, 6 pages(3) (3), 262 - 276, 英語[査読有り]研究論文（国際会議プロシーディングス）
- 2012年12月, 9th International Conference on CFD in the Minerals and Process Industries,, CD-ROM, 6 pagesEFFECTS OF NUMERICAL TREATMENT OF VISCOUS AND SURFACE TENSION FORCES ON PREDICTED INTERFACE MOTION[査読有り]研究論文（国際会議プロシーディングス）
- 2012年11月, Korea-Japan CFD Workshop 2012,, 8 pagesInterface Tracking Simulation of Bubbles in Clean and Contaminated Systems研究論文（国際会議プロシーディングス）
- 2012年09月, 6th Japanese-European Two-Phase Flow Group Meeting, USB memory, 10 pagesBUBBLE TRACKING SIMULATION OF BUBBLE-INDUCED PSEUDO TURBULENCE[査読有り]研究論文（国際会議プロシーディングス）
- A downscaled model of a steam separator is used to experimentally examine the effects of swirler shape on swirling annular flow. Several swirlers having different vane angles, different hub diameters, and different number of vanes are installed in the separator. Flow patterns, the ratio of the separated liquid flow rate to the total liquid flow rate, liquid film thickness, and pressure drop are measured. An improved swirler, which has a small hub diameter, the reduced number of vanes and modified vane shape, is also designed and its performance is examined. As a result, the following conclusions are obtained. (1)The reduction of the hub diameter and the number of vanes do not lower the separator performance. (2) The decrease in vane angle deteriorates the separation performance. This is due to the decrease in the centrifugal force, i.e., the decrease in the amount of droplet deposition. (3) The improved swirler effectively reduces the pressure drop while keeping high separator performance. (C) 2011 Elsevier B.V. All rights reserved.ELSEVIER SCIENCE SA, 2012年08月, Nuclear Engineering and Design, 249, 63 - 70, 英語[査読有り]研究論文（学術雑誌）
- 2012年06月, International Conference on Numerical Methods in Multiphase Flows, No. 7-1Numerical Simulation of Mass Transfer From a Taylor Bubble with Shape Oscillation研究論文（国際会議プロシーディングス）
- An immersed boundary-lattice Boltzmann method (IB-LBM) using a two-relaxation time model (TRT) is proposed. The collision operator in the lattice Boltzmann equation is modeled using two relaxation times. One of them is used to set the fluid viscosity and the other is for numerical stability and accuracy. Adirect-forcing method is utilized for treatment of immersed boundary. A multi-direct forcing method is also implemented to precisely satisfy the boundary conditions at the immersed boundary. Circular Couette flows between a stationary cylinder and a rotating cylinder are simulated for validation of the proposed method. The method is also validated through simulations of circular and spherical falling particles. Effects of the functional forms of the direct-forcing term and the smoothed-delta function, which interpolates the fluid velocity to the immersed boundary and distributes the forcing term to fixed Eulerian grid points, are also examined. As a result, the following conclusions are obtained: (1) the proposed method does not cause non-physical velocity distribution in circular Couette flows even at high relaxation times, whereas the single-relaxation time (SRT) model causes a large non-physical velocity distortion at a high relaxation time, (2) the multi-direct forcing reduces the errors in the velocity profile of a circular Couette flow at a high relaxation time, (3) the two-point delta function is better than the four-point delta function at low relaxation times, but worse at high relaxation times, (4) the functional form of the direct-forcing term does not affect predictions, and (5) circular and spherical particles falling in liquids are well predicted by using the proposed method both for two-dimensional and three-dimensional cases.2012年06月, Journal of Computational Multiphase Flows, 4(2) (2), 193 - 209, 英語[査読有り]研究論文（学術雑誌）
- 2012年06月, 2012 Japan-U.S. seminar on Two-Phase flow Dynamics, USB memory, 10 pagesEFFECTS OF SURFACTANT ON THE MOTION OF BUBBLES IN LINEAR SHEAR FLOWS研究論文（国際会議プロシーディングス）
- Effects of the total gas flow rate on the water level in a diffuser pipe for a membrane bioreactor, the gas flow rate from each aeration hole and the bubble diameter are investigated. The diffuser has evenly positioned five aeration holes on the top and a larger hole on the bottom for introducing the liquid into the pipe. The gas flow rate from each aeration hole is measured by capturing generated bubbles. The water level and gas velocity inside the diffuser are computed by processing video images. The bubble diameter is calculated using the gas flow rate and the bubble generation frequency measured from the video images. The conclusions obtained are as follows: (1) the gas flow rate from the aeration hole depends on the water level inside the diffuser and becomes constant for all the holes as the total gas flow rate increases since the high total gas flow rate make the water level uniform, which results in a constant gas pressure in the diffuser, (2) the onset of slugging in the diffuser is well correlated in terms of the local gas velocity and the Mishima-Ishii's slugging model, (3) the increase in the total gas flow rate decreases the water level, causing suppression of the onset of slugging, (4) the diameter of aeration hole strongly affects the gas flow rate from each aeration hole and water level, and (5) the Davidson-Schuler correlation gives reasonable estimations of the bubble diameter, provided that the influence of slugging is not significant.The Japan Society of Mechanical Engineers, 2012年04月, Journal of Fluid Science and Technology, Vol. 7, No. 1, pp. 78-88(1) (1), 78 - 88, 英語[査読有り]研究論文（学術雑誌）
- The purpose of this research was to provide further insights on bubble-induced agitation of heated bulk liquid. Fundamental studies on the bubble disturbance of a stratified thermal layer were carried out for a 6 mm sphere-volume equivalent diameter air bubble suspended in water (Eo=1.2). A video digital image system and thermochromic liquid crystals were used to reproduce the bubble movement as well as the wake drift of the liquid. A three-dimensional interface tracking simulation was used as a numerical tool. The results have revealed a long open wake region that is formed along the fluctuating bubble path. The amplitude of the bubble wake is much larger than that of the bubble path. In addition to longitudinal mixing, strong lateral mixing is also caused by the movement of vortices in the transverse direction. Details of numerical simulations have revealed a wake that tends to form a chain of Omega shaped vortex rings. These "rings" are connected to each other in the near wake region. The coherent effects further downstream lead to more complex vortex patterns in the far wake. The disintegration of the coherent chain of vortices due to bubble surface deformability is highly complex and not clear at this stage. A comparison with liquid crystal temperature response has revealed that the time scale of the mixing is much longer than the ascending bubble residence, approximately 8 s as obtained by several experimental runs. A transverse propagation of entrained cold water has given an estimate of the bubble-induced diffusion to be about 170 times larger than the diffusivity of momentum. The bulk liquid agitation cannot be correctly modeled without taking into account the bubble size dependent wake structure. (C) 2011 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2012年04月, CHEMICAL ENGINEERING SCIENCE, 72, 155 - 171, 英語[査読有り]研究論文（学術雑誌）
- 2012年03月, Journal of Computational Multiphase Flows, 3(4), 247-261Interface Tracking Simulation of Mass Transfer from a Dissolving Bubble[査読有り]研究論文（学術雑誌）
- Effects of soluble surfactant on the terminal velocity of a Taylor bubble rising through a vertical pipe are investigated using an interface tracking method. A level set method is utilized to track the interface. Transport of surfactant in the bulk liquid and at the interface is taken into account. The amount of adsorption and desorption is evaluated using the Frumkin and Levich model. The normal component of surface tension force is computed using a ghost fluid method, whereas the tangential component, i.e., the Marangoni force, is evaluated by making use of the continuum surface force model. Simulations of small air bubbles contaminated with soluble surfactant are carried out for validation. The Marangoni effects on the bubbles, i.e., the surface immobilization and the increase in drag coefficient, are well predicted. Then Taylor bubbles rising through vertical pipes filled with contaminated water at a low Morton number are simulated for various Eotvos numbers, various bulk surfactant concentrations and two different surfactants, i.e., 1-pentanol and Triton X-100. As a result, the following conclusions are obtained: (1) the reduction of surface tension near the bubble nose is the cause of the increase in terminal velocity, (2) the surfactant does not affect the terminal velocities of high Eotvos number bubbles since the bubbles at high Eotvos numbers are independent of surface tension, (3) the terminal velocity of a low Morton number Taylor bubble can be evaluated by making use of available correlations for clean Taylor bubbles, provided that the degree of contamination near the bubble nose is known and the Marangoni effect in the nose region is negligible, and (4) the Hatta number, which is the ratio of the adsorption velocity to the bubble velocity, is a primary factor governing the degree of contamination in the bubble nose region. (C) 2011 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2012年03月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 39, 78 - 87, 英語[査読有り]研究論文（学術雑誌）
- An experimental study on countercurrent flow limitation (CCFL) in vertical pipes is carried out. Effects of upper tank geometry and water levels in the upper and lower tanks on CCFL characteristics are investigated for air-water two-phase flows at room temperature and atmospheric pressure. The following conclusions are obtained: (1) CCFL characteristics for different pipe diameters are well correlated using the Kutateladze number if the tank geometry and the water levels are the same; (2) CCFL occurs at the junction between the pipe and the upper tank both for the rectangular and cylindrical tanks, and CCFL with the cylindrical tank occurs not only at the junction but also inside the pipe at high gas flow rates and small pipe diameters; (3) the flow rate of water entering into the vertical pipe at the junction to the rectangular upper tank is lower than that to the cylindrical tank because of the presence of low frequency first-mode sloshing in the rectangular tank; (4) increases in the water level in the upper tank and in the air volume in the lower tank increase water penetration into the pipe, and therefore, they mitigate the flow limitation.HINDAWI PUBLISHING CORPORATION, 2012年, SCIENCE AND TECHNOLOGY OF NUCLEAR INSTALLATIONS, No. NURETH14-057, CD-ROM (12 p, 英語[査読有り]研究論文（学術雑誌）
- An experimental study on countercurrent flow limitation (CCFL) in vertical pipes is carried out. Effects of upper tank geometry and water levels in the upper and lower tanks on CCFL characteristics are investigated for air-water two-phase flows at room temperature and atmospheric pressure. The following conclusions are obtained: (1) CCFL characteristics for different pipe diameters are well correlated using the Kutateladze number if the tank geometry and the water levels are the same; (2) CCFL occurs at the junction between the pipe and the upper tank both for the rectangular and cylindrical tanks, and CCFL with the cylindrical tank occurs not only at the junction but also inside the pipe at high gas flow rates and small pipe diameters; (3) the flow rate of water entering into the vertical pipe at the junction to the rectangular upper tank is lower than that to the cylindrical tank because of the presence of low frequency first-mode sloshing in the rectangular tank; (4) increases in the water level in the upper tank and in the air volume in the lower tank increase water penetration into the pipe, and therefore, they mitigate the flow limitation.HINDAWI PUBLISHING CORPORATION, 2012年, SCIENCE AND TECHNOLOGY OF NUCLEAR INSTALLATIONS, Vol. 2012, 754724, 10 pages, 英語[査読有り]研究論文（学術雑誌）
- Bubble tracking simulations of bubbly flows in a rectangular column are carried out to examine the applicability of a bubble tracking method to flows involving the bubble-induced pseudoturbulence. Experimental data of void fraction, mean liquid and bubble velocities, and liquid fluctuation velocity are obtained for validation. As a result, the following conclusions are obtained: (1) The order of liquid fluctuation velocity is similar to that of the bubble relative velocity, and the fluctuation is well scaled in terms of the bubble terminal velocity and the void fraction, as Risso & Ellingsen suggested [J. Fluid Mech., vol. 440, pp. 235-268, 2001]. (2) The Reynolds shear stress is not produced by the bubble motion when both void fraction and liquid velocity are uniform. (3) The Reynolds stress model proposed by Lopez de Bertodano et al. [Int. J. Multiphase Flow, vol. 20, pp. 805-818, 1994] underestimates the normal stress, whereas the bubble tracking simulation with a spatial resolution comparable to the bubble size gives better predictions than the model. (4) The fluctuation in liquid velocity induced by bubbles is partly resolved with a spatial resolution comparable to the bubble size, and subgrid turbulence models do not have much influence on predictions. (5) Eddy viscosity models would be a reasonable choice to capture the shear-induced turbulence in the near-wall region in bubbletracking simulations with an inevitably insufficient spatial resolution for large eddy simulation. © 2012 by Begell House, Inc.Begell house, 2012年, Multiphase Science and Technology, 24(3) (3), 197 - 222, 英語[査読有り]研究論文（学術雑誌）
- An immersed boundary method based on a finite difference lattice Boltzmann method (IB-FDLBM) is presented. The FDLBM solves the discrete Boltzmann equation including an additional collision term by using finite difference schemes. The additional term works as a negative viscosity in the macroscopic level and allows us to alter the fluid viscosity while keeping the other relevant parameters of the simulation fixed. The immersed boundary method employs a direct-forcing method, which utilizes external forces at Lagrangian points embedded on immersed boundaries to impose the no-slip boundary condition. Several benchmark simulations are carried out to validate the developed method, i.e., flows past a circular cylinder, a falling particle, and interaction between two falling particles. Couette flows between a stationary and a rotating cylinder are also simulated at various values of the relaxation time for collision. The main conclusions obtained are as follows: (1) steady flows past a stationary circular cylinder are well predicted, (2) the motions of particles falling through liquids predicted using IB-FDLBM quantitatively agree well with those obtained using immersed boundary methods based on the lattice Boltzmann equation (IB-LBM), (3) the developed method well predicts the interaction between two particles falling through a liquid, e.g., the drafting-kissing-tumbling motion, and (4) distortion of velocity fields in circular Couette flows at high relaxation times is removed by the additional collision term.一般社団法人 日本機械学会, 2011年12月, Journal of Fluid Science and Technology, 6(6), 1051-1064(6) (6), 1051 - 1064, 英語[査読有り]研究論文（学術雑誌）
- 2011年10月, Proc. of 1st International Symposium on Multiscale Multiphase Process Engineering (MMPE), CD-ROM (6 pages)Multiscale Simulation of Dispersed Flow in Bubble Columns研究論文（国際会議プロシーディングス）
- 2011年09月, International Conference on Advanced Technology in Experimental Mechanics 2011Mass transfer from a dissolving carbon dioxide bubble in glycerol-water solution研究論文（国際会議プロシーディングス）
- 2011年09月, International Conference on Advanced Technology in Experimental Mechanics 2011Assessment of turbulence models for bubbly flow in a bubble column研究論文（国際会議プロシーディングス）
- 2011年06月, 8th International Conference on COMPUTATIONAL FLUID DYNAMICS in the Oil & Gas, Metallurgical and Process IndustriesInterface tracking and multi-fluid simulation of bubbly flows in bubble columns研究論文（国際会議プロシーディングス）
- 2011年05月, 49th European Two-Phase Flow Group MeetingInterface tracking simulation of dissolving bubbles and of contaminated bubbles and drops研究論文（国際会議プロシーディングス）
- A scaling analysis based on the field equations for two phases and the jump conditions at the interface is carried out to deduce a balance of forces acting on a Taylor drop rising through stagnant liquid in a vertical pipe. The force balance is utilized to deduce a functional form of an empirical correlation of terminal velocity of the Taylor drop. Undetermined coefficients in the correlation are evaluated by making use of available correlations for two limiting cases, i.e. extremely high and low Reynolds number Taylor bubbles in large pipes. Terminal velocity data obtained by interface tracking simulations are also used to determine the coefficients. The proposed correlation expresses the Froude number Fr as a function of the drop Reynolds number Re-D, the Eotvos number Eo(D) and the viscosity ratio mu*. Comparisons between the correlation, simulations and experimental data confirm that the proposed correlation is applicable to Taylor drops under various conditions, i.e., 0.002 < Re-D < 4960, 4.8 < Eo(D) < 228, 0 <= mu* <= 70, 1 < N < 14700, -12 < log M < 4, and d/D < 1.6, where N is the inverse viscosity number, M the Morton number, d the sphere-volume equivalent drop diameter and D the pipe diameter. (C) 2010 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2011年04月, INTERNATIONAL JOURNAL OF MULTIPHASE FLOW, 37(3) (3), 241 - 251, 英語[査読有り]研究論文（学術雑誌）
- MASS TRANSFER FROM A BUBBLE IN A VERTICAL PIPEMass transfer from single carbon dioxide bubbles in a vertical pipe is measured using a stereoscopic image processing method to develop a mass transfer correlation applicable to a wide range of bubble and pipe diameters. The pipe diameters are 12.5, 18.2 and 25 0 mm and the bubble diameter ranges from 5 to 26 mm The ratio, lambda, of bubble diameter to pipe diameter is therefore varied from 0.2 to 1.8, which covers various bubble shapes such as spherical, oblate spheroidal, wobbling, cap, and Taylor bubbles. Measured Sherwood numbers, Sh, strongly depend on bubble shape, i.e., Sh of Taylor bubbles clearly differs from those of spheroidal and wobbling bubbles. Hence two Sherwood number correlations, which are functions of the Peclet number and the diameter ratio lambda, are deduced from the experimental data: one is for small bubbles (lambda < 0.6) and the other for Taylor bubbles (lambda > 0.6). The applicability of the proposed correlations for the prediction of bubble dissolution process is examined through comparisons between measured and predicted long-term bubble dissolution processes. The predictions are carried out by taking into account the presence of all the gas components in the system of concern, i.e. nitrogen, oxygen and carbon dioxide. As a result, good agreements for the dissolution processes for various bubble sizes and pipe diameters are obtained. It is also demonstrated that it is possible to evaluate an equilibrium bubble diameter and instantaneous volume concentration of carbon dioxide in a bubble using a simple model based on a conservation of gas components.AMER SOC MECHANICAL ENGINEERS, 2011年, PROCEEDINGS OF THE ASME/JSME 8TH THERMAL ENGINEERING JOINT CONFERENCE 2011, VOL 1 PTS A AND B, AJTEC2011-44089, 9 pages, 1363 - +, 英語[査読有り]研究論文（国際会議プロシーディングス）
- An interface tracking method for predicting motions of bubbles contaminated with soluble surfactants is presented. A level set method is utilized to track the interface. Transportations of surfactants in the bulk liquid and those at the interface are taken into account. The amount of adsorption and desorption is evaluated by using the Frumkin & Levich model. Simulations of bubbles contaminated with soluble surfactants are carried out, i.e., single air bubbles rising through stagnant water, Taylor bubbles in a vertical pipe filled with water, and a wobbling bubble in a vertical duct. As a result, the following conclusions are obtained: (1) the increase in drag coefficients of spherical bubbles due to the presence of surfactant, i.e. Marangoni effect, is well predicted, (2) surfactants mainly accumulate at the rear edge of a Taylor bubble and the Marangoni effect is very small in the nose region at high Eötvös and low Morton numbers, and therefore, the effects of surfactant on the bubble rising velocity are small in low viscosity systems, and (3) the surfactant concentration is low in the top region of a wobbling bubble, whereas it is high in the bottom region. The peak concentration appears at the side edge of the bubble and the location of the peak concentration moves with the bubble and wake movements. Copyright © 2011 by ASME.2011年, ASME-JSME-KSME 2011 Joint Fluids Engineering Conference, AJK 2011, 1(D) (D), 1641 - 1650, 英語[査読有り]研究論文（国際会議プロシーディングス）
- An interface tracking method for predicting bubble dissolution process is proposed. A non-diffusive scheme for advecting species concentrations is adopted to accurately compute the volume change due to mass transfer. The applicability of the proposed method is examined through several benchmark tests, i.e. mass transfer from a static bubble and that from free rising bubbles. Predicted species concentration distributions and mass transfer coefficients agree well with theoretical and empirical models. Dissolution of single carbon dioxide bubbles in a vertical pipe filled with water is also simulated. The bubbles consist only of carbon dioxide, and nitrogen and oxygen are initially dissolved in water. The volume change due to dissolution of carbon dioxide from the bubbles and evaporation of nitrogen and oxygen from water are well predicted. such as mass transfer from a stationary bubble, that from free rising bubbles, and single dissolving carbon dioxide bubbles in a vertical pipe filled with water. Copyright © 2011 by ASME.2011年, ASME-JSME-KSME 2011 Joint Fluids Engineering Conference, AJK 2011, 1(D) (D), 1593 - 1601, 英語[査読有り]研究論文（国際会議プロシーディングス）
- Three-dimensional interface tracking simulations were carried out to investigate the role of surface tension force in the process of vapor bubble lift-off from a hydrophilic heated surface in nucleate boiling. Since bubbles are frequently flattened along the heated surface in photographic experiments reported in literature, a bubble was assumed to be spheroidal in shape in the initial condition. The effect of phase change at the bubble interface was not taken into consideration for the sake of simplicity. In the present numerical simulations, the initially spheroidal bubble approached the spherical shape due to the surface tension force and was eventually lifted off the surface. The change in bubble shape induced local liquid flow directing toward the bubble base, that was the direct cause of the occurrence of the bubble lift-off. The dependence of the bubble migration velocity on several important parameters including the bubble size, surface tension coefficient and the density of surrounding liquid was also investigated. The change in bubble shape from flattened to more rounded causes the reduction of the surface energy, while the formation of local liquid flow leads to an increase in the kinetic energy. It was demonstrated that the bubble migration velocity after the lift-off can successfully be interpreted from the standpoint of energy conservation during the lift-off process. (c) 2010 Elsevier B.V. All rights reserved.ELSEVIER SCIENCE SA, 2010年12月, NUCLEAR ENGINEERING AND DESIGN, 240(12) (12), 3942 - 3948, 英語[査読有り]研究論文（学術雑誌）
- This paper investigates flow patterns encountered in two-phase swirling flow in a vertical pipe. Flows in the pipe are observed by using a high-speed video camera and the void fractions are calculated from the video images. Flows at various gas and liquid volumetric fluxes, JG and JL, are tested, i.e., 0.01 ≤ JG ≤ 8 m/s and 0.10 ≤ JL ≤ 0.80 m/s. When JL is higher than a certain value, no bubbles exist in a swirling flow but a gas column forms due to the centrifugal force even at low JG. The diameter of a gas column increases with JG at a constant JL. A tornado-shaped gas column forms at intermediate JG. At high JG, a liquid film including small bubbles and intermittent liquid slugs with many bubbles are observed. These flows are classified into four flow patterns based on the observations and the characteristics of void fraction. The statistical quantities of void fraction, i.e., the mean, standard deviation, skewness and kurtosis, are also computed. The validity of the flow pattern classification based on the statistical quantities has been confirmed.日本混相流学会, 2010年09月, 混相流, Vol. 24, No. 3, pp. 305-312(3) (3), 305 - 312, 日本語[査読有り]研究論文（学術雑誌）
- 2010年09月, 8th Japan-UK Seminar on Multiphase FlowsInterface Tracking simulation of mass transfer from a dissolving bubble in a vertical pipe研究論文（国際会議プロシーディングス）
- 2010年06月, 48th European Two-Phase Flow Group Meeting 2010, ETPFGM'10Agitation of thermally stratified water layer by free rising wobbling air bubble研究論文（国際会議プロシーディングス）
- 2010年05月, Proc. of 7th International Conference on Multiphase Flow, ICMF2010, 13 pagesTerminal velocity of a Taylor drop in a vertical pipe[査読有り]研究論文（国際会議プロシーディングス）
- 2010年05月, Proc. of International Conference on Multiphase Flow, ICMF2010, 11 pagesTerminal velocity of a single drop in a vertical pipe in clean and fully-contaminated systems[査読有り]研究論文（国際会議プロシーディングス）
- 2010年05月, Proceedings of 7th International Conference on Multiphase Flow, No. 2-1-2, 英語LES-Based Bubble Tracking Simulation of Turbulent Cavitation Flow in Injector Nozzle[査読有り]研究論文（国際会議プロシーディングス）
- 2010年03月, 混相流,, Vol. 23, No. 5, pp. 547-554, 日本語[査読有り]研究論文（学術雑誌）
- 神戸市立工業高等専門学校, 2010年03月, 神戸市立工業高等専門学校研究紀要, No. 48, pp. 45-50(48) (48), 45 - 50, 日本語非構造格子に基づく流体粒子数値予測手法に関する研究[査読有り]研究論文（学術雑誌）
- 2010年03月, 神戸高専研究紀要, 49, 日本語活性汚泥の減圧浮上に関する研究[査読有り]研究論文（学術雑誌）
- 神戸市立工業高等専門学校, 2010年03月, 神戸市立工業高等専門学校研究紀要, No. 48, pp. 39-44(48) (48), 39 - 44, 日本語T字形エルボによるキャビテーション壊食対策に関する研究[査読有り]研究論文（学術雑誌）
- Interface tracking simulations of single drops rising through a vertical pipe are carried out using three coordinate systems, i.e. cylindrical, general curvilinear and Cartesian coordinates, to investigate the effects of coordinate system and spatial resolution on the accuracy of predictions. Experiments of single drops in a vertical pipe are also conducted to obtain experimental data for comparisons with simulations. The drop shape observed are spheroidal and deformed spheroidal at low values of the diameter ratio, λ, of the sphere-volume equivalent diameter of a drop to the pipe diameter, whereas they take bullet-shapes at large λ. The conclusions obtained are as follows: (1) the effects of coordinate system on drop shape are small at low λ. At large λ, the effects are also small for drops in a low viscosity system, whereas non-physical shape distortion takes place when the Cartesian coordinates are used with low spatial resolution for drops in a high viscosity system, and (2) the drop terminal velocity and the velocity profile in the liquid film between a bullet-shaped drop and a pipe wall are well predicted using all the coordinate systems tested even at low spatial resolution.2010年03月, Journal of Computational Multiphase Flows, 2(1) (1), 47 - 57, 英語[査読有り]研究論文（学術雑誌）
- 2010年, The 8th Japan-UK Seminar on Multiphase FlowBubble Coalescence in a Slurry Bubble Column研究論文（国際会議プロシーディングス）
- スラリー気泡塔設計のための進化する数値計算手法の構築A numerical method for simulating gas-liquid-solid three-phase flows in a slurry bubble column reactor was developed based on a mixing model and an artificial neural network (ANN). Macroscopic variables required for the mathematical closure of the mixing model, such as the gas holdup, dispersion coefficient, and reaction rate, were evaluated by use of ANN, rather than empirical correlations. In contrast to the empirical correlations, ANN possesses the following advantages: (1) it is applicable to a wide range of flow conditions owing to its generalization ability, and (2) the accuracy is easily improved through a learning process using an up-to-date database. For the validation of ANN, gas holdups in an air-water bubble column were measured. Then simulations of gas-liquid-solid three-phase flows in a Fischer-Tropsch (FT) reactor were carried out using the combination of the mixing model and ANN to demonstrate its potential. The results obtained are as follows: (1) ANN accurately evaluates the gas holdups in the air-water bubble column, (2) the combination of the mixing model and ANN gives good predictions for CO conversions in FT synthesis, and (3) the accuracy of ANN is easily improved through the re-learning process. These results imply that the proposed method can be a framework for simulating multiphase flows in large industrial systems.SOC CHEMICAL ENG JAPAN, 2010年, 化学工学論文集, 36(1) (1), 17 - 24, 日本語[査読有り]研究論文（学術雑誌）
- An empirical correlation of the terminal velocity of a Taylor bubble in a vertical pipe is proposed. A fundamental functional form, Fr = f (Re D, EoD), of the correlation is deduced by carrying out a dimensional analysis based on the local instantaneous field equations and the jump conditions. Here Fr is the Froude number, ReD the bubble Reynolds number, and EoD the Eötvös number. In the two limiting cases, (EoD →∞ and ReD →∞) and (EoD →∞ and ReD → 0), the deduced functional form approaches those of the well-known Fr models. Coefficients appearing in the correlation are determined by making use of the limiting cases and available experimental data. Comparisons between the proposed Fr correlation and the experimental data show that the correlation gives a good estimation of terminal velocities of Taylor bubbles for a wide range of fluid properties and pipe diameters, i.e., 10-7 < ReD < 104, 4 < EoD < 3 × 103, 10-2 < N < 105, and -11 < log M < 10, where N is the inverse viscosity number and M the Morton number. © 2010 by Begell House, Inc.2010年, Multiphase Science and Technology, 22(3) (3), 197 - 210, 英語[査読有り]研究論文（国際会議プロシーディングス）
- Bubble coalescence and breakup models are implemented into the (N+2)-field model, which is a hybrid combination of a multi-fluid model and an interface tracking method to examine the applicability of available bubble coalescence and breakup models to poly-dispersed bubbly flows in bubble columns. The implemented models are a breakup model proposed by Luo & Svendsen and coalescence models proposed by Prince & Blanch and Wang et al. The validation of these models is carried out using experimental data for a coalescence-dominating flow, a breakup-dominating flow, a flow in-between, and a flow with low coalescence-breakup frequencies. Good predictions of void and bubble size distributions for the four kinds of poly-dispersed bubbly flows in bubble columns are obtained. It is also demonstrated that all the possible bubble coalescence mechanisms, i.e. bubble coalescence caused by bubble collisions due to turbulent eddies, the difference in bubble velocities and shear slip, and bubble coalescence due to wake entrainment, should be taken into account for obtaining good predictions.一般社団法人 日本機械学会, 2009年11月, 日本機械学会論文集(B編), 75-759, pp.2182-2189(759) (759), 2182 - 2189, 日本語[査読有り]研究論文（学術雑誌）
- Numerical simulations based on a three-dimensional two-way bubble tracking method are carried out to predict bubble motions in a square duct with an obstacle and in a two-by-three rod bundle with a grid spacer. Comparisons between measured and predicted bubble motions demonstrate that the two-way bubble tracking method gives good predictions for trajectories of small bubbles in the upstream side of the grid spacer in the rod bundle geometry. The predicted bubble trajectories clearly show that bubbles are apt to migrate toward the rod surface in the vicinity of the bottom of the grid spacer. Analysis of forces acting on the bubbles confirms that pressure gradient force induced by the presence of the spacer is the main cause of the bubble lateral migration toward the rod surface. Motions of steam bubbles at a nominal operating condition of a pressurized water reactor (PWR) are also predicted by using the bubble tracking method, which indicates that steam bubbles also migrate toward the rod surface at the upstream side of the spacer due to the spacer-induced pressure gradient force.The Japan Society of Mechanical Engineers, 2009年09月, Journal of Power and Energy Systems, Vol. 3, No. 2, pp. 393 - 404(2) (2), 393 - 404, 英語[査読有り]研究論文（学術雑誌）
- 2009年09月, 5th EuropeanTwo-Phase Flow Group Meeting, CD-ROM, 9 pages, 英語Dissolution of Carbon Dioxide Bubbles in Clean and Contaminated Systems[査読有り]研究論文（国際会議プロシーディングス）
- 2009年09月, 5th EuropeanTwo-Phase Flow Group Meeting, CD-ROM, 6 pages, 英語Dimensional Analysis of Terminal Velocityof a Taylor Bubble in a Vertical Pipe[査読有り]研究論文（国際会議プロシーディングス）
- A drag correlation for a fluid particle rising along the axis of a vertical pipe at low and intermediate Reynolds numbers, Re, is proposed by making use of available correlations and a numerical database accumulated by interface tracking simulations. The accuracy of the interface tracking method has been verified through comparisons between measured and predicted velocities of single drops in vertical pipes. Being similar to drag model for solid spheres proposed by Michaelicles, the developed drag correlation takes into account inertial and wall effects as their linear combination. The correlation gives good estimation of the drag coefficient for fluid particles rising through stagnant liquids in vertical pipes under the conditions of 0.13 <= Eo <= 30, -10.0 <= logM <= 2.0, 0.083 <= Re < 200, 0 <= kappa <= 10.0 and lambda <= 0.6, where Eo is the Eotvos number, M the Morton number, K the viscosity ratio and A the ratio of particle diameter to pipe diameter. (C) 2009 Elsevier Ltd. All rights reserved.PERGAMON-ELSEVIER SCIENCE LTD, 2009年06月, CHEMICAL ENGINEERING SCIENCE, 64(12) (12), 3019 - 3028, 英語[査読有り]研究論文（学術雑誌）
- 神戸市立工業高等専門学校, 2009年03月, 神戸市立工業高等専門学校研究紀要, No. 47, pp. 13-18(47) (47), 13 - 18, 日本語水平管内におけるマイクロバブルによる摩擦圧力損失低減に関する研究[査読有り]研究論文（学術雑誌）
- 神戸市立工業高等専門学校, 2009年03月, 神戸市立工業高等専門学校研究紀要, No. 47, pp. 19-24(47) (47), 19 - 24, 日本語鉛直管内旋回気液二相流の流動様式とその判別法[査読有り]研究論文（学術雑誌）
- The purpose of this study is to develop a practical numerical method for predicting multiphase flows in industrial systems of large scale. The practical method should be equipped with two functions. One is the ability of dealing with complex geometry of a given practical system without time-consuming grid generation procedures. The other is to simulate flows in large complex systems without using a high number of computational cells. Therefore, a numerical method based on the combination of a multifluid and interface tracking method and an immersed boundary method is proposed in this study. The hybrid combination of the multifluid and interface tracking methods enables us to simulate multiphase flows with various scales and various phases, and the immersed boundary method makes it possible to satisfy the two requirements. Several numerical simulations are carried out to demonstrate the potential of the proposed method. Comparisons between measured and predicted bubbly flows around single obstacles prove that the proposed method gives reasonable predictions for the interaction between bubbly flow and structures. Simulation of a bubble column with complex structure demonstrates its applicability to large industrial systems.Begell House Inc., 2009年, Multiphase Science and Technology, 21(1-2) (1-2), 141 - 155, 英語[査読有り]研究論文（学術雑誌）
- 2008年07月, Proc. 5th International Conference on Transport Phenomena in Multiphase Systems, CD-ROM, 英語Single Drops Rising Through Stagnant Liquids in a Vertical Pipe[査読有り]研究論文（国際会議プロシーディングス）
- T字形エルボにおけるキャビテーションの流動特性に関する研究Water is more liable to evaporate than oils since its saturated vapor pressure is higher than the ones of oils. Therefore it is important to prevent erosion due to cavity in hydraulic devices using water. In our previous study, we reported that erosion at elbow of pipeline can be reduced by using T-shape elbow. However its mechanism has not been clarified yet. In the present study, we observed cavitating flows in T-shape elbows with different shapes using a digital and high-speed video cameras. The amount of cavity was evaluated from taking images. As a result we confirmed that (1) circulating flow is steadily formed when the length of cavity-trapping region L ≥ 22.5 mm. If L ≤ 15 mm, circulating flow is formed intermittently, or it is not formed, (2) The amount of cavity in cavity-trapping region increases with L, and (3)there are two uses (a) T-shape elbows can be used to guide cavity into the region, in which circulating flow is formed, for defending the downstream of pipeline, and (b) Then the elbows can also be used to guide bubbles to the downstream so as to prevent erosion at the elbows.神戸市立工業高等専門学校, 2008年03月, 神戸市立工業高等専門学校研究紀要, No. 46, pp. 13-18, 13 - 18, 日本語[査読有り]研究論文（学術雑誌）
- 2008年, Proc. 6th Japan-Korea Symp. Nuclear Thermal Hydraulics and Safety, CD-ROM, pp.1-8, 英語A Numericsl Study on the Mechanism of Vapor Bubble Lift-off from a Flat Plate[査読有り]研究論文（国際会議プロシーディングス）
- The number of computational cells assigned to a small bubble is liable to be small in a volume tracking simulation of a poly-dispersed bubbly flow. The volume tracking method is, therefore, expected to at least give qualitatively reasonable predictions of bubble motion without assigning many cells. In this study, bubble motions in stagnant water and shear flows are simulated with low spatial resolutions using a volume tracking method proposed in our previous study. As a result, we confirm that (1) the method can yield correct characteristics of rising velocity and lateral motion of bubbles, and (2) the bubble volume and interface sharpness are kept well even after a bubble had experienced a large deformation in a strong shear flow.The Japan Society of Mechanical Engineers, 2007年12月, Journal of Fluid Science and Technology, Vol. 2, No. 2, pp. 490-501(2) (2), 490 - 501, 英語[査読有り]研究論文（学術雑誌）
- The interface tracking method based on a nonuniform subcell scheme (NSS) and a level set function (LS) is applied to three-dimensional (3-D) general curvilinear and cylindrical coordinate systems. The method for 3-D general curvilinear coordinates is based on the combination of NSS, LS, a finite volume method, and a multiblock method. The method for cylindrical coordinates is based on the combination of NSS, LS, the variable transformation proposed by Verzicco and Orlandi (1996), and a simple method for evaluating surface tension force at the singular point r 0. Simulations of single air and vapor bubbles in rod bundle geometry are carried out using the former method. Predicted bubble shapes and terminal velocities agree well with measured data and an available terminal velocity model. Bubbles and drops in pipes are simulated using the latter method. Terminal velocities of single drops confined in a pipe are well predicted. The bouncing motion of an air bubble in a laminar water flow in a vertical pipe is successfully simulated. Volumes of bubbles and drops are well conserved in all simulations.2007年, Multiphase Science and Technology, 19(2) (2), 121 - 140, 英語[査読有り]研究論文（学術雑誌）
- Numerical simulation of single drops in a vertical pipe using various coordinate systemsThe effects of the diameter ratio lambda (=d/D, where d and D are the diameters of a drop and a pipe, respectively), the Morton number M and the viscosity ratio kappa (= mu(d)/mu(c), where mu is the viscosity and the Subscripts d and c are the dispersed fluid particle and the continuous phase, respectively) on terminal velocities and shapes of single drops rising through stagnant liquids in a vertical pipe are investigated experimentally. Then, the drops in the pipe are simulated using a volume tracking method with various coordinate systems, i.e., three-dimensional (3D) cylindrical coordinates, 3D general curvilinear coordinates and 3D Cartesian coordinates. Predicted velocities and shapes of the drops using three coordinate systems are compared with the measured data to examine the effects of coordinate systems on the accuracy of prediction. As a result. (1) The velocity ratio V-T/V-T0 (V-T and V-T0 are the terminal velocity in a pipe and infinite liquid, respectively) decreases as lambda increases, and it depends not only on lambda but also on M and kappa, (2) Good predictions for the terminal velocities and shapes of drops Lire obtained not only with cylindrical coordinates and Curvilinear coordinates but also with Cartesian coordinates, provided that the spatial resolution is high, (3) When the spatial resolution is low, effects of coordinate systems on a drop shape are larger for Cartesian coordinate systems than for cylindrical coordinate and general Curvilinear coordinate systems, and (4) Errors in predicted drop velocities are not so large even with very low spatial resolution.AMER SOC MECHANICAL ENGINEERS, 2007年, FEDSM 2007: PROCEEDINGS OF THE 5TH JOINT ASME/JSME FLUIDS ENGINEERING SUMMER CONFERENCE, VOL 2, PTS A AND B, CD-ROM, FEDSM2007-37491, 291 - 296, 英語[査読有り]研究論文（国際会議プロシーディングス）
- 気泡運動の低空間分解能体積追跡計算の定性的妥当性The number of computational cells assigned to a small bubble is liable to be small in a volume tracking simulation of a poly-dispersed bubbly flow. The volume tracking method is, therefore, expected to at least give qualitatively reasonable predictions of bubble motion without assigning many cells. In this study, bubble motions in stagnant water and shear flows were simulated with low spatial resolutions using a volume tracking method proposed in our previous study. As a result, we confirmed that (1) the method could yield correct characteristics of rising velocity and lateral motion of bubbles, and (2) the bubble volume and interface sharpness were kept well even after a bubble had experienced a large deformation in a strong shear flow.一般社団法人日本機械学会, 2006年10月, 日本機械学会論文集, Vol. 72, No. 722, (B), pp. 241(722) (722), 2410 - 2417, 日本語[査読有り]研究論文（学術雑誌）
- 2006年09月, 混相流, 20巻3号, pp.244-251, 日本語[査読有り]研究論文（学術雑誌）
- 2006年09月, 4th Japanese-European Two-Phase Flow Group Meeting, CD-ROM (8 pages)Interface Tracking Simulation of Bubbles and Drops in Complex Geometries[査読有り]研究論文（国際会議プロシーディングス）
- 2006年07月, Computational Fluid Dynamics Journal, Vol. 15, Iss. 2, pp. 225-232, 英語A Volume Tracking Method Based on Non-Uniform Subcells and Continuum Surface Force Model Using a Local Level Set Function[査読有り]研究論文（学術雑誌）
- Computational Multi-Fluid Dynamics, CMFD, has come to play an important role in practical engineering and fundamental researches on multiphase flows. Since a multiphase flow is often comprised of interdependent elementary phenomena with length scales ranging from mesoscale to macroscale, a multipurpose CMFD should cover a wide range of scales by the combination of different numerical methods such as interface tracking, bubble tracking and averaging methods. This report reviews recent activities on CMFD for bubbly flows, in particular, an improvement of an interface tracking method based on volume tracking, two types of hybrid CMFD methods based on a combination of interface tracking and bubble tracking methods and that of interface tracking and averaging methods, and work in progress heading toward the hybrid combination of the three methods.2006年, Multiphase Science and Technology, 18(1) (1), 73 - 110, 英語[査読有り]研究論文（学術雑誌）
- 2005年09月, Proc. Int. Minisymposium on Challenger and Adavances in Flow SImulation and Modeling, p.276, 英語Toward Hybrid Simulation of Bubbly Flow (Invited Lecture)研究論文（国際会議プロシーディングス）
- A volume tracking method based on Advanced Subgrid Counting AlgorithmA volume tracking scheme, ASCA (Advanced Subgrid Counting Algorithm), which is easily extendable to three-dimensions and possesses a good volume conservation property, is proposed. To examine the potential of ASCA, several two-phase flow simulations are carried out. As a result, the following conclusions are obtained: (1) Predicted shapes and breakup characteristics of single fluid particles in simple shear flows agree well with available data even with a low spatial resolution, (2) A water drop impinging on water surface and an air bubble rising through a stagnant water are successfully simulated with little errors in volume conservation, (3) Predicted shapes and terminal velocities of single drops in stagnant liquids under wide ranges of the viscosity ratio and Morton number agree well with measured one, and (4) In spite of a low spatial resolution, vapor bubbles flowing in a strong shear flow are simulated with good volume conservation. The interface sharpness is well preserved even after a large deformation of bubbles.AMER SOC MECHANICAL ENGINEERS, 2005年, Proceedings of the ASME Fluids Engineering Division Summer Conference, Vol 2, CD-ROM, FEDSM2005-77396, 265 - 270, 英語[査読有り]研究論文（国際会議プロシーディングス）
- 2004年, Proc.6th US-Japan Seminar on Two-Phase Flow Dynamics, 1,1-16, 英語Modelong and Hybrid Simulation of Bubbly Flow(Keynote Lecture)[査読有り]研究論文（国際会議プロシーディングス）
- The validity of the Simple Counting Algorithm (SCA) for the volume tracking simulations of a single bubble in stagnant liquid is examined. As a result, it is confirmed that (i) the SCA can keep sharp interface, but (ii) the volume of bubble is increased. An improved SCA is proposed in order to reduce the volume error. Numerical simulations of a single air bubble rising in stagnant water are carried out using a course grid. i.e. only 6 cells to the bubble equivalent diameter. The volume of bubble is conserved well and a sharp interface is kept for as long as 12 seconds.Japan Society of Mechanical Engineers, 2004年, 日本機械学会論文集（Ｂ編）, 70(698) (698), 2538 - 2544, 英語[査読有り]研究論文（学術雑誌）
- 2003年07月, Proc. 4th ASMER-JSME Joint Fluids Engineering Conf., CD-ROM, 英語Evaluation of Volume Tracking Algorithms for Gas-Liquid Two-Phase Flows.[査読有り]研究論文（国際会議プロシーディングス）
- Numerical analysis of methanol droplet combustion in hot air flowNumerical analysis on evaporation and combustion phenomena of a deforming methanol droplet in a hot air flow was carried out in the present study. An interface tracking method, the Volume of Fluid (VOF) method, was employed to predict the droplet deformation. This approach holds a high potential to simulate the spray combustion including many deforming droplets moving freely in a gas flow using a Eulerian mesh. The advance of the gas-liquid interface in time was predicted by the CIP scheme to avoid numerical diffusion of gas-liquid interface. The numerical calculation of the evaporation and ignition of a single methanol droplet in a hot stagnant air was conducted to clarify the validity of the numerical method. Then the numerical simulation of the ignition and transient flame behavior of a deforming droplet in convective hot air flow was carried out, and the envelope flame developed around a fuel droplet was reasonably predicted. The effects of numerical scheme on calculated results were also examined.SOUTHEAST UNIV PRESS, 2003年, COMBUSTION SCIENCE AND TECHNOLOGY IN ASIA-PACIFIC AREA: TODAY AND TOMORROW, Paper No.071, 509 - 512, 英語[査読有り]研究論文（国際会議プロシーディングス）
- Numerical analysis of ignition and combustion of a deformable dropletA numerical analysis was conducted to predict the evaporation and combustion phenomenon of droplets with free moving interface. An interface tracking method, the Volume of Fluid (VOF) method, was assigned to predict the droplet motion, which has a high potential to simulate the spray phenomena including many deforming droplets moving freely with each velocity in the surrounding gases. A one-step global reaction of methanol vapor in air was assigned to predict chemical reaction, and the reaction rate was given by the Arrhenius law at the first step. The numerical simulation of the evaporation and ignition phenomena of a single methanol droplet in a hot stagnant air was conducted to clarify the validity of the numerical method. Then the numerical simulation of the ignition of a methanol droplet in hot air flow was carried out, and the envelope flame developed around a droplet in a hot airflow was predicted. The effects of evaporation on droplet deformation and the ones of deformation of a droplet on the combustion phenomena around a droplet were numerically studied.2002年, International Symposium on Combustion Abstracts of Works-in-Progress Posters, 397, 英語研究論文（国際会議プロシーディングス）
- 日本混相流学会, 2016年09月, 混相流 = Japanese journal of multiphase flow, 30(3) (3), 327 - 329, 日本語第41回レクチャーシリーズ「混相流の数値シミュレーション」開催報告
- B224 気水分離器における液相分離構造が旋回二相流動に及ぼす影響(OS3 軽水炉・新型炉・原子力安全(4))A downscaled model of a steam separator is used to understand characteristics of swirling flow from the viewpoint of the effects of flow restriction ring (FRR) which is installed in a drainage channel. Main conclusions obtained are as follows. (1) Carryunder performance, the weight percentage of entrained gas in the water leaving the separator, improves by narrowing the flow area in the drainage channel using the FRR because it becomes hard to exhaust the gas. On the other hand, (2) the ratio of the separated flow rate to the liquid flow rate may decrease because the liquid phase climbing over the 1st pick-off ring and scattered droplets increase. Therefore, (3) it is necessary to optimize the FRR design based on the gas-liquid separation performance and pressure-loss performance.一般社団法人日本機械学会, 2015年06月07日, 動力・エネルギー技術の最前線講演論文集 : シンポジウム, 2015(20) (20), 239 - 240, 日本語
- 日本流体力学会, 2014年12月, ながれ, 33(6) (6), 509 - 512, 日本語時空間フィルタ流速計による単一液滴に作用する界面せん断応力の評価[査読有り][招待有り]速報，短報，研究ノート等（学術雑誌）
- GS35 埋め込み境界を適用した格子ボルツマン法における流速の滑りの除去We verify the effectiveness of the two-relaxation-time (TRT) collision operator in reducing boundary slip computed by the immersed boundary-lattice Boltzmann method (IB-LBM). Rigorous analysis of the symmetric Poiseuille flows reveals that the relaxation time for the antisymmetric part controls the boundary velocity and the boundary slip velocity computed by the IB-LBM. In the simulation of the cylindrical Couette flows, the boundary slip caused by the use of the single relaxation time (SRT) model in the IB-LBM impairs the accuracy that the LBM intrinsically possesses. The IB-LBMs with the TRT exhibit second-order accuracy in space for both the direct forcing method and the implicit correction method.一般社団法人日本機械学会, 2014年10月25日, 流体工学部門講演会講演論文集, 2014, "GS35 - 1"-"GS35-2", 日本語
- 日本計算工学会, 2014年06月, 計算工学講演会論文集 Proceedings of the Conference on Computational Engineering and Science, 19, 6p, 日本語二緩和時間衝突則に基づく埋め込み境界法を適用した格子ボルツマン法の精度評価
- 一般社団法人日本機械学会, 2014年03月18日, 関西支部講演会講演論文集, 2014(89) (89), "8 - 11", 日本語807 散気管内外の流れに及ぼす液相粘度の影響(OS-3 流体数値計算と実験計測の展開および産業界への応用(2))
- 一般社団法人日本機械学会, 2014年03月18日, 関西支部講演会講演論文集, 2014(89) (89), "8 - 27", 日本語823 一様せん断流中単一気泡に働く揚力に関する研究 : 第2報 数値計算(GS-5 気泡)
- 一般社団法人日本機械学会, 2014年03月18日, 関西支部講演会講演論文集, 2014(89) (89), "8 - 26", 日本語822 一様せん断流中単一気泡に働く揚力に関する研究 : 第1報 実験(GS-5 気泡)
- B121 4x4ロッドバンドル内気液二相気泡流に関する研究(OS9 熱・流動)Although area-averaged void fractions in subchannels have been measured frequently in literatures, distribution of local void fraction in a rod bundle have been rarely reported. Little is also known on bubble motion in a rod bundle. Distributions of local void fraction and bubble velocity in a 4x4 rod bundle were, therefore, measured using a conductance probe in this study. Liquid velocity was also measured using LDV. Experimental results showed that the distributions of the void fractions in the inner and side subchannels depend not only on lift force acting on bubbles but also on geometrical interaction between the bubble and the rod wall. The relative velocity between bubbles and liquid phase in the subchannel possesses non-uniform distribution over the cross-section. This is caused by increase in drag force acting on bubbles due to bubble deformation induced by interaction between the bubbles and the rod wall.一般社団法人日本機械学会, 2012年06月20日, 動力・エネルギー技術の最前線講演論文集 : シンポジウム, 2012(17) (17), 61 - 62, 日本語
- 一般社団法人日本機械学会, 2011年03月19日, 関西支部講演会講演論文集, 2011(86) (86), "5 - 26", 日本語526 界面追跡法に基づく気液間物質輸送の数値計算(GS-5 二相流)
- 2009年12月16日, 微粒化シンポジウム講演論文集 = Symposium (ILASS-Japan) on Atomization, 18, 113 - 117, 日本語ノズル内キャビテーション乱流のLES
- G105 T字形エルボにおけるキャビテーションの流動(2)(GS1 混相流・振動流,一般セッション)Water is more liable to evaporate than oils since its saturated vapor pressure is higher than the ones of oils. Therefore it is important to prevent erosion due to cavitation bubbles in hydraulic devices using water. In our previous study, we reported that erosion at elbow of pipeline can be reduced by using T-shape elbow. However its mechanism has not been clarified yet. In the present study, we observed cavitating flows in T-shape elbows of different shapes using a digital and high-speed video cameras. The amount of cavity was evaluated from taking images. As a result we confirmed that (1) circulating flow is steadily formed when the length of cavity-trapping region L ≧ 22.5 mm. If L ≦ 15 mm, circulating flow is formed intermittently, or it is not formed, (2) The amount of cavitation bubbles in cavity-trapping region increases as L increases, and (3)there are two uses (a) T-shape elbows can be used to guide cavitation bubbles into the region, in which circulating flow is formed, for defending the downstream of pipeline. (b) Then the elbows can also be used to guide bubbles to the downstream so as to prevent erosion at the elbows.一般社団法人日本機械学会, 2007年11月17日, 流体工学部門講演会講演論文集, 2007, "G105 - 1"-"G105-4", 日本語
- 一般社団法人日本機械学会, 2005年11月17日, 計算力学講演会講演論文集, 2005(18) (18), 805 - 806, 日本語2504 表面張力評価へのカーン-ヒリヤード方程式の応用(OS-25A,OS-25 混相流および界面挙動の数値解析)
- 2005年, IEEE Trans.on Nuclear Science, 52(1),295-298, 英語Flow Visualization of Segregation Process in a Fluidized-Bed by Neutron Radiography[査読有り]その他
- 2004年, 日本混相流学会年会講演会講演論文集, 2004改良SCAによる流体粒子の界面追跡計算
- 一般社団法人日本機械学会, 2003年11月22日, 計算力学講演会講演論文集, 2003(16) (16), 111 - 112, 日本語208 簡素化界面再構築法に基づく界面追跡法
- 2003年, 日本混相流学会年会講演会講演論文集, 2003微粒化ノズル内キャビテーション流れと液体噴流の微粒化
- 2003年, 日本混相流学会年会講演会講演論文集, 2003液滴・噴霧燃焼の数値シミュレーション
- 2002年, 日本機械学会関西学生会卒業研究発表講演会講演前刷集, 2002界面追跡法の改良 新しい界面移流スキームの提案
- 分担執筆, 第１章, 化学工学会 関東支部, 2024年01月, 日本語, ISBN: 4866938943化学工学系流体シミュレーションの最前線 ―基礎・実践・将来展望― (最近の化学工学 72)
- 朝倉書店, 2020年04月, 日本語, ISBN: 4254237952流体力学 (機械工学基礎課程)
- その他, Chaps 1 and 2, World Scientific Publishing, 2018年05月, 英語, ISBN: 9789813229440Encyclopedia of Two-Phase Heat Transfer and Flow III: Macro and Micro Flow Boiling and Numerical Modeling Fundamentals学術書
- 第38回数値流体力学シンポジウム, 2024年12月, 日本語平板間静止液中単一気泡の運動に関する数値的研究口頭発表（一般）
- 第38回数値流体力学シンポジウム, 2024年12月, 日本語狭平板間静止液中を上昇する単一気泡の数値計算口頭発表（一般）
- 日本流体力学会年会2024, 2024年09月, 日本語狭平板間静止液中単一気泡の分裂に関する研究口頭発表（一般）
- 日本流体力学会年会2024, 2024年09月, 日本語鉛直細円柱群内を上昇する気泡の形状と終端速度
- 日本機械学会年次大会, 2024年09月, 日本語円管内旋回流による油水および固液分離に関する研究ポスター発表
- 混相流シンポジウム2024, 2024年09月, 日本語停滞水条件における3×3ロッドバンドル内での界面摩擦口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語鉛直円管内旋回環状流内体積率及び圧力損失の一次元三流体モデルによる予測口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語堆積層内気液二相流の流動特性に関する研究口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語気泡塔内ボイド率に及ぼす初期液位依存性に関する研究口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語気泡塔内非均質気泡流におけるドリフトフラックスパラメータの検討口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語Experimental analysis of bubble terminal velocity in fiber bundle口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語鉛直円管内テイラー気泡の液膜厚さに加振が及ぼす影響に関する研究ポスター発表
- 混相流シンポジウム2024, 2024年09月, 日本語平板間静止液中単一気泡の形状及び上昇速度に関する研究ポスター発表
- 混相流シンポジウム2024, 2024年09月, 日本語鉛直細円柱群内の気泡の運動に関する研究ポスター発表
- 混相流シンポジウム2024, 2024年09月, 日本語フラッディング状態での鉛直管内流動特性に対する管径の影響口頭発表（一般）
- 混相流シンポジウム2024, 2024年09月, 日本語微細流路内テイラー気泡の形状に及ぼす界面活性剤の影響ポスター発表
- 混相流シンポジウム2024, 2024年09月, 日本語平板間静止液中を上昇する単一気泡の形状及び抗力係数口頭発表（一般）
- 第67回理論応用力学講演会, 2024年09月, 日本語正方形マイクロチャネル内テイラー流の圧力損失と気泡速度に関する研究口頭発表（一般）
- 第37回数値流体力学シンポジウム, 2023年12月, 日本語脳動脈瘤内流れの血行力学パラメータの数値予測と機械学習による破裂リスク予測可能性の検討口頭発表（一般）
- IUTAM Symposium on Dynamics and Interface Phenomena of Bubbles and Droplets at Multiple Scales, 2023年12月, 英語Surface Coverage of Bubbles with Surfactant or Particles[招待有り]口頭発表（一般）
- 日本流体力学会年会2023, 2023年09月, 英語液滴に働く揚力の粘性比依存性に関する数値的研究口頭発表（一般）
- 混相流シンポジウム2023, 2023年08月, 日本語微細流路内気泡表面における界面活性剤吸着量に関する研究口頭発表（一般）
- 混相流シンポジウム2023, 2023年08月, 日本語微細流路を用いた粒子付着気泡の生成と抗力係数の評価に関する研究ポスター発表
- 混相流シンポジウム2023, 2023年08月, 日本語U字管内気液二相流のボイド率に及ぼすベンド曲率の影響ポスター発表
- 混相流シンポジウム2023, 2023年08月, 日本語球形気泡の抗力と表面渦度に及ぼす界面活性剤の影響口頭発表（一般）
- 混相流シンポジウム2023, 2023年08月, 日本語鉛直円管内旋回環状流における界面および壁面摩擦係数に及ぼす管径の影響に関する研究口頭発表（一般）
- ASME-JSME-KSME Fluids Engineering Division, AJKFED 2023, 2023年07月, 英語Effects of Fine Particles and Surfactant on Bubbly Flows in Narrow Rectangular Column口頭発表（一般）
- 11th International Conference on Multiphase Flow, ICMF-2023, 2023年04月, 英語Lift Correlations of Ellipsoidal Bubbles in Low and High Viscosity Liquids[招待有り]口頭発表（基調）
- 日本機械学会関西支部第98期定時総会講演会, 2023年03月, 日本語中空糸膜を有する塔内気泡流に関する研究口頭発表（一般）
- ANS Winter Meeting 2022, 2022年11月, 英語Effects of Upper and Lower Pipe Ends on Interfacial and Wall Friction Factors under Counter-Current Flow Limitation in a Vertical Pipe口頭発表（一般）
- ANS Winter Meeting 2022, 2022年11月, 英語Effects of Liquid Viscosity on Friction Factors of Gas-Liquid Swirling Annular Flows in a Vertical Pipe口頭発表（一般）
- 1st European-American-Japanese Two-Phase Flow Group Meeting (EAJTPFGM22), 2022年10月, 英語Fragmentation of Liquid Drop Falling through Miscible Liquid due to Rayleigh-Taylor Instabilityシンポジウム・ワークショップパネル（指名）
- 1st European-American-Japanese Two-Phase Flow Group Meeting (EAJTPFGM22), 2022年10月, 英語Bubbly Flows in a Column with Submerged Hollow Fiber Membranes口頭発表（一般）
- 日本流体力学会年会2022, 2022年09月, 日本語水流による管路内油性付着物の除去に及ぼす微細気泡の影響口頭発表（一般）
- 4th International Symposium on Multiscale Multiphase Process Engineering (MMPE), 2022年09月, 英語Pressure Drop of Single Bubble in Taylor Flows through Square Microchannels口頭発表（一般）
- 4th International Symposium on Multiscale Multiphase Process Engineering (MMPE), 2022年09月, 英語Bubbly Flows in a Vertical Pipe with Horizontal Oscillationポスター発表
- 混相流シンポジウム2022, 2022年08月, 日本語正方形マイクロチャネル内テイラー流における気泡の圧力損失に関する研究口頭発表（一般）
- 混相流シンポジウム2022, 2022年08月, 日本語狭平板間静止水中単一気泡の分裂に関する研究口頭発表（一般）
- 混相流シンポジウム2022, 2022年08月, 日本語鉛直円管内気液二相対向流の界面及び壁面摩擦係数に及ぼす管端形状の影響口頭発表（一般）
- 混相流シンポジウム2022, 2022年08月, 日本語鉛直円管内旋回環状流の界面及び壁面摩擦係数に及ぼす旋回強度及び液相粘度の影響口頭発表（一般）
- 日本流体力学会年会2021, 2021年09月, 日本語汚染系単一球形液滴の数値計算に関する研究口頭発表（一般）
- 混相流シンポジウム2021, 2021年08月, 日本語堆積層内における気液二相流に関する研究口頭発表（一般）
- 混相流シンポジウム2021, 2021年08月, 日本語振動を伴う鉛直円管内の気液二相流動特性に関する研究
- 混相流シンポジウム2021, 2021年08月, 日本語鉛直円管内気液二相旋回流の界面および壁面摩擦係数に及ぼす旋回強度の影響
- 混相流シンポジウム2021, 2021年08月, 日本語鉛直サーペンタイン管内気液二相流のボイド率及び液膜厚さ
- 混相流シンポジウム2021, 2021年08月, 日本語垂直U字管内気液二相流のボイド率に関する研究
- 混相流シンポジウム2021, 2021年08月, 日本語気泡流による中空糸膜の汚れ除去と揺動に関する研究
- 混相流シンポジウム2021, 2021年08月, 日本語気泡塔内平均及び局所ボイド率の初期液位依存性に関する研究
- 混相流シンポジウム2020, 2020年08月, 日本語鉛直円管内気液二相流に及ぼす管加速度の影響に関する研究口頭発表（一般）
- 混相流シンポジウム2020, 2020年08月, 日本語鉛直円管内旋回環状流の界面及び壁面摩擦係数に及ぼす液相粘度の影響口頭発表（一般）
- 混相流シンポジウム2020, 2020年08月, 日本語気泡群による液相中粒子の輸送能力評価口頭発表（一般）
- 混相流シンポジウム2020, 2020年08月, 日本語鉛直 U 字管内気液二相流の流動特性に関する研究口頭発表（一般）
- 混相流シンポジウム2020, 2020年08月, 日本語鉛直円管内アルカリ溶液中二酸化炭素気泡の化学吸収に関する研究口頭発表（一般）
- 日本機械学会関西支部第95期定時総会講演会, 2020年03月, 日本語格子ボルツマン法を用いた脳動脈瘤内血流の数値シミュレーション口頭発表（一般）
- 日本流体力学会年会2019, 2019年09月, 日本語, 日本流体力学会, 調布, 国内会議微細流路内テイラー流に及ぼす界面活性剤の影響口頭発表（一般）
- 混相流シンポジウム2019, 2019年08月, 日本語, 日本混相流学会, 福岡, 国内会議鉛直円管内気液対向流における界面及び壁面摩擦係数口頭発表（一般）
- 混相流シンポジウム2019, 2019年08月, 日本語, 日本混相流学会, 福岡, 国内会議旋回環状流における界面及び壁面摩擦係数口頭発表（一般）
- Annual European Rheology Conference 2019 (AERC 2019), 2019年04月, 英語, Portorož, Slovenia, 国際会議Three-Dimensional Flow Characteristics around a Bulge Structure in a Cavity Swept by a Viscoelastic Fluids口頭発表（一般）
- 第20回気液固分散工学サロン, 2019年03月, 日本語, 芝浦工業大学, 国内会議管路内気泡・液滴の運動に及ぼす界面活性剤の影響[招待有り]口頭発表（招待・特別）
- 第95期日本機械学会流体工学部門講演会, 2018年11月, 日本語, 日本機械学会, 室蘭, 国内会議単一液滴に作用する力の実験的分析ポスター発表
- 日本流体力学会年会2018, 2018年09月, 日本語, 大阪大学, 国内会議可溶性界面活性剤により汚染された球形液滴の数値シミュレーション口頭発表（一般）
- 混相流シンポジウム2018, 2018年08月, 日本語, 仙台, 国内会議上端シャープエッジ鉛直管における気液二相対向流の界面及び壁面摩擦係数口頭発表（一般）
- 混相流シンポジウム2018, 2018年08月, 日本語, 仙台, 国内会議鉛直円管内単一二酸化炭素気泡の物質移動に及ぼすアルコールと電解質の相互作用の影響口頭発表（一般）
- 混相流シンポジウム2018, 2018年08月, 日本語, 仙台, 国内会議U字管内気液二相流の圧力損失に及ぼす液相粘度の影響口頭発表（一般）
- 混相流シンポジウム2017, 2017年08月, 日本語, 東京, 国内会議気泡による微細粒子輸送に関する研究口頭発表（一般）
- 混相流シンポジウム2017, 2017年08月, 日本語, 東京, 国内会議鉛直円管内旋回気液二相流の圧力損失に関する研究口頭発表（一般）
- 混相流シンポジウム2017, 2017年08月, 日本語, 東京, 国内会議U字管内気液二相流の圧力損失に関する研究口頭発表（一般）
- International Conference on Numerical Methods in Multiphase Flow, 2017年06月, 英語, Tokyo, 国際会議Numerical Simulation of Contaminated Bubbles and Drops Using Level Set Method (Keynote Lecture)口頭発表（基調）
- 日本機械学会研究分科会RC270第3回分科会, 2016年12月, 日本語, 日本機械学会, 東京, 国内会議二相流数値計算検証のための実験計測シンポジウム・ワークショップパネル（指名）
- 第22回流動化・粒子プロセッシングシンポジウム, 2016年12月, 日本語, 国内会議気泡と液液のシミュレーションとモデリング（招待講演）口頭発表（一般）
- 第94期日本機械学会流体工学部門講演会, 2016年11月, 日本語, 国内会議時空間フィルタ流速計による界面活性剤の吸着を伴う液滴の測定口頭発表（一般）
- 混相流シンポジウム2016, 2016年08月, 日本語, 国内会議塔径及び初期液位が気泡塔内平均ボイド率に及ぼす影響口頭発表（一般）
- 混相流シンポジウム2016, 2016年08月, 日本語, 国内会議静止液中単一気泡の形状に及ぼす界面活性剤の影響に関する研究口頭発表（一般）
- 混相流シンポジウム2016, 2016年08月, 日本語, 国内会議静止液中単一気泡の形状におよぼす界面活性剤の影響に関する研究口頭発表（一般）
- 混相流シンポジウム2016, 2016年08月, 日本語, 国内会議鉛直円管内電解質水溶液中単一二酸化炭素気泡の物質移動に関する研究口頭発表（一般）
- 混相流シンポジウム2016, 2016年08月, 日本語, 国内会議U字管内気液二相流の圧力損失に関する研究口頭発表（一般）
- 日本機械学会関西支部第91期定時総会講演会, pp. 129-132, 2016年03月, 日本語, 大阪, 国内会議散気管内外の流れに散気孔の周方向角度が及ぼす影響口頭発表（一般）
- 第29回数値流体力学シンポジウム, 2015年12月, 日本語, 九州大学筑紫キャンパス, 国内会議埋め込み境界－熱流動格子ボルツマン法における熱流束境界条件の精度に及ぼす二緩和時間衝突則の影響口頭発表（一般）
- American Physical Society 68th Annual DFD Meeting, 2015年11月, 英語, Boston,Massachusetts, U.S.A., 国際会議Reduction of the Temperature Jump in the Immersed Boundary-Thermal Lattice Boltzmann Method口頭発表（一般）
- 日本機械学会研究分科会RC263第6回分科会, 2015年09月, 日本語, 日本機械学会, 東京, 国内会議境界適合型時空間フィルタ流速計と界面活性剤濃度評価の試みシンポジウム・ワークショップパネル（指名）
- 日本流体力学会年会2015, 2015年09月, 日本語, 日本流体力学会, 東京, 国内会議球形液滴界面への界面活性剤の吸着に関する研究口頭発表（一般）
- 混相流シンポジウム2015, 2015年08月, 日本語, 高知工科大学, 国内会議二緩和時間衝突則に基づく埋め込み境界法を適用した熱流動格子ボルツマン法口頭発表（一般）
- 混相流シンポジウム2015, 2015年08月, 日本語, 高知, 国内会議汚染系単一球形液滴に作用するマランゴニ応力の実験的評価口頭発表（一般）
- 混相流シンポジウム2015, 2015年08月, 日本語, 高知, 国内会議一様せん断流中単一気泡に働く揚力に関する研究口頭発表（一般）
- 混相流シンポジウム2015, 2015年08月, 日本語, 高知, 国内会議サブクール水中における蒸気凝縮音に関する研究口頭発表（一般）
- 第20回動力・エネルギー技術シンポジウム, 2015年06月, 日本語, 日本機械学会動力エネルギーシステム部門, 仙台, 国内会議気水分離器における液相分離構造が旋回二相流動に及ぼす影響口頭発表（一般）
- 日本原子力学会 春の大会 2015, 2015年03月, 日本語, 茨城, 国内会議冷却材喪失事故時の格納容器内温度シミュレーションモデルの検証；（2）空気及び蒸気吹込み実験口頭発表（一般）
- 日本機械学会流体工学部門講演会, 2014年10月, 日本語, 富山, 国内会議埋め込み境界を適用した格子ボルツマン法における流速の滑りの除去口頭発表（一般）
- 日本流体力学会年会2014, 2014年09月, 日本語, 日本流体力学会, 宮城県仙台市, 国内会議時空間フィルタ流速計による単一液滴に作用する界面せん断応力の評価口頭発表（一般）
- 理論応用力学講演会, 2014年09月, 日本語, 日本学術会議, 東京都, 国内会議時空間フィルタ流速計による界面近傍速度測定口頭発表（一般）
- 混相流シンポジウム2014, E151, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議非一様電場におけるオリフィスからの気泡成長口頭発表（一般）
- 混相流シンポジウム2014, A311, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議固体壁に付着した液滴に及ぼす界面活性剤の影響に関する研究口頭発表（一般）
- 混相流シンポジウム2014, B331, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議汚れ系における気泡の界面追跡計算と物質移動口頭発表（基調）
- 混相流シンポジウム2014, A311, USB Memory (2 pages), 2014年07月, 日本語, 国内会議スラリー気泡塔内平均ボイド率に及ぼす初期液位の影響口頭発表（一般）
- 混相流シンポジウム2014, B232, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議スラリー気泡塔における気泡合体に粒子径が及ぼす影響口頭発表（一般）
- 混相流シンポジウム2014, B223, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議U字ベンドにおける空気水系気液二相流の圧力損失に関する研究口頭発表（一般）
- 10th International Congress on Membranes and Membrane Processes (ICOM2014), 2014年07月, 英語, Suzhou, China, 国際会議Numerical Study of Concentrated Emulsion Flow Through a Membrane Pore in Oil/Water Separation口頭発表（一般）
- 日本機械学会関西支部定時総会・講演会2014, 2014年03月, 日本語, 大阪府立大学, 国内会議散気管内外の流れに及ぼす液相粘度の影響口頭発表（一般）
- 日本機械学会関西支部定時総会・講演会2014, 2014年03月, 日本語, 大阪府立大学, 国内会議一様せん断流中単一気泡に働く揚力に関する研究 (第2 報 数値計算)口頭発表（一般）
- 日本機械学会 関西支部第89期定時総会講演会, 2014年03月, 日本語, 日本機械学会関西支部, 大阪, 国内会議一様せん断流中単一気泡に働く揚力に関する研究(第 1 報 実験)口頭発表（一般）
- 膜シンポジウム2013, 2013年11月, 日本語, 日本膜学会, 京都府立医科大学 （京都）, 国内会議膜分離プロセスにおける流体挙動のモデリングとシミュレーション口頭発表（一般）
- 混相流シンポジウム2013, 2013年08月, 日本語, 日本混相流学会, 長野, 国内会議散気管内外の流れに液相粘性が及ぼす影響に関する研究口頭発表（一般）
- 混相流シンポジウム2013, 2013年08月, 日本語, 日本混相流学会, 長野, 国内会議鉛直円管内CO2気泡周りの濃度分布計測口頭発表（一般）
- 第62回理論応用力学講演会, 2013年03月, 日本語, 東京工業大学, 国内会議粘性力及び表面張力の評価法が界面運動の数値予測に及ぼす影響口頭発表（一般）
- 日本混相流学会年会講演会2012, pp. 298-299, 2012年08月, 日本語, 日本混相流学会, 柏, 国内会議一様せん断流中気泡に作用する揚力に及ぼす界面活性剤の影響口頭発表（一般）
- 日本混相流学会年会講演会2012, pp. 168-169, 2012年08月, 日本語, 日本混相流学会, 柏, 国内会議Two-Relaxation-Time によるIB-LBM の速度のすべりの除去口頭発表（一般）
- 日本混相流学会年会講演会2012, pp. 166-167, 2012年08月, 日本語, 日本混相流学会, 柏, 国内会議Immersed Boundary-Finite Difference Lattice Boltzmann Method Using Two Relaxation Times口頭発表（一般）
- 第17回動力・エネルギー技術シンポジウム, 2012年06月, 日本語, 日本機械学会動力・エネルギー部門, 福岡, 国内会議4x4ロッドバンドル内気液二相気泡流に関する研究口頭発表（一般）
- 日本混相流学会年会講演会2011, pp. 222-223, 2011年08月, 日本語, 京都, 国内会議膜分離活性汚泥法排水処理装置における散気管内の流れの基礎研究口頭発表（一般）
- 日本混相流学会年会講演会2011, pp. 216-217, 2011年08月, 日本語, 京都, 国内会議高シュミット数・低レイノルズ数気泡の溶解過程口頭発表（一般）
- 日本混相流学会年会2011, pp. 354-355, 2011年08月, 日本語, 日本混相流学会, 京都, 国内会議空間フィルタ流速計による気泡速度の測定口頭発表（一般）
- 日本混相流学会年会講演会2011, pp. 44-45, 2011年08月, 日本語, 日本混相流学会, 京都, 国内会議可溶性界面活性剤の吸脱着を伴う気泡の界面追跡計算口頭発表（一般）
- 日本混相流学会年会講演会2011, pp. 198-199, 2011年08月, 日本語, 京都, 国内会議スラリー系気泡塔における気泡合体に及ぼす粒子濃度及び粒子径の影響口頭発表（一般）
- 日本混相流学会年会講演会2011, pp. 146-147, 2011年08月, 日本語, 京都, 国内会議PWR加圧器サージ管の上端における気液対向流の落下水制限口頭発表（一般）
- 日本混相流学会年会2011, pp. 144-145, 2011年08月, 日本語, 日本混相流学会, 京都, 国内会議4x4ロッドバンドル内気液二相気泡流に関する研究口頭発表（一般）
- 日本機械学会関西支部 第86期定時総会講演会, 2011年03月, 日本語, 日本機械学会関西支部, 京都, 国内会議界面追跡法に基づく気液間物質輸送の数値計算口頭発表（一般）
- 第14回オーガナイズド混相流フォーラム2010, 2010年11月, 日本語, The Japan Society for Multiphase Flow, 札幌，北海道, 国内会議熱・物質輸送を考慮した界面追跡法の現状と課題[招待有り]口頭発表（招待・特別）
- 日本混相流年会2010, pp. 266-267, 2010年07月, 日本語, 浜松, 国内会議多流体モデルによるスラリー中気泡流の数値予測口頭発表（一般）
- 日本混相流学会年会講演会2010, 2010年07月, 日本語, 静岡, 国内会議汚泥減圧浮上濃縮装置の開発口頭発表（一般）
- 日本混相流学会年会講演会2010, pp. 226-227, 2010年07月, 日本語, 浜松, 国内会議鉛直円管内テイラー液滴の終端上昇速度相関式口頭発表（一般）
- 第18回微粒化シンポジウム, 2009年12月, 日本語, 日本液体微粒化学会, 九州大学, 国内会議ノズル内キャビテーション乱流のLES口頭発表（一般）
- 日本混相流学会年会講演会2009, 2009年08月, 日本語, 日本混相流学会, 熊本, 国内会議部分長燃料棒を有するロッドバンドル内の圧力損失口頭発表（一般）
- 日本混相流学会年会講演会2009, pp. 224-225, 2009年08月, 日本語, 国内会議水平管内マイクロバブル乱流の摩擦圧力損失特性に関する研究口頭発表（一般）
- 日本混相流学会年会講演会2009, 2009年08月, 日本語, 日本混相流学会, 熊本, 国内会議気泡による温度成層場の熱輸送促進に関する研究口頭発表（一般）
- 日本混相流学会年会講演会2009, pp. 152-153, 2009年08月, 日本語, 熊本, 国内会議気相体積率に基づく旋回気液二相流の流動様式判別口頭発表（一般）
- 日本混相流学会年会講演会2009, 2009年08月, 日本語, 日本混相流学会, 熊本, 国内会議鉛直円管内単一大液滴の終端上昇速度に関する研究口頭発表（一般）
- 日本高専学会第15回年会・講演会, 2009年08月, 日本語, 愛知, 国内会議サイホンを利用した汚泥減圧濃縮装置の開発口頭発表（一般）
- 日本混相流学会年会講演会2008, pp. 244-245, 2008年06月, 日本語, 会津, 国内会議T字形のエルボを流動するキャビテーションに関する研究口頭発表（一般）
- 日本機械学会流体工学部門講演会, CD-ROM (4 pages), 2007年11月, 日本語, 日本機械学会, 広島, 国内会議T字形エルボにおけるキャビテーションの流動口頭発表（一般）
- 日本混相流学会年会講演会2007, pp. 152-153, 2007年06月, 日本語, 日本混相流学会, 北海道, 国内会議鉛直管内旋回気液二相流の流動様式と旋回長さに関する研究口頭発表（一般）
- 日本流体力学会第20回数値流体力学シンポジウム, CO-ROM (6 pages), 2006年12月, 日本語, 日本流体力学会, 名古屋, 国内会議フェーズフィールド界面追跡法による二相流数値シミュレーション口頭発表（一般）
- 日本混相流学会年会講演会2006, 2006年08月, 日本語, 金沢, 国内会議一様せん断流中単一液滴に作用する揚力に関する研究口頭発表（一般）
- 日本計算工学会第11回計算工学講演会, pp. 415-416, 2006年06月, 日本語, 日本計算工学会, 大阪, 国内会議二相流シミュレーションの界面移流・再構成計算に対するカーン‐ヒリヤード方程式の適用性の検討口頭発表（一般）
- 日本機械学会第18回計算力学講演会, pp. 805-806, 2005年, 日本語, 日本機械学会, 筑波, 国内会議表面張力評価へのカーン‐ヒリヤード方程式の応用口頭発表（一般）
- 日本混相流学会年会講演会, 2004年, 日本語, 日本混相流学会, 岡山大学, 国内会議改良SCA法による液体粒子の界面追跡計算口頭発表（一般）
- 日本混相流学会年会講演会, 2003年, 日本語, 日本混相流学会, 大阪大学, 国内会議微粒化ノズル内キャビテーション流れと液体噴流の微粒化口頭発表（一般）
- 日本機械学会第16回計算力学講演会, 2003年, 日本語, 日本機械学会計算力学部門, 神戸大学工学部, 国内会議簡素化界面再構築法に基づく界面追跡法口頭発表（一般）
- 産学官技術フォーラム’03, 2003年, 日本語, 神戸市立工業高等専門学校, 神戸商工会議所, 神戸市産業振興センター, 国内会議界面追跡法による気液二相流数値予測の試みポスター発表
- 日本混相流学会年会講演会, 2003年, 日本語, 日本混相流学会, 大阪大学, 国内会議液滴・噴霧燃焼の数値シミュレーション口頭発表（一般）
■ 共同研究・競争的資金等の研究課題
- 日本学術振興会, 科学研究費助成事業, 基盤研究(B), 神戸大学, 2024年04月 - 2027年03月汚染系変形気泡に作用する揚力の表面渦度理論に基づく統一的理解と気泡流への展開
- 日本学術振興会, 科学研究費助成事業, 基盤研究(C), 神戸大学, 2023年04月 - 2026年03月, 研究分担者密な鉛直物体群内外における気泡流構造の解明と数値予測技術の開発
- 兵庫県, 成長産業育成のための研究開発支援事業, 可能性調査・研究, 森合精機・神戸大学, 2024年09月 - 2025年03月, その他洗浄機械における旋回流を利用した 油水分離促進技術の可能性調査
- 関西エネルギー・リサイクル科学研究振興財団, 研究助成, 神戸大学, 2024年04月 - 2025年03月, 研究代表者狭隘流路内を運動する気泡の実用的抗力モデルの開発
- JKA, 研究補助, 個別研究, 2023年04月 - 2024年03月, 研究代表者粒子に覆われ安定化された気泡の掩蔽状態よ揚力機構の解明
- 日本学術振興会, 科学研究費助成事業 基盤研究(C), 基盤研究(C), 神戸大学, 2020年04月01日 - 2023年03月31日微細流路内テイラー流を利用した界面活性剤分離器の基盤技術の確立本研究課題では，界面活性剤を含む微細流路内テイラー流の特性を解明するとともに，テイラー流を利用した高効率界面活性剤分離技術の基盤を確立することを目的として各要素研究を推進している．微細流路内テイラー気泡の形成に関して，これまでに界面活性剤の効果により気泡長さが増加する形成モードと減少する形成モードがあることを明らかにしてきた．2021年度は，気液体積流束の範囲を広くとって実験を行い，幅500マイクロメートルの微細流路内での形成モードの流量配分依存性について詳しく調べたところ，本条件の範囲においては，ほとんどの条件で増加のモードであることを明らかにした．その理由について，吸脱着および流れのタイムスケールからなる八田数と，マンランゴニ応力と粘性力のオーダーからなるマランゴニ数をもとにモード遷移の主要因を検討した．界面活性剤分離器開発について，種々の界面活性剤濃度で分離実験を行った．その結果，テイラー流として微細流路を通過させることによって，液相の初期界面活性剤濃度が臨界ミセル濃度前後の広い範囲で半分程度に低減できることを明らかにした．一方，本実験結果から，微細流路内の吸着状態を検討し，微細流路と分離槽における分離過程について検討した．また，微細流路通過後の分離槽における気泡挙動が分離器設計上重要であるため，気泡力学的検討を前年度に引き続いてさらに推進し，揚力や界面活性剤が及ぼす気泡運動への影響への理解を深化させた．
- カシオ科学振興財団, 研究協賛, 2021年12月 - 2022年11月, 研究代表者粒子付着気泡生成器開発とその介在物除去技術への展開
- 日本学術振興会, 科学研究費助成事業, 基盤研究(B), 関西大学, 2019年04月01日 - 2022年03月31日変形・変動を伴う流動界面への界面活性剤吸脱着現象の探求とモデル化多くの工業分野、自然界において気泡や液滴を含む流れが見られ、界面活性剤などの汚染物質の界面への吸着により気泡/液滴と周囲流体間の運動量/熱/物質移動が清浄系と異なることが知られている。静止した界面への界面活性剤の吸着脱離特性に比べ、流動を伴う界面への吸着脱離に関する知見は少なく、数値予測精度向上などを阻んでいる。本研究では、高精度速度分布計測に基づく界面上界面活性剤濃度評価手法を発展させ、液滴表面に作用する局所粘性応力、圧力を評価し、界面活性剤の吸着と流体力の関係を明らかにした。また、当該手法を一般座標系測定領域に拡張することにより、変形液滴界面の流体力学的特性と界面活性剤吸着状態を評価した。
- マツダ財団, マツダ研究助成, 2020年11月 - 2022年03月, 研究代表者強制振動管内気液二相流に及ぼす慣性力の影響の解明
- 科学研究費補助金／基盤研究(A), 2018年04月 - 2022年03月競争的資金
- (財)関西エネルギー・リサイクル科学研究振興財団, 研究助成：エネルギー・リサイクル分野, 2019年04月 - 2020年03月, 研究代表者海底ガス資源回収用垂直管に生じる振動管内気液二相流動に関する基礎研究
- 学術研究助成基金助成金／基盤研究(C), 2017年04月 - 2020年03月競争的資金
- 学術研究助成基金助成金／基盤研究(C), 2017年04月 - 2020年03月, 研究代表者競争的資金
- 日本学術振興会, 科学研究費助成事業, 基盤研究(C), 神戸大学, 2016年04月01日 - 2019年03月31日流れを伴う界面への界面活性剤の吸脱着現象の実験的解明とモデル化多くの工業分野、自然界において気泡や液滴を含む流れが見られ、界面活性剤などの汚染物質を含む場合、汚染物質の界面への吸着により気泡/液滴と周囲流体間の運動量/熱/物質移動が清浄系と異なることが知られている。静止した界面への界面活性剤の吸着脱離特性に比べ、流動を伴う界面への吸着脱離に関する知見は極めて少なく、数値予測精度向上などを阻んでいる。本研究では、高精度速度分布計測に基づく界面上界面活性剤濃度評価手法を構築するとともに、流動を伴う球形液滴界面への界面活性剤の吸着脱離特性を明らかにした。また、数値予測にしばしば使用されている静的な吸着脱離モデルが流動界面には利用できない場合があることを示した。
- 科学研究費一部基金／基盤研究(B), 2015年04月 - 2018年03月競争的資金
- 伊藤忠兵衛基金, 学術研究事業助成, 2015年04月 - 2016年03月, 研究代表者気泡を用いたろ過膜洗浄の計測及び数値シミュレーョンに基づく洗浄速度モデリングに関する研究
- 学術研究助成基金助成金／若手研究(B), 2013年04月 - 2015年03月, 研究代表者競争的資金
- (財)関西エネルギー・リサイクル科学研究振興財団, 研究助成：エネルギー・リサイクル分野, 2012年04月 - 2013年03月, 研究代表者天然ガス液体燃料製造工程において触媒粒子が合成ガス気泡合一に及ぼす影響の解明
- 科学研究費補助金／若手研究(B), 2011年04月 - 2013年03月, 研究代表者競争的資金
- (財)ひょうご科学技術協会, 学術研究助成事業, 2011年04月 - 2012年03月, 研究代表者水質汚染環境下気泡の運動に対する汚染物質の影響のモデル化
- (財)油空圧機器技術振興財団, 研究助成事業, 2009年04月 - 2010年03月, 研究代表者配管内キャビテーション損傷対策のための高安全特殊エルボの開発
研究シーズ
■ 研究シーズ- 複数の相が織りなす複雑な流れ（混相流）の研究シーズカテゴリ：環境・農学, エネルギー, ものづくり技術（機械・電気電子・化学工業）研究キーワード：マイクロ混相流, 微細流路, 気泡, 液滴, 界面活性剤研究内容：マイクロンサイズの微細な流路（マイクロチャネル）を流れるマイクロ気泡や粒子の運動特性、工業機器の各種配管系に現れる気液二相流の流動特性、気泡・液滴に作用する抗力・揚力・壁力、界面活性剤の吸着に起因する界面力や粒子との相互作用など、複数の相（気相・液相・固相）が織りなす複雑な流れを理解し、工学応用へ貢献することを目指して、実験と数値シミュレーションを活用した研究を行っています。