Use of grid stretching functions in DiveMESH
 

Hi! When I am using reef3d,I do not understand some commands of divemesh,such as,B 101,B 102,B 103.Could you reply it?

The commands you name are relvant if you are using a stretched grid fro your simulations.

For eg.
B 1 0.05
B 10 0.0 10.0 0.0 0.05 0.0 1.0
Will result in a two-dimensional domain that is 10 m long in the x-direction, 1 cell wide (2D) and 1 m high.
In this domain, the grid size is dx=dy=dz= 0.05 m in all directions, with a total of 4000 grids (see fig1)

On the other hand, if you want finer meshes in a part of the domain compared to the rest, you can apply a stretched mesh. For eg.

B 2 200 1 20
B 10 0.0 10.0 0.0 0.05 0.0 1.0

B 101 5
B 111 2.0
B 114 5.0

B 103 5
B 111 2.0
B 116 0.5

will create a mesh with 4000 cells but with a focus point for x-grids at x=5.0 m, and z=0.5 m. The stretching factor is 2.0. This will result in a mesh as seen in Fig 2.
The file DIVEMesh_Log/DIVEMesh-Non-Uniform-Mesh.txt will show you the max and min sizes of the grids and the location of the finest grid.
This example uses a simple sinh function based stretching. Other options can be explored in the B 101/2/3 commands to define different types of grid stretching.

Hope that helps.

PS. stretching factor of 3.0 is used purely for demonstration. In practical cases, this might be too high. You need check the aspect ratios of the grids in the DIVEMesh-Non-Uniform-Mesh.txt file and is recommended to keep the maximum cell ratios to around 1.05 or lower to ensure simulation stability.

Ok! Thanks for your reply.

Hii,

Is it necessary to use B 1 or B 2 in control.txt ??

I have used the following codes to generate the grid ---

############################
B 10 0.0 5.0 0.0 0.06 0.0 0.3 //domain

B 101 11 // cell based
B 127 0.0025 0.005 3.7725 0.545 1.03
B 102 11
B 128 0.0025 0.0025 0.03 0.01 1.0 // uniform grid
B 103 11
B 129 0.0025 0.005 0.15 0.14 1.03

#######################

In B 127/128/129 - we are giving min/max cell size, the refinement region and cell ratio, so this would work without B1 and B2, correct ?

Paul

Yes it is essential to define
B 2
in your input.

As you can see the inputs for your commands do not have an idea of what the maximum number of cells is. This is obtained from B 2.
Without this information, the stretched grids cannot be defined.

Thank you for the clarification. I will correct my code.

Hi! I have one question.How to control the maximum cell ratios?

Depending on what you want to control, you need to select the grid generation method.
If you have a very particular mesh in mind, then simply use B 101/2/3 10 and provide files x-spacing, y-spacing, z-spacing define the grids in these directions.
The contents should range from 0 to 1 in these files.

Unable to generate the required geometry and fluid domain
 

Hello Arun
I am trying to create a cylinder in z-direction at location x = 17.0 and z = 0.2 and radius 0.1. However, the entire fluid domain is removed down below the cylinder. I am attaching the file below.

B 2 1000 1 10
B 10 0.0 40.0 0.0 0.04 0.0 1.0

S 32 17.0 0.2 0.1

// S 61 11.0 17.0 0.0 0.04 0.0 0.3
// S 10 17.0 19.0 0.0 0.04 0.0 0.3
// S 61 19.0 22.0 0.0 0.04 0.3 0.0

// B 101 5
// B 111 10.0
// B 114 17.0

// B 103 5
// B 113 10.0
// B 116 0.3

B 101 5
B 111 10
B 114 17.0

B 103 5
B 113 10.0
B 116 0.2

M 10 4
M 20 2

Is there a more up-to-date DiveMesh user guide that has explanations for the meaning of B 101 and the other parameters?

Can you please send a link to that user guide if it exists?

Yes, there is an explanation of the B 101 functions in the User Guide. The best way to see what they are doing is to take a simple case (e.g. the second CFD tutorial case) and try them out a bit. In general, for refinement around objects, B 101 11 gives the best result.