> อัปเดต 28 กันยายน 2026: ทั้งสี่กริดจบ6000แล้ว มีเพียง N32 ผ่านทุก case gate; N30/N36/N44 ไม่ผ่าน residual จึงงด qualified sensitivity/GCI ดูรายงานสุดท้าย `ms1101-family-final-2026-09-28.md` ข้อความด้านล่างเป็นบันทึกสถานะในวันที่27

# ผลคัดกรองกริดชั้นใกล้ผิว MS1101: protocol v3

วันที่ 27 กันยายน 2026 · สถานะ: ผ่าน construction screen 4 ระดับ; solve จบ N30/N32 และกำลังรัน N36/N44

## เหตุที่ต้องลงทะเบียน protocol ใหม่

ชุดแรก N=24/32/44 ไม่ได้สามกริดที่ผ่านครบ เพราะ N=24 มีเซลล์ที่ตรวจพบจากผิวเป็นศูนย์ สร้างชั้นไม่ได้ และ `checkMesh` ไม่ผ่าน (max skewness 4.517) หลังลงทะเบียน v2 ล่วงหน้า ได้ทดลอง N=28 โดยคงเกณฑ์เดิม แต่แม้มีการสร้างชั้นบน 96.02% ของหน้า `checkMesh` ยังพบ skewness สูงสุด 4.737525193 และหน้าที่ skew สูง 4 หน้า กริดนี้จึงไม่ผ่านเช่นกัน

ตาม protocol v2 เราหยุดแล้วบันทึกกรณีที่ไม่ผ่าน จากนั้นลงทะเบียน [protocol v3](ms1101-layer-family-protocol-v3-2026-09-27.md) ก่อนสร้างกริดชุดถัดไป ไม่มีการปรับเกณฑ์ผ่าน/ไม่ผ่านย้อนหลัง

## วิธีคัดกรองที่ตรึงไว้

สร้างกริดใหม่จากเคส MS1101 ต้นทางเดียวกัน โดยตรึงรูปทรง โดเมนลูกบาศก์ด้านละ 0.84 m, surface/rotor/wake refinement 5/3/2, วิธีขยับผิว และ mesh quality controls เปลี่ยนเฉพาะจำนวนเซลล์พื้นฐาน N เป็น 30, 32, 36, 44 ต่อด้าน ระยะพื้นฐานคือ `h=0.84/N`; ความสูงชั้นแรกใช้ `t₁=50 µm×32/N`; ขอ 8 ชั้นด้วยอัตราขยาย 1.2 และ displacement motion solver

กริดต้องผ่านทั้ง `checkMesh`, ไม่เกิน 3,000,000 เซลล์ และสร้างชั้นอย่างน้อยบางส่วนบนไม่น้อยกว่า 95% ของหน้าผิวที่ร้องขอ หากกรณีใดไม่ผ่านจะเก็บผลไว้และไม่เปลี่ยนเกณฑ์

## ผล mesh construction

| N | เซลล์ | หน้าที่สร้างชั้น / หน้าที่ร้องขอ | ความครอบคลุมตามจำนวนหน้า | ชั้นเฉลี่ยจาก 8 | checkMesh | gate |
|---:|---:|---:|---:|---:|---|---|
| 24 · v1 | 100,682 | 0 / 553 | 0.00% | ไม่มี | ไม่ผ่าน; skewness 4.517 | ไม่ผ่าน |
| 28 · v2 | 286,141 | 11,519 / 11,996 | 96.02% | 4.58 | ไม่ผ่าน; skewness 4.738 | ไม่ผ่าน |
| 30 · v3 | 340,769 | 13,089 / 13,662 | 95.81% | 4.74 | ผ่าน | ผ่าน |
| 32 · v3 | 384,364 | 14,818 / 15,425 | 96.06% | 4.81 | ผ่าน | ผ่าน |
| 36 · v3 | 529,897 | 19,042 / 19,731 | 96.51% | 5.15 | ผ่าน | ผ่าน |
| 44 · v3 | 921,354 | 28,340 / 29,302 | 96.72% | 5.64 | ผ่าน | ผ่าน |

ความครอบคลุมนี้นับจำนวนหน้า ไม่ได้ถ่วงน้ำหนักตามพื้นที่ และหน้าใดที่ได้ชั้นบางส่วนก็ถูกนับ ไม่ได้แปลว่าผิวบริเวณนั้นมีครบ 8 ชั้น ดังนั้น N=30/32/36/44 ผ่านเพียงขั้นสร้างกริด ไม่ใช่ผล y⁺ หรือแรงที่ลู่เข้าแล้ว

ภาพสรุปแสดงจำนวนเซลล์และ coverage พร้อมผล gate ของทุก protocol:

![กริด N=24 และ N=28 ไม่ผ่าน checkMesh; กริด N=30, 32, 36 และ 44 ผ่านเฉพาะขั้น construction](/diagrams/aerodynamics/ms1101-global-layer-family.png)

## ขั้น solve และข้อจำกัดการสรุป

ผลระหว่างดำเนินการ: N30 จบที่ 6,000 แต่ residual ผ่าน 0/100 รอบท้าย (สูงสุด 5.029×10⁻⁵) จึงไม่ผ่าน case gate; N32 จบและผ่านทุก case gate โดย residual ผ่าน 100/100 รอบ (สูงสุด 2.458×10⁻⁶) N36/N44 ยังรันอยู่ อ่านค่าแรง สถานะ checkpoint และข้อจำกัดใน [บันทึกการรันและตรวจรับ](ms1101-layer-family-execution-notes-2026-09-27.md) ยังไม่มี qualified grid sensitivity หรือ GCI ของทั้งชุด

ตาม protocol v3 กริดทั้งสี่กำลังแก้สมการด้วย OpenFOAM 2512 `simpleFoam`/steady SIMPLE, Open MPI 4.1.6 จำนวน 8 ranks จาก fields เริ่มต้นที่ iteration 0 ถึง 6,000 โดยใช้เกณฑ์เดียวกันทุกกริด ต้องตรวจ `End` ที่ iteration 6000, initial residual ทั้งหกสมการต่ำกว่า 1e-5 ครบทุก iteration ในช่วงท้าย 100 รอบ, thrust และ resisting torque มี half-window mean change <1% และ peak-to-peak <2%; y⁺ ต้องเป็นข้อมูลที่ iteration 6000 หาก invocation หนึ่งชนเพดานเวลา ให้ต่อจาก checkpoint ที่เขียนครบด้วย `startFrom latestTime` โดยคง input ทางตัวแก้สมการและเกณฑ์ทั้งหมด; logs, checkpoints, เวลาเริ่ม/จบ และ SHA-256 แยกต่อช่วงตาม [protocol addendum](ms1101-layer-family-protocol-v3-checkpoint-continuation-2026-09-27.md) จนกว่าจะถึงรอบ 6000 ผล solver จะถูกสรุปในรายงานนี้หลังจบครบทุกกรณี

แม้ครบสามหรือสี่ solver runs แล้ว จะคำนวณ Richardson/GCI ต่อเมื่อแรงขับและแรงบิดแสดงลำดับ monotonic และข้อมูลรองรับ asymptotic range โดยคำนึงถึงอัตราส่วน h ที่ไม่เท่ากัน หากไม่รองรับจะรายงาน grid sensitivity และงด GCI การศึกษานี้ไม่แทน physical validation กับ ENOLA ซึ่งยังต้องใช้ metadata สภาวะทดสอบ uncertainty รายจุด และหลักฐานยืนยัน specimen

## แหล่งอ้างอิง

- [NASA NPARC: Examining Spatial (Grid) Convergence](https://www.grc.nasa.gov/www/wind/valid/tutorial/spatconv.html)
- [Celik et al. (2008), Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications](https://doi.org/10.1115/1.2960953)
- [OpenCFD: Layer addition in snappyHexMesh](https://doc.openfoam.com/2212/tools/pre-processing/mesh/generation/snappyhexmesh/layers/)