Browse code

Обновил README описанием некоторых тестов. Проверен тест test-kt-split-2

Dmitrii Kashin authored on 26/05/2015 01:46:07
Showing 2 changed files
... ...
@@ -1,3 +1,5 @@
1
+Описание репозитория:
2
+
1 3
 В данном репозитории содержится программа на языке Common Lisp,
2 4
 реализующая функционал для численного моделирования многофазных
3 5
 течений.
... ...
@@ -11,3 +13,25 @@
11 11
 следовательно были исключены из рассмотрения компоненты скорости,
12 12
 связанные с этой осью.
13 13
 
14
+
15
+Описание тестов:
16
+
17
+(test-kt-constant)
18
+Тест на постоянном поле.
19
+
20
+(test-kt-split-1)
21
+Тест на скачке при phi~1
22
+Остановка при t=0.005
23
+
24
+(test-kt-split-1/convergence)
25
+Тест на скачке при phi~1
26
+Остановка при t=0.005
27
+Остановка при невязке.
28
+
29
+(test-kt-split-1v2)
30
+Тест на скачке при phi~1
31
+Скорости искусственно заданы побольше, чтобы наблюдать движение скачка.
32
+Не наблюдаю.
33
+
34
+(test-kt-split-2)
35
+Тест на скачке при phi=0.9
... ...
@@ -196,8 +196,8 @@
196 196
       (loop for i from (min-x-index grid) to (max-x-index grid) do
197 197
 	   (loop for j from (min-y-index grid) to (min-y-index grid) do
198 198
 		(let ((pv (conservative->primitive (cell grid (index i j)))))
199
-		  (incf (ug pv) 100)
200
-		  (incf (ul pv) 100)
199
+		  (incf (ug pv) 500)
200
+		  (incf (ul pv) 500)
201 201
 		  (setf (cell grid (index i j)) (primitive->conservative pv)))))
202 202
       (format t "~%INCREMENT DEBUG:~%")
203 203
       (print-content t (cell grid (index 0 0)))
... ...
@@ -207,7 +207,7 @@
207 207
       (gridbg "4 (boundaries)" grid)
208 208
       (gridbg "start-grid" grid)
209 209
       (let ((newgrid (march grid
210
-			    :finish-step 1
210
+			    :finish-step 100
211 211
 			    :output-every-step 1
212 212
 			    :cu cu
213 213
 			    )))
... ...
@@ -216,7 +216,7 @@
216 216
 	)
217 217
       )))
218 218
 
219
-(defun test-kt-split-2 (&key (cells-nx 100) (cells-ny 1) (cu 1) (output-dir #P"~/prog/clfpv/test-kt-split-2/"))
219
+(defun test-kt-split-2 (&key (cells-nx 100) (cells-ny 1) (cu 0.001) (output-dir #P"~/prog/clfpv/test-kt-split-2/"))
220 220
   "Тест на стоячем скачке уплотнения при phi = 0.9. Этот тест нужен,
221 221
 чтобы проверить, как ведёт себя метод при расчёте не чисто газового
222 222
 течения."
... ...
@@ -268,8 +268,10 @@
268 268
       (gridbg "4 (boundaries)" grid)
269 269
       (gridbg "start-grid" grid)
270 270
       (let ((newgrid (march grid
271
-			    :finish-step 1000
272
-			    :output-every-step 50
271
+			    :finish-time 0.005
272
+			    :finish-convergence 1.0d-6
273
+			    :output-every-time 0.0001
274
+			    :output-every-step 100
273 275
 			    :cu cu
274 276
 			    )))
275 277
 	(gridbg "final-grid" newgrid)