WORKSHOP:-
SEM prototype: A lighter fluid penetrating heavier liquid.

This workshop example shows how to set up and run the simulation of the lighter fluid (oil) stream penetrating, from the slot at the top wall, the domain filled with heavier (water) liquid.

The task is to simulate the velocity field and concentrations distributions.

The geometry and expected velocity vectors are shown in the figure below.

The inlet oil velocity is 2.0 m. The water density is 1000.0 and the oil density is 700, , kg/m3. Two-dimensional domain has 1 x 1 meters in Y and Z directions.

Accessing PHOENICS-VR

From the system level:

To enter the PHOENICS-VR environment, click on the PHOENICS icon on the desktop, or click on Start, programs, PHOENICS, PHOENICS.

From the commander level:

To enter the PHOENICS-VR environment, click on the 'Run vre' icon in the left column.

In PHOENICS-VR environment,

Start with an 'empty' case - click on 'File' then on 'Start New Case', then on 'Core', then click on 'OK'; to confirm the resetting.

To enter VR Editor:

This is the default mode of operation .

Within VR Editor.

Set the domain size:

Leave the default sizes.

Activate the solution settings for velocities and mass fraction of water:

Click on 'Menu' and set the 'Oil penetrating into water' as the Title.

Click on 'Models'

Leave 'Solution for velocity and pressure' to ON.

Click on 'Page down' and click 'Settings' of 'Extra variables'.

Set 'SOLVE' C1 (mass fraction of water)

and click on 'Apply'

Click on 'Previous panel'

Click on 'Properties'

Set 'Use property tables' OFF

Click on 'Density' CONSTANT, from the list of options select'Linear in scalar' and click 'OK'

Set A=700 and B=1000-700=300

Change 'Storage' ON

Change 'Viscosity' 1.e-6

Click on 'Sources'

Set 'Gravitational Forces' ON and click 'OK' to confirm.

Set 'Gravitational acceleration' -9.81 in Z direction.

Set 'Reference density' 1000.

Click 'Top menu' and click 'OK'

Adjust the view point:

Click on 'Move right' toggle to view in X-direction.

Click 'Reset' to line up.

Create the OIL-IN object:

Click on 'Settings', 'New' and 'New Object'.

Change name to OIL-IN.

Click on 'Size' and set SIZE of object as:

Xsize: 1.0

Ysize: 0.2

Zsize: 0.0

Click on 'Place' and set POSITION of object as:

Xpos: 0.0

Ypos: 0.4

Zpos: 1.0

Click on 'General'.

Define Type: Inlet.

Click on 'Attributes'

Set velocity in Z-direction to -2.0 m/s.

Set 'Inlet density' at 'Domain fluid'.

Set C1 'Inlet value' 0.0

Click on 'OK' to close the Attributes menu,

and on 'OK' in the Object Specification Dialogue Box.

Create the FREE-N object:

Click on 'Settings', 'New' and 'New Object'.

Change name to FREE-N.

Click on 'Size' and set SIZE of object as:

Xsize: 1.0

Ysize: 0.0

Zsize: 1.0

Click on 'General'.

Define Type: Outlet.

Click on 'Attributes' and set C1 'External value' USER-SET.

Set 'User-set' 1.0

Click on 'OK' twice to exit the Object Dialogue Box.

Create the FREE-S object:

Click on 'Settings', 'New' and 'New Object'.

Change name to FREE-S.

Click on 'Size' and set SIZE of object as:

Xsize: 1.0

Ysize: 0.0

Zsize: 1.0

Click on 'Place' and set POSITION of object as:

Xpos: 0.0

Ypos: 1.0

Zpos: 0.0

Click on 'General'.

Define Type: Outlet.

Click on 'Attributes' and set C1 'External value' USER-SET.

Set 'User-set' 1.0

Click on 'OK' twice to exit the Object Dialogue Box.

Create the FREE-LOW object:

Click on 'Settings', 'New' and 'New Object'.

Change name to FREE-LOW.

Click on 'Size' and set SIZE of object as:

Xsize: 1.0

Ysize: 1.0

Zsize: 0.0

Click on 'General'.

Define Type: Outlet.

Click on 'Attributes' and set C1 'External value' USER-SET.

Set 'User-set' 1.0

Click on 'OK' twice to exit the Object Dialogue Box.

 

Set the spatial grid:

Click on the 'Mesh toggle' button. The default mesh (shown as orange lines) will appear on the screen.

Click anywhere on the image, and the 'Gridmesh settings' dialog box will appear.

Change the Number of cells to 10 cells in Y and 10 cells in Z.

Click on 'OK'; to close the dialog box. The mesh will now be shown as blue and orange lines. The orange lines are 'regions', defined by the bounding boxes of the objects.

Click on 'Mesh toggle' again to turn off the mesh display.

Set the remaining solution-output-control parameters:

Click on'Menu'.

Click on 'Numerics'.

Set the number of iterations to 300.

Click on 'Relaxation control'.

Switch 'Automatic convergence control' to OFF

Set relaxations for V1, W1 and C1 0.1, in the fields provided for 'DTFALS'.

Click on 'Previous panel' to get up.

Click on 'Output'.

Change the monitoring point by setting

IYMON to 5,

IZMON to 5.

Click on 'Previous panel' to get up.

Click on 'Top menu' to return to the top menu.

Click on 'OK' to exit the Main Menu.

Running the Solver.

In the PHOENICS-VR environment, click on 'Run', 'Solver'(Earth), and click on 'Yes' to confirm the action.

 

Using the VR Viewer.

In the PHOENICS-VR environment, click on 'Run', 'Post processor',then GUI Post processor (VR Viewer) , and view the velocity vectors, contours of water concentration and mixture density.

The calculated water concentration contour is shown below:

Saving the results.

In the PHOENICS-VR environment, click on 'Save as a case', make a new folder called 'OILIN' (e.g.) and save as 'CASE1' (e.g.).