Authorship：Last author,　Corresponding author   Publishing type：Research paper (scientific journal)   Publisher：Begell House  

Takuma Kogawa;Yousuke Nakamura;Soma Taniguchi;Hisashi Nakajima;Junnosuke Okajima;Atsuki Komiya;Shigenao Maruyama;Atsushi Sakurai
Recently, selective heating techniques by laser therapy have been proposed. Among them, laser therapy adopting metallic nanoparticles has garnered attention. By injecting metallic nanoparticles inside a tissue, high absorption can be achieved by the local surface plasmonic resonance at the surface of the nanoparticle. In particular, gold nanorods exhibit significantly high resonance against laser and their effectiveness has been investigated experimentally. However, the modeling of laser therapy by gold nanorods is limited, and a quantitative numerical simulation is required to further investigate the possibility of effective heating by gold nanorods. To investigate the effectiveness of selective heating by laser therapy with gold nanorods, a three-dimensional calculation of laser therapy is performed in this study. To model the thermal conduction and propagation of laser inside a tissue, the Pennes bio-heat transfer equation is applied. The radiative heat transfer is calculated using the radiation element method. The dispersion effect of gold nanorod injection is considered. Comparing the thermal damage between gold nanorod injection and noninjection cases, the performance of selective heating by gold nanorods is evaluated. Furthermore, the effect of gold nanorod injection location on the thermal damage is discussed. Numerical calculation results indicate that gold nanorod injection enables selective heating against the target as well as noninvasive heating on the tissue surface. Furthermore, the heating depth can be controlled by changing the injection location appropriately.

DOI： 10.1615/heattransres.2021036835

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：American Chemical Society (ACS)  

Atsushi Sakurai;Kyohei Yada;Tetsushi Simomura;Shenghong Ju;Makoto Kashiwagi;Hideyuki Okada;Tadaaki Nagao;Koji Tsuda;Junichiro Shiomi
We computationally designed an ultranarrow-band wavelength-selective thermal radiator via a materials informatics method alternating between Bayesian optimization and thermal electromagnetic field calculation. For a given target infrared wavelength, the optimal structure was efficiently identified from over 8 billion candidates of multilayers consisting of multiple components (Si, Ge, and SiO2). The resulting optimized structure is an aperiodic multilayered metamaterial exhibiting high and sharp emissivity with a Q-factor of 273. The designed metamaterials were then fabricated, and reasonable experimental realization of the optimal performance was achieved with a Q-factor of 188, which is significantly higher than those of structures empirically designed and fabricated in the past. This is the first demonstration of the experimental realization of metamaterials designed by Bayesian optimization. The results facilitate the machine-learning-based design of metamaterials and advance our understanding of the narrow-band thermal emission mechanism of aperiodic multilayered metamaterials.

DOI： 10.1021/acscentsci.8b00802

PubMed

researchmap

A. Sakurai;Y. Matsuno
In this study, a tungsten-SiO2-based metal–insulator–metal-structured metasurface for the thermal emitter of the thermophotovoltaic system was designed and fabricated. The proposed emitter was fabricated by applying the photolithography method. The fabricated emitter has high emissivity in the visible to near-infrared region and shows excellent wavelength selectivity. This spectral emissivity tendency agreed well with the result calculated by the finite-difference time-domain method. Additionally, the underlying mechanism of its emission was scrutinized. Study of the fabrication process and theoretical mechanisms of the emission, clarified in this research, will be fundamental to design the wavelength-selective thermal emitter.

DOI： 10.3390/mi10020157

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Elsevier Ltd  

Takuma Kogawa;Lin Chen;Junnosuke Okajima;Atsushi Sakurai;Atsuki Komiya;Shigenao Maruyama
To use and understand the large-scale natural convection in the building and systems, heat transfer analysis of the natural convection with radiation is important. The radiation effect is influenced by concentration of the participating media in the considered system. In this study, an effects of concentration of participating media on turbulent natural convection inside a cubic cavity were investigated. The turbulence effect was modeled by large eddy simulation. Further, the radiative heat transfer was calculated by the radiation element method using the ray emission model. In addition, the gas radiation was modeled by the full-spectrum k-distribution method. By varying the partial pressures of H2O and CO2, the effect of the participating media were evaluated. From the calculation, gas radiation effects affected the flow instability and restrained the temperature stratification when the concentration of the participating media was higher. However, the convective heat transfer have not been affected by the significant flow instability by the radiation effects.

DOI： 10.1016/j.applthermaleng.2017.11.135

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Taylor and Francis Ltd.  

Atsushi Sakurai;Ryo Kanbayashi;Koji Matsubara
We performed a direct numerical simulation (DNS) of the turbulent natural convection between two parallel plates at a Rayleigh number of Ra = 8.0 × 106, focusing on the turbulent natural convection affected by radiation in an optically thick fluid (τ ≥ 10). When the effects of the radiation were considered, the flow structure and temperature distribution in the channel changed as the optical thickness of the fluid increased. The effects of the radiation on the turbulent natural convection were clearly explained by the turbulence statistics from the DNS results.

DOI： 10.1080/10407782.2017.1412713

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：OPTICAL SOC AMER  

Yuki Matsuno;Atsushi Sakurai
Metasurfaces with metal-insulator-metal structures have attracted significant attention because of their high optical performance. Their spectral emissivity/absorptivity can be modulated with geometrically controlled electromagnetic resonances. However, practical fabrication of these materials in the meter scale remains a challenge. In the present study, an aluminum-and ceria-based metasurface that exhibits nearly perfect emission/absorption in the infrared region was designed. Furthermore, we succeeded in fabricating the proposed metasurface in the meter scale via photolithography and wet etching. This study enhances understanding of underlying resonance mechanisms with electromagnetic simulations and equivalent circuit model, and will facilitate the experimental design of high-efficiency largearea metasurfaces. (C) 2017 Optical Society of America

DOI： 10.1364/OME.7.000618

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：ELSEVIER SCIENCE BV  

Yuki Matsuno;Atsushi Sakurai
The purpose of this study is to clarify the optical characteristics of a Tungsten SiO2-based film-couple metamaterial with a fishnet-shaped grating, and we aim to show that this structure could potentially be used a solar selective absorber for a solar thermophotovoltaic system. The proposed film-coupled metamateris absorber combined with a fishnet-shaped grating shows significant enhancement in its absorption in the spectral region compared with a flat Tungsten surface, and it keeps its spectral emission low in the infrare region, thereby reducing radiative heat loss. The underlying mechanisms of the proposed absorber are discusse through a 3D full-wave electromagnetic simulation, the results of which are compared with that of theoretic equations. Furthermore, its spectral absorption under oblique incident light at the transverse magnetic am transverse electric waves is scrutinized. The underlying absorption mechanisms and relations between eac optical resonance are discussed in this paper and will prove to be fundamental not only in the design of selective absorbers but also in other wavelength-selective thermal radiation controlling devices.

DOI： 10.1016/j.optcom.2016.10.044

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Takuma Kogawa;Junnosuke Okajima;Astushi Sakurai;Atsuki Komiya;Shigenao Maruyama
The turbulence structure, boundary layer, and heat transfer of turbulent natural convection, including radiation effects inside the cubic cavity, were investigated in this study. Large eddy simulation was conducted, and Vreman model was adopted for the dynamic subgrid-scale model. To calculate the radiative heat flux efficiently, a coupled calculation method, using the radiation element method by ray emission model, was constructed. To separate the radiation effects of the gas and surface radiations, four calculation conditions have been analyzed; non-radiation, gas radiation, surface radiation, and combined radiation. Observing the vortices structure using the Q value revealed that the surface radiation effect was more dominant for the flow instability than the gas radiation effect. An evaluation of the boundary layers for both the gas and surface radiation effects showed that the flow circulation inside the cubic cavity was enhanced. The surface radiation was dominant in the generation of the shear stress by the turbulent flow. The total heat transfer, which includes the convective and radiative heat transfers, have been investigated. The surface radiation affected on the radiative heat transfer significantly as compared to the gas radiation. The radiation effects changed the radiative heat transfers, while the convective heat transfers of all the calculation conditions were similar. (C) 2016 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.ijheatmasstransfer.2016.08.059

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Atsushi Sakurai;Tomoaki Kawamata
We numerically investigate a solar-selective absorber with tungsten core-shell nano particle arrays embedded in an SiO2 layer. The 3D full-wave finite-difference time-domain (FDTD) simulations are performed to investigate the geometric effects of different types of solar-selective absorbers. Consequently, broadband light absorption was achieved with either a tungsten nanoparticle array or a tungsten core-shell nanoparticle array because of the strong electric field enhancement in the gap between the core nanoparticles. The solar performance of the proposed structure is shown for high-efficiency solar light absorption. This study enhances understanding of the light absorption mechanism of metallic nanoparticle/dielectric composite and facilitates the design of high-efficiency solar-selective absorbers. (C) 2016 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.jqsrt.2016.08.005

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：AMER INST AERONAUTICS ASTRONAUTICS  

Bo Zhao;Atsushi Sakurai;Zhuomin M. Zhang
Due to the interference effect, the reflectance and transmittance of an arbitrary linearly polarized wave cannot be fully described by the reflectance and transmittance of transverse electric and transverse magnetic waves. In this work, the polarization dependence of reflectance and transmittance for anisotropic metamaterials is examined. The explicit mathematical relationship between the Fresnel coefficients and the reflectance (or transmittance) for a given polarization angle is derived. A method is presented that can be used to obtain the reflectance (or transmittance) for any polarization angle from three other polarizations. Several different anisotropic metamaterials are numerically investigated with the finite-difference time-domain technique to elucidate the polarization dependence of reflectance or transmittance. The reflectance and transmittance obtained from the three-polarization-angle method are in agreement with those from the rigorous numerical simulation.

DOI： 10.2514/1.T4587

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Atsushi Sakurai;Bo Zhao;Zhuomin M. Zhang
This work deals with the spectral radiative properties of metamaterials made of a 2D metallic pattern on a dielectric thin film that separates the top periodic structure from a metal ground plane. Two distinct patterns are considered: one consists of disconnected double-rectangle gold pattern and the other is made of L-shape pattern. Both structures exhibit dual-band emission or absorption peaks in the infrared region. Unlike the disconnected rectangular pattern for which the normal emittance is independent of the polarization angle, the L-shape structure shows strong polarization dependent due to near-field coupling. The electromagnetic fields at the resonant frequencies are analyzed using the finite-difference time-domain technique to elucidate the physical mechanisms accounting for the different behaviors. For both structures, the resonance mechanisms can be explained by magnetic polaritons. Inductor-capacitor circuit models are developed to quantitatively predict the resonance frequencies. This work helps the understanding of polarization dependence in metamaterials and may facilitate their applications in biosensing and infrared spectroscopic system. (C) 2014 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.jqsrt.2014.11.018

Web of Science

researchmap

Language：English   Publishing type：Research paper (international conference proceedings)   Publisher：Elsevier BV  

K. Matsubara;H. Sakai;Y. Kazuma;A. Sakurai;T. Kodama;N. Gokon;H.S. Cho;K. Yoshida
This study describes the flow and thermal field of a two-tower fluidized receiver, aimed to be incorporated into a beam-down reflector system. An experimental visualization and numerical simulation were made to accumulate phenomenological knowledge, in order to complete the design of a demonstration model. Visualization of the cold particle bed with no irradiation revealed that global circulation occurs between the two towers by the aeration with different line velocities. The numerical simulation revealed that this global circulation enhances the transport of sensible heat from the irradiated layer in the high pressure tower to the low pressure tower. The global circulation in the two-tower fluidized receiver has the potential to be extended for use in a thermal receiver and direct storage system. Unexpectedly, local circulations occur on the high and low pressure sides, and are stronger than the global circulation. These local circulations contribute to the thermal mixing in each tower. The bottom distributor on the low pressure side can cause an uprising flow which leads to the local circulation. The design of the distributor is suggested to be elaborated based on these findings. (C) 2015 The Authors. Published by Elsevier Ltd.

DOI： 10.1016/j.egypro.2015.03.057

Web of Science

Scopus

researchmap

Language：English   Publisher：Japan Institute of Energy  

Koji Matsubara;Yuki Kazuma;Atsushi Sakurai;Tatsuya Kodama;Nobuyuki Gokon;Hyun-Seok Cho;Kazuo Yoshida;Yoshinori Nagase
This study proposes a t wo-tower t ype fluidized receiver for receiving concentrated s olar light at high temperatures is proposed. This fluidized receiver is coupled with the Miyazaki beam-down reflector system. The visualization and the numerical simulation was made for the laboratory model to show the fluid-dynamics and the heat transport as well as to elaborate the existing demonstration model. The experimental visualization of the cold model demonstrated that in each tower, the particles are moved by the aerated organized flow. The numerical simulation was made using the thermal input cited from the 5.0kWth sun simulator. The numerical simulation demonstrated that the inner particles are heated above 950 °C after 3.0 minutes of the aeration start time. The increase in temperature is caused by the intrusion of the heated particles from the high-pressure side tower to the low-pressure side tower. It is thus suggested that collective circulation occurs between high- and low-pressure side towers and this can be utilized for a direct thermal storage system.

DOI： 10.3775/jie.94.1323

Scopus

CiNii Article

researchmap

Other Link： https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-25289357/

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Atsushi Sakurai;Bo Zhao;Zhuomin M. Zhang
Metamaterial thermal emitters and absorbers have been widely studied for different geometric patterns by exciting a variety of electromagnetic resonances. A resistor-inductor-capacitor (RLC) circuit model is developed to describe the magnetic resonances (i.e. magnetic polaritons) inside the structures. The RLC circuit model allows the prediction of not only the resonance frequency, but also the full width at half maximum and quality factor for various geometric patterns. The parameters predicted by the RLC model are compared with the finite-difference time-domain simulation. The magnetic field distribution and the power dissipation density profile are also used to justify the RLC circuit model. The geometric effects on the resonance characteristics are elucidated in the wire (or strip), cross, and square patterned metamaterial in the infrared region. This study will facilitate the design of metamaterial absorbers and emitters based on magnetic polaritons. (C) 2014 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.jqsrt.2014.07.024

Web of Science

researchmap

Language：English   Publishing type：Research paper (international conference proceedings)   Publisher：Elsevier BV  

K. Matsubara;Y. Kazuma;A. Sakurai;S. Suzuki;L. Soon-Jae;T. Kodama;N. Gokon;C. Hyun Seok;K. Yoshida
This study proposes a novel fluidized receiver for absorbing concentrated solar light at high temperatures. Previously, a tubular receiver and a volumetric receiver were developed to make high-temperature air for a solar gas turbine system. The aim was to combine these elements with a tower reflector; however, it was challenging to install these heavy receivers on the top of the tower. Currently, a fluidized receiver prototype is tested by a 3 kWh solar simulator in preparation for a field test at the Miyazaki beam-down reflector system. The fluid dynamics of the prototype receiver is numerically investigated. The currently treated receiver is an inner-circulating fluidized bed spouted by concentric gas streams with high and low velocities in the center and outer annulus, respectively. The draft tube is submerged in the particles to organize particle circulation. Concentrated light irradiates the particles through a quartz window at the top of the receiver container. Such a fluidized bed was first adopted by Kodama et al. for thermochemical reactions; however, it is currently pursued for its potential as a high-temperature receiver aimed at concentrated solar power generation. Experiments of the prototype receiver (inner diameter = 45 mm) demonstrated that the inner particles are heated to a temperature greater than 900 degrees C and that an increase of the central gas velocity removes the excess temperature near the particle bed surface. A numerical computation suggests that the large-scale circulation of particles leads to the activation of thermal mixing. The currently proposed receiver is thus expected to attenuate re-radiation losses likely to occur in a conventional volumetric porous receiver. The scale-up of the receiver is being considered by the numerical computation for a field test in the Miyazaki 100 kWh beam-down reflector system. (C) 2013 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license

DOI： 10.1016/j.egypro.2014.03.048

Web of Science

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Atsushi Sakurai;Hiroya Tanikawa;Makoto Yamada
The purpose of this study is to computationally design a wide-angle cermet-based solar selective absorber for high temperature applications by using a characteristic matrix method and a genetic algorithm. The present study investigates a solar selective absorber with tungsten-silica (W-SiO2) cermet. Multilayer structures of 1, 2, 3, and 4 layers and a wide range of metal volume fractions are optimized. The predicted radiative properties show good solar performance, i.e., thermal emittances, especially beyond 2 mu m, are quite low, in contrast, solar absorptance levels are successfully high with wide angular range, so that solar photons are effectively absorbed and infrared radiative heat loss can be decreased. (C) 2013 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.jqsrt.2013.03.004

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Atsushi Sakurai;Koji Matsubara;Kenji Takakuwa;Ryo Kanbayashi
The purpose of this study is to clarify the radiation effects on mixed turbulent convection in a horizontal channel. The present study provides turbulence statistics using direct numerical simulation (DNS) in an optically thin medium. When the radiation effect is considered, the flow structure and the temperature distribution in the channel change with an increase in the optical thickness of the fluid. The radiation effect changes the distributions of the temperature fluctuation intensity and the turbulent heat flux. These radiation effects on mixed convection can be clearly explained by the turbulence statistics obtained from the DNS results. (C) 2012 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.ijheatmasstransfer.2012.01.006

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：ASME-AMER SOC MECHANICAL ENG  

Koji Matsubara;Atsushi Sakurai;Takahiro Miura;Takuya Kawabata
The near-wall turbulent heat transport of the three orthogonal directions was directly solved for the Prandtl numbers ranging from 0.025 to 5.0 to validate the algebraic models of the turbulent heat flux. Two kinds of thermal situations were considered in the low Reynolds number turbulent flow: (a) the case with a uniform heat flux in the spanwise direction (UHF) and (b) the case with the mean spanwise temperature gradient (STG). Among the turbulent heat flux models tested, the model of Rogers preferably predicted over the treated range of the Prandtl numbers, but it failed to reproduce the low Prandtl number effects very accurately. This paper revealed that the coefficient of the Rotta model can be modified to include the low Prandtl number effects by means of the correlation between the exact coefficient suggested by DNS and the Prandtl number. [DOI: 10.1115/1.4005077]

DOI： 10.1115/1.4005077

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：TAYLOR & FRANCIS INC  

Koji Matsubara;Katsuya Nishiwaki;Okiteru Fukuda;Togo Yamaguchi;Atsushi Sakurai
A newly developed electric furnace was used to pyrolyze benzene, methane, or acetylene diluted by a nonreactive gas, and the carbon aggregates were separated from the hydrogen-rich gases. The temperature was changed up to 2200K for various levels of initial concentration of hydrocarbon and residence time, and how the experimental parameters affected the mean diameter of carbon aggregates was examined. The temperature increase in experiments led to reduction of the mean diameter of the carbon aggregates, whereas increase in the initial concentration of hydrocarbon or the residence time resulted in slightly lifting it. The mean diameter removed of temperature effects correlated well with the initial carbon atom mol concentration of hydrocarbon, N-C, and this parameter was used to develop the experimental equations. This article presents three kinds of experimental equations, and the prediction accuracy was confirmed to be improved when accounting for increasing the experimental parameters.

DOI： 10.1080/00102202.2012.691582

Web of Science

researchmap

Language：Japanese   Publishing type：Research paper (scientific journal)  

Takahiro Miura;Koji Matsubara;Atsushi Sakurai
Direct numerical simulation was performed for air turbulent flows and heat transfer in a channel with square ribs attached to one wall. The spatially advancing cases and the fully advanced cases were both treated for the bulk Reynolds number of about 4,560. The rib height, H, was kept constant at 0.2 times the channel half-width, whereas the rib pitch was changed from Pi = 2H to 16 H. The heat transfer performance was evaluated using the ratio between the mean Nusselt number of the ribbed case and that of the smooth case requiring the same pumping power. This ratio was maximized at Pi/H = 9 for the fully advanced cases, but at Pi/H = 2 or 4 for the advancing cases. The mechanisms leading to this difference were discussed through turbulence statistics and their visualization. The thin thermal boundary layer in the advancing region was shown to increase heat transfer on the extended surface of Pi/H = 2 and 4, resulting the higher heat transfer performance than the case of sparse ribs. © 2012 The Japan Society of Mechanical Engineers.

DOI： 10.1299/kikaib.78.607

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：JAPAN SOC MECHANICAL ENGINEERS  

Takahiro Miura;Koji Matsubara;Atsushi Sakurai
Heat transfer and fluid flow in a single-rib mounting channel were investigated by directly solving Navier-Stokes and energy equations. The flow and thermal fields were considered to be fully developed at the inlet of the channel, and the simulation was made for spatial advancement of turbulent heat transfer. The Reynolds number based on the friction velocity at the inlet and the channel half width was 150. The Prandtl number was 0.71. The budgets for turbulent heat fluxes and temperature variance at various sections were presented and were investigated, which would be useful for testing and developing turbulence models. Near a circular vortex in front of the rib, pressure diffusion terms made an important contribution. Remarkable production terms were generated near a reattachment point. Production and dissipation terms were not dominant in front of and above the rib, and a time scale ratio exceeded 2.0 in the region.

DOI： 10.1299/jtst.7.120

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：AMER INST PHYSICS  

Koji Matsubara;Atsushi Sakurai;Kenya Yamazaki;Makoto Takeda
Direct numerical simulation was performed for the spatially advancing turbulent flow in a two-dimensional curved channel. The fully developed turbulence was introduced into the curved channel through a straight part, and observation was made for how the near-wall turbulence advanced along curved walls. The frictional Reynolds number, Re(tau 0), was fixed at 150, whereas the radius ratio of the curved channel, alpha, was changed at three steps; 0.8, 0.92, and 0.975. Five computational tests were performed, which include simulations with changing the spanwise length of the channel or the grid arrangement. The maximum number of grid points was 43,610,112 (=936 x 91 x 512) used in the case of alpha = 0.92. Numerical results for the curved flow were nearly free from the grid arrangement or the spanwise length of the channel, and changes of the mean velocity qualitatively agreed with the experiment by Kobayashi [Trans. Jpn. Soc. Mech. Eng., Ser. B 57, 4064 (1991)], which both supported credibility of the simulation. The spanwise pre-multiplied power spectrum and the spanwise two-point correlation captured the peaks corresponding to the micro-scale structures near the straight wall, which continuously grew along the concave wall and advanced to exhibit the wave length of organized flows. Therefore, the micro-scale structures such as quasi-stream vortices and the ejection near the outer wall of the straight channel were suggested to trigger the onset of initial seeds of the organized structure, which grew into the large scale vortices expanding almost over the channel width.(C) 2011 American Institute of Physics. [doi: 10.1063/1.3584126]

DOI： 10.1063/1.3584126

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：TAYLOR & FRANCIS INC  

Koji Matsubara;Takahiro Miura;Atsushi Sakurai;Kenya Yamazaki;Makoto Takeda
Direct numerical simulation was performed for the spatially advancing turbulent flow and heat transfer in a two-dimensional curved channel equating the radius ratio to 0.92 or 0.8. The frictional Reynolds number was fixed at 150, whereas the Prandtl number was set at 0.71. According to the numerical result, the remarkable enhancement of heat transfer occurred on the outer wall, suggesting the organized vortex activated the heat transfer. The budgets of Reynolds stresses clarified that the onset and growth of the organized flow was assisted by the direct energy transfer from the mean flow.

DOI： 10.1080/10407782.2011.588562

Web of Science

researchmap

Language：Japanese   Publisher：Japan Society of Mechanical Engineers  

Takashi Uchida;Koji Matsubara;Takahiro Miura;Suguru Tachikawa;Atsushi Sakurai
Direct numerical simulation was performed for a spatially advancing turbulent flow and heat transfer in a two-dimensional curved channel. A rib was placed on either an inner-or outer-wall of the straight part nearby curved-channel inlet for controlling the flow and thermal fields. Outer-wall was heated at the constant temperature and inner-wall was cooled similarly. In the simulation, fully developed flow and temperature computed by the straight channel was given to the inlet of curved channel domain. The frictional Reynolds number was assigned at 150, and the Prandtl number was given 0.71. When a rib was attached to the inner-wall, heat transfer in the upstream region of curved channel was enhanced effectively. On the other hand, placing a rib on the outer-wall caused an enhancement of heat transfer in the relatively downstream region of curved part, suggesting active heat transport of the large-scale vortices. Consequently, the inner-rib case affected more greatly on heat transfer than the outer-rib case. Mean velocity vector and power spectrum analysis implied that development of coherent structures near the outer-wall were controlled by a rib. It was found that inner-rib restrained large-scale vortex, while outer-rib promoted growth and power of large-scale vortex. © 2011 The Japan Society of Mechanical Engineers.

DOI： 10.1299/kikaib.77.2170

Scopus

CiNii Article

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：JAPAN SOC MECHANICAL ENGINEERS  

Atsushi Sakurai;Ryo Kanbayashi;Koji Matsubara;Shigenao Maruyama
Thermal radiation in a turbulent natural convection plays an important role in a wide area of engineering and nature. The purposes of this study are to investigate the effects of turbulent fluctuation on radiative heat transfer, and to evaluate radiative heat transfer models applied to turbulent natural convection. The present radiative heat transfer analysis of a turbulent natural convection using direct numerical simulation (DNS) provides a useful fundamental data for the complete coupling simulation in the future.

DOI： 10.1299/jtst.6.449

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：ASME-AMER SOC MECHANICAL ENG  

Atsushi Sakurai;Shigenao Maruyama;Koji Matsubara;Takahiro Miura;Masud Behnia
The purpose of this paper is to consider a possibility of the independent column approximation for solving the radiative heat fluxes in a 3D turbulent channel flow. This simulation method is the simplest extension of the plane-parallel radiative heat transfer. The test case of the temperature profile was obtained from the direct numerical simulation. We demonstrate the comparison between the 3D radiative transfer simulation and the independent column approximation with an inhomogeneous temperature field and optical properties. The above mentioned results show the trivial discrepancies between the 3D simulation and the independent column approximation. The required processing time for the independent column approximation is much faster than the 3D radiative transfer simulation due to the simple algorithm. Although the independent column simulation is restricted to simple configurations such as channel flow in this paper, wide application areas are expected due to the computational efficiency.

DOI： 10.1115/1.4000240

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Koji Matsubara;Akihiko Matsui;Takahiro Miura;Atsushi Sakurai;Hitoshi Suto;Koji Kawai
Direct numerical simulation was performed for a spatially advancing turbulent flow and heat transfer in a two-dimensional curved channel, where one wall was heated to a constant temperature and the other wall was cooled to a different constant temperature. In the simulation, fully developed flow and temperature from the straight-channel driver was passed through the inlet of the curved-channel domain. The frictional Reynolds number was assigned 150, and the Prandtl number was given 0.71. Since the flow field was examined in the previous paper, the thermal features are mainly targeted in this paper. The turbulent heat flux showed trends consistent with a growing process of large-scale vortices. In the curved part, the wall-normal component of the turbulent heat flux was twice as large as the counterpart in the straight part, suggesting active heat transport of large-scale vortices. In the inner side of the same section, temperature fluctuation was abnormally large compared with the modest fluctuation of the wall-normal velocity. This was caused by the combined effect of the large-scale motion of the vortices and the wide variation of the mean temperature
in such a temperature distribution, large-scale ejection of the hot fluid near the outer wall, which is transported into the near inner-wall region, should have a large impact on the thermal boundary layer near the inner wall. Wave number decomposition was conducted for various statistics, which showed that the contribution of the large-scale vortex to the total turbulent heat flux normal to the wall reached roughly 80% inside the channel 135° downstream from the curved-channel inlet. © 2009 Wiley Periodicals, Inc.

DOI： 10.1002/htj.20275

Scopus

researchmap

Language：Japanese   Publishing type：Research paper (scientific journal)   Publisher：Japan Society of Mechanical Engineers  

Atsushi Sakurai;Koji Matsubara
Radiative Properties of high temperature water vapor play an important role for radiative heat transfer analysis in combustion engineering, fire safety science, etc. Absorption spectrum of visible to infrared light inherently depends on the wavelength of electromagnetic waves. Ignoring the spectral dependence makes radiative heat transfer analyses wrong results. To overcome the problem, various nongray gas models have been proposed. However, there are a few studies for summarizing these characteristics in prediction. The goal of this study is to give valuable information for engineering applications. In the present study, radiative properties and simulation characteristics of high-temperature water vapor using the Planck-mean approximation, the statistical narrow band model, and the full-spectrum k-distribution model are assessed rigorously by comparing with the Line-by-Line method. The recent HITRAN/HITEMP databases are employed to calculate Line-by-Line spectrum.

DOI： 10.1299/kikaib.76.771_1681

Scopus

CiNii Article

CiNii Books

researchmap

Other Link： http://id.ndl.go.jp/bib/10916168

Language：English   Publishing type：Research paper (scientific journal)   Publisher：TAYLOR & FRANCIS INC  

Atsushi Sakurai;Subhash C. Mishra;Shigenao Maruyama
This article deals with the implementation of the radiation element method (REM) with the lattice Boltzmann method (LBM) to solve a combined mode transient conduction-radiation problem. Radiative information computed using the REM is provided to the LBM solver. The planar conducting-radiating participating medium is contained between diffuse gray boundaries, and the system may contain a volumetric heat generation source. Temperature and heat flux distributions in the medium are studied for different values of parameters such as the extinction coefficient, the scattering albedo, the conduction-radiation parameter, the emissivity of the boundaries, and the heat generation rate. To check the accuracy of the results, the problem is also solved using the finite-volume method (FVM) in conjunction with the LBM. In this case, the data for radiation field are calculated using the FVM. The REM has been found to be compatible with the LBM, and in all the cases, results of the LBM-REM and the LBM-FVM have been found to provide an excellent comparison.

DOI： 10.1080/10407780903583008

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：JAPAN SOC MECHANICAL ENGINEERS  

Takahiro Miura;Koji Matsubara;Atsushi Sakurai
Heat transfer and fluid flow in a single-rib mounting channel were investigated by directly solving Navier-Stokes and energy equations. Flow and thermal fields were considered to be fully developed at the inlet of the channel, and the simulation was made for spatial advancement of turbulent heat transfer. Keeping the frictional Reynolds number, Re-tau 0, at 150, the rib height ratio was changed in four steps from H/delta = 0.05 to H/delta = 0.4. Computational results were confirmed to be nearly independent of grid meshes. In addition, numerical accuracy was confirmed through close agreement between computed mean pressure and the experiment by Yao et al. (1995). The numerical results revealed that the highest value of the mean Nusslet number was as large as 1.3 times the smooth surface consuming the same pumping power, and the local enhancement of heat transfer was correlated with the turbulence increase near the rib front and the reattachment point. According to the Reynolds stress budgets for H/delta = 0.2, there were mechanisms to induce powerful fluctuations: (1) Streamwise fluctuation was increased through production by flow deceleration in the upstream of the rib; (2) Redistribution to wall-normal and spanwise fluctuations was fortified by the fluid splattering to the rib front. Therefore, excellent performance of heat transfer was concluded to occur due to flow structures, which induce the strong disturbance near the rib front triggering smooth transition of the separated shear layer.

DOI： 10.1299/jtst.5.135

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：TAYLOR & FRANCIS INC  

Atsushi Sakurai;Shigenao Maruyama;Koji Matsubara
The purpose of this study is to develop the radiation element method by ray emission method (REM2) code appropriate for coupling with the bioheat transfer equation, and to clarify the photothermal effect of various parameters inside biological tissues. First, the REM2, which involves the air-tissue interface effect, is validated with the existing literature. In order to clarify the effects of optical and thermophysical properties of biological tissues, a one-dimensional tissue model of CW light transport and bioheat transfer is employed. The present study provides nondimensional results, which are obtained by varying refractive index, extinction coefficient, scattering albedo, blood perfusion parameter, and conduction-radiation parameter, show valuable guidance for understanding the coupled light and bioheat transport in tissues.

DOI： 10.1080/10407782.2010.516698

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

A. Safavinejad;S. H. Mansouri;A. Sakurai;S. Maruyama
In this study, a combinatorial optimization methodology has been presented for determining the optimal number and location of equally powered heaters over some parts of the boundary, called the heater surface, to satisfy the desired heat flux and temperature profiles over the design surface while keeping the total heaters power constant but floating the number of heaters. In a typical enclosure, candidate locations were numerous for placing the heaters. The optimal number and location could be found by checking among all the possible combinations of heater power ranges and locations on the heater surface. The possibility of checking only a small portion of the total search space was increasingly desirable for finding an overall optimal solution. Micro-genetic algorithm was a candidate method which displayed a significant potential in achieving that task. Micro-genetic algorithm was used to minimize an objective function which was expressed by the sum of square errors between estimated and desired heat fluxes on the design surface. Radiation element method by ray emission model (REM(2)) Was used to calculate the radiative heat flux on the design surface. It enabled us to handle the effects of specular surfaces and blockage radiation due to enclosure geometry. The capabilities of this methodology were demonstrated by finding the optimal number and position of heaters in two irregular enclosures. The effects of refractory surface characteristics (i.e., diffuse and/or specular) on the optimal solution have been studied in detail. The results show that the refractory surface characteristics have profound effects on the optimal number and location of heaters. (C) 2008 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.applthermaleng.2008.05.025

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：JAPAN SOC MECHANICAL ENGINEERS  

Atsushi Sakurai;Isami Nitta;Shigenao Maruyama;Junnosuke Okajima;Koji Matsubara
The objective of this study is to develop a new simulator for laser therapy. In order to minimize undesired effects on the laser therapy, it is important to precisely predict for heat and photon transport inside a biological tissue. In the present study, the photon transport is solved with the Radiation Element Method by Ray Emission Model (REM2), and then coupled with bioheat transfer equation. The dimensionless results become valuable guidance for developing a laser therapy system.

DOI： 10.1299/jtst.4.314

Web of Science

researchmap

Language：Japanese   Publishing type：Research paper (scientific journal)   Publisher：Japan Society of Mechanical Engineers  

Koji Matsubara;Akihiko Matsui;Takahiro Miura;Atsushi Sakurai;Hitoshi Suto;Koji Kawai
Direct numerical simulation was performed for a spatially advancing turbulent flow and heat transfer in a two-dimensional curved channel, where one wall was heated at the constant temperature and another wall was cooled similarly. In the simulation, fully developed flow and temperature computed by the straight-channel driver was given to the inlet of curved channel domain. The frictional Reynolds number was assigned 150, and the Prandtl number was given 0.71. Since the flow field was examined in the previous paper, thermal features are mainly targeted in the present paper. Turbulent heat flux showed trends consistent with growing process of large-scale vortices. In the curved part, the wall-normal component of the turbulent heat flux was twice as large as the counterpart in the straight part, suggesting active heat transport of large-scale vortices. In the inner side of the same part, temperature fluctuation was abnormally large when compared with modest fluctuation of wall-normal velocity. This was caused by the combined effect of the large-scale motion of the vortices and the monotonous variation of the mean temperature
in such a temperature distribution, large-scale ejection of the hot fluid near the outer, which is transported into the near inner-wall region, should made large impact on thermal boundary layer near the inner wall. Wave number decomposition was made for various statistics, and such effort showed that contribution of the large-scale vortex to the total turbulent heat flux normal to the wall was roughly 80 percent at the channel center of 130° downstream from the curved-channel inlet.

DOI： 10.1299/kikaib.75.754_1336

Scopus

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Atsushi Sakurai;Shigenao Maruyama;Atsuki Komiya;Koji Miyazaki

CiNii Article

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Shigenao Maruyama;Hirotaka Nakai;Atsushi Sakurai;Atsuki Komiya
Simplifications of the model for nongray radiative heat transfer analysis in participating media comprised of polydisperse water droplets are presented. Databases of the radiative properties for a water droplet over a wide range of wavelengths and diameters are constructed using rigorous Mie theory. The accuracy of the radiative properties obtained from the database interpolation is validated by comparing them with those obtained from the Mie calculations. The radiative properties of polydisperse water droplets are compared with those of monodisperse water droplets with equivalent mean diameters. Nongray radiative heat transfer in the anisotropic scattering fog layer, including direct and diffuse solar irradiations and infrared sky flux, is analyzed using REM2. The radiative heat fluxes within the fog layer containing polydisperse water droplets are compared with those in the layer containing monodisperse water droplets. Through numerical simulation of the radiative heat transfer, polydisperse water droplets can be approximated by using the Sauter diameter, a technique that can be useful in several research fields, such as engineering and atmospheric science. Although this approximation is valid in the case of pure radiative transfer problems, the Sauter diameter is reconfirmed to be the appropriate diameter for approximating problems in radiative heat transfer, although volume-length mean diameter shows better accordance in some cases. The CPU time for nongray radiative heat transfer analysis with a fog model is evaluated. It is proved that the CPU time is decreased by using the databases and the approximation method for polydisperse particulate media. (C) 2007 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.jqsrt.2007.05.006

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：JAPAN SOC MECHANICAL ENGINEERS  

Ali Safavinejad;Shigenao Maruyama;Seyed Hossein Mansouri;Atsushi Sakurai
This study presents an optimization methodology for finding the heaters setting that produce a desired heat flux and temperature distribution over a region of the enclosure surface, called the design surface. Radiation element method by ray emission model (REM2) was used to calculate the radiative heat flux on the design surface, which enable us to handle reflecting surfaces. Micro-genetic algorithm was used to minimize an objective function, which was expressed by the sum of square error between estimated and desired heat fluxes on the design surface. The novel features of this methodology were demonstrated by finding the optimal heaters setting of a three dimensional enclosure for three aspect ratio. The effects of refractory properties have been studied as well.

DOI： 10.1299/jtst.3.179

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：TAYLOR & FRANCIS INC  

Shigenao Maruyama;Atsushi Sakurai;Atsuki Komiya
A new algorithm, the discrete ordinates radiation element method (DOREM), for modeling radiative heat transfer in inhomogeneous three-dimensional participating media is described. The DOREM uses advantages of the both the radiation element method (REM) and the discrete ordinates method. Benchmark comparisons are conducted against several radiation models. The DOREM successfully implements radiative heat transfer simulations precisely, since false scattering never occurs. The DOREM has advantages of computational speed against Monte Carlo, and the CPU time and the memory size of the DOREM are 82 times faster and 767 times smaller at the maximum than that of the REM.

DOI： 10.1080/10407790600878726

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

A Sakurai;S Maruyama;S Sakai;T Nishikawa
The effects of three-dimensional (3D) radiative heat transfer in cloud fields are investigated by using the radiation element method by ray emission model (REM2). We have performed a 3D radiative heat transfer simulation in two fair-weather cumuli that are subjected to solar and sky irradiation. One of the features in this simulation is that both shortwave radiation and longwave radiation are considered. As a result, the effect of shortwave radiation is not so substantial as compared with that of longwave radiation in terms of the heat generation rate within the cloud. Also, the radiative interaction does not strongly affect the relative differences of the average heat generation rate between the single cloud model and the double cloud model. However, the distributions of heat generation rate at the thin layer beneath the cloud surface are influenced by radiative interaction and shadowing effect. Hence, these effects caused by 3D radiative heat transfer are important for horizontal inhomogeneous cloud fields. (C) 2004 Elsevier Ltd. All rights reserved.

DOI： 10.1016/j.jqsrt.2004.08.013

Web of Science

researchmap

Language：English   Publishing type：Research paper (scientific journal)  

Atsushi Sakurai;Tae-Ho Song;Shigenao Maruyama;Hyun Keol Kim
This paper compares and scrutinizes the characteristics of radiation element method by ray emission model (REM2) and discrete ordinates interpolation method (DOIM). A three-dimensional radiative heat transfer problem is solved within a rectangular enclosure containing absorbing, emitting, nonscattering medium with temperatures given as a sine function. The calculation results are validated with the zonal results and the analytical solution. Generally, the calculation results agree well with the exact solution. However, when the optically thickness is much larger than unity, the zonal and REM2 results overestimate the heat source strength and wall heat flux, and when it is small, DOIM results underestimate them. The reasons and remedies are carefully discussed.

DOI： 10.1299/jsmeb.48.259

Scopus

researchmap

Language：Japanese Book type：Scholarly book

researchmap

Authorship：Corresponding author  

Metal-Insulator-Metal(MIM)構造を用いた熱輻射メタマテリアルの設計・作製方法に関する研究と、機械学習による超狭帯域熱輻射メタマテリアルの実現に関する研究について紹介する。

researchmap

Authorship：Lead author,　Last author,　Corresponding author  

researchmap

researchmap

researchmap

Language：English   Publisher：PERGAMON-ELSEVIER SCIENCE LTD  

Zhuomin Zhang;Shigenao Maruyama;Atsushi Sakurai;M. Pinar Menguc

DOI： 10.1016/j.jqsrt.2013.06.012

Web of Science

researchmap

Event date： 2021

researchmap

Event date： 2020.10

Language：Japanese  

In a semiconductor light emitting device such as an LED, when a voltage is applied to increase the chemical potential of photons, a thermodynamically nonequilibrium state is generated, and light emission may be different from Planck's law black body radiation. In this study, we focus on Thermophotonic power generation system composed of LED and PV cell, and discuss the mechanism of non-equilibrium thermal radiation via fluctuational electromagnetic simulation.

researchmap

Event date： 2020

researchmap

Event date： 2020

researchmap

Event date： 2020

researchmap

Event date： 2019.12

researchmap

Event date： 2019.11

researchmap

Event date： 2019.11

researchmap

Event date： 2019.11

researchmap

Event date： 2019.10

researchmap

Event date： 2019.8

The concentrated solar power is utilized for industrial thermal usage, electricity generation and fuel production. To increase the solar thermal temperature is essential for enhancing the thermal efficiency of the energy conversion and for widening application. The synthetic oil, widely used for solar thermal storage, pylorizes at 400 °C, and this cannot be used for very high temperatures. The authors used quartz sand as thermal storage medium in the fluidized bed receiver, and achieved the thermal receiver temperatures beyond 1000 °C in the previous report. The research work was further made on using the particle steady flow for high-temperature solar thermal storage. The steady-flow-type particle receiver was designed through visualization of cold model and numerical simulation. The experimental apparatus was fabricated and tested using beam-down sun simulator of 1.0 kWth. The apparatus comprises of the internal-circulating fluidized bed receiver placed between a loop seal with an inlet port and another loop seal with a weir and an exit port of the particles. The tested particles are quartz sand of 0.1 - 0.2 mm and fused brown alumina of 0.09 - 0.1 mm. The experimental set up produces hot particles of quartz sand at 428°C and those of fused brown alumina at 491°C. The highly efficient steady-flow-type particle receiver created steady particle flow smoothly. This apparatus takes in the concentrated light through the quartz window to continuously discharge the particles at high temperatures beyond the criterion in synthesis oil (400°) for solar thermal storage.

researchmap

Event date： 2019.8

researchmap

Event date： 2019.7

researchmap

Event date： 2019.4

researchmap

Event date： 2019.1

researchmap

Event date： 2019

Language：Japanese  

In this paper, we conduct TPV power generation experiment using wavelength selective emitter. This experiments were conducted in a vacuum chamber to ignore the effect of convection. We use emitter is a MIM structure due to metal-dielectric-metal. This emitter shows high emissivity at shorter wavelength than 1.75 um, which is the high sensitivity wavelength range of GaSb PV cell. In the experiment, blackbody and selective radiant material are used for the emitter, and the power generation experiment is performed by directly radiating the emitter to the PV cell by heating the emitter directly with the heater. The results showed the usefulness of selective emitter.

researchmap

Event date： 2019

researchmap

Event date： 2019

Language：Japanese  

In the present study, we propose a wavelength selective emitter by using Graphene Metasurface. Graphene Metasurface consists of graphene ribbon array, potassium bromide thin film and silver substrate. Geometric parameters of Graphene Metasurface is determined based on equivalent circuit model that agrees well with the result of electromagnetic field analysis (Rigorous Coupled-Wave Analysis). The proposed emitter causes impedance matching depending on the conductivity of the graphene ribbon at very narrow wavelength range. Therefore, the proposed emitter shows emissivity of 100%, narrowband peaks with high quality factor and active controllable switching, for various wavelengths. This study may facilitate the initial design of wavelength selective emitter by Graphene Metasurface.

researchmap

Event date： 2019

researchmap

Event date： 2019

researchmap

Event date： 2019

researchmap

Event date： 2019

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

Language：Japanese  

In this paper we created ultra-narrowband wavelength selective emitter with Q-factor exceeding 180. The devices are designed by Bayesian optimization. The mechanisms of radiative transfer emission are clarified by using Rigorous Coupled-Wave Analysis. The optimized structure consists of aperiodic multiple layers of alternating materials with varying refractive index. The obtained structure was fabricated by sputtering and measured by fourier transform infrared spectroscopy. As a result, we demonstrated high and sharp emissivity in the infrared region.

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

researchmap

Event date： 2018

The technology for drying at a low temperature is required in the drying industry because flammable solvents are ignited and the heat damages are generated in a solute by a high temperature. Recently, the wavelength selective heater with a band pass filter has been developed to prevent from the ignition. However, the big heat loss occurs at the filter. In this study, the wavelength selective emitter with the metamaterial of a metal-insulator-metal (MIM) is developed to prevent from the big heat loss. The MIM emitter has the metamaterial of Au-Al2O3-Au and can emit the radiation with the normal emissivity peak over 0.8 at the wavelength of 6.60 µm. The experiments for evaporating toluene are carried out in the cases of using the MIM emitter in a vacuum and a black paint with a band pass filter in an atmospheric pressure and a vacuum. Toluene has an absorption band at the wavelength of 6.68 µm. It has been found out that the heat required for evaporating one gram of toluene decreases 39.7 % and 21.4 % in the case of using the MIM emitter in a vacuum than using the black paint with a band pass filter in an atmospheric pressure and a vacuum, respectively. Numerical analyses are clarified that the MIM emitter has non-negligible heat loss and the thermal management allows to decrease the heat loss moreover. It is concluded that the MIM emitter can contribute to decrease big heat loss at the wavelength selective heater with the filter.

researchmap

Event date： 2017.11

Language：English  

researchmap

Event date： 2017

Language：English  

Koji Matsubara;Atsushi Sakurai;Takahiro Suzuki;Selvan Bellan;Nobuyuki Gokon;Cho Hyun-Seok;Tatsuya Kodama

researchmap

Event date： 2017

A radiator is a thermal control device which radiates waste heat from the spacecraft to space and cools equipment. Since an infrared space telescope observes weak infrared signals, it is important to shift thermal radiation wavelength to longer values than the observation signals and to reduce the radiation intensity by cooling the equipment. Especially, the next generation infrared space telescope will be cooled at cryogenic temperatures so as to enhance the observation performance. However, a conventional radiator material, such as black paint, shows lower emissivity at cryogenic temperatures. That is why a new radiator with high emissivity at cryogenic temperatures is required. Therefore, we focused on a structure consisting of periodic array of metallic elements separated from a ground plane by a dielectric spacer layer in order to realize a high emittance radiator at cryogenic temperatures. In this work, we reported the performance of the IR emitter with periodic array. It was designed by the FDTD method. The periodic structure is a circle. The metal layer was made by aluminum and the dielectric layer was made by SiO2. The prototype IR emitters were fabricated. Spectral reflectance of the prototype emitter was measured by FT-IR. Moreover, spectral emittance was measured by the blackbody comparison method, and it was compared with spectral absorptance calculated by spectral reflectance.

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2017

researchmap

Event date： 2016

researchmap

Event date： 2016

researchmap

Event date： 2016

researchmap

Event date： 2016

researchmap

Event date： 2016

researchmap

Event date： 2016

researchmap

Event date： 2016

Language：Japanese  

In the present study, we numerically demonstrate wavelength-selective thermal radiation emitter made of graphene-covered Ag grating and it filled with silicon. Underlying physical mechanism is clarified by rigorous coupled-wave analysis. We find that the emission peaks associated with magnetic polariton (MP), whose resonance wavelength can be dynamically tuned by varying graphene chemical potential. This study could facilitate the new tunable thermal radiation sources in the infrared region.

researchmap

Event date： 2016

researchmap

Event date： 2015

researchmap

Event date： 2015

researchmap

Event date： 2015

researchmap

Event date： 2015

researchmap

Event date： 2015

researchmap

Event date： 2015

researchmap

Event date： 2015

researchmap

Event date： 2015

Language：Japanese  

Technologies for drying a material with a flammable solvent at a low temperature are required. It is reasonable to evaporate a flammable solvent using radiation at the absorption band of a flammable solvent and to eject the flammable vapor with a cold air. The meta material with Au/Al2O3/Au was fabricated by using an sputtering equipment, an atomic layer deposition equipment, and an electron beam lithography equipment. The reflectances of the meta material structure were measured by using an FT-IR equipment with an integral sphere. It is clarified that the emissivities of the meta material structure are larger than 0.8 at the absorption band of toluene.

researchmap

Event date： 2015

researchmap

Event date： 2014.11

Language：Japanese  

We numerically demonstrate wavelength-selective absorption in core-shell deposited on top of a continuous Au film by excitation of magnetic resonance. The resonance of magnetic polariton at wavelength 653 nm is confirmed by electromagnetic Geld distribution. The influences of Au core diameter, SiO_2 shell thickness on the resonance wavelength are studied with 3D finite-differential time-domain (FDTD) method. The fundamental study will facilitate the development of a core-shell metamaterial for energy and sensing applications.

researchmap

Event date： 2014.10

Language：English  

Makoto Yamada;Atsushi Sakurai

researchmap

Event date： 2014.6

Language：English  

Atsushi Sakurai;Bo Zhao;Zhuomin M. Zhang

researchmap

Event date： 2014.1

Atsushi Sakurai;Bo Zhao;Zhuomin M. Zhang

researchmap

Event date： 2014

researchmap

Event date： 2014

researchmap

Event date： 2014

researchmap

Event date： 2013.11

Language：Japanese  

Direct numerical simulation was performed for a spatially advancing turbulent flow in a two-dimensional curved channel. The radius ratio of the curved part, α , was set at 0.92, the same as Kobayashi et al.'s experiment. And the friction Reynolds number, Re_<τ0>, was set at 550. Numerically solved mean velocity field showed trends consistent with the experiment. As a result, the validity of the simulation was confirmed. Instantaneous flow field implied that micro-scale structures near the outer wall were related to the birth and growth of large-scale vortices.

researchmap

Event date： 2013.8

Language：Japanese  

researchmap

Event date： 2013.6

Language：English  

Atsushi Sakurai;Hiroya Tanikawa;Makoto Yamada

researchmap

Event date： 2013.3

Language：Japanese  

researchmap

Event date： 2013.3

Language：Japanese  

researchmap

Event date： 2013.3

Language：Japanese  

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2013

researchmap

Event date： 2012.11

Language：Japanese  

It's necessary to further improve technology of energy-saving with increasing energy consumption. And there is a current situation that low-temperature exhaust heat has been discarded from a factory. So we have started a study of the heat pipe for utilizing this exhaust heat. It's necessary to carry out heat transportation of long distance and high density for utilizing the exhaust heat. Therefore the present study had been used the loop thermosyphon. But quantities of the water which are enclosed in loop thermosyphon will affect heat transportation. So the present study has been practiced the experiment which had been changed quantities of the water and thermal dose. And we have researched the effect of heat transportation. As the result, heat transportation is improved with increasing thermal dose.

researchmap

Event date： 2012.9

Language：Japanese  

An experimental study was performed for freeze-concentration of 1 % water-solution of sodium chloride in a cylindrical container, where a heat transfer pipe was set vertically. Experimental results denoted that temperature gradient in the ice promote solute rejection. When the ice is given a temperature gradient for long time , solute is eliminated little by little from the ice . By varying the temperature of the inner wall, the difference of concentration rate, distribution coefficient and concentration of ice was measured.

researchmap

Event date： 2012.3

Language：Japanese  

researchmap

Event date： 2012.3

Language：Japanese  

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2012

Language：Japanese  

Double-walled reactor tubes with molten-salt thermal storage were developed for solar air receiver to generate electricity and for solar reforming to produce solar fuel. A mixture of Na_2CO_3 and MgO was loaded as a thermal storage medium into the outer annulus of that. Performances of the reactor tubes were studied as an air receiver and as a reformer tube in the present study.

researchmap

Event date： 2012

researchmap

Event date： 2012

Language：English  

Atsushi Sakurai;So Sakuma;Gokon Nobuyuki;Koji Matsubara;Kodama Tatsuya

researchmap

Event date： 2012

researchmap

Event date： 2012

researchmap

Event date： 2011.11

Language：English  

A.Sakurai;Y.Sato;S.Maruyama;J.Okajima;A.Komiya

researchmap

Event date： 2011.11

Language：English  

N.Yamada;A.Sakurai;A.Komiya;S.Maruyama

researchmap

Event date： 2011.11

Language：English  

So Sakuma;Atsushi Sakurai;Kyohei Ogino;Seung-Jae Lee;Koji Matsubara;Nobuyuki Gokon;Tatsuya Kodama

researchmap

Event date： 2011.11

Language：English  

researchmap

Event date： 2011.10

Language：Japanese  

Direct numerical simulation was performed for a spatially advancing turbulent flow in a two-dimensional curved channel. The radius ratio of the curved part, α, was set at 0.92, the same as Kobayashi et al.'s experiment. And the friction Reynolds number, Re_<τ0>, was set at 550 and 300. Numerically solved mean velocity field showed trends consistent with the experiment. As a result, the validity of the simulation was confirmed. Instantaneous flow field implied that micro-scale structures near the outer wall were related to the birth and growth of large-scale vortices.

researchmap

Event date： 2011.9

Language：English  

Direct numerical simulation was performed for a spatially advanced turbulent flow in a two-dimensional curved channel. In the past research, the influence of radius ratio was not discussed. Therefore, the radius ratio of the curved part, a=0.92, the same as Kobayashi et al.'s experiment, and 0.8 were calculated. Then, we carried out simulations of the temperature condition where two walls were sustained at different temperatures. Prandtl number was assumed to be 0.71. Analogy between the friction factor and the Stanton number was confirmed through computed cases. Turbulence statistics including the turbulent heat flux budgets were presented and discussed.

researchmap

Event date： 2011.9

Language：English  

Nobuyuki Gokon;Ryuta Ono;Tsuyoshi Hatamachi;Liuyun Li;Hee Joon Kim;Atsushi Sakurai;Koji Matsubara;Masayuki Yagi;Tatsuya Kodama

researchmap

Event date： 2011.9

Language：English  

Atsushi Sakurai;So Sakuma;Nobuyuki Gokon;Koji Matsubara;Tatsuya Kodama

researchmap

Event date： 2011.7

Language：English  

Takahiro Miura;Koji Matsubara;Atsushi Sakurai

researchmap

Event date： 2011.6

Language：English  

Keisuke Sato;Koji Matsubara;Atsushi Sakurai

researchmap

Event date： 2011.6

Language：English  

O. Fukuda;T. Yamaguchi;K. Matsubara;A. Sakurai;Y. Aoki;M. Maki;F. Taguchi

researchmap

Event date： 2011.6

Language：English  

T. Miura;K. Matsubara;A. Sakurai

researchmap

Event date： 2011.6

Language：English  

S.Hanashima;A.Sakurai;K.Matsubara

researchmap

Event date： 2011.6

Language：English  

S. Sakuma;A. Sakurai;K. Matsubara;N. Gokon;T. Kodama

researchmap

Event date： 2011.3

Language：Japanese  

researchmap

Event date： 2011.1

Language：Japanese  

researchmap

Event date： 2011

researchmap

Event date： 2011

Importance of turbulence and radiation interaction (TRI) has been investigated in a turbulent channel flow by using direct numerical simulation (DNS) to clarify detailed turbulent flow structure and heat transfer mechanisms. To investigate the effect of correlation functions between gas absorption and temperature fluctuation, the two cases of correlation are tested. Consequently, the TRI effect can be clearly observed when the correlation is positive. This fact provides the evidence that radiative intensity is enhanced by the turbulent fluctuation. The DNS results suggest the significance in the fundamental aspect of TRI. Furthermore, effects of frictional Reynolds number, Reô , are investigated. Comparing with the case of Reô =150, the location of the enhancement peaks of Reô =300 shifts toward the walls. It is found that the relative importance of the TRI correspond to the structure of temperature fluctuation intensity originated from the differences of the Reτ. Copyright © 2011 by ASME.

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

Language：Japanese  

researchmap

Event date： 2011

Language：Japanese  

researchmap

Event date： 2011

Language：Japanese  

researchmap

Event date： 2011

Language：English  

Takashi Uchida;Koji Matsubara;Takahiro Miura;Atsushi Sakurai

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2011

Direct numerical simulations were performed for turbulent flows in the spatial developing region of a ribbed channel, where square ribs were attached to one wall. The rib pitch-to-height ratio was changed at three steps; Pi/H = 2, 4 and 7. The bulk Reynolds number was given 4,560 through simulations. For the same parameters, simulations were also made for periodic fully developed cases. Although the friction factor decreased as flow goes downstream in the case of Pi/H = 2, it increased after unsteadiness in the case of Pi/H = 7 due to the turbulence enhancement. For Pi/H = 7, the developing flow near the rough wall even in the second cavity had a similar structure to the fully developed flow, whereas developed structures extended further away from the ribbed wall. The vortex structures were revealed by instantaneous fields and the two-point correlation of vorticities. Streamwise and spanwise vortices appeared near the rough wall. The streamwise vortex was smaller than that in a smooth channel, especially in the streamwise direction, and the inclination angle of the vortex was larger. Although the flow near the ribs was more isotropic than the smooth wall, coherent structures existed and produced the high correlation between the streamwise and wall-normal velocity fluctuations.

researchmap

Event date： 2011

researchmap

Event date： 2011

researchmap

Event date： 2010.11

Language：English  

Atsushi Sakurai;Koji Matsubara;Shigenao Maruyama

researchmap

Event date： 2010.11

Language：English  

Ryo Kanbayashi;Atsushi Sakurai;Koji Matsubara

researchmap

Event date： 2010.11

Language：English  

Shouta Hanashima;Atsushi Sakurai;Koji Matsubara

researchmap

Event date： 2010.10

Language：Japanese  

This paper is two dimensional numerical simulations of solid-liquid flow in straight micro-channel. It is shown that mixing effect of particle by applying AC electric field. Mixing is one of the important operations of μ-TAS. However, flow in micro-channel is laminar flow because of low Reynolds number. Therefore, mixing in micro-channel is difficult. Electric potential is applied at four channel wall electrode. Particle mixing is encouraged by interaction of forced convection and electro-osmotic flow. Mixing characteristics of particle is depends on frequency of AC electric field.

researchmap

Event date： 2010.10

Language：Japanese  

An experimental study was performed for freeze-concentration of 1% or 10% water-solution of sodium chloride in a cylindrical container, where a heat transfer pipe was set vertically. Experimental results denoted that concentration in case of non-insulation is better than the case of insulation for outside of a container. In case of keeping the temperature at inner wall constant, distribution factor shown better results as preset temperature is raised. In terms of a relation of freezing speed and coefficient of distribution, high temperature at inner wall prevent ice surface from figuration of mushy region. Though measurement of density distribution, it was suggested that natural convection enhances the concentration.

researchmap

Event date： 2010.10

Language：Japanese  

Direct numerical simulation was performed for the coexistent convection in a channel affected by the heat radiation. Four kinds of numerical condition were targeted by changing the Grashof number, Gr, and the optical thickness, T; the pure forced convection (Gr=0, τ=0), only buoyancy (Gr/10^6=1.3, τ=0) and the buoyancy with radiation (Gr/10^6=1.3, τ=0.01), (Gr/10^6=1.3, τ=0.1). According to the numerical result, buoyancy effects created the large-scale organized vortex flow all over the channel, which enhanced dramatically heat transfer. The effect of the radiation weakens the large-scale organized vortex and reduces the heat transfer by the buoyancy.

researchmap

Event date： 2010.3

Language：Japanese  

researchmap

Event date： 2010.3

Language：Japanese  

researchmap

Event date： 2010

Language：Japanese  

researchmap

Event date： 2010

Language：English  

Takahiro Miura;Koji Matsubara;Yoshiyuki Nagai;Atsushi Sakurai
Direct numerical simulations were performed for air flows and related heat transfer in a developing region of the ribbed channel, where five ribs were attached to one wall. Three rib pitch-to-height ratios, Pi/H = 7, 4 and 2, were considered. The Reynolds number based on the bulk velocity and the channel width was 4,560. For the same parameters, simulations were also made for periodic fully developed cases. Although the largest rib pitch, Pi/H = 7, was the best choice for a heat transfer enhancement ratio in the fully developed cases, the enhancement ratio was larger for the smaller values of rib pitch, Pi/H = 4 or Pi/H = 2 in the initial region. Preference of smaller rib spacing in the initial region was found to come from: (I) the pressure drop was kept low due to the modest meandering of the stream lines; (2) downstream ribs were placed immediately behind the upstream ribs being exposed to the thin thermal layer; (3) the turbulent heat transport was active because of the fluctuation based on the Kelvin-Helmholtz type instability in Pi/H = 4 though fluctuation was weak in Pi/H = 2. Therefore, flow structures positively contributed the heat transfer enhancement of Pi/H = 4 and 2 in the developing region, and the smaller values of the rib pitch was recommended to be selected in such situations.

researchmap

Event date： 2010

Language：Japanese  

Direct numerical simulation was performed for a spatially advancing turbulent flow in a two-dimensional curved channel. The radius ratio of the curved part, α, was set at 0.92, the same as Kobayashi et al.'s experiment. However, the friction Reynolds number, Re_<τ0>, was assigned at 300, which was roughly half the experiment. Numerically solved mean velocity showed trends consistent with the experiment in spite of the Reynolds number difference. Numerical data of mean velocity field illustrated that large scale vortices grew along curved channel. Instantaneous flow field and power spectrum analysis implied that micro structures near the outer wall were related to the birth and formation of large scale vortices.

researchmap

Event date： 2010

researchmap

Event date： 2010

researchmap

Event date： 2010

researchmap

Event date： 2010

researchmap

Event date： 2010

researchmap

Event date： 2010

researchmap

Event date： 2010

researchmap

Event date： 2010

Language：Japanese  

researchmap

Event date： 2010

Language：Japanese  

researchmap

Event date： 2009.10

Language：Japanese  

researchmap

Event date： 2009.4

Language：English  

Atsushi Sakurai;Subhash C. Mishra;Shigenao Maruyama

researchmap

Event date： 2009

researchmap

Event date： 2009

Language：Japanese  

researchmap

Event date： 2009

Language：Japanese  

researchmap

Event date： 2009

Language：Japanese  

Direct numerical simulation was performed for a spatially developing turbulent flow in a channel with a rectangular rib. Inflow condition was a fully developed channel flow for a Reynolds number of 150, based on the friction velocity and the channel half-width. The budgets of transport equations of Reynolds stress are indicated. In the upstream of the rib, the streamwise velocity fluctuations increase by the deceleration of the mean streamwise velocity. Near the front of the rib, very large splatting effect observes. Near the upper part of the rib, the production term by the deceleration of the mean streamwise velocity is dominant.

researchmap

Event date： 2009

researchmap

Event date： 2009

researchmap

Event date： 2009

Language：English  

K. Matsubara;H. Suto;T. Miura;A. Sakurai
Simulation was made for non-isothermal air jet of the Reynolds number set 1200, 10,000, 100,000 and 1,000,000 to reveal coherent structures in its developed regime. For the case of low Reynolds number, visualization was made by dye contamination and PTV, and experimental results supported the numerical validity. For the cases of higher Reynolds number, numerical result was confirmed to agree with existing literatures. Coherent structures are examined through three-dimensional mapping of vortex tube, instantaneous distribution of velocity vectors and temperature, two-point correlation of velocity and temperature, and PDF of vorticity vector angle. Conceptual model of coherent structures in developed jet was thus proposed; (1) There appeared frequently hairpin shape structures in developed jet, and they enhanced momentum transport and scalar mixing: (2) Hairpin-shape structures are dominant for low Reynolds number, but incompleteness of structures grows as the Reynolds number increases.

researchmap

Event date： 2009

researchmap

Event date： 2009

researchmap

Event date： 2009

researchmap

Event date： 2009

researchmap

Event date： 2009

researchmap

Event date： 2007.4

Language：English  

Atsushi Sakurai;Shigenao Maruyama;Koji Miyazaki;Atsuki Komiya

researchmap

Event date： 2007.4

Language：English  

Ali Safavinejad;Seyed Hossein Mansouri;Shigenao Maruyama;Atsushi Sakurai

researchmap

Event date： 2007.4

Language：English  

Stephanie Veran;Michel Quintard;Shigenao Maruyama;Atsushi Sakurai

researchmap

Event date： 2007

researchmap

Event date： 2006.4

Language：English  

Atsushi Sakurai;Masud Behnia;Shigenao Maruyama;Atsuki Komiya;Seigo Sakai

researchmap

Event date： 2006

researchmap

Event date： 2006

researchmap

Event date： 2006

researchmap

Event date： 2005

researchmap

Event date： 2004

researchmap

Event date： 2004

Language：Japanese  

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

Presentation type：Oral presentation (general)  

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

Applicant：国立研究開発法人物質・材料研究機構

Application no：特願2018-087363  Date applied：2018.4

Announcement no：特開2019-192608  Date announced：2019.10

J-GLOBAL

researchmap

Applicant：国立研究開発法人物質・材料研究機構

Application no：特願2018-087363  Date applied：2018.4

Announcement no：特開2019-192608  Date announced：2019.10

Patent/Registration no：特許第7066120号  Date registered：2022.5 

J-GLOBAL

researchmap

Applicant：日本碍子株式会社

Application no：JP2017037742  Date applied：2017.10

Publication no：WO2018-079386  Date published：2018.5

J-GLOBAL

researchmap

Applicant：北海道大学、新潟大学、日本碍子株式会社

Application no：特願2016-207571  Date applied：2016.10

researchmap

Applicant：北海道大学、新潟大学、日本碍子株式会社

Application no：特願2015-175068  Date applied：2015.9

Announcement no：特開2017-050254  Date announced：2017.3

Patent/Registration no：特許6692046  Date registered：2020.4 

researchmap

Applicant：国立大学法人北海道大学

Application no：特願2015-175068  Date applied：2015.9

Announcement no：特開2017-050254  Date announced：2017.3

J-GLOBAL

researchmap

Applicant：国立大学法人北海道大学

Application no：特願2015-175068  Date applied：2015.9

Announcement no：特開2017-050254  Date announced：2017.3

Patent/Registration no：特許第6692046号  Date registered：2020.4 

J-GLOBAL

researchmap

Applicant：国立大学法人 新潟大学

Application no：特願2013-222867  Date applied：2013.10

Announcement no：特開2015-086232  Date announced：2015.5

Patent/Registration no：特許第6232923号  Date registered：2017.11 

J-GLOBAL

researchmap

Applicant：国立大学法人 新潟大学

Application no：特願2013-222867  Date applied：2013.10

Announcement no：特開2015-086232  Date announced：2015.5

J-GLOBAL

researchmap

Applicant：国立大学法人 新潟大学

Application no：特願2014-534366  Date applied：2013.9

Patent/Registration no：特許第6165743号  Date registered：2017.6 

J-GLOBAL

researchmap

Applicant：国立大学法人 新潟大学

Application no：JP2013073693  Date applied：2013.9

Publication no：WO2014-038553  Date published：2014.3

J-GLOBAL

researchmap

Applicant：国立大学法人 新潟大学

Application no：JP2013052211  Date applied：2013.1

Publication no：WO2013-115316  Date published：2013.8

J-GLOBAL

researchmap

Applicant：国立大学法人 新潟大学

Application no：特願2013-556494  Date applied：2013.1

Patent/Registration no：特許第5986589号  Date registered：2016.8 

J-GLOBAL

researchmap

Application no：特願2018-082171 

Patent/Registration no：特許6977943  Date registered：2021.11 

researchmap

Application no：特願2018-087363 

Patent/Registration no：特許7066120  Date registered：2022.5 

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

Award type：International academic award (Japan or overseas)  Country：Japan

researchmap

researchmap

Country：Japan

researchmap

researchmap

researchmap

Award type：Award from Japanese society, conference, symposium, etc.  Country：Japan

researchmap

Country：Japan

researchmap

Award type：Award from Japanese society, conference, symposium, etc.  Country：Japan

researchmap

researchmap

researchmap