Researcher profile

Shu Takagi

Shu Takagi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

An Eulerian-based immersed boundary method for particle suspensions with implicit lubrication model

We describe an immersed boundary method in which the fluid-structure coupling is achieved in an Eulerian framework. The method is an improved extension of the immersed boundary method originally developed by Kajishima et al. [1], which accounts for the inertia of the fictitious fluid inside the particle volume and is thus able to reproduce the behaviour of particles both in the case of neutrally-buoyant objects and in the presence of density difference between the particles and the fluid. The method is capable to handle the presence of multiple suspended objects, i.e., a suspension, by including a soft-sphere normal collision model, while the lubrication correction typically added to similar immersed boundary methods in order to capture the sub-grid unresolved lubrication force is here treated implicitly, i.e., naturally obtained without any explicit expression, thus no additional computation is required. We show that our methodology can successfully reproduce the rheology of a particle suspension in a shear flow up to a dense regime (with a maximum particle volume fraction around 46%) without any additional correction force. The applicability of this methodology is also tested in a turbulent pressure-driven duct flow at high Reynolds number in the presence of non-negligible inertia and non-uniform shear-rate, showing good agreement with experimental measurements.

preprint2022arXiv

Gas Adsorption and Diffusion Behaviors in Interfacial Systems Composed of a Polymer of Intrinsic Microporosity and Amorphous Silica: A Molecular Simulation Study

We investigate the adsorption and diffusion behaviors of CO2, CH4, and N2 in interfacial systems composed of a polymer of intrinsic microporosity (PIM-1) and amorphous silica using grand canonical Monte Carlo (GCMC) and molecular dynamics (MD) simulations. We build model systems of mixed matrix membranes (MMMs) with PIM-1 chains sandwiched between silica surfaces. Gas adsorption analysis using GCMC simulations shows that gas molecules are preferentially adsorbed in microcavities distributed near silica surfaces, resulting in an increase in the solubility coefficients of CO2, CH4, and N2 compared to bulk PIM-1. In contrast, diffusion coefficients obtained from MD simulations and then calibrated using the dual-mode sorption model show different tendencies depending on gas species: CO2 diffusivity decreases in MMMs compared to PIM-1, whereas CH4 and N2 diffusivities increase. These differences are attributed to competing effects of silica surfaces: the emergence of larger pores as a result of chain packing disruption, which enhances gas diffusion, and a quadrupole-dipole interaction between gas molecules and silica surface hydroxyl groups, which retards gas diffusion. The former has a greater impact on CH4 and N2 diffusivities, whereas the latter has a greater impact on CO2 diffusivity due to the strong quadrupole-dipole interaction between CO2 and surface hydroxyls. These findings add to our understanding of gas adsorption and diffusion behaviors in the vicinity of PIM-1/silica interfaces, which are unobtainable in experimental studies.

preprint2021arXiv

Molecular Insights into the Mechanical Properties of Polymer-Fullerene Bulk Heterojunctions for Organic Photovoltaic Applications

We investigate the mechanical properties of $π$-conjugated polymeric materials composed of regioregular poly(3-hexylthiophene) (P3HT) and fullerene C$_{60}$ using coarse-grained molecular dynamics simulations. Specifically, we perform tensile simulations of P3HT:C$_{60}$ composites with varied degrees of polymerization and C$_{60}$ mass fractions to obtain their stress-strain responses. Decomposition of stress tensor into kinetic energy and virial contributions indicates that the tensile moduli of the pure P3HT samples are greatly dependent on non-bonded interactions and on bonded interactions associated with bond-stretching, while the addition of C$_{60}$ leads to an increase in the tensile modulus originating from enhanced non-bonded interactions associated with C$_{60}$. Additionally, the tensile strength of the P3HT:C$_{60}$ samples correlates well with molecular chain entanglements, which are characterized by the average number of kinks per chain obtained from primitive path analysis. We also find that the upper and lower yield points characterizing strain softening become more pronounced with increasing C$_{60}$ mass fraction. Persistent homology analysis indicates that the emergence of the yield points correlates well with the coalescence of microvoids in the course of tensile deformation, resulting in the generation of larger voids. These results provide a fundamental understanding of the molecular determinants of the mechanical properties of $π$-conjugated polymer-fullerene composites, which can also help to interpret and predict the mechanical properties of other polymer composites containing fullerene.

preprint2020arXiv

On viscid-inviscid interactions of a pair of bubbles rising near the wall

Series of experiments on turbulent bubbly channel flows observed bubble clusters near the wall which can change large-scale flow structures. To gain insights into clustering mechanisms, we study the interaction of a pair of spherical bubbles rising in a vertical channel through combined experiments and modeling. Experimental imaging identifies that pairwise bubbles of 1.0 mm diameter take two preferred configurations depending on their mutual distance: side-by-side positions for a short distance ($S<5$) and nearly inline, oblique positions for a long distance ($S>5$), where $S$ is the mutual distance normalized by the bubble radius. In the model, we formulate the motions of pairwise bubbles rising at $Re=O(100)$. Analytical drag and lift, and semi-empirical, spatio-temporal stochastic forcing are employed to represent the mean acceleration and the fluctuation due to turbulent agitation, respectively. The model is validated against the experiment through comparing Lagrangian statistics of the bubbles. Simulations using this model identify two distinct timescales of interaction dynamics which elucidate the preferred configurations. For pairs initially in-line, the trailing bubble rapidly escapes from the viscous wake of the leading bubble to take the oblique position. Outside of the wake, the trailing bubble travels on a curve-line path with a slower velocity driven by potential interaction and horizontally approaches the leading bubble to become side-by-side. Moreover, statistical analysis identifies that the combination of the wake and the agitation can significantly accelerate the side-by-side clustering of in-line pairs. These results indicate positive contributions of liquid viscosity and turbulence to the formation of bubble clusters.

preprint2020arXiv

Relating the thermal properties of a micro pulsating heat pipe to the internal flow characteristics via experiments, image recognition of flow patterns and heat transfer simulations

We investigate the relationship between the thermal properties of a micro pulsating heat pipe (MPHP) and the internal flow characteristics. The MPHP consists of an eleven-turn closed-loop of a meandering square microchannel with a hydraulic diameter of $350\ {}μ{\rm m}$ engraved on a silicon substrate. The MPHP charged with Fluorinert FC-72 tends to exhibit higher effective thermal conductivities for the coolant temperature of $T_{\rm c} = 40\ {}^\circ\mathrm{C}$ compared to $T_{\rm c} = 20\ {}^\circ\mathrm{C}$, and provides the highest effective thermal conductivity of about $700\ {}{\rm W/(m{\cdot}K)}$ for $T_{\rm c} = 40\ {}^\circ\mathrm{C}$ and a filling ratio of 48%. Interestingly, we observe two different self-oscillation modes having different thermal conductivities, even for identical heat input rates. This tendency indicates a hysteresis of the effective thermal conductivity, which originates from the difference in the heat input rates at which the MPHP falls into and recovers from dryout. Subsequently, semantic segmentation-based image recognition is applied to the recorded flow images to identify the flow characteristics, successfully extracting four different flow patterns involving liquid slugs, liquid films, dry walls, and rapid-boiling regions. The image recognition results indicate that high effective thermal conductivities of the MPHP relate to stable self-oscillations with large amplitudes and high frequencies, along with long and thin liquid films beneficial for latent heat transfer. Finally, we perform numerical simulations of latent/sensible heat transfer via vapor plugs and of sensible heat transfer via liquid slugs using the extracted flow patterns as inputs. We find that latent heat transfer via liquid films accounts for a considerable portion of the overall heat transfer, while the sensible heat transfer via liquid slugs is much less significant.

preprint2008arXiv

The lipid bilayer at the mesoscale: a physical continuum model

We study a continuum model of the lipid bilayer based on minimizing the free energy of a mixture of water and lipid molecules. This paper extends previous work by Blom & Peletier (2004) in the following ways. (a) It formulates a more physical model of the hydrophobic effect to facilitate connections with microscale simulations. (b) It clarifies the meaning of the model parameters. (c) It outlines a method for determining parameter values so that physically-realistic bilayer density profiles can be obtained, for example for use in macroscale simulations. Points (a)-(c) suggest that the model has potential to robustly connect some micro- and macroscale levels of multiscale blood flow simulations. The mathematical modelling in point (a) is based upon a consideration of the underlying physics of inter-molecular forces. The governing equations thus obtained are minimized by gradient flows via a novel numerical approach; this enables point (b). The numerical results are shown to behave physically in terms of the effect of background concentration, in contrast to the earlier model which is shown here to not display the expected behaviour. A &#34;short-tail&#34; approximation of the lipid molecules also gives an analytical tool which yields critical values of some parameters under certain conditions. Point (c) involves the first quantitative comparison of the numerical data with physical experimental results.