CST MWS同軸連接器自適應(yīng)網(wǎng)格優(yōu)化發(fā)現(xiàn)的問(wèn)題
閑暇之余,按照《微波工作室用戶全書(shū)》卷二的同軸連接器的例子做仿真,結(jié)果發(fā)現(xiàn)了一個(gè)非常有意思的地方,采用自適應(yīng)網(wǎng)格加密(Energy based&Expert based)后,查看結(jié)果的Mesh Adaptation Logfile,結(jié)果發(fā)現(xiàn)Pass2計(jì)算速度居然比Pass1快,很是驚訝!因?yàn)樽赃m應(yīng)加密后,Pass2的網(wǎng)格數(shù)是Pass1的23倍,理論上應(yīng)該P(yáng)ass2的仿真時(shí)間更長(zhǎng)啊。
Energy based: Mesh Adaptation Logfile
Expert based: Mesh Adaptation Logfile
不知道大家有沒(méi)有研究過(guò)這方面的情況?
我猜的是,加密后的情況只計(jì)算了加密的部分,未加密的地方?jīng)]有進(jìn)行計(jì)算。
我覺(jué)得在第二次加密的過(guò)程中可能出現(xiàn)網(wǎng)絡(luò)的最小分辨率變大,網(wǎng)格總數(shù)比第一次的網(wǎng)格總數(shù)多。而且計(jì)算總的時(shí)間與網(wǎng)格數(shù)目以及網(wǎng)格的分辨率是有關(guān)系的。所以出現(xiàn)上述的情況應(yīng)該是合理的。
而且除此之外,第二次計(jì)算時(shí)可能不需要重復(fù)第一次的某些初始化計(jì)算,因此可以減少這部分的時(shí)間。
在優(yōu)化中,也會(huì)出現(xiàn)第二次優(yōu)化的時(shí)間比第一次短的情形,雖然網(wǎng)格總數(shù)目沒(méi)有變化。
你說(shuō)的最小分辨率變大的情況不成立,因?yàn)锳daptive mesh refinement是朝著加密方向優(yōu)化的,因而網(wǎng)格密度只會(huì)增加或者不變,絕不會(huì)減小,下圖為證
Energy based Global Mesh
Expert based Global Mesh
No refinement Global Mesh
初始化計(jì)算這個(gè)概念太籠統(tǒng)了,缺乏根據(jù)和說(shuō)服力
要是這樣就好了,如果網(wǎng)格數(shù)表較多的情況下進(jìn)行優(yōu)化,也不會(huì)占用太多時(shí)間;不過(guò)CST有這么智能嗎?可以分著算嗎?希望有算法方面的高手給詳解
以下是個(gè)人理解,僅供參考。
首先,1樓的情況可能只能通過(guò)CST官方解答了。只看幫助文件《The Simulation Method》好像并沒(méi)有明確指出電磁場(chǎng)參數(shù)只在單個(gè)網(wǎng)格內(nèi)計(jì)算還是會(huì)受到相鄰網(wǎng)格計(jì)算數(shù)值的影響,不了解具體的算法原理,不能確定。因此不知道如果網(wǎng)格大小保持不變的話,Solver是否會(huì)保存上一次的計(jì)算結(jié)果。
第二,每次adaptive mesh,都要重新進(jìn)行初始化計(jì)算,所以2樓所說(shuō)的后半部分是不正確的。當(dāng)然,2樓的前半部分也持保留意見(jiàn),這個(gè)應(yīng)該不可能出現(xiàn)。
第三,樓主的這個(gè)模型solver的計(jì)算時(shí)間很短,很難判斷是哪個(gè)“環(huán)節(jié)”的時(shí)間被縮短了,excitation duration、matrix calculation time、solver setup time、solver loop time、solver postprocessing time,這些參數(shù)哪些變化了哪些沒(méi)有變樓主給出的數(shù)據(jù)體現(xiàn)不出來(lái)。
第四,并不能說(shuō)網(wǎng)格數(shù)增加仿真時(shí)間就必定延長(zhǎng)。有沒(méi)有可能在“大網(wǎng)格”條件下模型內(nèi)的能量因?yàn)?ldquo;諧振”而耗散得慢,在“小網(wǎng)格”條件下計(jì)算沒(méi)有出現(xiàn)“諧振”因而能量耗散得快?我想沒(méi)有人可以很確定得回答是或者不是吧。
建議:重新運(yùn)行,在transient solver parameter對(duì)話框里勾選"Store result data in cache",看看能不能把每一次的solver logfile保存下來(lái)。如果不能的話,有可能需要通過(guò)VBA把solver logfile保存起來(lái),先要確定是哪個(gè)環(huán)節(jié)的時(shí)間被縮短了。
收兵回營(yíng)!
有一點(diǎn)我沒(méi)搞明白,為什么大網(wǎng)格下諧振小網(wǎng)格下不是諧振。 一個(gè)微波元器件是否諧振應(yīng)該是其自身的屬性 應(yīng)該不是由網(wǎng)格決定的吧
網(wǎng)格數(shù)決定了模型離散化導(dǎo)致的誤差大小,按照你的說(shuō)法,應(yīng)該是大網(wǎng)格下的誤差太大了。你通過(guò)自適應(yīng)網(wǎng)格優(yōu)化看delta S變化就可以看出誤差的大小
的確不是由網(wǎng)格決定的。
但是對(duì)于這個(gè)問(wèn)題的理解是這樣的,不是說(shuō)網(wǎng)格決定了結(jié)果,而是結(jié)果就在那里,你用更加準(zhǔn)確的網(wǎng)格可以得到正確的諧振,而不夠密的網(wǎng)格可能對(duì)模型剖分不夠精確,得不到該得的結(jié)果。
正所謂,你加密或者不加密,微波元器件諧振就在那里,不偏不倚。