| 查看: 5227 | 回复: 116 | |||
[交流]
激光反射吸收热源模型！！！！代码重量爆料！！！！！求解答！！！！！！！！！
|
|||
|
代码有点错误，我现在是想用最小二乘法拟合自由表面的曲线，首先轮循所有自由表面网格，计算自由表面网格数量N，开辟数值存储自由表面上网格的重心（优良拟合三次曲线方程），gt,GS,max三个函数是用来求三次曲线方程y=a+bx+cx2+dx*x*x四个系数。然后把求出的四个系数保存到C_UDMI中，用于与离散的激光光线求交。代码编译之后有错误，可能思路上都已经与udf编程思想向偏离了，求各路大神指点，重谢！！！！！！！！！！！至少400金币！！！！呵呵 #include "udf.h" #include "sg.h" #include "sg_mphase.h" #include "flow.h" #include <math.h> double gt(double x[],double y[]); double GS(double a[4][4],double b[4],double x[4]); double max(double array[4]); DEFINE_ADJUST(store_gradient,domain) { int phase_domain_index = 0; real xc[ND_ND],xd[ND_ND],x[4],n1,n2,thita_l,thita_g; Thread *t; Thread **pt; cell_t c; Domain *pDomain = DOMAIN_SUB_DOMAIN(domain,phase_domain_index); { Alloc_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL); Scalar_Reconstruction(pDomain, SV_VOF,-1,SV_VOF_RG,NULL); Scalar_Derivatives(pDomain,SV_VOF,-1,SV_VOF_G,SV_VOF_RG,Vof_Deriv_Accumulate); } mp_thread_loop_c(t,domain,pt) { if (FLUID_THREAD_P(t)) { int N=0; Thread *ppt = pt[phase_domain_index]; begin_c_loop (c,t) { C_CENTROID(xd,c,t); if(C_VOF(c, ppt)>0.1&&C_VOF(c, ppt)<0.9&&fabs(xd[0]-0)>0.0005) N++; } end_c_loop (c,t) double a[N],b[N],c[N],d[N],e[N],f[N],g[N]; int i; for(i=0; i<N; i++) { a=1; } int j=0; begin_c_loop (c,t) { C_CENTROID(xc,c,t); if(C_VOF(c, ppt)>0.1&&C_VOF(c, ppt)<0.9&&fabs(xc[0]-0)>0.0005) { b[j]=xc[0]; e[j]=xc[1]; c[j]=b[j]*b[j]; d[j]=b[j]*b[j]*b[j]; j++; } } end_c_loop (c,t) f[0]=gt(a,a);f[1]=gt(a,b);f[2]=gt(a,c);f[3]=gt(a,d);f[4]=gt(b,b);f[5]=gt(b,c);f[6]=gt(b,d); f[7]=gt(c,c);f[8]=gt(c,d);f[9]=gt(d,d);g[0]=gt(a,e);g[1]=gt(b,e);g[2]=gt(c,e);g[3]=gt(d,e); double array[4][4]={{f[0],f[1],f[2],f[3]},{f[1],f[4],f[5],f[6]},{f[2],f[5],f[7],f[8]},{f[3],f[6],f[8],f[9]}}; GS(array,g,x); C_UDMI(c,t,0) = x[0]; C_UDMI(c,t,1) = x[1]; C_UDMI(c,t,2) = x[2]; C_UDMI(c,t,3) = x[3]; } } Free_Storage_Vars(pDomain,SV_VOF_RG,SV_VOF_G,SV_NULL); } double gt(double x[],double y[]) { double sum=0; int i; for(i=0;i<N;i++) sum=x*y+sum; return sum; } double GS(double a[4][4],double b[4],double x[4]) { double c[4]={0}; double x0[4]={0}; int i,j,k; double sum=0; for(k=1;;k++) { for(i=0;i<4;i++) { for(j=0;j<4;j++) { sum=a[j]*x0[j]+sum; } x=x0+(b-sum)/a; c=fabs(x-x0); x0=x; sum=0; } r=max(c); if(r<0.0001) { break; } } return 0; } double max(double array[4]) { double a=array[0]; int i; for(i=0;i<4;i++) { if(a<array) a=array; } return a; } |
回复此楼» 猜你喜欢
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有159人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有27人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有44人回复
» 抢金币啦！回帖就可以得到:
-
南京大学智能驱动与感知材料实验室诚招博士/硕士/科研助理/访问学生
+1/279
-
西湖大学工学院智能微型机器人实验室博士后招聘启事
+4/228
-
散金币，求好运，祝面上顺利！
+2/112
-
山东师范大学有机化学专业胡忠燕老师课题组招收2026届硕士研究生以及调剂生
+1/40
-
欢迎交流咨询
+1/38
-
湘潭大学“过程强化与绿色化工”创新团队补招2026年秋入学博士生
+1/35
-
上海交通大学电催化方向博后招聘-有锂电、高分子、燃料电池背景者优先
+1/31
-
华南师范大学（211）申请审核制博士招生- 电子、自动化、机械、生物学、物理相关专业
+1/25
-
欢迎生物与医药、药学、化学等相关专业的同学
+1/20
-
重庆文理学院化学与环境工程学院 2026 届硕士研究生调剂通知
+1/17
-
北京某研究院生物化学与分子生物学相关专业申请考核制博士招生（春季）
+1/14
-
南昌航空大学代威力教授招收光催化和光热催化研究方向博士生
+1/10
-
香港大学David Srolovitz教授课题组招聘「香江学者」博士后
+1/7
-
福建师范大学化学与材料学院杜克钊团队博士招生
+1/7
-
首次招收资格，招收工科调剂！
+1/3
-
国家杰青低维材料与器件力学团队2026年招收博士研究生
+1/3
-
2026年上海大学王亮课题组招收电催化或光催化研究背景的博士生
+2/2
-
Top-88悉尼科技大学数据科学/AI 招收2027年入学 校奖 博士生1到2名(国际和本地学生）
+1/2
-
香港科技大学（广州） 招收温室气体监测与减排方向的博士、博后、硕士、RA
+1/1
-
静电纺丝：工艺原理、形貌调控与典型应用
+1/1
11楼2013-12-02 11:06:04
12楼2013-12-02 11:06:09
20楼2013-12-02 11:26:38
22楼2013-12-02 11:37:53
37楼2013-12-03 10:27:32
42楼2013-12-03 11:03:56
43楼2013-12-03 11:04:57
48楼2013-12-03 11:23:48
简单回复
湮没黑洞2楼
2013-12-02 10:21
回复
焊接网格划分(金币+3): 谢谢参与
051031113楼
2013-12-02 10:25
回复
焊接网格划分(金币+3): 谢谢参与
风雨飘然4楼
2013-12-02 10:29
回复
焊接网格划分(金币+3): 谢谢参与
2013-12-02 10:30
回复
焊接网格划分(金币+3): 谢谢参与
haixiawu6楼
2013-12-02 10:30
回复
焊接网格划分(金币+3): 谢谢参与
tzynew7楼
2013-12-02 10:33
回复
焊接网格划分(金币+3): 谢谢参与
slz20088楼
2013-12-02 10:34
回复
焊接网格划分(金币+3): 谢谢参与
2013-12-02 10:48
回复
焊接网格划分(金币+3): 谢谢参与
wlxiao10楼
2013-12-02 10:57
回复
焊接网格划分(金币+3): 谢谢参与
峰锋风逢烽13楼
2013-12-02 11:08
回复
焊接网格划分(金币+3): 谢谢参与
doctorbang14楼
2013-12-02 11:14
回复
焊接网格划分(金币+3): 谢谢参与
wullww15楼
2013-12-02 11:18
回复
焊接网格划分(金币+3): 谢谢参与
aboluo19833116楼
2013-12-02 11:23
回复
焊接网格划分(金币+3): 谢谢参与
pzchen17楼
2013-12-02 11:24
回复
焊接网格划分(金币+3): 谢谢参与
xianzhi201018楼
2013-12-02 11:24
回复
焊接网格划分(金币+3): 谢谢参与
fujiating19楼
2013-12-02 11:24
回复
焊接网格划分(金币+3): 谢谢参与
wendytxy21楼
2013-12-02 11:33
回复
焊接网格划分(金币+3): 谢谢参与
abczhang23楼
2013-12-02 11:41
回复
焊接网格划分(金币+3): 谢谢参与
Spinage24楼
2013-12-02 11:41
回复
焊接网格划分(金币+3): 谢谢参与
顶
chlq314101a25楼
2013-12-02 11:43
回复
焊接网格划分(金币+3): 谢谢参与
宋全伟26楼
2013-12-02 11:44
回复
焊接网格划分(金币+3): 谢谢参与
， [ 发自手机版 http://muchong.com/3g ]
会飞的乌龟--27楼
2013-12-02 11:44
回复
焊接网格划分(金币+3): 谢谢参与
wander051528楼
2013-12-02 11:45
回复
焊接网格划分(金币+3): 谢谢参与
yuzhichao29楼
2013-12-02 11:50
回复
焊接网格划分(金币+3): 谢谢参与
祝福
Sammy230楼
2013-12-02 11:52
回复
焊接网格划分(金币+3): 谢谢参与
doublehaohao31楼
2013-12-02 11:58
回复
焊接网格划分(金币+3): 谢谢参与
尘轩无梦32楼
2013-12-02 11:59
回复
焊接网格划分(金币+3): 谢谢参与
、 [ 发自手机版 http://muchong.com/3g ]
lgward33楼
2013-12-02 12:01
回复
焊接网格划分(金币+3): 谢谢参与
huiminchlq34楼
2013-12-02 12:01
回复
焊接网格划分(金币+3): 谢谢参与
xiejf35楼
2013-12-02 12:09
回复
焊接网格划分(金币+1): 谢谢参与
silence136楼
2013-12-02 12:25
回复
ts_chenchao38楼
2013-12-03 11:00
回复
焊接网格划分(金币+1): 谢谢参与
天涯201239楼
2013-12-03 11:02
回复
焊接网格划分(金币+1): 谢谢参与
neu23440楼
2013-12-03 11:02
回复
焊接网格划分(金币+1): 谢谢参与
miaojiabing41楼
2013-12-03 11:03
回复
焊接网格划分(金币+1): 谢谢参与
gyfwon44楼
2013-12-03 11:09
回复
焊接网格划分(金币+1): 谢谢参与
gander45楼
2013-12-03 11:09
回复
焊接网格划分(金币+1): 谢谢参与
v
liu12333846楼
2013-12-03 11:11
回复
焊接网格划分(金币+1): 谢谢参与
lefantian12247楼
2013-12-03 11:18
回复
焊接网格划分(金币+1): 谢谢参与
zhoupeng8749楼
2013-12-03 11:32
回复
焊接网格划分(金币+1): 谢谢参与
·
base0803030850楼
2013-12-03 11:34
回复
焊接网格划分(金币+1): 谢谢参与
