[sourav8016]

Dear Sir,

I have been trying to install sedfoam. I am using openfoam6 in ubuntu 16.04. While running ./Allwmake command, I get the following errors

Foam::sedIncompressibleTurbulenceModel<TransportMo del>:Primef() const’:
../TurbulenceModels/turbulenceModels/lnInclude/sedIncompressibleTurbulenceModel.C:128:12: error: ‘New’ is not a member of ‘Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >’
return tmp<surfaceScalarField>::New
^
/home/sourav/OpenFOAM/OpenFOAM-6/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt64Opt/sedFoam_MULES.o' failed
make: *** [Make/linux64GccDPInt64Opt/sedFoam_MULES.o] Error 1

Please help me out. Thanks in advance.

[sourav8016]

Dear Sir,

I have been trying to install sedfoam. I am using openfoam6 in ubuntu 16.04. While running ./Allwmake command, I get the following errors

Foam::sedIncompressibleTurbulenceModel<TransportMo del>:Primef() const’:
../TurbulenceModels/turbulenceModels/lnInclude/sedIncompressibleTurbulenceModel.C:128:12: error: ‘New’ is not a member of ‘Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> >’
return tmp<surfaceScalarField>::New
^
/home/sourav/OpenFOAM/OpenFOAM-6/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt64Opt/sedFoam_MULES.o' failed
make: *** [Make/linux64GccDPInt64Opt/sedFoam_MULES.o] Error 1

Please help me out. Thanks in advance.

[cyss38]

Sorry, we are now only supporting the versions of OpenFoam.com.

I have to update the readme and the doc to clarify this point


Cyrille

[sourav8016]

Originally Posted by cyss38

Sorry, we are now only supporting the versions of OpenFoam.com.

I have to update the readme and the doc to clarify this point


Cyrille

Is sedfoam compatible with OpenFOAM-v2012 ?

[cyss38]

Originally Posted by sourav8016

Is sedfoam compatible with OpenFOAM-v2012 ?

yes, sedfoam is compatible with OpenFOAM-v2012

[qi.yang@polimi.it]

Originally Posted by julienC

Dear all,

The third release of sedFoam, a two-phase flow model for sediment transport applications, is available on github:
https://github.com/SedFoam/sedfoam

and an open-access article describing the model capabilities is available here:
https://www.geosci-model-dev.net/10/...-4367-2017.pdf

A short introduction and the description of the tutorials can be found here:
http://servforge.legi.grenoble-inp.fr/pub/soft-sedfoam/

This thread is to inform the community about the recent evolutions of the code as well as to gather feedbacks, questions and comments from the community. Don't hesitate!

All the best,

Julien Chauchat

Dear Julien,


Thanks for your updating. But I always met a problem when I was compiling the solver. Could you please give me some starting tutorials?

[cyss38]

Which version of Openfoam do you use?
As said before, we now support openfoam.com (and not openfoam.org)...

[qi.yang@polimi.it]

Originally Posted by cyss38

Which version of Openfoam do you use?
As said before, we now support openfoam.com (and not openfoam.org)...

Hi Thanks for your reply. In the end, I compiled the solver successfully.

[qi.yang@polimi.it]

Dear all,

Please help me see the issue of leading wrong fields set (refer to this threadthe inlet proile is not numbers while using funkySetFieldsDict). Thanks.

[qi.yang@polimi.it]

already solved.

[qi.yang@polimi.it]

Dear all,

"After I finished the calculation, should I use another tool to plot the scour characteristics, instead of paraview? Now I can only plot the distribution of solid concentration to roughly see the scouring process. Do you have any idea? Thanks a lot." solved

[cyss38]

Originally Posted by qi.yang@polimi.it

Dear all,

"After I finished the calculation, should I use another tool to plot the scour characteristics, instead of paraview? Now I can only plot the distribution of solid concentration to roughly see the scouring process. Do you have any idea? Thanks a lot." solved

Sorry for the response time...
Yes You can use Python to postprocess and display your results.

In the Tutorials/Py folder, there are some examples that are explained in the sedfoam documentation https://sedfoam.github.io.

Note that Sedfoam has evolved a lot in recent months (thanks to the developers of OpenFoam.com for their help and advice)

[qi.yang@polimi.it]

Originally Posted by cyss38

Sorry for the response time...
Yes You can use Python to postprocess and display your results.

In the Tutorials/Py folder, there are some examples that are explained in the sedfoam documentation https://sedfoam.github.io.

Note that Sedfoam has evolved a lot in recent months (thanks to the developers of OpenFoam.com for their help and advice)

Thanks a lot. I found I could deal with it via paraview directly.
I have another question that always makes me confused that the flow velocity in the 3DScour tutorial case is not specified since cyclic BC was used in the inlet and outlet in the 1D case. Thus, if I require exploring the effect of different flow velocities, how can it be achieved?

[cyss38]

Originally Posted by qi.yang@polimi.it

I have another question that always makes me confused that the flow velocity in the 3DScour tutorial case is not specified since cyclic BC was used in the inlet and outlet in the 1D case. Thus, if I require exploring the effect of different flow velocities, how can it be achieved?

As inlet boundary condition, the velocity of both phases, the sediment concentration, the TKE and the TKE dissipation are set based on one-dimensional simulation results using funkySetFields.


The associated one-dimensional case can be found in 1D subdirectory and the fluidfoam python package is used in order to generate the 1d_profil files which are necessary for funkySetFields (see Allrun.pre and Allrun).

So, if you want to explore the effect of different flow velocities you have to change the velocity in the 1D case. In this 1D case, the flow is driven by a pressure gradient, that you can change via the variable gradPMean of the constant/forceProperties file (of the 1D case) : https://github.com/SedFoam/sedfoam/b...Properties#L19

[qi.yang@polimi.it]

Originally Posted by cyss38

As inlet boundary condition, the velocity of both phases, the sediment concentration, the TKE and the TKE dissipation are set based on one-dimensional simulation results using funkySetFields.


The associated one-dimensional case can be found in 1D subdirectory and the fluidfoam python package is used in order to generate the 1d_profil files which are necessary for funkySetFields (see Allrun.pre and Allrun).

So, if you want to explore the effect of different flow velocities you have to change the velocity in the 1D case. In this 1D case, the flow is driven by a pressure gradient, that you can change via the variable gradPMean of the constant/forceProperties file (of the 1D case) : https://github.com/SedFoam/sedfoam/b...Properties#L19

I got the idea. Thanks. So that's to say, I can roughly calculate the inlet velocity by the expression (1/2*rho*V^2=P) if I specified a mean pressure in the dictionary. Is it correct?

[mathiant]

Originally Posted by qi.yang@polimi.it

I got the idea. Thanks. So that's to say, I can roughly calculate the inlet velocity by the expression (1/2*rho*V^2=P) if I specified a mean pressure in the dictionary. Is it correct?

You can not directly target the mean velocity in a boundary layer by changing the pressure gradient gradPMEAN in constant/forceProperties.

Wath you can do is to change gradPMEAN to match a given bottom friction velocity "utau" following: gradPMEAN = utau^2/h

Regards,


Antoine

[qi.yang@polimi.it]

Originally Posted by mathiant

You can not directly target the mean velocity in a boundary layer by changing the pressure gradient gradPMEAN in constant/forceProperties.

Wath you can do is to change gradPMEAN to match a given bottom friction velocity "utau" following: gradPMEAN = utau^2/h

Regards,


Antoine

Thank you, Antoine. I must know the relationship regarding mean velocity. Where can I find the relevant literature or equations?

[ms.hashempour]

Dear friends

I can not resume my running with "mpirun -np 4 sedfoam -parallel"
while I change my controldict for continuing based on the latest time of run.
the error is:
"mpirun was unable to find the specified executable file, and therefore
did not launch the job. This error was first reported for process
rank 0; it may have occurred for other processes as well."


hope to guide me

[cyss38]

Originally Posted by ms.hashempour

Dear friends

I can not resume my running with "mpirun -np 4 sedfoam -parallel"
while I change my controldict for continuing based on the latest time of run.
the error is:
"mpirun was unable to find the specified executable file, and therefore
did not launch the job. This error was first reported for process
rank 0; it may have occurred for other processes as well."


hope to guide me

The name of the binary is sedFoam_rbgh and not sedfoam.
So, just type :
"mpirun -np 4 sedFoam_rbgh -parallel"

[qi.yang@polimi.it]

Hi,

I am simulating a scouring problem in terms of my own structure. However, I always encountered some convergence issues and the results seem crazy. Due to the fact that I adopted ansys workbench to create geometry and mesh, also I specified symmetry BC on two sides of the domain, instead of cyclic used in the tutorial case. I am not sure that boundary condition will result in the wrong solution (see attachment, alpha.a at time 1 s). Alternatively, I tried to change cyclic to symmetry based on the original 3DScour case in the tutorial I also met the same problem but after I adjusted the cell number in the vertical direction, which improved the solution. I get stuck here. Do you have any ideas or comments? Appreciated.