Source author record

Yosuke Hasegawa

Yosuke Hasegawa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn physics.comp-ph Biological Physics Computational Engineering, Finance, and Science Machine Learning math.NA Numerical Analysis

Catalog footprint

What is connected

8works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2026arXiv

A volume penalization method for solving conjugate scalar transport with interfacial jump conditions

Conjugate scalar transport with interfacial jump conditions on complex interfacial geometries is common in thermal and chemical processes, while its accurate and efficient simulations are still quite challenging. In the present study, a novel treatment of a two-phase interface in the volume penalization method, a kind of immersed boundary method, for solving conjugate scalar transport with general interfacial boundary conditions is developed. We first propose an interfacial treatment for solving an advection-diffusion equation with a Neumann boundary condition, and then extend it to general conjugate scalar transport with both interfacial flux and scalar jumps. A one-dimensional diffusion problem is solved to verify the present scheme and demonstrate the advantage of the present scheme in improving accuracy and unifying the governing equations in the two phases with an additional source term representing the local jump condition of the interfacial scalar flux. Then, the present scheme is further applied to fluid-solid coupled scalar diffusion and advection-diffusion problems with the scalar and its flux jumps across the interface. The simulation results of the present scheme generally show good agreement with reference results obtained by body-fitted mesh simulations with average relative deviations less than 3.0%.

preprint2026arXiv

Volume penalization method for simulating flows around a rotating solid with multiple reference frame and sliding mesh

Despite the significant role of turbomachinery in fluid-based energy transfer, precise simulation of rotating solid objects with complex geometry is a challenging task. In the present study, the volume penalization method (VPM) is combined with multiple reference frame (MRF) and sliding mesh (SLM), respectively, so as to develop immersed-boundary approaches for simulating flows around a rotating solid. The level-set function is adopted to represent arbitrary geometries embedded in Cartesian grids. The VPM body-forcing terms in the momentum equation are proposed for MRF and SLM, respectively, so as to build unified governing equations for both fluid and solid regions. The flows around a rotating cuboid under various rotating speeds are simulated by the present schemes, namely, VPM with MRF, and VPM with SLM, and compared to corresponding simulations by the body-fitted method (BFM). The results suggest the relative deviations of predicted pressure drop and torque between the present VPM and BFM are around 5%, demonstrating the validity of the present VPM.

preprint2022arXiv

Multi-objective optimization of actuation waveform for high-precision drop-on-demand inkjet printing

Drop-on-demand (DOD) inkjet printing has been considered as one of promising technologies for the fabrication of advanced functional materials. For a DOD printer, high-precision dispensing techniques for achieving satellite-free smaller droplets, have long been desired for patterning thin-film structures. The present study considers the inlet velocity of a liquid chamber located upstream of a dispensing nozzle as a control variable and aims to optimize its waveform using a sample-efficient Bayesian optimization algorithm. Firstly, the droplet dispensing dynamics are numerically reproduced by using an open-source OpenFOAM solver, interFoam, and the results are passed on to another code based on pyFoam. Then, the parameters characterizing the actuation waveform driving a DOD printer are determined by the Bayesian optimization (BO) algorithm so as to maximize a prescribed multi-objective function expressed as the sum of two factors, i.e., the size of a primary droplet and the presence of satellite droplets. The results show that the present BO algorithm can successfully find high-precision dispensing waveforms within 150 simulations. Specifically, satellite droplets can be effectively eliminated and the droplet diameter can be significantly reduced to 24.9% of the nozzle diameter by applying the optimal waveform.

preprint2022arXiv

Optimization of a Moving Sensor Trajectory for Observing a Point Scalar Source in Turbulent Flow

We propose a strategy for optimizing a sensor trajectory in order to estimate the time dependence of a localized scalar source in turbulent channel flow. The approach leverages the view of the adjoint scalar field as the sensitivity of measurement to a possible source. A cost functional is constructed so that the optimal sensor trajectory maintains a high sensitivity and low temporal variation in the measured signal, for a given source location. This naturally leads to the adjoint-of-adjoint equation based on which the sensor trajectory is iteratively optimized. It is shown that the estimation performance based on the measurement obtained by a sensor moving along the optimal trajectory is drastically improved from that achieved with a stationary sensor. It is also shown that the ratio of the fluctuation and the mean of the sensitivity for a given sensor trajectory can be used as a diagnostic tool to evaluate the resultant performance. Based on this finding, we propose a new cost functional which only includes the ratio without any adjustable parameters, and demonstrate its effectiveness in predicting the time dependence of scalar release from the source.

preprint2021arXiv

Multiscale Parareal Algorithm for Long-Time Mesoscopic Simulations of Microvascular Blood Flow in Zebrafish

Various biological processes such as transport of oxygen and nutrients, thrombus formation, vascular angiogenesis and remodeling are related to cellular/subcellular level biological processes, where mesoscopic simulations resolving detailed cell dynamics provide a key to understanding and identifying the cellular basis of disease. To break this bottleneck and achieve a biologically meaningful timescale, we propose a multiscale parareal algorithm in which a continuum-based solver supervises a mesoscopic simulation in the time-domain. Using an iterative prediction-correction strategy, the parallel-in-time mesoscopic simulation supervised by its continuum-based counterpart can converge fast. The effectiveness of the proposed method is first verified in a time-dependent flow with a sinusoidal flowrate through a Y-shaped bifurcation channel. Physical quantities of interest including velocity, wall shear stress and flowrate are computed to compare against those of reference solutions, showing a less than 1% relative error on flowrate in the Newtonian flow and a less than 3\% relative error in the non-Newtonian blood flow. The proposed method is then applied to a large-scale mesoscopic simulation of microvessel blood flow in a zebrafish hindbrain for temporal acceleration. The time-dependent blood flow from heartbeats in this realistic vascular network of zebrafish hindbrain is simulated using dissipative particle dynamics as the mesoscopic model, which is supervised by a one-dimensional blood flow model (continuum-based model) in multiple temporal sub-domains. The computational analysis shows that the resulting microvessel blood flow converges to the reference solution after only two iterations. The proposed method is suitable for long-time mesoscopic simulations with complex fluids and geometries.

preprint2015arXiv

Does the choice of the forcing term affect flow statistics in DNS of turbulent channel flow?

We seek possible statistical consequences of the way a forcing term is added to the Navier--Stokes equations in the Direct Numerical Simulation (DNS) of incompressible channel flow. Simulations driven by constant flow rate, constant pressure gradient and constant power input are used to build large databases, and in particular to store the complete temporal trace of the wall-shear stress for later analysis. As these approaches correspond to different dynamical systems, it can in principle be envisaged that these differences are reflect by certain statistics of the turbulent flow field. The instantaneous realizations of the flow in the various simulations are obviously different, but, as expected, the usual one-point, one-time statistics do not show any appreciable difference. However, the PDF for the fluctuations of the streamwise component of wall friction reveals that the simulation with constant flow rate presents lower probabilities for extreme events of large positive friction. The low probability value of such events explains their negligible contribution to the commonly computed statistics; however, the very existence of a difference in the PDF demonstrates that the forcing term is not entirely uninfluential. Other statistics for wall-based quantities (the two components of friction and pressure) are examined; in particular spatio-temporal autocorrelations show small differences at large temporal separations, where unfortunately the residual statistical uncertainty is still of the same order of the observed difference. Hence we suggest that the specific choice of the forcing term does not produce important statistical consequences, unless one is interested in the strongest events of high wall friction, that are underestimated by a simulation run at constant flow rate.

preprint2014arXiv

Numerical simulation of turbulent duct flows with constant power input

The numerical simulation of a flow through a duct requires an externally specified forcing that makes the fluid flow against viscous friction. To this aim, it is customary to enforce a constant value for either the flow rate (CFR) or the pressure gradient (CPG). When comparing a laminar duct flow before and after a geometrical modification that induces a change of the viscous drag, both approaches (CFR and CPG) lead to a change of the power input across the comparison. Similarly, when carrying out the (DNS and LES) numerical simulation of unsteady turbulent flows, the power input is not constant over time. Carrying out a simulation at constant power input (CPI) is thus a further physically sound option, that becomes particularly appealing in the context of flow control, where a comparison between control-on and control-off conditions has to be made. We describe how to carry out a CPI simulation, and start with defining a new power-related Reynolds number, whose velocity scale is the bulk flow that can be attained with a given pumping power in the laminar regime. Under the CPI condition, we derive a relation that is equivalent to the Fukagata--Iwamoto--Kasagi relation valid for CFR (and to its extension valid for CPG), that presents the additional advantage of natively including the required control power. The implementation of the CPI approach is then exemplified in the standard case of a plane turbulent channel flow, and then further applied to a flow control case, where the spanwise-oscillating wall is used for skin friction drag reduction. For this low-Reynolds number flow, using 90% of the available power for the pumping system and the remaining 10% for the control system is found to be the optimum share that yields the largest increase of the flow rate above the reference case, where 100% of the power goes to the pump.

preprint2013arXiv

Prediction of turbulence control for arbitrary periodic spanwise wall movement

In order to generalize the well-known spanwise-oscillating-wall technique for drag reduction, non-sinusoidal oscillations of a solid wall are considered as a means to alter the skin-friction drag in a turbulent channel flow. A series of Direct Numerical Simulations is conducted to evaluate the control performance of nine different temporal waveforms, in addition to the usual sinusoid, systematically changing the wave amplitude and the period for each waveform. The turbulent average spanwise motion is found to coincide with the laminar Stokes solution that is constructed, for the generic waveform, through harmonic superposition. This allows us to define and compute, for each waveform, a new penetration depth of the Stokes layer which correlates with the amount of turbulent drag reduction, and eventually to predict both turbulent drag reduction and net energy saving rate for arbitrary waveforms. Among the waveforms considered, the maximum net energy saving rate is obtained by the sinusoidal wave at its optimal amplitude and period. However, the sinusoid is not the best waveform at every point in the parameter space. Our predictive tool offers simple guidelines to design waveforms that outperform the sinusoid for given (suboptimal) amplitude and period of oscillation. This is potentially interesting in view of applications, where physical limitations often preclude the actuator to reach its optimal operating conditions.