Giải pháp điều kiện biên Dirichlet-Neumann trở nên không ổn định - Phương pháp hiệu chỉnh áp suất


12

Tôi đang mô phỏng dòng chảy không thể nén trên một hình trụ ở số Reynold là 500. Tôi đang giải phương trình stokes điều hướng bằng phương pháp hiệu chỉnh áp suất. Giải pháp của tôi trở nên không ổn định sau một thời gian nhất định (khoảng 5s).

Tôi đã thử tinh chỉnh lưới của mình, stepize (0,05) (đảm bảo CFL của tôi <1, mặc dù tôi đang sử dụng các phương thức ngầm)

Điều kiện biên, lưới và kết quả không ổn định của tôi được hiển thị trong các hình đính kèm. Tên miền lớn hơn khoảng 25 lần so với đường kính xi lanh.

Tôi đã thử mô phỏng vấn đề này lưới O (trở nên không ổn định gần như ngay lập tức).

Các liên kết sau đây chứa các hình ảnh của các điều kiện biên và kết quả.

Điều kiện biên

Sự bất ổn

Tôi sẽ biết ơn nếu bất cứ ai có thể chia sẻ suy nghĩ / kinh nghiệm của họ về vấn đề này. Cảm ơn nhiều.

đã chỉnh sửa:

Xin lỗi vì lỗi đánh máy:

Tôi đang sử dụng các điều kiện biên sau: Ranh giới Neumann

bạnn-np= =0;

trên ranh giới Dirichlet

bạn= =bạnx= =1

đã chỉnh sửa:

tôi đã áp dụng các điều kiện biên vận tốc trên các nút xung quanh ranh giới dirichlet. Ngoài ra, nút góc trên bên phải và dưới cùng bên phải là ranh giới dirichlet với vận tốc 1.

Sau đó, tôi nhìn sâu hơn vào các kết quả mô phỏng, tôi nhận thấy rằng sự không ổn định bắt đầu len vào ở ngã ba dòng chảy / dòng chảy.


2
Làm thế nào, cụ thể, bạn đang thực hiện các điều kiện biên của bạn? Điều này có thể làm cho tất cả sự khác biệt trong một mô phỏng như thế này.
Kyle Mandli

0-np= =0nbạn= =x(bạnx,0,0)= =0

2
Phương pháp mà bạn sử dụng là gì? FEM? Với sự ổn định? Bạn đã cố gắng để giảm số Reynold?
Dr_Sam

Câu trả lời:


1

Tôi đã tìm ra vấn đề. Tôi đã phải tăng kích thước của miền hơn nữa để loại bỏ các hiệu ứng biên. Hơn nữa, tôi đã phải giảm số CFL xuống khoảng 0,5-1,0

Tôi nghĩ rằng số CFL cần phải giảm hơn nữa cho số reynold cao hơn.

Ban đầu, tôi nghĩ rằng tôi đã giảm kích thước bước đủ, nhưng thực tế không phải vậy.


2
bạnnnbạn

3
Thay vì trả lời câu hỏi của riêng bạn, bạn nên chỉnh sửa câu hỏi ban đầu để bao gồm thông tin bổ sung. Điều này giúp dễ dàng có tất cả thông tin trong một địa điểm và do đó để trả lời câu hỏi của bạn.
Christian Clason

2
Một nhận xét về suy nghĩ của bạn - số CFL có lẽ không cần phải giảm cho số Reynold cao hơn. Max Gunzberger trong cuốn sách FEM của mình cho Viscous Incomp Flows lưu ý rằng bán kính hội tụ của phương pháp Newton bị thu hẹp khi tăng số Reynold và giảm CFL hạn chế dấu thời gian, có thể được hiểu là tăng thời gian chính quy phép lặp Newton thuần túy.
Jesse Chan

Không phải một ranh giới Neumann cho vận tốc trên hai ranh giới ngang sẽ phù hợp hơn? Tôi đoán là khi bạn đang áp đặt Dirichlet, ranh giới vẫn không còn xa.
Discrete_Reynolds
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.