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

Centripetal forces on particles in incompressible fluid

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

Like Tree1Likes
  • 1 Post By vikramaditya91

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 13, 2012, 11:15
Default Centripetal forces on particles in incompressible fluid
  #1
New Member
 
Sam
Join Date: Dec 2012
Posts: 1
Rep Power: 0
sam1984_b is on a distinguished road
Hi all,

I am trying to simulate very small oil droplets in water, in a rotating reference frame. I would like to use icoUncoupledKinematicParcelFoam, as the problem is an incompressible one and I would like to uncouple the particles from the flow. However, I am struggling to implement the centripetal/Coriolis forces on the particle.

As far as I can see, my options for particle forces include NonInertialFrame and SRF, neither of which seem to work.

When I put in a NonInertialFrame in the kinematicCloudProperties, as

particleForces
{ nonInertialFrame
{omega (0 0 180);
centreOfRotation (0 0 0);
}}


It compiles and runs but the particles stay perfectly still! It's always worrying when I can replace that code with nonInertialFrame {blah blah;} and get the same result.

When I try and use SRF I get the error message:
request for SRFModel SRFProperties from objectRegistry region0 failed
available objects of type SRFModel are 0 ( )

So I guess the solver doesn't know to supply the SRF data. Please could someone advise, I would be happy fixing either of the above, or putting in a custom code to define the particle forces.

Thanks in advance,

Sam
sam1984_b is offline   Reply With Quote

Old   July 29, 2016, 04:43
Default I have the same problem as above
  #2
New Member
 
eu sou cfd
Join Date: Jun 2012
Location: Brazil
Posts: 18
Rep Power: 13
vikramaditya91 is on a distinguished road
Hey,

I have 2 questions about the post above:

1) I am stuck in the exact same problem of not being able to provide an SRF force. I get this error when I try to do some Lagrangian particle tracking using the icoUncoupledKineticParcelFoam on a case which has flow field results from an SRFPimpleFoam solver.
From what I understand in the error, it seems that objectRegistry does not have the SRFProperties file loaded in it. I tried to go to the place where it actually loads the file and I came across this particular piece of code located in $FOAM_SRC/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF in OpenFOAM3.0.1

Code:
template<class CloudType>
void Foam::SRFForce<CloudType>::cacheFields(const bool store)
{
    if (store)
    {
        const typename SRF::SRFModel& model = this->mesh().template
            lookupObject<SRF::SRFModel>("SRFProperties");
        srfPtr_ = &model;
    }
    else
    {
        srfPtr_ = NULL;
    }
}
Do you have any suggestions on how I can let the SRFForce file know that I already have a SRFProperties file in the constant/ directory.


https://drive.google.com/file/d/0B7O...ew?usp=sharing

I have attached a modified tutorial case (2DRotor), on which you can recreate the problem by first running "SRFPimpleFoam" which will output results of velocity and pressure flow fields. Then run "icoUncoupledKineticParcelFoam" which gives the error that
Code:
request for SRFModel SRFProperties from objectRegistry region0 failed
available objects of type SRFModel are 0 ( )

2) My second question is about what Sam posted above. He said he wanted to use the NonIntertialFrame forces as well. I do not get the reasoning behind using the NonInterialFram forces because the SRFForce makes a correction of exactly that, while also taking densities of the particle and fluid also into account with this line of code:
Code:
    // Coriolis and centrifugal acceleration terms
    value.Su() =
        mass*(1.0 - p.rhoc()/p.rho())
       *(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega)));
Regards

Last edited by vikramaditya91; July 29, 2016 at 05:54.
vikramaditya91 is offline   Reply With Quote

Old   July 29, 2016, 06:15
Default Quick Update
  #3
New Member
 
eu sou cfd
Join Date: Jun 2012
Location: Brazil
Posts: 18
Rep Power: 13
vikramaditya91 is on a distinguished road
This error doesnt seem to appear now after I made these changes. It seems that the SRFProperties file wasnt even read earlier. Adding these lines to the createFields.H seems to make the error go away. I have not tested this just yet, will edit this post soon

Code:
Info<< "Reading field Urel\n" << endl;
volVectorField Urel
(
    IOobject
    (
        "Urel",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::AUTO_WRITE
    ),
    mesh
);

Info<< "Creating SRF model\n" << endl;
autoPtr<SRF::SRFModel> SRF(SRF::SRFModel::New(Urel));
vikramaditya91 is offline   Reply With Quote

Old   October 26, 2016, 06:16
Default icoUncoupledKineticParcelFoam is not for rotating frames of reference
  #4
New Member
 
eu sou cfd
Join Date: Jun 2012
Location: Brazil
Posts: 18
Rep Power: 13
vikramaditya91 is on a distinguished road
Hey, it is me again.

I posted a bug on the OpenFOAM Foundation's bug tracking facility. The issue is listed here http://bugs.openfoam.org/view.php?id=2265 .
The site manager Henry Weller, said that the icoUncoupledKineticParcelFoam is not really meant for single rotating frames of reference, so it was never meant to work with the SRFPimpleFoam fluid velocity. However, the new and edited (fully uncoupled) Lagrangian particle tracker is listed in
https://github.com/vikramaditya91/Op...cParcelSRFFoam
in case you are interested.

The only major change I made was to read Urel instead of U for particle tracking. Please let me know if there are still bugs or if it can be improved further.
cesarvenier likes this.
vikramaditya91 is offline   Reply With Quote

Old   December 14, 2018, 07:16
Default still doesn't work
  #5
New Member
 
Join Date: Nov 2017
Location: Germany
Posts: 11
Rep Power: 8
elomri is on a distinguished road
I tried the modifications suggested in vikramaditya91 post but i still get the same error as before. I'm using the 4.1 release of OF.


How doeas your KinematicProperties looks like ? vikramaditya91
elomri is offline   Reply With Quote

Reply

Tags
lagrangian particle, noninertialframe, srf

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
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 18:44
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
How to choose the mean diameter value for dispersed fluid? creddy_trddc CFX 1 October 30, 2011 05:30
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
Particles in slow dense fluid newbie Main CFD Forum 6 December 6, 2001 10:43


All times are GMT -4. The time now is 10:40.