[转帖]问与答

hyacinth

问：如何计算物体在水波作用下运动的问题？大家用Morison公式结合波的方程对该问题求解，一般方法是用波的方程计算流体的速度和压力，再用Morison公式计算物体受到力。Morison公式为：f=Cd*密度*d*u*u+Cm*d*d*pi/4*du/dt,公式中的速度u和速度导数du/dt是怎么求解的?
答：如果要想求得流场内的速度分布，你首先要做的就是选择合适的波浪方程，然后结合里面放置结构物后的边界条件，利用数值方法计算出波浪方程中的速度分布，然后可以利用墨尔森公式计算受力情况，但是墨尔森公式适合于小杆件的计算，有一些假设条件需要注意！

问：求解无量纲化的方程和有量纲的方程有什么不同，那个比较好？
无量纲的方程求解可以把计算域模型化，也便于分析结果。但是很难说是否无量纲化的方程比较好,根据情况选择合适的方程.其实可以把结果计算出来以后再化成无量纲形式的结果。

问：有限差分和有限体积法在应用上差别
答：有限差分方法必须是均匀网格，而有限体积方法可以是非均匀网格，甚至是无结构网格。

问：求解大型稀疏矩阵一般都用什么方法，求解速度如何？
答：1.预处理双稳定共轭梯度算法,计算速度和稳定性都还可以!
2.GMRES。
3采用LU分解，应该可行！
4.综合来讲， Bio- Conjugated Gradient Method （BiCG） 是性能最好的，一般的 比较的迭代揭发，Matlab 都有 ， 网上也有原代码下载。