首页 > 最新文献

Computers & Fluids最新文献

英文 中文
Invariant control strategies for active flow control using graph neural networks 基于图神经网络的主动流控制的不变量控制策略
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-09-30 DOI: 10.1016/j.compfluid.2025.106854
Marius Kurz , Rohan Kaushik , Marcel Blind , Patrick Kopper , Anna Schwarz , Felix Rodach , Andrea Beck
Reinforcement learning (RL) has recently gained traction for active flow control tasks, with initial applications exploring drag mitigation via flow field augmentation around a two-dimensional cylinder. RL has since been extended to more complex turbulent flows and has shown significant potential in learning complex control strategies. However, such applications remain computationally challenging owing to its sample inefficiency and associated simulation costs. This fact is worsened by the lack of generalization capabilities of these trained policy networks, often being implicitly tied to the input configurations of their training conditions. In this work, we propose the use of graph neural networks (GNNs) to address this particular limitation, effectively increasing the range of applicability and getting more value out of the upfront RL training cost. GNNs can naturally process unstructured, three-dimensional flow data, preserving spatial relationships without the constraints of a Cartesian grid. Additionally, they incorporate rotational, reflectional, and permutation invariance into the learned control policies, thus improving generalization and thereby removing the shortcomings of commonly used convolutional neural networks (CNNs) or multilayer perceptron (MLP) architectures. To demonstrate the effectiveness of this approach, we revisit the well-established two-dimensional cylinder benchmark problem for active flow control. The RL training is implemented using Relexi, a high-performance RL framework, with flow simulations conducted in parallel using the high-order discontinuous Galerkin framework FLEXI. Our results show that GNN-based control policies achieve comparable performance to existing methods while benefiting from improved generalization properties. This work establishes GNNs as a promising architecture for RL-based flow control and highlights the capabilities of Relexi and FLEXI for large-scale RL applications in fluid dynamics.
强化学习(RL)最近在主动流动控制任务中获得了关注,最初的应用是通过在二维圆柱体周围增加流场来减少阻力。RL已经扩展到更复杂的湍流,并在学习复杂的控制策略方面显示出巨大的潜力。然而,由于其样本效率低下和相关的模拟成本,此类应用在计算上仍然具有挑战性。由于这些训练过的策略网络缺乏泛化能力,这一事实变得更加糟糕,这些策略网络通常隐含地与它们的训练条件的输入配置联系在一起。在这项工作中,我们建议使用图神经网络(gnn)来解决这一特殊限制,有效地增加适用性范围,并从前期RL训练成本中获得更多价值。gnn可以自然地处理非结构化的三维流量数据,在没有笛卡尔网格约束的情况下保持空间关系。此外,他们将旋转,反射和排列不变性纳入学习的控制策略,从而提高泛化,从而消除常用的卷积神经网络(cnn)或多层感知器(MLP)架构的缺点。为了证明这种方法的有效性,我们重新审视了主动流动控制的二维圆柱体基准问题。RL训练使用高性能RL框架Relexi实现,流模拟使用高阶不连续Galerkin框架FLEXI并行进行。我们的研究结果表明,基于gnn的控制策略可以获得与现有方法相当的性能,同时受益于改进的泛化特性。这项工作建立了gnn作为基于RL的流控制的有前途的架构,并强调了Relexi和FLEXI在流体动力学中大规模RL应用的能力。
{"title":"Invariant control strategies for active flow control using graph neural networks","authors":"Marius Kurz ,&nbsp;Rohan Kaushik ,&nbsp;Marcel Blind ,&nbsp;Patrick Kopper ,&nbsp;Anna Schwarz ,&nbsp;Felix Rodach ,&nbsp;Andrea Beck","doi":"10.1016/j.compfluid.2025.106854","DOIUrl":"10.1016/j.compfluid.2025.106854","url":null,"abstract":"<div><div>Reinforcement learning (RL) has recently gained traction for active flow control tasks, with initial applications exploring drag mitigation via flow field augmentation around a two-dimensional cylinder. RL has since been extended to more complex turbulent flows and has shown significant potential in learning complex control strategies. However, such applications remain computationally challenging owing to its sample inefficiency and associated simulation costs. This fact is worsened by the lack of generalization capabilities of these trained policy networks, often being implicitly tied to the input configurations of their training conditions. In this work, we propose the use of graph neural networks (GNNs) to address this particular limitation, effectively increasing the range of applicability and getting more <em>value</em> out of the upfront RL training cost. GNNs can naturally process unstructured, three-dimensional flow data, preserving spatial relationships without the constraints of a Cartesian grid. Additionally, they incorporate rotational, reflectional, and permutation invariance into the learned control policies, thus improving generalization and thereby removing the shortcomings of commonly used convolutional neural networks (CNNs) or multilayer perceptron (MLP) architectures. To demonstrate the effectiveness of this approach, we revisit the well-established two-dimensional cylinder benchmark problem for active flow control. The RL training is implemented using Relexi, a high-performance RL framework, with flow simulations conducted in parallel using the high-order discontinuous Galerkin framework FLEXI. Our results show that GNN-based control policies achieve comparable performance to existing methods while benefiting from improved generalization properties. This work establishes GNNs as a promising architecture for RL-based flow control and highlights the capabilities of Relexi and FLEXI for large-scale RL applications in fluid dynamics.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106854"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145218983","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A novel interface reconstruction method based on B-Spline parametric surfaces: Application to free-falling natural particle in thermo-buoyant flows 一种基于b样条参数曲面的界面重建新方法:在热浮力流动中自由落体自然粒子中的应用
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-09-24 DOI: 10.1016/j.compfluid.2025.106843
Lorenzo Vallisa , Delphine Laboureur , Maria Teresa Scelzo , Silvania Lopes , Michel De Paepe
This paper provides an innovative formulation based on B-Spline surfaces to explicitly represent liquid–solid interfaces in multiphase problems. B-Spline surfaces constitute a robust framework for representing complex geometries. Their ability to provide continuous access to surface information across the entire parametric domain offers a substantial advantage in managing boundary point distributions, particularly when compared to the traditional connectivity matrix employed in hybrid front-tracking methods and immersed boundary approaches. This work begins by introducing a definition for a closed parametric surface that fully encompasses the boundary of a solid object. Building on this foundation, it presents a parsing algorithm designed for systematic organization. The algorithm processes an unordered set of boundary points, representing the rigid body’s boundary, and arranges them into a structure aligned with the defined closed parametric surface. Leveraging the advanced geometric capabilities of B-Splines, we introduce a novel procedure for the computation of the indicator function, whose computational advantages over established methods such as the Poisson equation is then proved in a benchmark test. The presented geometrical tools are then integrated into a direct-forcing immersed boundary method. The solver is enhanced with a extension field technique, to address numerical instabilities caused by spurious oscillations arising near the immersed boundary. Additionally, it is extended to incorporate the effects of thermobuoyant flow dynamics. The solver is rigorously validated through both isothermal and non-isothermal test cases and subsequently applied to investigate the comparative free-falling behavior of a spherical object and a natural particle with an equivalent spherical volume. The results demonstrate that, despite the morphology of the natural particle exhibits only a slight deviation from that of a spherical object, its asymmetries significantly influence its trajectory.
本文提出了一种基于b样条曲面的创新公式来明确表示多相问题中的液固界面。b样条曲面构成了表示复杂几何形状的健壮框架。它们能够连续访问整个参数域的表面信息,这在管理边界点分布方面具有很大的优势,特别是与混合前沿跟踪方法和浸入边界方法中使用的传统连接矩阵相比。这项工作首先介绍了一个完全包含固体物体边界的封闭参数曲面的定义。在此基础上,提出了一种面向系统组织的解析算法。该算法对代表刚体边界的一组无序边界点进行处理,并将其排列成与所定义的封闭参数曲面对齐的结构。利用b样条的先进几何能力,我们引入了一种新的程序来计算指标函数,其计算优势优于既定的方法,如泊松方程,然后在基准测试中得到证明。然后将所提出的几何工具集成到直接强迫浸入边界法中。采用扩展场技术对求解器进行了改进,以解决在浸入边界附近产生的伪振荡引起的数值不稳定性问题。此外,它被扩展到包括热浮力流动动力学的影响。通过等温和非等温测试案例对求解器进行了严格验证,并随后应用于研究具有等效球形体积的球形物体和自然粒子的自由落体行为的比较。结果表明,尽管天然粒子的形态与球形物体的形态只有轻微的偏差,但其不对称性显著影响其轨迹。
{"title":"A novel interface reconstruction method based on B-Spline parametric surfaces: Application to free-falling natural particle in thermo-buoyant flows","authors":"Lorenzo Vallisa ,&nbsp;Delphine Laboureur ,&nbsp;Maria Teresa Scelzo ,&nbsp;Silvania Lopes ,&nbsp;Michel De Paepe","doi":"10.1016/j.compfluid.2025.106843","DOIUrl":"10.1016/j.compfluid.2025.106843","url":null,"abstract":"<div><div>This paper provides an innovative formulation based on B-Spline surfaces to explicitly represent liquid–solid interfaces in multiphase problems. B-Spline surfaces constitute a robust framework for representing complex geometries. Their ability to provide continuous access to surface information across the entire parametric domain offers a substantial advantage in managing boundary point distributions, particularly when compared to the traditional connectivity matrix employed in hybrid front-tracking methods and immersed boundary approaches. This work begins by introducing a definition for a closed parametric surface that fully encompasses the boundary of a solid object. Building on this foundation, it presents a parsing algorithm designed for systematic organization. The algorithm processes an unordered set of boundary points, representing the rigid body’s boundary, and arranges them into a structure aligned with the defined closed parametric surface. Leveraging the advanced geometric capabilities of B-Splines, we introduce a novel procedure for the computation of the indicator function, whose computational advantages over established methods such as the Poisson equation is then proved in a benchmark test. The presented geometrical tools are then integrated into a direct-forcing immersed boundary method. The solver is enhanced with a <em>extension field</em> technique, to address numerical instabilities caused by spurious oscillations arising near the immersed boundary. Additionally, it is extended to incorporate the effects of thermobuoyant flow dynamics. The solver is rigorously validated through both isothermal and non-isothermal test cases and subsequently applied to investigate the comparative free-falling behavior of a spherical object and a natural particle with an equivalent spherical volume. The results demonstrate that, despite the morphology of the natural particle exhibits only a slight deviation from that of a spherical object, its asymmetries significantly influence its trajectory.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106843"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145263970","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Applying a compact porous media model to numerically derive resistance coefficients for lattice structures 应用致密多孔介质模型数值推导晶格结构的阻力系数
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-08-16 DOI: 10.1016/j.compfluid.2025.106798
Kyle Pittman , Jacob Riglin , Jay Chen , Cesar Dominguez , Marvin Davis , Rami Batrice
Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×106 and 2.90×106 [1/(Pa*m*s)] for the permeability coefficients, 6.37×101 and 5.44×101 [m2/kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.
增材制造允许探索各种几何形状,以达到特定的工程标准。晶格是一种具有独特性质的几何结构,包括周期性重复的结构,根据Darcy-Forchheimer方程,流动可以通过它们被表示为多孔介质。这些方程的系数通常是通过实验推导出来的,但这项工作证明了用CFD数值推导它们的能力。利用Ansys Fluent软件采用三维稳态reynolds -average Navier-Stokes和k-ω剪切应力输运湍流模型进行了模拟。研究了三种晶格几何形状,并推导了阻力系数。通过外部公布的类似几何形状的数据验证了该方法的有效性,并使用网格收敛指数方法计算了所有模拟的网格收敛性。壁面粗糙度被证明对结果有不可忽略的影响,并且粗糙度值被考虑为主要焦点八面体几何,其中光滑壁面和粗糙壁面系数都得到了推导。在1.0 [m/s]条件下,八面体孔隙率系数为2.89×106和2.90×106 [1/(Pa*m*s)],惯性阻力系数为6.37×101和5.44×101 [m2/kg],光滑壁面和粗糙壁面的最大压降分别为5116.7 [Pa]和4429.5 [Pa]。导出的数值方法能够快速探索和优化新的晶格设计,以适应不同的工程应用。
{"title":"Applying a compact porous media model to numerically derive resistance coefficients for lattice structures","authors":"Kyle Pittman ,&nbsp;Jacob Riglin ,&nbsp;Jay Chen ,&nbsp;Cesar Dominguez ,&nbsp;Marvin Davis ,&nbsp;Rami Batrice","doi":"10.1016/j.compfluid.2025.106798","DOIUrl":"10.1016/j.compfluid.2025.106798","url":null,"abstract":"<div><div>Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10<sup>6</sup> and 2.90×10<sup>6</sup> [1/(Pa*m*s)] for the permeability coefficients, 6.37×10<sup>1</sup> and 5.44×10<sup>1</sup> [m<sup>2</sup>/kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106798"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145264408","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
An Eulerian–Lagrangian approach to simulate turbidity currents 用欧拉-拉格朗日方法模拟浊度流
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-09-25 DOI: 10.1016/j.compfluid.2025.106832
Boyang Chen , Zhen Liu , Bruño Fraga
We present an Eulerian–Lagrangian four-way coupled algorithm within a large-eddy simulation framework to simulate turbidity currents. Our approach preserves the particle-driven, dispersed nature of such currents with significantly reduced reliance on semi-empirical parametrisation. We are capable of reproducing key processes such as the entrainment of fluid within the particle-laden current and the settlement and re-suspension of solid particles. Particle interactions are handled using a soft-sphere collision model. Our results are successfully validated versus experimental results. We investigate a lock-exchange setup in a numerical flume, comparing predicted front velocities and deposition profiles with experimental measurements. Furthermore, we analyse the differences between particle-driven and gravity-driven currents (simulated via an Eulerian–Eulerian approach), focusing on propagation velocity, transition to turbulence and the generation of coherent structures in the shear layer.
We use our model to examine the evolution and driving mechanisms of turbidity currents. We describe in detail how Kelvin–Helmholtz singularities evolve into a well-defined current head and sediment trail, while also accounting for the mechanical effect of lifting the sliding gate that initially separates the laden and unladen liquids. The dynamics of particle settling and re-suspension and their correlation with bed shear stress and the current’s reach are predicted, showing good agreement with experimental data. The local friction velocity at the flume’s bed peaks at the current head during early development and is later redistributed via a re-suspension mechanism linked to instantaneous turbulent structures. Finally, an energy budget analysis reveals that turbulent kinetic energy dissipation is primarily due to the settling of solid particles.
我们提出了一个欧拉-拉格朗日四路耦合算法在大涡模拟框架内模拟浊度流。我们的方法保留了这种电流的粒子驱动,分散的性质,大大减少了对半经验参数化的依赖。我们能够重现关键的过程,如在含颗粒的电流中流体的夹带以及固体颗粒的沉降和再悬浮。粒子相互作用使用软球碰撞模型处理。我们的结果与实验结果进行了验证。我们研究了一个锁交换设置在一个数值水槽,比较预测的前速度和沉积剖面与实验测量。此外,我们分析了粒子驱动和重力驱动电流之间的差异(通过欧拉-欧拉方法模拟),重点关注传播速度,向湍流的过渡以及剪切层中相干结构的产生。我们使用我们的模型来研究浊度流的演变和驱动机制。我们详细描述了开尔文-亥姆霍兹奇点是如何演变成一个定义明确的水流头和沉积物轨迹的,同时也考虑了提起最初分离有载液体和无载液体的滑动闸门的机械效应。预测了颗粒沉降和再悬浮的动力学过程及其与床层剪切应力和水流长度的关系,与实验数据吻合较好。在早期发展过程中,水槽底部的局部摩擦速度在水流头部达到峰值,随后通过与瞬时湍流结构相关的再悬浮机制重新分配。最后,能量收支分析表明,湍流动能耗散主要是由于固体颗粒的沉降。
{"title":"An Eulerian–Lagrangian approach to simulate turbidity currents","authors":"Boyang Chen ,&nbsp;Zhen Liu ,&nbsp;Bruño Fraga","doi":"10.1016/j.compfluid.2025.106832","DOIUrl":"10.1016/j.compfluid.2025.106832","url":null,"abstract":"<div><div>We present an Eulerian–Lagrangian four-way coupled algorithm within a large-eddy simulation framework to simulate turbidity currents. Our approach preserves the particle-driven, dispersed nature of such currents with significantly reduced reliance on semi-empirical parametrisation. We are capable of reproducing key processes such as the entrainment of fluid within the particle-laden current and the settlement and re-suspension of solid particles. Particle interactions are handled using a soft-sphere collision model. Our results are successfully validated versus experimental results. We investigate a lock-exchange setup in a numerical flume, comparing predicted front velocities and deposition profiles with experimental measurements. Furthermore, we analyse the differences between particle-driven and gravity-driven currents (simulated via an Eulerian–Eulerian approach), focusing on propagation velocity, transition to turbulence and the generation of coherent structures in the shear layer.</div><div>We use our model to examine the evolution and driving mechanisms of turbidity currents. We describe in detail how Kelvin–Helmholtz singularities evolve into a well-defined current head and sediment trail, while also accounting for the mechanical effect of lifting the sliding gate that initially separates the laden and unladen liquids. The dynamics of particle settling and re-suspension and their correlation with bed shear stress and the current’s reach are predicted, showing good agreement with experimental data. The local friction velocity at the flume’s bed peaks at the current head during early development and is later redistributed via a re-suspension mechanism linked to instantaneous turbulent structures. Finally, an energy budget analysis reveals that turbulent kinetic energy dissipation is primarily due to the settling of solid particles.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106832"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145264458","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Comparison of entropy stable collocation high-order DG methods for compressible turbulent flows 可压缩湍流的熵稳定配置高阶DG方法比较
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-10-16 DOI: 10.1016/j.compfluid.2025.106874
Anna Schwarz , Daniel Kempf , Jens Keim , Patrick Kopper , Christian Rohde , Andrea Beck
High-order methods are well-suited for the numerical simulation of complex compressible turbulent flows, but require additional stabilization techniques to capture instabilities arising from the underlying non-linear hyperbolic equations. This paper provides a detailed comparison of the effectiveness of entropy stable discontinuous Galerkin methods for the stabilization of compressible (wall-bounded) turbulent flows. For this investigation, an entropy stable discontinuous Galerkin spectral element method is applied on Gauss–Legendre and Gauss–Lobatto nodes. In the compressible regime, an additional stabilization technique for shock capturing based on a convex blending of a low-order finite volume with the high-order discontinuous Galerkin operator is utilized. The present investigation provides a systematic study from convergence tests, to the Taylor–Green vortex and finally to a more intricate turbulent wall-bounded 3D diffuser flow, encompassing both weakly compressible and compressible regimes. The comparison demonstrates that the DGSEM on Gauss–Lobatto nodes is generally less accurate for an equal amount of degrees of freedom. Conversely, it is faster than the DGSEM on Gauss–Legendre nodes due to a less severe time step restriction and simpler numerical operator. A performance comparison reveals that the DGSEM on Gauss–Lobatto nodes generally outperforms the DGSEM on Gauss nodes for under-resolved turbulence in the subsonic regime on a periodic domain. Conversely, the opposite effect can be observed for wall-bounded flows as well as the supersonic regime, the latter depending of course on the chosen shock-capturing scheme. To the author’s knowledge, this is the first time for which a comparison of entropy stable DGSEM on Gauss–Lobatto and Gauss–Legendre has been performed for compressible, wall-bounded turbulent flows with separation.
高阶方法非常适合于复杂的可压缩湍流的数值模拟,但需要额外的稳定技术来捕捉由潜在的非线性双曲方程引起的不稳定性。本文详细比较了熵稳定不连续伽辽金方法在可压缩(有壁)湍流稳定中的有效性。在高斯-勒让德节点和高斯-洛巴托节点上应用熵稳定的不连续伽辽金谱元方法。在可压缩状态下,采用了基于低阶有限体积与高阶不连续伽辽金算子的凸混合的附加稳定技术来捕获激波。本研究提供了一个系统的研究,从收敛测试,到泰勒-格林涡,最后到一个更复杂的湍流壁界三维扩散流,包括弱可压缩和可压缩两种状态。比较表明,在相同的自由度下,高斯-洛巴托节点上的DGSEM通常精度较低。相反,由于时间步长限制较轻和数值运算符较简单,它比高斯-勒让德节点上的DGSEM更快。性能比较表明,在周期域上,对于亚音速区域的欠分辨湍流,高斯- lobatto节点上的DGSEM通常优于高斯节点上的DGSEM。相反,在有壁流动和超音速状态下,可以观察到相反的效果,后者当然取决于所选择的激波捕获方案。据作者所知,这是第一次在高斯-洛巴托和高斯-勒让德的熵稳定DGSEM上对可压缩的、有壁界的分离湍流进行了比较。
{"title":"Comparison of entropy stable collocation high-order DG methods for compressible turbulent flows","authors":"Anna Schwarz ,&nbsp;Daniel Kempf ,&nbsp;Jens Keim ,&nbsp;Patrick Kopper ,&nbsp;Christian Rohde ,&nbsp;Andrea Beck","doi":"10.1016/j.compfluid.2025.106874","DOIUrl":"10.1016/j.compfluid.2025.106874","url":null,"abstract":"<div><div>High-order methods are well-suited for the numerical simulation of complex compressible turbulent flows, but require additional stabilization techniques to capture instabilities arising from the underlying non-linear hyperbolic equations. This paper provides a detailed comparison of the effectiveness of entropy stable discontinuous Galerkin methods for the stabilization of compressible (wall-bounded) turbulent flows. For this investigation, an entropy stable discontinuous Galerkin spectral element method is applied on Gauss–Legendre and Gauss–Lobatto nodes. In the compressible regime, an additional stabilization technique for shock capturing based on a convex blending of a low-order finite volume with the high-order discontinuous Galerkin operator is utilized. The present investigation provides a systematic study from convergence tests, to the Taylor–Green vortex and finally to a more intricate turbulent wall-bounded 3D diffuser flow, encompassing both weakly compressible and compressible regimes. The comparison demonstrates that the DGSEM on Gauss–Lobatto nodes is generally less accurate for an equal amount of degrees of freedom. Conversely, it is faster than the DGSEM on Gauss–Legendre nodes due to a less severe time step restriction and simpler numerical operator. A performance comparison reveals that the DGSEM on Gauss–Lobatto nodes generally outperforms the DGSEM on Gauss nodes for under-resolved turbulence in the subsonic regime on a periodic domain. Conversely, the opposite effect can be observed for wall-bounded flows as well as the supersonic regime, the latter depending of course on the chosen shock-capturing scheme. To the author’s knowledge, this is the first time for which a comparison of entropy stable DGSEM on Gauss–Lobatto and Gauss–Legendre has been performed for compressible, wall-bounded turbulent flows with separation.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106874"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145359265","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Physics appropriate interface capturing reconstruction approach for viscous compressible multicomponent flows 粘性可压缩多组分流动的物理适宜界面捕获重建方法
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-10-09 DOI: 10.1016/j.compfluid.2025.106858
Amareshwara Sainadh Chamarthi
The paper proposes a physically consistent numerical discretization approach for simulating viscous compressible multicomponent flows. It has two main contributions. First, a contact discontinuity (and material interface) detector is developed. In those regions of contact discontinuities, the THINC (Tangent of Hyperbola for INterface Capturing) approach is used for reconstructing appropriate variables (phasic densities). For other flow regions, the variables are reconstructed using the Monotonicity-preserving (MP) scheme (or Weighted essentially non-oscillatory scheme (WENO)). For reconstruction in the characteristic space, the THINC approach is used only for the contact (or entropy) wave and volume fractions. For the reconstruction of primitive variables, the THINC approach is used for phasic densities and volume fractions only, offering an effective solution for reducing dissipation errors near contact discontinuities. The numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods. The second contribution is the development of an algorithm that uses a central reconstruction scheme for the tangential velocities, as they are continuous across material interfaces in viscous flows. In this regard, the Ducros sensor (a shock detector that cannot detect material interfaces) is employed to compute the tangential velocities using a central scheme across material interfaces. Using the central scheme does not produce any oscillations at the material interface. The proposed approach is thoroughly validated with several benchmark test cases for compressible multicomponent flows, highlighting its advantages. The physics appropriate approach also shown to prevent spurious vortices, despite being formally second-order accurate for nonlinear problems, on a coarser mesh than a genuinely high-order accurate method.
本文提出了一种物理一致的数值离散方法来模拟粘性可压缩多分量流动。它有两个主要贡献。首先,研制了接触不连续(和材料界面)探测器。在这些接触不连续区域,使用THINC(界面捕获双曲线切线)方法重建适当的变量(相密度)。对于其他流动区域,使用保持单调性(MP)格式(或加权本质非振荡格式(WENO))重建变量。对于特征空间的重建,THINC方法仅用于接触(或熵)波和体积分数。对于原始变量的重建,THINC方法仅用于相密度和体积分数,为减少接触不连续点附近的耗散误差提供了有效的解决方案。基准试验的数值结果表明,与现有方法相比,该方法能较好地捕捉材料界面。第二个贡献是开发了一种算法,该算法使用切向速度的中心重建方案,因为它们在粘性流动中跨越材料界面是连续的。在这方面,Ducros传感器(一种不能检测材料界面的冲击探测器)被用于使用跨材料界面的中心方案来计算切向速度。采用中心方案不会在材料界面处产生任何振荡。通过对可压缩多组分流进行基准测试,充分验证了该方法的优越性。物理上适当的方法也被证明可以防止虚假的漩涡,尽管在非线性问题上是正式的二阶精度,但在更粗糙的网格上比真正的高阶精度方法。
{"title":"Physics appropriate interface capturing reconstruction approach for viscous compressible multicomponent flows","authors":"Amareshwara Sainadh Chamarthi","doi":"10.1016/j.compfluid.2025.106858","DOIUrl":"10.1016/j.compfluid.2025.106858","url":null,"abstract":"<div><div>The paper proposes a physically consistent numerical discretization approach for simulating viscous compressible multicomponent flows. It has two main contributions. First, a contact discontinuity (and material interface) detector is developed. In those regions of contact discontinuities, the THINC (Tangent of Hyperbola for INterface Capturing) approach is used for reconstructing appropriate variables (phasic densities). For other flow regions, the variables are reconstructed using the Monotonicity-preserving (MP) scheme (or Weighted essentially non-oscillatory scheme (WENO)). For reconstruction in the characteristic space, the THINC approach is used only for the contact (or entropy) wave and volume fractions. For the reconstruction of primitive variables, the THINC approach is used for phasic densities and volume fractions only, offering an effective solution for reducing dissipation errors near contact discontinuities. The numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods. The second contribution is the development of an algorithm that uses a central reconstruction scheme for the tangential velocities, as they are continuous across material interfaces in viscous flows. In this regard, the Ducros sensor (a shock detector that cannot detect material interfaces) is employed to compute the tangential velocities using a central scheme across material interfaces. Using the central scheme does not produce any oscillations at the material interface. The proposed approach is thoroughly validated with several benchmark test cases for compressible multicomponent flows, highlighting its advantages. The physics appropriate approach also shown to prevent spurious vortices, despite being formally second-order accurate for nonlinear problems, on a coarser mesh than a genuinely high-order accurate method.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106858"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145264410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Obtaining converged flow solutions using quantum annealing 利用量子退火获得收敛流解
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-10-08 DOI: 10.1016/j.compfluid.2025.106856
Takahito Asaga , Yuichi Kuya
This paper proposes numerical methods to obtain converged flow solutions by quantum annealing. The proposed quantum annealing methods are developed for lattice gas automata (LGA) and finite difference methods (FDMs). The quadratic unconstrained binary optimization (QUBO) model for LGA consists of the cost functions for the steady-state flow condition, collision law condition, boundary condition, and flow field condition. In contrast, the QUBO model for FDMs is built directly from the discretized governing equations expressed in a binary form. In the numerical experiments of channel flows, both proposed methods successfully extract the converged velocity profiles from a large number of flow state combinations by quantum annealing. The obtained solutions closely match those obtained by the conventional or analytical approach. Since, due to the difference in characteristics between LGA and FDMs, FDMs can reduce the scale of combinatorial optimization problems more efficiently than LGA, the proposed FDM-based method obtains more accurate solutions than the proposed LGA-based method.
本文提出了用量子退火法求收敛流解的数值方法。提出了晶格气体自动机(LGA)和有限差分方法(fdm)的量子退火方法。LGA的二次无约束二元优化(QUBO)模型由稳态流条件、碰撞律条件、边界条件和流场条件的代价函数组成。相比之下,fdm的QUBO模型是直接从以二进制形式表示的离散控制方程中建立的。在通道流动的数值实验中,两种方法都成功地通过量子退火从大量流态组合中提取了收敛速度分布。所得到的解与传统方法或分析方法得到的解非常接近。由于LGA和fdm的特性不同,fdm比LGA更能有效地减少组合优化问题的规模,因此基于fdm的方法比基于LGA的方法得到更精确的解。
{"title":"Obtaining converged flow solutions using quantum annealing","authors":"Takahito Asaga ,&nbsp;Yuichi Kuya","doi":"10.1016/j.compfluid.2025.106856","DOIUrl":"10.1016/j.compfluid.2025.106856","url":null,"abstract":"<div><div>This paper proposes numerical methods to obtain converged flow solutions by quantum annealing. The proposed quantum annealing methods are developed for lattice gas automata (LGA) and finite difference methods (FDMs). The quadratic unconstrained binary optimization (QUBO) model for LGA consists of the cost functions for the steady-state flow condition, collision law condition, boundary condition, and flow field condition. In contrast, the QUBO model for FDMs is built directly from the discretized governing equations expressed in a binary form. In the numerical experiments of channel flows, both proposed methods successfully extract the converged velocity profiles from a large number of flow state combinations by quantum annealing. The obtained solutions closely match those obtained by the conventional or analytical approach. Since, due to the difference in characteristics between LGA and FDMs, FDMs can reduce the scale of combinatorial optimization problems more efficiently than LGA, the proposed FDM-based method obtains more accurate solutions than the proposed LGA-based method.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106856"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145325681","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
A kinetic scheme based on positivity preservation with exact shock capture 一种基于精确冲击捕获的正性保持的动力学方案
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-10-08 DOI: 10.1016/j.compfluid.2025.106857
Shashi Shekhar Roy , S.V. Raghurama Rao
In this paper, we present a kinetic model with flexible velocities that satisfy positivity preservation conditions for the Euler equations. Our 1D kinetic model consists of two velocities and employs both the asymmetrical and symmetrical models. Switching between the two models is governed by our formulation of kinetic relative entropy, together with an additional criterion that ensures a robust and accurate scheme yielding entropic results. In 2D, we introduce a novel three-velocity kinetic model, defined to ensure a locally one-dimensional formulation for the resulting macroscopic normal flux. For first order accuracy, we also obtain a limit on the time step that ensures positivity preservation. The resulting numerical scheme captures grid-aligned steady shocks exactly. Several benchmark compressible flow test cases are solved in 1D and 2D to demonstrate the efficacy of the proposed solver.
本文给出了满足欧拉方程正保持条件的具有柔性速度的动力学模型。我们的一维动力学模型由两种速度组成,采用不对称和对称两种模型。在两个模型之间的切换由我们的动力学相对熵公式控制,以及一个额外的标准,以确保产生熵结果的鲁棒和准确的方案。在二维中,我们引入了一种新的三速度动力学模型,该模型的定义是为了确保所得宏观法向通量的局部一维公式。对于一阶精度,我们还得到了保证正性保持的时间步长极限。所得到的数值格式准确地捕获了网格对齐的稳态冲击。通过对一维和二维可压缩流测试用例的求解,验证了该算法的有效性。
{"title":"A kinetic scheme based on positivity preservation with exact shock capture","authors":"Shashi Shekhar Roy ,&nbsp;S.V. Raghurama Rao","doi":"10.1016/j.compfluid.2025.106857","DOIUrl":"10.1016/j.compfluid.2025.106857","url":null,"abstract":"<div><div>In this paper, we present a kinetic model with flexible velocities that satisfy positivity preservation conditions for the Euler equations. Our 1D kinetic model consists of two velocities and employs both the asymmetrical and symmetrical models. Switching between the two models is governed by our formulation of kinetic relative entropy, together with an additional criterion that ensures a robust and accurate scheme yielding entropic results. In 2D, we introduce a novel three-velocity kinetic model, defined to ensure a locally one-dimensional formulation for the resulting macroscopic normal flux. For first order accuracy, we also obtain a limit on the time step that ensures positivity preservation. The resulting numerical scheme captures grid-aligned steady shocks exactly. Several benchmark compressible flow test cases are solved in 1D and 2D to demonstrate the efficacy of the proposed solver.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106857"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145325684","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Transfer function of low-pass filters on unstructured grids 非结构网格上低通滤波器的传递函数
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-10-13 DOI: 10.1016/j.compfluid.2025.106873
Hiroyuki Asada, Kanako Maruyama, Soshi Kawai
Low-pass filters designed on unstructured grids are investigated in terms of a transfer function in the wavenumber space. The transfer functions on unstructured grids are derived, and the properties of the low-pass filters for removing high-wavenumber components and inducing phase errors are investigated through the derived transfer functions. The transfer function reveals that the low-pass filters on unstructured grids can achieve the property that higher-wavenumber components are removed more by adjusting a filter coefficient to a small value, whereas large filter coefficients induce unfavorable amplifications of high-wavenumber components. The presence of phase errors induced by the low-pass filters on triangle unstructured cells is also found by the transfer function. Furthermore, the transfer function shows that the numerical methods for evaluating the gradients that appear in the filter formulation affect the characteristics of the low-pass filters and that the simplest central scheme can have an advantage in terms of retaining numerical robustness by removing high-wavenumber components compared to the edge-normal augmentation scheme. The numerical experiments of inviscid Taylor–Green vortex and shock-vortex interaction are also conducted with the low-pass filter coupled with the non-dissipative kinetic energy and entropy preserving (KEEP) scheme on unstructured grids, demonstrating the validity of the present transfer function of the low-pass filter.
从波数空间传递函数的角度研究了非结构网格低通滤波器的设计。推导了非结构栅格上的传递函数,并通过推导的传递函数研究了低通滤波器去除高波数分量和诱导相位误差的特性。传递函数表明,非结构网格上的低通滤波器通过调整滤波器系数,可以实现高波数分量被去除的特性,而大的滤波器系数会导致高波数分量的不利放大。通过传递函数也发现了低通滤波器在三角形非结构单元上引起的相位误差。此外,传递函数表明,用于评估滤波器公式中出现的梯度的数值方法会影响低通滤波器的特性,并且与边缘法向增强方案相比,最简单的中心方案在通过去除高波数分量来保持数值鲁棒性方面具有优势。利用低通滤波器与非结构网格上的非耗散动能和熵保持(KEEP)方案耦合进行了无粘Taylor-Green涡和激波-涡相互作用的数值实验,验证了所提出的低通滤波器传递函数的有效性。
{"title":"Transfer function of low-pass filters on unstructured grids","authors":"Hiroyuki Asada,&nbsp;Kanako Maruyama,&nbsp;Soshi Kawai","doi":"10.1016/j.compfluid.2025.106873","DOIUrl":"10.1016/j.compfluid.2025.106873","url":null,"abstract":"<div><div>Low-pass filters designed on unstructured grids are investigated in terms of a transfer function in the wavenumber space. The transfer functions on unstructured grids are derived, and the properties of the low-pass filters for removing high-wavenumber components and inducing phase errors are investigated through the derived transfer functions. The transfer function reveals that the low-pass filters on unstructured grids can achieve the property that higher-wavenumber components are removed more by adjusting a filter coefficient to a small value, whereas large filter coefficients induce unfavorable amplifications of high-wavenumber components. The presence of phase errors induced by the low-pass filters on triangle unstructured cells is also found by the transfer function. Furthermore, the transfer function shows that the numerical methods for evaluating the gradients that appear in the filter formulation affect the characteristics of the low-pass filters and that the simplest central scheme can have an advantage in terms of retaining numerical robustness by removing high-wavenumber components compared to the edge-normal augmentation scheme. The numerical experiments of inviscid Taylor–Green vortex and shock-vortex interaction are also conducted with the low-pass filter coupled with the non-dissipative kinetic energy and entropy preserving (KEEP) scheme on unstructured grids, demonstrating the validity of the present transfer function of the low-pass filter.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106873"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145325679","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
Convergence acceleration algorithms for non-overlapping domain decomposition in near-wall turbulence modeling 近壁湍流建模中非重叠区域分解的收敛加速算法
IF 3 3区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2025-12-15 Epub Date: 2025-10-03 DOI: 10.1016/j.compfluid.2025.106855
N. Smirnova , S. Utyuzhnikov , V. Titarev , M. Petrov
In turbulence modeling the resolution of near-wall boundary layer requires most of computing time. The near-wall domain decomposition (NDD) approach proved to be efficient in tackling this problem. It represents a trade-off between computing time and accuracy. In this method, the computational domain is divided into two non-overlapping subdomains: the inner layer and region outside. The interface boundary conditions of Robin type are set by transferring the boundary conditions from the wall to the interface boundary. In contrast to the exact NDD, in the approximate NDD a simplified system of equations is solved in the near-wall subdomain. In the current paper a variant of the exact NDD is proposed, that uses an operator corresponding to the approximate NDD approach as a preconditioner. To improve the efficiency of NDD methods the GMRES method is applied. The efficacy of NDD algorithms are compared against for the low-Reynolds-number model.
在湍流模拟中,近壁边界层的分辨率占用了大量的计算时间。近壁区域分解(NDD)方法被证明是解决这一问题的有效方法。它代表了计算时间和准确性之间的权衡。该方法将计算域划分为两个不重叠的子域:内层和外部区域。Robin型界面边界条件是通过将边界条件从壁面转移到界面边界来设置的。与精确NDD相比,近似NDD是在近壁子域中求解一个简化的方程组。本文提出了精确NDD的一种变体,它使用与近似NDD方法相对应的算子作为前置条件。为了提高NDD方法的效率,引入了GMRES方法。比较了NDD算法在低雷诺数模型下的有效性。
{"title":"Convergence acceleration algorithms for non-overlapping domain decomposition in near-wall turbulence modeling","authors":"N. Smirnova ,&nbsp;S. Utyuzhnikov ,&nbsp;V. Titarev ,&nbsp;M. Petrov","doi":"10.1016/j.compfluid.2025.106855","DOIUrl":"10.1016/j.compfluid.2025.106855","url":null,"abstract":"<div><div>In turbulence modeling the resolution of near-wall boundary layer requires most of computing time. The near-wall domain decomposition (NDD) approach proved to be efficient in tackling this problem. It represents a trade-off between computing time and accuracy. In this method, the computational domain is divided into two non-overlapping subdomains: the inner layer and region outside. The interface boundary conditions of Robin type are set by transferring the boundary conditions from the wall to the interface boundary. In contrast to the exact NDD, in the approximate NDD a simplified system of equations is solved in the near-wall subdomain. In the current paper a variant of the exact NDD is proposed, that uses an operator corresponding to the approximate NDD approach as a preconditioner. To improve the efficiency of NDD methods the GMRES method is applied. The efficacy of NDD algorithms are compared against for the low-Reynolds-number model.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"303 ","pages":"Article 106855"},"PeriodicalIF":3.0,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145263969","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
引用次数: 0
期刊
Computers & Fluids
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1