研究者紹介システム

林 公祐
ハヤシ コウスケ
大学院工学研究科 機械工学専攻
准教授
機械工学関係
Last Updated :2022/06/23

研究者情報

所属

  • 【主配置】

    大学院工学研究科 機械工学専攻
  • 【配置】

    工学部 機械工学科, 未来医工学研究開発センター

学位

  • 博士(工学), 神戸大学
  • 複数の相が織りなす複雑な流れ(混相流)の研究

授業科目

ジャンル

  • 科学・技術 / 工学

コメントテーマ

  • 流体工学
  • 数値流体力学
  • 二酸化炭素溶解過程

研究活動

研究キーワード

  • 数値流体力学
  • 混相流
  • 流体工学

研究分野

  • ものづくり技術(機械・電気電子・化学工学) / 流体工学

受賞

  • 2021年10月 神戸大学工学部, 工学部優秀教育賞

    林 公祐

  • 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 Bundle

    HOSOKAWA SHIGEO, HAYASHI KOSUKE, TOMIYAMA AKIO

    The “best” poster based on the structure of the presentation (visual quality), イタリア共和国

    国際学会・会議・シンポジウム等の賞

  • 2013年03月 日本学術会議IUTAM分科会, 第62回理論応用力学講演会 優秀講演賞, 粘性力及び表面張力の評価法が界面運動の数値予測に及ぼす影響

    林 公祐

    36歳未満の若手研究者による優秀な講演に対して授与される

    国内学会・会議・シンポジウム等の賞

  • 2010年07月 日本混相流学会, 日本混相流学会 奨励賞, 円管内気泡及び液滴の終端上昇速度に関する研究

    林 公祐

    「奨励賞」は、本会会誌または他の学術誌、あるいは本会主催・共催の国際会議論文集に掲載された論文または技術報告(原則として3年以内)を対象とし、際立った独創性、新規性、将来性などが期待され40歳以下の本会会員に贈呈する。

    国内外の国際的学術賞

  • 2004年09月 日本混相流学会, 日本混相流学会年会講演会2004 学生優秀講演賞, 改良SCAによる流体粒子の界面追跡計算

    林 公祐

    国内学会・会議・シンポジウム等の賞

  • 2004年02月 日本機械学会, 第18回計算力学講演会 学生優秀講演賞, 表面張力評価へのカーン‐ヒリヤード方程式の応用

    林 公祐

    国内学会・会議・シンポジウム等の賞

論文

  • Devy Setiorini Sa'adiyah, Yuki Matsuo, Michael Schlüter, Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    Elsevier BV, 2021年12月, Chemical Engineering Science, 245, 116852 - 116852, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Void Fraction and Interfacial Friction in Vertical Circular Pipes with the Square Top End under Flooding Conditions

    Toshiya Takaki, Michio Murase, Kosuke Hayashi, Akio Tomiyama

    2021年10月, Nuclear Technology, 英語

    [査読有り]

  • Susumu Osaki, Kosuke Hayashi, Hidehito Kimura, Takeshi Seta, Takashi Sasayama, Akio Tomiyama

    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月25日, Fluids, 6 (10), 338 - 338, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Hendrik Hessenkemper, Dirk Lucas, Dominique Legendre, Akio Tomiyama

    Elsevier BV, 2021年09月, International Journal of Multiphase Flow, 142, 103653 - 103653, 英語

    [査読有り]

    研究論文(学術雑誌)

  • 高木 俊弥, 合田 頼人, 栗本 遼, 林 公祐, 村瀬 道雄, 冨山 明男

    The Japanese Society for Multiphase Flow, 2021年08月02日, 混相流, 35 (3), 463 - 472, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Hidehito Kimura, Susumu Osaki, Kosuke Hayashi, Masaaki Taniguchi, Yuichi Fujita, Takeshi Seta, Akio Tomiyama, Takashi Sasayama, Eiji Kohmura

    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月15日, 混相流, 35 (1), 85 - 92, 日本語

    [査読有り][招待有り]

    研究論文(学術雑誌)

  • Raito Goda, Ryo Kurimoto, Kosuke Hayashi, Michio Murase, Akio Tomiyama

    Elsevier BV, 2021年03月, Nuclear Engineering and Design, 373, 111020 - 111020, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Toshiya Takaki, Raito Goda, Kosuke Hayashi, Michio Murase, Akio Tomiyama

    Elsevier BV, 2021年01月, Nuclear Engineering and Design, 371, 110951 - 110951, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Devy Setiorini Sa'adiyah, Kosuke Hayashi, Ryo Kurimoto, Akio Tomiyama

    Wiley, 2021年01月, Chemie Ingenieur Technik, 93 (1-2), 247 - 259, 英語

    [査読有り][招待有り]

    研究論文(学術雑誌)

  • Yuegui Li, Shigeo Hosokawa, Kosuke Hayashi, Akio Tomiyama, Naoki Shibata, Yasunari Maeda

    Begell House, 2021年, Multiphase Science and Technology, 33 (3), 17 - 31, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Masatoshi Kawabata, Kosuke Hayashi, Ryo Kurimoto, Akio Tomiyama

    Begell House, 2022年01月, Multiphase Science and Technology, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    Begell House, 2021年, Multiphase Science and Technology, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Masatoshi Kawabata, Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    Begell House, 2021年, Multiphase Science and Technology, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Nobuya Kobayashi, Shun Adomi, Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    Begell House, 2021年, Multiphase Science and Technology, 33 (2), 87 - 101, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Dirk Lucas, D. Legendre, Akio Tomiyama

    Begell House, 2021年, Multiphase Science and Technology, 33 (2), 69 - 85, 英語

    [査読有り]

    研究論文(学術雑誌)

  • 村瀬 道雄, 佐野 直樹, 高木 俊弥, 合田 頼人, 林 公祐, 冨山 明男

    The Japanese Society for Multiphase Flow, 2020年11月16日, 混相流, 35 (1), 159 - 167, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Junichiro Kazi, Naoyuki Yoshida, Akio Tomiyama

    Elsevier BV, 2020年10月, International Journal of Multiphase Flow, 131, 103403 - 103403, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Dominique Legendre, Akio Tomiyama

    Elsevier BV, 2020年08月, International Journal of Multiphase Flow, 129, 103350 - 103350, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Numerical Simulations of Flows in Cerebral Aneurysms Using the Lattice Boltzmann Method

    Susumu Osaki, Kosuke Hayashi, Hidehito Kimura, Eiji Kohmura, Akio Tomiyama

    2020年07月, ASME 2020 Fluids Engineering Division Summer Meeting, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Masaaki Hashida, Kosuke Hayashi, Akio Tomiyama

    Elsevier BV, 2020年06月, International Journal of Multiphase Flow, 127, 103270 - 103270, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Douglas Martins Rocha, Fabio Toshio Kanizawa, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama, Gherhardt Ribatski

    Springer Science and Business Media LLC, 2020年04月24日, Flow, Turbulence and Combustion, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Douglas Martins Rocha, Fabio Toshio Kanizawa, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama, Gherhardt Ribatski

    Springer Science and Business Media LLC, 2020年04月, Journal of the Brazilian Society of Mechanical Sciences and Engineering, 42 (4), 英語

    [査読有り]

    研究論文(学術雑誌)

  • Yohei Hori, Dieter Bothe, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    Elsevier BV, 2020年03月, International Journal of Multiphase Flow, 124, 103207 - 103207, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    Begell House, 2020年, Multiphase Science and Technology, 32 (3), 221 - 236, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, J. Kazi, N. Yoshida, Akio Tomiyama

    Begell House, 2020年, Multiphase Science and Technology, 32 (4), 325 - 340, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Application of Spatiotemporal Filter Velocimetry to Evaluation of Drag Force Acting on Single Drop

    Gaku Shigekane, Shigeo Hosokawa, Kosuke Hayashi, Akio Tomiyama

    2019年11月, 11th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF)

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • GODA RAITO, HAYASHI KOSUKE, MURASE MICHIO, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2019年11月, Nuclear Engineering and Design, 353, 110223, 英語

    [査読有り]

    研究論文(学術雑誌)

  • CHEN JUNMING, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2019年09月, Multiphase Science and Technology, 31 (3), 215 - 234, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Evaluation of Drag Force Acting on Single Contaminated Drop based on Velocity Distribution

    Gaku Shigekane, Shigeo Hosokawa, Kosuke Hayashi, Akio Tomiyama

    2019年07月, AJK-Fluids, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Yohei Hori, Yutaka Hirota, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    Elsevier BV, 2019年06月, International Journal of Heat and Mass Transfer, 136, 521 - 530, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Viscous Stress and Pressure Acting on Interface of Contaminated Drop

    HOSOKAWA SHIGEO, SHIGEKANE GAKU, HAYASHI KOSUKE, TOMIYAMA AKIO

    2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • OSAKI SUSUMU, HAYASHI KOSUKE, KIMURA HIDEHITO, SETA TAKESHI, KOHMURA EIJI, TOMIYAMA AKIO

    2019年05月, Computers & Mathematics with Applications, 78, 2746 - 2760, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Effects of Surfactant on Bubble Breakup at T-Junction and Characteristics of Taylor Flows in Mini and Micro Channels

    HAYASHI KOSUKE, SUZUKI DAIKI, KURIMOTO RYO, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Effects of Initial Liquid Height on Gas Holdup in an Air-Water Bubble Column

    KOBAYASHI NOBUYA, HASHIDA MASAAKI, HAYASHI KOSUKE, MARASHDEH QUSSAI, TOMIYAMA AKIO

    2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Combined Effects of Surfactant and Electrolyte on Mass Transfer from Single Carbon-Dioxide Bubbles in Vertical Pipes

    HORI YOHEI, BOTHE DIETER, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    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

    ROCHA DOUGLAS, KANIZAWA FABIO, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO, RIBATSKI GHERHARDT

    2019年05月, 10th International Conference on Multiphase Flow (ICMF2019), 9 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2019年03月, Multiphase Science and Technology, 31 (1), 61 - 71, 英語

    [査読有り]

    研究論文(学術雑誌)

  • HASHIDA MASAAKI, HAYASHI KOSUKE, TOMIYAMA AKIO

    2019年02月, Internation Journal of Multiphase Flow, 111, 285 - 293, 英語

    [査読有り]

    研究論文(学術雑誌)

  • KAWABATA MASATOSHI, TSUKAGUCHI YUICHI, KANG JYUNGYU, HAYASHI KOSUKE, TOMIYAMA AKIO

    2019年02月, ISIJ International, 59 (2), 209 - 215, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Hidehito Kimura, Kosuke Hayashi, Masaaki Taniguchi, Kohkichi Hosoda, Atsushi Fujita, Takeshi Seta, Akio Tomiyama, Eiji Kohmura

    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, 英語, 国際誌

    [査読有り]

    研究論文(学術雑誌)

  • Hidehito Kimura, Masaaki Taniguchi, Kosuke Hayashi, Yosuke Fujimoto, Youichi Fujita, Takashi Sasayama, Akio Tomiyama, Eiji Kohmura

    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, 英語, 国際誌

    [査読有り]

    研究論文(学術雑誌)

  • KURIMOTO RYO, HAYASHI KOSUKE, MINAGAWA HISATO, TOMIYAMA AKIO

    2018年12月, International Journal of Heat and Fluid Flow, 74, 28 - 35, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Shohei Aoyama, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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月01日, Experimental Thermal and Fluid Science, 96, 460 - 469, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Evaluation of Adsorption of Surfactant at a Moving Interface of a Single Spherical Drop

    HOSOKAWA SHIGEO, HAYASHI KOSUKE, TOMIYAMA AKIO

    Elsevier, 2018年09月, Experimental Thermal and Fluid Science, 96, 397 - 405, 英語

    [査読有り]

    研究論文(学術雑誌)

  • AL ISSA SULEIMAN, MURASE MICHIO, TOMIYAMA AKIO, HAYASHI KOSUKE, MACIAN-JUAN RAFAEL

    2018年08月, Nuclear Science and Engineering, 193, 147 - 159, 英語

    [査読有り][招待有り]

    研究論文(学術雑誌)

  • Diagnostics of flow about a contaminated single drop using spatiotemporal filter velocimetry

    HOSOKAWA SHIGEO, SHIGEKANE GAKU, HAYASHI KOSUKE, TOMIYAMA AKIO

    2018年07月, 19th Symposium on Applications of Laser and Imaging Techniques to Fluid Mechanics, Lisbon 2016, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • H. Funahashi, K. Vierow Kirkland, K. Hayashi, S. Hosokawa, A. Tomiyama

    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月15日, Nuclear Engineering and Design, 330, 97 - 105, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Wall and Interfacial Shear Stress of Swirling Annular Flow

    FUNAHASHI HAYATO, 林 公祐, 細川 茂雄, 冨山 明男

    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

    HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2018年04月, 8th European-Japanese Twp-Phase Flow Group Meeting, EJTPFGM-8, 1, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Takashi Seta, Kosuke Hayashi, Akio Tomiyama

    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月10日, International Journal for Numerical Methods in Fluids, 86 (7), 454 - 490, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Raito Goda, Kosuke Hayashi, Karen Vierow Kirkland, Michio Murase, Akio Tomiyama

    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月01日, Nuclear Engineering and Design, 328, 182 - 187, 英語

    [査読有り]

    研究論文(学術雑誌)

  • HAYASHI KOSUKE, TOMIYAMA AKIO

    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), 59 - 66, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Shigeo Hosokawa, Yuya Masukura, Kosuke Hayashi, Akio Tomiyama

    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月01日, International Journal of Multiphase Flow, 97, 157 - 167, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Yohei Hori, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Evaluation of Adsorption of Surfactant at a Drop Interface using Spatiotemporal Filter Velocimetry

    HOSOKAWA SHIGEO, OKAMOTO RYOTA, HAYASHI KOSUKE, TOMIYAMA AKIO

    2017年12月, 10th International Symposium on Measurement Techniques for Multiphase Flow (ISMTMF2017), 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Ryo Sato, Tatsuya Miyayoshi, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Spatiotemporal Flow Structure of Counter-Current Two-Phase Annular Flows

    GODA RAITO, KATSUYA MORI, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2017年11月, ANS Winter Meeting and Nuclear Technology Expo, 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • S. Aoyama, K. Hayashi, S. Hosokawa, D. Lucas, A. Tomiyama

    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), 315 - 320, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • H. Funahashi, K. Hayashi, S. Hosokawa, A. Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • On CCFL in Collider Test Facility: CCFL Results at PWR Hot-Leg Pipe Geometry with 190 mm Diameter

    AL ISSA SULEIMAN, MURASE MICHIO, TOMIYAMA AKIO, HAYASHI KOSUKE, MACIAN-JUAN RAFAEL

    2017年09月, 17th International Topical Meeting on Nuclear Reactor Thermal Hydraulics(NURETH-17), 14 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • 楠木 貴世志, 山本 泰功, 村瀬 道雄, 林 公祐, 細川 茂雄, 冨山 明男

    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 (1), 37 - 46, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • 鉛直管での気液対向流制限に対する流体物性の影響

    村瀬 道雄, 楠木 貴世志, 山本 泰功, 合田 頼人, 林 公祐, 細川 茂雄, 冨山 明男

    2017年06月, 混相流, 31 (2), 152 - 161, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Terminal Velocities of Single Bubble Rising through Stagnant Liquids and Slurry in a Hele-Shaw Cell

    HASHIDA MASAAKI, HAYASHI KOSUKE, TOMIYAMA AKIO

    2017年06月, 9th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamic (ExHFT-9), 9 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Pressure Drop in Air-Water Two-Phase Flows in Horizontal U-Bends

    HAYASHI KOSUKE, SUZUKI DAIKI, MIYACHI NAOTO, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2017年06月, 9th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamic (ExHFT-9), 8 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Evaluation of Adsorption of Surfactant at a Moving Interface of a Single Spherical Drop

    HOSOKAWA SHIGEO, OKAMOTO RYOTA, HAYASHI KOSUKE, TOMIYAMA AKIO

    2017年06月, 9th World Conference on Experimental Heat Transfer, Fluid Mechanics and Thermodynamic (ExHFT-9), 10 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Jiro Aoki, Yohei Hori, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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), 1991 - 2001, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Effects of electrolyte on mass transfer from single carbon-dioxide bubbles in vertical pipes

    HORI YOHEI, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2017年05月, Third International Symposium on Multiscale Multiphase Process Engineering (MMPE2017), 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Experimental Investigation of Temperature Distribution in Condensing Steam-Air Flow in a Circular Tube

    MORI KATSUYA, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO, UTANOHARA YOICHI, YAMAMOTO YASUNORI, MURASE MICHIO

    Experiments 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 m

    2017年04月, 2017 International Congress on Advances in Nuclear Power Plants (ICAPP 2017), 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Shohei Sasaki, Kengo Uchida, Kosuke Hayashi, Akio Tomiyama

    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), 37 - 46, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • S. Aoyama, K. Hayashi, S. Hosokawa, D. Lucas, A. Tomiyama

    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), 392 - 401, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • H. Funahashi, K. Hayashi, S. Hosokawa, A. Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • A Semi-Empirical Correlation for Counter-Current Flow Limitation in a Vertical Pipe

    MORI KATSUYA, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2016年11月, ANS Conference 2016 Winter Meeting and Nuclear Technology Expo, 3 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Shohei Sasaki, Kengo Uchida, Kosuke Hayashi, Akio Tomiyama

    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), 824 - 830, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Numerical Simulation of Surfactant Concentration at an Interface of Single Spherical Drop

    HAYASHI KOSUKE, MASUKURA YUYA, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2016年05月, 9th International Conference on Multiphase Flow, No. 262, 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Mass Transfer from Single Carbon Dioxide Bubbles in Water-Alchol Solutions in a Verical Pipe

    HORI YOHEI, AOKI JIRO, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2016年05月, 9th International Conference on Multiphase Flow, No. 245, 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Experimental Evaluation of Surfacant Concentration at an Interface of Single Spherical Drop

    HOSOKAWA SHIGEO, MASUKURA YUYA, HAYASHI KOSUKE, TOMIYAMA AKIO

    2016年05月, 9th International Conference on Multiphase Flow, No. 262, 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Effects of Azimuthal Angle of Aeration Holes on Flows inside and outside a Bubble Diffuser Pipe

    SATO RYO, MIYAYOSHI TATSUYA, HAYASHI KOSUKE, TOMIYAMA AKIO

    2016年05月, 9th International Conference on Multiphase Flow, No. 168, 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Takayoshi Kusunoki, Takahiro Nozue, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama, Michio Murase

    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), 486 - 495, 英語

    [査読有り]

    研究論文(学術雑誌)

  • SASAKI SHOHEI, HAYASHI KOSUKE, TOMIYAMA AKIO

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • S. Aoyama, K. Hayashi, S. Hosokawa, A. Tomiyama

    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), 343 - 350, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Two-Phase Swirling Flow in a Gas-Liquid Separator

    FUNAHASHI HAYATO, 林 公祐, 細川 茂雄, 冨山 明男

    2015年11月, ANS Conference 2015 Winter Meeting and Nuclear Technology Expo, 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Jiro Aoki, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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), 1955 - 1964, 英語

    [査読有り][招待有り]

    研究論文(学術雑誌)

  • Evaluation of Interfacial Shear Stress of Single Drop Using Spatiotemporal Filter Velocimetry

    HOSOKAWA Shigeo, MASUKURA Yuya, HAYASHI Kosuke, TOMIYAMA Akio

    2015年10月, 7th European-Japanese Two-Phase Flow Group Meeting, 3 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Ryo Sato, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Effects of Liquid Height on Gas Holdup in a Rectangular Bubble Column

    SASAKI SHOHEI, 林 公祐, 冨山 明男

    2015年09月, 9th International Symposium on Measurement Techniques for Multiphase Flow, 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Motion of Contaminated Drops Falling onto Flat Plate

    林 公祐, KITANO MASAAKI, DHAR ABHINAV, 冨山 明男

    2015年07月, ASME-JSME-KSME Joint Fluids Engineering Conference 2015 (AJK2015-FED), AJK2015 - 21048, 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Takayoshi Kusunoki, Michio Murase, Yuki Fujii, Takahiro Nozue, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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), 887 - 896, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Implicit Temperature-Correction-based Immersed-Boundary Thermal Lattice-Boltzmann Method with Two Relaxation Times

    SETA TAKESHI, HAYASHI KOSUKE, TOMIYAMA AKIO

    2015年05月, 23th International Conference on Nuclear Engineering (ICONE23), 10 pages, ICONE23 - 1039, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Jiro Aoki, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Abhinav Dhar, Naoki Shimada, Kosuke Hayashi, Akio Tomiyama

    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月01日, Journal of Computational Multiphase Flows, 7 (1), 15 - 32, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Shimpei Ojima, Shohei Sasaki, Kosuke Hayashi, Akio Tomiyama

    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), 181 - 189, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Masaya Tateyama, Kazutaka Ikeuchi, Shigeo Hosokawa, Akio Tomiyama, Makoto Hirotsu, Nobuhiro Takeda

    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), 67 - 75, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Shimpei Ojima, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • CCFL of Steam-Condensate Two-Phase Flow in a U-Tube

    NOZUE TAKAHIRO, FUJII YUKI, MURASE MICHIO, 林 公祐, 細川 茂雄, 冨山 明男

    2014年11月, ANS Conference 2014 Winter Meeting and Nuclear Technology Expo, CD-ROM, 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Effects of Surfactant on Mass Transfer from Single Carbon Dioxide Bubbles in Vertical Pipes

    AOKI JIRO, HAYASHI KOSUKE, HOSOKAWA SHIGEO, TOMIYAMA AKIO

    2014年09月, Proc. 2nd International Symposium on Multiscale Multiphase Process Engineering, pp. 94 - 99, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Effects of Height-To-Diameter Ratio on Gas Holdup in a Slurry Bubble Column

    SASAKI SHOHEI, OJIMA SHIMPEI, HAYASHI KOSUKE, 冨山 明男

    2014年09月, Proc. 2nd International Symposium on Multiscale Multiphase Process Engineering, pp. 230 - 235, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • An Interface Capturing Scheme for Free-Surface Flows in a Flow Channel Consisting of Solid Obstacles

    SHIMADA NAOKI, DHAR ABHINAV, HAYASHI KOSUKE, 冨山 明男

    2014年09月, Proc. 2nd International Symposium on Multiscale Multiphase Process Engineering, pp. 106 - 115, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Bubble Growth in Microgravity under the Action of Electric Forces: Experiments and Numerical Simulation

    DI MARCO PAOLO, KURIMOTO RYO, SACCONE GIACOMO, HAYASHI KOSUKE, TOMIYAMA AKIO

    2014年08月, Proc. 15th International Heat Transfer Conference, IHTC-15, IHTC15 - 08960, 11 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Numerical Simulations of Clean and Contaminated Bubbles in Linear Shear Flows

    林 公祐, AOYAMA SHOHEI, 冨山 明男

    2014年07月, International Conference on Numerical Methods for Multiphase Flow 2014, ICNMMF2014, 2 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Evaluation of Interfacial Shear Stress of Single Drop using Spatiotemporal Filter Velocimetry

    HOSOKAWA SHIGEO, KITAHATA KYOSUKE, HAYASHI KOSUKE, TOMIYAMA AKIO

    2014年07月, 17th International Symposium on Applications of Laser Techniques to Fluid Mechamics (LISBON 2014), 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • K. Hayashi, A. Tomiyama

    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月01日, Journal of Computational Multiphase Flows, 6 (2), 111 - 126, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Void distribution and bubble motion in bubbly flows in a 4×4 rod bundle. Part II: numerical simulation

    林 公祐, 細川 茂雄, 冨山 明男

    2014年02月, Journal of Nuclear Science and Technology, 51 (5), 580 - 589, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Shigeo Hosokawa, Kosuke Hayashi, Akio Tomiyama

    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), 220 - 230, 英語

    [査読有り]

    研究論文(学術雑誌)

  • SETA, Takeshi, ROJAS, Roberto, 林 公祐, 冨山 明男

    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), 023307, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Shogo Hosoda, Gretar Tryggvason, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Shimpei Ojima, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Abhinav Dhar, Naoki Shimada, Kosuke Hayashi, Akio Tomiyama

    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), 230 - 240, 英語

    [査読有り]

    研究論文(学術雑誌)

  • MASS TRANSFER FROM SINGLE CARBON DIOXIDE BUBBLES IN CONTAMINATED WATER

    Jiro Aoki, Kosuke Hayashi, Shogo Hosoda, Shigeo Hosokawa, Akio Tomiyama

    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 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, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Yusuke Kagawa, Toru Ishigami, Kosuke Hayashi, Hiromi Fuse, Yasushi Mino, Hideto Matsuyama

    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), 7985 - 7992, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Yusuke Kagawa, Toru Ishigami, Kosuke Hayashi, Hiromi Fuse, Yasushi Mino, Hideto Matsuyama

    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), 7985 - 7992, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Mohd Hatta bin Mohd Akbar, Kosuke Hayashi, Dirk Lucas, Akio Tomiyama

    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-Finite Difference Lattice Boltzmann Method Using Two Relaxation Times

    ROJAS, Roberto, SETA, Takeshi, 林 公祐, 冨山 明男

    JSME, 2013年09月, Journal of Fluid Science and Technology, 8 (3), 262 - 276, 英語

    [査読有り]

    研究論文(学術雑誌)

  • 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

    OJIMA Shinpei, HOSOKAWA Shigeo, HAYASHI Kosuke, TOMIYAMA Akio

    2013年09月, 9th UK-Japan Seminar on Multi-Phase Flow, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Paolo Di Marco, Ryo Kurimoto, Giacomo Saccone, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Effects of Inlet Condition on Flow Structure of Bubbly Flow in Rectangular Column

    Mohd Hatta, bin Moh, Akbar, 林 公祐, 冨山 明男, LUCAS, Dirk

    2013年06月, 8th World Conference on Experimental Heat Transfer, Fluid Mechanics, and Thermodynamics (ExHFT8), USB memory, 8 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Void Distribution and Bubble Motion in Bubbly Flows in a 4 X 4 Rod Bundle

    HOSOKAWA Shigeo, HAYASHI Kosuke, TOMIYAMA Akio

    2013年05月, 15th International Topical Meeting on Nuclear Reactor Thermal - Hydraulics, NURETH-15, NURETH15 - 244, USB memory, 12 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • NUMERICAL SIMULATION OF BUBBLY FLOWS IN A 4 X 4 ROD BUNDLE USING A MULTI-FLUID MODEL

    林 公祐, 細川 茂雄, 冨山 明男

    2013年05月, 15th International Topical Meeting on Nuclear Reactor Thermal - Hydraulics, NURETH-15, NURETH15 - 245, USB memory, 12 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Mass Transfer from Single Carbon Dioxide Bubbles in Glycerol-Water Solution

    青木 二郎, 細田 将吾, 林 公祐, 冨山 明男, 細川 茂雄

    2013年05月, 8th International Conference on Multiphase Flow, ICMF2013, ICMF2013 - 229, USB memory, 9 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Effects of Shape Oscillation on Mass Transfer from a Taylor Bubble

    細田 将吾, 林 公祐, 細川 茂雄, 冨山 明男

    2013年05月, 8th International Conference on Multiphase Flow, ICMF2013, ICMF2013 - 286, USB memory, 8 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Effects of Hydrophilic Particles on Bubble Coalescence in Slurry Bubble Column

    小嶋 真平, 林 公祐, 細川 茂雄, 冨山 明男

    2013年05月, 8th International Conference on Multiphase Flow, ICMF2013, ICMF2013 - 157, USB memory, 13 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Roberto Rojas, Kosuke Hayashi, Takeshi Seta, Akio Tomiyama

    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月01日, Journal of Computational Multiphase Flows, 5 (1), 27 - 41, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Two-Relaxation-Time を用いたIB-LBMによる熱流動解析

    瀬田 剛, 内田 洋助, ROJAS Roberto, 林 公祐, 冨山 明男

    2012年12月, 計算数理工学論文集, Vol. 12, pp. 37-42, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Numerical analysis of thermal-hydraulics by immersed boundary-lattice Boltzmann method using two-relaxation-time

    瀬田 剛, 内田 洋助, ROJAS Roberto, 林 公祐, 冨山 明男

    2012年12月, ICOME2012/JASCOME2012

    研究論文(国際会議プロシーディングス)

  • IMMERSED BOUNDARY-FINITE DIFFERENCE LATTICE BOLTZMANN METHOD USING TWO RELAXATION TIMES

    ROJAS, Roberto, 瀬田 剛, 林 公祐, 冨山 明男

    2012年12月, 9th International Conference on CFD in the Minerals and Process Industries,, CD-ROM, 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • EFFECTS OF NUMERICAL TREATMENT OF VISCOUS AND SURFACE TENSION FORCES ON PREDICTED INTERFACE MOTION

    林 公祐, 冨山 明男

    2012年12月, 9th International Conference on CFD in the Minerals and Process Industries,, CD-ROM, 6 pages

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Interface Tracking Simulation of Bubbles in Clean and Contaminated Systems

    林 公祐, 冨山 明男

    2012年11月, Korea-Japan CFD Workshop 2012,, 8 pages

    研究論文(国際会議プロシーディングス)

  • BUBBLE TRACKING SIMULATION OF BUBBLE-INDUCED PSEUDO TURBULENCE

    Mohd Hatta, bin Moh, Akbar, 林 公祐, 細川 茂雄, 冨山 明男

    2012年09月, 6th Japanese-European Two-Phase Flow Group Meeting, USB memory, 10 pages

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • MATSUBAYASHI, Toshiki, KATONO, Kenichi, HAYASHI, Kosuke, TOMIYAMA, Akio

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Roberto Rojas, Takeshi Seta, Akio Tomiyama

    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月01日, Journal of Computational Multiphase Flows, 4 (2), 193 - 209, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Numerical Simulation of Mass Transfer From a Taylor Bubble with Shape Oscillation

    HOSODA, Shogo, HAYASHI, Kosuke, HOSOKAWA, Shigeo, TOMIYAMA, Akio

    2012年06月, International Conference on Numerical Methods in Multiphase Flows, No. 7-1

    研究論文(国際会議プロシーディングス)

  • EFFECTS OF SURFACTANT ON THE MOTION OF BUBBLES IN LINEAR SHEAR FLOWS

    HAYASHI, Kosuke, KURIMOTO, Ryo, HOSOKAWA, Shigeo, TOMIYAMA, Akio

    2012年06月, 2012 Japan-U.S. seminar on Two-Phase flow Dynamics, USB memory, 10 pages

    研究論文(国際会議プロシーディングス)

  • Study on Flows inside and outside an Air Diffuser for Membrane Bioreactor

    KIRA, Fumihiro, FURUNO, Shinsuke, HAYASHI, Kosuke, SAMPEI, Tomoyuki, TOMIYAMA, Akio

    2012年04月, Journal of Fluid Science and Technology, Vol. 7, No. 1, pp. 78-88

    [査読有り]

    研究論文(学術雑誌)

  • Iztok Zun, Matjaz Perpar, Jurij Gregorc, Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Interface Tracking Simulation of Mass Transfer from a Dissolving Bubble

    HAYASHI Kosuke, TOMIYAMA Akio

    2012年03月, Journal of Computational Multiphase Flows, 3(4), 247-261

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Taiga Doi, Takashi Futatsugi, Michio Murase, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Taiga Doi, Takashi Futatsugi, Michio Murase, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Mohd Hatta Bin Mohd Akbar, Kosuke Hayashi, Shigeo Hosokawa, Akio Tomiyama

    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), 197 - 222, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Immersed Boundary-Finite Difference Lattice Boltzmann Method for Liquid-Solid Two-Phase Flows

    ROJAS Roberto, SETA Takeshi, HAYASHI Kosuke, TOMIYAMA Akio

    2011年12月, Journal of Fluid Science and Technology, 6(6), 1051-1064

    [査読有り]

    研究論文(学術雑誌)

  • Multiscale Simulation of Dispersed Flow in Bubble Columns

    TOMIYAMA Akio, OJIMA Shinpei, HAYASHI Kosuke, HOSOKAWA Shigeo

    2011年10月, Proc. of 1st International Symposium on Multiscale Multiphase Process Engineering (MMPE), CD-ROM (6 pages)

    研究論文(国際会議プロシーディングス)

  • Mass transfer from a dissolving carbon dioxide bubble in glycerol-water solution

    HOSODA, Shogo, HAYASHI, Kosuke, TOMIYAMA, Akio

    2011年09月, International Conference on Advanced Technology in Experimental Mechanics 2011

    研究論文(国際会議プロシーディングス)

  • Assessment of turbulence models for bubbly flow in a bubble column

    SATO, Hiroyasu, Mohd Hatta Bin Mohd, Akbar, HAYASHI, Kosuke, TOMIYAMA Akio

    2011年09月, International Conference on Advanced Technology in Experimental Mechanics 2011

    研究論文(国際会議プロシーディングス)

  • Interface tracking and multi-fluid simulation of bubbly flows in bubble columns

    TOMIYAMA Akio, HAYASHI Kosuke

    2011年06月, 8th International Conference on COMPUTATIONAL FLUID DYNAMICS in the Oil & Gas, Metallurgical and Process Industries

    研究論文(国際会議プロシーディングス)

  • Interface tracking simulation of dissolving bubbles and of contaminated bubbles and drops

    TOMIYAMA Akio, HAYASHI Kosuke

    2011年05月, 49th European Two-Phase Flow Group Meeting

    研究論文(国際会議プロシーディングス)

  • K. Hayashi, R. Kurimoto, A. Tomiyama

    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), 241 - 251, 英語

    [査読有り]

    研究論文(学術雑誌)

  • MASS TRANSFER FROM A BUBBLE IN A VERTICAL PIPE

    Shogo Hosoda, Ryosuke Sakata, Kosuke Hayashi, Akio Tomiyama

    Mass 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 - +, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Ryo Kurimoto, Kosuke Hayashi, Akio Tomiyama

    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), 1641 - 1650, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Kosuke Hayashi, Akio Tomiyama

    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), 1593 - 1601, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Tomio Okawa, Kosuke Hayashi, Akio Tomiyama

    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), 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), 305 - 312, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Interface Tracking simulation of mass transfer from a dissolving bubble in a vertical pipe

    HAYASHI Kosuke, TOMIYAMA Akio

    2010年09月, 8th Japan-UK Seminar on Multiphase Flows

    研究論文(国際会議プロシーディングス)

  • Agitation of thermally stratified water layer by free rising wobbling air bubble

    ZUN Iztok, PERPAR Matjaz, GREGORC Jurij, HAYASHI Kosuke, TOMIYAMA Akio

    2010年06月, 48th European Two-Phase Flow Group Meeting 2010, ETPFGM'10

    研究論文(国際会議プロシーディングス)

  • Terminal velocity of a Taylor drop in a vertical pipe

    HAYASHI Kosuke, Ryo Kurimoto, TOMIYAMA Akio

    2010年05月, Proc. of 7th International Conference on Multiphase Flow, ICMF2010, 13 pages

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Terminal velocity of a single drop in a vertical pipe in clean and fully-contaminated systems

    KURIMOTO Ryo, HAYASHI Kosuke, TOMIYAMA Akio

    2010年05月, Proc. of International Conference on Multiphase Flow, ICMF2010, 11 pages

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • LES-Based Bubble Tracking Simulation of Turbulent Cavitation Flow in Injector Nozzle

    SOU Akira, GOSHIMA Daisuke, KINUGASA Tomoyuki, HAYASHI Kosuke, TOMIYAMA Akio

    2010年05月, Proceedings of 7th International Conference on Multiphase Flow, No. 2-1-2, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Kosuke Hayashi, Ryo Kurimoto, Akio Tomiyama

    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月01日, Journal of Computational Multiphase Flows, 2 (1), 47 - 57, 英語

    [査読有り]

    研究論文(学術雑誌)

  • 太田 基希, 細川 茂雄, 林 公祐, 冨山 明男

    2010年03月, 混相流,, Vol. 23, No. 5, pp. 547-554, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • 非構造格子に基づく流体粒子数値予測手法に関する研究

    長 正和, 林 公祐, 赤対 秀明

    2010年03月, 神戸市立工業高等専門学校研究紀要, No. 48, pp. 45-50, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • 活性汚泥の減圧浮上に関する研究

    奥 友晃, 宮明 大輝, 澤井 正和, 林 公祐, 鈴木 隆起, 赤対 秀明

    2010年03月, 神戸高専研究紀要, 49, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • T字形エルボによるキャビテーション壊食対策に関する研究

    苗木 伸也, 柴田 高宏, 林 公祐, 鈴木 隆起, 赤対 秀明

    2010年03月, 神戸市立工業高等専門学校研究紀要, No. 48, pp. 39-44, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Bubble Coalescence in a Slurry Bubble Column

    OBATA Satoshi, OJIMA Shinpei, HAYASHI Kosuke, TOMIYAMA Akio

    2010年, The 8th Japan-UK Seminar on Multiphase Flow

    研究論文(国際会議プロシーディングス)

  • スラリー気泡塔設計のための進化する数値計算手法の構築

    田中 正博, 林 公祐, 冨山 明男

    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), 17 - 24, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Ryo Kurimoto, Akio Tomiyama

    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), 197 - 210, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • 多流体モデルによる気泡塔内多分散気泡流の数値計算

    田中 正博, 林 公祐, 冨山 明男

    2009年11月, 日本機械学会論文集(B編), 75-759, pp.2182-2189, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • ZHANG Zheng, HOSOKAWA Shigeo, HAYASHI Kosuke, TOMIYAMA Akio

    2009年09月, Journal of Power and Energy Systems, Vol. 3, No. 2, pp. 393 - 404, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Dissolution of Carbon Dioxide Bubbles in Clean and Contaminated Systems

    TOMIYAMA, Akio, ABE, Satoru, HAYASHI, Kosuke, HOSOKAWA, Shigeo

    2009年09月, 5th EuropeanTwo-Phase Flow Group Meeting, CD-ROM, 9 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Dimensional Analysis of Terminal Velocityof a Taylor Bubble in a Vertical Pipe

    HAYASHI Kosuke, KURIMOTO Ryo, TOMIYAMA Akio

    2009年09月, 5th EuropeanTwo-Phase Flow Group Meeting, CD-ROM, 6 pages, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • K. Hayashi, A. Tomiyama

    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), 3019 - 3028, 英語

    [査読有り]

    研究論文(学術雑誌)

  • 水平管内におけるマイクロバブルによる摩擦圧力損失低減に関する研究

    星野 剛, 青田 晃一, 林 公祐, 赤対 秀明

    2009年03月, 神戸市立工業高等専門学校研究紀要, No. 47, pp. 13-18, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • 鉛直管内旋回気液二相流の流動様式とその判別法

    高垣 慎太郎, 本間 直樹, 林 公祐, 赤対 秀明

    2009年03月, 神戸市立工業高等専門学校研究紀要, No. 47, pp. 19-24, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Masahiro Tanaka, Kosuke Hayashi, Akio Tomiyama

    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), 141 - 155, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Single Drops Rising Through Stagnant Liquids in a Vertical Pipe

    HAYASHI Kosuke, YOSHIKAWA Yasuki, SOU Akira, TOMIYAMA Akio

    2008年07月, Proc. 5th International Conference on Transport Phenomena in Multiphase Systems, CD-ROM, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • T字形エルボにおけるキャビテーションの流動特性に関する研究

    松岡 隆史, 宇崎 裕真, 林 公祐, 赤対 秀明, 菊谷 功

    2008年03月, 神戸市立工業高等専門学校研究紀要, No. 46, pp. 13-18, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • A Numericsl Study on the Mechanism of Vapor Bubble Lift-off from a Flat Plate

    OKAWA Tomio, HAYASHI Kosuke, TOMIYAMA Akio

    2008年, Proc. 6th Japan-Korea Symp. Nuclear Thermal Hydraulics and Safety, CD-ROM, pp.1-8, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • HAYASHI Kosuke, SOU Akira, TOMIYAMA Akio

    2007年12月, Journal of Fluid Science and Technology, Vol. 2, No. 2, pp. 490-501, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Akio Tomiyama

    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), 121 - 140, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Numerical simulation of single drops in a vertical pipe using various coordinate systems

    Keishun Hirokawa, Kosuke Hayashi, Akira Sou, Akio Tomiyama, Naoki Takada

    The 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), 2410 - 2417, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • 林 公祐, 高田 尚樹, 冨山 明男

    2006年09月, 混相流, 20巻3号, pp.244-251, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Interface Tracking Simulation of Bubbles and Drops in Complex Geometries

    林 公祐, 冨山 明男

    2006年09月, 4th Japanese-European Two-Phase Flow Group Meeting, CD-ROM (8 pages)

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • A Volume Tracking Method Based on Non-Uniform Subcells and Continuum Surface Force Model Using a Local Level Set Function

    HAYASHI Kosuke, SOU Akira, TOMIYAMA Akio

    2006年07月, Computational Fluid Dynamics Journal, Vol. 15, Iss. 2, pp. 225-232, 英語

    [査読有り]

    研究論文(学術雑誌)

  • A. Tomiyama, K. Sakoda, K. Hayashi, A. Sou, N. Shimada, S. Hosokawa

    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), 73 - 110, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Toward Hybrid Simulation of Bubbly Flow (Invited Lecture)

    TOMIYAMA Akio, SAKODA Kenichi, HAYASHI Kosuke

    2005年09月, Proc. Int. Minisymposium on Challenger and Adavances in Flow SImulation and Modeling, p.276, 英語

    研究論文(国際会議プロシーディングス)

  • A volume tracking method based on Advanced Subgrid Counting Algorithm

    Kosuke Hayashi, Win Myint, Akira Sou, Akio Tomiyama

    A 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, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Modelong and Hybrid Simulation of Bubbly Flow(Keynote Lecture)

    TOMIYAMA Akio, SAKODA Kenichi, HAYASHI Kosuke, SOU Akira, SHIMADA Naoki, HOSOKAWA Shigeo

    2004年, Proc.6th US-Japan Seminar on Two-Phase Flow Dynamics, 1,1-16, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • 宋 明良, 林 公祐, 冨山 明男

    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年, 日本機械学会論文集(B編), 70 (698), 2538 - 2544, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Evaluation of Volume Tracking Algorithms for Gas-Liquid Two-Phase Flows.

    SOU Akira, HAYASHI Kosuke, NAKAJIMA Tsuyoshi

    2003年07月, Proc. 4th ASMER-JSME Joint Fluids Engineering Conf., CD-ROM, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Numerical analysis of methanol droplet combustion in hot air flow

    A Sou, K Hayashi, JY Chen

    Numerical 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, 英語

    [査読有り]

    研究論文(国際会議プロシーディングス)

  • Kosuke Hayashi, Yuya Motoki, Matheus J. A. van der Linden, Niels G. Deen, Shigeo Hosokawa, Akio Tomiyama

    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月25日, Fluids, 7 (2), 55 - 55, 英語

    [査読有り]

    研究論文(学術雑誌)

  • 佐野 直樹, 合田 頼人, 林 公祐, 村瀬 道雄, 冨山 明男

    The Japanese Society for Multiphase Flow, 2021年09月30日, 混相流, 35 (4), 533 - 542, 日本語

    [査読有り]

    研究論文(学術雑誌)

  • Hidehito Kimura, Kosuke Hayashi, Susumu Osaki, Ayaka Shibano, Yuichi Fujita, Hiroaki Nagashima, Akio Tomiyama, Takashi Sasayama

    Elsevier BV, 2022年06月, World Neurosurgery, 162, 42 - 42, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Roberto Rojas, Verónica Sotomayor, Tomohiro Takaki, Kosuke Hayashi, Akio Tomiyama

    Elsevier BV, 2022年05月, Computers & Mathematics with Applications, 114, 180 - 187, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Kosuke Hayashi, Junming Chen, Ryo Kurimoto, Akio Tomiyama

    Begell House, 2022年03月, Multiphase Science and Technology, 英語

    [査読有り]

    研究論文(学術雑誌)

  • Fragmentation of Drops Falling through a Miscible Liquid with and without Drop-Drop Interactions

    Rina Nakajima, Kosuke Hayashi

    2022年07月, International Journal of Multiphase Flow, 英語

    [査読有り]

    研究論文(学術雑誌)

MISC

  • 時空間フィルタ流速計による単一液滴に作用する界面せん断応力の評価

    細川 茂雄, 北畑 恭助, 林 公祐, 冨山 明男

    日本流体力学会, 2014年12月, ながれ, 33 (6), 509 - 512, 日本語

    [査読有り][招待有り]

    速報,短報,研究ノート等(学術雑誌)

  • GS35 埋め込み境界を適用した格子ボルツマン法における流速の滑りの除去

    瀬田 剛, ROJAS Roberto, 林 公祐, 冨山 明男

    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", 日本語

  • 二緩和時間衝突則に基づく埋め込み境界法を適用した格子ボルツマン法の精度評価

    瀬田 剛, Rojas Roberto, 林 公祐

    日本計算工学会, 2014年06月, 計算工学講演会論文集 Proceedings of the Conference on Computational Engineering and Science, 19, 6p, 日本語

  • 807 散気管内外の流れに及ぼす液相粘度の影響(OS-3 流体数値計算と実験計測の展開および産業界への応用(2))

    三平 智之, 林 公祐, 細川 茂雄, 冨山 明男, ゾー ランム

    一般社団法人日本機械学会, 2014年03月18日, 関西支部講演会講演論文集, 2014 (89), "8 - 11", 日本語

  • 823 一様せん断流中単一気泡に働く揚力に関する研究 : 第2報 数値計算(GS-5 気泡)

    浜田 剛士, 松永 晃, 林 公祐, 細川 茂雄, 冨山 明男

    一般社団法人日本機械学会, 2014年03月18日, 関西支部講演会講演論文集, 2014 (89), "8 - 27", 日本語

  • 822 一様せん断流中単一気泡に働く揚力に関する研究 : 第1報 実験(GS-5 気泡)

    松永 晃, 濱田 剛士, 細川 茂雄, 林 公祐, 冨山 明男

    一般社団法人日本機械学会, 2014年03月18日, 関西支部講演会講演論文集, 2014 (89), "8 - 26", 日本語

  • 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), 61 - 62, 日本語

  • 526 界面追跡法に基づく気液間物質輸送の数値計算(GS-5 二相流)

    坂田 了介, 細田 将吾, 林 公祐, 冨山 明男

    一般社団法人日本機械学会, 2011年03月19日, 関西支部講演会講演論文集, 2011 (86), "5 - 26", 日本語

  • ノズル内キャビテーション乱流のLES

    五島 代介, 宋 明良, 林 公祐, 冨山 明男

    2009年12月16日, 微粒化シンポジウム講演論文集 = Symposium (ILASS-Japan) on Atomization, 18, 113 - 117, 日本語

  • Flow Visualization of Segregation Process in a Fluidized-Bed by Neutron Radiography

    FURUI Shuji, UMEKAWA Hisashi, HAYASHI Kousuke, OZAWA Mamoru, TAKENAKA Nobuyuki

    2005年, IEEE Trans.on Nuclear Science, 52(1),295-298, 英語

    [査読有り]

    その他

書籍等出版物

  • 流体力学

    冨山 明男, 梶島 岳夫, 加藤 健司, 宋 明良, 高比良 裕之, 林 公祐, 細川 茂雄

    共著, 3章,5章, 朝倉書店, 2020年04月

  • Encyclopedia of Two-Phase Heat Transfer and Flow III: Macro and Micro Flow Boiling and Numerical Modeling Fundamentals

    林 公祐

    その他, World Scientific Publishing, 2018年05月, 英語

    学術書

講演・口頭発表等

  • 汚染系単一球形液滴の数値計算に関する研究

    元木 裕哉, 林 公祐, 細川 茂雄, 冨山 明男

    日本流体力学会年会2021, 2021年09月21日, 日本語

    口頭発表(一般)

  • 堆積層内における気液二相流に関する研究

    田中 裕太, 栗本 遼, 林 公祐, 村瀬 道雄, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

    口頭発表(一般)

  • 振動を伴う鉛直円管内の気液二相流動特性に関する研究

    加藤 大樹, 芳田 直征, 林 公祐, 栗本 遼, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

  • 鉛直円管内気液二相旋回流の界面および壁面摩擦係数に及ぼす旋回強度の影響

    今泉 敦哉, 栗本 遼, 林 公祐, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

  • 鉛直サーペンタイン管内気液二相流のボイド率及び液膜厚さ

    福間 淳平, 加次 淳一郎, 栗本 遼, 林 公祐, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

  • 垂直U字管内気液二相流のボイド率に関する研究

    安井 翔亮, 林 公祐, 栗本 遼, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

  • 気泡流による中空糸膜の汚れ除去と揺動に関する研究

    高屋 慎, 栗本 遼, 林 公祐, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

  • 気泡塔内平均及び局所ボイド率の初期液位依存性に関する研究

    安富 舜, 小林 暢也, 栗本 遼, 林 公祐, 冨山 明男

    混相流シンポジウム2021, 2021年08月22日, 日本語

  • 鉛直円管内気液二相流に及ぼす管加速度の影響に関する研究

    芳田 直征, 林 公祐, 栗本 遼, 冨山 明男

    混相流シンポジウム2020, 2020年08月23日, 日本語

    口頭発表(一般)

  • 鉛直円管内旋回環状流の界面及び壁面摩擦係数に及ぼす液相粘度の影響

    矢原 喜裕, 今泉 敦哉, 栗本 遼, 林 公祐, 冨山 明男

    混相流シンポジウム2020, 2020年08月21日, 日本語

    口頭発表(一般)

  • 気泡群による液相中粒子の輸送能力評価

    川端 雅俊, 藤原 京平, 林 公祐, 冨山 明男

    混相流シンポジウム2020, 2020年08月21日, 日本語

    口頭発表(一般)

  • 鉛直 U 字管内気液二相流の流動特性に関する研究

    森 仁志, 福間 淳平, 加次 淳一郎, 林 公祐, 冨山 明男

    混相流シンポジウム2020, 2020年08月21日, 日本語

    口頭発表(一般)

  • 鉛直円管内アルカリ溶液中二酸化炭素気泡の化学吸収に関する研究

    松尾 優樹, Devy Setiorini Sa’adiyah, 栗本 遼, 林 公祐, 冨山 明男

    混相流シンポジウム2020, 2020年08月21日, 日本語

    口頭発表(一般)

  • 格子ボルツマン法を用いた脳動脈瘤内血流の数値シミュレーション

    大崎 侑, 林 公祐, 木村 英仁, 甲村 英二, 冨山 明男

    日本機械学会関西支部第95期定時総会講演会, 2020年03月11日, 日本語

    口頭発表(一般)

  • 微細流路内テイラー流に及ぼす界面活性剤の影響

    奥居 柊哉, 林 公祐, 芳田 直征, 細川 茂雄, 冨山 明男

    日本流体力学会年会2019, 2019年09月15日, 日本語, 日本流体力学会, 調布, 国内会議

    口頭発表(一般)

  • 鉛直円管内気液対向流における界面及び壁面摩擦係数

    合田 頼人, 島村 長幸, 林 公祐, 細川 茂雄, 村瀬 道雄, 冨山 明男

    混相流シンポジウム2019, 2019年08月07日, 日本語, 日本混相流学会, 福岡, 国内会議

    口頭発表(一般)

  • 旋回環状流における界面及び壁面摩擦係数

    矢原 喜裕, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2019, 2019年08月05日, 日本語, 日本混相流学会, 福岡, 国内会議

    口頭発表(一般)

  • Three-Dimensional Flow Characteristics around a Bulge Structure in a Cavity Swept by a Viscoelastic Fluids

    SUZUKI Hiroshi, KAWATA Masaki, HIDEMA Ruri, HOSOKAWA Shigeo, HAYASHI Kosuke, TOMIYAMA Akio

    Annual European Rheology Conference 2019 (AERC 2019), 2019年04月, 英語, Portorož, Slovenia, 国際会議

    口頭発表(一般)

  • 管路内気泡・液滴の運動に及ぼす界面活性剤の影響

    林 公祐

    第20回気液固分散工学サロン, 2019年03月, 日本語, 芝浦工業大学, 国内会議

    [招待有り]

    口頭発表(招待・特別)

  • 単一液滴に作用する力の実験的分析

    重兼 楽, 細川 茂雄, 林 公祐, 冨山 明男

    第95期日本機械学会流体工学部門講演会, 2018年11月, 日本語, 日本機械学会, 室蘭, 国内会議

    ポスター発表

  • 可溶性界面活性剤により汚染された球形液滴の数値シミュレーション

    林 公祐, 酒井 航平, 細川 茂雄, 冨山 明男

    日本流体力学会年会2018, 2018年09月, 日本語, 大阪大学, 国内会議

    口頭発表(一般)

  • 上端シャープエッジ鉛直管における気液二相対向流の界面及び壁面摩擦係数

    島村 長幸, 合田 頼人, 林 公祐, 細川 茂雄, 村瀬 道雄, 西田 浩二, 冨山 明男

    混相流シンポジウム2018, 2018年08月, 日本語, 仙台, 国内会議

    口頭発表(一般)

  • 鉛直円管内単一二酸化炭素気泡の物質移動に及ぼすアルコールと電解質の相互作用の影響

    廣田 雄嵩, 堀 陽平, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2018, 2018年08月, 日本語, 仙台, 国内会議

    口頭発表(一般)

  • U字管内気液二相流の圧力損失に及ぼす液相粘度の影響

    林 公祐, 鈴木 大貴, 細川 茂雄, 冨山 明男

    混相流シンポジウム2018, 2018年08月, 日本語, 仙台, 国内会議

    口頭発表(一般)

  • 気泡による微細粒子輸送に関する研究

    姜 俊奎, 林 公祐, 冨山 明男, 田口 謙治, 塚口 友一

    混相流シンポジウム2017, 2017年08月, 日本語, 東京, 国内会議

    口頭発表(一般)

  • 鉛直円管内旋回気液二相流の圧力損失に関する研究

    船橋 駿斗, 木下 博仁, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2017, 2017年08月, 日本語, 東京, 国内会議

    口頭発表(一般)

  • U字管内気液二相流の圧力損失に関する研究

    鈴木 大貴, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2017, 2017年08月, 日本語, 東京, 国内会議

    口頭発表(一般)

  • Numerical Simulation of Contaminated Bubbles and Drops Using Level Set Method (Keynote Lecture)

    HAYASHI KOSUKE

    International Conference on Numerical Methods in Multiphase Flow, 2017年06月, 英語, Tokyo, 国際会議

    口頭発表(基調)

  • 二相流数値計算検証のための実験計測

    細川 茂雄, 林 公祐, 冨山 明男

    日本機械学会研究分科会RC270第3回分科会, 2016年12月, 日本語, 日本機械学会, 東京, 国内会議

    シンポジウム・ワークショップパネル(指名)

  • 気泡と液液のシミュレーションとモデリング(招待講演)

    林 公祐

    第22回流動化・粒子プロセッシングシンポジウム, 2016年12月, 日本語, 国内会議

    口頭発表(一般)

  • 時空間フィルタ流速計による界面活性剤の吸着を伴う液滴の測定

    細川 茂雄, 益倉 侑冶, 林 公祐, 冨山 明男

    第94期日本機械学会流体工学部門講演会, 2016年11月, 日本語, 国内会議

    口頭発表(一般)

  • 塔径及び初期液位が気泡塔内平均ボイド率に及ぼす影響

    内田 賢吾, 橋田 昌明, 佐々木 翔平, 林 公祐, 冨山 明男

    混相流シンポジウム2016, 2016年08月, 日本語, 国内会議

    口頭発表(一般)

  • 静止液中単一気泡の形状に及ぼす界面活性剤の影響に関する研究

    竹川 稔彦, 青山 昌平, 林 公祐, 冨山 明男

    混相流シンポジウム2016, 2016年08月, 日本語, 国内会議

    口頭発表(一般)

  • 静止液中単一気泡の形状におよぼす界面活性剤の影響に関する研究

    竹川 稔彦, 青山 昌平, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2016, 2016年08月, 日本語, 国内会議

    口頭発表(一般)

  • 鉛直円管内電解質水溶液中単一二酸化炭素気泡の物質移動に関する研究

    堀 陽平, 青木 二郎, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2016, 2016年08月, 日本語, 国内会議

    口頭発表(一般)

  • U字管内気液二相流の圧力損失に関する研究

    宮地 祥人, 林 公祐, 細川 茂雄, 冨山 明男, 舘山 将也, 廣津 誠, 竹田 信宏

    混相流シンポジウム2016, 2016年08月, 日本語, 国内会議

    口頭発表(一般)

  • 散気管内外の流れに散気孔の周方向角度が及ぼす影響

    佐藤 稜, 宮吉 達也, 林 公祐, 冨山 明男

    日本機械学会関西支部第91期定時総会講演会, pp. 129-132, 2016年03月, 日本語, 大阪, 国内会議

    口頭発表(一般)

  • 埋め込み境界-熱流動格子ボルツマン法における熱流束境界条件の精度に及ぼす二緩和時間衝突則の影響

    瀬田 剛, 林 公祐, 冨山 明男

    第29回数値流体力学シンポジウム, 2015年12月, 日本語, 九州大学筑紫キャンパス, 国内会議

    口頭発表(一般)

  • Reduction of the Temperature Jump in the Immersed Boundary-Thermal Lattice Boltzmann Method

    SETA TAKESHI, HAYASHI KOSUKE, TOMIYAMA AKIO

    American Physical Society 68th Annual DFD Meeting, 2015年11月, 英語, Boston,Massachusetts, U.S.A., 国際会議

    口頭発表(一般)

  • 境界適合型時空間フィルタ流速計と界面活性剤濃度評価の試み

    細川 茂雄, 林 公祐, 冨山 明男

    日本機械学会研究分科会RC263第6回分科会, 2015年09月, 日本語, 日本機械学会, 東京, 国内会議

    シンポジウム・ワークショップパネル(指名)

  • 球形液滴界面への界面活性剤の吸着に関する研究

    益倉 侑冶, 細川 茂雄, 林 公祐, 冨山 明男

    日本流体力学会年会2015, 2015年09月, 日本語, 日本流体力学会, 東京, 国内会議

    口頭発表(一般)

  • 二緩和時間衝突則に基づく埋め込み境界法を適用した熱流動格子ボルツマン法

    瀬田 剛, 林 公祐, 冨山 明男

    混相流シンポジウム2015, 2015年08月, 日本語, 高知工科大学, 国内会議

    口頭発表(一般)

  • 汚染系単一球形液滴に作用するマランゴニ応力の実験的評価

    益倉 侑冶, 細川 茂雄, 林 公祐, 冨山 明男

    混相流シンポジウム2015, 2015年08月, 日本語, 高知, 国内会議

    口頭発表(一般)

  • 一様せん断流中単一気泡に働く揚力に関する研究

    青山 昌平, 竹川 稔彦, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2015, 2015年08月, 日本語, 高知, 国内会議

    口頭発表(一般)

  • サブクール水中における蒸気凝縮音に関する研究

    池内 和孝, 舘山 将也, 林 公祐, 細川 茂雄, 冨山 明男, 廣津 誠, 竹田 信宏

    混相流シンポジウム2015, 2015年08月, 日本語, 高知, 国内会議

    口頭発表(一般)

  • 気水分離器における液相分離構造が旋回二相流動に及ぼす影響

    上遠野 健一, 船橋 駿斗, 細川 茂雄, 林 公祐, 冨山 明男

    第20回動力・エネルギー技術シンポジウム, 2015年06月, 日本語, 日本機械学会動力エネルギーシステム部門, 仙台, 国内会議

    口頭発表(一般)

  • 冷却材喪失事故時の格納容器内温度シミュレーションモデルの検証;(2)空気及び蒸気吹込み実験

    野末 貴大, 藤井 雄基, 林 公祐, 細川 茂雄, 冨山 明男, 歌野原 陽一, 山本 泰功, 村瀬 道雄

    日本原子力学会 春の大会 2015, 2015年03月, 日本語, 茨城, 国内会議

    口頭発表(一般)

  • 埋め込み境界を適用した格子ボルツマン法における流速の滑りの除去

    瀬田 剛, ROJAS ROBERT, 林 公祐, 冨山 明男

    日本機械学会流体工学部門講演会, 2014年10月, 日本語, 富山, 国内会議

    口頭発表(一般)

  • 時空間フィルタ流速計による単一液滴に作用する界面せん断応力の評価

    細川 茂雄, 北畑 恭助, 林 公祐, 冨山 明男

    日本流体力学会年会2014, 2014年09月, 日本語, 日本流体力学会, 宮城県仙台市, 国内会議

    口頭発表(一般)

  • 時空間フィルタ流速計による界面近傍速度測定

    細川 茂雄, 北畑 恭助, 林 公祐, 冨山 明男

    理論応用力学講演会, 2014年09月, 日本語, 日本学術会議, 東京都, 国内会議

    口頭発表(一般)

  • 非一様電場におけるオリフィスからの気泡成長

    栗本 遼, DI MARCO PAOLO, 南川 久人, 林 公祐, 冨山 明男

    混相流シンポジウム2014, E151, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議

    口頭発表(一般)

  • 固体壁に付着した液滴に及ぼす界面活性剤の影響に関する研究

    北野 正章, ABHINAV DHAR, 林 公祐, 冨山 明男

    混相流シンポジウム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月, 日本語, 札幌, 国内会議

    口頭発表(一般)

  • U字ベンドにおける空気水系気液二相流の圧力損失に関する研究

    小林 慎吾, 林 公祐, 細川 茂雄, 冨山 明男, 廣津 誠, 舘山 将也, 竹田 信宏

    混相流シンポジウム2014, B223, USB Memory (2 pages), 2014年07月, 日本語, 札幌, 国内会議

    口頭発表(一般)

  • Numerical Study of Concentrated Emulsion Flow Through a Membrane Pore in Oil/Water Separation

    Yusuke KAGAWA, Toru ISHIGAMI, Hiromi FUSE, Yasushi MINO, Kosuke HAYASHI, Hideto MATSUYAMA

    10th International Congress on Membranes and Membrane Processes (ICOM2014), 2014年07月, 英語, Suzhou, China, 国際会議

    口頭発表(一般)

  • 散気管内外の流れに及ぼす液相粘度の影響

    三平 智之, 林 公祐, 細川 茂雄, 冨山 明男, ラン ムゾー

    日本機械学会関西支部定時総会・講演会2014, 2014年03月, 日本語, 大阪府立大学, 国内会議

    口頭発表(一般)

  • 一様せん断流中単一気泡に働く揚力に関する研究 (第2 報 数値計算)

    濵田 剛士, 松永 晃, 林 公祐, 細川 茂雄, 冨山 明男

    日本機械学会関西支部定時総会・講演会2014, 2014年03月, 日本語, 大阪府立大学, 国内会議

    口頭発表(一般)

  • 一様せん断流中単一気泡に働く揚力に関する研究(第 1 報 実験)

    松永 晃, 濱田 剛士, 細川 茂雄, 林 公祐, 冨山 明男

    日本機械学会 関西支部第89期定時総会講演会, 2014年03月, 日本語, 日本機械学会関西支部, 大阪, 国内会議

    口頭発表(一般)

  • 膜分離プロセスにおける流体挙動のモデリングとシミュレーション

    石神 徹, 布施 ひろみ, 香川 裕輔, 浅尾 慎一, 林 公祐, 三野 泰志, 松山 秀人

    膜シンポジウム2013, 2013年11月, 日本語, 日本膜学会, 京都府立医科大学 (京都), 国内会議

    口頭発表(一般)

  • 散気管内外の流れに液相粘性が及ぼす影響に関する研究

    林 公祐, 三平 智之, 冨山 明男, 古野 真介, 吉良 文博

    混相流シンポジウム2013, 2013年08月, 日本語, 日本混相流学会, 長野, 国内会議

    口頭発表(一般)

  • 鉛直円管内CO2気泡周りの濃度分布計測

    細田 将吾, 青木 二郎, 林 公祐, 細川 茂雄, 冨山 明男

    混相流シンポジウム2013, 2013年08月, 日本語, 日本混相流学会, 長野, 国内会議

    口頭発表(一般)

  • 粘性力及び表面張力の評価法が界面運動の数値予測に及ぼす影響

    林 公祐, 冨山 明男

    第62回理論応用力学講演会, 2013年03月, 日本語, 東京工業大学, 国内会議

    口頭発表(一般)

  • 一様せん断流中気泡に作用する揚力に及ぼす界面活性剤の影響

    林 公祐, 松永 晃, 濵田 剛士, 冨山 明男

    日本混相流学会年会講演会2012, pp. 298-299, 2012年08月, 日本語, 日本混相流学会, 柏, 国内会議

    口頭発表(一般)

  • Two-Relaxation-Time によるIB-LBM の速度のすべりの除去

    瀬田 剛, ROJAS, Roberto, 林 公祐, 冨山 明男

    日本混相流学会年会講演会2012, pp. 168-169, 2012年08月, 日本語, 日本混相流学会, 柏, 国内会議

    口頭発表(一般)

  • Immersed Boundary-Finite Difference Lattice Boltzmann Method Using Two Relaxation Times

    ROJAS, Roberto, 瀬田 剛, 林 公祐, 冨山 明男

    日本混相流学会年会講演会2012, pp. 166-167, 2012年08月, 日本語, 日本混相流学会, 柏, 国内会議

    口頭発表(一般)

  • 4x4ロッドバンドル内気液二相気泡流に関する研究

    細川 茂雄, 村松 祐哉, 多田 遼, 林 公祐, 冨山 明男

    第17回動力・エネルギー技術シンポジウム, 2012年06月, 日本語, 日本機械学会動力・エネルギー部門, 福岡, 国内会議

    口頭発表(一般)

  • 膜分離活性汚泥法排水処理装置における散気管内の流れの基礎研究

    吉良 文博, 川岸 朋樹, 古野 真介, 三平 智之, 林 公祐, 冨山 明男

    日本混相流学会年会講演会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月, 日本語, 京都, 国内会議

    口頭発表(一般)

  • PWR加圧器サージ管の上端における気液対向流の落下水制限

    柳 千裕, 二木 高志, 林 公祐, 冨山 明男, 木下 郁男, 村瀬 道雄

    日本混相流学会年会講演会2011, pp. 146-147, 2011年08月, 日本語, 京都, 国内会議

    口頭発表(一般)

  • 4x4ロッドバンドル内気液二相気泡流に関する研究

    村松 祐哉, 細川 茂雄, 林 公祐, 冨山 明男

    日本混相流学会年会2011, pp. 144-145, 2011年08月, 日本語, 日本混相流学会, 京都, 国内会議

    口頭発表(一般)

  • 界面追跡法に基づく気液間物質輸送の数値計算

    坂田 了介, 細田 将吾, 林 公祐, 冨山 明男

    日本機械学会関西支部 第86期定時総会講演会, 2011年03月, 日本語, 日本機械学会関西支部, 京都, 国内会議

    口頭発表(一般)

  • 熱・物質輸送を考慮した界面追跡法の現状と課題

    林 公祐, 冨山 明男, ZUN Iztok

    第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月, 日本語, 浜松, 国内会議

    口頭発表(一般)

  • ノズル内キャビテーション乱流のLES

    五島 代介, 宋 明良, 林 公祐, 冨山 明男

    第18回微粒化シンポジウム, 2009年12月, 日本語, 日本液体微粒化学会, 九州大学, 国内会議

    口頭発表(一般)

  • 部分長燃料棒を有するロッドバンドル内の圧力損失

    太田 基希, 林 公祐, 細川 茂雄, 冨山 明男

    日本混相流学会年会講演会2009, 2009年08月, 日本語, 日本混相流学会, 熊本, 国内会議

    口頭発表(一般)

  • 水平管内マイクロバブル乱流の摩擦圧力損失特性に関する研究

    星野 剛, 青田 晃一, 林 公祐, 赤対 秀明

    日本混相流学会年会講演会2009, pp. 224-225, 2009年08月, 日本語, 国内会議

    口頭発表(一般)

  • 気泡による温度成層場の熱輸送促進に関する研究

    林 公祐, PERPAR Matjaz, ZUN Iztok, 冨山 明男

    日本混相流学会年会講演会2009, 2009年08月, 日本語, 日本混相流学会, 熊本, 国内会議

    口頭発表(一般)

  • 気相体積率に基づく旋回気液二相流の流動様式判別

    山下 恭平, 高垣 慎太郎, 本間 直樹, 林 公祐, 赤対 秀明

    日本混相流学会年会講演会2009, pp. 152-153, 2009年08月, 日本語, 熊本, 国内会議

    口頭発表(一般)

  • 鉛直円管内単一大液滴の終端上昇速度に関する研究

    吉河 泰輝, 栗本 遼, 林 公祐, 冨山 明男

    日本混相流学会年会講演会2009, 2009年08月, 日本語, 日本混相流学会, 熊本, 国内会議

    口頭発表(一般)

  • サイホンを利用した汚泥減圧濃縮装置の開発

    奥 友晃, 島田 光浩, 林 公祐, 澤井 正和, 望月 修一, 赤対 秀明

    日本高専学会第15回年会・講演会, 2009年08月, 日本語, 愛知, 国内会議

    口頭発表(一般)

  • T字形のエルボを流動するキャビテーションに関する研究

    苗木 伸也, 松岡 隆史, 林 公祐, 赤対 秀明, 菊谷 功, 苗木 伸也

    日本混相流学会年会講演会2008, pp. 244-245, 2008年06月, 日本語, 会津, 国内会議

    口頭発表(一般)

  • T字形エルボにおけるキャビテーションの流動

    松岡 隆史, 宇崎 裕真, 林 公祐, 赤対 秀明, 掛札 理美, 日詰 誠, 菊谷 功

    日本機械学会流体工学部門講演会, CD-ROM (4 pages), 2007年11月, 日本語, 日本機械学会, 広島, 国内会議

    口頭発表(一般)

  • 鉛直管内旋回気液二相流の流動様式と旋回長さに関する研究

    高垣 慎太郎, 竹田 博昭, 林 公祐, 赤対 秀明

    日本混相流学会年会講演会2007, pp. 152-153, 2007年06月, 日本語, 日本混相流学会, 北海道, 国内会議

    口頭発表(一般)

  • フェーズフィールド界面追跡法による二相流数値シミュレーション

    高田 尚樹, 廣川 景俊, 林 公祐, 冨山 明男

    日本流体力学会第20回数値流体力学シンポジウム, CO-ROM (6 pages), 2006年12月, 日本語, 日本流体力学会, 名古屋, 国内会議

    口頭発表(一般)

  • 一様せん断流中単一液滴に作用する揚力に関する研究

    小川 耕平, Win Myint, 林 公祐, 細川 茂雄, 冨山 明男

    日本混相流学会年会講演会2006, 2006年08月, 日本語, 金沢, 国内会議

    口頭発表(一般)

  • 二相流シミュレーションの界面移流・再構成計算に対するカーン‐ヒリヤード方程式の適用性の検討

    高田 尚樹, 廣川 景俊, 林 公祐, 冨山 明男

    日本計算工学会第11回計算工学講演会, pp. 415-416, 2006年06月, 日本語, 日本計算工学会, 大阪, 国内会議

    口頭発表(一般)

  • 表面張力評価へのカーン‐ヒリヤード方程式の応用

    林 公祐, 高田 尚樹, 冨山 明男

    日本機械学会第18回計算力学講演会, pp. 805-806, 2005年, 日本語, 日本機械学会, 筑波, 国内会議

    口頭発表(一般)

  • 改良SCA法による液体粒子の界面追跡計算

    林 公祐, Win Myint, 宋 明良, 冨山 明男

    日本混相流学会年会講演会, 2004年, 日本語, 日本混相流学会, 岡山大学, 国内会議

    口頭発表(一般)

  • 微粒化ノズル内キャビテーション流れと液体噴流の微粒化

    宋 明良, 林 公祐, 新田 真一, 冨山 明男

    日本混相流学会年会講演会, 2003年, 日本語, 日本混相流学会, 大阪大学, 国内会議

    口頭発表(一般)

  • 簡素化界面再構築法に基づく界面追跡法

    宋 明良, 林 公祐, 冨山 明男

    日本機械学会第16回計算力学講演会, 2003年, 日本語, 日本機械学会計算力学部門, 神戸大学工学部, 国内会議

    口頭発表(一般)

  • 界面追跡法による気液二相流数値予測の試み

    林 公祐, 宋 明良, 冨山 明男

    産学官技術フォーラム’03, 2003年, 日本語, 神戸市立工業高等専門学校, 神戸商工会議所, 神戸市産業振興センター, 国内会議

    ポスター発表

  • 液滴・噴霧燃焼の数値シミュレーション

    宋 明良, 林 公祐, 國本 雅也, 冨山 明男

    日本混相流学会年会講演会, 2003年, 日本語, 日本混相流学会, 大阪大学, 国内会議

    口頭発表(一般)

所属学協会

  • 日本流体力学会

    2018年05月 - 現在
  • 日本高専学会

    2009年07月 - 2010年06月
  • 日本混相流学会

  • 日本機械学会

共同研究・競争的資金等の研究課題