infistream

6.6.2.2 Volume dan Pusat Massa Elemen

Volume dan pusat massa elemen. Setelah menghitung normal, luas, dan pusat wajah, mungkin untuk mengevaluasi metrik sel. Seperti halnya dengan wajah, ide dasar untuk sel polihedral adalah memecahkannya menjadi jumlah elemen tetrahedra. Operasi yang menghitung volume dan pusat massa elemen kemudian didefinisikan dalam file “$FOAM_SRC/OpenFOAM/meshes/primitiveMesh/primitiveMeshCellCentresAndVols.C” menggunakan fungsi yang ditunjukkan dalam Listing 6.9.

Listing 6.9 Fungsi yang digunakan untuk menghitung volume dan pusat massa elemen

Fungsi ini menggunakan empat argumen di mana dua pertama mewakili, masing-masing, pusat dan luas wajah. Dua argumen yang tersisa mengembalikan objek yang berisi pusat dan volume sel. Seperti yang ditunjukkan dalam Listing 6.10, loop “for” didefinisikan dua kali. Redundansi ini terkait dengan alamat LDU yang digunakan dalam OpenFOAM®, yang akan diperkenalkan dan dijelaskan dalam Bab 7.

Listing 6.10 Menghitung pusat perkiraan sel

Untuk menghitung pusat massa dan volume sel, langkah pertama adalah evaluasi , pusat massa sel yang disimbolkan dalam Listing 6.10 sebagai cEst, menggunakan hubungan pertama dalam Persamaan (6.25). Setelah nilai-nilai  diperoleh, perhitungan volume piramida dan pusat massa sel dilakukan dengan menerapkan Persamaan (6.25) menggunakan kode yang ditampilkan dalam Listing 6.11.

Listing 6.11 Menghitung pusat sel dan volume

Dalam kode di atas, fCtrs sesuai dengan xCE sementara “cEst – fCtrs” sesuai dengan vektor jarak dGf yang ditampilkan dalam Gambar 6.21. Nilai akhir diperoleh, seperti yang ditunjukkan dalam Listing 6.12, dengan membagi pusat massa sel oleh volume sel dan kemudian membagi volume sel dengan 3.

Listing 6.12 Nilai akhir pusat massa sel dan volume

Struktur data mesh untuk uFVM dan OpenFOAM® akan dijelaskan secara rinci dalam Bab 7.