In simulating viscous incompressible SPH fluids, incompressibility and viscosity are typically solved in two separate stages. However, the interference between pressure and shear forces could cause the missing of beha...
详细信息
In simulating viscous incompressible SPH fluids, incompressibility and viscosity are typically solved in two separate stages. However, the interference between pressure and shear forces could cause the missing of behaviors that include preservation of sharp surface details and remarkable viscous behaviors such as buckling and rope coiling. To alleviate this problem, we introduce for the first time the semi-implicit method for pressure linked equations (simple) into SPH to solve incompressible fluids with a broad range viscosity. We propose to link incompressibility and viscosity solvers, and impose incompressibility and viscosity constraints iteratively to gradually remove the interference between pressure and shear forces. We will also discuss how to solve the particle deficiency problem for both incompressibility and viscosity solvers. Our method is stable at simulating incompressible fluids whose viscosity can range from zero to an extremely high value. Compared to state-of-the-art methods, our method not only produces realistic viscous behaviors, but is also better at preserving sharp surface details.
In 1998, Hu and Liu developed a strongly polynomial algorithm for solving the inverse arborescence problem that aims at minimally modifying a given cost-function on the edge-set of a digraph D so that an input spannin...
详细信息
In 1998, Hu and Liu developed a strongly polynomial algorithm for solving the inverse arborescence problem that aims at minimally modifying a given cost-function on the edge-set of a digraph D so that an input spanning arborescence of D becomes a cheapest one. In this note, we develop a conceptually simpler algorithm along with a new min-max formula for the minimum modification of the cost-function. The approach is based on a link to a min-max theorem and a simple (two-phase greedy) algorithm by the first author from 1979 concerning the primal optimization problem of finding a cheapest subgraph of a digraph that covers an intersecting family along with the corresponding dual optimization problem, as well. (C) 2021 The Author(s). Published by Elsevier B.V.
A comparative assessment between simple and its variant simple-C is conducted based on two-dimensional (2-D) incompressible flows, using a cell-centered finite-volume Li-formulation on a collocated grid. The simple-C ...
详细信息
A comparative assessment between simple and its variant simple-C is conducted based on two-dimensional (2-D) incompressible flows, using a cell-centered finite-volume Li-formulation on a collocated grid. The simple-C (simple Compressibility) scheme additionally combines the concept of artificial compressibility (AC) with the pressure Poisson equation, provoking the diagonal dominance of influence coefficients. An improved nonlinear momentum interpolation scheme is employed at the cell face in discretizing the continuity equation to suppress pressure oscillations. The pseudo-time step Delta t(i) remains the same in both schemes to conserve an analogous scaling with momentum and scalar nodal influence coefficients. Numerical experiments in reference to buoyancy-driven cavity flow dictate that both contrivances execute a residual smoothing enhancement, facilitating an avoidance of the velocity/pressure under-relaxation (UR). However, compared with the simple approach, included benefits of the simple-C method are the use of larger Courant numbers, enhanced robustness and convergence. Excellent consistency is obtained between results available in the literature and numerical solutions obtained by both simple and simple-C solvers. The segregated simple algorithm is finally reformulated in conjunction with a new slope/flux limiter function to predict fluid flow at all speeds. Numerical results show that the compressible variant of simple replicates correct shock speed, well-resolved shock front, contact discontinuity and rarefaction waves when compared with analytical solutions. Compressible laminar flows are computed to further support the accuracy and robustness of proposed algorithm. (C) 2020 The Author(s). Published by Elsevier B.V. on behalf of International Association for Mathematics and Computers in Simulation (IMACS).
The segregated simple algorithm is developed on a curvilinear collocated grid using a cell-centered finite-volume D-formulation to speculate all-speed fluid flows. The occurrence of velocity-pressure decoupling is unr...
详细信息
The segregated simple algorithm is developed on a curvilinear collocated grid using a cell-centered finite-volume D-formulation to speculate all-speed fluid flows. The occurrence of velocity-pressure decoupling is unraveled by invoking a consistently formulated dual-dissipation scheme. A modified flux-difference scheme calculates viscous fluxes at all speeds. A local Mach number based indicator is ini-tiated to impose a constraint on the actual state of compressibility, related to the density-updating in the corrector step of simple algorithm. The positivity in influence coefficients' diagonal-dominance with the corrector step is guaranteed by the current contrivance, preventing the need for an under-relaxation (UR) factor. For a selected low Mach number flow, comparisons reveal the relevance of incompressible approach, sprung from the unique compressible solver. (c) 2022 THE AUTHORS. Published by Elsevier BV on behalf of Faculty of Engineering, Ain Shams Uni-versity. This is an open access article under the CC BY license (http://***/licenses/by/4.0/).
This paper deals with a parallel implementation of the simple (Semi-Implicit Method for Pressure-Linked Equations) algorithm to numerically solve the Navier-Stokes system of equations for viscous incompressible flows....
详细信息
This paper deals with a parallel implementation of the simple (Semi-Implicit Method for Pressure-Linked Equations) algorithm to numerically solve the Navier-Stokes system of equations for viscous incompressible flows. A mechanism of interprocess communication using a mesh decomposition with virtual cells and an algebraic multigrid method is proposed. A method of distributed matrix storage and an algorithm for matrix-vector operations reducing the number of interprocess communications are presented. The results of a series of numerical experiments on structured and unstructured grids (including a problem of external aerodynamics) are presented. Based on the results obtained, an analysis of the influence of multigrid solver settings on the overall efficiency of the algorithm is made. It is shown that the parallel algorithm based on an algebraic multigrid technique proposed for the simple method makes it possible to efficiently calculate problems on hundreds of processors.
This study proposes a simple control algorithm for LCL-filter-based distribution static compensator (DSTATCOM) for power quality improvement. The proposed algorithm is able to operate DSTATCOM in both current control ...
详细信息
This study proposes a simple control algorithm for LCL-filter-based distribution static compensator (DSTATCOM) for power quality improvement. The proposed algorithm is able to operate DSTATCOM in both current control mode (CCM) and voltage control mode (VCM) by controlling point of common coupling (PCC) voltage. However, in the conventional method, controlling both filter current and PCC voltage is required to operate DSTATCOM in CCM and VCM, respectively. The achieving of CCM and VCM operations with only PCC voltage controlling method reduces the computational burden on the controller since it eliminates the necessity of filter current sensors. For smooth reference tracking, the interfacing LCL-filter parameters are designed based on the switching dynamics of hysteresis controller (HC) and harmonic order to be compensated by DSTATCOM. The total interfacing inductance value required in LCL DSTATCOM is low when compared with conventional DSTATCOM (L-filter-based DSTATCOM);therefore, the required dc-link voltage is reduced. An adaptive dc-link voltage control has been implemented to improve LCL-DSTATCOM performance in terms of switching voltage stress and efficiency. To show the efficacy of the proposed algorithm, simulation and experimental studies have been performed.
Background Existing algorithms to predict the location of an accessory pathway (AP) in Wolff-Parkinson-White Syndrome (WPW) have good sensitivity and specificity but complex with various accuracy and inter-observer ag...
详细信息
Background Existing algorithms to predict the location of an accessory pathway (AP) in Wolff-Parkinson-White Syndrome (WPW) have good sensitivity and specificity but complex with various accuracy and inter-observer agreement rates. A simple algorithm with high accuracy and inter-observer agreement rates is needed. Methods This was a cross-sectional and retrospective diagnostic study. The data were collected by total population sampling from January 2015 to January 2017. Forty-seven patients were included in the study. Data collected were pre-ablation 12-lead ECGs and ablation reports. These ECGs were evaluated by two independent observers using the simplified algorithm and compared with ablation results. Results The algorithm had a sensitivity of 45% on the left free wall, 80% on septal, 92% on the right free wall, and the specificity of 96% on the left free wall, 69% on the septal, 85% on the right free wall for AP prediction. The positive predictive value was 90% on the left free wall, 55% on the septal, and 67% on the right free wall APs. The negative predictive value was 70% on the left free wall, 88% on the septal, and 97% on right free wall AP. The positive likelihood ratio was 11.23 on the left free wall, 2.23 on septal and 6.57 on right free wall APs, and the negative likelihood ratio was 0.57 on left free wall APs, 0.28 on septal, and 0.09 on the right free wall APs. algorithm accuracy varied from 73-87%. Inter-observer agreement calculation was kappa 0.93 for left free wall AP, 0.78 for septal AP, and 0.74 for right free wall AP. Conclusion This simple algorithm has a remarkable accuracy and inter-observer agreement;therefore, it may prove to be helpful even to non-electrophysiologists and has the potential to be integrated into clinical practice
This paper introduces a new pressure-velocity coupling algorithm based on the simpleC method. The new approach considers the neighbour velocity corrections of simpleC as a Taylor series expansion, introducing a first-...
详细信息
This paper introduces a new pressure-velocity coupling algorithm based on the simpleC method. The new approach considers the neighbour velocity corrections of simpleC as a Taylor series expansion, introducing a first-order term to increase the accuracy of the approximation. The new term includes a velocity correction gradient which is assumed to be a scalar matrix constrained by means of a mass conservation equation. The stability of the method is analyzed via a Fourier decomposition of the error showing a better convergence rate than simple and simpleC for high relaxation factors. The new method is tested in two incompressible laminar flow problems. Then, the analysis is extended to a turbulent flow case. In all cases, the conclusions of the stability analysis are verified. The current proposal sets a theoretical baseline for further improvements of simple-based algorithms. (C) 2019 Elsevier Ltd. All rights reserved.
The simple algorithm for the pressure-velocity coupling has been adopted by many commercial and non-commercial CFD codes. The simple algorithm requires all normal velocities at the boundaries before solving the pressu...
详细信息
The simple algorithm for the pressure-velocity coupling has been adopted by many commercial and non-commercial CFD codes. The simple algorithm requires all normal velocities at the boundaries before solving the pressure correction equation. These velocities must imply an overall mass balance. In this paper, a new procedure to recalculate the outflow velocity for the Neumann type outflow boundary for forced convection with variable density is proposed. As an example of the utilization of this procedure, the outflow normal velocity of a circular duct was computed and the results were compared with the results using the new modified procedure.
This paper describes the details of the parallel implementation of the simple algorithm for numerical solution of the Navier-Stokes system of equations on arbitrary unstructured grids. The iteration schemes for the se...
详细信息
This paper describes the details of the parallel implementation of the simple algorithm for numerical solution of the Navier-Stokes system of equations on arbitrary unstructured grids. The iteration schemes for the serial and parallel versions of the simple algorithm are implemented. In the description of the parallel implementation, special attention is paid to computational data exchange among processors under the condition of the grid model decomposition using fictitious cells. We discuss the specific features for the storage of distributed matrices and implementation of vector-matrix operations in parallel mode. It is shown that the proposed way of matrix storage reduces the number of interprocessor exchanges. A series of numerical experiments illustrates the effect of the multigrid SLAE solver tuning on the general efficiency of the algorithm;the tuning involves the types of the cycles used (V, W, and F), the number of iterations of a smoothing operator, and the number of cells for coarsening. Two ways (direct and indirect) of efficiency evaluation for parallelization of the numerical algorithm are demonstrated. The paper presents the results of solving some internal and external flow problems with the evaluation of parallelization efficiency by two algorithms. It is shown that the proposed parallel implementation enables efficient computations for the problems on a thousand processors. Based on the results obtained, some general recommendations are made for the optimal tuning of the multigrid solver, as well as for selecting the optimal number of cells per processor.
暂无评论