📄 二维、结构网格单相流simple 算法的c语言实现.htm
字号:
lang=EN-US>.lib文件。</SPAN></SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt">”</SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">%matlab%\extern\include</SPAN><SPAN
lang=EN-US style="FONT-SIZE: 9pt">”</SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">有产生<SPAN
lang=EN-US>.lib文件的.def文件。使用</SPAN></SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt">”</SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">%vc%\bin\lib.exe /def[:deffilename]
/out[:libfilename]</SPAN><SPAN lang=EN-US style="FONT-SIZE: 9pt">”</SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">将<SPAN
lang=EN-US>.def编译为.lib文件。%vc%是你的Visual C++的安装主目录。需要产生的.lib文件包括libeng.lib,
libmat.lib与libmx.lib。它们所使用的.def文件分别是libeng.def, libmat.def与libmx.def。不同版本的Visual
C++使用的.lib文件好像并不兼容。</SPAN></SPAN></P>
<P class=MsoNormal style="MARGIN-LEFT: 18pt; TEXT-INDENT: -18pt"><SPAN
lang=EN-US style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">4)<SPAN
style="FONT-WEIGHT: normal; FONT-SIZE: 7pt; FONT-STYLE: normal; FONT-FAMILY: Times New Roman; FONT-VARIANT: normal">
</SPAN></SPAN><SPAN style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">将<SPAN
lang=EN-US>libeng.lib,
libmat.lib与libmx.lib加入到你的项目中。VC7与VC6的设置见下面的图：</SPAN></SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US><!--[if gte vml 1]><v:shapetype
id=_x0000_t75 coordsize = "21600,21600" o:preferrelative = "t" o:spt = "75"
filled = "f" stroked = "f" path = " m@4@5 l@4@11@9@11@9@5 xe"><v:stroke
joinstyle = "miter"></v:stroke><v:formulas><v:f eqn =
"if lineDrawn pixelLineWidth 0 "></v:f><v:f eqn = "sum @0 1 0 "></v:f><v:f eqn =
"sum 0 0 @1 "></v:f><v:f eqn = "prod @2 1 2 "></v:f><v:f eqn =
"prod @3 21600 pixelWidth "></v:f><v:f eqn =
"prod @3 21600 pixelHeight "></v:f><v:f eqn = "sum @0 0 1 "></v:f><v:f eqn =
"prod @6 1 2 "></v:f><v:f eqn = "prod @7 21600 pixelWidth "></v:f><v:f eqn =
"sum @8 21600 0 "></v:f><v:f eqn = "prod @7 21600 pixelHeight "></v:f><v:f eqn =
"sum @10 21600 0 "></v:f></v:formulas><v:path o:extrusionok = "f"
gradientshapeok = "t" o:connecttype = "rect"></v:path><o:lock aspectratio="t"
v:ext="edit"></o:lock></v:shapetype><v:shape id=_x0000_s1025
style="WIDTH: 415.5pt; HEIGHT: 290.25pt" type = "#_x0000_t75" coordsize =
"21600,21600"><v:imagedata o:title="" src =
"index.files/image001.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=554 height=387
src="index.files/image002.jpg" v:shapes="_x0000_s1025"><![endif]></SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">
图23. 在VC7中添加matlab engine库文件</SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US><!--[if gte vml 1]><v:shape id=_x0000_s1026
style="WIDTH: 414.75pt; HEIGHT: 270pt" type = "#_x0000_t75" coordsize =
"21600,21600"><v:imagedata o:title="" src =
"index.files/image003.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=553 height=360
src="index.files/image004.jpg" v:shapes="_x0000_s1026"><![endif]></SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">
图24. 在VC6中添加matlab engine库文件</SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: center" align=center><SPAN
lang=EN-US></SPAN> </P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">另外，关于制作matlab 5.3下动画制作的问题</SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">可以</SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">用下面的方法</SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">mpgwrite :matlab user contributed
function<BR>功能，把一系列matlab的图像帧转化为mpg文件<BR>平台matlab5
:win32,solaris,sgi,hpux,linux<BR></SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">具体可以参考</SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体"> </SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体"><A
href="http://fluid.stanford.edu/~fringer/movies/gif_howto/gif_howto.html">http://fluid.stanford.edu/~fringer/movies/gif_howto/gif_howto.html</A><BR><A
href="http://www.rpi.edu/~holmes/NumDiffEqs/Demos/matlab.movie.info.pdf">http://www.rpi.edu/~holmes/NumDiffEqs/Demos/matlab.movie.info.pdf</A></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体"><BR>需要从www.mathworks.com下载一个文件</SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">,我已经上载到 <A
href="http://www.flotu.org/~qzhen/cfd/simplec/mpgwrite.tar">http://</A></SPAN><A
href="http://www.flotu.org/~qzhen/cfd/simplec/mpgwrite.tar"><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">www.flotu.org/~qzhen/</SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">cfd/simplec/</SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">mpgwrite</SPAN><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">.</SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">tar</SPAN></A><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体"> </SPAN><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体"> 2M</SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: center" align=center><B><SPAN lang=EN-US
style="FONT-SIZE: 12pt; FONT-FAMILY: 宋体"></SPAN></B> </P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><B><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">4.1 例子exam1：二维稳态传热</SPAN></B></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">边壁上的温度<SPAN lang=EN-US><SPAN
style="POSITION: relative; TOP: 5pt"><!--[if gte vml 1]><v:shape id=_x0000_s1028
style="WIDTH: 60.75pt; HEIGHT: 15.75pt" type = "#_x0000_t75" coordsize =
"21600,21600"> <v:imagedata o:title="" src =
"index.files/image005.wmz"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=81 height=21
src="index.files/image006.gif" v:shapes="_x0000_s1028"><![endif]></SPAN><!--[if gte mso 9]><xml>
<o:OLEObject Type="Embed" ProgID="Equation.DSMT4" ShapeID="_x0000_i1043" DrawAspect="Content" ObjectID="_1100594551">
</o:OLEObject>
</xml><![endif]-->。
</SPAN></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">扩散系数<SPAN lang=EN-US><SPAN
style="POSITION: relative; TOP: 2pt"><!--[if gte vml 1]><v:shape id=_x0000_s1029
style="WIDTH: 27pt; HEIGHT: 12.75pt" type = "#_x0000_t75" coordsize =
"21600,21600"> <v:imagedata o:title="" src =
"index.files/image007.wmz"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=36 height=17
src="index.files/image008.gif" v:shapes="_x0000_s1029"><![endif]></SPAN><!--[if gte mso 9]><xml>
<o:OLEObject Type="Embed" ProgID="Equation.DSMT4" ShapeID="_x0000_i1044" DrawAspect="Content" ObjectID="_1100594552">
</o:OLEObject>
</xml><![endif]-->。扩散系数在user_gamsor( )中设定。</SPAN></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">求解域<SPAN lang=EN-US>x方向宽度＝1.0
。求解域y方向宽度＝2.0 。</SPAN></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left> </P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">计算结果：</SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN lang=EN-US><!--[if gte vml 1]><v:shape id=_x0000_s1030
style="WIDTH: 415.5pt; HEIGHT: 327pt" type = "#_x0000_t75" coordsize =
"21600,21600"><v:imagedata o:title="" src =
"index.files/image009.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=554 height=436
src="index.files/image010.jpg" v:shapes="_x0000_s1030"><![endif]></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: center" align=center><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">图<SPAN lang=EN-US>25.
exam1的计算结果－二维温度场的分布</SPAN></SPAN></P>
<P class=MsoNormal><B><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体"></SPAN></B> </P>
<P class=MsoNormal><B><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">4.2 例子exam10：湍流中被动标量扩散（基于时均的 k
–ε方程）</SPAN></B></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体"></SPAN> </P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">exam10中没有考虑温度对于密度以及粘度的影响，所以可以把温度当成只是一个被动标量。</SPAN></P>
<P class=MsoNormal><SPAN style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">入口左半边的速度为<SPAN
lang=EN-US>5m/s，温度为100摄氏度；右半边的速度为10m/s，温度为400摄氏度。</SPAN></SPAN></P>
<P class=MsoNormal><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">计算结果：</SPAN></P>
<P class=MsoNormal> </P>
<P class=MsoNormal><SPAN lang=EN-US><!--[if gte vml 1]><v:shape id=_x0000_s1053
style="WIDTH: 341.25pt; HEIGHT: 477pt" type = "#_x0000_t75" coordsize =
"21600,21600" o:borderleftcolor = "this" o:bordertopcolor = "this"
o:borderrightcolor = "this" o:borderbottomcolor = "this"><v:imagedata o:title=""
src = "index.files/image011.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=455 height=636
src="index.files/image012.jpg" v:shapes="_x0000_s1053"><![endif]><!--[if gte vml 1]><v:shape
id=_x0000_s1054 style="WIDTH: 70.5pt; HEIGHT: 370.5pt" type = "#_x0000_t75"
coordsize = "21600,21600" o:borderleftcolor = "this" o:bordertopcolor = "this"
o:borderrightcolor = "this" o:borderbottomcolor = "this"><v:imagedata o:title=""
src = "index.files/image013.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=94 height=494
src="index.files/image014.jpg" v:shapes="_x0000_s1054"><![endif]></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: center" align=center><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体"></SPAN> </P>
<P class=MsoNormal style="TEXT-ALIGN: center" align=center><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">图<SPAN lang=EN-US>26.
exam10的计算结果－温度与速度的分布</SPAN></SPAN></P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN lang=EN-US
style="BORDER-RIGHT: windowtext 1pt solid; PADDING-RIGHT: 0cm; BORDER-TOP: windowtext 1pt solid; PADDING-LEFT: 0cm; FONT-SIZE: 9pt; PADDING-BOTTOM: 0cm; BORDER-LEFT: windowtext 1pt solid; PADDING-TOP: 0cm; BORDER-BOTTOM: windowtext 1pt solid; FONT-FAMILY: 新宋体"></SPAN> </P>
<P class=MsoNormal><B><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">4.3 例子exam6：常粘度系数湍流的动态模拟</SPAN></B></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体"></SPAN> </P>
<P class=MsoNormal><SPAN style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">一股热流体（速度<SPAN
lang=EN-US>100m/s，温度500摄氏度）从右下角注入系统（壁面温度为300摄氏度），出口的流速分布设为均匀。</SPAN></SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体"></SPAN> </P>
<P class=MsoNormal><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">微分方程包括时均的连续性方程、动量、能量的输运方程。所用的粘度系数是<SPAN
lang=EN-US>1。热传导系数的量级与粘度系数相仿。</SPAN></SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体"></SPAN> </P>
<P class=MsoNormal style="TEXT-ALIGN: left" align=left><SPAN lang=EN-US><!--[if gte vml 1]><v:shape id=_x0000_s1055
style="WIDTH: 141.75pt; HEIGHT: 553.5pt" type = "#_x0000_t75" coordsize =
"21600,21600" o:borderleftcolor = "this" o:bordertopcolor = "this"
o:borderrightcolor = "this" o:borderbottomcolor = "this"><v:imagedata o:title=""
src = "index.files/image015.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=189 height=738
src="index.files/image016.jpg" v:shapes="_x0000_s1055"><![endif]><!--[if gte vml 1]><v:shape
id=_x0000_s1056 style="WIDTH: 99pt; HEIGHT: 528.75pt" type = "#_x0000_t75"
coordsize = "21600,21600" o:borderleftcolor = "this" o:bordertopcolor = "this"
o:borderrightcolor = "this" o:borderbottomcolor = "this"><v:imagedata o:title=""
src = "index.files/image017.png"></v:imagedata></v:shape><![endif]--><![if !vml]><img border=0 width=132 height=705
src="index.files/image018.jpg" v:shapes="_x0000_s1056"><![endif]></SPAN></P>
<P class=MsoNormal><SPAN lang=EN-US
style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体"></SPAN> </P>
<P class=MsoNormal style="TEXT-INDENT: 22.5pt"><SPAN
style="FONT-SIZE: 9pt; FONT-FAMILY: 宋体">图<SPAN lang=EN-US>27.
exam6的计算结果－温度与速度的分布</SPAN></SPAN></P>
<P><SPAN style="FONT-SIZE: 9pt; FONT-FAMILY: 新宋体">动画文件可以从<SPAN lang=EN-US><A
style="COLOR: blue; TEXT-DECORATION: underline; text-underline: single"
href="http://www.flotu.org/~qzhen/cfd/simplec/exam6.mpg">http://www.flotu.org/~qzhen/cfd/simplec/exam6.mpg</A>
下载。</SPAN></SPAN></P>
<P> </P></BODY></HTML>
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -