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

The characteristic of two-phase immiscible flow in porous media is controlled by different kinds of force, such as interface surface tension (capillary force) and fluid viscosity (viscous force). Based on those two forces, three typical displacement patterns in porous media are divided, named as capillary fingering, viscous fingering and stable displacement. However, the impact of inertial force in displacement pattern, which generate due to the flow direction or velocity change of fluid, is always neglected. This research used direct numerical simulation (DNS) to study the two-phase immiscible pattern with a wide range of capillary number (Ca) and ratio of Reynolds number (Re). The Ca is a dimensionless value which represents the ratio of viscous force and capillary force. While Re represents the ratio of inertial force and viscous force. The impact of forces on displacement patterns were determined based on the quantitative analyses of the saturation distribution as functions of Ca and Re. The result shows that inertial effects have minimal influence on flow conditions at low capillary numbers. However, at high capillary numbers, capillary forces become less significant, and inertial effects strongly influence flow conditions. These findings contribute the different insight of fluid displacement patterns, controlled by balance of inertial, capillary and viscous forces, has a noticeable influence on recovery or storage efficiency in subsurface process.

DOI： 10.26855/jamc.2024.12.006

researchmap

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

In fluid dynamics, the generation and bursting of surface bubble liquid films in surfactant-laden environments involve complex phenomena, one of which is marginal pinching at the film's foot—a crucial yet inadequately understood aspect due to experimental limitations. Considering its profound impact on liquid film drainage and lifetime, we utilize high-resolution numerical simulations incorporating a weakly compressible scheme and adaptive mesh refinement to dissect the marginal pinching dynamics with unprecedented detail. Our approach elucidates the pinching dynamics and tracks the evolution of film thickness during the critical late-stage drainage process. By leveraging detailed geometric data from pinched regions, we significantly refine existing drainage models, enhancing their predictive accuracy regarding rupture thickness and film lifetime across various viscosities, surfactant concentrations, and bubble sizes. This refined model demonstrates robust alignment with our simulation results. Furthermore, we establish a quantifiable relationship between the prefactor governing reverse flow induced by the Marangoni effect and surfactant concentration. The methodologies and findings of this study provide foundational knowledge that paves the way for optimized industrial processes and an enhanced understanding of natural phenomena.

DOI： 10.1063/5.0235150

researchmap

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

DOI： 10.1016/j.jcp.2024.113292

researchmap

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

Haines jumps and meniscus reconfiguration are commonly observed in immiscible two-phase flows through porous media. The experiment found that the fluid-fluid interface moves rapidly during these two phenomena. This leads to hydraulic change and entropy production events, which have important implications for relative permeability of two-phase flow in porous media. However, owing to the energy-transfer process remains to be elucidated, their impact on transport of heat and mass in porous media has not yet been well predicted. By analyzing the exchange process between surface energy, kinetic energy, and viscous dissipation, this study provides a better understanding of meniscus dynamics and whether these effects are significant on a larger scale. In this study, Direct Numerical Simulation (DNS) combined with level-set method was used to simulate the process of Haines jumps and meniscus reconfiguration. A governing equation combined with density-scale continuum surface force (CSF) model was derived for energy conservation, and simulations were performed considering high-density and high-viscosity ratios. The rate of surface-energy change was evaluated using CSF model. A capillary tube was designed to verify the rate of surface-energy change in simulation. Comparison of simulation and theoretical data indicates that the density-scale CSF model can satisfactorily evaluate the rate of change in surface energy. Besides, to compare energy transfer process in Haines jumps and meniscus reconfiguration, two simulation cases were designed: (I) Haines jump and meniscus reconfiguration occur individually, and (II) Haines jump and meniscus reconfiguration occur simultaneously. We investigated the energy exchange process and compared the effects of the inertial force and viscous dissipation between cases I and II. Based on the simulation results, the interface releases surface free energy to the system during the Haines jumps and meniscus reconfiguration. In the case of only Haines jumps, a portion of the released surface energy is absorbed by another fluid-fluid interface, thus reducing the surface-energy release rate in the system. This energy can be dissipated promptly, resulting in a rapid viscous-dissipation growth; but the kinetic-energy change is negligible. However, when meniscus reconfiguration occurs, released surface energy is difficult to reabsorb and cannot dissipate quickly. Thus, a significant amount of released energy is converted into kinetic energy, which is subsequently damped by the fluid. This significant kinetic-energy oscillation generates a large inertial force and spreads over a wider sphere of influence in the porous media. Energy transfer process shows a more detailed fundamental understanding of mass transfer in porous media, such as flow perturbations, entropy production and permeability in porous media, and it helps a rigorous process for upscaling models used for thermodynamic, heat and mass transfer in porous media.

DOI： 10.1016/j.ijheatmasstransfer.2024.125749

researchmap

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

Investigating Rayleigh–Bénard convection, its onset, and its effect on mass transfer are crucial for understanding the effective dissolution processes between carbon-dioxide (CO2)-saturated and CO2-free brines in CO2 geological storage. However, most previous investigations on the onset of convection and mass transfer properties were conducted by 2D and 3D numerical simulations. Only a few 3D experimental studies of Rayleigh–Bénard convection were reported. In this study, we proposed a non-destructive experimental approach to visualize the three-dimensional (3D) fingering structure of Rayleigh–Bénard convection using a micro-focused X-ray computed tomography (X-ray CT) scanner. The experimental conditions covered the Rayleigh number (Ra) range between 3400 and 23000, which also includes the Ra value of the actual condition of gas reservoir field. Results show that fingers and coalescence dynamics significantly developed upon increasing Ra. At higher Ra, the finger concentration experienced greater dilution due to the combined effects of unsteady finger motion and enhanced dispersion strength. The maximum finger extension velocity is linearly proportional to Ra. The convective mass flux (Sh) correlated with Ra by a power-law relationship. The findings in this study provide a better understanding of the influence of Ra on mass transport properties and can be used as a reference for developing the convective mass transfer model of Rayleigh–Bénard convection in 3D porous media.

DOI： 10.1016/j.advwatres.2024.104666

researchmap

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

DOI： 10.1029/2023WR034849

researchmap

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

DOI： 10.1063/5.0139528

CiNii Research

researchmap

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

CiNii Research

J-GLOBAL

researchmap

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

DOI： 10.3390/app13031955

researchmap

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

Drying process of porous medium is of interest in engineering applications. However, pore-scale studies of wettability effects on the drying process of three-dimensional (3D) porous media remain limited. X-ray microtomography was used to investigate the wettability effects on the drying process of 3D porous medium from a pore-scale perspective. Results show higher wettability increases capillary pressure, which enhances the rearrangement of liquid from large pores to small pores. We also found that liquid clusters are unlikely to disappear for a long time confirming that higher wettability supports the flow through liquid film. This study also provides new insight that capillary flow by liquid film is associated with branch clusters. These clusters maintain the gas–liquid interfacial area that acts as the evaporation surface. As a result, drying is intensive in a water-wet porous medium with a longer constant rate period in which drying is controlled by external diffusion. In a neutral-wet porous medium, an early drop in drying rate occurs. Thus, drying is mainly a falling rate period limited by vapor diffusion through the newly emerging dried region within the medium.

DOI： 10.1016/j.icheatmasstransfer.2022.106527

researchmap

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

X-ray micro-tomography was used to comprehensively study the pore-scale drying process and characterize the liquid film region and the invasion front. Our finding reveals that the decrease in saturation is not uniform across the pore size. The capillary rearrangement, which refers to the restructuring of the liquid phase due to the capillary pressure difference, drives the liquid to flow from large pores to small pores. As a result, the decrease in liquid saturation occurs preferentially within the larger pores. A sharp decrease in liquid saturation occurred at the invasion front, where saturation ranged from 1 to around 0.15. The obtained invasion width that characterized the extent of roughness and disorder of the invasion front obeyed the scaling law of the percolation model. Therefore, this region is significantly limited by the interaction between capillary and gravity. However, with more multi-pore connected clusters observed within the porous medium, the critical exponent of the cumulative cluster size distribution obtained from the cluster region deviates from the universal power-law behavior. This result indicates the presence of liquid film within the cluster region, where the viscous force becomes comparable to the capillary force. We also found that the total gas-liquid-specific interfacial area has no effect on the drying rate. This result provides new evidence that the liquid film region is vapor saturated.

DOI： 10.1016/j.ijheatmasstransfer.2022.123299

researchmap

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

An in-situ dual-surfactant formulation is developed by using a novel cationic surfactant adsorbing the in-situ acid components in crude oil. As a result, various oil deformation patterns are first revealed and applied to control the viscous instabilities for oil recovery application. To this end, four dynamic oil displacement tests are conducted via camera snapshotting or X-ray microtomography in two and three dimensions, without and with porous media, in single to multipore-throat structures as well as through pore-scale and macroscale visualizations. As suggested in the results, displacing the oil phase using a 20 mmol/L cationic surfactant is the most effective method to suppress interfacial instability. The preferential flow path is frequently diverted by the repeated change of oil shape at a low flow rate. By contrast, emulsification plays a key role when the flow rate is high. It depends on which is dominant between the oil deformation and shearing velocity. From a deep insight into the pore spaces, oil deformation patterns are classified into oil spreading, shrinking, splitting, film detachment, and bridge breakup. The most significant oil deformation in the tortuous pore spaces occurs after 10 min of aging. Based on this time, we propose a new injection–alternate–stop scenario to suppress the viscous fingering and maximize the oil recovery efficiency to 80% in a large-scale reservoir.

DOI： 10.1021/acs.iecr.2c01735

researchmap

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

Microemulsion flooding has recently gained attention as an effective chemical-based method for enhanced oil recovery (EOR). The high performance of this technique is attributed to the classical mechanism: plugging effect and divergence of preferential flow path. In this study, miscible behavior was identified as a novel mechanism by which oil recovery was controlled by dissolution process instead of capillary forces. To this end, a new type of solvent-based microemulsion phase was prepared by mixing solvent, water, and surfactant as a flooding fluid. A series of experiments with direct water flooding, microemulsion flooding, and water-alternate-emulsion flooding for EOR were carried out for strong water-wet and oil-wet artificial rock sample. X-ray microtomography was used to visualize the oil configurations in three-dimensional pore spaces. An advanced workflow image processing method was developed for the first time to capture the miscible behavior in a pore-scale view. As suggested in the results, microemulsion flooding has a miscible displacement front where the oil phase is gradually solubilized. Therefore, capillary forces can be eliminated because the emulsion-oil interface is invisible but has a compositional gradient, leading to higher oil recovery. Compared with the immiscible waterflooding process, the capillary forces remained dominant, which hinders additional oil recovery. Based on EOR scenarios, we found that emulsion flooding could effectively solubilize larger well-connected oil clusters instead of small separate ganglia or singlets. Therefore, considering the miscible behavior, emulsion flooding can be applied to secondary oil recovery regardless of the reservoirs’ wettability, or to tertiary oil recovery in weakly/strongly oil-wet reservoirs with a well-connected oil phase after waterflooding.

DOI： 10.1016/j.petrol.2022.110788

researchmap

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

When viscosity and density contrast exist in the vertical miscible displacement in porous media between two fluids, the interplay between the viscous force and gravity determines the interface stability. Two stability criteria are derived to determine the interface stability. Hill's and Dumore's stability criteria are used to determine the interface stability of the sharp and diffused interface, respectively. In this study, we visualized the crossover between unstable displacement and stable displacement for a vertical displacement in porous media using microfocused x-ray computed tomography. The experiments were divided into four possible configurations: (1) unconditionally stable (gravitationally stable-viscously stable), (2) unconditionally unstable (gravitationally unstable-viscously unstable), (3) conditionally stable (gravitationally stable-viscously unstable), and (4) conditionally stable (gravitationally unstable, viscously stable). The structure of the displacement interface was visualized for the critical velocity ratio (V/Vc) in the range of 0.5–11.9. In the conditionally stable configurations 3 and 4, a crossover between stable and unstable displacements was observed. We found that Dumore's stability criterion is more appropriate for predicting interface stability than Hill's stability criterion. Viscous fingering occurs in configuration 3 when V/Vc is higher than Dumore's critical velocity, whereas gravity fingering occurs in configuration 4 when V/Vc is lower than Dumore's critical velocity. Similar events in two-dimensional experiments, such as tip-splitting, shielding, and coalescence, were also observed three-dimensionally. The significant changes in the mixing length and sweep efficiency signify the crossover between the stable and unstable displacements.

DOI： 10.1063/5.0090387

researchmap

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

Crude oil often contains long-chain fatty acid components. The fatty acid usually exists as a dissociative anionic surfactant in the oil phase. This study introduced a novel cationic surfactant to absorb the in situ anionic surfactant. Thus, a dual surfactant system was formed, and surfactant aggregates were produced near the water–oil interface. Due to the viscoelastic and fragile aggregates, spontaneous oil deformation and splitting were discovered as a novel mechanism for chemical enhanced oil recovery (EOR). First, we used a water phase doped with stearyltrimethylammonium chloride to displace the heavy oil with palmitic acid in the Hele-Shaw cell. The fingering patterns were observed in the water flooding (WF) and chemical flooding (CF) systems under different flow rates. After that, the experiments were transferred to the wettability-altered micromodel to simulate real oil recovery applications. The Hele-Shaw cell experiments show that viscous fingering was intensive at a low flow rate whereas it was significantly suppressed at a high flow rate in the CF system. The micromodel experiments show that oil recovery is higher in the CF system than in the WF system. The spontaneous cocurrent and counter-current oil blebbing was observed in the dynamic CF process. As the blebbing developed to a certain degree, the oil phase split into tiny droplets that could easily pass through the pore spaces. Additionally, the wettability was also altered from oil-wet to water-wet. As a result, oil recovery was improved significantly. The residual oil was classified into cluster, ganglia, and singlet based on the circularity and Euler number. We found the large clusters dominated in the WF system, whereas the small size of ganglia and singlets were the main forms in the CF system. The new chemical used can be applied for large-scale industrial EOR fields due to its good performance and low cost.

DOI： 10.1021/acs.energyfuels.2c00772

researchmap

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

Hydrodynamic fingering induced by gel formation shares common features with growing biofilms, bacterial colonies, and the instability of a confined chemical garden. Fluid displacement with gel formation is also essential in various engineering applications, including CO2 leakage remediation from storage reservoirs and enhanced oil recovery. We conducted Hele-Shaw cell displacement experiments for a miscible fluid system using skim milk and aqueous citric acid solution. This study aimed to investigate the effects of gel film formation on the fingering instability of a miscible fluid system and develop a mathematical model of the sequential growth of gel film formation at the fingertip. We found that the gel film formation thickens with time, resulting in instability at the interface. A distinctive fingering pattern, resembling tentacles, appears where miscibility is suppressed, and the growth of the finger is localized at the fingertip. The finger width remains constant with increasing flow rate, whereas the number of fingers increases linearly before the fingers merge. The gap width significantly limits the finger width. Finally, a mathematical model of sequential film thickness growth for a bubble-like fingertip structure was developed. This model is based upon the interplay between the diffusion of citric acid through the existing gel film formation and elongation of the fingertip. The model provides an understanding of the fundamental mechanism of the growth of the bubble-like fingertip.

DOI： 10.3390/app12105043

researchmap

Language：English   Publishing type：Research paper (scientific journal)   Publisher：Academic Press Inc.  

DOI： 10.1016/j.jcp.2021.110605

Scopus

researchmap

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

researchmap

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

researchmap

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

researchmap

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

researchmap

Publisher：The Japanese Society for Multiphase Flow  

DOI： 10.3811/jjmf.2021.004

J-GLOBAL

researchmap

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

DOI： 10.1016/j.jiec.2021.03.046

researchmap

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

DOI： 10.1063/5.0044756

Scopus

researchmap

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

researchmap

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

researchmap

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

researchmap

Language：Japanese   Publisher：Japanese Society of Tribologists  

Simulation of gas-liquid two-phase flows include multi-scale and multi-physics problems. We have developed a weakly compressible computational scheme for incompressible flows with a purely explicit time integration and successfully implemented on GPU (Graphics Processing Unit) providing very high computational performance. AMR (Adaptive Mesh Refinement) method is available for the scheme and high-resolution meshes are assigned dynamically near the region of the gas-liquid interface. The AMR method greatly enhances the computational efficiency and makes it possible to compute the dynamics of thin liquid film such as a soap bubble which is a typical multi-scale problem. At the surface, viscoelasticity becomes more important and transportation of surfactant should be taken into consideration. It has been confirmed that the Marangoni effect stabilizes and avoids breaking-up a thin film covering a bubble on the water surface against the gravity in the simulation.

DOI： 10.18914/tribologist.65.11_678

CiNii Books

CiNii Research

J-GLOBAL

researchmap

Other Link： https://ndlsearch.ndl.go.jp/books/R000000004-I030778045

Language：Japanese   Publishing type：Doctoral thesis  

researchmap

Language：Japanese  

J-GLOBAL

researchmap

Authorship：Lead author   Language：English   Publishing type：Research paper (scientific journal)   Publisher：Elsevier BV  

A weakly compressible scheme for gas–liquid two-phase flows with low Mach number is proposed for a fully-explicit time integration. The scheme is based on the fractional-step method with the Euler equation solver and time integration for viscous, surface tension and gravity terms. For the efficient solution of the Euler equation, we introduce the directional-splitting method and the method of characteristics which is applicable to semi-Lagrangian methods. A reduction of the speed of sound is also taken into consideration to enlarge the time step intervals. The primitive flow variables are time-integrated on the node-center configuration. To capture the gas–liquid interface, the conservative phase-field equation is solved by using the multi-moment method. It is computed as a conservative form for the volume-averaged value defined at the cell center, in which the flow velocities are given on the node. We discard the exact momentum conservation, however the total mass of gas and liquid is conserved by introducing a phase-field variable. Several benchmark problems are examined and the results of the explicit scheme are in good agreement with the experimental and computational references.

DOI： 10.1016/j.jcp.2018.10.019

researchmap

Authorship：Lead author   Language：Japanese   Publisher：THE JAPANESE SOCIETY FOR MULTIPHASE FLOW  

<p>A weakly compressible scheme for low-Mach number gas-liquid two-phase flows have been developed for a full-explicit time integration to avoid solving Poisson equation for large-scale two-phase flow simulations. To describe the gas-liquid interface the conservative Allen-Cahn equation, which is one of the phase-field models, is introduced and combined with the continuum equation. The accuracy of numerical results of two-phase flow strongly depends on the mesh resolution near the interface. The AMR(Adaptive Mesh Refinement) method greatly reduces the computational cost to assign high-resolution mesh to the region around the moving interface. We have developed a GPU-code of the tree-based AMR. We choose the node-centered configuration and the interpolation can be closed within a leaf in refining the leaf. We successfully have carried out two-phase flow simulations including very thin liquid film efficiently.</p>

DOI： 10.3811/jjmf.2019.008

CiNii Books

CiNii Research

J-GLOBAL

researchmap

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

researchmap

Authorship：Lead author   Language：Japanese   Publisher：JAPAN SOCIETY FOR COMPUTATIONAL ENGINEERING AND SCIENCE  

A phase-field model describes the gas-liquid interface of two-phase flows. We solve the conservative Allen-Cahn equation combined with the continuum equation on a two-dimensional computational domain with a given velocity field. The accuracy of numerical result strongly depends on the mesh resolution around the interface. The AMR (Adaptive Mesh Refinement) method greatly reduces the computational cost, since it is possible to assign high-resolution mesh to the region around the moving interface. We have developed a code to solve the equation in a manner of the tree-based AMR with multi-moment methods, the conservative IDO and CIP-CSL schemes. To reduce the implementation difficulties of AMR method, we introduce the fractional step method and the directional splitting method. In a benchmark test of the single vortex problem, the AMR computation with 5-level refinement for the interface achieves 9.26-times speed up and 1/12.3 mesh reduction to compare with the computation on a uniform mesh.

DOI： 10.11421/jsces.2018.20180005

CiNii Research

J-GLOBAL

researchmap

Language：Japanese  

CiNii Research

J-GLOBAL

researchmap

Other Link： https://ndlsearch.ndl.go.jp/books/R000000004-I031500280

J-GLOBAL

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Language：Japanese  

J-GLOBAL

researchmap

Language：Japanese  

CiNii Research

J-GLOBAL

researchmap

Language：Japanese  

identifier:oai:t2r2.star.titech.ac.jp:50588402

CiNii Research

researchmap

Language：Japanese  

CiNii Books

CiNii Research

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Language：Japanese  

CiNii Research

J-GLOBAL

researchmap

Language：Japanese  

CiNii Books

CiNii Research

J-GLOBAL

researchmap

Language：Japanese   Publisher：The Japan Society of Mechanical Engineers  

<p>A fully explicit gas-liquid two-phase flow solver without solving Poisson equation based on weakly compressible assumption with interface-adapted AMR method is proposed to achieve high resolution simulation. Solving independent hyperbolic pressure evolution equation derived from isothermal Navier-Stokes equation allows us to achieve a fully explicit time integration. We have developed a GPU-code of the tree-based AMR method, which can greatly reduce the computational cost to assign high-resolution mesh to the region around the moving interfaces. The memory pool data structure having extra allocated space is prepared to prevent frequent memory allocation and deallocation to enhance the computational performance. The results of benchmarks have good agreements with that of incompressible solver and experiments. The soap bubble simulation using the finest mesh width (= 0.07825mm) on single GPU are demonstrated and the thin liquid film can be reproduced by achieving high-resolution two-phase flow simulation with explicit scheme and AMR method.</p>

DOI： 10.1299/jsmecmd.2019.32.173

CiNii Research

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Language：Japanese  

CiNii Research

J-GLOBAL

researchmap

Language：Japanese   Publisher：THE JAPANESE SOCIETY FOR MULTIPHASE FLOW  

<p>Incompressible flows have much lower flow speeds to compare with the sound speed and called low Mach number flows. So far, we have simulated incompressible flows by using semi-implicit solvers, in which we iteratively compute the linear equation derived from the pressure Poisson equation. In the case of incompressible gas-liquid two-phase flows, the sparse matrix would be very stiff with non-zero elements changing largely from the gas density to the liquid density at the interface. Even if we use a sophisticated Krylov sub-space iteration solver with the multi-grid preconditioner, the convergence becomes poor for large-scale problems. We have developed a weakly compressible gas-liquid two-phase solver based on fully-explicit time integration of compressible Navier-Stokes equation. For several benchmark problems, the results of the weakly compressible solver are in good agreement with those of the semi-implicit incompressible solver. Since the weakly compressible solver uses a local stencil access, it is possible to implement AMR (Adaptive Mesh Refinement) method with less programming cost. Fine meshes are assigned near the gas-liquid interface and we efficiently compute a liquid film generated by a spoon and a growing soap bubble with the AMR two-phase flow solver.</p>

DOI： 10.3811/jjmf.2018.T011

CiNii Books

CiNii Research

J-GLOBAL

researchmap

Language：Japanese   Publisher：The Japan Society of Mechanical Engineers  

<p>A weakly compressible scheme for low-Mach number gas-liquid two-phase flows have been developed for a fullexplicit time integration to avoid solving Poisson equation. We introduce the fractional splitting method and directional splitting method to solve Euler equation which an efficient semi-Lagrangian scheme is applicable to. To describe the gas-liquid interface the conservative Allen-Cahn equation, which is one of the phase-field models, combined with the continuum equation is introduces. The accuracy of numerical results of two-phase flow strongly depends on the mesh resolution near the interface. The AMR (Adaptive Mesh Refinement) method greatly reduces the computational cost, since it is possible to assign high-resolution mesh to the region around the moving interface. We have developed a code to solve the equation in a manner of the tree-based AMR. Fractional step method and the directional splitting method helps to reduce the implementation difficulties of AMR method. We successfully carried out some violent two-phase flow simulation and flow including very thin liquid film by using the method which is combined weakly compressible scheme with AMR method.</p>

DOI： 10.1299/jsmecmd.2018.31.222

CiNii Research

J-GLOBAL

researchmap

Language：Japanese   Publisher：The Japan Society of Mechanical Engineers  

<p>In distributed implementations of stencil AMR applications, the inter-node communication time often represents a major performance bottleneck. Thus minimizing communication is an objective as important as maintaining a good load balance. We propose a new domain partitioning method for block-based AMR applications based on the multi-phase-field (MPF) model. The MPF model for polycrystalline growth minimizes the interfacial energy and forms a convex shape for each crystal grain. In our method, each phase of the MPF model represents a computational sub-domain. We show that the MPF partitioning can successfully keep the load balance and reduce the communication cost.</p>

DOI： 10.1299/jsmecmd.2018.31.213

J-GLOBAL

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Language：Japanese  

CiNii Research

J-GLOBAL

researchmap

Language：Japanese   Publisher：The Japan Society of Mechanical Engineers  

DOI： 10.1299/jsmecmd.2017.30.034

CiNii Research

J-GLOBAL

researchmap

Language：Japanese   Publisher：The Japan Society of Mechanical Engineers  

DOI： 10.1299/jsmecmd.2017.30.070

J-GLOBAL

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Language：Japanese   Publisher：The Japan Society of Mechanical Engineers  

DOI： 10.1299/jsmecmd.2016.29.089

J-GLOBAL

researchmap

J-GLOBAL

researchmap

Event date： 2018

Language：Japanese  

<p>A weakly compressible scheme for low-Mach number gas-liquid two-phase flows have been developed for a fullexplicit time integration to avoid solving Poisson equation. We introduce the fractional splitting method and directional splitting method to solve Euler equation which an efficient semi-Lagrangian scheme is applicable to. To describe the gas-liquid interface the conservative Allen-Cahn equation, which is one of the phase-field models, combined with the continuum equation is introduces. The accuracy of numerical results of two-phase flow strongly depends on the mesh resolution near the interface. The AMR (Adaptive Mesh Refinement) method greatly reduces the computational cost, since it is possible to assign high-resolution mesh to the region around the moving interface. We have developed a code to solve the equation in a manner of the tree-based AMR. Fractional step method and the directional splitting method helps to reduce the implementation difficulties of AMR method. We successfully carried out some violent two-phase flow simulation and flow including very thin liquid film by using the method which is combined weakly compressible scheme with AMR method.</p>

researchmap

Event date： 2018

Language：Japanese  

A phase-field model describes the gas-liquid interface of two-phase flows. We solve the conservative Allen-Cahn equation combined with the continuum equation on a two-dimensional computational domain with a given velocity field. The accuracy of numerical result strongly depends on the mesh resolution around the interface. The AMR (Adaptive Mesh Refinement) method greatly reduces the computational cost, since it is possible to assign high-resolution mesh to the region around the moving interface. We have developed a code to solve the equation in a manner of the tree-based AMR with multi-moment methods, the conservative IDO and CIP-CSL schemes. To reduce the implementation difficulties of AMR method, we introduce the fractional step method and the directional splitting method. In a benchmark test of the single vortex problem, the AMR computation with 5-level refinement for the interface achieves 9.26-times speed up and 1/12.3 mesh reduction to compare with the computation on a uniform mesh.

researchmap

Event date： 2018

Language：Japanese  

<p>In distributed implementations of stencil AMR applications, the inter-node communication time often represents a major performance bottleneck. Thus minimizing communication is an objective as important as maintaining a good load balance. We propose a new domain partitioning method for block-based AMR applications based on the multi-phase-field (MPF) model. The MPF model for polycrystalline growth minimizes the interfacial energy and forms a convex shape for each crystal grain. In our method, each phase of the MPF model represents a computational sub-domain. We show that the MPF partitioning can successfully keep the load balance and reduce the communication cost.</p>

researchmap

Event date： 2017.5

Language：Japanese  

researchmap

Event date： 2017

Language：Japanese  

researchmap

Event date： 2017

Language：Japanese  

researchmap

Event date： 2016

Language：Japanese  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

Presentation type：Oral presentation (general)  

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap

researchmap