Q144：FS，求解流体方程（逻辑总结）

最新推荐文章于 2021-07-13 11:34:12 发布
原创 最新推荐文章于 2021-07-13 11:34:12 发布 · 589 阅读 · 1 ·
CC 4.0 BY-SA版权
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
文章标签：

#流体 #模拟 #求解流体方程 #逻辑

本文针对《FluidSimulationforComputerGraphics》的基础部分进行了学习总结，主要从非技术角度阐述了流体方程求解的基本逻辑，包括精简版和完整版两个部分。

看完了《Fluid Simulation for Computer Graphics》的第一部分（基础部分）。
真的只能算是“看完”，着实搞不懂其中的数学原理and物理原理。

本文简单理清一下“求解流体方程的逻辑”（不含数学原理、物理原理）
分两部分：其一、精简逻辑版本；其二，笔记完整版本。

其一、精简逻辑版本

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述

其二，笔记完整版本

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

这里写图片描述这里写图片描述

NS方程求解-PointNet和升维思想（效果很差）
forrestguang的博客
04-14 3176
PointNet求解NS方程 参考文献：Kashefi, A. , and T. Mukerji . “Physics-Informed PointNet: A Deep Learning Solver for Steady-State Incompressible Flows and Thermal Fields on Multiple Sets of Irregular Geometries.” (2022). PointNet网络的特点在于中间通过pool提取过一次全局特征，这个特点导致一个问题。
流体动力学控制方程（详细推导）
diyhoo的编程小屋
07-06 4万+
流体流动要受到物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律。控制方程（governing equations）是这些守恒定律的数学描述。本篇文章将对三大流体力学方程进行推导，以加深对其理解。
计算流体力学三大基本方程：连续性方程、动量方程（纳维-斯托克斯方程）和能量方程推导
热门推荐
BertramZhang的博客
01-31 28万+
方法取自安德森《计算流体力学基础》 连续性方程 固定流体微元内质量变化率=流体从笛卡尔坐标三个方向流出量 因此可得： 质量变化率: 则： 连续性方程： 用散度表示则可得到： 对于不可压缩流体，其密度为一常数，因此可以得到： 动量方程（纳维-斯托克斯方程） 根据牛顿第二定律可以得出：F=ma；
定点运算
Rihaong' Blog
09-05 3103
定点运算： 移位 加、减 乘 除 一、移位 1.算数移位规则 【不同码制机器数移位后的空位添补规则】 【Example】 2.算数移位和逻辑移位 （1）算数移位：有符号数的移位 （2）逻辑移位：无符号数的移位 二、加法 、减法 1.补码加法、减法运算基本公式 可以把符号位和数位同等处理。 无论操作数正负，在做补码加减法时，只需将符号位和数值部分一起参加运算...
CFD求解NS方程代码
12-26
求解NS方程会碰到一个问题就是CFL问题，该代码研究时间参数设置对结果影响
[物理学与PDEs]第2章习题10 一维理想流体力学方程组的 Lagrange 形式
weixin_34290631的博客
04-12 186
试证明: 一维理想流体力学方程组的 Lagrange 形式 (5. 22)-(5. 24) 也可写成如下形式 $$\beex \bea \cfrac{\p \tau}{\p t}-\cfrac{\p u}{\p x}&=0,\\ \cfrac{\p u}{\p t}+\cfrac{\p p}{\p x}&=F,\\ \cfrac{\p }{\p t}\sex{e+\cfrac{u^...
流体力学基础学习
Mr.Tan的博客
09-06 5859
前言 笔者因为需要稍微看了一下流体力学相关的知识，是b站中一个技术邻的up主传的视频（流体力学基础）。老师讲的很好，课件做的很用心。这里是笔者在学习中觉得重要的一些截图，如有侵权马上删除。由于笔者在学习过程中，b站的资源内容出现了点问题，暂时就先到这里。后面会补充 图片如下 01 02 03 04 05 06 07 08 09 10 ...
【MATLAB自适应网格】：Gray–Scott方程求解精度的提升策略
[【MATLAB自适应网格】：Gray–Scott方程求解精度的提升策略](https://media.springernature.com/lw1200/springer-static/image/art%3A10.1007%2Fs00211-022-01296-x/MediaObjects/211_2022_1296_Fig9_HTML.png) ...
【MATLAB数值分析】：Gray–Scott方程求解的实战技巧与案例剖析
：Gray–Scott方程求解的实战技巧与案例剖析](https://media.springernature.com/lw1200/springer-static/image/art%3A10.1038%2Fs41598-023-39741-y/MediaObjects/41598_2023_39741_Fig1_HTML.png) # 1. Gray–...
【Broyden方法全解析】：提升非线性方程组求解效率的7大实用技巧
本文全面介绍了Broyden方法，一种用于求解非线性方程组的有效迭代算法。文章首先概述了Broyden方法的概念和理论基础，包括其定义、数学原理以及与其他数值方法的比较。接着，文章深入探讨了该方法在实践中的技巧，...
Matlab神功揭秘：非线性方程求解的终极指南（2023版）
[Matlab神功揭秘：非线性方程求解的终极指南（2023版）](https://media.springernature.com/lw1200/springer-static/image/art%3A10.1007%2Fs40305-023-00453-9/MediaObjects/40305_2023_453_Figa_HTML.png) ...
高斯展开法求解薛定谔方程：从入门到精通（完整学习路径）
[高斯展开法求解薛定谔方程：从入门到精通（完整学习路径）](https://media.springernature.com/lw1200/springer-static/image/art%3A10.1007%2Fs42967-023-00250-4/MediaObjects/42967_2023_250_Fig18_HTML.png) ...
因子分析（FA）算法简述
iceberg7012的博客
10-16 7万+
文章目录前言一、什么是因子分析？1.1 因子分析应用背景1.2 因子分析算法的基本步骤1.3 因子分析算法的数学解释1.3.1 因子模型1.3.2 因子载荷矩阵的求解二、因子分析的应用实例三、主成分分析（PCA）与因子分析（FA）的联系与区别总结 前言 在学习数据降维时，了解到因子分析（FA）算法是其中的一种方式，因此，在这里对因子分析算法做一个简要的归纳、梳理，后续会对数据降维的几种方式做个总结，感兴趣的朋友，可以持续关注。 一、什么是因子分析？ 因子分析法是指： 研究从变量群中提取共性因子的统计技术，
流体力学控制方程的总结
diyhoo的编程小屋
07-07 1万+
一、粘性流动的纳维-斯托克斯（Navier-Stokes）方程       粘性流动是包括摩擦、热传导和质量扩散等输运现象的流动，这些输运现象是耗散性的，它们使流体的熵增加。但是质量扩散只有当流动中不同化学组分之间存在浓度梯度时才发生。目前没有考虑质量扩散。1.连续性方程守恒形式：2.动量方程守恒形式：X方向：Y方向：Z方向：3.能量方程守恒形式：二、无粘流欧拉（Euler）方程无粘流的定义是忽略...
简谐振动基本方程的推导
the193thDoctor的博客
02-01 3万+
——本博客为博主原创，如有雷同，纯属巧合，未经允许，不得抄袭 简谐振动，作为自然界极为常见的运动方式，具有简洁而数学化的美妙，而对于简谐运动的公式：      的推导过程，一直是比较复杂的一个问题。博主通过自己的思考，得出了一种比较简洁的方法。 首先，我们可以构建这样一个模型（如下图所示），在一根无限长的数轴上，有一个运动的质点。若t为运动时间，当t=0时，质点具有初速度v0，不妨假设速度的方
计算不可压缩流体- NS方程求解算法
lizhengjiang的专栏
03-20 1万+
在上一篇介绍差分格式可以看到求解不可压缩NS方程，主要是求解速度的动量方程和压力的波松方程。本文介绍一些历史中的算法，下一文介绍目前常用的算法。 一  Marker and Celler （MAC) 算法 在自然差分格式中给出PPE如下：    laplace ( p) = laplace ( body_force) - rho * div ( u * grad_u)    （1） 。
流体力学方程的秘密（1）
無名黑洞
01-12 2301
NS方程组没有简单的求解方法。一个专门的学科叫做计算流体力学就是寻求使用计算机来解决此类问题的。困难的原因在于NS方程是非线性微分方程。比如前一节中的移流项就是非线性项（它意味着速度的变化不仅跟时间有关，还和空间有关，所以它的解无法构成线性空间）。对于非线性微分方程没有通用解。在某些特例的情况下，NS方程组可以化简为线性微分方程组，这个是比较好解的。比如斯托克斯蠕动流，这种流体假定了很低的雷诺数，流体大致上只受到粘性系数影响，而忽略了空间的影响。
超短强脉冲激光与材料相互作用双温模型的研究
midiio的博客
07-13 2398
激光与材料相互作用的问题一直受到人们的青睐。目前已经建立起各种各样的理论模型(如热传导模型、流体动力学模型，分子动力学模型、双温模型等) 大多数只是在特定的范围内有用。 经典的热传导模型主要用在长脉冲烧蚀的领域( 如纳秒级脉冲激光) ，当激光脉冲达到飞秒级时，由于飞秒激光烧蚀靶材过程中电子和晶格之间能量的转移时间短于脉冲激光辐照靶材的时间，烧蚀机理则有明显不同。对于分析超短脉冲激光烧蚀过程中电子和离子亚系统之间相互作用需采用双温模型。双温模型在激光和金属材料作用的领域里，取得了许多成果，得到了广泛的应..
亥姆霍兹流体谐振器MATLAB
最新发布
04-28
### MATLAB 中亥姆霍兹流体谐振器的模拟与实现 #### 背景介绍 亥姆霍兹谐振器是一种基于声学原理设计的装置，其共振频率可以通过几何参数计算得出。在 MATLAB 中，可以利用数值方法来模拟亥姆霍兹谐振器的行为，包括求解共振频率以及分析压力场分布。 #### 数学模型 亥姆霍兹谐振器的核心方程为： \[ f_r = \frac{c}{2\pi} \sqrt{\frac{A}{V L}} \] 其中： - \( f_r \) 是共振频率， - \( c \) 是声音的速度[^1]， - \( A \) 是颈部横截面积， - \( V \) 是腔体体积， - \( L \) 是颈部长度。 通过该公式可计算理论上的共振频率，并将其用于验证仿真结果。 #### MATLAB 实现代码示例 以下是 MATLAB 的代码示例，展示如何模拟亥姆霍兹谐振器并绘制相关图形： ```matlab % 参数定义 c = 343; % 声音速度 (m/s) A = 0.01^2 * pi; % 颈部横截面积 (m^2)，假设直径为 1 cm 圆形开口 L = 0.05; % 颈部长度 (m) V = 0.001; % 腔体体积 (m^3) % 计算共振频率 fr = (c / (2*pi)) * sqrt(A / (V*L)); disp(['亥姆霍兹谐振器的共振频率为: ', num2str(fr), ' Hz']); % 时间域信号生成 fs = 10000; % 采样率 (Hz) t = 0:1/fs:1; % 时间向量 (秒) f_signal = fr; % 输入信号频率等于共振频率 signal = sin(2*pi*f_signal*t); % 正弦波输入信号 % 绘制时间域信号 figure; plot(t, signal); title('输入正弦波信号'); xlabel('时间 (s)'); ylabel('幅值'); % FFT 分析以观察频谱特性 NFFT = length(signal); % FFT 点数 Y = fft(signal, NFFT)/NFFT; frequencies = fs/2*linspace(0,1,NFFT/2+1); % 取单边幅度谱 magnitude = 2*abs(Y(1:NFFT/2+1)); figure; plot(frequencies, magnitude); title('频谱图'); xlabel('频率 (Hz)'); ylabel('幅值'); grid on; % 找到峰值频率 [~, idx_max] = max(magnitude); peak_frequency = frequencies(idx_max); disp(['检测到的主要频率成分: ', num2str(peak_frequency), ' Hz']); ``` #### 结果解释 上述代码实现了以下几个功能： 1. **共振频率计算**：根据给定的物理参数，计算出亥姆霍兹谐振器的理论共振频率。 2. **时间域信号生成**：创建了一个与共振频率一致的正弦波作为输入信号。 3. **频谱分析**：通过对输入信号进行快速傅里叶变换（FFT），展示了频域中的主要频率分量，并确认其接近于理论共振频率。 #### 进一步扩展 为了更全面地研究亥姆霍兹谐振器的行为，还可以引入有限元法（FEM）或边界元法（BEM）。这些高级技术能够处理复杂的几何形状和材料属性，从而提供更加精确的结果。 ---
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值