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

Thermophysical transport model to include shear rate

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 29, 2016, 14:19
Default Thermophysical transport model to include shear rate
  #1
New Member
 
Janez Turk
Join Date: Jul 2010
Posts: 6
Rep Power: 12
Janez Turk is on a distinguished road
Greetings!

I'm currently working on new thermophysical model with its own equation of state and transport properties. I get stuck when I try to implement shear rate dependence of viscosity. Is this feasible? Where would you calculate shear rate and pass it to viscosity function? Goal is to capture shear thinning behaviour of polymers with multiphase flow.

This is where ("/src/thermophysicalModels/specie/transport/crossWLF/crossWLFTransportI.H") viscosity is evaluated. It takes p and T.
Code:
template<class Thermo>
inline Foam::scalar Foam::crossWLFTransport<Thermo>::mu
(
    const scalar p,
    const scalar T
) const
{
    return D1_*::exp(-(A1_*(T-D2_))/(A2_+T-D2_));
}
This is file /src/thermophysicalModels/basic/rhoThermo/heRhoThermo.H, where I suppose mu is called.

Code:
template<class BasicPsiThermo, class MixtureType>
void Foam::heRhoThermo<BasicPsiThermo, MixtureType>::calculate()
{
    const scalarField& hCells = this->he().internalField();
    const scalarField& pCells = this->p_.internalField();
    //const vectorField& UCells = this->U_.internalField(); ofcourse it doesn't work like this

    scalarField& TCells = this->T_.internalField();
    scalarField& psiCells = this->psi_.internalField();
    scalarField& rhoCells = this->rho_.internalField();
    scalarField& muCells = this->mu_.internalField();
    scalarField& alphaCells = this->alpha_.internalField();

    forAll(TCells, celli)
    {
        const typename MixtureType::thermoType& mixture_ =
            this->cellMixture(celli);

        TCells[celli] = mixture_.THE
        (
            hCells[celli],
            pCells[celli],
            TCells[celli]
        );

        psiCells[celli] = mixture_.psi(pCells[celli], TCells[celli]);
        rhoCells[celli] = mixture_.rho(pCells[celli], TCells[celli]);

        muCells[celli] = mixture_.mu(pCells[celli], TCells[celli]); // is this the line where mu is called?
        alphaCells[celli] = mixture_.alphah(pCells[celli], TCells[celli]);
    }
Thank you!
Janez Turk is offline   Reply With Quote

Old   March 6, 2018, 04:44
Default
  #2
Member
 
Alberto
Join Date: Sep 2013
Posts: 36
Rep Power: 9
malv83 is on a distinguished road
Hi Janez,

Were you able to implement this modification?

Thanks
malv83 is offline   Reply With Quote

Old   March 18, 2018, 10:41
Default
  #3
New Member
 
Janez Turk
Join Date: Jul 2010
Posts: 6
Rep Power: 12
Janez Turk is on a distinguished road
Hello sir,

yes we managed, source code can be downloaded:

https://github.com/krebeljk/openInjMoldSim

Regards,
Janez Turk
Janez Turk is offline   Reply With Quote

Old   March 24, 2018, 12:30
Default
  #4
Senior Member
 
Derek Mitchell
Join Date: Mar 2014
Location: UK, Reading
Posts: 163
Rep Power: 9
derekm is on a distinguished road
Is there a way to create a new thermomodel without having to modify the solver to reference it? I'm trying to create a modified version of heSolidThermo but cant see a way round modifying the solvers to reference a new version solidThermo as well
__________________
A CHEERING BAND OF FRIENDLY ELVES CARRY THE CONQUERING ADVENTURER OFF INTO THE SUNSET
derekm is offline   Reply With Quote

Reply

Tags
multiphase shear, multiphase shear thinning, shear, shear thinning, thermophysical

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
How to determine the mean filtered strain rate for the Smagorinsky model levinperson Main CFD Forum 5 November 30, 2015 14:00
Superlinear speedup in OpenFOAM 13 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 06:36
Problems in compiling paraview in Suse 10.3 platform chiven OpenFOAM Installation 3 December 1, 2009 08:21
OpenFOAM15 paraFoam bug koen OpenFOAM Bugs 19 June 30, 2009 11:46
Species Transport model Vs discrete phase Karthick FLUENT 0 April 27, 2004 03:23


All times are GMT -4. The time now is 02:26.