where the partial density of pure b

where the partial density of pure brine and the dissolved CO2 mass fraction are both functions of temperature, pressure and salt mass fraction, calculated by Battistelli et al. [38] and Spycher and Pruess [5], respectively. The partial density of dissolved CO2 in Eq. (6) is calculated from:
q2 ¼
MCO2
Vu
ð7Þ
where Vu is a function of temperature and independent of pressure, given by García [39]. In this work, T is a constant, accordingly Vu = 34.92  106 m3/mol is used. MCO2 is equal to 44.0  103 kg/mol. When the two values are substituted into Eq. (7), q2 = 1260 kg/m3 is obtained. In the above, Eq. (6) shows that the density of aqueous mixture increases when CO2 dissolves into the brine at the reservoir conditions. Eqs. (2)–(7) constitute the fundamental governing equations for
numerical simulations of the flow in porous media for solubility trapping, which are a set of coupled nonlinear equations and could be solved by Newton–Raphson iteration.
2.2. Numerical method
Eqs. (2) and (3) are discretized in time using an implicit finite difference scheme and in space using an integral finite difference method. The discretization method, the spatial and temporal resolutions used in the simulations are described as follows. For the discretization, the accumulation terms are discretized by introducing appropriate volume averages:
@ @t
Z
Vn
/qXidVn ¼ /tþDt
n qtþDt
n XtþDt
i;n  /t
nqt
nXti
;n
 

Vn
Dt ð8Þ
where t + Dt represents the new time step. The fluxes of convective and diffusive terms are approximated as a discrete sum of averages over all surfaces of meshes, Z
Sn
ðqXiqÞ  ndSn ¼
X
m
qtþDt
nm XtþDt
i;nm qtþDt
nm Snm ð9Þ
Z
Sn
ð/DqrXiÞ  ndSn ¼
X
m
/DqtþDt
nm rXtþDt
i;nm Snm ð10Þ
where these flux terms are treated as fully implicit, given by the values at the new time step. The Darcy’s law is discretized in the following way [40], qtþDt
nm ¼ 
K
l
 
nm
PtþDt
n  PtþDt
m
Dnm
þ qtþDt
nm g cos a
!
ð11Þ
where a is the intersection angle between gravitational acceleration and the line segment from mesh m to n with rotation direction from g to the line segment clockwise as indicated in Fig. 2. In Eq. (11), cosa can be written as follows:
cos a ¼
Z2  Z1
Dnm
ð12Þ
where Z1, Z2 are the respective grid block elevations. The Darcy’s law between mesh m and n is then written as
qtþDt
nm ¼ 
K
l
 
nm
PtþDt
n  PtþDt
m
Dnm
þ qtþDt
nm g
Z2  Z1
Dnm
!
ð13Þ
where qnm, Xi,nm, (K/l)nm are the harmonic averages between meshes m and n, given by, qnm ¼
D1qm þ D2qn
Dnm
ð14Þ
Xi;nm ¼
D1Xi;m þ D2Xi;n
Dnm
ð15Þ
K
l
 
nm
¼
D1ðK=lÞm þ D2ðK=lÞn
Dnm
ð16Þ
where Dnm is given as:
Dnm ¼ D1 þ D2 ð17Þ
2.3. The consideration of resolutions
Evolution and interaction of flow and transport over a broad range of space and time scales are essential features of the solubility trapping. Therefore temporal and spatial resolutions have to be adequate to resolve the essential flow features. In order to capture the early stage diffusive behaviours and obtain the onset time of convection accurately, a small time step of 0.01 s is used. At the same time, due to the long-term behaviours of this problem, an automatic time step adjustment scheme is used to track the process accurately and effectively, by changing the time steps according
to the variations of solutions between adjacent time steps. The spatial resolution should be able to resolve the concentration gradients that drive the process. The diffusive boundary layer depends on the formation and fluid properties, where the gradients can be rather large, not only because the thickness of diffusive boundary layer is quite thin, but also because some regions near the top boundary layer are compressed due to the upwelling water with low CO2 concentration after the onset of convection. For the vertical discretization, the grid block sizes near the top boundary
are chosen such that they are sufficiently smaller than the thickness of the diffusive boundary layer obtained by linear
stability analysis, while a larger grid resolution is used in other regions. The penetration depth at the onset of instability is given as dc ¼ 34l/D=ðKDqgÞ ¼ 12 mm [3]. A vertical grid is set to be 1 mm near the top boundary, which is able to resolve the gradients of dissolved CO2 concentrations. Away from the top boundary, vertical grid spacing was stretched from 1 mm to 1 cm. Horizontal grid resolution is set up to be smaller than the size of the critical wavelength. The critical wavelength is given as kc ¼ 96:23 l/D=ðKDqgÞ ¼ 34 mm [14]. A horizontal grid resolution of 1 cm is used to capture the convective transport and the resulting solubility trapping. In order to obtain a better understanding on how the grid resolution affects the numerical solutions, grid dependence is examined for the 2D simulations with the physical domain size of 1m 1 m. According to the discretization method discussed above, the model is partitioned into 100  134 points. The domain

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

ความหนาแน่นของน้ำเกลือบริสุทธิ์บางส่วนและละลาย CO2 โดยรวมเศษส่วนเป็นทั้งฟังก์ชันของอุณหภูมิ ความดัน และมวลเศษเกลือ ตาม Battistelli et al. [38] และ Spycher และ Pruess [5], ตามลำดับ ความหนาแน่นบางส่วนของ CO2 ละลายใน Eq. (6) คำนวณได้จาก:ไตรมาสที่ 2 ¼MCO2วูð7Þที่วูเป็นฟังก์ชันของอุณหภูมิและเป็นอิสระจากแรงกดดัน โดย García [39] ในงานนี้ T เป็นค่าคง ตามวู = 34.92 ใช้ 10 6 m3/โมล MCO2 มีค่าเท่ากับ 44.0 10 กก. 3 โมล เมื่อแทนค่าสองค่าเป็น Eq. (7), q2 = 1260 กิโลกรัม/m3 จะรับ ในข้างต้น Eq. (6) รายการที่เพิ่มความหนาแน่นของส่วนผสมอควีเมื่อ CO2 ละลายลงในน้ำเกลือที่สภาพอ่างเก็บน้ำ Eqs (2)–(7) เป็นสมการควบคุมพื้นฐานสำหรับจำลองตัวเลขของการไหลในสื่อสำหรับละลายดัก porous ซึ่งเป็นกลุ่มของสมการไม่เชิงเส้นที่ควบคู่ และสามารถแก้ไขได้ โดยการวนซ้ำนิวตัน – Raphson2.2. วิธีตัวเลขEqs (2) และ (3) มี discretized เวลาใช้โครงร่างแตกต่างแน่นอนนัย และ ในพื้นที่โดยใช้วิธีผลต่างจำกัดเป็นการ วิธีการ discretization แก้ปัญหาชั่วคราว และพื้นที่ที่ใช้ในแบบจำลองมีอธิบายไว้ดังนี้ สำหรับ discretization เงื่อนไขการสะสมจะ discretized โดยการแนะนำค่าเฉลี่ยของระดับเสียงที่เหมาะสม:@ @tZวีเอ็น/tþDt /qXidVn ¼n qtþDtn XtþDtฉัน n /tnqtnXti; n วีเอ็นDt ð8Þที่ t + Dt แทนขั้นตอนเวลาใหม่ Fluxes เงื่อนไขด้วยการพา และ diffusive จะเลียนแบบเป็นผลรวมที่แยกกันของค่าเฉลี่ยมากกว่าพื้นผิวทั้งหมดของตาข่าย ZSnðqXiqÞ ndSn ¼XmqtþDtnm XtþDtฉัน nm qtþDtnm Snm ð9ÞZSnð/DqrXiÞ ndSn ¼Xm/ DqtþDtnm rXtþDtฉัน nm Snm ð10Þซึ่งเงื่อนไขเหล่านี้ไหลถือว่าทั้งหมดเป็นนัย กำหนดค่าในขั้นตอนเวลาใหม่ กฎหมายของ Darcy เป็น discretized ใน [40], qtþDtnm ¼Kl nmPtþDtn PtþDtmDnmþ qtþDtnm g cos เป็น!ð11Þซึ่งการจะแยกมุมระหว่างอัตราเร่งของความโน้มถ่วงและบรรทัดเซ็กเมนต์จากตาข่าย m กับ n มีทิศทางการหมุนจาก g บรรทัดเซ็กเมนต์ตามเข็มนาฬิกาตามที่ระบุใน Fig. 2 ใน Eq. (11), cosa เขียนได้ดังนี้:cos a ¼Z2 Z1Dnmð12Þที่ Z1, Z2 เป็น elevations บล็อกตารางที่เกี่ยวข้อง แล้วมีเขียนกฎหมายของ Darcy ระหว่างตาข่าย m และ n เป็นqtþDtnm ¼Kl nmPtþDtn PtþDtmDnmþ qtþDtnm gZ2 Z1Dnm!ð13Þที่ qnm, Xi, nm, nm (K/l) มีค่าเฉลี่ยมีค่าระหว่างตาข่าย m และ n กำหนด โดย qnm ¼D1qm þ D2qnDnmð14Þซีอานซีกวน nm ¼D1Xi; m þ D2Xi; nDnmð15ÞKl nm¼D1ðK = lÞm þ D2ðK = lÞnDnmð16Þที่ Dnm ได้เป็น:Þง 1 Dnm ¼ D2 ð17Þ2.3. การพิจารณาการแก้ปัญหาวิวัฒนาการและการโต้ตอบของไหลและการขนส่งช่วงกว้างของระดับพื้นที่และเวลาคือ คุณลักษณะสำคัญของดักละลาย ดังนั้น วิธีแก้ไขชั่วคราว และพื้นที่มีให้เพียงพอเพื่อแก้ไขลักษณะการทำงานของขั้นตอนที่สำคัญ สั่งในอากัปกิริยาการ diffusive ขั้นต้น และได้รับการพาเวลาเริ่มมีอาการได้อย่างถูกต้อง ขั้นเล็กเวลา 0.01 s จะใช้ ในเวลาเดียวกัน เนื่องจากพฤติกรรมระยะยาวของปัญหา แผนงานปรับปรุงขั้นตอนการเวลาอัตโนมัติจะใช้ในการติดตามกระบวนการได้อย่างถูกต้อง และมี ประสิทธิภาพ โดยการเปลี่ยนขั้นตอนเวลาตามกับรูปแบบของโซลูชั่นระหว่างขั้นตอนเวลาติด การแก้ปัญหาพื้นที่ควรผิดไดรฟ์นั้นการไล่ระดับสีเข้มข้น ชั้นขอบ diffusive ขึ้นอยู่กับการก่อตัวและน้ำมันคุณสมบัติ ที่การไล่ระดับสีจะค่อนข้างใหญ่ ไม่เท่านั้นเนื่องจากความหนาของชั้นขอบเขต diffusive ค่อนข้างบาง แต่เนื่องจากมีรวมบางภูมิภาคใกล้ชั้นขอบด้านบนเนื่องจากน้ำความเข้มข้น CO2 ต่ำหลังจากเริ่มมีอาการของการพา upwelling สำหรับ discretization แนวตั้ง เส้นบล็อกขนาดใกล้ขอบด้านบนจะเลือกที่จะพอมีขนาดเล็กกว่าความหนาของชั้นขอบเขต diffusive ได้ โดยเส้นวิเคราะห์เสถียรภาพ ในขณะที่ใช้ความละเอียดเส้นใหญ่ในภูมิภาคอื่น ๆ เจาะความลึกที่เริ่มมีอาการของความไม่แน่นอนให้เป็น dc ¼ 34l/D = ðKDqgÞ ¼ 12 มม. [3] เส้นตารางแนวตั้งถูกตั้งค่าให้เป็น 1 มม.ใกล้ขอบด้านบน สามารถแก้ไขไล่ระดับสีของความเข้มข้นของ CO2 ที่ละลาย ได้ยืดเส้นตารางแนวตั้งระยะห่างจากขอบด้านบน จาก 1 มม.ถึง 1 ซม.แนวนอน ความละเอียดตารางค่าจะเล็กกว่าขนาดของความยาวคลื่นสำคัญ ความยาวคลื่นสำคัญได้เป็น kc ¼ 96:23 l/D = ðKDqgÞ ¼มม. 34 [14] เส้นตารางแนวนอนความละเอียดของ 1 ซม.จะใช้การขนส่งด้วยการพาและดักละลายได้ เพื่อให้ได้เข้าใจในโซลูชั่นเลขกระทบความละเอียดตาราง ตารางพึ่งพาเป็นการตรวจสอบสำหรับแบบจำลอง 2D กับโดเมนที่มีอยู่จริงขนาด 1 เมตร 1 เมตร ตามวิธี discretization ที่กล่าวถึงข้างต้น รูปแบบแบ่งออกเป็น 100 คะแนน 134 โดเมน

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

ที่ความหนาแน่นของบางส่วนของน้ำเกลือบริสุทธิ์และส่วนมวล CO2 ละลายมีฟังก์ชั่นทั้งสองอุณหภูมิความดันและมวลส่วนเกลือคำนวณโดย Battistelli และคณะ [38] และ Spycher และ PRUESS [5] ตามลำดับ ความหนาแน่นของบางส่วนของ CO2 ละลายในสมการ (6) จะถูกคำนวณจาก:
q2 ¼
MCO2
Vu
ð7Þ
ที่ Vu เป็นหน้าที่ของอุณหภูมิและเป็นอิสระจากความดันที่กำหนดโดยGarcía [39] ในงานนี้ T เป็นค่าคงที่ตาม Vu = 34.92? 10? 6 m3 / mol ถูกนำมาใช้ MCO2 เท่ากับ 44.0? 10 3 กก. / mol เมื่อทั้งสองค่าถูกเปลี่ยนตัวลงไปในสมการ (7), q2 = 1260 kg / m3 จะได้รับ ในข้างต้นสม (6) แสดงให้เห็นว่ามีความหนาแน่นของสารผสมน้ำเพิ่มขึ้นเมื่อ CO2 ละลายลงไปในน้ำเกลือที่มีเงื่อนไขอ่างเก็บน้ำ สม (2) - (7) เป็นสมการปกครองพื้นฐานสำหรับการ
จำลองเชิงตัวเลขของการไหลในสื่อที่มีรูพรุนสำหรับดักละลายซึ่งเป็นชุดของสมการเชิงเส้นคู่และสามารถแก้ไขได้โดยการย้ำนิวตันราฟสัน.
2.2 วิธีการเชิงตัวเลข
สมการ (2) และ (3) มีการ discretized ในเวลาโดยใช้ความแตกต่างแน่นอนนัยโครงการและอยู่ในพื้นที่โดยใช้วิธีการที่แตกต่างกันแน่นอนหนึ่ง วิธีเนื่อง, ความละเอียดพื้นที่และเวลาที่ใช้ในการจำลองมีรายละเอียดดังนี้ สำหรับ discretization เงื่อนไขการสะสมกำลัง discretized โดยการนำค่าเฉลี่ยปริมาณที่เหมาะสม:
@ @t
Z
Vn
/ qXidVn ¼ / tþDt
n qtþDt
n XtþDt
ฉัน; n? / T
NQT
nXti
; n
? ?
?
Vn
Dt ð8Þ
ที่ t + Dt หมายถึงขั้นตอนเวลาใหม่ ฟลักซ์ของคำไหลเวียนและ diffusive ห้วงจะเป็นผลรวมที่ไม่ต่อเนื่องของค่าเฉลี่ยกว่าทุกพื้นผิวของตาข่าย, Z
Sn
ðqXiqÞ? NDSN ¼
X
เมตร
qtþDt
นาโนเมตรXtþDt
ฉัน; นาโนเมตรqtþDt
นาโนเมตรสนมð9Þ
Z
Sn
D / DqrXiÞ? NDSN ¼
X
เมตร
/ DqtþDt
นาโนเมตรrXtþDt
ฉัน; นาโนเมตรสนมð10Þ
ที่แง่ฟลักซ์เหล่านี้จะถือว่าเป็นนัยอย่างเต็มที่กำหนดโดยค่าที่ขั้นตอนเวลาใหม่ กฎหมายดาร์ซีเป็น discretized ในทางต่อไปนี้ [40], qtþDt
นาโนเมตร¼?
K
L
? ?
นาโนเมตร
PtþDt
n? PtþDt
เมตร
DNM
þqtþDt
นาโนกรัม cos
!
ð11Þ
ที่เป็นมุมสี่แยกระหว่างเร่งแรงโน้มถ่วงและส่วนของเส้นจากตาข่ายเมตรถึง n กับทิศทางการหมุนจาก g เพื่อส่วนของเส้นเข็มนาฬิกาตามที่ระบุไว้ในรูป 2. ในสมการ (11), โคซาสามารถเขียนได้ดังนี้
cos ¼
Z2? Z1
DNM
ð12Þ
ที่ Z1, Z2 เป็นเอนไซม์บล็อกตารางตามลำดับ กฎหมายดาร์ซีระหว่างตาข่ายเมตรและ n เป็นลายลักษณ์อักษรแล้วเป็น
qtþDt
นาโนเมตร¼?
K
L
? ?
นาโนเมตร
PtþDt
n? PtþDt
เมตร
DNM
þqtþDt
นาโนกรัม
Z2? Z1
DNM
!
ð13Þ
ที่ Qnm, Xi, นาโนเมตร (K / ลิตร) นาโนเมตรเป็นค่าเฉลี่ยประสานระหว่างตาข่ายเมตรและ n ที่ได้รับโดย Qnm ¼
D1qm þ D2qn
DNM
ð14Þ
Xi; นาโนเมตร¼
D1Xi; M þ D2Xi; n
DNM
ð15Þ
K
L
? ?
นาโนเมตร
¼
D1ðK = lÞmþD2ðK = lÞn
DNM
ð16Þ
DNM ที่จะได้รับเป็น:
DNM ¼ D1 D2 þð17Þ
2.3 พิจารณามติ
วิวัฒนาการและการมีปฏิสัมพันธ์ของการไหลและการขนส่งในช่วงกว้างของเครื่องชั่งพื้นที่และเวลาเป็นคุณสมบัติที่สำคัญของการวางกับดักสามารถในการละลาย ดังนั้นมติเวลาและพื้นที่จะต้องมีเพียงพอที่จะแก้ไขคุณสมบัติการไหลที่จำเป็น เพื่อที่จะจับภาพช่วงเริ่มต้นพฤติกรรม diffusive และได้รับเวลาที่เริ่มมีอาการของการพาความร้อนถูกต้องขั้นตอนเวลาเล็ก ๆ ของ 0.01 จะใช้ ในเวลาเดียวกันเนื่องจากพฤติกรรมในระยะยาวของปัญหานี้ขั้นตอนเวลาอัตโนมัติรูปแบบการปรับใช้ในการติดตามกระบวนการอย่างถูกต้องและมีประสิทธิภาพโดยการเปลี่ยนขั้นตอนเวลาตาม
การเปลี่ยนแปลงของการแก้ปัญหาระหว่างขั้นตอนเวลาที่อยู่ติดกัน ความละเอียดเชิงพื้นที่ควรจะสามารถที่จะแก้ปัญหาการไล่ระดับสีเข้มข้นที่ผลักดันกระบวนการ บริเวณชั้น diffusive ขึ้นอยู่กับการก่อตัวและคุณสมบัติของของเหลวที่ไล่ระดับสีจะค่อนข้างใหญ่ไม่ใช่เพียงเพราะความหนาของชั้นขอบเขต diffusive ผอมมาก แต่ยังเป็นเพราะพื้นที่บางส่วนที่อยู่ใกล้บริเวณชั้นบนสุดจะถูกบีบอัดเนื่องจากน้ำเต็มตื่น ที่มีความเข้มข้นของ CO2 ต่ำหลังจากที่เริ่มมีอาการของการพาความร้อน สำหรับเนื่องแนวตั้งขนาดของบล็อกตารางใกล้พรมแดนด้านบน
ได้รับการแต่งตั้งดังกล่าวว่าพวกเขาจะพอมีขนาดเล็กกว่าความหนาของชั้นขอบเขต diffusive ได้โดยการเชิงเส้น
การวิเคราะห์เสถียรภาพในขณะที่ความละเอียดตารางขนาดใหญ่ถูกนำมาใช้ในภูมิภาคอื่น ๆ เจาะลึกที่เริ่มมีอาการของความไม่แน่นอนจะได้รับเป็น dc ¼ 34l / D = ðKDqgÞ¼ 12 มม [3] ตารางแนวตั้งถูกตั้งค่าให้เป็น 1 มิลลิเมตรซึ่งอยู่ใกล้กับพรมแดนด้านบนซึ่งจะทำให้สามารถที่จะแก้ปัญหาการไล่ระดับสีของความเข้มข้นของ CO2 ที่ละลายในน้ำ ห่างจากพรมแดนด้านบนระยะห่างตารางแนวนอนเหยียดยาวตั้งแต่วันที่ 1 มมถึง 1 ซม. ความละเอียดตารางแนวนอนจะตั้งขึ้นจะมีขนาดเล็กกว่าขนาดของความยาวคลื่นที่สำคัญ ความยาวคลื่นที่สำคัญจะได้รับเป็น kc ¼ 96:23 ลิตร / D = ðKDqgÞ¼ 34 มิลลิเมตร [14] ความละเอียดแนวนอนของตาราง 1 ซม. ที่ใช้ในการจับภาพการขนส่งไหลเวียนและการละลายดักผล เพื่อที่จะได้รับความเข้าใจที่ดีขึ้นเกี่ยวกับวิธีการที่มีความละเอียดตารางมีผลกระทบต่อการแก้ปัญหาเชิงตัวเลขการพึ่งพาตารางการตรวจสอบสำหรับการจำลองแบบ 2 มิติที่มีขนาดทางกายภาพของโดเมน 1m? 1 เมตร ตามวิธีการ discretization ที่กล่าวข้างต้นรุ่นที่มีการแบ่งออกเป็น 100? 134 จุด โดเมน

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

ที่ความหนาแน่นของน้ำบริสุทธิ์บางส่วนและละลาย CO2 เศษส่วนมวลเป็นฟังก์ชันของอุณหภูมิ ความดัน และเศษส่วนมวลเกลือ คำนวณโดย battistelli et al . [ 38 ] และ spycher และ พรูส [ 5 ] ตามลำดับ ส่วนความหนาแน่นของปริมาณ CO2 ในอีคิว ( 6 ) คำนวณจาก :
2


¼ mco2 วูð 7 Þ
ที่วูเป็นฟังก์ชันของอุณหภูมิและอิสระของความดันให้โดยกาโอ การ์ซีอา [ 39 ]ในงานนี้ ไม่ได้เป็นค่าคงที่ดังนั้น Vu = 34.92  10  6 m3 / โมลที่ใช้คือ mco2 เท่ากับ 44.0  10  3 กิโลกรัมต่อโมล เมื่อทั้งสองค่าทดแทนในอีคิว ( 7 ) , Q2 = 1260 kg / m3 จะได้รับ ในข้างต้น อีคิว ( 6 ) พบว่า ความหนาแน่นของสารละลายผสมเพิ่มขึ้นเมื่อ CO2 ละลายในน้ำเกลือที่สภาพอ่างเก็บน้ำ EQS . ( 2 ) และ ( 7 ) ถือเป็นพื้นฐานสมการควบคุมสำหรับ
การจำลองเชิงตัวเลขสำหรับการไหลในวัสดุพรุนสำหรับดักการละลายซึ่งเป็นชุดของสมการไม่เชิงเส้นคู่ และสามารถแก้ไขได้โดยนิวตัน ราฟสันและทำซ้ำ .
2.2 .
EQS วิธีเชิงตัวเลข ( 2 ) และ ( 3 ) แบบจุดในเวลา ใช้ผลต่างสืบเนื่องโดยปริยาย และในพื้นที่โดยใช้วิธีผลต่างอันตะปริพันธ์ . ส่วนค่าโดยพื้นที่และเวลาที่ใช้ในการจำลองมติอธิบายไว้ดังนี้ สำหรับค่า , การสะสมเงื่อนไขแบบจุดโดยการนำค่าเฉลี่ยปริมาณที่เหมาะสม :
@ @ T
z

/ qxidvn VN ¼ / T þ DT
n Qt þ DT
n XT þ DT
i ; N  / T
nqt nxti

; N



   จากแฟรช ð 8 Þ
ที่ t DT เป็นก้าวครั้งใหม่ต่อการพาและการแพร่ของเทอมโดยประมาณเป็นแบบผลรวมของค่าเฉลี่ยกว่าพื้นผิวทั้งหมดของตาข่าย , Z
SN
ð qxiq Þ  ndsn ¼
x
m
Qt þ DT
nm XT þ DT
i ; nm Qt þ DT
nm SNM ð 9 Þ
z

ð / สินี dqrxi Þ  ndsn ¼
x
m
/ dqt þ DT
nm rxt þ DT
i ; nm SNM ð 10 Þ
ที่เงื่อนไขเหล่านี้จะถือว่าเป็นค่าเต็มที่เปิดเผย ให้คุณค่า ในขั้นตอนเวลาใหม่ กฎของดาร์ซีเป็นแบบจุดในวิธีดังต่อไปนี้ [ 40 ]Qt þ DT

L
k  nm ¼
 

þ DT
n nm PT PT  þ DT
m

þ dnm Qt þ DT
nm G เพราะเป็น
! Þ

ð 11 ซึ่งเป็นจุดตัดระหว่างความเร่งโน้มถ่วงมุมและส่วนของเส้นตรงจากตาข่าย M N กับทิศทางการหมุนตามเข็มนาฬิกาจาก G ไปเส้นตรงตามที่ระบุไว้ในรูปที่ 2 ในอีคิว ( 11 ) , โคซาสามารถเขียนได้ดังนี้

เพราะเป็น¼กขึ้น  Z1


ที่ 12 dnm ðÞ Z1 กขึ้นตามตารางบล็อก , มีความสูงกฎของดาร์ซีระหว่างตาข่าย m และ n เป็นแล้วเขียน Qt þ DT


L
k  nm ¼
 

þ DT
n nm PT PT  þ DT
m

þ dnm Qt þ DT
nm g
กขึ้น  dnm Z1

! !

ที่ 13 Þð qnm Xi , NM ( K / L ) nm เป็นค่าเฉลี่ยฮาร์มอนิก ระหว่าง ตาข่าย และ ให้ โดย qnm ¼
d1qm þ d2qn

ð dnm 14 Þ
Xi ; nm ¼
d1xi ; M þ d2xi ; N

dnm 15 ðÞ
K
L
 
nm
¼
D1 ð K = L Þ M þ D2 ð K = L Þ n

ð dnm 16 Þ
ที่ dnm ให้ :
dnm ¼þ D1 D2 ð 17 Þ
2.3การพิจารณาของมติ
วิวัฒนาการและปฏิสัมพันธ์ของการไหลและการขนส่งในช่วงกว้างของพื้นที่และเวลาเป็นเครื่องลักษณะสําคัญของการละลายและการ ดังนั้นเวลาพื้นที่และมติต้องเพียงพอที่จะแก้ไขคุณสมบัติการไหลที่สําคัญ เพื่อที่จะจับช่วงแรกกระจายพฤติกรรมและได้รับการโจมตีเวลาของการพาได้อย่างถูกต้องขั้นตอนขนาดเล็กที่เวลา 0.01 S ใช้ ทั้งนี้ พฤติกรรมในระยะยาวปัญหานี้โดยอัตโนมัติขั้นตอนการปรับเวลาแบบใช้ในการติดตามกระบวนการที่ถูกต้อง และมีประสิทธิภาพ โดยการเปลี่ยนเวลา ขั้นตอนตาม
กับรูปแบบของโซลูชั่นระหว่างขั้นตอนเวลาที่อยู่ติดกันความละเอียดเชิงพื้นที่สามารถแก้ไขอัตราส่วนของความเข้มข้นที่ทำให้กระบวนการ ชั้นขอบเขตการแพร่ขึ้นอยู่กับรูปแบบและคุณสมบัติของของไหลที่มีการไล่ระดับสีได้ค่อนข้างมาก เพราะไม่เพียงความหนาของชั้นขอบเขตการแพร่ค่อนข้างบางแต่ยังเพราะบางพื้นที่ใกล้ boundary layer ด้านบนจะถูกบีบอัดเนื่องจากการกระบวนการน้ำผุดน้ำกับ CO2 ความเข้มข้นต่ำหลังจากการโจมตีของการพาความร้อน สำหรับค่าแนวตั้งตารางบล็อกขนาดใกล้
ขอบด้านบนเลือกเช่นที่พวกเขามีเพียงพอมีขนาดเล็กกว่าความหนาของชั้นขอบเขตการแพร่ที่ได้จากการวิเคราะห์เสถียรภาพเชิงเส้น
,ในขณะที่ความละเอียดตารางขนาดใหญ่ใช้ในภูมิภาคอื่น ๆ ความลึกของการสอดใส่ที่เริ่มขาดเสถียรภาพให้ DC ¼ 34l / D = ð kdqg Þ¼ 12 มม. [ 3 ] ตารางตามแนวตั้งถูกตั้งค่าเป็น 1 มม. ตรงบริเวณด้านบน ซึ่งสามารถแก้ปัญหาการไล่สีของปริมาณคาร์บอนไดออกไซด์เข้มข้น ห่างจากขอบด้านบนตารางแนวตั้งถูกยืดระยะตั้งแต่ 1 มม. ถึง 1 เซนติเมตรความละเอียดแนวนอนมีการตั้งค่าตารางมีขนาดเล็กกว่าขนาดของความยาวคลื่นที่สําคัญ คลื่นวิกฤตให้ KC ¼ 96:23 L / D = ð kdqg Þ¼ 34 มม. [ 14 ] ความละเอียดแนวนอนของตาราง 1 ซม. ใช้ในการจับและการขนส่งโดยการเกิดการละลายและการ เพื่อที่จะได้รับความเข้าใจที่ดีขึ้นเกี่ยวกับวิธีการตารางความละเอียดต่อการแก้ปัญหาเชิงตัวเลขตารางการตรวจสอบสำหรับ 2D จำลองกับทางเมนขนาด 1m  1 เมตร ตามวิธีที่กล่าวข้างต้นค่า นางแบบก็แบ่งเป็น 100  134 คะแนน โดเมน