关注
分享
陆嵩
中国科学院数学与系统科学研究院 计算数学与科学工程计算研究所 科学与工程计算国家重点实验室
State Key Laboratory of Scientific and Engineering Computing,
Institute of Computational Mathematics and Scientific/Engineering Computing,
Academy of Mathematics and Systems Science,
Chinese Academy of Sciences
Email: lusong@lsec.cc.ac.cn
展开
-
关于向量值函数方程变分的一点注记
一般来说，我们喜欢把梯度定义为列向量，一列一列地来，但是上述的关于向量值函数的梯度，却是对每个分量求梯度，然后按行排的，这是反直觉的，所以导致了这个地方的。有人说，自己做的都是标量解函数的偏微分方程变分，当遇到解函数是向量值函数时，就不知道有限元空间怎么取，变分也不知道怎么做了。再 一个对于矩阵的散度和梯度算子，是支持分块矩阵的操作的。有了以上的准备，我们举个简单的例子来看向量值解的方程的变分，头上加转置变回原来的定义方式，就得到了最开始的上述的表达。的地方，是可以随意进出梯度和散度算子的，即。...原创 2022-08-03 04:51:48 · 947 阅读 · 2 评论 -
模型降阶方法之 DEIM
为什么叫插值，为什么叫greedy，可以看得出来，我们选的就是投影之后，距离最大的一个点，然后把它固定。也就是说，我们对于$\mathbf{F}$选行，我们只要计算里面必要分量作用即可，也就是可以把选行操作放到$\mathbf{F}$里面去。POD方法固有的局限性，使得无法很好地处理非线性项，所以必须有一个合适的方法处理非线性项。这里看起来，似乎还需要把非线性项施加于每个分量算出来，其实并不然，这里的。同理，对于稳态相关的问题，我们需要计算其雅克比，非线性项的雅克比矩阵为，...原创 2022-07-24 02:52:18 · 3160 阅读 · 4 评论 -
模型降阶方法之 POD
我们想要找到解所在的超平面，那么事实上，我们只要找到这个超平面上足够多个“典型”的点，那么事实上，我们就获得了平面的信息。在更高维的空间中，一般的ODEs的解的分布区域也不是充满整个空间的，它是位于一个特定的平面或者流形上的。这揭示一种共性，即虽然是解空间是一个三维的空间，但是他们的解都位于一个二维的超平面里面，也就是说，用两个变量，就可以把三维的向量值解表达出来了。那么问题来了，上述的两个例子很简单，我们用肉眼都能观察到方程之间的关系，从而进行降阶，对于更一般的问题，我们应该怎么处理呢？...原创 2022-07-23 10:39:27 · 6064 阅读 · 0 评论 -
曲面曲率驱动的斑图传播
前几天 PRL 杂志发布了一篇文章，是关于图灵斑图的。大概意思是说，对于在曲面的上斑图演化，受到曲面的形状的影响。这个 Letter 首次报告了在平面上的静态图灵斑图在曲面上是会传播的，而我们以往的研究里，都是假设图灵斑图在曲面上也是静态的。这篇文章研究了反对称曲面上的反应扩散系统。结果表明，曲面的对称性和斑图的样子都会影响斑图的传播。作者研究了曲面的曲率对斑图传播影响的一个机制。一些研究已经表明了，当图灵斑图放到曲面上时，不稳定性条件是如何变化的，并且图案的位置如何受到曲面曲率的不均匀性的影响，这被称为原创 2022-06-12 13:12:38 · 546 阅读 · 0 评论 -
曲面抛物方程的质量集中法
曲面抛物方程的质量集中法该篇内容来自文章 Xiao X, Feng X, Yuan J. The lumped mass finite element method for surface parabolic problems: error estimates and maximum principle[J]. Computers & Mathematics with Applications, 2018, 76(3): 488-507.部分内容来自张林波老师 slides。特此说明。文原创 2022-05-29 10:08:04 · 1050 阅读 · 0 评论 -
特征值问题误差分析：Babuska–Osborn紧算子谱逼近理论
特征值问题误差分析：Babuska–Osborn紧算子谱逼近理论文章目录特征值问题误差分析：Babuska–Osborn紧算子谱逼近理论从一个简单的例子开始一维例子P1 元逼近一些观测和思考一些基本定义变分特征值问题特征值问题特征值刻画误差的刻画Babuska-Osborn 理论收敛性通常的框架收敛阶对称变分特征值问题上午好，各位。非常荣幸站在这里。且，今天我想给大家分享一些关于特征值问题的误差分析。关于我们的讨论班，我提议用说英文。就是，关于我们的讨论班，我提一个建议。我觉得大家可以英文来讲。事实上原创 2022-05-18 22:57:37 · 562 阅读 · 0 评论 -
紧算子的谱逼近理论及其在拉普拉斯特征值问题上的应用
紧算子的谱逼近理论及其在拉普拉斯特征值问题上的应用文章目录紧算子的谱逼近理论及其在拉普拉斯特征值问题上的应用紧算子的谱逼近理论预备知识收敛阶对称变分特征值问题拉普拉斯特征值问题简述T=TVT=T_VT=TV 时的分析T=THT=T_HT=TH 时的分析紧算子的谱逼近理论预备知识V1V_1V1 和 V2V_2V2 是复希尔伯特空间，我们需要找 λ∈C\lambda \in \mathbb{C}λ∈C 和 u∈V1u \in V_{1}u∈V1 且 u≠0u\neq 0u=0，使得，a(原创 2022-05-18 22:53:58 · 828 阅读 · 0 评论 -
密度泛函理论平面波基组展开
密度泛函理论：平面波基组展开文章目录密度泛函理论：平面波基组展开基本理论与概念原胞绝热近似、 HF 方法DFT 与 KS 方程自洽迭代交换关联势赝势DFT 应用赝势平面波 DFT 方法各项能量表达Fourier 展开动量空间总能量及发散问题动量空间 KS 方程截断能平面波 DFT 方法概述参考书基本理论与概念密度泛函理论 (Density Functional Theory, DFT) 是一种研究多电子体系电子结构的量子力学方法。密度泛函理论在物理、 化学、材料上都有广泛的应用, 特别是用来研究分子和原创 2022-05-18 22:44:05 · 1906 阅读 · 0 评论 -
tmux 的简单使用记录
tmux 的简单使用记录安装apt install tmux新建一个 tmux 会话tmux new -s abacus退出 tmux 会话窗口并杀死进程exitdetach 退出 tmux 但保持后台Ctrl+B + d查看有几个 tmux 会话窗口tmux ls杀死会话tmux kill-session -t abacusattach 进入后台运行的会话窗口tmux a -t abacus增加 windows 窗口Ctrl+B + c创建一个新窗口，状原创 2022-05-18 22:22:13 · 191 阅读 · 0 评论 -
小议移动有限元方法（转载）
小议移动网格有限元方法（转载）转载一篇提到移动网格有限元的文章。网格是偏微分方程数值解法的基础，网格体系的好坏直接影响计算结果的精度，甚至影响计算的成败。网格方法的研究经历了从结构化到非结构化，从单一网格到混合网格的过程。经过几十年的发展，这些网格方法已经很好地用于各种问题的计算，并不断出现新的针对不同情况的网格生成技术，而且形成了一些好的网格生成软件。近三十年来，自适应网格方法（主要有移动网格方法和局部细化或粗化的网格方法）一直引起了国际学术界和各类应用部门的高度重视，并且成为网格方法研究的热点问原创 2022-05-18 22:10:48 · 278 阅读 · 0 评论 -
电子结构计算简介
电子结构计算简介文章目录电子结构计算简介问题背景基本概念密度泛函理论方法波函数方法最基础的 Hartree-Fock 算法薛定谔方程Hartree-Fock 方法软件包求解薛定谔方程是定义在 R3N\mathbb{R}^{3 N}R3N 高维线性特征值问题. 除极少数体系外, 无法直接求解. 因此, 人们需要寻找与其等价或简化的可计算模型. 这些等价或简化的可计算模型中最典型的是 Hartree-Fock 方程和 Kohn-Sham 方程. 求解这两种方程，其对应了波函数方法和密度泛函理论方法。前原创 2022-05-04 23:08:28 · 2619 阅读 · 0 评论 -
水平集和符号距离函数
水平集和符号距离函数零水平集定义： 对于一个函数 ϕ(x⃗)：Rn→R\phi(\vec x)：{\mathbf{R}^n}\rightarrow \mathbf{R}ϕ(x)：Rn→R（其中x⃗∈Rn\vec x \in {\mathbf{R}^n}x∈Rn ，下同），取其值域为零部分对应的定义域：Γ={x⃗∣ϕ(x⃗)=0}(1)\Gamma = \{ \vec x|\phi (\vec x) = 0\} \tag{1}Γ={x∣ϕ(x)=0}(1)这里， Γ∈Rn−1\Gamma \in原创 2021-12-19 15:27:33 · 3229 阅读 · 1 评论 -
C++/C 使用过程中的经验总结荟萃琐碎汇总（持续更新）
同一个算法，不同的人写出来，效率上可能会差1000倍。写程序的同时，要注意性能。文章目录玄学命名空间左值非const结构体 set 和 mapSTL 容器选择关于-O3的过度优化内存上的性能优化关于 set 和 mapconst_cast的妙用关于函数传递数组与指针匿名表达式Debug 模式和 Release 模式全局变量Ubuntu 修改C头文件（链接库、环境变量）目录类中vector的初始化vector在使用迭代器遍历过程中，不能使用erase (),否则出错关于类型重命令explicit关键字s.原创 2021-12-12 13:34:26 · 2346 阅读 · 1 评论 -
图灵与图灵斑图
图灵与图灵斑图文章目录图灵与图灵斑图图灵简介图灵成就与贡献战争中的英雄计算机科学之父人工智能之父图灵斑图同性恋、迫害与平反怪人和运动天赋英镑和电影图灵斑图（Turing Pattern）自然和物理现象形成机理形而上的理解反应扩散系统形成条件图灵斑图的数值方法之参数化有限元方法简介演化表面上反应扩散方程的推导曲面梯度演化表面上的反应扩散系统模型说明曲面有限元方法变分形式演化曲面有限元方法时间离散肿瘤增长模型建立数值方法数值结果网格自适应和重新划分网格图灵斑图的 Trace 方法图灵简介图灵成就与贡献原创 2021-12-05 17:01:18 · 3268 阅读 · 0 评论 -
计算数学领域一些期刊（持续更新）
计算数学领域一些期刊（持续更新）前言前些日子某某同学告诉我，XX 期刊在计算数学领域排第四，我其实挺纳闷的，这么具体的一个数值是从哪来的呢。就挺好玩的认知。有的同学也总喜欢问，发的是不是 SCI，是几区的，有时候都不知道怎么回答。就好像有个声音突然有一天在我耳边问我，“你是不是个人？” 我总不能一本正经地回答，“我确实是个人” 吧，感觉怪怪的。可以说是，外行看分区，内行看内涵。只能说，这个问题本身没意义，很多二、三区的期刊也比一区的好，很多非 SCI 的新期刊，也比所谓的 SCI 要好。也有遇到生物原创 2021-11-28 14:29:13 · 10735 阅读 · 7 评论 -
有限元分析中的常识（持续更新）
有限元分析中的常识（持续更新）介绍一些学过有限元的人都不得不超级熟练掌握的基本常识。文章目录有限元分析中的常识（持续更新）通用符号基本不等式Cauchy-Schwarz 不等式Hölder's 不等式Young's 不等式Poincaré 不等式离散版本的 Gronwall 不等式Sobolev 不等式嵌入定理嵌入定理紧嵌入定理Cea 引理证明技巧Aubin-Nitsche 对偶技巧Lax-Milgram 定理通用符号∥⋅∥m,p\|\cdot\|_{m, p}∥⋅∥m,p 表示 Sobol原创 2021-11-20 19:45:55 · 2625 阅读 · 0 评论 -
梯度流是个什么玩意儿
梯度流是个什么玩意儿举个最最最最简单的例子。考虑一个经典的极小化问题：minxf(x)\min _{\mathbf{x}} f(\mathbf{x})xminf(x)其中 x∈Rn,f:Rn→R\mathrm{x} \in \mathbb{R}^{n}, f: \mathbb{R}^{n} \rightarrow \mathbb{R}x∈Rn,f:Rn→R 是一个二次可微的函数。x∗\mathrm{x}^{*}x∗ 是上述问题的解的必要条件是∇f(x∗)=0\nabla f\left原创 2021-11-15 10:40:26 · 5628 阅读 · 3 评论 -
肿瘤增长数学模型
曲面偏微分方程的一个实际应用——肿瘤增长文章目录曲面偏微分方程的一个实际应用——肿瘤增长简介演化表面上反应扩散方程的推导曲面梯度演化表面上的反应扩散系统模型说明曲面有限元方法变分形式演化曲面有限元方法时间离散肿瘤增长模型建立数值方法数值结果网格自适应和重新划分网格简介基本想法是利用表面有限元数值方法求解演化曲面上的反应扩散方程。不断增长的生物表面上的模式形成，是通过求解曲面上的反应扩散方程得到的，可以参考图灵的经典文章 1。这种方程表现出空间均匀结构的扩散驱动不稳定性，导致空间不均匀的图案。这项工作的原创 2021-11-06 00:00:00 · 3314 阅读 · 3 评论 -
常见的使用水平集函数的零水平集表示的曲面汇总（3D）（MATLAB 画三维图）
常见的使用水平集函数的零水平集表示的曲面汇总（3D）文章目录常见的使用水平集函数的零水平集表示的曲面汇总（3D）Donut甜甜圈Atom原子Popcorn爆米花Tooth牙齿Heart爱心Apple苹果Orthocircle正交圆Swiss cheese block瑞士奶酪块Heart心形Flower花代码感谢周同学和张师兄的公式提供。Donut甜甜圈水平集函数：ϕ(x,y,z)=(x2+y2−0.6)2+z2−0.09\phi(x, y, z)=\left(\sqrt{x^{2}+y^{2原创 2021-07-25 21:27:02 · 2667 阅读 · 8 评论 -
【小波分析】二、小波分析基础知识
【小波分析】二、小波分析基础知识文章目录【小波分析】二、小波分析基础知识基本记号信号空间内积和模（范数）傅里叶变换及其性质正变换（分解）逆变换（重构）傅里叶变换的性质小波基傅里叶变换优缺点傅里叶变换缺点的问题和本质空间问题数值计算偏差问题不可保性质的问题问题的本质和根源小波基小波定义小波的波形特征小波变换小波变换定义小波变换的基本性质Parseval恒等式（也称 能量守恒、内积恒等式 、Plancherel 定理）保长度逆变换其他性质小波变换的一些分析吸收小波吸收小波引入吸收小波定义吸收小波性质二进小波原创 2021-05-05 22:07:27 · 2177 阅读 · 2 评论 -
【小波分析】一、小波分析入门基础介绍
【小波分析】一、小波分析入门基础介绍文章目录【小波分析】一、小波分析入门基础介绍引入从线性代数谈起到高等数学历史和发展傅里叶分析哈尔基加窗傅里叶变换小波分析参考和致谢引入从线性代数谈起我们在线性代数当中学过点、向量、矩阵，线性变换等一些概念，我们一起回顾一下，从而引出函数空间，为小波分析的介绍做准备。我们知道在一个线性空间下，比如说三维的线性空间，有一个坐标系，我们喜欢叫 xxx 轴、yyy 轴、zzz 轴，在这个三维空间中的每一个点，都有一个坐标表达（点的表示），比如说 (x,y,z)(x,y,原创 2021-05-04 11:46:53 · 15941 阅读 · 8 评论 -
[有限元方法基础理论] 质量集中有限元（谱元）
[有限元方法基础理论] 质量集中有限元（谱元）集中质量有限元，也称谱元，要求具有对角质量矩阵，且对角线上的值等于行和。可以通过用数值积分公式以及其他方法来构造谱元。不妨考虑一个初边值问题，ut−Δu=f in Ω,t>0u=0 on ∂Ω,t>0, with u(⋅,0)=v in Ω\begin{aligned}&u_{t}-\Delta u=f \quad \text { in } \Ome原创 2021-01-12 22:49:50 · 2270 阅读 · 0 评论 -
曲面偏微分方程：参数化有限元方法
曲面偏微分方程：参数化有限元方法前面介绍的P\mathbf{P}P和Pd\mathbf{Pd}Pd，及其对应的同参等值延伸及投影自然延伸，对于曲面上有限元方法的理论分析起着至关重要的作用。利普希茨参数曲面上的FEM利普希茨参数化曲面因为同参映射P\mathbf{P}P的双利普希茨性，存在一个常数LLL，使得 (64) L−1∣x1−x2∣≤∣x~1−x~2∣≤L∣x1...原创 2019-11-06 23:40:21 · 1656 阅读 · 0 评论 -
曲面扰动理论之C2曲面的H2扩张
曲面扰动理论之C2曲面的H2扩张我们在后续提到的窄带（narrow band）方法和迹方法（trace）的分析中，对解要求：u~∈H2(γ)\widetilde{u} \in H^{2}(\gamma)u∈H2(γ)∥u∥H2(N(δ))≲δ12∥u~∥H2(γ)\|u\|_{H^{2}(\mathcal{N}(\delta))} \lesssim \delta^{\frac{1}{2...原创 2019-10-17 00:27:04 · 392 阅读 · 0 评论 -
曲面逼近和等参trace有限元
曲面逼近和等参trace有限元用水平集方法来表示曲面在处理演化曲面问题上用得比较多，即Γ\GammaΓ是水平集函数的一个水平集：Γ={x∈Ω:ϕ(x)=0}\Gamma=\{x \in \Omega : \phi(x)=0\}Γ={x∈Ω:ϕ(x)=0}如果用ϕh\phi_hϕh来表示ϕ\phiϕ的插值，即ϕh=Ik(ϕ)\phi_{h}=I^{k}(\phi)ϕh=Ik(ϕ)，我们...原创 2019-04-28 20:40:38 · 937 阅读 · 4 评论 -
曲面上变分的基础知识（方程组）
曲面上变分的基础知识（方程组）定理：曲面上变分的格林公式（向量形式）定理描述如下，为了方便，我们省略了严格的数学书写：∫ΓdivΓ(GP)⋅gds=−∫Γtr(G∇Γg)ds\int_{\Gamma} \operatorname{div}_{\Gamma}(\mathbf{GP} ) \cdot \mathbf{g} d s=-\int_{\Gamma} \operatorname{t...原创 2019-06-02 20:51:45 · 1257 阅读 · 0 评论 -
曲面上的微积分（一）
曲面上的微积分（一）参数化方法、迹方法和窄带方法是三种主要方法去求解Laplace-Beltrami问题。首先，我会介绍一些关于曲面上微积分的知识。参数化曲面定义曲面的参数化表示，χi:Vi→Ui∩γ⊂Rn+1\chi_{i} : \mathcal{V}_{i} \rightarrow \mathcal{U}_{i} \cap \gamma \subset \mathbb{R}^{n+1}...原创 2019-09-11 22:24:30 · 2425 阅读 · 2 评论 -
曲面上的微积分（二）
曲面上的微积分（二）从上一次讨论结束，回去之后，我花了一点时间走了一遍证明。心中依然遗留了两个问题。对于Laplace-Beltrami方程的弱形式，∫γφ~Δγv~=−∫γ∇γφ~⋅∇γv~\int_{\gamma} \tilde{\varphi} \Delta_{\gamma} \tilde{v}=-\int_{\gamma} \nabla_{\gamma} \tilde{\varph...原创 2019-09-18 21:46:07 · 552 阅读 · 0 评论 -
曲面扰动理论之C1曲面的扰动
曲面扰动理论之C1曲面的扰动定义曲面χΓ\chi_\GammaχΓ，定义uΓ∈H#1(Γ)u_{\Gamma} \in H_{\#}^{1}(\Gamma)uΓ∈H#1(Γ)，它是如下弱形式的解：(33)∫Γ∇ΓuΓ⋅∇Γv=∫ΓfΓv∀v∈H#1(Γ)(33) \quad \int_{\Gamma} \nabla_{\Gamma} u_{\Gamma} \cdot \nabla_{\...原创 2019-09-26 20:38:33 · 465 阅读 · 0 评论 -
Surfice PDE常用的微分几何工具
曲面偏微分常用工具稳定界面曲面上梯度的定义曲面上的梯度，无非就是梯度在曲面切面上的一个投影，我们关心得是投影算子P\mathbf{P}P，容易想到，向量v\mathbf{v}v的投影为：Pv=v−(n⋅v)n=v−nnTv=(I−nnT)v\mathbf{Pv}=\mathbf{v}-(\mathbf{n} \cdot \mathbf{v})\mathbf{n}=\mathbf{v}-\m...原创 2019-04-09 22:18:04 · 743 阅读 · 0 评论 -
曲面扰动理论之C2曲面的扰动
曲面扰动理论之C2曲面的扰动依然假设邻域：(48)N={x∈Rn+1:∣d(x)∣<12K∞}(48) \quad \mathcal{N}=\left\{\mathrm{x} \in \mathbb{R}^{n+1}:|d(\mathrm{x})|<\frac{1}{2 K_{\infty}}\right\}(48)N={x∈Rn+1:∣d(x)∣<2K∞1}假设有...原创 2019-10-10 21:19:40 · 280 阅读 · 0 评论