|
| float | extent = 30.0 |
| float | gridDelta = 0.5 |
| tuple | bounds = (-extent, extent, -extent, extent) |
| tuple | boundaryCons |
| | substrate = vls.Domain(bounds, boundaryCons, gridDelta) |
| tuple | origin = (0.0, 0.0) |
| tuple | planeNormal = (0.0, 1.0) |
| | mesh = vls.Mesh() |
| | minNodeDistFactor |
| | trench = vls.Domain(bounds, boundaryCons, gridDelta) |
| tuple | minCorner = (-extent / 6.0, -25.0) |
| tuple | maxCorner = (extent / 6.0, 1.0) |
| | surfaceMeshBox = vls.ToSurfaceMesh(trench, mesh, minNodeDistFactor=0.02) |
| | surfaceMeshTrench = vls.ToSurfaceMesh(substrate, mesh, minNodeDistFactor=0.02) |
| | substrateFE = vls.Domain(substrate) |
| | newLayerFE = vls.Domain(substrateFE) |
| | substrateRK2 = vls.Domain(substrate) |
| | newLayerRK2 = vls.Domain(substrateRK2) |
| | substrateRK = vls.Domain(substrate) |
| | newLayerRK = vls.Domain(substrateRK) |
| | velocities = velocityField() |
| float | totalSimulationTime = 16.5 |
| float | outputInterval = 0.5 |
| | advectionKernelFE = vls.Advect() |
| | passedTimeFE |
| | advectionKernelRK2 = vls.Advect() |
| | passedTimeRK2 |
| | advectionKernelRK = vls.Advect() |
| | passedTimeRK |
| | writer = vls.WriteVisualizationMesh() |
| | multiMeshKernel = vls.ToMultiSurfaceMesh(minNodeDistFactor=0.02) |
| | multiMesh = vls.Mesh() |