CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Modeling free-surface river connection with interIsoFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 17, 2020, 11:14
Question Modeling free-surface river connection with interIsoFoam
  #1
Member
 
Grivalszki Péter
Join Date: Mar 2019
Location: Budapest, Hungary
Posts: 32
Rep Power: 4
GrivalszkiP is on a distinguished road
Hi Foamers,

Im trying to simulate a free surface river connection using interIsoFoam.
  • Geometry: fictitious from the internet. Two inlet, one outlet. Quite complex.
  • Mesh: with snappyHexMesh, checkMesh shows that everything is OK.
  • Boundary: separated inlet for water and air phases. Values are also fictitious. variableHeightFlowRateInletVelocity have been tried, failure.
  • Turbulence model: k-epsilon.

Unfortunately, it crashes after ~0.01 s. Where is my mistake?

Boundaries:

U
Code:
    in1_w
    {
        type            	fixedValue;	
		value				uniform (0.5 0 0);

    }

    in1_a
    {
        type            	fixedValue;
		value				uniform (0.5 0 0);

    }

    in2_w
    {
        type            	fixedValue;
		value				uniform (0 -0.6 0);

    }

    in2_a
    {
        type            	fixedValue;
		value				uniform (0 -0.6 0);

    }

    out
    {
        type            zeroGradient;

    }
	
	wall
	{
        type            fixedValue;
		value			uniform (0 0 0);
	}

	RiversComplete
	{
        type            fixedValue;
		value			uniform (0 0 0);
	}

	atm
	{
        type            slip;
	}
p_rgh
Code:
    atm
    {
        type            totalPressure;
        p0              uniform 0;
    }

    ".*"
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }
alpha
Code:
    in1_w
    {
        type            fixedValue;
		value           uniform 1; 
    }
	
    in1_a
    {
        type            fixedValue;
		value           uniform 0; 
    }

    in2_w
    {
        type            fixedValue;
		value           uniform 1; 
    }
	
    in2_a
    {
        type            fixedValue;
		value           uniform 0; 
    }

    out
    {
        type            zeroGradient;
		value			uniform 0;
    }

    atm
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    wall
    {
        type            zeroGradient;
    }

    RiversComplete
    {
        type            zeroGradient;
    }
k

Code:
    in1_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in1_a
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_a
    {
        type            fixedValue;
        value           $internalField;
    }

    out
    {
        type            zeroGradient;
    }
	
	wall
    {
        type            kqRWallFunction;
        value           $internalField;
    }

    atm
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

	RiversComplete
    {
        type            kqRWallFunction;
        value           $internalField;
    }
epsilon

Code:
    in1_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in1_a
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_w
    {
        type            fixedValue;
        value           $internalField;
    }
	
    in2_a
    {
        type            fixedValue;
        value           $internalField;
    }
	
    out
    {
        type            zeroGradient;
    }

    wall
    {
        type            epsilonWallFunction;
        value           $internalField;
    }

    atm
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }

    RiversComplete
    {
        type            epsilonWallFunction;
        value           $internalField;
    }
nut
Code:
    in1_w
    {
        type            calculated;
		value			uniform 0;
    }
	
    in1_a
    {
        type            calculated;
		value			uniform 0;
    }
	
    in2_w
    {
        type            calculated;
		value			uniform 0;
    }
	
    in2_a
    {
        type            calculated;
		value			uniform 0;
    }

    wall
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
	
	RiversComplete
    {
        type            nutkWallFunction;
        value           uniform 0;
    }


    out
    {
        type            zeroGradient;
    }
	
    atm
    {
        type            calculated;
		value			uniform 0;
    }
GrivalszkiP is offline   Reply With Quote

Reply

Tags
channel water flow, free surface flow, interisofoam, multiphase flow, river

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SU2 7.0.7 Built on CentOS 7, parallel computation pyscript mpi exit error? EternalSeekerX SU2 3 October 9, 2020 18:28
Problem in free surface modeling of water-oil mixture scythe FLUENT 4 November 12, 2013 11:47
Which code to go for free surface wave modeling.. Ted Chu Main CFD Forum 7 December 29, 2008 15:42
Hydrofoil modeling near free surface in fluent Mojtaba Lajevardi FLUENT 2 June 7, 2008 04:59
Modeling Free Surface Flows Elliot Schwartz Main CFD Forum 5 August 25, 1998 21:03


All times are GMT -4. The time now is 00:28.