浸没边界法空间精度相关的论文的阅读笔记

于 2024-09-15 19:32:36 发布
阅读量723 收藏 10
点赞数 6
文章标签： 笔记
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hijackedbycsdn/article/details/142275926
版权

Convergence proof of the velocity field for a stokes flow immersed boundary method

https://doi.org/10.1002/cpa.20233

研究对象的选取

在这里插入图片描述

他这里为什么能够选取一个周期性边界的流场啊？为什么不是狄利克雷边界或者诺伊曼边界？

方形流场的边界值

为什么要定义方形流场的边界值为 -pi 到 pi 而不是别的？

浸没边界上的坐标系的定义方式

在这里插入图片描述

这个浸没边界上的坐标系的定义方式，这个记号，之前我也没见过

周期性边界

在这里插入图片描述

他定义周期性边界或许是为了在这个公式中能够得到边界上的积分为 0？

还是说我想错了，并不是先希望它积分为 0 再定的周期性边界条件，而是先有定义周期性边界条件的惯例，然后再自然地得出这些量的积分为 0？

在这里插入图片描述

然后就是这个 delta 函数的利用很灵性

就是，它本身是一个域内的 delta，那么就可以用对这个域的积分来消去它

就像是解引用一样，获得了那个 delta 修饰的，限定了范围的那个地方上的值

在这里插入图片描述

但是我不知道为什么它想得到这个式子

或许他在写 NS 方程的时候，在写 f-g 的时候，就是故意要写一个常数 g，为了之后可以说明 f 中不包含常数项？

或者说这是确保自己的公式中的某个物理量不含常数项，所必要的手段？就是直接把这个物理量 u 拆成 u* + u0 认为 u0 是常数，然后经过公式推导得到 u0 的解析式或者怎么滴

因为 u0 是常数，所以不用害怕推导，性质良好，都可以从积分中提出来，所以肯定方便推导？

速度约束

在这里插入图片描述

凭什么给速度添加这个约束啊，我搞不懂，你怎么知道自己可以这样设置约束呢？万一这个约束并不物理呢？

离散格式

在这里插入图片描述

为什么 N h Nh Nh 和 M Δ θ M \Delta \theta MΔθ 会都是 2 π 2\pi 2π 啊，我不理解，这是直接设的随便一个值，还是说他会有物理意义？

我感觉是有物理意义，但是我理解不了，为什么是 2 π 2\pi 2π 比如如果你说是周长的话，那么正方形的周长不应该是 4 π 4\pi 4π 吗

然后浸没边界的坐标为什么也是 2 π 2\pi 2π 跟方形边界类似？为啥啊

后面的格林函数法似乎是泛函分析相关的文章

Accuracy analysis of immersed boundary method using method of manufactured solutions

https://doi.org/10.1007/s10483-010-1353-x

Method of Manufactured Solutions 制造解法

在这里插入图片描述

这个我还是第一次听说

原来做解析解如此简单，只需要自己随便写一个解析解然后凑源项就好了

把源项设置好了，我们就可以认为带了这个源项的 NS 的解就是你用来推这个源项的解析解

就，暂时认为 NS 方程的解是具有唯一性的

代入压强梯度到 NS

在这里插入图片描述

挺神奇的……代入一个阶跃函数之后，居然可以直接转化成一个对 delta 修饰之后的项求积分

我知道这应该是等价的，就是，从物理上，都是对某个区域

但是就是很神奇

这个神奇的另外一个侧面就是，明明是阶跃的压强，但是代入 NS 之后就直接转换为了对 delta 修饰的 F 的求积分

就，这两个物理量之间是怎么转化的他也不说

that is all

好吧，这个文章不咋滴，2010 年的文章仅仅是把前人的方法验证一遍精度就没有了

Stability And Convergence Of Immersed Boundary Computations

https://citeseerx.ist.psu.edu/document?repid=rep1&type=pdf&doi=8c6d600532e8ac22d3d4a48bad1a8012293049b8

在这里插入图片描述

把速度做成一个，速度随着时间衰减，随着 x 改变相位，的 exp 乘以一个只与 y 有关的速度 u(y)

看了好久才知道单纯在讨论纤维

似乎是单纯判断纤维的稳定性

On the order of accuracy of the immersed boundary method: Higher order convergence rates for sufficiently smooth problems

https://doi.org/10.1016/j.jcp.2005.02.011

混合投影方法

在这里插入图片描述

居然可以投影到一个近似无散的场

居然可以多次投影，称之为混合投影

不用算平方根的 delta 函数

在这里插入图片描述

不用算平方根的 delta 函数可能更快，好合理

时滞的压力梯度

在这里插入图片描述

近似的压力投影

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

他用混合投影是因为可以消除压力振荡……

emmmm 我没做过弹性体，所以不知道这种缺点，所以攻击都不知道要攻击，坏了……

评估精度的用例

在这里插入图片描述

我不知道我要不要像他们那样做一个薄壳……

或者说其实我这个圆柱或者圆球也就是薄壳而已

我也不知道我要不要做弹性的 IB 固体来验证，应该不用吧

etc

人家 Charles S. Peskin 在 JCP 提出新格式，也没有从理论上证明收敛精度，就是做了数值实验证明了收敛精度而已

好吧，这些都是空间精度的，我想看时间精度的

hijackedbycsdn
关注
A Survey on Metaverse: Fundamentals, Security, and Privacy
c_cpp_csharp的专栏
11-08 1343
元宇宙作为下一代互联网的一个不断发展的范例，旨在为人类建立一个完全沉浸式、超时空和自我维持的虚拟共享空间，供人们玩耍、工作和社交。在扩展现实、人工智能和区块链等新兴技术的最新进展推动下，元宇宙正从科幻小说走向即将到来的现实。然而，元宇宙的严重隐私侵犯和安全漏洞（继承自底层技术或在新的数字生态中出现）可能会阻碍其广泛部署。同时，由于元宇宙的固有特性，如沉浸式现实主义、超时空性、可持续性和异质性，在元宇宙安全供应中可能会出现一系列基本挑战（例如，可扩展性和互操作性）。
字典查询（部分）
myth_HG的专栏
05-28 14万+
#include #include #include using namespace std; const int enMax=46; const int chMax=256; const int attrMax = 66; const int spaceMax = 100; const int inputMax = 30; char atrrSum[13][attrMax]= {"n.","
用于模拟颗粒流的直接强迫浸没边界法 An immersed boundary method with direct forcing 笔记
hijackedbycsdn的博客
02-26 1180
使用固定且均匀的计算网格来计算悬浮刚性颗粒周围的不可压缩粘性流主要思想是将 Peskin 的正则化 delta 函数方法 [Acta Numerica 11 (2002) 1] 纳入流固相互作用力的直接公式中，以便允许欧拉和拉格朗日表示之间的平滑转换，同时避免时间步长的严格限制。该技术是在有限差分和分数步环境中实现的。我们提出了一种改进的浸入边界方法，直接公式化了流固相互作用力。Peskin [23] 的正则化 delta 函数用于任意拉格朗日位置和离散欧拉位置​​之间的关联。
用C语言写一个查单词的小demo
m0_47068593的博客
07-03 4万+
用C语言写一个查单词的小demo #include <stdio.h> #include <string.h> #include <stdlib.h> int main(int argc, const char *argv[]) { // 用只读的方式打开我们的单词库 FILE *fp = fopen("dict.txt","r"); if (NULL == fp){ perror("fopen"); return -1
网络编程_8(项目附件)
weixin_38717634的博客
01-13 3万+
dict.txt abandonment n.放弃 abbreviation n.缩写 abeyance n.缓办，中止 abide v.遵守 ability n.能力 able adj.有能力的，能干的 abnormal adj.反常的，变态的 aboard adv.船(车)上 abolish v.废除，取消 abolition n.废除，取消 abortion n.流产 abortive adj.无效果的，失败的 about prep.关于，大约 ab
简单英译汉SQL脚本
热门推荐
shencomd的专栏
10-19 19万+
create table t_dictionary(word varchar(200), trans varchar(1000)); insert into t_dictionary values('abandon','v.抛弃，放弃'); insert into t_dictionary values('abandonment','n.放弃'); insert into t_dictionary
一种基于过滤与反卷积的新型高阶浸没边界法 (2013年)
06-01
传统浸没边界法在边界附近只有一阶精度,而高精度的改进方法都需要额外引入跳跃条件,因此不具备普适性.文中设计了一种基于过滤和反卷积的新型算法,既在一定程度上提高了精度,又避免了以往方法中引入额外跳跃条件的...
一种稳定的浸没边界-格子波尔兹曼方法 (2013年)
05-13
基于强制边界条件的浸没边界-格子波尔兹曼方法(Immersed boundary-lattice Boltz mann method , IB- LBM)提出了一种稳定的IB- LBM。在该方法中,采用了具有良好稳定性的多松弛时间(Multi-relaxation time , MRT )模型...
浸入边界法及其应用
03-10
浸入边界法是一种数值计算方法，主要用于模拟粘性不可压缩流场中浸入边界的运动问题。这种方法特别适用于处理几何外形复杂的弹性膜结构在流场中的动态问题，例如在生物流体动力学中模拟血液流动与心脏活动等。该方法...
WENO格式与虚拟单元浸入边界法在笛卡尔网格中的应用 (2013年)
05-26
高精度有限差分WENO格式在结构网格上处理具有复杂几何外形绕流问题时较困难，而虚拟单元浸入边界法却是一种较新颖且对网格的要求较低的方法，适用于复杂几何外形边界的处理.为此，在笛卡尔网格上采用WENO格式以求解...
STM32学习笔记（二、初识stm32单片机）
tommorrowwill的博客
09-10 777
首先stm32是意法半导体公司（ST）使用ARM公司的Cortex-M为核心生产的32位的单片机。其中，ST---意法半导体公司，即SOC厂商。M---为Microelectronics的缩写，即微型处理器。32---表示控制器为32位的。103---表示F系列的子系列。
C++学习笔记----6、内存管理（五）---- 智能指针（3）
weixin_71738303的博客
09-11 1152
与指向特定类型的原始指针可以转化为不同类型的指针一样，shared_ptr保存特定的类型可以转化为一个另一种类型的shared_ptr。转化shared_ptr的函数是const_pointer_cast()，dynamic_pointer_cast(),static_pointer_cast()和reinterpret_pointer_cast()。前面简要提过，当拥有共享属主的智能指针，例如shared_ptr不在活动范围或者被重置时，只有它是最后指向的智能指针时才能释放其指向的资源。
UEFI学习笔记（五）：EDK II PCD的概念、类型、使用
qq_44189622的博客
09-11 824
如果在BIOS里面有一些模块是binary方式集成进来的而这些binary又需要用到PCD（用于Binary Release），那么这些Binary集成的要用到的PCD就必须要设置为。PCD的值存在memory里面，下次启动时，上次更改的值丢失了，每次启动都是从default值开始。是存在VPD空间的（在FLASH上，只读），一般是出厂配置。如果platform是从源码build出来的，没有binary在里面的时候，PCD用的都是。作用域在一个模块中（模块级的），可以在Binary Level进行修改。
OJ在线判题系统项目笔记
今天不coding的博客
09-11 768
OJ在线判题系统项目笔记
经验笔记：Web 浏览器存储手段
qq_45831414的博客
09-10 874
以上就是使用 IndexedDB 的基本步骤。IndexedDB 提供了强大的数据存储能力，并且可以通过索引来优化数据检索性能。在实际应用中，可能还需要处理更复杂的场景，比如数据迁移、错误处理等。务必确保你的代码能够妥善处理各种异常情况。
Netty笔记03-组件Channel
weixin_46425661的博客
09-12 715
本文主要讲解Netty中的组件Channel、ChannelFuture和CloseFuture
rk3568系统buildroot开发笔记
Reasonss的博客
09-11 1240
构建RK3568系统期间遇到的各种问题,我主要分为编译问题和运行问题
17、电科院FTU检测标准学习笔记-录波性能
最新发布
sixtome
09-14 524
作者简介：本人从事电力系统多年，岗位包含研发，测试，工程等，具有丰富的经验在配电自动化验收测试以及电科院测试中，本人全程参与，积累了不少现场的经验。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值