[求助]如何在UDF中得到压力梯度的散度？

zhouxytj

即求  ▽.▽p[DISABLELBCODE]

zhouxytj

我的微分方程中没有三阶偏导，只是在一般的传输方程中多了一个源项，即压力梯度的
散度C▽.▽p (C为常数) 。因为我是刚接触UDF，所以不知道还有更好的方法。
下面的方法是否可行：
由于对方程积分后，源项为∫(C▽.▽p) dV =∫( d(C▽p,i)/dx,i ) dV =∑( (C▽p,i).(dA,i) )
在实际添加UDF源项时的定义为 ∫(C▽.▽p) dV / △V， 即 ∑( (C▽p,i).(dA,i) ) / △V。
即遍历每个cell的所有face，得到其face上的▽p及面积矢量，把▽p在i方向的分量乘以面积在i方
向的分量，再把三个方向的乘积相加（包含该cell的所有面），再除以cell的体积。
可以这样做吗？
谢谢！[DISABLELBCODE]

kaiser

从理论上来讲似乎没有什么错误

mingjing

沙子和眼睛的故事
她和他相识了，邂逅在初夏暖暖的沙滩上。她被他粗犷豪放的性格渲染着，那宽宽的肩膀像一座山，让她无依的心，无条件的停靠在他的港湾里。
当他有力的手握紧她时，她没有拒绝。而是任他牵着自己的手，相伴而行。排排海浪戏谑的卷过，吞噬着歪歪扭扭浅浅的两行足迹。她回眸浅笑，因为这段美好会深深的印无数个甜蜜的日子里，她由封闭的情感中走了出来。开始接纳着海的温存，浪花的洗礼。开始在他无微的呵护中编织着一个个旖旎的梦。
有一次，他托着她灿烂的笑脸说“此生我愿成为你心海的一粒沙，经风雨的磨砺后蜕变成金”。她话语很少，没有向他表白什么。只是在心里告诉他，“此生愿为菩萨玉瓶里的一滴甘露，让他无烦无恼，无伤无痛，永永远远滋润在他那片寂寞的沙洲”。
走累的时候，他喜欢背起她。看到他洒落的汗水，她心痛的挣扎着想逃离他的臂弯。但是他把她箍的很紧。那时她就想，她要做他一辈子的影子，跟随他今生来世。他背着她，在心里告诉自己。她是他甜蜜的负担。每一个时针的转动，都是他为她跳动的幸福音符。
就这样，他们决定结婚了。结婚前的那天，他告诉她，让她去海边等他。他要给她一份惊喜。她懂他，心里甜蜜的幸福着。他是要让大海和千年巨石，作为他们爱的见证。可是，那天的他失约了。她从上午等到中午。电话打了无数次，无人接听。她发疯似的到处找他。却在朋友打来的电话里得知，他在医院了。在医院的病床前，她看到了他因车祸而被碾压的面目全非的冰凉的尸体，手里还紧紧握着一大束红色的玫瑰。
以后的日子里，郑州治白癜风哪家好，她不知什么是白天，什么是黑夜。父母轮流守着她。终于在父母极度困乏的时候。她逃离了他们的视线。独自走在沙滩上，麻木的她任海风吹乱秀发。
一粒沙子随风飘进了她的眼睛，她泪眼滂沱。“是你吗?，为什麽你在心海跑出来，成为我眼睛里的一粒尘呢?。让我今生为你痛，为你流泪呢?”她心里嘶喊着，没有去揉她的眼睛。她要留它在眼眸里，她这样想着，慢慢的向大海深处走去••••••。