[Jose Manuel Olmos]

Hello everyone,


I'm new with OpenFOAM and I would like to know if it is possible to simulate the flux inside a three-dimensional pipe of a fluid containing a solute. The fluid is introduced at the inlet of the pipe at a constant velocity and with a constant concentration of a solute. I want to simulate the distribution of velocities of the fluid and concentrations of the solute along the inside of the pipe.


I used a geometry file (in stl format) of the pipe to create the internal mesh by using the "snappyHexMesh" utility.


https://ibb.co/FmSRh5b


However, I have problems for running the simulation. The diffusion of the solute only takes place from the inlet of the pipe instead of at any point of the pipe (I mean, in sincronization with the flux of the fluid). The following images refer to the distribution of velocities and concentrations at the same time of simulation


https://ibb.co/bQLCf2v
https://ibb.co/7XSj25K


I modified the icoFoam solver to include the effect of the diffusion of the solute, according to the instructions provided in the web "https://openfoamwiki.net/index.php/How_to_add_transport_of_N_solutes_to_icoFoam". I attach the files for this new solver.


Here are the boundary conditions for the concentration of the solute:
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object C_azoxystrobin;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 -3 0 0 1 0 0];

internalField uniform 0;

boundaryField
{
CAD_inlet
{
type fixedValue;
value uniform 0.01;
}
CAD_outlet
{
type zeroGradient;
}
CAD_wall
{
type zeroGradient;
}
}


// ************************************************** *********************** //


I hope my explanations about what I want to do are clear. I would be very grateful if anyone could tell me if I'm taking any mistake (or maybe the case I want to simulate is not possible with OpenFOAM)


Thank you in advance,
Jose Manuel