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

How to remove viscosity term from the governing equation in interFoam solver

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 17, 2022, 15:56
Default How to remove viscosity term from the governing equation in interFoam solver
  #1
New Member
 
Basar H.Ibrahim
Join Date: Sep 2021
Posts: 7
Rep Power: 5
Pasar is on a distinguished road
Hello All

After I modified interFoam solver to newInterFoam solver,
I want to remove the viscosity term from the governing equation,
And I know that the library that controls the transport properties is src/transportModels/immiscibleIncompressibleTwoPhaseMixture,
And this library includes two other libraries which are

#include "incompressibleTwoPhaseMixture.H"
#include "interfaceProperties.H"

So, I modified these libraries and copied them into my user directory as well,
and I made wmake and my new solver works very well.

I am trying to modify these libraries but I don't know which part should I remove

My question
1- is this way is true to remove the viscosity term?
2- which libraries should I modify?

By the way, I have made a lot of tests by modifying the libraries but I still got errors.

Best Regards
Pasar is offline   Reply With Quote

Old   February 18, 2022, 06:23
Default
  #2
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 798
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
1/ To set viscosity to zero, it suffices to set the molecular diffusion constant equal zero in the constant/transportProperties file;

2/ It is likely that the solver fails to converge in the absence of diffusion;

Best, Domenico.
dlahaye is offline   Reply With Quote

Old   February 18, 2022, 06:38
Default
  #3
New Member
 
Basar H.Ibrahim
Join Date: Sep 2021
Posts: 7
Rep Power: 5
Pasar is on a distinguished road
Quote:
Originally Posted by dlahaye View Post
1/ To set viscosity to zero, it suffices to set the molecular diffusion constant equal zero in the constant/transportProperties file;

2/ It is likely that the solver fails to converge in the absence of diffusion;

Best, Domenico.
I did these steps, but the solver didn't work.
If the solver works correctly it is not what I want
What I want to do is to remove from the source code.
In general I want to make new solver with out viscosity.

thankyou for your reply dlahaye
Best regards.
Pasar is offline   Reply With Quote

Old   February 18, 2022, 06:56
Default
  #4
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 798
Blog Entries: 1
Rep Power: 17
dlahaye is on a distinguished road
Removing the diffusion term from the source code is a very humble endaveour.

I seriously doubt it will converge. Diffusion is often required to obtain converge. Stated differently, to obtain convergence is the absence of diffusion, other discretization techniques and solvers might be required.

The book by LeVeque https://www.cambridge.org/core/books...D52EAD6909E2B9 among many others might be a good reference.
dlahaye is offline   Reply With Quote

Reply

Tags
remove, solver development, viscosity.


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
Error SIGSEGV using VOF and UDF JERC_UTFSM Fluent UDF and Scheme Programming 14 November 8, 2021 00:17
Question about interFoam Solver Kahnbein.Kai OpenFOAM Running, Solving & CFD 2 August 26, 2019 16:36
How can use MacCormack & Baldwin artificial viscosity in term of generalized nonortho Nam Tran Ha CFD Freelancers 0 December 21, 2016 04:31
how to troduce a source term in ICOFOAM solver's equation sawyer86 OpenFOAM Running, Solving & CFD 0 July 20, 2012 12:15
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45


All times are GMT -4. The time now is 06:48.