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

FSI solver

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

Like Tree3Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 14, 2016, 00:25
Default
  #21
Member
 
Viraj Belekar
Join Date: Jun 2016
Posts: 68
Rep Power: 9
viraj20feb is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Quick answer: Please run the following commands:
Code:
cd $WM_PROJECT_USER_DIR
cd FluidStructureInteraction/src/
./Allwmake > log.make 2>&1
gzip < log.make > log.make.gz
Then please attach the file "log.make.gz" to your next post. I ask this because the PDF file did not have any clear indication of any problem at all.
Hi Mr. Bruno,

Please find the attached file.
Attached Files
File Type: gz log.make.gz (2.2 KB, 2 views)
viraj20feb is offline   Reply With Quote

Old   September 14, 2016, 16:11
Default
  #22
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Quick answer: The problem was that there were a few steps that failed, where I had written how to download the updated "constitutiveModel" files.
Please follow these steps instead, which hopefully will fix the problem:
  1. Download the attached file "src_fluidStructureInteraction_stressModels_consti tutiveModel.tar.gz".
  2. In the terminal (command line) where you have foam-extend 3.1 ready to be used, please run:
    Code:
    cd $WM_PROJECT_USER_DIR/FluidStructureInteraction
    
    pwd
  3. The last command above should have told you about a path. Please place the file that you downloaded from here, into that folder (directory).
  4. Then run the following command:
    Code:
    tar -xzf src_fluidStructureInteraction_stressModels_constitutiveModel.tar.gz
  5. Then let's try building the toolbox again, by running:
    Code:
    cd src
    ./Allwmake > log.make 2>&1
    gzip < log.make > log.make.gz
Once it's over, please attach the new "log.make.gz" file once more, to see if everything went well or not.
wyldckat is offline   Reply With Quote

Old   September 15, 2016, 00:33
Default
  #23
Member
 
Viraj Belekar
Join Date: Jun 2016
Posts: 68
Rep Power: 9
viraj20feb is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Quick answer: The problem was that there were a few steps that failed, where I had written how to download the updated "constitutiveModel" files.
Please follow these steps instead, which hopefully will fix the problem:
  1. Download the attached file "src_fluidStructureInteraction_stressModels_consti tutiveModel.tar.gz".
  2. In the terminal (command line) where you have foam-extend 3.1 ready to be used, please run:
    Code:
    cd $WM_PROJECT_USER_DIR/FluidStructureInteraction
    
    pwd
  3. The last command above should have told you about a path. Please place the file that you downloaded from here, into that folder (directory).
  4. Then run the following command:
    Code:
    tar -xzf src_fluidStructureInteraction_stressModels_constitutiveModel.tar.gz
  5. Then let's try building the toolbox again, by running:
    Code:
    cd src
    ./Allwmake > log.make 2>&1
    gzip < log.make > log.make.gz
Once it's over, please attach the new "log.make.gz" file once more, to see if everything went well or not.

Hi wyldckat,

Please find the attached file. There already existed an src folder at this location: cd $WM_PROJECT_USER_DIR/FluidStructureInteraction

from the previous installation. I did exactly what you said.
Attached Files
File Type: gz log.make.gz (2.2 KB, 1 views)
viraj20feb is offline   Reply With Quote

Old   September 16, 2016, 09:06
Default
  #24
Member
 
Viraj Belekar
Join Date: Jun 2016
Posts: 68
Rep Power: 9
viraj20feb is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Quick answer: The problem was that there were a few steps that failed, where I had written how to download the updated "constitutiveModel" files.
Please follow these steps instead, which hopefully will fix the problem:
  1. Download the attached file "src_fluidStructureInteraction_stressModels_consti tutiveModel.tar.gz".
  2. In the terminal (command line) where you have foam-extend 3.1 ready to be used, please run:
    Code:
    cd $WM_PROJECT_USER_DIR/FluidStructureInteraction
    
    pwd
  3. The last command above should have told you about a path. Please place the file that you downloaded from here, into that folder (directory).
  4. Then run the following command:
    Code:
    tar -xzf src_fluidStructureInteraction_stressModels_constitutiveModel.tar.gz
  5. Then let's try building the toolbox again, by running:
    Code:
    cd src
    ./Allwmake > log.make 2>&1
    gzip < log.make > log.make.gz
Once it's over, please attach the new "log.make.gz" file once more, to see if everything went well or not.

Hi wyldckat,

Please find the attached file. I did exactly what you told me
Attached Files
File Type: gz log.make.gz (2.1 KB, 1 views)
viraj20feb is offline   Reply With Quote

Old   September 17, 2016, 07:06
Default
  #25
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Hi viraj20feb,

Sorry for the late reply, but I didn't manage to find enough free time to look into this sooner.

Something strange is going on with your version of foam-extend 3.1. I need more details, namely:
  1. Which exact installation instructions did you follow for installing foam-extend 3.1?
  2. Where and how exactly did you download foam-extend 3.1?
  3. What do the following commands give you?
    Code:
    foam
    git log -1
Best regards,
Bruno
__________________
wyldckat is offline   Reply With Quote

Old   September 19, 2016, 00:43
Default
  #26
Member
 
Viraj Belekar
Join Date: Jun 2016
Posts: 68
Rep Power: 9
viraj20feb is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Hi viraj20feb,

Sorry for the late reply, but I didn't manage to find enough free time to look into this sooner.

Something strange is going on with your version of foam-extend 3.1. I need more details, namely:
  1. Which exact installation instructions did you follow for installing foam-extend 3.1?
  2. Where and how exactly did you download foam-extend 3.1?
  3. What do the following commands give you?
    Code:
    foam
    git log -1
Best regards,
Bruno

Hi Bruno,

Following are the answers of your question:

1. I am working as Research Assistant in a project and foam-extend 3.1 was already installed in the system alloted to me. So I don't know about the installation instructions, sorry!
2. Same as 1.
3.
Code:
foam
If I do not source foam -extend 3.1, this command gives me the output :
Code:
foam: command not found
After sourcing, the output is:
Code:
/foam/foam-extend-3.1$
In both the cases, the command
Code:
git log -1
gives the output:
Code:
fatal: Not a git repository (or any of the parent directories): .git
Thanks.
viraj20feb is offline   Reply With Quote

Old   September 26, 2016, 18:09
Default
  #27
Retired Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,974
Blog Entries: 45
Rep Power: 128
wyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to allwyldckat is a name known to all
Quick answer: My apologies for taking so long to give an answer. Since there aren't enough details about the specific snapshot of foam-extend 3.1, it becomes considerably harder to diagnose which version it is.
Running the following commands might help to diagnose, but it's not 100% guaranteed it will work:
Code:
foam
cat .build
For example, for me, it gives me this:
Code:
3.1-1dd681f6e943
wyldckat is offline   Reply With Quote

Old   March 12, 2019, 10:05
Default
  #28
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7
Hgholami is on a distinguished road
Dear Samkhaniani
Did you simulate a non-Newtonian flow with flexible pipe?
I want to simulate non-newtonian fluid with fsiFoam. I was modified icoFsiElasticNonLinULSolidFoam in createField and fluidSolver and was seen correct. But In this Solver I can't do that. Can you favor me? Thanks
Hgholami is offline   Reply With Quote

Old   March 13, 2019, 11:16
Default
  #29
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
nimasam is on a distinguished road
could you please tell me what problems are you dealing with?



Best regards
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)
nimasam is offline   Reply With Quote

Old   March 14, 2019, 08:32
Default
  #30
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7
Hgholami is on a distinguished road
Dear Samkhaniani
I uses fsiFoam as the solver for FSI problem. As you know, this solver has fluidSolver such as icoFluid, PisoFluid and ConsistIcoFluid for solving Newtonian fluid. I want to solve non-Newtonian fluid with fsiFoam.
I modified simply some code of icoFsiElasticNonLinULSolidFoam in createField and fluidSolver to solve now-Newtonian fluid. Only I changed reading nu in solvers. I do that change in this solver but have error. the icoNonFluid.C is
Code:
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | foam-extend: Open Source CFD
   \\    /   O peration     | Version:     4.0
    \\  /    A nd           | Web:         http://www.foam-extend.org
     \\/     M anipulation  | For copyright notice see file Copyright
-------------------------------------------------------------------------------
License
    This file is part of foam-extend.

    foam-extend is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by the
    Free Software Foundation, either version 3 of the License, or (at your
    option) any later version.

    foam-extend is distributed in the hope that it will be useful, but
    WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.

\*---------------------------------------------------------------------------*/

#include "icoNonFluid.H"
#include "singlePhaseTransportModel.H"
#include "volFields.H"
#include "fvm.H"
#include "fvc.H"
#include "fvMatrices.H"
#include "addToRunTimeSelectionTable.H"
#include "findRefCell.H"
#include "adjustPhi.H"
#include "fluidSolidInterface.H"
#include "fixedGradientFvPatchFields.H"

#include "EulerDdtScheme.H"
#include "CrankNicolsonDdtScheme.H"
#include "backwardDdtScheme.H"

#include "elasticSlipWallVelocityFvPatchVectorField.H"
#include "elasticWallVelocityFvPatchVectorField.H"
#include "elasticWallPressureFvPatchScalarField.H"

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

namespace Foam
{
namespace fluidSolvers
{

// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //

defineTypeNameAndDebug(icoNonFluid, 0);
addToRunTimeSelectionTable(fluidSolver, icoNonFluid, dictionary);


// * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //

icoNonFluid::icoNonFluid(const fvMesh& mesh)
:
    fluidSolver(this->typeName, mesh),
    U_
    (
        IOobject
        (
            "U",
            runTime().timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    ),
    p_
    (
        IOobject
        (
            "p",
            runTime().timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    ),

    gradp_(fvc::grad(p_)),
    gradU_(fvc::grad(U_)),

    phi_
    (
        IOobject
        (
            "phi",
            runTime().timeName(),
            mesh,
            IOobject::READ_IF_PRESENT,
            IOobject::AUTO_WRITE
        ),
        fvc::interpolate(U_) & mesh.Sf()
    ),
    singlePhaseTransportModel fluid(U, phi); //the new code
    transportProperties_
    (
        IOobject
        (
            "transportProperties",
            runTime().constant(),
            mesh,
            IOobject::MUST_READ,
            IOobject::NO_WRITE
        )
    ),
//    nu_(transportProperties_.lookup("nu")),
     rho_( transportProperties.lookup("rho"))

{}

// * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //

const volVectorField& icoNonFluid::U() const
{
    return U_;
}


const volScalarField& icoNonFluid::p() const
{
    return p_;
}


const volScalarField& icoNonFluid::nu() const
{
    return nu_;
}

//- Patch viscous force (N/m2)
tmp<vectorField> icoNonFluid::patchViscousForce(const label patchID) const
{
    tmp<vectorField> tvF
    (
        new vectorField(mesh().boundary()[patchID].size(), vector::zero)
    );

    tvF() = rho_.value()*nu().value()*U().boundaryField()[patchID].snGrad();

    vectorField n = mesh().boundary()[patchID].nf();
    tvF() -= n*(n&tvF());

    return tvF;
}

//- Patch pressure force (N/m2)
tmp<scalarField> icoNonFluid::patchPressureForce(const label patchID) const
{
    tmp<scalarField> tpF
    (
        new scalarField(mesh().boundary()[patchID].size(), 0)
    );

    tpF() = rho_.value()*p().boundaryField()[patchID];

    return tpF;
}

//- Patch viscous force (N/m2)
tmp<vectorField> icoNonFluid::faceZoneViscousForce
(
    const label zoneID,
    const label patchID
) const
{
    vectorField pVF = patchViscousForce(patchID);

    tmp<vectorField> tvF
    (
        new vectorField(mesh().faceZones()[zoneID].size(), vector::zero)
    );
    vectorField& vF = tvF();

    const label patchStart =
        mesh().boundaryMesh()[patchID].start();

    forAll(pVF, i)
    {
        vF[mesh().faceZones()[zoneID].whichFace(patchStart + i)] =
            pVF[i];
    }

    // Parallel data exchange: collect pressure field on all processors
    reduce(vF, sumOp<vectorField>());


    return tvF;
}

//- Patch pressure force (N/m2)
tmp<scalarField> icoNonFluid::faceZonePressureForce
(
    const label zoneID,
    const label patchID
) const
{
    scalarField pPF = patchPressureForce(patchID);

    tmp<scalarField> tpF
    (
        new scalarField(mesh().faceZones()[zoneID].size(), 0)
    );
    scalarField& pF = tpF();

    const label patchStart =
        mesh().boundaryMesh()[patchID].start();

    forAll(pPF, i)
    {
        pF[mesh().faceZones()[zoneID].whichFace(patchStart + i)] =
            pPF[i];
    }

    // Parallel data exchange: collect pressure field on all processors
    reduce(pF, sumOp<scalarField>());

    return tpF;
}

tmp<scalarField> icoNonFluid::faceZoneMuEff
(
    const label zoneID,
    const label patchID
) const
{
    tmp<scalarField> tMuEff
    (
        new scalarField
        (
            mesh().faceZones()[zoneID].size(),
            rho_.value()*nu_.value()
        )
    );

    return tMuEff;
}

void icoNonFluid::evolve()
{
    Info << "Evolving fluid solver: " << this->type() << endl;

    const fvMesh& mesh = fluidSolver::mesh();

    int nCorr(readInt(fluidProperties().lookup("nCorrectors")));

    int nNonOrthCorr =
        readInt(fluidProperties().lookup("nNonOrthogonalCorrectors"));

    int nOuterCorr =
        readInt(fluidProperties().lookup("nOuterCorrectors"));

    // Prepare for the pressure solution
    label pRefCell = 0;
    scalar pRefValue = 0.0;
    setRefCell(p_, fluidProperties(), pRefCell, pRefValue);

    for (int oCorr = 0; oCorr < nOuterCorr; oCorr++)
    {
        if(mesh.moving())
        {
            // Make the fluxes relative
            phi_ -= fvc::meshPhi(U_);
        }

        // CourantNo
        {
          scalar CoNum = 0.0;
          scalar meanCoNum = 0.0;
          scalar velMag = 0.0;

          if (mesh.nInternalFaces())
          {
              surfaceScalarField SfUfbyDelta =
                  mesh.surfaceInterpolation::deltaCoeffs()*mag(phi_);

              CoNum = max(SfUfbyDelta/mesh.magSf())
                .value()*runTime().deltaT().value();

              meanCoNum = (sum(SfUfbyDelta)/sum(mesh.magSf()))
                .value()*runTime().deltaT().value();

              velMag = max(mag(phi_)/mesh.magSf()).value();
          }

          Info<< "Courant Number mean: " << meanCoNum
              << " max: " << CoNum
              << " velocity magnitude: " << velMag << endl;
        }

        fvVectorMatrix UEqn
        (
            fvm::ddt(U_)
          + fvm::div(phi_, U_)
          - fvm::laplacian(nu_, U_)
        );

        solve(UEqn == -gradp_);    

        // --- PISO loop

        volScalarField rAU = 1.0/UEqn.A();
        surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU));

        for (int corr=0; corr<nCorr; corr++)
        {
            U_ = rAU*UEqn.H();
            phi_ = (fvc::interpolate(U_) & mesh.Sf());
//             + fvc::ddtPhiCorr(rUA, U_, phi_);

#           include "updateRobinFsiInterface.H"

            adjustPhi(phi_, U_, p_);

            for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
            {
                fvScalarMatrix pEqn
                (
                    fvm::laplacian
                    (
                        rAUf, p_, "laplacian((1|A(U)),p)"
                    )
                 == fvc::div(phi_)
                   // fvm::laplacian(rAUf, p_) == fvc::div(phi_)
                );

                pEqn.setReference(pRefCell, pRefValue);
                pEqn.solve();

                gradp_ = fvc::grad(p_);

                if (nonOrth == nNonOrthCorr)
                {
                    phi_ -= pEqn.flux();
                }
            }

            // Continuity error
            {
                volScalarField contErr = fvc::div(phi_);

                scalar sumLocalContErr = runTime().deltaT().value()*
                    mag(contErr)().weightedAverage(mesh.V()).value();

                scalar globalContErr = runTime().deltaT().value()*
                    contErr.weightedAverage(mesh.V()).value();

                Info<< "time step continuity errors : sum local = "
                    << sumLocalContErr << ", global = "
                    << globalContErr << endl;
            }

            U_ -= rAU*gradp_;
            U_.correctBoundaryConditions();

            gradU_ = fvc::grad(U_);
        }
    }
}


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

} // End namespace fluidSolvers
} // End namespace Foam

// ************************************************************************* //
and the icoNonFluid.H is
Code:
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | foam-extend: Open Source CFD
   \\    /   O peration     | Version:     4.0
    \\  /    A nd           | Web:         http://www.foam-extend.org
     \\/     M anipulation  | For copyright notice see file Copyright
-------------------------------------------------------------------------------
License
    This file is part of foam-extend.

    foam-extend is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by the
    Free Software Foundation, either version 3 of the License, or (at your
    option) any later version.

    foam-extend is distributed in the hope that it will be useful, but
    WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
    General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.

Class
    icoNonFluid

Description
    icoFoam fluid solver

Author
    Hrvoje Jasak, Wikki Ltd.  All rights reserved. modified

SourceFiles
    icoNonFluid.C

\*---------------------------------------------------------------------------*/

#ifndef icoNonFluid_H
#define icoNonFluid_H

#include "fluidSolver.H"
#include "volFields.H"
#include "surfaceFields.H"

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

namespace Foam
{
namespace fluidSolvers
{

/*---------------------------------------------------------------------------*\
                        Class icoNonFluid Declaration
\*---------------------------------------------------------------------------*/

class icoNonFluid
:
    public fluidSolver
{
    // Private data

        //- Velocity field
        volVectorField U_;

        //- Pressure field
        volScalarField p_;


      //- viscousity field
        volScalarField nu_;

        //- Pressure gradient
        volVectorField gradp_;

        //- Velocity gradient
        volTensorField gradU_;

        //- Flux field
        surfaceScalarField phi_;

        // Transport properties dictionary
        IOdictionary transportProperties_;

        //- Density
        dimensionedScalar rho_;


//         //-
//         scalarField nGradUn_;

    // Private Member Functions

        //- Disallow default bitwise copy construct
        icoNonFluid(const icoNonFluid&);

        //- Disallow default bitwise assignment
        void operator=(const icoNonFluid&);


public:

    //- Runtime type information
    TypeName("icoNonFluid");

    // Constructors

        //- Construct from components
        icoNonFluid(const fvMesh& mesh);

//         //- Construct from components
//         icoNonFluid
//         (
//             const word& type,
//             const fvMesh& mesh
//         );

    // Destructor

        virtual ~icoNonFluid()
        {}


    // Member Functions

        // Access

            //- Return velocity field
            virtual const volVectorField& U() const;

            //- Return velocity field
            volVectorField& U()
            {
                return U_;
            }

            //- Return pressure field
            virtual const volScalarField& p() const;

            //- Return pressure field
            volScalarField& p()
            {
                return p_;
            }

            // return viscousity field
            virtual const volScalarField& nu() const;
            //- Return viscousity field
            volScalarField& nu()
            {
                return nu_;
            }
            //- Return pressure gradient
            volVectorField& gradp()
            {
                return gradp_;
            }

            //- Return pressure gradient
            volTensorField& gradU()
            {
                return gradU_;
            }
            //- Return flux field
            surfaceScalarField& phi()
            {
                return phi_;
            }

            //- Density
            const dimensionedScalar& rho()
            {
                return rho_;
            }

            //- Patch viscous force (N/m2)
            virtual tmp<vectorField> patchViscousForce
            (
                const label patchID
            ) const;

            //- Patch pressure force (N/m2)
            virtual tmp<scalarField> patchPressureForce
            (
                const label patchID
            ) const;

            //- Patch viscous force (N/m2)
            virtual tmp<vectorField> faceZoneViscousForce
            (
                const label zoneID,
                const label patchID
            ) const;

            //- Patch pressure force (N/m2)
            virtual tmp<scalarField> faceZonePressureForce
            (
                const label zoneID,
                const label patchID
            ) const;

            //- Face zone effective dynamic viscosity
            virtual tmp<scalarField> faceZoneMuEff
            (
                const label zoneID,
                const label patchID
            ) const;

        // Edit

            //- Evolve the fluid solver
            virtual void evolve();
};


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

} // End namespace fluidSolvers
} // End namespace Foam

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

#endif

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

Last edited by wyldckat; March 16, 2019 at 18:34. Reason: [QUOTE]->[CODE]
Hgholami is offline   Reply With Quote

Old   March 14, 2019, 08:36
Default
  #31
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7
Hgholami is on a distinguished road
first error is from "singlePhaseTransportModel fluid(U, phi);" type is not direct base of "icoNonFluid"
Hgholami is offline   Reply With Quote

Old   March 14, 2019, 13:00
Default
  #32
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
nimasam is on a distinguished road
no need to write new solver, use pisoFlow as your fluid solver, then you can use non-Newtonian fluid
Hgholami likes this.
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)
nimasam is offline   Reply With Quote

Old   March 16, 2019, 23:31
Default
  #33
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7
Hgholami is on a distinguished road
Dear Samkhaniani
Thanks for your comment. I used the pisoFluid, and worked with non-Newtonian fluid. But the problem have diverging. I will try to change fvscheme for the converging. thanks again
Hgholami is offline   Reply With Quote

Old   March 20, 2019, 23:31
Default
  #34
Senior Member
 
Hojatollah Gholami
Join Date: Jan 2019
Posts: 171
Rep Power: 7
Hgholami is on a distinguished road
Dear Samkhaniani
I changed fvscheme to converge the solution, but still I unsuccessful.
I also modify the icoFluid to add add nu() as non-Newtonian fluid in that.
the problem is in here "part of icoFluid"
Quote:
tmp<scalarField> icoNonFluid::faceZoneMuEff
(
const label zoneID,
const label patchID
) const
{

tmp<scalarField> tMuEff
(
new scalarField
(
mesh().faceZones()[zoneID].size(),
rho_.value()*nu_.value() //Here
)
);

return tMuEff;
}
void icoNonFluid::evolve()
I want to change "nu_.value" with
Quote:
fluid_.nu().boundaryField()[patchID]
I add non-Newtonian fluid as fluid_(U_, phi_), then to call that I used fluid_nu.
the compiler give this error.
Quote:
fluidSolvers/icoNonFluid/icoNonFluid.C: In member function ‘virtual Foam::tmp<Foam::Field<double> > Foam::fluidSolvers::icoNonFluid::faceZoneMuEff(Foa m::label, Foam::label) const’:
fluidSolvers/icoNonFluid/icoNonFluid.C:285:9: error: no matching function for call to ‘Foam::Field<double>::Field(Foam::label, Foam::tmp<Foam::Field<double> >)’
)
^
Do you know what is the problem?
Thanks
Hgholami is offline   Reply With Quote

Old   March 26, 2019, 09:16
Default
  #35
Senior Member
 
Nima Samkhaniani
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,266
Blog Entries: 1
Rep Power: 24
nimasam is on a distinguished road
it seems that you used the wrong constructor to make your object!
maybe you need to do something like this (im not sure )

Code:
tmp<scalarField> icoNonFluid::faceZoneMuEff
(
    const label zoneID,
    const label patchID
) const
{
 
    tmp<scalarField> tMuEff
      (
        new scalarField
        (
            mesh().faceZones()[zoneID].size(),
            rho_.value()   //Here
        )
    );

tMuff=tMuff*fluid_.nu().boundaryField()[patchID];
     return tMuEff;
}
void icoNonFluid::evolve()
Hgholami likes this.
__________________
My Personal Website (http://nimasamkhaniani.ir/)
Telegram channel (https://t.me/cfd_foam)
nimasam is offline   Reply With Quote

Reply

Tags
flexible pipe, fsi

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
different results between serial solver and parallel solver wlt_1985 FLUENT 11 October 12, 2018 08:23
error in two way fsi kmgraju CFX 1 May 2, 2011 02:32
Working directory via command line Luiz CFX 4 March 6, 2011 20:02
FSI Negative Element Volume Fatal Error Aaron CFX 2 July 1, 2008 20:48
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08


All times are GMT -4. The time now is 23:46.