Fluid Simulation for Computer Graphics第二版，8.2节。
在求解自由表面流体时，简单地把计算域划分为液体格和气体格，会导致artifact.如果能更精确地追踪表面，可以修改一下在水-气界面附近计算压强差的方式。
根据Gibou et al. 2002[……]
作者： wmdcstdio
Ubuntu 20.04 Server软路由配置折腾记录
记录一下最近的一次折腾，使用一台安装了Ubuntu 20.04 Server的工控机搭建软路由，并成功配置家庭NAS的内网访问。计算机网络没好好学，全靠本能debug，感谢lcy同志的全程技术科普（
Real-time exploration of regular volume data by adaptive reconstruction of isosurfaces读书笔记
本文讨论使用Octree（八叉树）重建表面的问题。
3 Continuous isosurfaces
讨论连续性问题。
如果在Octree上储存数据，那么在不同大小的格子交界的地方，数据可能会不连续。
本文采用的方法：在大小格子交界处，让大格子在小格子上sample数值：
[……]
关于马拉戈尼效应(Marangoni Effect)
本文基于维基百科内容写作：
https://en.wikipedia.org/wiki/Marangoni_effect
https://en.wikipedia.org/wiki/Marangoni_number
示例视频（需科学上网）：
OpenGL屏幕像素坐标问题
传给glutMouseFunc()的(x,y)是左上至右下的坐标：x表示离屏幕左边缘的像素数，y表示离屏幕上边缘的像素数，如图：
而glReadPixels()采用的坐标系统则与此不同，它是从左下到右上的。也就是说，x表示离屏幕左边缘的像素数，y表示离屏幕下边缘的像素数，如图：
[……]
读书笔记：A Model for Soap Film Dynamics with Evolving Thickness
3 SOAP FILM DYNAMICS WITH VARYING THICKNESS
3.1 Definitions and Setting
流体在薄膜上的流动遵循带有表面张力项的无粘NS方程：
这个x上面一点就是速度的意思。(1)右端第一项是表面张力。里面的u代表气-液界面[……]
Efficient elasticity for character skinning with contact and collisions部分翻译（碰撞检测部分）
7 约束和碰撞
我们使用点约束（point constraints），这样一方面施加软体约束（例如bone attachments），一方面解决物体之间的碰撞/自碰撞。特别地，我们在规则网格（lattices）上放一些代理点（proxy points），在这些点上记录力，并用三次插值法把这些力[……]
读书笔记：An Introduction to Physically Based Modeling: Rigid Body Simulation I—Unconstrained Rigid Body Dynamics(SIGGRAPH’97 COURSE NOTES)
1 模拟基础
用“状态向量”表示一个质点的状态：
如果三维，就是六个数。这个Y可以扩展到n个质点：
我们先暂且考虑一个质点。F(t)是该质点t时刻所受合力。假设该点有质量m，那么状态向量的变化率就是：
2 刚体基本概念
2.1 位置和方向
质点没有方向[……]
Peskin2002：The immersed boundary method阅读笔记
2 运动方程
文中推导的是不可压缩弹性固体的方程。目标是让弹性方程看起来尽量像流体方程一样。
(q,r,s)是固体上的曲线坐标。X(q,r,s,t)是该固体质点在t时刻的位置。M(q,r,s)是质量密度。由能量密度函数E[X]定义固体材料性质。弹性力就是它的Frechet导数：
[……]
The Material Point Method for Simulating Continuum Materials读书笔记
3 简介
MPM是结合PIC和FLIP而发展的一种算法。MPM对Lagrangian mesh连通性没有要求。
和PIC/FLIP类似，MPM算法在背景Eularian网格的辅助下，隐式处理自碰撞和破碎。与传统的Lagrangian方法（例如FEM算固体）Eularian方法（流体）相比[……]