Drones near buildings: wind that changes with place and time

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

ลมที่ไหลเข้าหาอาคารอาจเร่งเมื่ออ้อมมุม และเกิดบริเวณความเร็วกับทิศทางเปลี่ยนหลังอาคาร คำว่า 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