Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid

Smoothed particle hydrodynamics method applied to oral region: A narrative review

Traditionally, simulation studies in dentistry have relied on the finite element method (FEM). However, the smoothed particle hydrodynamics (SPH) method, which represents objects as particle collections without the use of meshes, has gained recent attention. Despite its application in dentistry, there is currently a lack of comprehensive literature summarizing the specific applications of the SPH method in the oral region. This review aims to provide a summary of studies that have utilized the SPH method in dentistry, focusing on its applications in analyzing large deformations, such as dental ceramic collisions, soft material analysis (e.g., denture adhesive), and virtual training simulations for dental treatments. By combining the advantages of the SPH and FEM methods, more accurate simulations can be achieved, and further applications of the SPH method in dentistry are anticipated.

Moving Particle Semi-Implicit and Finite Element Method Coupled Analysis for Brain Shift Estimation

Neuronavigation is a computer-assisted technique for presenting three-dimensional images of a patient’s brain to facilitate immediate and precise lesion localization by surgeons. Neuronavigation systems use preoperative medical images of patients. In neurosurgery, when the dura mater and arachnoid membrane are incised and the cerebrospinal fluid (CSF) drains out, the brain loses the CSF buoyancy and deforms in the direction of gravity, which is referred to as brain shift. This brain shift yields inaccurate neuronavigation. To reduce this inaccuracy, an intraoperative brain shift should be estimated. This paper proposes a dynamic simulation method for brain-shift estimation combining the moving-particle semi-implicit (MPS) method and the finite element method (FEM). The CSF was modeled using fluid particles, whereas the brain parenchyma was modeled using finite elements (FEs). Node particles were attached to the surface nodes of the brain parenchyma in the FE model. The interaction between the CSF and brain parenchyma was simulated using the repulsive force between the fluid particles and node particles. Validation experiments were performed using a gelatin block. The gelatin block was dipped into silicone oil, which was then gradually removed; the block deformation owing to the buoyancy loss was measured. The experimental deformation data were compared with the results of the MPS-FEM coupled analysis. The mean absolute error (MAE) between the simulated deformation and the average across the four experiments was 0.26 mm, while the mean absolute percentage error (MAPE) was 27.7%. Brain-shift simulations were performed using the MPS-FEM coupled analysis, and the computational cost was evaluated.

Development of an MPS Code for Corium Behavior Analysis: 3D Alloy Melting

The moving particle semi-implicit (MPS) method as a Lagrangian method is attracting increasing attention in severe accident analysis. In this paper, we developed an MPS code for the corium behavior analysis with several additional models added: an improved heat transfer model to improve the calculation between different materials, an enthalpy-based viscosity model to realize a smooth transition of viscosity at the solid-liquid interface, and a surface tension model for better simulation of surface shape. Validation of the developed simulation approach is carried out on a classical water column collapse example. The development of the heat transfer model is validated by the example of a one-dimensional semi-infinite plate. A comprehensive example of the melting of “Wood’s alloy” is carried out to verify the capacity of MPS method in the simulation of melting and expansion. The simulation results are in good agreement with the experimental results, which indicates that MPS method promises well in the field of severe accidents.

Splitting of a three-dimensional liquid plug at an airway bifurcation

Employing the moving particles' semi-implicit (MPS) method, this study presents a numerical framework for solving the Navier-Stokes equations for the propagation and the split of a liquid plug through a three-dimensional air-filled bifurcating tube, where the inner surface is coated by a thin fluid film, and surface tension acts on the air-liquid interface. The detailed derivation of a modified MPS method to handle the air-liquid interface of liquid plugs is presented. When the front air-liquid interface of the plug splits at the bifurcation, the interface deforms quickly and causes large wall shear stress. We observe that the presence of a transverse gravitational force causes asymmetries in plug splitting, which becomes more pronounced as the capillary number decreases or the Bond number increases. We also observe that there exists a critical capillary number below which the plug does not split into two daughter tubes but propagates into the lower daughter tube only. In order to deliver the plug into the upper daughter tube, the driving pressure to push the plug is required to overcome the hydrostatic pressure due to gravity. These tendencies agree with our previous experimental and theoretical studies.

Review of the State-of-Art of MPS Method in Ocean Engineering

When dealing with the complex deformation of free surface such as wave breaking, traditional mesh-based Computational Fluid Dynamics (CFD) methods often face problems arising alongside grid distortion and re-meshing. Therefore, the meshless method became robust for treating large displaced free surface and other boundaries caused by moving structures. The particle method, which is an important branch of meshless method, is mainly divided into the Smoothed Particle Hydrodynamics (SPH) and Moving Particle Semi-implicit (MPS) methods. Different from the SPH method, which involves continuity and treat density as a variable when building kernel functions, the kernel function in the MPS method is a weight function which treats density as a constant, and the spatial derivatives are discretized by establishing the gradient operator and Laplace operator separately. In other words, the first- or second-order continuity of the kernel functions in the MPS method is not a necessity as in SPH, though it might be desirable. At present, the MPS method has been successfully applied to various violent-free surface flow problems in ocean engineering and diverse applications have been comprehensively demonstrated in a number of review papers. This work will focus on algorithm developments of the MPS method and to provide all perspectives in terms of numerical algorithms along with their pros and cons.

Numerical Computation of Sloshing-Induced Force in Complex Ship Tanks under the Excitation of Ship Rolling Motion Based on the MPS Method

Liquid sloshing in ship tanks would have a direct effect on ship dynamic stability, and thus is of great importance for navigation safety. To calculate the sloshing-induced force on real ship tanks, this paper presents an approach to numerically simulate the liquid-sloshing in complex tanks using the moving particle semi-implicit (MPS) method. The sloshing-induced force is numerically calculated and used to investigate the effect of different excitation conditions in which the realistic ship motions under different loading conditions have been taken into account. Simulation results show that the maximum sloshing-induced force is much bigger than the corresponding static one. Meanwhile, both the rolling angle and period have significant effects on liquid sloshing.

A phantom and in vivo simulation of coronary flow to calculate fractional flow reserve using a mesh-free model

Moving particle semi-implicit (MPS) method is a mesh-free method to perform computational fluid dynamics (CFD). The purpose of this study was to calculate the simulated fractional flow reserve (sFFR) using a coronary stenosis model, and to validate the MPS-derived sFFR against invasive FFR using clinical coronary CT data. Coronary flow simulation included 21 stenosis models with stenosis ranging 30-70%. Patient coronary analysis was performed in 76 consecutive patients (100 vessels) who underwent coronary CT angiography and subsequent invasive FFR between November 2016 and March 2020. Accuracy of sFFR and CT angiography for diagnosis of invasive FFR ≤ 0.80 was compared. Quantitative morphological stenosis data of CT angiography were also obtained. Area stenosis showed a good correlation to sFFR (R2 = 0.996, p < 0.001) in coronary stenosis models. In the patient study, the mean FFR value was 0.82 ± 0.10, and 37 out of 100 vessels showed FFR ≤ 0.80. FFR and sFFR values showed a good correlation (R2 = 0.59, p < 0.001) with a slight underestimation of sFFR as compared with FFR (mean difference - 0.015 ± 0.096, p = 0.12). The sensitivity, specificity, positive predictive value, and negative predictive value of sFFR to predict FFR ≤ 0.80 was 86%, 89%, 82%, 92%, respectively. The accuracy to predict FFR ≤ 0.80 using sFFR was greater than using diameter stenosis and minimum lumen area (88% vs. 74%, p = 0.008). CFD using the MPS method showed feasible results validated against invasive FFR. The accuracy to predict significant stenosis was higher than morphological stenosis.

Weighted Moving Square-Based Solver for Unsteady Incompressible Laminar Flow Simulations

For computational fluid dynamics simulations, grid systems are generally used in the Eulerian frame for both structured and unstructured grids and solvers designed for the chosen grid systems. In contrast to the grid-based method, in which the connection information with neighboring grids must be maintained, gridless methods do not require an underlying connectivity in the form of control volumes or elements. Hence, gridless methods are feasible and robust for the problems with moving boundary and/or complicated boundary shapes. In this study, a Eulerian gridless solver is proposed, and its application for simulating two-dimensional unsteady viscous flows in low Reynolds number regions is investigated. The solver utilizes the weighted moving square method to obtain the spatial derivatives of the governing equations and solve the pressure Poisson equation iteratively. Simple remedies to satisfy the boundary conditions in the finite difference method are applied. The fractional step method with the second-order Adams–Bashforth method is used for time integration. Some benchmark problems were solved using the developed solver, and the results were compared with those of other experimental and computational methods. Good agreement with the results of other methods confirmed the validity of the proposed solver.

Application of moving particle semi-implicit (MPS) method on retro-oil fluid using three-dimensional vitreous cavity models from magnetic resonance imaging

Silicone oil (SO) is a safe and widely used intraocular tamponade agent for treating complicated vitreoretinal diseases, such as retinal detachments (RRDs) with inferior proliferative vitreoretinopathy (PVR). However, as the human vitreous cavity is irregularly shaped, it is difficult to predict the area of the inferior retina covered with SO and the retro-oil fluid currents in each patient. Here, we performed fluid simulation analysis using the moving particle semi-implicit method on the oil cover rates and absolute velocity gradient of retro-oil fluid to the retina using vitreous cavity models derived from magnetic resonance imaging of patients to determine the appropriate amount of SO and postoperative position to achieve a sufficient tamponade effect on the inferior retina. In all seven vitreous cavity models tested, the inferior quadrant of the retina was completely covered by SO in more positions and the absolute velocity gradient of the retro-oil fluid in contact with the retinal wall caused by eye and head movements was lower when the vitreous cavity was filled with 95% SO and 5% retro-oil fluid versus 80% SO and 20% retro-oil fluid. Taken together, these findings have clinical implications for the treatment of complicated RRDs with inferior PVR requiring SO tamponade.

Evaluation of Discharging Surplus Soils for Relative Stirred Deep Mixing Methods by MPS-CAE Analysis

Most of the ground in Japan is soft, leading to great damage in the event of liquefaction. Various ground-improvement measures are being taken to suppress such damage. However, it is difficult to carry out ground-improvement work while checking the internal conditions of the ground during the construction. Therefore, a visible and measurable evaluation of the performance of the ground-improvement work was conducted in this study. The authors performed a simulation analysis of the relative stirred deep mixing method (RS-DMM), a kind of ground-improvement method, using a computer-aided engineering (CAE) analysis based on particle-based methods (PBMs). In the RS-DMM, the “displacement reduction type (DRT)” suppresses displacement during construction. Both the DRT and the normal type (NT) were simulated, and a visible and measurable evaluation was performed on the internal conditions during each construction, the quality of the improved body, and the displacement reduction performance. As an example of these results, it was possible to visually evaluate the discharge of surplus soil by the spiral rod attached to the stirring wing of the DRT. In addition, the authors succeeded in quantitatively showing that more surplus soil was discharged when the stirring wing of the DRT was used than when the stirring wing of the NT was used.

Numerical Investigation of Shallow Liquid Sloshing in a Baffled Tank and the Associated Damping Effect by BM-MPS Method

Understanding the damping mechanism of baffles is helpful to make more reasonable use of them in suppressing liquid sloshing. In this study, the damping effect and mechanism of vertical baffles in shallow liquid sloshing under a rotational excitation are investigated by an improved particle method. By incorporation of a background mesh scheme and a modified pressure gradient model, the accuracy of impact pressure during sloshing is significantly enhanced. Combined with the advantages of the particle method, the present numerical method is a wonderful tool for the investigation of liquid sloshing issues. Through the analysis of impact pressure, the influences of baffle height and baffle position on the damping mechanism are discussed. The results show that the damping effect of vertical baffles increases with the increase of the elevation of baffle top and decreases with the increase of the elevation of the baffle bottom. Moreover, the resonance characteristics of sloshing are altered when static water is divided into two parts by the vertical baffle. The dominant damping mechanism of vertical baffles depends on the configurations.

Experimental and numerical studies of sediment removal in double bottom ballast tanks

The sediments inside the ballast tanks might compromise the effectiveness of the ballast water exchange procedure. In the present study, the sediment removal from double bottom ballast tanks due to flow-through ballast water exchange is investigated by simplified experimental and numerical models. Two flow rates and two sediment densities were considered to verify their influence on the process. Besides, the effectiveness of two tank configuration modifications to improve the sediments removal was investigated. The first is the introduction of a flow deflector inside the tank and the second is the water inflow from the bottom. A fully Lagrangian particle-based method is adopted to simulate the sediment removal process with complex tank geometry. A new boundary condition is proposed to model the mesh of the bottom inlet in the bottom inflow configuration. Despite the simplifications of the numerical model, the computational results are in good agreement with the experimental ones. The results show that the flow-through ballast water exchange method might not be effective to remove the sediments entrapped between the bottom stiffeners, and the proposed modifications, despite challenging implementation, improve the sediment removal.

Simulating Groups and the IntraGroup Medium: The Surprisingly Complex and Rich Middle Ground between Clusters and Galaxies

Galaxy groups are more than an intermediate scale between clusters and halos hosting individual galaxies, they are crucial laboratories capable of testing a range of astrophysics from how galaxies form and evolve to large scale structure (LSS) statistics for cosmology. Cosmological hydrodynamic simulations of groups on various scales offer an unparalleled testing ground for astrophysical theories. Widely used cosmological simulations with ∼(100 Mpc)3 volumes contain statistical samples of groups that provide important tests of galaxy evolution influenced by environmental processes. Larger volumes capable of reproducing LSS while following the redistribution of baryons by cooling and feedback are the essential tools necessary to constrain cosmological parameters. Higher resolution simulations can currently model satellite interactions, the processing of cool (T≈104−5 K) multi-phase gas, and non-thermal physics including turbulence, magnetic fields and cosmic ray transport. We review simulation results regarding the gas and stellar contents of groups, cooling flows and the relation to the central galaxy, the formation and processing of multi-phase gas, satellite interactions with the intragroup medium, and the impact of groups for cosmological parameter estimation. Cosmological simulations provide evolutionarily consistent predictions of these observationally difficult-to-define objects, and have untapped potential to accurately model their gaseous, stellar and dark matter distributions.

Numerical Modelling of Flow-Debris Interaction during Extreme Hydrodynamic Events with DualSPHysics-CHRONO

Floods can transport debris of a very wide range of dimensions, from cohesive sediments to large floating debris, such as trees and cars. The latter increases the risk associated with floods by, for example, obstructing the flow or damaging structures due to impact. The transport of this type of debris and their interaction with structures are often studied experimentally in the context of tsunamis and flash floods. Numerical studies on this problem are rare, therefore the present study focuses on the numerical modelling of the flow-debris interaction. This is achieved by simulating multiple laboratory experiments, available in the literature, of a single buoyant container transported by a dam-break flow in order to validate the chosen numerical approach. The numerical simulations are carried using the open source DualSPHysics model based on the smoothed particle hydrodynamics method coupled with the multi-physics engine CHRONO, which handles the container–bottom interactions. The trajectory, as well as the velocity of the centroid of the container, were tracked throughout the simulation and compared with the same quantities measured in the laboratory. The agreement between the model and the experiment results is quantitatively assessed using the normalised root mean squared error and it is shown that the model is accurate in reproducing the floating container trajectory and velocity.

Two-Stage Water Jet Landing Point Prediction Model for Intelligent Water Shooting Robot

In this paper, an intelligent water shooting robot system for situations of carrier shake and target movement is designed, which uses a 2 DOF (degree of freedom) robot as an actuator, a photoelectric camera to detect and track the desired target, and a gyroscope to keep the robot’s body stable when it is mounted on the motion carriers. Particularly, for the accurate shooting of the designed system, an online tuning model of the water jet landing point based on the back-propagation algorithm was proposed. The model has two stages. In the first stage, the polyfit function of Matlab is used to fit a model that satisfies the law of jet motion in ideal conditions without interference. In the second stage, the model uses the back-propagation algorithm to update the parameters online according to the visual feedback of the landing point position. The model established by this method can dynamically eliminate the interference of external factors and realize precise on-target shooting. The simulation results show that the model can dynamically adjust the parameters according to the state relationship between the landing point and the desired target, which keeps the predicted pitch angle error within 0.1°. In the test on the actual platform, when the landing point is 0.5 m away from the position of the desired target, the model only needs 0.3 s to adjust the water jet to hit the target. Compared to the state-of-the-art method, GA-BP (genetic algorithm-back-propagation), the proposed method’s predicted pitch angle error is within 0.1 degree with 1/4 model parameters, while costing 1/7 forward propagation time and 1/200 back-propagation calculation time.

Three-Dimensional Fluid Dynamical Features of Coronary Plaque Rupture Provoking Acute Coronary Syndrome

AIM

Coronary plaque rupture is the main cause of acute coronary syndrome (ACS), but the role of blood flow features around plaque rupture for ACS is still unknown. The present study aimed to assess the relationship between the geometric configuration of ruptured plaque and ACS occurrence using computational fluid dynamics (CFD) by moving particle method in patients with coronary artery disease.

METHODS

In this study, 45 patients with coronary artery disease who underwent three-dimensional intravascular ultrasound (IVUS) and had a coronary ruptured plaque (24 plaques with provoked ACS, 21 without) were included. To compare the difference in blood flow profile around ruptured plaque between the patients with and without ACS, the IVUS images were analyzed via the novel CFD analysis.

RESULTS

There were no significant differences in localized flow profile around ruptured plaque between the two groups when the initial particle velocity was 10.0 cm/s corresponded to a higher coronary flow velocity at ventricular diastole. However, when it was 1.0 cm/s corresponded to lower coronary flow velocity at ventricular systole, particles with lower velocity (0 ≤ V ≤ 5 cm/s) were more prevalent around ACS-PR (p=0.035), whereas particles with higher velocity (10 ≤ V ≤ 20 cm/s) were more often detected in silent plaque ruptures (p=0.018).

CONCLUSIONS

Three-dimensional IVUS revealed that coronary plaque rupture was a complex one with a wide variety of its stereoscopic configuration, leading to various patterns of the local coronary flow profile. A novel CFD analysis suggested that the local flow was more stagnant around ACS-provoked ruptures than in silent ones.

Numerical and Experimental Investigation of Rectangular Liquid-Containing Structures under Seismic Excitation

Liquid-containing structures are used for various municipal and industrial applications. The functionality of these structures in seismic regions is crucial. The main purpose of this study is to investigate the behavior of liquid under seismic excitations using numerical modelling. For this purpose, experimental and numerical studies are conducted. In the experimental tests, a ground-supported rectangular tank is excited on a shaking table. The tests are videotaped from two directions and subsequently analyzed frame-by-frame. Four different orientations are tested to investigate the effect of bilateral excitation. In the numerical simulations, the same tank is modeled in OpenFOAM—a computational fluid dynamics program—and the same excitations are applied. The results from the numerical and the experimental studies are compared, and reliability of the numerical model is discussed. Furthermore, using the numerical model, the pressure on the roof of the tank is obtained at various locations and examined for different excitations.