{"title":"Multi-Objective Optimization and Experimental Research of Ship Form Based on Improved Bare-Bones Multi-Objective Particle Swarm Optimization Algorithm","authors":"Jie Liu, Baoji Zhang, Yuyang Lai, Liqiao Fang","doi":"10.1002/fld.5346","DOIUrl":"https://doi.org/10.1002/fld.5346","url":null,"abstract":"<div>\u0000 \u0000 <p>Ship form optimization poses a complex and high-dimensional engineering challenge. Therefore, when conducting multi-objective optimization research of ship forms, traditional intelligent optimization algorithms are prone to falling into local optima solution and difficult to converge. In order to effectively improve the diversity and convergence performance of the algorithm, this paper improves the bare-bones multi-objective particle swarm optimization (BBMOPSO) algorithm by dynamically adjusting the local and global search step sizes, and verifies the algorithm's reliability through standard function testing. Then, a multi-objective optimization design framework with high efficiency and high integration is constructed. Taking DTMB 5512 as the research case, Free Form Deformation (FFD) method is used for hull deformation, and the proposed algorithm is used for multi-objective optimization of resistance performance and motion response. Ship model tests were conducted on the DTMB 5512's original hull. And the numerical simulations were compared with the ship model tests. Finally, under the constructed multi-objective optimization design framework, satisfactory solutions were obtained through the improved algorithm, which confirms the effectiveness and practicality of the improved algorithm. The results show that the algorithm improved in this paper can provide some theoretical basis and technical support for green ship design and low-carbon shipping.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 3","pages":"267-282"},"PeriodicalIF":1.7,"publicationDate":"2024-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143117678","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Daniel B. V. Santos, Rafael Vidal, Prashant Valluri, Gustavo R. Anjos
{"title":"Computation of Deformable Interface Two-Phase Flows: A Semi-Lagrangian Finite Element Approach","authors":"Daniel B. V. Santos, Rafael Vidal, Prashant Valluri, Gustavo R. Anjos","doi":"10.1002/fld.5342","DOIUrl":"https://doi.org/10.1002/fld.5342","url":null,"abstract":"<div>\u0000 \u0000 <p>This work aims at presenting a new computational approach to study two and three dimensional two-phase flows and two dimensional coalescence phenomenon using direct numerical simulation. The flows are modeled by the incompressible Navier–Stokes equations, which are approximated by the finite element method. The Galerkin formulation is used to discretize the Navier–Stokes equations in the spatial domain and the semi-Lagrangian method is used to discretize the material derivative. In order to satisfy the Ladyzhenskaya–Babuška–Brezzi condition, high-order stable pair of elements are used, with pressure and velocity fields being calculated on different degrees of freedom in the unstructured mesh nodes. The interface is modeled by an unfitted adaptive moving mesh, where interface nodes are tracked in a Lagrangian fashion and moved with the velocity solution of the fluid motion equations. The surface tension is computed using the interface curvature and the gradient of a Heaviside function, and added in the momentum equations as a body force. In order to avoid undesired spurious modes at the interface due to high property ratios, a smooth transition between fluid properties is defined on the interface region. Several benchmark tests have been carried out to validate the proposed approach, and the obtained results have demonstrated agreement with analytical solutions and numerical results reported in the literature. A coalescence model is also proposed based on geometric criteria and results show interesting dynamics.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 3","pages":"233-243"},"PeriodicalIF":1.7,"publicationDate":"2024-10-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143115014","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Collocation Finite Element Method for the Fractional Fokker–Planck Equation","authors":"Hatice Karabenli, Alaattin Esen, Yusuf Uçar","doi":"10.1002/fld.5343","DOIUrl":"https://doi.org/10.1002/fld.5343","url":null,"abstract":"<p>In this study, the approximate results of the fractional Fokker–Planck equations have been investigated. First, finite element schemes have been obtained using collocation finite element method based on the trigonometric quintic B-spline basis functions. Then, the present method is tested on two fundamental problems having appropriate initial conditions. The newly obtained numerical results contained the error norms <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msub>\u0000 <mi>L</mi>\u0000 <mn>2</mn>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ {L}_2 $$</annotation>\u0000 </semantics></math> and <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msub>\u0000 <mi>L</mi>\u0000 <mi>∞</mi>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ {L}_{infty } $$</annotation>\u0000 </semantics></math> for various temporal and spatial steps are compared with the exact ones and other solutions. More accurate results have been obtained for large numbers of spatial and temporal elements.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 3","pages":"224-232"},"PeriodicalIF":1.7,"publicationDate":"2024-10-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5343","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143114762","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Weizhao Li, Aditya Pandare, Hong Luo, Jozsef Bakosi, Jacob Waltz
{"title":"Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods","authors":"Weizhao Li, Aditya Pandare, Hong Luo, Jozsef Bakosi, Jacob Waltz","doi":"10.1002/fld.5340","DOIUrl":"https://doi.org/10.1002/fld.5340","url":null,"abstract":"<p>A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>m</mi>\u0000 <mo>≥</mo>\u0000 <mn>2</mn>\u0000 </mrow>\u0000 <annotation>$$ mge 2 $$</annotation>\u0000 </semantics></math>) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 2","pages":"188-209"},"PeriodicalIF":1.7,"publicationDate":"2024-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5340","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143113182","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Novel lattice Boltzmann method for simulation of strongly shear thinning viscoelastic fluids","authors":"Richard Kellnberger, Tomasz Jüngst, Stephan Gekle","doi":"10.1002/fld.5335","DOIUrl":"https://doi.org/10.1002/fld.5335","url":null,"abstract":"<p>The simulation of viscoelastic liquids using the Lattice–Boltzmann method (LBM) in full three dimensions remains a formidable numerical challenge. In particular the simulation of strongly shear-thinning fluids, where the ratio between the high-shear and low-shear viscosities is large, is often prevented by stability problems. Here we present a novel approach to overcome this issue. The central idea is to artificially increase the solvent viscosity which allows the method to benefit from the very good stability properties of the LBM. To compensate for this additional viscous stress, the polymer stress is reduced by the same amount. We apply this novel method to simulate two realistic cell carrier fluids, methyl cellulose and alginate solutions, of which the latter exhibits a viscosity ratio exceeding 10,000.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 2","pages":"164-187"},"PeriodicalIF":1.7,"publicationDate":"2024-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5335","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143113181","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A finite volume scheme employing the multipoint flux approximation with diamond stencil for the diffusive-viscous wave equation on general polyhedral meshes","authors":"Di Yang, Zhiming Gao, Wenjing Yan","doi":"10.1002/fld.5336","DOIUrl":"https://doi.org/10.1002/fld.5336","url":null,"abstract":"<p>Based on three-dimensional seismic wave, simulations have become a pivotal aspect of seismic exploration. The diffusive-viscous wave equation, initially proposed by Goloshubin et al., is frequently utilized to describe seismic wave propagation in fluid-saturated media. However, obtaining numerical solutions for this equation has become an urgent issue in recent years. In this study, we present a cell-centered finite volume scheme utilizing a multipoint flux approximation that employs a “diamond stencil” on general polyhedral meshes to address the diffusive-viscous wave equation. Numerical tests exhibit that this new scheme attains optimal convergence, and its effectiveness is demonstrated through simulating vibrations induced by an earthquake source.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 2","pages":"151-163"},"PeriodicalIF":1.7,"publicationDate":"2024-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143113242","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Optimizing Winglet Cant Angle for Enhanced Aircraft Wing Performance Using CFD Simulation and Hybrid ANN-GA","authors":"Vidhit Mandia, Vipul Sharma, Yash Chandra, Gaurav Kumar, Raj Kumar Singh","doi":"10.1002/fld.5341","DOIUrl":"https://doi.org/10.1002/fld.5341","url":null,"abstract":"<div>\u0000 \u0000 <p>Winglets are an extended angled or vertical projected at the wing tips used to reduce the drag encountered during the flight of an aircraft. The main aim of this research was to study the effects of winglets on NACA 4412 airfoil at 15° angle of attack. The simulation was done on the basis of the aerodynamic properties such as lift (CL), drag (CD), and lift/drag (CL/CD) ratio for both with and without the winglets at various cant angles. The designing was carried out in ANSYS Design Modeler for both with and without winglet. Further, the meshing part was again carried out in ANSYS Mesh. K-Epsilon (two equation) turbulence model is used for the simulation at the inlet speed of 100 m/s, since it is the most common model used to simulate the mean flow characteristics for high turbulent conditions. Further, the cant angle has been optimized to get the maximum coefficient of lift using Nelder Mead, Genetic Algorithm, and Genetic Algorithm with ANN optimization techniques.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 3","pages":"211-223"},"PeriodicalIF":1.7,"publicationDate":"2024-10-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143113241","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Iason Tsetoglou, Mélody Cailler, Pierre Bénard, Ghislain Lartigue, Vincent Moureau, Julien Réveillon
{"title":"A volume-of-solid implicit volume penalty method for moving-body flows","authors":"Iason Tsetoglou, Mélody Cailler, Pierre Bénard, Ghislain Lartigue, Vincent Moureau, Julien Réveillon","doi":"10.1002/fld.5334","DOIUrl":"https://doi.org/10.1002/fld.5334","url":null,"abstract":"<p>An original Immersed Boundary Method for solving moving body flows is proposed. This method couples (i) a Lagrangian Volume-of-Solid description of the solid object avoiding conservation issues and (ii) a robust implicit volume penalty forcing embedded in a low-Mach number projection method to account for the solid's impact on the fluid dynamics. A new composite velocity field is introduced to describe both solid and fluid domains in a single set of governing equations. The accuracy of the method has been assessed on several academic cases, involving stationary or moving bodies and with different mesh resolutions. The predicted forces on the solid are in excellent agreement with body-fitted reference cases. The system of equations is also proven to be fully mass conservative. Application of the method on a two-dimensional vertical axis turbine case shows a <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>30</mn>\u0000 <mo>%</mo>\u0000 </mrow>\u0000 <annotation>$$ 30% $$</annotation>\u0000 </semantics></math> reduction in computational cost compared to a body-fitted method.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 2","pages":"117-150"},"PeriodicalIF":1.7,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5334","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143117229","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Semi-implicit Lagrangian Voronoi approximation for the incompressible Navier–Stokes equations","authors":"Ondřej Kincl, Ilya Peshkov, Walter Boscheri","doi":"10.1002/fld.5339","DOIUrl":"10.1002/fld.5339","url":null,"abstract":"<p>We introduce semi-implicit Lagrangian Voronoi approximation (SILVA), a novel numerical method for the solution of the incompressible Euler and Navier–Stokes equations, which combines the efficiency of semi-implicit time marching schemes with the robustness of time-dependent Voronoi tessellations. In SILVA, the numerical solution is stored at particles, which move with the fluid velocity and also play the role of the generators of the computational mesh. The Voronoi mesh is rapidly regenerated at each time step, allowing large deformations with topology changes. As opposed to the reconnection-based Arbitrary-Lagrangian-Eulerian schemes, we need no remapping stage. A semi-implicit scheme is devised in the context of moving Voronoi meshes to project the velocity field onto a divergence-free manifold. We validate SILVA by illustrative benchmarks, including viscous, inviscid, and multi-phase flows. Compared to its closest competitor, the Incompressible Smoothed Particle Hydrodynamics method, SILVA offers a sparser stiffness matrix and facilitates the implementation of no-slip and free-slip boundary conditions.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 1","pages":"88-115"},"PeriodicalIF":1.7,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5339","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142252123","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A new non-equilibrium modification of the \u0000 \u0000 \u0000 k\u0000 −\u0000 ω\u0000 \u0000 $$ k-omega $$\u0000 turbulence model for supersonic turbulent flows with transverse jet","authors":"Altynshash Naimanova, Assel Beketaeva","doi":"10.1002/fld.5337","DOIUrl":"10.1002/fld.5337","url":null,"abstract":"<p>The goal of this research is to propose a new modification of a non-equilibrium effect in the <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>k</mi>\u0000 <mo>−</mo>\u0000 <mi>ω</mi>\u0000 </mrow>\u0000 <annotation>$$ k-omega $$</annotation>\u0000 </semantics></math> turbulence model to better predict high-speed turbulent flows. For that, the two local compressibility coefficients are included in the balance production/dissipation terms in a specific dissipation rate equation. The specific dissipation rate reacts to changes in the local Mach number and density through these local coefficients. The developed model is applied to the numerical simulation of the spatial supersonic turbulent airflow with round hydrogen injection. In that, the effects of the proposed turbulence model on the flow field behavior (shock wave and vortex formations, shock wave/boundary layer interaction, and mixture layer) are studied via the solution of three-dimensional Favre-averaged Navier–Stokes equations with a third-order Essentially Non-Oscillatory scheme. A series of numerical experiments are performed, in which an allowable range of local constants by comparing results with experimental data is obtained. The non-equilibrium modification by simultaneous decrease of the turbulence kinetic energy and increase of the specific dissipation rate gives a good agreement of the hydrogen depth penetration with experimental data. Also, the numerical experiment of the supersonic airflow with a nitrogen jet shows wall pressure distribution is consistent well with experimental data.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 1","pages":"69-87"},"PeriodicalIF":1.7,"publicationDate":"2024-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142252122","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}