ssmodule.vb

一个lbm的程序

Module ssModule
    'field description
    Public field(,) As BasicCell
    Public rho(,), u(,), v(,), eps(,), mass(,) As Double
    Public flag(,) As Byte
    'used for calculation
    Public ei(,) As Int16 = {{0, 0}, {1, 0}, {0, 1}, {-1, 0}, {0, -1}, {1, 1}, {-1, 1}, {-1, -1}, {1, -1}}
    Public eiu(,) As Int16 = {{0, 0}, {1, 0}, {0, 1}, {-1, 0}, {0, -1}, {1 / Math.Sqrt(2), 1 / Math.Sqrt(2)}, {-1 / Math.Sqrt(2), 1 / Math.Sqrt(2)}, {-1 / Math.Sqrt(2), -1 / Math.Sqrt(2)}, {1 / Math.Sqrt(2), -1 / Math.Sqrt(2)}}
    Public wi() As Double = {4.0 / 9, 1.0 / 9, 1.0 / 9, 1.0 / 9, 1.0 / 9, 1.0 / 36, 1.0 / 36, 1.0 / 36, 1.0 / 36}
    ''Public sita() As double = {0, 0, Math.PI / 2, Math.PI, 3 * Math.PI / 2, Math.PI / 4, 3 * Math.PI / 4, 5 * Math.PI / 4, 7 * Math.PI / 4}
    Public s() As Double = {0, 0.5, 0.5, 0.5, 0.5, 0.25, 0.25, 0.25, 0.25}
    Public oppk() As Double = {0, 3, 4, 1, 2, 7, 8, 5, 6}
    Public gw() As Double = {0, 1, 1 / Math.Sqrt(2), 1, 1 / Math.Sqrt(2), 1, 1 / Math.Sqrt(2), 1, 1 / Math.Sqrt(2)}
    'array for interface cells
    Public itfArray As New List(Of Int32())
    'parameters
    Public w, r0, ra, beta, g As Double
    Public m, n, steps, interval As Int32
    'const variables
    Public Const kapa = 0.001
    Public Const q = 8, CONST_FLUID = 9, CONST_EMPTY = 0, CONST_WALL = 1, CONST_NONE = 2, CONST_IOLET = 3
    Public Const CONST_TMPS = 5, CONST_MEDIUM = 6, CONST_SAW = 7, CONST_INTERFACE = 8
    Public Const strNextLine = Chr(13) + Chr(10)
    'needed for GUI
    Public strWrite As String
    Public version As String
    Public DataDescription As String
    Public showon As Boolean
    Public candraw As Boolean
    Public stepnum As Int32
End Module