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

[isoAdvector] Issues with modification of IsoAdvector Code

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

Like Tree6Likes
  • 1 Post By roenby
  • 1 Post By roenby
  • 3 Post By roenby
  • 1 Post By roenby

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 16, 2017, 17:13
Default Issues with modification of IsoAdvector Code
  #1
Member
 
Raunak Bardia
Join Date: Jan 2015
Posts: 32
Rep Power: 11
raunakbardia is on a distinguished road
Hello,

I came across the release of IsoAdvector code and have recently started testing this code.

I had a question about the UEqn.H file in the code. In this file, there is a reference to UEqnSolved.H.

Quote:
Code:
...
        phi = uFactor*phi0;
        U = uFactor*U0;
    }
}
else
{
    #include "UEqnSolved.H"
}
I am unable to find this header file anywhere in the entire directory of isoAdvector and am wondering where this is invoked from.

I tried to change the isoAdvector code for a test case and rightfully, an error is thrown by the build due to its inability to find the UEqnSolved.H file.

Quote:
In file included from kappaFlow.C:109:0:
UEqn.H:38:28: fatal error: UEqnSolved.H: No such file or directory
#include "UEqnSolved.H"
^
compilation terminated.
make: *** [Make/linux64GccDPInt32Opt/kappaFlow.o] Error 1
Given that I cannot find this header file in my original isoAdvector download either is interesting; because I have been able to compile the original code and use it for some other test cases.

If someone involved in the development of isoAdvector code can answer this question, it will really speed up my testing of the code.

Thank you.
raunakbardia is offline   Reply With Quote

Old   August 23, 2017, 04:05
Default
  #2
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29
akidess will become famous soon enough
My guess is you are working on a broken version.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   August 23, 2017, 15:45
Default
  #3
Member
 
Raunak Bardia
Join Date: Jan 2015
Posts: 32
Rep Power: 11
raunakbardia is on a distinguished road
This is what I first thought. But, even on the git repository of isoAdvector I see the same thing. The file, UEqn.H includes a header file UEqnSolved.H, if we don't use a prescribed velocity.
This header file does not exist in the repository, nor in the entire OpenFOAM v1612+ directory, which is the version I am using.

Work Around:
Because I want to solve for U & P in my dynamic test cases, I simply copied the original interFoam header file for velocity and the solver compiled with the isoAdvector advection scheme without a problem.

However, I am still at a loss on how the UEqn.H that comes from IsoAdvector git repository, works?

Thank you.
raunakbardia is offline   Reply With Quote

Old   August 24, 2017, 03:47
Default
  #4
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29
akidess will become famous soon enough
There is no include for UEqnSolved.H in the foam-extend tree, and neither in the v1712 release, and neither in previous versions of the github repository you are checking. Thus I still believe it's a bug in the latest repository commit.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   August 24, 2017, 17:57
Default
  #5
Member
 
Raunak Bardia
Join Date: Jan 2015
Posts: 32
Rep Power: 11
raunakbardia is on a distinguished road
Thank you for your response.
I did not realize that the new release has got the interIsoFoam solver built-in.
Until I shift to this version of OpenFOAM, which is dependent on my system admin, I will look for the older repositories of the code.
raunakbardia is offline   Reply With Quote

Old   August 25, 2017, 04:28
Default
  #6
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 29
akidess will become famous soon enough
If you don't want to wait just get the app from the repository develop.openfoam.com and compile it yourself.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   September 11, 2017, 19:31
Default
  #7
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 92
Rep Power: 20
roenby will become famous soon enough
Hi Raunak

If you install isoAdvector using the Allwmake script from github.com/isoAdvector it will copy the interFoam solver from OpenFOAM-x.y.z/applications/solvers/multiphase/interFoam to the isoAdvector installation directory and modify the files relating to the interface advection step before compiling it into the interFlow solver. In other words the UEqn.H file included in the github.com/isoAdvector source code is never used. I will remove it to avoid confusion.

Have a look in the README.md file for installation and usage. Also feel free to report an issue in the github repo. This may save yourself and others a lot of time :-)

Best,
Johan
raunakbardia likes this.
roenby is offline   Reply With Quote

Old   September 11, 2017, 20:25
Default
  #8
Member
 
Raunak Bardia
Join Date: Jan 2015
Posts: 32
Rep Power: 11
raunakbardia is on a distinguished road
Thank you, Johan.

I started to do that a while back but forgot to update this thread about that.

You have further clarified it for all of us.
raunakbardia is offline   Reply With Quote

Old   February 27, 2020, 08:02
Default installing IsoAdvector
  #9
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 21
Rep Power: 6
Arghavani is on a distinguished road
Hello everyone,

I cannot get how we can have isoAdvector while the OpenFoam 5.x is installed?

and what would happen to the tutorials in interFoam folder, all will be deleted?

kind regards,
Arghavan
Arghavani is offline   Reply With Quote

Old   February 27, 2020, 08:15
Default
  #10
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 92
Rep Power: 20
roenby will become famous soon enough
Hi Arghavani


Go to https://github.com/isoAdvector/isoAdvector and follow the installation instructions in the README file.


Installation will not touch or delete interFoam (the new solver is called interFlow).



Best,
Johan
Aabadani likes this.
roenby is offline   Reply With Quote

Old   February 27, 2020, 08:25
Default
  #11
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 21
Rep Power: 6
Arghavani is on a distinguished road
Hi Johan,

Thank you for your reply.
I did it once but It damaged my OpenFoam and I had to reinstall it.
But I will ask my boss to help me at this time.

kind regards,
Arghavan
Arghavani is offline   Reply With Quote

Old   February 27, 2020, 08:54
Default
  #12
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 92
Rep Power: 20
roenby will become famous soon enough
That sounds terrible and should not be possible.
I would be very interested to get more details on what you did.

Do you remember to which directory you cloned isoAdvector?
Do you know what your $FOAM_USER_APPBIN and $FOAM_USER_LIBBIN variables are set to?



A general advice: Once you have installed your OpenFOAM version and confirmed that the installation works, you should always remove write access for anyone to the installation. You can for instance do this by going to the installation directory and type:


chmod -R ugo-w


Which means "change permissions (-R)ecursively removing write access (-w) to for User, Group and Others (ugo)".


Then if you at a later time accidentally try to write to the installation directory you will get a warning (note that you can still accidentally remove files with rm -f). If you want to be even safer (and you have permissions), you could transfer ownership of the installation files to the root using chown.
roenby is offline   Reply With Quote

Old   February 27, 2020, 09:36
Default
  #13
New Member
 
Deutschland
Join Date: Jun 2019
Posts: 21
Rep Power: 6
Arghavani is on a distinguished road
thank you for additional explanation.
actually I did it before Chrismas holiday and I only remember that I downloaded the IsoAdvector package and I coyed to my OpenFOAM folder and then since I needed a master version of the OpenFOAM I downloaded that one and in the middle of the installation I got an error and after that, I couldn't run the OpenFoam. I don't remember more details. I have never noticed these things you mentioned (Do you remember to which directory you cloned isoAdvector?
Do you know what your $FOAM_USER_APPBIN and $FOAM_USER_LIBBIN variables are set to?) and I don't have enough knowledge with the installation part. since that time I was busy with another project, I stopped dealing with that problem and now I want to continue it and then it is my problem. I will ask my expert boss to help me. and Thank you again.

Best,
Arghavan
Arghavani is offline   Reply With Quote

Old   May 26, 2020, 22:48
Smile Adding source terms to alpha equation and solve using isoAdvector
  #14
New Member
 
Jobin Raju
Join Date: Nov 2019
Location: Germany
Posts: 3
Rep Power: 6
jobin2600 is on a distinguished road
Hello Everyone

I recently found this excellent method for surface reconstruction. Since I am dealing with mass transfer, I need to add a source term to the alpha equation. I see that MULES take the source terms Su, Sp from alphaSuSp.H . But couldn't figure out if isoAdvector handles those terms as well. I had installed isoAdvector from Github on OpenFOAM 6(and it works perfect).

How can I add the source terms to the alpha equation and make use of isoAdvector to solve the equation ?

Any help is highly appreciated

Best regards

Jobin Raju
jobin2600 is offline   Reply With Quote

Old   May 27, 2020, 06:23
Default
  #15
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 92
Rep Power: 20
roenby will become famous soon enough
Hi Jobin Raju


Adding source terms for compressibility was done by Henning Scheufler here:
https://github.com/DLR-RY/VoFLibrary


In particular, in alphaEqn.H of the compressibleInterFlow solver we now have:


#include "alphaSuSp.H"
advector->advect(Sp,(Su + divU*min(alpha1(), scalar(1)))());

To see how the advect function of isoAdvector handles the source terms see here.


Note that hopefully this work will be integrated in the upcoming OpenFOAM-v2006.


Kind regards,
Johan
roenby is offline   Reply With Quote

Old   September 3, 2020, 23:41
Smile Understanding the theory behind interface reconstruction - reference material
  #16
New Member
 
Jobin Raju
Join Date: Nov 2019
Location: Germany
Posts: 3
Rep Power: 6
jobin2600 is on a distinguished road
Hi Johan

Is there any reference material or publications available from which I can understand the theory behind the interface reconstruction part(iso-method) in the isoAdvector library. I went through your papers and couldn't find the mathematical details of the algorithm used.

Any help is highly appreciated.

Best,
Jobin
jobin2600 is offline   Reply With Quote

Old   September 4, 2020, 06:36
Default
  #17
Member
 
Johan Roenby
Join Date: May 2011
Location: Denmark
Posts: 92
Rep Power: 20
roenby will become famous soon enough
The reconstruction step of isoAdvector is for instance described here (section 3.1):https://doi.org/10.1016/j.jcp.2019.01.009


Or here (also section 3.1):
https://royalsocietypublishing.org/d...98/rsos.160405
jobin2600 likes this.
roenby is offline   Reply With Quote

Reply

Tags
isoadvector

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
Convergence issues for Flat plate with sharp edge rajnarayang FLUENT 3 June 20, 2017 13:02
[ANSYS Meshing] Multizone issues (on my project) crenaudo ANSYS Meshing & Geometry 8 April 13, 2016 03:59
Multigrid Stability Issues ThomasHermann SU2 1 November 5, 2014 17:18
[General] Some Paraview Issues I can not solve MR_Chicho ParaView 1 September 24, 2012 06:03
compressible modification of nearwall turbulence Quain Tchew Main CFD Forum 0 March 4, 2002 02:29


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