ในพลศาสตร์ของไหลกฎของสโตกส์ให้แรงเสียดทาน – หรือที่เรียกว่าแรงต้าน – ที่กระทำต่อวัตถุทรงกลม ที่เคลื่อนที่ด้วย เลขเรย์โนลด์ ต่ำมาก ในของไหลหนืด^{[ 1 ]} กฎนี้ได้รับการพัฒนาโดยจอร์จ กาเบรียล สโตกส์ในปี พ.ศ. 2494 โดยการแก้ขีด จำกัด การไหลของสโตกส์สำหรับเลขเรย์โนลด์ต่ำของสมการนาเวียร์-สโตกส์^{[ 2 ]}

แรงหนืดบนทรงกลมขนาดเล็กที่เคลื่อนที่ผ่านของเหลวหนืดกำหนดโดย: ^{[ 3 ] [ 4 ]}

${\displaystyle {\vec {F}__{\rm {d}}=-6\pi \mu R{\vec {v}}}$

โดยที่ (ในหน่วย SI ):

• ${\displaystyle {\vec {F}}_{\rm {d}}}$คือแรงเสียดทาน – หรือที่รู้จักกันในชื่อแรงต้านของสโตกส์ – ที่กระทำต่อพื้นผิวสัมผัสระหว่างของเหลวกับอนุภาค ( นิวตัน , กิโลกรัม ms⁻² ⁾ ;
• μ (ผู้เขียนบางท่านใช้สัญลักษณ์η ) คือความหนืดไดนามิก ( ปาสคาลวินาที, กก.ม. ^⁻¹วินาที^⁻¹ );
• Rคือรัศมีของวัตถุทรงกลม (เมตร)
• ${\displaystyle {\vec {v}}}$คือเวกเตอร์ความเร็วของวัตถุ ไม่ใช่ความเร็วการไหลสัมพัทธ์กับวัตถุ (เมตรต่อวินาที) โปรดสังเกตเครื่องหมายลบในสมการ แรงต้านอากาศชี้ไปในทิศทางตรงกันข้ามกับความเร็วสัมพัทธ์: แรงต้านอากาศต้านการเคลื่อนที่

กฎของสโตกส์ตั้งสมมติฐานต่อไปนี้เกี่ยวกับพฤติกรรมของอนุภาคในของเหลว:

• การไหลแบบลามินาร์
• ไม่มีผลกระทบจากแรงเฉื่อย ( เลขเรย์โนลด์ เป็นศูนย์ )
• อนุภาคทรงกลม
• วัสดุที่เป็นเนื้อเดียวกัน (มีองค์ประกอบสม่ำเสมอ)
• พื้นผิวเรียบ
• อนุภาคไม่รบกวนซึ่งกันและกัน

ขึ้นอยู่กับความแม่นยำที่ต้องการ การที่ไม่เป็นไปตามข้อสมมติเหล่านี้อาจจำเป็นหรือไม่จำเป็นต้องใช้แบบจำลองที่ซับซ้อนกว่า ตัวอย่างเช่น สำหรับความคลาดเคลื่อน 10% ความเร็วจะต้องถูกจำกัดไว้ที่ค่า Re < 1

สำหรับโมเลกุลนั้นกฎของสโตกส์ใช้ในการกำหนดรัศมีสโตกส์และเส้นผ่านศูนย์กลางของ โมเลกุล

หน่วยCGSของความหนืดจลน์ได้รับการตั้งชื่อว่า "สโตกส์" ตามผลงานของเขา

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

ความสำคัญของกฎของสโตกส์แสดงให้เห็นได้จากข้อเท็จจริงที่ว่ากฎนี้มีบทบาทสำคัญในการวิจัยซึ่งนำไปสู่รางวัลโนเบลอย่างน้อยสามรางวัล^{[ 5 ]}

กฎของสโตกส์มีความสำคัญต่อการทำความเข้าใจการว่ายน้ำของจุลินทรีย์และอสุจิรวมถึงการตกตะกอนของอนุภาคขนาดเล็กและสิ่งมีชีวิตในน้ำภายใต้แรงโน้มถ่วง^{[ 5 ]}

ในอากาศ ทฤษฎีเดียวกันนี้สามารถนำมาใช้อธิบายได้ว่าทำไมหยดน้ำขนาดเล็ก (หรือผลึกน้ำแข็ง) จึงสามารถลอยอยู่ในอากาศได้ (เป็นเมฆ) จนกระทั่งมีขนาดใหญ่ถึงขนาดวิกฤตและเริ่มตกลงมาเป็นฝน (หรือหิมะและลูกเห็บ) ^{[ 6 ]}สามารถใช้สมการในลักษณะเดียวกันนี้ในการตกตะกอนของอนุภาคละเอียดในน้ำหรือของเหลวอื่นๆ ได้

ที่ความเร็วปลายทาง (หรือความเร็วการตก)แรงส่วนเกินFeเนื่องมาจากความแตกต่างระหว่างน้ำหนักและแรงลอยตัวของทรงกลม (ซึ่งเกิดจากแรงโน้มถ่วง ทั้งคู่ ^{[ 7 ]} ) จะได้รับดังนี้

${\displaystyle F_{e}=(\rho _{p}-\rho _{f})\,g\,{\frac {4}{3}}\pi \,R^{3},}$

โดยที่ (ในหน่วย SI ):

• ρpคือความหนาแน่นมวลของทรงกลม [kg ^m³ ]
• ρfคือความหนาแน่นมวลของของเหลว [กก. ^ลบ.ม. ]
• gคือความเร่งโน้มถ่วง [m/ ^s² ]

เมื่อกำหนดสมดุลแรงF = F และแก้หาความเร็วvจะได้ความเร็วปลายv โปรดทราบว่าเนื่องจากแรงส่วนเกินเพิ่มขึ้นเป็นR ³และแรงต้านของสโตกส์เพิ่มขึ้นเป็นRความเร็วปลายจึงเพิ่มขึ้นเป็นR ²และดังนั้นจึงแปรผันอย่างมากกับขนาดอนุภาคดังแสดงด้านล่าง หากอนุภาคได้รับน้ำหนักของตัวเองเท่านั้นขณะตกในของเหลวหนืด ความเร็วปลายจะเกิดขึ้นเมื่อผลรวมของแรงเสียดทานและแรงลอยตัวบนอนุภาคเนื่องจากของเหลวสมดุลกับแรงโน้มถ่วง พอดี ความเร็ว v [m/s] นี้กำหนดโดย: ^{[ 7 ]}

$${\displaystyle v={\frac {2}{9}}{\frac {\rho _{p}-\rho _{f}}{\mu }}g\,R^{2}\quad {\begin{cases}\rho _{p}>\rho _{f}&\implies {\vec {v}}{\text{ ในแนวตั้งลง}}\\\rho _{p}<\rho _{f}&\implies {\vec {v}}{\text{ ในแนวตั้งขึ้น}}\end{cases}}}$$

โดยที่ (ในหน่วย SI):

• gคือความแรงของสนามโน้มถ่วง [m/ ^s² ]
• Rคือรัศมีของอนุภาคทรงกลม [เมตร]
• ρpคือความหนาแน่นมวลของอนุภาค [กก. ^ลบ.ม. ]
• ρfคือความหนาแน่นมวลของของเหลว [กก. ^ลบ.ม. ]
• μคือค่าความหนืดไดนามิก [กก./(ม.•วินาที)]

ในการไหลแบบสโตกส์ ที่ เลขเรย์โนลด์ต่ำมากเทอมเร่งความเร็วการพาความร้อนในสมการนาเวียร์-สโตกส์จะถูกละเลย จากนั้นสมการการไหลจะกลายเป็น สำหรับการไหลคงที่ที่ไม่สามารถอัดได้ : ^{[ 8 ]}

${\displaystyle {\begin{aligned}&\nabla p=\mu \,\nabla ^{2}\mathbf {u} =-\mu \,\nabla \times \mathbf {\boldสัญลักษณ์ {\omega }} ,\\[2pt]&\nabla \cdot \mathbf {u} =0,\end{aligned}}}$

ที่ไหน:

• pคือความดันของของเหลว (หน่วยเป็น Pa)
• uคือความเร็วการไหล (หน่วยเป็นเมตร/วินาที) และ
• ωคือค่าความหมุน (ในหน่วย s⁻¹ ⁾ซึ่งกำหนดโดย  ${\displaystyle {\boldsymbol {\omega }}=\nabla \times \mathbf {u} .}$

โดยใช้เอกลักษณ์แคลคูลัสเวกเตอร์ บางอย่าง สมการเหล่านี้สามารถแสดงให้เห็นว่าส่งผลให้เกิดสมการลาปลาสสำหรับความดันและแต่ละส่วนประกอบของเวกเตอร์การหมุนวน: ^{[ 8 ]}

${\displaystyle \nabla ^{2}{\boldสัญลักษณ์ {\omega }}=0}$   และ   ${\displaystyle \nabla ^{2}p=0.}$

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

สำหรับกรณีของทรงกลมในกระแสการไหลแบบสม่ำเสมอในระยะไกลการใช้ระบบพิกัดทรงกระบอก( r , φ , z ) จะเป็นประโยชน์ แกนzผ่านศูนย์กลางของทรงกลมและวางตัวในทิศทางเดียวกับทิศทางการไหลเฉลี่ย ในขณะที่rคือรัศมีที่วัดตั้งฉากกับแกนz จุดกำเนิดอยู่ที่ศูนย์กลางของทรงกลม เนื่องจากการไหลมี ความสมมาตรตาม แกน zจึงไม่ขึ้นอยู่กับมุมอะซิมุธφ

ในระบบพิกัดทรงกระบอกนี้ การไหลที่ไม่สามารถอัดได้สามารถอธิบายได้ด้วยฟังก์ชันกระแสสโตกส์ψซึ่งขึ้นอยู่กับrและz : ^{[ 9 ] [ 10 ]}

${\displaystyle u_{z}={\frac {1}{r}}{\frac {\partial \psi }{\partial r}},\qquad u_{r}=-{\frac {1}{r}}{\frac {\partial \psi }{\partial z}},}$

โดยที่u และu คือส่วนประกอบความเร็วการไหลใน ทิศทาง rและzตามลำดับ ส่วนประกอบความเร็วเชิงมุมใน ทิศทาง φเท่ากับศูนย์ในกรณีสมมาตรตามแกนนี้ ฟลักซ์ปริมาตรผ่านท่อที่ล้อมรอบด้วยพื้นผิวที่มีค่าคงที่ψเท่ากับ2 πψและคงที่^{[ 9 ]}

สำหรับกรณีการไหลแบบสมมาตรตามแกนนี้ ส่วนประกอบที่ไม่เป็นศูนย์เพียงอย่างเดียวของเวกเตอร์ความหมุนωคือส่วนประกอบ เชิงมุม φ ω ^{[ 11 ] [ 12 ]}

${\displaystyle \omega _{\varphi }={\frac {\partial u_{r}}{\partial z}}-{\frac {\partial u_{z}}{\partial r}}=-{\frac {\partial }{\partial r}}\left({\frac {1}{r}}{\frac {\partial \psi }{\partial r}}\right)-{\frac {1}{r}}\,{\frac {\partial ^{2}\psi }{\partial z^{2}}}.}$

ตัวดำเนินการลาปลาซที่ใช้กับความหมุนω จะกลายเป็นในระบบพิกัดทรงกระบอกที่มีสมมาตรแกนนี้: ^{[ 12 ]}

${\displaystyle \nabla ^{2}\omega _{\varphi }={\frac {1}{r}}{\frac {\partial }{\partial r}}\left(r\,{\frac {\partial \omega _{\varphi }}{\partial r}}\right)+{\frac {\partial ^{2}\omega _{\varphi }}{\partial z^{2}}}-{\frac {\omega _{\varphi }}{r^{2}}}=0.}$

จากสมการสองสมการก่อนหน้านี้ และด้วยเงื่อนไขขอบเขตที่เหมาะสม สำหรับความเร็วการไหลสม่ำเสมอระยะไกลuใน ทิศทาง zและทรงกลมรัศมีRพบว่าคำตอบคือ^{[ 13 ]}

${\displaystyle \psi (r,z)=-{\frac {1}{2}}\,u\,r^{2}\,\left[1-{\frac {3}{2}}{\frac {R}{\sqrt {r^{2}+z^{2}}}}+{\frac {1}{2}}\left({\frac {R}{\sqrt {r^{2}+z^{2}}}}\right)^{3}\;\right].}$

การหาค่าความเร็วในระบบพิกัดทรงกระบอกและส่วนประกอบต่างๆ มีดังนี้:

${\displaystyle {\begin{aligned}u_{r}(r,z)&={\frac {3Rrzu}{4{\sqrt {r^{2}+z^{2}}}}}\left(\left({\frac {R}{r^{2}+z^{2}}}\right)^{2}-{\frac {1}{r^{2}+z^{2}}}\right)\\[4pt]u_{z}(r,z)&=u+{\frac {3Ru}{4{\sqrt {r^{2}+z^{2}}}}}\left({\frac {2R^{2}+3r^{2}}{3(r^{2}+z^{2})}}-\left({\frac {rR}{r^{2}+z^{2}}}\right)^{2}-2\right)\end{aligned}}}$

วิธีแก้ปัญหาเรื่องความหมุนวนในระบบพิกัดทรงกระบอกมีดังนี้:

${\displaystyle \omega _{\varphi }(r,z)=-{\frac {3Ru}{2}}\cdot {\frac {r}{{\sqrt {r^{2}+z^{2}}}^{3}}}}$

การหาค่าความดันในระบบพิกัดทรงกระบอกมีวิธีการดังนี้:

${\displaystyle p(r,z)=-{\frac {3\mu Ru}{2}}\cdot {\frac {z}{{\sqrt {r^{2}+z^{2}}}^{3}}}}$

การหาค่าความดันในระบบพิกัดทรงกลมมีสูตรดังนี้:

${\displaystyle p(r,\theta )=-{\frac {3\mu Ru}{2}}\cdot {\frac {\cos \theta }{r^{2}}}}$

สูตรของความดันเรียกอีกอย่างว่าศักย์ไดโพลซึ่งคล้ายคลึงกับแนวคิดในไฟฟ้าสถิต

สูตรทั่วไปที่ครอบคลุมกว่า โดยใช้เวกเตอร์ความเร็วระยะไกลแบบใดก็ได้ในพิกัดคาร์ทีเซียนมีดังนี้: ${\displaystyle \mathbf {u} _{\infty }}$${\displaystyle \mathbf {x} =(x,y,z)^{T}}$$${\displaystyle {\begin{aligned}\mathbf {u} (\mathbf {x} )&=\underbrace {\underbrace {{\frac {R^{3}}{4}}\cdot \left({\frac {3\left(\mathbf {u} _{\infty }\cdot \mathbf {x} \right)\cdot \mathbf {x} }{\|\mathbf {x} \|^{5}}}-{\frac {\mathbf {u} _{\infty }}{\|\mathbf {x} \|^{3}}}\right)} _{{\text{conservative: curl=0,}}\ \nabla ^{2}\mathbf {u} =0}+\underbrace {\mathbf {u} _{\infty }} _{\text{far-field}}} _{\text{Terms of Boundary-Condition}}\;\underbrace {-{\frac {3R}{4}}\cdot \left({\frac {\mathbf {u} _{\infty }}{\|\mathbf {x} \|}}+{\frac {\left(\mathbf {u} _{\infty }\cdot \mathbf {x} \right)\cdot \mathbf {x} }{\|\mathbf {x} \|^{3}}}\right)} _{{\text{non-conservative: curl}}={\boldsymbol {\omega }}(\mathbf {x} ),\ \mu \nabla ^{2}\mathbf {u} =\nabla p}\\[8pt]&=\left[{\frac {3R^{3}}{4}}{\frac {\mathbf {x\otimes \mathbf {x} } }{\|\mathbf {x} \|^{5}}}-{\frac {R^{3}}{4}}{\frac {\mathbf {I} }{\|\mathbf {x} \|^{3}}}-{\frac {3R}{4}}{\frac {\mathbf {x} \otimes \mathbf {x} }{\|\mathbf {x} \|^{3}}}-{\frac {3R}{4}}{\frac {\mathbf {I} }{\|\mathbf {x} \|}}+\mathbf {I} \right]\cdot \mathbf {u} _{\infty }\end{aligned}}}$$

${\displaystyle {\boldsymbol {\omega }}(\mathbf {x} )=-{\frac {3R}{2}}\cdot {\frac {\mathbf {u} _{\infty }\times \mathbf {x} }{\|\mathbf {x} \|^{3}}}}$

${\displaystyle p\left(\mathbf {x} \right)=-{\frac {3\mu R}{2}}\cdot {\frac {\mathbf {u} _{\infty }\cdot \mathbf {x} }{\|\mathbf {x} \|^{3}}}}$

ในสูตรนี้ เทอม ที่ไม่คงตัวแสดงถึงสิ่งที่เรียกว่าStokeslet Stokeslet คือฟังก์ชันกรีนของสมการการไหลของ Stokes ส่วนเทอมที่คงตัวนั้นเท่ากับสนามเกรเดียนต์ไดโพลสูตรของความหมุนวนนั้นคล้ายคลึงกับกฎของ Biot–Savart ในแม่เหล็กไฟฟ้า

อีกทางเลือกหนึ่งที่กระชับกว่าคือ สามารถกำหนดสูตรสนามความเร็วได้ดังนี้:

${\displaystyle \mathbf {u} (\mathbf {x} )=\left[\mathbf {I} +\mathrm {H} \left({\frac {R^{3}}{4}}{\frac {1}{\|\mathbf {x} \|}}\right)-\mathrm {S} \left({\frac {3R}{4}}\|\mathbf {x} \|\right)\right]\cdot \mathbf {u} _{\infty },\quad \|\mathbf {x} \|\geq R}$,

โดยที่คือตัวดำเนินการเชิงอนุพันธ์ของเมทริกซ์เฮสเซียน และคือตัวดำเนินการเชิงอนุพันธ์ที่ประกอบขึ้นจากผลต่างของเมทริกซ์ลาปลาเซียนและเมทริกซ์เฮสเซียน ด้วยวิธีนี้จึงเห็นได้ชัดเจนว่าคำตอบประกอบขึ้นจากอนุพันธ์ของศักยภาพคูลอมบ์ ( ) และศักยภาพไบฮาร์มอนิก ( ) ตัวดำเนินการเชิงอนุพันธ์ที่ใช้กับนอร์มเวกเตอร์จะสร้างสโตกส์เล็ต ${\displaystyle \mathrm {H} =\nabla \otimes \nabla }$${\displaystyle \mathrm {S} =\mathbf {I} \nabla ^{2}-\mathrm {H} }$${\displaystyle 1/\|\mathbf {x} \|}$${\displaystyle \|\mathbf {x} \|}$${\displaystyle \mathrm {S} }$${\displaystyle \|\mathbf {x} \|}$

สูตรต่อไปนี้อธิบายเทนเซอร์ความเค้นหนืดสำหรับกรณีพิเศษของการไหลแบบสโตกส์ ซึ่งจำเป็นในการคำนวณแรงที่กระทำต่ออนุภาค ในระบบพิกัดคาร์ทีเซียนเวกเตอร์เกรเดียนต์จะเหมือนกับเมทริกซ์จาโคเบียนเมทริกซ์Iแทน เมท ริก ซ์เอกลักษณ์${\displaystyle \nabla \mathbf {u} }$

${\displaystyle {\boldsymbol {\sigma }}=-p\cdot \mathbf {I} +\mu \cdot \left((\nabla \mathbf {u} )+(\nabla \mathbf {u} )^{T}\right)}$

แรงที่กระทำต่อทรงกลมสามารถคำนวณได้จากการอินทิเกรตเทนเซอร์ความเค้นบนพื้นผิวของทรงกลม โดยที่e แทนเวกเตอร์หน่วยรัศมีของพิกัดทรงกลม :

${\displaystyle {\begin{aligned}\mathbf {F} &=\iint _{\partial V}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\subset \!\supset \;{\boldsymbol {\sigma }}\cdot {\text{d}}\mathbf {S} \\[4pt]&=\int _{0}^{\pi }\int _{0}^{2\pi }{\boldsymbol {\sigma }}\cdot \mathbf {e_{r}} \cdot R^{2}\sin \theta {\text{d}}\varphi {\text{d}}\theta \\[4pt]&=\int _{0}^{\pi }\int _{0}^{2\pi }{\frac {3\mu \cdot \mathbf {u} _{\infty }}{2R}}\cdot R^{2}\sin \theta {\text{d}}\varphi {\text{d}}\theta \\[4pt]&=6\pi \mu R\cdot \mathbf {u} _{\infty }\end{aligned}}}$

${\displaystyle {\begin{aligned}\mathbf {u} (\mathbf {x} )&=-\;R^{3}\cdot {\frac {{\boldsymbol {\omega }}_{R}\times \mathbf {x} }{\|\mathbf {x} \|^{3}}}\\[8pt]{\boldsymbol {\omega }}(\mathbf {x} )&={\frac {R^{3}\cdot {\boldsymbol {\omega }}_{R}}{\|\mathbf {x} \|^{3}}}-{\frac {3R^{3}\cdot ({\boldsymbol {\omega }}_{R}\cdot \mathbf {x} )\cdot \mathbf {x} }{\|\mathbf {x} \|^{5}}}\\[8pt]p(\mathbf {x} )&=0\\[8pt]{\boldsymbol {\sigma }}&=-p\cdot \mathbf {I} +\mu \cdot \left((\nabla \mathbf {u} )+(\nabla \mathbf {u} )^{T}\right)\\[8pt]\mathbf {T} &=\iint _{\partial V}\!\!\!\!\!\!\!\!\!\!\!\!\!\!\subset \!\supset \mathbf {x} \times \left({\boldsymbol {\sigma }}\cdot {\text{d}}{\boldsymbol {S}}\right)\\&=\int _{0}^{\pi }\int _{0}^{2\pi }(R\cdot \mathbf {e_{r}} )\times \left({\boldsymbol {\sigma }}\cdot \mathbf {e_{r}} \cdot R^{2}\sin \theta {\text{d}}\varphi {\text{d}}\theta \right)\\&=8\pi \mu R^{3}\cdot {\boldsymbol {\omega }}_{R}\end{aligned}}}$

แม้ว่าของเหลวจะอยู่นิ่งและทรงกลมจะเคลื่อนที่ด้วยความเร็วค่าหนึ่ง แต่เมื่อพิจารณาจากกรอบอ้างอิงของทรงกลมแล้ว ทรงกลมจะอยู่นิ่งและของเหลวจะไหลในทิศทางตรงกันข้ามกับการเคลื่อนที่ของทรงกลม

• ความสัมพันธ์ของไอน์สไตน์ (ทฤษฎีจลน์)
• กฎทางวิทยาศาสตร์ที่ตั้งชื่อตามบุคคล
• วิถีกระสุน
• สมการลาก
• การวัดความหนืด
• เส้นผ่านศูนย์กลางทรงกลมเทียบเท่า
• การสะสมตัว (ธรณีวิทยา)
• หมายเลขสโตกส์ – ตัวกำหนดผลกระทบเพิ่มเติมของความปั่นป่วนต่อความเร็วการตกขั้นสุดท้ายของอนุภาคในของเหลว^{[ 14 ]}

• Batchelor, GK (1967). บทนำสู่พลศาสตร์ของไหล . สำนักพิมพ์มหาวิทยาลัยเคมบริดจ์. ISBN 0-521-66396-2.
• แลมบ์, เอช. (1994). อุทกพลศาสตร์ (ฉบับที่ 6). สำนักพิมพ์มหาวิทยาลัยเคมบริดจ์. ISBN 978-0-521-45868-9.หนังสือเล่มนี้ตีพิมพ์ครั้งแรกในปี 1879 และฉบับแก้ไขเพิ่มเติมครั้งที่ 6 ตีพิมพ์ครั้งแรกในปี 1932