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

Parcel angular velocity and Magnus Force

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 11, 2017, 11:08
Default Parcel angular velocity and Magnus Force
  #1
Member
 
Joaquín Neira
Join Date: Oct 2017
Posts: 38
Rep Power: 5
cojua8 is on a distinguished road
Hello,

I am trying to create a new Magnus force for particles.

for this purpose I need to access the angular velocity in CollidingParcelI.H:


Code:
template<class ParcelType>
inline Foam::vector Foam::CollidingParcel<ParcelType>::omega() const
{
    return angularMomentum_/this->momentOfInertia();
}
And the code of the force is this:

Code:
template<class CloudType>
Foam::forceSuSp Foam::MagnusForce<CloudType>::calcCoupled
(
    const typename CloudType::parcelType& p,
    const scalar dt,
    const scalar mass,
    const scalar Re,
    const scalar muc
    const vector omega
) const
{
    forceSuSp value(Zero, 0.0);

    vector curlUc =
        curlUcInterp().interpolate(p.coordinates(), p.currentTetIndices());

    scalar Cl = this->Cl(p, curlUc, Re, muc);

    value.Su() =  Cl*0.75*mass*(omega^(p.U()-p.Uc()));
    
    Info << "Calculate Lift force" << endl;

    return value;
}
the compiler compiles the new force, but when running the solver there ir no "Calculate Lift Force" output, which means the solver doesn't get into the force code.

When I delete all the "omega's" the code works correctly (but obviously it doesnt work as needed), so I guess the omega call is not okay

Last edited by cojua8; December 11, 2017 at 11:09. Reason: more clear title
cojua8 is offline   Reply With Quote

Reply

Tags
angular velocity, magnus

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 calculate swirl velocity or angular velocity in a in-cylinder flow analysis us tgsankar Main CFD Forum 2 April 26, 2017 16:55
Accessing parcel angular velocity from particel force class slint OpenFOAM Programming & Development 1 September 27, 2016 12:46
Relation between Swirl flow and angular velocity padmanathan FLUENT 0 May 4, 2011 07:18
Angular velocity Haritz Main CFD Forum 0 March 19, 2006 13:59
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 09:11


All times are GMT -4. The time now is 09:00.