当前位置：首页 > news >正文

大规模CFD仿真计算中，SIMPLE或者PISO算法中加速压力场方程迭代求解

news 来源：原创 2025/6/20 15:12:33

文章目录

在SIMSOL或PISO算法中加速压力场方程的迭代求解是提高CFD计算效率的关键。以下从算法优化、数值技巧和并行计算等方面总结加速策略：
- - **1. 压力方程求解器的选择与优化**
  - **2. 算法层面的加速**
  - **3. 离散格式与网格优化**
  - **4. 并行计算与硬件加速**
  - **5. 代码级优化**
  - **6. 其他技巧**
  - **实际应用建议**

在SIMSOL或PISO算法中加速压力场方程的迭代求解是提高CFD计算效率的关键。以下从算法优化、数值技巧和并行计算等方面总结加速策略：

1. 压力方程求解器的选择与优化

代数多重网格法（AMG）
- 压力方程（泊松型）的求解通常占主要计算时间，AMG对这类椭圆型方程具有最优的收敛速度（接近O(n)）。
- 优先选用基于AMG的求解器（如Hypre、PETSc或开源库AMGCL）。
- 调整AMG参数：粗化策略（Ruge-Stuben或PMIS）、平滑迭代次数（通常1-2次Gauss-Seidel）、循环类型（V-cycle比W-cycle更快）。
Krylov子空间方法预处理
- 结合AMG预处理的CG或GMRES方法，适用于非对称矩阵（如瞬态PISO）。
- 示例：PCG + AMG 或 Flexible GMRES + AMG。
几何多重网格（GMG）
- 如果网格结构规则，GMG比AMG更高效（但适应性较差）。

2. 算法层面的加速

SIMPLE家族的变体选择
- SIMPLEC：通过近似抵消速度-压力耦合项，减少压力修正的欠松弛需求（可增大松弛因子）。
- PISO：瞬态问题中通过多次修正步（通常2-3次）显式处理耦合，减少迭代次数。
- SIMPLE-R：通过重构压力方程减少误差传递。
松弛因子优化
- SIMPLE中压力松弛因子（αₚ）通常取0.1-0.3，但可动态调整：
  - 初始阶段用较小值（如0.1），稳定后逐步增大（如0.3）。
  - 根据残差变化自适应调整（需谨慎避免发散）。
初始猜测优化
- 利用前一时间步的解或外插法提供更好的初始猜测，减少迭代次数。

3. 离散格式与网格优化

Rhie-Chow插值改进
- 确保动量插值避免压力振荡，减少非物理波动导致的额外迭代。
梯度重构方法
- 采用最小二乘法或Green-Gauss高阶重构，提高压力梯度计算精度。
网格质量
- 避免高长宽比或扭曲网格，否则会导致压力方程条件数恶化。
- 局部加密关键区域（如边界层），但需保持过渡平滑。

4. 并行计算与硬件加速

区域分解与MPI并行
- 将计算域划分为多个子区域，使用Schur补或重叠区域法处理边界交换。
- 确保压力求解器支持分布式计算（如Hypre）。
GPU加速
- 使用CUDA实现的AMG求解器（如AmgX）或定制Krylov求解器，可显著提升速度。
- 注意：内存带宽是瓶颈，需优化数据局部性。
混合并行（MPI+OpenMP）
- 在多核节点上结合进程级和线程级并行。

5. 代码级优化

矩阵存储格式
- 使用压缩稀疏行（CSR）或块格式（Block-CSR）存储压力矩阵，减少内存访问开销。
残差收敛准则
- 设置合理的相对残差容差（如1e-5），避免过度求解。
Jacobian矩阵重用
- 在非线性迭代中，若矩阵变化较小，可复用前几步的矩阵或预处理子。

6. 其他技巧

时间步长策略
- 瞬态模拟中，采用自适应时间步长（如CFL条件控制），避免过小步长导致冗余计算。
耦合求解尝试
- 对某些问题，可尝试完全耦合算法（如Monolithic方法），但需更多内存。

实际应用建议

工具选择：
- OpenFOAM用户可切换至GAMG（几何代数多重网格）或调用外部库（如Hypre）。
- 商业软件（如Fluent）中启用AMG并调整参数。
性能分析：
- 使用Profiler（如VTune、gprof）定位压力求解阶段的瓶颈。
验证稳定性：
- 加速后需检查质量守恒性和残差收敛性，避免虚假收敛。

通过综合上述方法，通常可将压力场求解时间缩短50%以上，尤其在大规模问题时AMG和并行优化效果显著。

大规模CFD仿真计算中，SIMPLE或者PISO算法中加速压力场方程迭代求解

文章目录在SIMSOL或PISO算法中加速压力场方程的迭代求解是提高CFD计算效率的关键。以下从算法优化、数值技巧和并行计算等方面总结加速策略：**1. 压力方程求解器的选择与优化****2. 算法层面的加速****3. 离散格式与网格优化****4. 并行计算与硬件加速****5. 代码级…...

编程日记 2025/6/20 15:12:33

【C#】 lock 关键字

在 C# 里，lock 关键字就是对 Monitor.Enter/Exit 的简写。它的作用是保证“同一时刻只有一个线程能进入被保护的代码块”，从而避免多个线程同时修改同一个共享状态导致竞态条件（race condition）。一、结合Jog 的例子讲解 // Mot…...

编程日记 2025/5/31 9:18:29

前端脚手架开发指南：提高开发效率的核心操作

前端脚手架通过自动化的方式可以提高开发效率并减少重复工作，而最强大的脚手架并不是现成的那些工具而是属于你自己团队量身定制的脚手架！本篇文章将带你了解脚手架开发的基本技巧，帮助你掌握如何构建适合自己需求的工具，并带着你…...

编程日记 2025/6/2 22:02:05

职坐标AIoT技能培训课程实战解析

职坐标AIoT技能培训课程以人工智能与物联网技术深度融合为核心，构建了“理论实战行业应用”三位一体的教学体系。课程体系覆盖Python编程基础、传感器数据采集、边缘计算开发、云端服务部署及智能硬件开发全链路，通过分层递进的知识模块帮助学员建立系统…...

编程日记 2025/5/31 18:15:25

Yocto Project 快速构建

此文为机器辅助翻译，仅供个人学习使用，如有翻译不当之处欢迎指正 1 Yocto 项目快速构建 1.1 欢迎！ 本简短文档将引导您完成使用 Yocto Project 进行典型镜像构建的流程，并介绍如何为特定硬件配置构建环境。您将使用 Yocto Proj…...

编程日记 2025/5/26 18:29:17

git相关配置

git相关配置欢迎使用Markdown编辑器修改Git默认编辑器为vimgit配置默认用户名和密码： 欢迎使用Markdown编辑器修改Git默认编辑器为vim #方法1：直接执行 git config --global core.editor vim#方法2：修改git的配置文件.git/config文件&am…...

编程日记 2025/5/28 8:09:16

ci/cd全流程实操

本次采用架构，gitlab + jenkins + 镜像仓库+ k8s 准备工作一、gitlab部署拉取镜像部署环境： macbook m2中docker部署gitlab （m2平台架构问题，这里只能用yrzr/gitlab-ce-arm64v8 这个容器镜像） docker pull yrzr/gitlab-ce-arm64v8 在 Docker 里，–privileged=tr…...

编程日记 2025/5/31 7:29:09

Python中in和is关键字详解和使用

在 Python 中，in 和 is 是两个常用但含义不同的关键字，初学者很容易混淆它们的用法。下面是关于它们的详细解释、注意事项及常见示例。一、关键字 in：成员运算符 1. 功能用于判断某个元素是否存在于序列（如列表、元组、字符串…...

编程日记 2025/6/1 17:44:15

ACM模式用Scanner和System.out超时的解决方案和原理

Hi~！这里是奋斗的明志，很荣幸您能阅读我的文章，诚请评论指点，欢迎欢迎 ~~ 🌱🌱个人主页：奋斗的明志 🌱🌱所属专栏：笔试强训 📚本系列文章为个人学…...

编程日记 2025/5/29 20:26:18

微服务中服务降级和异常的区别

在Java中，服务降级和异常处理是两个相关但不同的概念。它们的主要区别如下： 1. 服务降级（Service Degradation）: 定义：服务降级是指在系统中某个服务或功能出现问题时，通过采取某些策略来降低服务的质量或…...

编程日记 2025/5/27 7:34:27

MYSQL创建索引的原则

创建索引的原则包括： 表中的数据量超过10万以上时考虑创建索引。选择查询频繁的字段作为索引，如查询条件、排序字段或分组字段。尽量使用复合索引，覆盖SQL的返回值。如果字段区分度不高，可以将其放在组合索引的后面。对于…...

编程日记 2025/5/27 6:48:12

29、魔法微前端——React 19 模块化架构

一、时空结界分割术（模块化架构设计） 1. 次元切割协议 // 主应用入口const HogwartsMain () > {const [subApps] useState({potion: React.lazy(() > import(./PotionShop)),library: React.lazy(() > import(./LibraryApp)),quidditch: R…...

编程日记 2025/6/1 9:46:49

【PmHub后端篇】PmHub 中缓存与数据库一致性的实现方案及分析

在软件开发项目中，缓存的使用十分普遍。缓存作为一种存储机制，能够暂时保存数据，从而加速数据的读取和访问。然而，当数据同时存在于缓存和数据库中时，如何保证两者的数据一致性成为了一个关键问题。在 PmHub 项目中&am…...

编程日记 2025/5/29 9:43:24

Verilog HDL 语言整理

Verilog HDL 语言 Verilog HDL 简介硬件描述语言Hardware Description Language是一种用形式化方法即文本形式来描述和设计数字电路和数字系统的高级模块化语言 Verilog HDL（Hardware Description Language）是一种硬件描述语言，用于建模…...

编程日记 2025/5/26 10:17:54

[250516] OpenAI 升级 ChatGPT：GPT-4.1 及 Mini 版上线！

目录 ChatGPT 迎来重要更新：GPT-4.1 和 GPT-4.1 mini 正式上线用户如何访问新模型？技术亮点与用户体验优化 ChatGPT 迎来重要更新：GPT-4.1 和 GPT-4.1 mini 正式上线 OpenAI 宣布在 ChatGPT 平台正式推出其最新的 AI 模型 GPT-4.1 和 GPT-4.…...

编程日记 2025/5/29 14:49:53

R语言学习--Day03--数据清洗技巧

在一般情况下，我们都是在数据分析的需求前提下去选择使用R语言。而实际上，数据分析里，百分之八十的工作，都是在数据清洗。并不只是我们平时会提到的异常值处理或者是整合格式，更多会涉及到将各种各样的数据整合&#x…...

编程日记 2025/5/28 8:41:19

文件系统交互实现

关于之前的搭建看QT控件文件系统的实现-CSDN博客，接下来是对本程序的功能完善，我想着是这样设计的，打开一个目录以后，鼠标选中一个项可以是目录，也可以是文件，右键可以出现一个菜单选择操作，比如…...

编程日记 2025/6/1 0:41:00

SqlHelper 实现类，支持多数据库，提供异步操作、自动重试、事务、存储过程、分页、缓存等功能。

/// <summary> /// SqlHelper 实现类，支持多数据库，提供异步操作、自动重试、事务、存储过程、分页、缓存等功能。 /// </summary> public class SqlHelper : IDbHelper {private readonly IDbConnectionFactory _connectionFactory;private…...

编程日记 2025/5/28 14:37:53

DevExpressWinForms-RichEditControl-基础应用

RichEditControl-基础应用在企业级WinForms应用开发中，富文本编辑与文档处理是常见需求。DevExpress WinForms的RichEditControl作为一款功能强大的富文本编辑控件，提供了媲美Microsoft Word的文档处理能力，支持复杂格式编辑、打印导出、界…...

编程日记 2025/6/20 15:11:25

Elasticsearch 索引副本数

作者：来自 Elastic Kofi Bartlett 解释如何配置 number_of_replicas、它的影响以及最佳实践。更多阅读：Elasticsearch 中的一些重要概念: cluster, node, index, document, shards 及 replica 想获得 Elastic 认证？查看下一期 Elasticsearc…...

编程日记 2025/6/20 15:08:45

RabbitMQ 扇形交换器工作原理详解

目录一、扇形交换器简介二、扇形交换器工作原理2.1 消息广播机制2.2 路由键的忽略三、代码示例3.1 生产者代码3.2 消费者代码四、实际应用场景4.1 日志收集系统4.2 实时通知系统4.3 事件驱动架构五、总结在 RabbitMQ 的众多交换器类型中，扇形交换器（Fanout Exchange）是一种…...

编程日记 2025/6/20 15:11:24

IDEA中springboot项目中连接docker

具体内容如下： 1、在Linux中安装docker 使用安装命令： apt-get install docker.io 还有一个是更新软件并安装docker： sudo apt-get update sudo apt-get install docker-ce docker-ce-cli containerd.io 运行docker systemctl start …...

编程日记 2025/6/20 15:09:44

arxiv等开源外文书数据的获取方式

一、一些基本说明开放API接口文档：https://info.arxiv.org/help/api/user-manual.html#2-api-quickstart研究领域分类说明文档：https://arxiv.org/category_taxonomy 二、基于url接口方式检索并获取数据本质是get方式，在url中传检索参数…...

编程日记 2025/6/20 15:08:46

ChatGPT再升级！

近日，OpenAI 正式发布 GPT-4.1 和轻量级版本 GPT-4.1mini，并已全面上线 ChatGPT 平台，迅速引发全球 AI 圈热议，标志着 ChatGPT 在智能化和效率上再登新高峰。 GPT-4.1 是为编程与任务处理优化的高性能模型。相较前作 GPT-4o&#…...

编程日记 2025/5/29 17:21:51

23、电网数据管理与智能分析 - 负载预测模拟 - /能源管理组件/grid-data-smart-analysis

76个工业组件库示例汇总电网数据管理与智能分析组件 1. 组件概述本组件旨在模拟一个城市配电网的运行状态，重点关注数据管理、可视化以及基于模拟数据的智能分析，特别是负载预测功能。用户可以通过界面交互式地探索电网拓扑、查看节点状态、控制时间…...

编程日记 2025/5/27 15:02:43

#跟着若城学鸿蒙# web篇-获取定位

前言在业务中，某些网页上需要获取用户的地理位置，然后按照用户搜索的兴趣点与用户的距离远近进行排序，这就需要h5能够获取到用户的位置。由于 web 组件基于Chromium M114 版本开发，前端就可以使用navigator.geolocation.getC…...

编程日记 2025/5/28 7:28:42

前端批量下载文件打包为zip

多文件需要一次性下载为zip文件这是近期遇到的一个需求，本身是多文件上传的，下载时单个下载太慢又繁杂，用户希望能一次性批量下载，就选择了jszip import axios from "axios" import JSZip from "jszip" im…...

编程日记 2025/6/4 10:16:46

Vue百日学习计划Day9-15天详细计划-Gemini版

重要提示： 番茄时钟： 每个番茄钟为25分钟学习，之后休息5分钟。每完成4个番茄钟，进行一次15-30分钟的长休息。灵活性： JavaScript 的概念较多，尤其是 this、原型链、闭包和异步编程，可能需要更多…...

编程日记 2025/5/24 12:00:03

MySQL8.x新特性：与mysql5.x的版本区别

MySQL8.x新特性 1.与mysql5.x的区别：MySQL8.x新特性：与mysql5.x的版本区别-CSDN博客 2.窗口函数（Window Functions）：MySQL8.x新特性：窗口函数（Window Functions）-CSDN博客引言 …...

编程日记 2025/5/30 17:57:37

RabbitMQ 消息模式实战：从简单队列到复杂路由（三）

精准投递：路由模式路由模式详解路由模式是 RabbitMQ 中一种功能强大且灵活的消息传递模式，它在发布订阅模式的基础上，引入了路由键（Routing Key）的概念，实现了消息的精准路由和分发。在路由模式中&…...

编程日记 2025/6/1 8:08:54

STM32 定时器主从模式配置解析

STM32 定时器主从模式配置解析下面这两行代码是配置STM32定时器主从模式的关键设置代码功能解析 TIM_SelectInputTrigger(TIM3, TIM_TS_TI2FP2); // 选择从模式输出的触发源 TIM_SelectSlaveMode(TIM3, TIM_SlaveMode_Reset); // 选择从模式1. TIM_SelectInputTrigger(T…...

编程日记 2025/5/31 9:18:58

覆盖最小子串代码来自b站左程云 class Solution {public String minWindow(String str, String tar) {char[] s str.toCharArray();char[] t tar.toCharArray();int[] cnt new int[256];for (char cha : t) { cnt[cha]--;}int len Integer.MAX_VALUE;int debt t.length…...

编程日记 2025/5/27 6:33:49

Perl语言深度考查：从文本处理到正则表达式的全面掌握

阅读原文前言：为什么Perl依然值得学习？ "这个脚本用Perl写只需要5分钟！"——在当今Python大行其道的时代，你依然能在不少企业的运维部门听到这样的对话。Perl作为一门有着30多年历史的语言，凭借其强大的文…...

编程日记 2025/6/1 16:30:06

idea中Lombok失效的解决方案

Lombok 是一个 Java 库，旨在通过注解简化 Java 代码的编写，减少样板代码，提高开发效率。它通过自动生成常见的代码（如 getter、setter、构造函数等）来减少开发者的手动编码工作。一般Lombok失效有四步排查方案&#…...

编程日记 2025/5/27 8:03:25

【LeetCode 热题 100】动态规划系列

📁 70. 爬楼梯状态标识：爬到第i层楼梯时，有多少种方法。状态转移方程：dp[i] dp[i-1] dp[i-2]，表示从走一步和走两步的方式。初始化：dp[1] 1 , dp[2] 2。返回值：dp[n]，即走到…...

编程日记 2025/5/29 22:10:18

刷leetcodehot100返航版--双指针5/16

for (int i 0, j 0; i < n; i ) { while (j < i && check(i, j)) j ; // 具体问题的逻辑 } 常见问题分类： (1) 对于一个序列，用两个指针维护一段区间 (2) 对于两个序列，维护某种次序，比如归并排序中…...

编程日记 2025/5/27 6:45:10

DAY24元组和OS模块

元组元组的特点： 有序，可以重复，这一点和列表一样元组中的元素不能修改，这一点非常重要，深度学习场景中很多参数、形状定义好了确保后续不能被修改。很多流行的 ML/DL 库（如 TensorFlow, PyTorch, Num…...

编程日记 2025/5/27 6:33:39

CSS：三大特性

文章目录一、层叠性二、继承性三、优先级一、层叠性二、继承性可以在MDN网站上查看属性是否可以被继承例如color 三、优先级...

编程日记 2025/5/22 8:16:04

Cross-Site Scripting（XSS）

1. XSS介绍跨站脚本攻击（Cross-Site Scripting）简称XSS，人们经常将跨站脚本攻击（Cross Site Scripting）缩写为CSS，但这会与层叠样式表（Cascading Style Sheets，CSS）的缩…...

编程日记 2025/5/27 5:26:51

掌握HTML文件上传：从基础到高级技巧

HTML中input标签的上传文件功能详解一、基础概念 1. 文件上传的基本原理在Web开发中，文件上传是指将本地计算机中的文件（如图片、文档、视频等）传输到服务器的过程。HTML中的<input type"file">标签是实现这一功能的基础…...

编程日记 2025/5/29 9:20:57

WebRTC中的几个Channel

一、我指的是谁？ 以视频为例，常见的有：MediaChannel、VideoMediaChannel、WebRtcVideoChannel、BaseChannel、VideoChannel，那么，为什么要这么多Channel，只写一个叫做SuperChannel行不行（很多程…...

编程日记 2025/5/29 17:02:00

【设计模式】- 行为型模式1

模板方法模式定义了一个操作中的算法骨架，将算法的一些步骤推迟到子类，使得子类可以不改变该算法结构的情况下重定义该算法的某些步骤【主要角色】： 抽象类：给出一个算法的轮廓和骨架（包括一个模板方法和若干基…...

编程日记 2025/5/18 8:45:57

容器化-k8s-使用和部署

一、K8s 使用 1、基本概念集群：由 master 节点和多个 slaver 节点组成，是 K8s 的运行基础。节点：可以是物理机或虚拟机，是 K8s 集群的工作单元，运行容器化应用。Pod： K8s 中最小的部署单元，一个 Pod 可以包含一个或多个紧密相关的容器，这些容器共享网络和存储资源。…...

编程日记 2025/5/28 7:56:32

黑马k8s（九）

1.Pod-生命周期概述 2.Pod生命周期-创建和终止 3.Pod生命周期-初始化容器...

编程日记 2025/5/28 14:47:32

Android trace中CPU的RenderThread与GPU

Android trace中CPU的RenderThread与GPU RenderThread是系统的GPU绘制线程，GPU渲染就是通常所谓的硬件加速，如果应用关闭硬件加速，就没有了RenderThread，只有UI Thread，即Android主线程。 Android GPU渲染SurfaceFlin…...

编程日记 2025/5/18 10:17:14

测试工程师如何学会Kubernetes(k8s)容器知识

Kubernetes(K8s)作为云原生时代的关键技术之一，对于运维工程师、开发工程师以及测试工程师来说，都是一门需要掌握的重要技术。作为一名软件测试工程师，学习Kubernetes是一个有助于提升自动化测试、容器化测试以及云原生应用测试能力的重要过程…...

编程日记 2025/5/29 9:51:45

接触感知钳位电路分析

以下是NG板接触感知电路的原理图。两极分别为P3和P4S，电压值P4S < P3。电路结构分两部分，第一部分对输入电压进行分压钳位。后级电路使用LM113比较器芯片进行电压比较，输出ST接触感知信号。钳位电路输出特性分析输出电压变化趋势&a…...

编程日记 2025/5/18 9:08:21

码蹄集——圆包含

MT1181 圆包含输入2个圆的圆心的坐标值（x，y）和半径，判断断一个圆是否完全包含另一个圆，输出YES或者NO。另：内切不算做完全包含。格式输入格式：输入整型，空格分隔。每行输入一组…...

编程日记 2025/6/1 12:40:17

ConcurrentSkipListMap的深入学习

目录 1、介绍 1.1、线程安全 1.2、有序性 1.3、跳表数据结构 1.4、API 提供的功能 1.5、高效性 1.6、应用场景 2、数据结构 2.1、跳表（Skip List） 2.2、节点类型： 1.Node 2.Index 3.HeadIndex 2.3、特点 3、选择层级 3.1、随…...

编程日记 2025/5/30 0:08:41

ProfibusDP主站转modbusTCP网关接DP从站网关通讯案例

ProfibusDP主站转modbusTCP网关接DP从站网关通讯案例在工业自动化领域，Profibus DP和Modbus TCP是两种常见的通信协议。Profibus DP广泛应用于过程自动化、工厂自动化等场景，而Modbus TCP则常见于楼宇自动化、能源管理等领域。由于设备和系统之间往往存…...

编程日记 2025/5/22 8:21:53