Fig.1 shows the schematic view of t

Fig.1 shows the schematic view of the cyclone
separator used in this simulation. This is the simplified
cyclone model corresponding to our previous study[1, 3].
It has a swirl chamber with a radius of 72 mm and a
height of 340 mm. The air flow enters through the top of
the swirl chamber in the tangential direction and is
exhausted through the outlet pipe with a radius of 43.2
Received Sep. 24, 2005
Kazuyoshi MATSUZAKI: Ph. D. 182 Journal of Thermal Science, Vol.15, No.2, 2006
mm, which is coaxially located below the swirl chamber.
As shown in Fig.1, the outlet pipe is installed on the
bottom of the swirl chamber (TYPE-A) or inside the
swirl chamber by 85, 113 and 170 mm from the bottom
of the swirl chamber (TYPE-B, C and D, respectively) in
order to investigate the influence of the position of the
outlet pipe of the swirl chamber on the swirling flows in
the swirl chamber. The length of outlet pipe from the
bottom of the swirl chamber is fixed at 680 mm. The
inlet mean velocity Vin, which is the bulk velocity in the
inlet duct, was about 7.8 m/s corresponding to our
experimental study. The Reynolds number based on the
radius of the swirl chamber and the inlet mean velocity
was about 37,000. The cylindrical coordinate system as
shown in Fig.1 is used and the origin is located at the
center point on the bottom of the swirl chamber.
Numerical Analysis Method
The governing equations are the Navier-Stokes equation
and the equation of continuity for incompressible viscous
flow. They are transformed into cylindrical coordinates
and the spatial filter operation is applied to them. The
finite-difference method is used for the discretization.
The variable is introduced in order to simplify the
discretization at the singular point[9]. A fractional step
method is used for the solution of the pressure. A
Smagorinsky model with van Driest wall damping
function is applied to subgrid scale (SGS) stresses. In
this computation, the Smogorinsky constant Cs is set at
0.10. All spatial terms except for the convective terms
are discretized with 2nd-order central difference scheme.
The 3rd-order modified upstream-biased finite difference
scheme is applied to the convective terms[10]. In the time
integration, the 3rd-order Adams-Bashforth method is
applied to the terms including convection and eddy
viscosity. The 2nd-order Crank-Nicolson method is used
for viscous terms.
A staggered grid system is used for variable array.
Grid spacing in the r and z directions is generated more
finely close to the wall and by a function of hyperbolic
tangent. Uniform grid spacing is adopted in the θ
direction. Fig.2 shows the computational grid for TYPE-B.
As seen in this figure, the whole computational region is
separated into two parts. Block-1and Block-2 is the
region of the swirl chamber and the outlet pipe,
respectively. Grid points for Block-1 are 64×65×128 and
those for Block-2 are 32×65×64 in the r, θ and z
directions, respectively. Non-dimensional time increment
Δt is set at 5×10-4. After the fully developed flow fields
are achieved, time history data of the 200,000 time steps
is taken in order to calculate the time average flow fields.
r
z
θ
Block-1
Block-2
Fig.2 The computational grid for TYPE-B
Neumann conditions are imposed on pressure at the
solid walls. Non-slip conditions are used for velocities on
the walls. The uniform flow is given for the inlet
conditions. Velocities on the outlet boundary are obtained
by solving Orr-Sommerfeld equations. Pressure condition
on the outlet boundary is given by the pressure equation
obtained from substituting the Orr-Sommerfeld equation
into Navier-Stokes equation in the r direction[11].
The motion of a small rigid sphere is described by
the equations of translational motion and rotational
motion[12,13]. The motion of individual particle is treated
by a Lagrangian method. The motions of fluid and
particles are computed simultaneously with the one-way
method. The fluid velocity at a particle position is
calculated by linear interpolation of the velocities at 8
grid points surrounding the particle and the velocity and rotational velocity of particles were integrated by the
3rd-order Adams-Bashforth method and the particle
position by the 2nd-order Crank-Nicolson method.
Initially, 100,000 particles are distributed randomly in
the inlet part of the swirl chamber. Particles are
lycopodium with a density of 700 kg/m3
. The diameters
dp are set at 5 μm in this simulation. Inter-particle
collisions are not considered because the solid volume
fraction is very small. Impulsive equations are solved
when particles collide with the walls [13]. A fully
developed state is used as the initial condition of the
fluid flow. The initial velocity and rotational velocity of
particles are given by the same values as those of the
fluid at each particle position.
The performance for the particle separation of the
cyclone separator is estimated by the collection rate. The
collection rate of particles is calculated by the following
equation.
100
total
ou ×= N
N t η (%) (1)
where Nout is the number of particles which flow out
from the outlet pipe and Ntotal is the total number of
particles mixed into the swirl chamber.

ผลลัพธ์ (ไทย) 1:

ภาพแสดงมุมมองแผนผังตัวอย่างของพายุตัวแยกที่ใช้ในการจำลองนี้ เป็นการง่ายรุ่นไซโคลนที่สอดคล้องกับการศึกษาของเราก่อนหน้านี้ [1, 3]มีหอหมุนและรัศมี 72 มม.และสูง 340 มม. ป้อนอากาศไหลผ่านด้านบนของหมุนในทิศทางที่ tangential chamber และเป็นหมดผ่านท่อร้านและรัศมี 43.2รับกันยายน 24, 2005Kazuyoshi MATSUZAKI: สมุด Ph. D. 182 ของร้อนวิทยาศาสตร์ Vol.15, No.2, 2006มม. ซึ่ง coaxially อยู่ด้านล่างหอหมุนดังแสดงในภาพ ท่อร้านติดตั้งด้านล่างของหอหมุน (ชนิด A) หรือภายในswirl หอ 85, 113 และ 170 มม.จากด้านล่างของหอหมุน (ประเภท B, C และ D ตามลำดับ) ในสั่งการตรวจสอบอิทธิพลของตำแหน่งท่อร้านหอหมุนในการหมุนรอบของไหลในหอหมุน ความยาวของท่อร้านจากด้านล่างของหอหมุนได้คงที่ 680 mm.ทางเข้าของวิน ซึ่งเป็นความเร็วขนาดใหญ่ความเร็วเฉลี่ยในการทางเข้าของท่อ ถูกประมาณ 7.8 m/s ที่สอดคล้องกับของเราศึกษาทดลอง หมายเลขเรย์โนลด์สตามรัศมีของหอหมุนและความเร็วเฉลี่ยทางเข้าของมีประมาณ 37000 ระบบพิกัดทรงกระบอกเป็นแสดงในภาพถูกใช้ และจุดเริ่มต้นอยู่ที่การจุดศูนย์กลางด้านล่างของหอหมุนวิธีการวิเคราะห์เชิงตัวเลข สมการควบคุมสมการของ Navier-สโตกส์และสมการของความต่อเนื่องสำหรับ incompressible ข้นขั้นตอนการ พวกเขาจะเปลี่ยนเป็นพิกัดทรงกระบอกและการอนุรักษ์พื้นที่ตัวกรองไปใช้ ที่วิธีผลต่างจำกัดใช้สำหรับการ discretizationตัวแปรเป็นที่รู้จักเพื่อทำการdiscretization จุด singular [9] ขั้นตอนที่เป็นเศษส่วนใช้วิธีการแก้ไขปัญหาของความดัน Aรุ่น Smagorinsky กับ Driest แวนผนังลดเสียงรบกวนมีใช้ฟังก์ชันกับเส้นตารางย่อยกระบวนสเกล (SGS) เครียด ในนี้คำนวณ Cs คง Smogorinsky ตั้งอยู่ที่0.10 เงื่อนไขทุกพื้นที่ยกเว้นเงื่อนไขด้วยการพามี discretized กับแผน 2 สั่งต่างกลางสั่ง 3 แก้ไขภาพขั้นต้นน้ำความแตกต่างแน่นอนประโยชน์ที่ใช้เงื่อนไขด้วยการพา [10] ในเวลารวม วิธี Adams Bashforth 3 ลำดับคือกับเงื่อนไขรวมถึงการพาและเอ็ดดี้ความหนืด ใช้วิธี Crank Nicolson 2 สั่งสำหรับเงื่อนไขความหนืดระบบกริดที่เหลื่อมกันจะใช้สำหรับตัวแปรแถวลำดับระยะห่างเส้นในทิศทาง r และ z ถูกสร้างขึ้นปิดผนัง และฟังก์ชันไฮเพอร์โบลิของประณีตแทนเจนต์ ระยะห่างระหว่างเส้นตารางเป็นรูปแบบเป็นนำในθทิศทางการ Fig.2 แสดงตารางการคำนวณสำหรับชนิดบีเห็นในรูปนี้ ภูมิภาคทั้งหมดคำนวณเป็นแบ่งออกเป็นสองส่วน 1and บล็อกบล็อก 2ภูมิภาคหอหมุนและท่อร้านตามลำดับ ตารางสำหรับบล็อก 1 คะแนน 64 × 65 × 128 และที่บล็อก 2 มี 32 × 65 × 64 r θ และ zทิศทาง ตามลำดับ เวลาไม่มีมิติเพิ่มขึ้นΔt ตั้งอยู่ที่ 5 × 10-4 หลังจากฟิลด์ขั้นตอนพัฒนาอย่างเต็มมีข้อมูลประวัติทำได้ เวลา 200000 เวลาขั้นตอนจะใช้เพื่อคำนวณเขตข้อมูลกระแสเฉลี่ยเวลาrzΘบล็อก-1บล็อก 2Fig.2 ตารางคำนวณสำหรับประเภท Bเงื่อนไข Neumann มีเก็บความดันที่การผนังทึบ เงื่อนไขการลื่นใช้ในตะกอนในผนัง การไหลสม่ำเสมอให้สำหรับทางเข้าของเงื่อนไขการ ตะกอนบนขอบเขตของร้านจะได้รับโดยแก้สมการ Orr-Sommerfeld เงื่อนไขความดันในร้านขอบเขตถูกกำหนด โดยสมการของแรงดันได้จากการแทนที่สมการของ Orr Sommerfeldลงในสมการ Navier-สโตกส์ในทิศทาง r [11]มีอธิบายการเคลื่อนที่ของทรงกลมแข็งขนาดเล็กโดยสมการ ของการเคลื่อนไหว translational และในการหมุนเคลื่อนไหว [12,13] การเคลื่อนไหวของอนุภาคแต่ละตัวจะถือว่าโดยวิธี Lagrangian เคลื่อนไหวของของเหลว และอนุภาคจะคำนวณพร้อมกันกับทางเดียววิธีการ คือความเร็วของเหลวที่ตำแหน่งที่อนุภาคคำนวณ โดยสอดแทรกเชิงเส้นของตะกอนที่ 8เส้นรอบ ๆ อนุภาค และความเร็ว และความเร็วในการหมุนของอนุภาคถูกรวมเข้าโดยการวิธีสั่ง 3 Adams Bashforth และอนุภาคตำแหน่ง โดยวิธี Crank Nicolson สั่ง 2เริ่มแรก 100000 อนุภาคกระจายแบบสุ่มในส่วนทางเข้าของของท่อหมุน มีอนุภาคlycopodium มีความหนาแน่นของ 700 kg/m3. สมมาตรdp จะตั้งที่ μm 5 ในการจำลองนี้ ระหว่างอนุภาคจะไม่มีพิจารณาตามเนื่องจากปริมาณของแข็งเศษมีขนาดเล็กมาก มีแก้ไขสมการ impulsiveเมื่ออนุภาคชนกับกำแพง [13] A ทั้งหมดรัฐที่พัฒนาแล้วใช้เป็นเงื่อนไขเริ่มต้นของการกระแสของของเหลว ความเร็วเริ่มต้นและความเร็วในการหมุนอนุภาคจะกำหนดค่าเดียวกันเป็นผู้น้ำมันในแต่ละตำแหน่งของอนุภาคประสิทธิภาพการแยกอนุภาคของตัวคาดพายุแบ่งตามอัตราการเรียกเก็บเงิน ที่มีคำนวณอัตราการเรียกเก็บเงินของอนุภาค โดยต่อไปนี้สมการ100ผลรวมซื้อ ou = NΗ t N (%) (1)Nout อยู่จำนวนอนุภาคที่ไหลออกจากร้านท่อ Ntotal เป็นจำนวนอนุภาคที่ผสมเข้าไปในหอหมุน

ผลลัพธ์ (ไทย) 2:

รูปที่ 1 แสดงให้เห็นถึงมุมมองแผนผังของพายุไซโคลน
คั่นที่ใช้ในการจำลองนี้ นี่คือง่าย
แบบพายุไซโคลนที่สอดคล้องกับการศึกษาก่อนหน้านี้ [1, 3].
แต่ก็มีห้องหมุนที่มีรัศมีของ 72 มิลลิเมตรและ
ความสูง 340 มม การไหลของอากาศผ่านเข้าสู่ด้านบนของ
ห้องหมุนในทิศทางวงและ
หมดผ่านท่อระบายมีรัศมี 43.2
ได้รับ 24 กันยายน 2005
Kazuyoshi Matsuzaki: Ph. D. 182 วารสารวิทยาศาสตร์ความร้อน, ปีที่ 15, ฉบับที่ 2, 2006
มิลลิเมตรซึ่งตั้งอยู่ด้านล่างห้อง coaxially หมุน.
ดังแสดงในรูปที่ 1, ท่อระบายติดตั้งอยู่บน
ด้านล่างของห้องหมุน (TYPE-A) หรือภายใน
ห้องหมุนได้ถึง 85, 113 และ 170 มมจากด้านล่าง
ของห้องหมุน (TYPE-B, C และ D ตามลำดับ) ใน
การสั่งซื้อเพื่อตรวจสอบอิทธิพลของตำแหน่งของ
ท่อระบายของห้องหมุนในการหมุนไหลใน
ห้องหมุน ความยาวของท่อระบายจาก
ด้านล่างของห้องหมุนคงที่ 680 มิลลิเมตร
เข้าหมายถึงความเร็วเหล้าซึ่งเป็นความเร็วที่เป็นกลุ่มใน
ท่อทางเข้าเป็นประมาณ 7.8 เมตร / วินาทีของเราสอดคล้องกับ
การศึกษาทดลอง จำนวนนาดส์ขึ้นอยู่กับ
รัศมีของห้องหมุนและทางเข้าหมายถึงความเร็ว
เป็นเรื่อง 37,000 ทรงกระบอกระบบพิกัดเป็น
ที่แสดงในรูปที่ 1 จะใช้และที่มาตั้งอยู่ที่
จุดกึ่งกลางที่ด้านล่างของห้องหมุน.
วิธีการวิเคราะห์เชิงตัวเลข
สมการปกครองเป็นสมการ Navier-Stokes
และสมการของความต่อเนื่องสำหรับความหนืดอัดไม่
ไหล . พวกเขากลายเป็นพิกัดทรงกระบอก
และการดำเนินการกรองอวกาศถูกนำไปใช้กับพวกเขา
วิธีการ จำกัด แตกต่างกันจะใช้สำหรับการ discretization.
ตัวแปรที่ถูกนำมาใช้ในการสั่งซื้อเพื่อให้ง่ายต่อ
เนื่องที่จุดเอกพจน์ [9] ขั้นตอนบางส่วน
วิธีการที่ใช้สำหรับการแก้ปัญหาของความดัน
รูปแบบ Smagorinsky กับผนังตู้ Driest หมาด ๆ
ฟังก์ชั่นที่ใช้กับ subgrid ขนาด (SGS) ความเครียด ใน
การคำนวณนี้ Smogorinsky คง Cs ตั้งไว้ที่
0.10 ทุกแง่เชิงพื้นที่ยกเว้นเงื่อนไขไหลเวียน
จะ discretized กับ 2 คำสั่งรูปแบบที่แตกต่างกันกลาง.
3 สั่งแก้ไขต้นน้ำลำเอียงต่างสืบ
โครงการจะนำไปใช้เงื่อนไขการไหลเวียน [10] ในเวลาที่
บูรณาการ 3 วิธีการสั่งซื้ออดัมส์-Bashforth ถูก
นำไปใช้กับเงื่อนไขรวมทั้งการพาความร้อนและการไหลวน
ความหนืด 2 สั่งซื้อวิธีการหมุน-Nicolson ถูกนำมาใช้
สำหรับคำหนืด.
ระบบกริดเซจะใช้สำหรับตัวแปรอาร์เรย์.
ระยะห่างกริดในทิศทาง R z และถูกสร้างขึ้นมากขึ้น
อย่างประณีตใกล้กับผนังและฟังก์ชั่นของการผ่อนชำระ
สัมผัส ระยะห่างตารางชุดถูกนำมาใช้ในθ
ทิศทาง รูปที่ 2 แสดงให้เห็นตารางการคำนวณสำหรับ TYPE-B.
เท่าที่เห็นในรูปนี้ภูมิภาคคำนวณทั้งหมดจะถูก
แบ่งออกเป็นสองส่วน บล็อกบล็อก 1 และที่ 2 เป็น
พื้นที่ของห้องหมุนและท่อระบาย,
ตามลำดับ จุดกริดหา Block-1 เป็น 64 × 65 × 128 และ
ผู้ที่หา Block-2 มี 32 × 65 × 64 R, θ z และ
ทิศทางตามลำดับ เวลาที่ไม่ได้เพิ่มมิติ
Δtตั้งไว้ที่ 5 × 10-4 หลังจากที่กระแสการไหลของการพัฒนาอย่างเต็มที่
จะประสบความสำเร็จ, ข้อมูลประวัติเวลาในขั้นตอนเวลา 200,000
ถูกนำมาเพื่อที่จะคำนวณระยะเวลาที่กระแสการไหลเฉลี่ย.
อาร์
ซี
θ
บล็อก 1
บล็อก-2
รูปที่ 2 ตารางการคำนวณสำหรับ TYPE-B
เงื่อนไขนอยมันน์ ที่กำหนดไว้ในความดันที่
ผนังทึบ เงื่อนไขกันลื่นใช้สำหรับความเร็วบน
ผนัง การไหลสม่ำเสมอจะได้รับสำหรับเข้า
เงื่อนไข ความเร็วในขอบเขตเต้าเสียบจะได้รับ
โดยการแก้สมการออร์-Sommerfeld สภาพความดัน
ที่ชายแดนเต้าเสียบจะได้รับจากสมการความดัน
ที่ได้จากสมการแทนออร์-Sommerfeld
เป็นสมการ Navier-Stokes ในทิศทาง R [11].
การเคลื่อนที่ของทรงกลมขนาดเล็กและแข็งอธิบายโดย
สมการการเคลื่อนที่การแปลและการหมุน
การเคลื่อนไหว [12,13] การเคลื่อนที่ของอนุภาคของแต่ละบุคคลได้รับการปฏิบัติ
ด้วยวิธีการลากรองจ์ การเคลื่อนไหวของของเหลวและ
อนุภาคจะคำนวณพร้อมกันกับทางเดียว
วิธี ความเร็วของไหลที่ตำแหน่งของอนุภาคจะ
คำนวณโดยสอดแทรกเชิงเส้นของความเร็วที่ 8
จุดรอบอนุภาคและตารางความเร็วและความเร็วการหมุนของอนุภาคที่ถูกบูรณาการโดย
3 วิธีการสั่งซื้ออดัมส์-Bashforth และอนุภาค
ตำแหน่งโดย 2 สั่ง วิธีการหมุน-Nicolson.
ในขั้นต้น 100,000 อนุภาคมีการกระจายแบบสุ่มใน
ส่วนที่ทางเข้าของห้องหมุน อนุภาคมี
Lycopodium ที่มีความหนาแน่น 700
กก. / m3 เส้นผ่าศูนย์กลาง
DP จะถูกตั้งไว้ที่ 5 ไมครอนในการจำลองนี้ ระหว่างอนุภาค
ชนกันจะไม่ถือว่าเป็นเพราะปริมาณของแข็ง
ส่วนมีขนาดเล็กมาก สมห่ามจะแก้ไข
เมื่ออนุภาคชนกับผนัง [13] อย่างเต็มที่
รัฐพัฒนาจะใช้เป็นเงื่อนไขเริ่มต้นของ
การไหลของของไหล ความเร็วเริ่มต้นและความเร็วการหมุนของ
อนุภาคที่จะได้รับจากค่าเช่นเดียวกับบรรดาของ
ของเหลวที่ตำแหน่งของอนุภาคแต่ละ.
ประสิทธิภาพการทำงานสำหรับการแยกอนุภาค
คั่นพายุไซโคลนประมาณจากอัตราการจัดเก็บ
อัตราการจัดเก็บของอนุภาคมีการคำนวณโดยต่อไปนี้
สม.
100
ทั้งหมด
อู× = ไม่มี
? n T η (%) (1)
ที่ Nout คือจำนวนของอนุภาคที่ไหลออกมา
จากท่อระบายและ Ntotal เป็นจำนวนรวมของ
อนุภาคผสม เข้าไปในห้องหมุน

ผลลัพธ์ (ไทย) 3:

วิว”แสดงแผนผังของไซโคลน
คั่นที่ใช้ในการจำลองนี้ นี้เป็นแบบง่าย
ไซโคลนแบบที่สอดคล้องกับการศึกษาก่อนหน้าของเรา [ 1 , 3 ] .
มันมีห้องหมุนด้วยรัศมี 72 mm และ
ความสูง 340 mm อากาศไหลเข้ามาทางด้านบนของห้อง
หมุนในทิศทางแนวและ
หมดผ่านร้านท่อด้วยรัศมี 90
ได้รับ 24 กันยายน ,2005
คาซูโยชิ มัตสึซากิ : วิทยานิพนธ์ดุษฎีบัณฑิต . 182 วารสารเชิงวิทยาศาสตร์ vol.15 2 , 2006
มิลลิเมตร ซึ่งเป็น coaxially ตั้งอยู่ด้านล่างหมุน ห้อง
ดัง” เต้าเสียบท่อที่ติดตั้งบน
ด้านล่างของห้องหมุน ( - A ) หรือภายใน
หมุนห้อง 85 113 และ 170 อืมจากด้านล่าง
ของหมุน ( ห้อง Type C และ D ตามลำดับ )
เพื่อศึกษาอิทธิพลของตำแหน่งของ
ร้านท่อของห้องบนหมุนหมุนหมุนไหล
Chamber ความยาวของท่อขาออกจาก
ด้านล่างของหมุนห้องที่ 680 มิลลิเมตร
หมายถึงความเร็วเข้าวิน ซึ่งมีความเร็วมากใน
เข้าท่อ ประมาณ 7.8 M / S ที่สอดคล้องกับการศึกษาของเรา

เลขเรย์โนลด์ตาม
รัศมีของการหมุนและหอการค้าปากน้ำหมายถึงความเร็ว
ประมาณ 37 , 000 . ระบบพิกัดทรงกระบอกเป็น
แสดงใน”คือการใช้และที่มาตั้งอยู่ที่
ศูนย์จุดที่ด้านล่างของห้องหมุน

วิธีวิเคราะห์เชิงตัวเลขของสมการเป็นสมการ
navier Stokes และสมการความต่อเนื่องสำหรับอัดหนืด
ไหลพวกเขาจะกลายเป็นทรงกระบอกและพิกัด
งานกรองเชิงพื้นที่ที่ใช้กับพวกเขา
ผล เป็นวิธีที่ใช้สำหรับค่า .
ตัวแปรที่เป็นแนะนำเพื่อลดความซับซ้อนในที่จุดเอกพจน์ค่า
[ 9 ] วิธีขั้นเศษส่วน
ใช้โซลูชั่นของความดัน เป็นรูปแบบ smagorinsky
กับรถตู้แห้งหมาดๆ
ผนังฟังก์ชันที่ใช้กับขนาด subgrid ( SGS ) ความเครียด ใน
การคำนวณนี้ smogorinsky คงที่ CS ตั้งอยู่ที่
0.10 . ทุกพื้นที่ยกเว้นเงื่อนไขข้อตกลงการออกแบบ
แบบจุดสั่งซื้อ 2 กลางแตกต่างกับโครงการ เพื่อแก้ไข
3
ลำเอียงเหนือขอบเขตความแตกต่างนี้ใช้กับเงื่อนไขโดย [ 10 ] ในการสั่ง 3 เวลา

วิธี bashforth อดัมส์ใช้กับเงื่อนไข รวมถึงการพาและเอ็ดดี้
ความหนืด 2 วิธีที่ใช้เพื่อเหวี่ยงนิโคลสัน
คำศัพท์หนืด
ระบบกริดเป็นโงนเงนใช้ตัวแปรอาร์เรย์ .
ตารางระยะห่างใน R และ ซีเส้นทางถูกสร้างขึ้นเพิ่มเติม
ให้ชิดกับผนัง และฟังก์ชันไฮเพอร์โบลิกแทนเจนต์
. ช่องว่าง : เครื่องแบบเป็นบุตรบุญธรรมในθ
ทิศทาง fig.2 แสดงตารางการคำนวณสำหรับ type-b.
เท่าที่เห็นในรูปนี้ ภูมิภาค การคำนวณทั้งหมดจะ
แยกออกเป็นสองส่วน block-1and block-2 เป็น
ภูมิภาคของหมุนและหอการค้าร้านท่อ
ตามลำดับ ตารางคะแนน block-1 เป็น 64 × 65 × 128
ผู้ block-2 เป็น 32 × 65 × 64 ใน R , θและ Z
เส้นทางตามลำดับ เพิ่มมิติΔ T
เวลาไม่ตั้งไว้ที่ 5 × 10-4 . หลังจากพัฒนาเต็มที่เขต
ความไหล ,เวลาข้อมูลประวัติของ 200000 เวลาขั้นตอน
ถ่ายเพื่อคำนวณค่าเฉลี่ยของเวลาการไหลของเขตข้อมูล .
r
z



θ block-1 block-2 fig.2 ตารางการคำนวณสำหรับ Type-B
นอยมันน์เงื่อนไขกำหนดความดันที่
ผนังทึบ เงื่อนไขการจัดส่ง ไม่ใช้ความเร็วบน
ผนัง การไหลสม่ำเสมอให้เงื่อนไขเข้า

ความเร็วในขอบเขตที่ได้รับ
ร้านออร์ ซอมเมอร์เฟลด์ โดยการแก้สมการ แรงดันสภาพ
บนเต้าเสียบ ขอบเขตจะได้รับโดยแรงดันที่ได้จากสมการ

เข้าไปแทน ออร์ ซอมเมอร์เฟลด์สมการสมการใน navier สโต R ทิศทาง [ 11 ] .
การเคลื่อนไหวของเล็กแข็งทรงกลมถูกอธิบายด้วย

และสมการของการเคลื่อนไหวสำหรับการหมุนเคลื่อนไหว [ 12 , 13 ‘ ] การเคลื่อนที่ของอนุภาคแต่ละตัวจะถือว่า
โดยระบบวิธี การเคลื่อนที่ของของไหลและอนุภาคจะคำนวณพร้อมกัน
ด้วยวิธี ?

ความเร็วของอนุภาคที่คำนวณโดยการประมาณตำแหน่ง
เชิงเส้นของความเร็วที่ 8
ตารางคะแนนรอบอนุภาคและความเร็วและความเร็วการหมุนของอนุภาคเป็นแบบบูรณาการโดย
สั่งซื้อ 3 วิธี และอนุภาค
bashforth อดัมส์ตำแหน่งโดยวิธีการสั่งซื้อ 2 ข้อเหวี่ยง นิโคลสัน .
ตอนแรก อนุภาคมีการกระจายแบบสุ่มใน 100000
ปากน้ำส่วนหนึ่งของ Swirl Chamber อนุภาค
ยอดจำหน่ายที่มีความหนาแน่นของ 700 kg / m3
ขนาด
DP ตั้งค่าไว้ที่ 5 μ M ในแบบจำลองนี้ ระหว่างอนุภาค
ไม่ถือว่าเป็นเพราะปริมาณของแข็ง
เศษส่วนมีขนาดเล็กมาก สมการห่ามจะหมดไป
เมื่ออนุภาคชนกับผนัง [ 13 ] อย่างเต็มที่
พัฒนารัฐจะใช้เป็นเงื่อนไขเริ่มต้นของ
การไหลของของไหล ความเร็วเริ่มต้นและความเร็วการหมุนของอนุภาคจะถูกกำหนดโดยค่า

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

รวม 100 = N

อุ๊× n t η ( % ) ( 1 )
ที่ออกมาคือจำนวนของอนุภาคซึ่งไหลออกมาจากร้านท่อ และ ntotal

คือจำนวนของอนุภาคผสมลงใน Swirl Chamber