8月的最后一天了
之前的目標(biāo)和承諾也都沒(méi)有實(shí)現(xiàn)
哎每天得過(guò)且過(guò)啊
開(kāi)始共工作吧
把100%轉(zhuǎn)速各個(gè)工況的堵塞系數(shù)都處理掉 10:03
很簡(jiǎn)單啊把cst復(fù)制過(guò)去然后執(zhí)行一下cse此疹，五分鐘吧（如果沒(méi)有意外的話還是用了7分鐘誒贝淤，不過(guò)做好了 10:10
接下來(lái)把80%轉(zhuǎn)速的也處理掉吧處理80%轉(zhuǎn)速各工況的堵塞系數(shù) 10:10
預(yù)計(jì)時(shí)間就一些赔桌，10分鐘吧已經(jīng)處理好了甘耿，比想象中快 10:19
不過(guò)有個(gè)問(wèn)題砾赔，閾值要不要調(diào)旦事？
看了下10.5kgs的，好像沒(méi)啥問(wèn)題
但是看8.4kgs的族吻，發(fā)現(xiàn)如果閾值不變的話，已經(jīng)有部分虧損區(qū)域到75%葉高以下的位置了珠增，這就很蛋疼
那么接下來(lái)看一下80%轉(zhuǎn)速各個(gè)工況的閾值看了一下超歌，其實(shí)5作為閾值是挺合適的 10:22
但是有個(gè)問(wèn)題，在8.4kgs工況蒂教，堵塞區(qū)域?qū)嶋H上下降到比較低的位置巍举，比如低于75%葉高，但是與此同時(shí)較高葉高位置的堵塞值確實(shí)小于5的凝垛，就是說(shuō)有比較大的“island”懊悯。這種情況如何處理呢蜓谋？比較費(fèi)解
我現(xiàn)在說(shuō)不來(lái)怎么搞，但是可以把圖畫(huà)出來(lái)看一下就知道了畫(huà)堵塞的圖炭分，先處理軸向坐標(biāo) 10:25
就取95%葉高的坐標(biāo)好了
一來(lái)之前處理過(guò)
二來(lái)也比較折衷桃焕，在葉頂和75%葉高中間（其實(shí)也沒(méi)有很“中間”
三來(lái)跟I2O的圖能夠?qū)ι献鴺?biāo)啊
所以就不重新生成坐標(biāo)了收集完100N和080N的數(shù)據(jù) 10:33
處理完坐標(biāo)和數(shù)據(jù)，準(zhǔn)備開(kāi)始做dat文件 10:40
忘了除以面積了捧毛，計(jì)算下出口面積观堂，就0.75的好了 10:42
處理完了，量級(jí)合理 10:48
出口處B≈0.1呀忧，和Khalid博士論文中高速風(fēng)扇出口B值相當(dāng)繼續(xù)做dat文件师痕，出lay 10:48
光壁堵塞的后處理處理完了 10:57
這個(gè)小時(shí)之內(nèi)把堵塞的描述完成掉嘗試把壁面的堵塞和泄漏堵塞區(qū)分開(kāi)，好像不好做到 11:08
想把壁面堵塞和泄漏造成的堵塞區(qū)域分開(kāi)
先是把靠近端壁的部分去掉了而账，只要到98%葉高以下的
然后提高閾值胰坟。提高閾值確實(shí)能去掉很大一部分泄漏流造成的堵塞區(qū)域
但是還是存在一些高閾值的區(qū)域，所以這個(gè)衡量沒(méi)法太準(zhǔn)確泞辐。
先放一放笔横，是在不能定量說(shuō)明，后期有機(jī)匣處理結(jié)果的時(shí)候定性說(shuō)明也好铛碑。上午全力學(xué)了一個(gè)小時(shí)狠裹，休息下吧 11:10
休息了倆小時(shí) 13:14
這么看來(lái)真的學(xué)習(xí)時(shí)間太短了
上午學(xué)了一小時(shí)，中午休息兩小時(shí)汽烦，我真的不夠努力啊
下午敢學(xué)四小時(shí)嗎涛菠？
有點(diǎn)夠嗆，先學(xué)兩小時(shí)再說(shuō)開(kāi)始下午的學(xué)習(xí) 13:15
就把堵塞的東西寫(xiě)掉啊
怎么寫(xiě)呢撇吞，單獨(dú)這樣寫(xiě)嗎俗冻，感覺(jué)有點(diǎn)單薄
能不能和前面的流場(chǎng)分析結(jié)合在一起
如果要的話我要在哪里寫(xiě)呢
其實(shí)不太好我，我前面一個(gè)cp分布牍颈，熵分布迄薄，激波結(jié)構(gòu)分析，三維泄漏流分析
堵塞放在哪里比較合適呢煮岁？
照理說(shuō)應(yīng)該和泄漏流發(fā)展放在一起
但是泄漏流發(fā)展里面以來(lái)面是斜的讥蔽，二來(lái)視角是斜的，曲線也疊加不進(jìn)去
要不就分開(kāi)寫(xiě)好了画机， 沒(méi)關(guān)系的么
要不要做堵塞的非定常分布冶伞？
好像沒(méi)什么必要
出堵塞我主要還是想著加上機(jī)匣處理之后可以做對(duì)比，放在前面其實(shí)沒(méi)有特別分析的點(diǎn)
要說(shuō)沒(méi)有也不好這么說(shuō)步氏，100N和080N的分布特點(diǎn)就不太一樣响禽，不過(guò)先不管這個(gè)了
先寫(xiě)起來(lái)開(kāi)始分析光壁堵塞 13:27
有個(gè)疑問(wèn)，我出口堵塞量下降，會(huì)不會(huì)是因?yàn)槲姨幚矸椒ǖ膯?wèn)題
因?yàn)榇_實(shí)可以觀察到在虧損區(qū)域的里面存在著大大小小的island芋类，有些其實(shí)是挺大
但是我沒(méi)有把這些算進(jìn)去隆嗅，單純是用閾值來(lái)計(jì)算的，所以會(huì)不會(huì)是這里出現(xiàn)了問(wèn)題侯繁？
這個(gè)值對(duì)最終結(jié)構(gòu)有多大影響胖喳？
我要確認(rèn)一下重新考察堵塞量 13:36
先觀察各個(gè)工況出口的虧損區(qū)域
然后在對(duì)工況C重新處理，閾值設(shè)低巫击，把island的區(qū)域包含進(jìn)去禀晓，然后來(lái)計(jì)算。觀察各個(gè)工況的虧損情況 13:37
把各工況的視角和顯示配置都調(diào)了 13:46
不可否認(rèn)坝锰，兩個(gè)轉(zhuǎn)速工況C的分布都是那種粹懒，有很大的Island
試試看把缺的那塊補(bǔ)上然后執(zhí)行cse看看感覺(jué)不太對(duì)，我就應(yīng)該提出這個(gè)情況啊顷级，這其實(shí)是流動(dòng)的一個(gè)變化 13:47
近失速工況下堵塞并不是完全堵塞整個(gè)葉高范圍凫乖，實(shí)際上在較低葉高的堵塞（出口）反而更厲害
這是為什么呢
好困啊我要睡會(huì)兒睡會(huì)兒 13:50
睡醒了，剛好15分鐘 14:05
接著處理 14:06
我覺(jué)得我可以說(shuō)明啊弓颈，為什么出口堵塞就一定是變大呢
試試就是虧損區(qū)域雖然范圍跟大了帽芽，到更低的區(qū)域了，但是堵塞面積實(shí)際上變小了-
看了軸向速度 14:18
軸向速度分布看起來(lái)翔冀，其實(shí)工況C中低閾值的區(qū)域反而速度更低呢导街，對(duì)比一下看：
9.8kgs_BA0.75_0-100.png
這圖看起來(lái)就超明顯了，左邊是速度軸向速度分布纤子，右邊是密流梯度
密流梯度里那上部那塊值比較低的區(qū)域搬瑰，在速度分布里依然是低速區(qū)域
甚至是速度更低的
這是什么原因呢？是我梯度的算法不對(duì)嗎控硼？應(yīng)該分別對(duì)r和θ求梯度泽论？
我覺(jué)得沒(méi)必要深究這塊，但是現(xiàn)在要趕緊把這個(gè)處理掉 想辦法把那塊缺的補(bǔ)上 14:28
考慮到軸向位置比較靠前的面卡乾，閾值又不能定得太低
因?yàn)榍懊娴拿嬉磴玻潛p區(qū)域是整個(gè)在75%葉高以上的
如果閾值太低，虧損區(qū)域會(huì)明顯增大
我現(xiàn)在看一下前兩兩個(gè)工況幔妨，如果前兩個(gè)工況設(shè)一個(gè)閾值OK鹦赎，只是工況C不行的話， 考慮手動(dòng)來(lái)設(shè)置確認(rèn)完畢误堡，其他兩個(gè)工況很好 14:42
把100和80轉(zhuǎn)速的AB工況都看了一下钙姊，用5做閾值很好，很明確地指示出了堵塞區(qū)域現(xiàn)在找不到一個(gè)合適的通用的辦法 14:54
先不要一個(gè)勁看了埂伦，沒(méi)有意義，我可以先手動(dòng)看一下0.75出口截面在調(diào)小閾值之后得到的Ab跟之前區(qū)別大不大思恐，如果根本區(qū)別都不大沾谜，就沒(méi)有必要糾結(jié)了事實(shí)證明差老多了 14:58
按照原來(lái)的處理方法膊毁，出口Ab是5e-4量級(jí)，現(xiàn)在是8e-4量級(jí)……
B則直接從0.119上升到0.182
所以還是要改可能確實(shí)梯度做的有問(wèn)題 15:12
剛才試著用兩個(gè)r和θ兩個(gè)方向梯度相加的新梯度來(lái)作為閾值
很好用基跑，原來(lái)上面比較低的值囊括在內(nèi)了也不是特別好用 15:36
考慮到前面的面的話婚温，判定的虧損區(qū)域似乎會(huì)大一些
暫且不管這個(gè)，用這個(gè)方法輸出一條線來(lái)看看
閾值還是5媳否，設(shè)成10的話栅螟，下沿會(huì)高一些，但是的island又會(huì)出現(xiàn)-
出了一組Ck≥5的但是結(jié)果不太好 15:45
現(xiàn)在不止是出口的B提高了篱竭，40%前的B也顯著提高力图，然后40%這個(gè)地方變成一個(gè)明顯的拐點(diǎn)
又出了一個(gè)Ck≥10的，和5的結(jié)果相近
圖片1.png
這說(shuō)明閾值定義方法對(duì)虧損區(qū)域的定義影響很大 去泡個(gè)咖啡 15:52
碰到問(wèn)題之后推進(jìn)的速度會(huì)慢很多回來(lái)了掺逼，接著干活吧 16:18
把用新的準(zhǔn)則把另兩個(gè)工況也處理下試試發(fā)現(xiàn)問(wèn)題了吃媒，前部的虧損區(qū)域絕對(duì)高估了 16:30
如果說(shuō)后部的堵塞，兩種梯度的計(jì)算方法得到的結(jié)果相近吕喘，那么前部的結(jié)果差別是在太大了
就不放圖了赘那，之前的方法很明顯，進(jìn)口處的堵塞是很低的氯质，僅在上端壁附近有一些附面層的堵塞
而后一種方法募舟， 基本無(wú)法準(zhǔn)確辨別主流區(qū)，很大一片都是有較高值的
甚至在在前緣上游都顯示混亂的虧損區(qū)
實(shí)際上這是不可能的啊闻察，還沒(méi)進(jìn)通道怎么會(huì)有大面積堵塞呢
所以我覺(jué)得可能是激波的存在使得梯度計(jì)算錯(cuò)誤估計(jì)了虧損區(qū)域拱礁。現(xiàn)在還是用以前的梯度方法把后部的結(jié)果處理一下 16:51
其實(shí)應(yīng)該挺快的，畢竟只是一個(gè)工況的一半
重新找一個(gè)合適的閾值蜓陌，把后半部分替換掉就好了或者用組合的方法觅彰，反正就是判斷的虧損區(qū)域嘛 16:57
后半部分靠近出口的位置用新的梯度定義？確定處理方法 17:04
用后一種梯度處理方法
然后將閾值提高到20
這樣會(huì)去掉一些上部的island
但是我可以不回去钮热，因?yàn)檫@些island范圍比較高填抬，不會(huì)和下面的主流混起來(lái)
確定了85%葉高以上的區(qū)域-
處理出來(lái)效果不錯(cuò) 17:25
圖片2.png
把0.45前面的（含0.45）用老的梯度方法，0.50后面的（含0.50）用新的梯度方法隧期，得到組合的結(jié)果
這個(gè)結(jié)果看起來(lái)就合理得多了
事實(shí)證明虧損區(qū)域的判定十分重要 現(xiàn)在對(duì)比一下11.5kgs用兩種方法處理出來(lái)的效果 17:31
在11.5kgs中激波位置更靠近下游飒责，所以第二種梯度方法有問(wèn)題的范圍更大
BA0.45跟梯度1方法就很不一樣，發(fā)生區(qū)別的區(qū)域很明顯是激波影響下速度梯度較大的區(qū)域-
結(jié)論是：工況AB非但不需要用第二種梯度方法仆潮，而且用了是錯(cuò)的 17:45
存在兩個(gè)問(wèn)題：- 一是在弦長(zhǎng)前半部分宏蛉，由于激波的影響，第二種梯度方法判斷的虧損區(qū)域過(guò)大性置，與激波形態(tài)有關(guān)
- 二是即使在弦長(zhǎng)后半部分拾并，不再受到激波影響了，這時(shí)虧損區(qū)域中并不存在Island，根本用不著補(bǔ)充Iso Clip 4嗅义，補(bǔ)充了反而是錯(cuò)的屏歹，因?yàn)楹竺娴奶潛p區(qū)域部分高于85%葉高的
現(xiàn)在只需要把80%轉(zhuǎn)速工況C處理掉就行了 17:49
好了80%的也處理好了 18:30
本來(lái)還想去踢球，算了吧還是把堵塞的東西弄弄好有點(diǎn)餓了去稍微買(mǎi)點(diǎn)東西墊墊肚子 18:30
吃了個(gè)牛肉烤餅歇了一會(huì)兒 19:18
開(kāi)始重新出堵塞系數(shù)的圖 19:18
把數(shù)據(jù)一替換就行了做好了新的堵塞系數(shù) 19:26
接下來(lái)把描述寫(xiě)一下之碗，換成新的梯度
