|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
x
最近学习Lee&Liu的2010年文章‘Lattice Boltzmann simulations of micron-scale drop impact on dry surfaces’的HCZ的拓展模型,处理边界的时候,边界湿润性需要采用三次自由能处理。
现在遇到一点问题,主要在边界处的差分处理,请教下论坛有经验的各位。
文章提到,计算化学势的二阶导数时,为保证守恒,因此需要[math]n\cdot \bigtriangledown \mu =0[/math]
但是在计算参数C的二阶导数时,由于需要体现接触角,和化学势的计算有所不同,需要用到[math]n\cdot \bigtriangledown C=\frac{\phi }{\kappa }(C-C^{2})[/math]
问题就是,这个式子是只在计算二阶导数是使用还是计算一阶导数就使用这个。
如果是计算一阶导数需要这个,那么,中心差分可以理解,偏差分该怎么处理?
困扰好久了,试了好多都搞不定啊啊啊。。。
还请做过的大大们指点下啊~~~
|
|