โดรนใกล้อาคาร: ลมที่เปลี่ยนตามตำแหน่งและเวลา

ลมจุดเดียวไม่ได้แทนทั้งอาคาร

ลมที่ไหลเข้าหาอาคารอาจเร่งเมื่ออ้อมมุม และเกิดบริเวณความเร็วกับทิศทางเปลี่ยนหลังอาคาร คำว่า wake หมายถึงบริเวณการไหลที่ถูกรบกวนด้านท้ายสิ่งกีดขวาง ในงานสำรวจอาคารหรือค้นหาผู้ประสบภัย โดรนจึงอาจได้รับแรงต่างกันแม้อยู่ห่างจากตึกเท่ากัน งาน NASA ที่จำลองอากาศยานใกล้อาคาร ศึกษาลมและปฏิสัมพันธ์กับ wake ขณะที่งานอุโมงค์ลมเรื่องอาคารหลายหลัง วัดตำแหน่งที่ความเร็วและความปั่นป่วนเปลี่ยน

จบบทนี้ คุณจะอธิบายได้ว่าทำไมต้องบันทึก ทิศลม ตำแหน่ง เวลา และแรงบนโดรน ไม่ใช้เพียงค่าความเร็วลมที่สนามบินหรือจุดเดียว ภาพใต้ชื่อบทเป็นภาพแนวคิดของนักศึกษาไทย ไม่ได้แสดงผลวัดลมจริงของตึกในภาพ

เริ่มจากตัวเลขที่คำนวณมือได้

ความดันพลวัต q = ½ρU² ช่วยดูว่าความเร็วลมมีผลต่อสเกลแรงอย่างไร โดย ρ คือความหนาแน่นอากาศ (kg/m³), U คือความเร็วอ้างอิง (m/s) และ q มีหน่วย Pa ถ้ากำหนดเพื่อฝึกว่า ρ = 1.225 kg/m³ ลมเพิ่มจาก 8 เป็น 12 m/s จะได้ q จาก 39.2 เป็น 88.2 Pa หรือ 2.25 เท่า ค่านี้ ไม่ใช่แรงจริงบนโดรน เพราะยังไม่มีพื้นที่อ้างอิง รูปร่าง ท่าทาง ลมที่แต่ละใบพัดได้รับ และการตอบสนองของตัวควบคุม

ในสนามลมรอบตึก U เปลี่ยนทั้งตำแหน่งและเวลา สำหรับผู้เริ่มเรียนให้ดูกราฟ U(t) เป็นความเร็วลมตามเวลา จุดที่กราฟขึ้นลงเร็วคือสิ่งที่ค่าเฉลี่ยอาจซ่อน หากทำ CFD ต้องระบุว่าจะใช้การคำนวณแบบค่ากลางหรือแบบเปลี่ยนตามเวลา โดยผูกกับคำถาม เช่น ต้องการแรงเฉลี่ยหรือแรงกระชาก งานของ NASA แสดงตัวอย่างโดรนสี่ใบพัดใน wake อาคารพร้อมกราฟค่าสัมประสิทธิ์แรงขับที่เปลี่ยนตามเวลาในรายงาน UAM รูปที่ 10 แต่ไม่ได้ให้เกณฑ์ความปลอดภัยสำหรับโดรนเล็กทุกแบบ

จาก Navier–Stokes ถึงแรงบนโดรน

สมการการไหลอธิบายว่าความเร็ว ความดัน และความหนืดเปลี่ยนอย่างไร เมื่อคำนวณสนามลมแล้วจึงหาแรงจากความดันและแรงเฉือนบนผิวโดรน: F = ∫_ผิว(−p n + τ·n)dA โดย p คือความดัน, n คือเวกเตอร์ตั้งฉากออกจากผิว, τ คือความเค้นหนืด และ dA คือพื้นที่ผิวเล็ก ๆ เครื่องหมายของแรงขึ้นกับนิยาม n และว่ากำลังนับแรงที่อากาศกระทำต่อวัตถุหรือกลับกัน ผู้เรียนเริ่มจากภาพทิศแรงก่อน แล้วค่อยเปิดบทสมการหลัก

การจำลองหนึ่งตึกยังต้องกำหนดรูปร่างตึก พื้น ผิวอาคาร โปรไฟล์ลมขาเข้า ขนาดโดเมน กริด และวิธีแทนใบพัด ผลอาจเปลี่ยนเมื่อเปลี่ยนค่าเหล่านี้ จึงต้องตรวจประวัติแรงและความไวต่อกริดตามNASA CFD Analysis Process

กิจกรรมเริ่มบนคอมพิวเตอร์

วาดแผนผังตึกหนึ่งหลังและกำหนดจุด A ด้านหน้าลม, B ที่มุม, C ด้านหลัง วางตารางสำหรับ U เฉลี่ย, U สูงสุดในช่วงที่กำหนด, ทิศลม และแรง/โมเมนต์บนลำ จากนั้นคำนวณ q ของค่าลม 8 และ 12 m/s ตามตัวอย่าง แล้วเปรียบเทียบกับชุด CFD ขนาดเล็กด้านล่าง ซึ่งใช้เงื่อนไขต่างจากลมเมืองจริง

ห้องทดลอง: แยกรอบแก้สมการออกจากเวลาจริง

ดาวน์โหลดชุด CFD พร้อมผลและวิธีรัน ใช้ Python ตรวจผลที่ให้มาบน Windows/Linux ได้ก่อนติดตั้ง solver: รัน python verify.py และ python compare.py --output my-comparison ในโฟลเดอร์ที่แตก ZIP ถ้าต้องการรันใหม่ อ่านขั้นตอน OpenFOAM 2512 สำหรับ Ubuntu 24.04/WSL ใน README

เริ่มจากก้อนสี่เหลี่ยมกว้างและสูง 0.4 m ในช่องลม 4 × 2 × 1.6 m พื้นกับก้อนหยุดนิ่งและเป็น no-slip ด้านข้างกับด้านบนเป็น symmetry ลมขาเข้าเป็นค่าคงที่ 0.2 หรือ 0.4 m/s ของไหลมี ν = 0.01 m²/s ทำให้ Re = UH/ν = 8 หรือ 16 จึงเป็นการไหลหนืดแบบ laminar เพื่อฝึกคำนวณ ไม่มีโดรนอยู่ในโดเมน และแรงที่รายงานเป็นแรงบนก้อนอาคาร

ซ้าย: SIMPLE หาคำตอบคงตัว โดยแกนนอนคือ iteration ขวา: PIMPLE ติดตามการเปลี่ยนตามเวลาหลังเริ่มการไหล โดยแกนนอนคือวินาที เป็นผลคำนวณจริงจากคนละวิธีในชุดฝึก เปิดกราฟขนาดเต็ม

กริดสองระดับมี 12,736 และ 101,888 เซลล์ ตามลำดับ เมื่อคำนวณแบบเปลี่ยนตามเวลา จะเพิ่ม Δt หรือก้าวเวลา: กรณีลม 0.2 m/s บนกริดหยาบ รันถึง t = 4 s ด้วย Δt = 0.02 และ 0.01 s ได้ Fx = 0.01698087 และ 0.01700083 N ต่าง 0.117% ของค่าจากก้าวเวลาละเอียด

ข้อสังเกต: ผลต่างของก้าวเวลาน้อยไม่ได้ลบผลต่างของกริด และ t = 4 s ก็ยังไม่ใช่หลักฐานว่าการไหลคงตัวแล้ว การเปรียบเทียบต้องใช้ตำแหน่ง กริด และเวลาจริงเดียวกัน ตามหลัก NASA เรื่องความคลาดเคลื่อนจากกริดและก้าวเวลา ยังไม่มีการจำลองความปั่นป่วนหรือลมกระโชกในชุดนี้

ลองทำ: เปิด results/urban-dt002/forces.csv และ urban-dt001/forces.csv ใช้แถวที่ step = 4 ทั้งสองไฟล์ แล้วคำนวณ 100 × |Fx₁ − Fx₂| / |Fx₂| ตรวจได้ประมาณ 0.117% ส่วนไฟล์ urban-v1-r1/forces.csv ใช้ step เป็น iteration ห้ามตีความว่าเป็นเวลา 4 s

การนำไปศึกษาอาคารจริงต้องเพิ่มโปรไฟล์ลมขาเข้า ความขรุขระ การไหลปั่นป่วน การตรวจขนาดโดเมน และข้อมูลวัดที่ตรงบริบท ส่วนแรงบนโดรนต้องมีรูปทรงหรือแบบจำลองโดรนเพิ่ม ไม่ใช่นำแรงบนก้อนนี้ไปใช้แทน

ต่อยอดที่ การตรวจและส่งต่อแรงสู่ DroneSim: แรงบนอาคารและค่าความเร็วสามจุดยังใช้เป็นแรงลำหรือสนามลมตลอดวิถีไม่ได้ ลองอ่านเหตุผลในรายงานความพร้อมก่อนเลือกข้อมูลเข้า simulator

แบบฝึกพร้อมแนวคำตอบ

โจทย์: ลมตัวอย่างเพิ่มจาก 8 เป็น 16 m/s ที่ ρ เดิม q เพิ่มกี่เท่า? ตัวเลขนี้ใช้สรุปแรงบิดบนโดรนใกล้ตึกได้เลยหรือไม่?

แนวคำตอบ: q เพิ่ม 4 เท่า เพราะ U เพิ่ม 2 เท่า แต่ยังบอกแรงบิดไม่ได้ ต้องมีการกระจายลมบนผิวและแต่ละใบพัด ระยะแขนจากจุดอ้างอิง และการเปลี่ยนตามเวลา

ต่อไป: ผลของใบพัดต่อลำตัวและสิ่งบรรทุก

บรรณานุกรม

• NASA — Simulated Tiltrotor Operation Near a Building
• NASA — Wakes from Arrays of Buildings
• NASA — Accessibility for UAM Vehicles, Figure 10
• NASA — CFD Analysis Process
• NASA — Examining Spatial Convergence
• OpenCFD — Forces

ตรวจแหล่งต้นทาง 26 กันยายน 2569