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

embedded/2024/9/19 21:15:20/ 标签: 笔记

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 提出新格式,也没有从理论上证明收敛精度,就是做了数值实验证明了收敛精度而已

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

On the Immersed Boundary Method with Time-Filter-SAV for Solving Fluid–Structure Interaction Problem

https://doi.org/10.1007/s10915-024-02591-5

SAV

在这里插入图片描述

在对流项前面乘一个 1

转成 q/q 的形式,然后算 dq/dt

然后这个 dq/dt 就可以离散化

之后的证明部分完全看不懂……

ANALYSIS OF THE PARAREAL TIME-PARALLEL TIME-INTEGRATION METHOD

https://doi.org/10.1137/05064607X

感觉这其中提到的

multiple shooting method

multigrid-in-time method

parabolic multigrid method

space-time multigrid

time-parallel multigrid

multigrid waveform relaxation methods

都是数值分析的方法啊……不知道怎么直接应用到 NS


http://www.ppmy.cn/embedded/113928.html

相关文章

高级算法设计与分析 学习笔记6 B树

B树定义 一个块里面存了1000个数和1001个指针,指针指向的那个块里面的数据大小介于指针旁边的两个数之间 标准定义: B树上的操作 查找B树 创建B树 分割节点 都是选择正中间的那个,以免一直分裂。 插入数字 在插入的路上就会检查节点需不需要…

Testbench编写与Vivado Simulator的基本操作

Testbench编写与Vivado Simulator的基本操作 Testbench编写 Testbench 是一种用Verilog或者systemVerilog语言编写的程序或模块,编写testbench的主要目的是为了对使用硬件描述语言(HDL)设计的电路UUT(unit under test)进行仿真验证&#xf…

一招教你解决excel表格打印预览时候表格线条显示不全的问题

1、如图,我们在制作好excel表格后再需要打印时候,点击打印预览会出现以下情况: 最下边的表格线条显示不全,这样即使打印出来或者导出为pdf,文件中依然显示不全,这时候我们只需要在excel表格中轻轻设置一下就…

CleanMyMac X 4.15.6正式版 mac直装破解版

你知道 CleanMyMac是什么吗?它的字面意思为“清理我的Mac”,作为软件,那就是一款 Mac清理工具 ,Mac OS X 系统下知名系统清理软件,是数以万计的Mac用户的选择。它可以流畅地与系统性能相结合,只需…

JVM 运行时数据区域

目录 前言 程序计数器 java虚拟机栈 本地方法栈 java堆 方法区 运行时常量池 前言 首先, java程序在被加载在内存中运行的时候, 会把他自己管理的内存划分为若干个不同的数据区域, 就比如你是一个你是一个快递员, 一堆快递过来需要你分拣, 这个时候, 你就需要根据投放的目…

三、(JS)JS中常见的表单事件

一、onfocus、onblur事件 这个很容易理解&#xff0c;就不解释啦。 <!DOCTYPE html> <html lang"en"><head><meta charset"UTF-8"><meta name"viewport" content"widthdevice-width, initial-scale1.0"&…

关于 Goroutines 和并发控制的 Golang 难题

下面是一道关于 Goroutines 和并发控制的 Golang 难题&#xff0c;它涉及到 Go 的并发编程模型、Goroutines、通道&#xff08;Channels&#xff09;以及 sync.WaitGroup 的使用&#xff1a; 问题描述&#xff1a; 你有一个需要并发执行的任务&#xff0c;其中有 100 个 URL …

基于Spring Boot的学生社区故障维修预约系统的设计与实现(开题报告)

毕业论文(设计)开题报告 基于Spring Boot的学生社区故障维修预约系统设计与实现 姓 名 学 院 数学与数据科学学院 专业班级 信息与计算科学202 学 号 202021314223 校内指导教师 职称/职务 副教授 校外指导教师 职称/职务 技术经理 起始时间 2023年9月 教务部制 一、开…

【LLM多模态】文生视频评测基准VBench

note VBench的16个维度自动化评估指标代码实践&#xff08;待完成&#xff09;16个维度的prompt举例人类偏好标注&#xff1a;计算VBench评估结果与人类偏好之间的相关性、用于DPO微调 文章目录 note一、相关背景二、VBench评测基准概述&#xff1a;论文如何解决这个问题&…

Go语言现代web开发13 方法和接口

方法 As you probably noticed, there are no classes in the Go programming language. But we can mimic this by declaring functions on types. The type which declares functions is called the receiver argument and the function declared on the type is called the…

JS基础之【对象详解 -- 对象的属性与方法、遍历对象与内置对象】

&#x1f680; 个人简介&#xff1a;某大型国企高级前端开发工程师&#xff0c;7年研发经验&#xff0c;信息系统项目管理师、CSDN优质创作者、阿里云专家博主&#xff0c;华为云云享专家&#xff0c;分享前端后端相关技术与工作常见问题~ &#x1f49f; 作 者&#xff1a;码…

WPF 的TreeView的TreeViewItem下动态生成TreeViewItem

树形结构仅部分需要动态生成TreeViewItem的可以参考本文。 xaml页面 <TreeView MinWidth"220" ><TreeViewItem Header"功能列表" ItemsSource"{Binding Functions}"><TreeViewItem.ItemTemplate><HierarchicalDataTempla…

【设计模式】UML类图和六大设计原则

前言 在实践中经常看到工厂模式、观察者模式等字眼&#xff0c;渐觉设计模式的重要性&#xff0c;于是开刷设计模式。 本文讲述了UML类图以及设计模式的六大原则 参考资料&#xff1a; 课程视频&#xff1a;黑马程序员Java设计模式 一、UML类图 1. 类和接口的表示方式 如…

JavaScript在数据可视化领域的探索与实践

目录 引言 JavaScript可视化库概览 D3.js基础入门 1. 引入D3.js 2. 绘制简单的条形图 3. 添加轴 交互性与动画 实际应用场景 结论 引言 在数据驱动决策日益重要的今天&#xff0c;数据可视化成为连接数据与洞察的桥梁。JavaScript&#xff0c;作为前端开发的主力军&am…

Python识别文字中的省市区

一、库安装 pip install cpca 注意&#xff1a;目前 cpca 模块仅支持Python3及以上版本。 二、基本使用 通过两行代码就能实现最基本的省市区提取&#xff1a; import cpcalocation_str ["新疆古阿贾克斯就打开房间啊开始","河北省石家庄市动物园",&…

CSP-J 计算机网络

文章目录 前言计算机网络的定义计算机网络的发展计算机网络的主要功能计算机网络的分类按网络地理范围分类按网络拓扑结构分类 OSI模型与TCP/IP模型OSI模型TCP/IP模型OSI模型与TCP/IP模型的网络协议及功能 IP地址域名1. **通用顶级域名&#xff08;gTLD&#xff0c;Generic Top…

掌握RESTful API设计:构建高效、可扩展的Web服务

设计RESTful API时&#xff0c;遵循最佳实践可以提高API的可维护性、可扩展性和用户体验。以下是一些关键的最佳实践&#xff1a; 使用HTTP方法正确&#xff1a;GET用于获取资源。POST用于创建新资源。PUT或PATCH用于更新现有资源。DELETE用于删除资源。无状态&#xff1a;每个…

iptables限制网速

1、使用hashlimit来限速 #从eth0网卡进入INPUT链数据&#xff0c;使用模块hashlimit 限制网速为100kb/s或2mb/s,超过限制的数据包会被DROP。OUTPUT链同理&#xff0c;mode为srcip&#xff0c;有4个mode选项: srcip&#xff08;默认匹配每个源地址IP&#xff0c;配置指定源地址…

如何使用命令安装android的.aab包

在Android设备上安装 .aab (Android App Bundle) 文件不是直接使用标准命令行工具的过程&#xff0c;通常需要通过Google Play或其他发布平台来安装它们。但是&#xff0c;如果你想在本地测试 .aab 文件&#xff0c;可以通过以下步骤进行操作&#xff1a; 方法 1&#xff1a;使…

【电脑组装】✈️从配置拼装到安装系统组装自己的台式电脑

目录 &#x1f378;前言 &#x1f37b;一、台式电脑基本组成 &#x1f37a;二、组装 &#x1f379;三、安装系统 &#x1f44b;四、系统设置 &#x1f440;五、章末 &#x1f378;前言 小伙伴们大家好&#xff0c;上篇文章分享了在平时开发的时候遇到的一种项目整合情况&…