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

How to use the dimensionless N-S equations in OpenFOAM

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 7, 2016, 04:40
Smile How to use the dimensionless N-S equations in OpenFOAM
  #1
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
Dear All,
I am trying to use a set of dimensionless equations in OpenFOAM, as in the post :
How does one turn off the dimensional checking
and in
http://www.cfd-online.com/Forums/ope...quation-2.html
I have run a case with dimension forms without any problem.
Attached please find the new solver I want to use. Everytime I run the case with the new solver, there exists an error as below. I have copied the controlDict out to my working directory and changed the dimensionSet from 1 to 0, but the same error shows up.
\*--------------------------------------------------------------------------------------------------*/
Build : 3.0.1-d8a290b55d28
Exec : myIcoFoam
Date: Dec 10 2016
Time: 17:53:57
Host:
PID:
Case:
nPorcs: 1
sigFpe: Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChcking: Monitoring run-time modified files using timeStampMaster
allowSystemOperation : Allowing user-supplied system call operations

//************************************************** ************************************************//
Create Time

Create mesh for time=0

PISO: Operating solver in PISO mode

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi

Starting time loop

cellSource functionObjects_allvolAverage:All():
total cells =103462
total volume = 0.0400007

Time=0.0125

Courant Number mean:0 max:0

-->FOAM FATAL IO ERROR:
incompatible dimensions for operation
[U[0 0 -1 0 0 0 0]]+ [U[0 -1 0 0 0 0 0]]

From function checkMethod (const fcMatrix<Type>&, const fvMatrix<Type>&)
in file /usr/local/OpenFOAM/OpenFOAM-3.0.1/src/finiteVolume/lnInclude/fcMatrix.C at line 1295.


FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
# Foam::IOerror::abort() at ??:?
#2 void Foam::checkMethod<Foam::Vector<double>>(Foam::fvMa trix<Foam::Vector<double>> const&, Foam::fvMatrix<Foam::Vector<double>> const&, char const*) at ??:?
#3 ? at ??:?
#4 ? at ??:?
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? at ??:?

Aborted

I wonder how to make a change to the dimensionedType.C.
I wonder if any one could give me some help.
Thank you in advance.
Best Regards,
Bill
Attached Files
File Type: zip applications.zip (169.8 KB, 17 views)
File Type: txt transportProperties.txt (1.0 KB, 12 views)
File Type: zip 0.zip (1.1 KB, 9 views)
File Type: zip system.zip (4.7 KB, 7 views)

Last edited by 6863523; December 12, 2016 at 16:13.
6863523 is offline   Reply With Quote

Old   December 7, 2016, 05:29
Default
  #2
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
Hi,
I copied the dimensionedType.C out, as what I have done for controlDict; Even, I deleted the dimension checking syntax, there is still the same error, which means the dimensionedType.C was not changed at all. Or the problem is caused by something else?
Please help, I have stuck for several days.
Thank you in advance.
Best Regards,
Bill
6863523 is offline   Reply With Quote

Old   December 9, 2016, 06:04
Default
  #3
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
Any suggestion would appreciated.
Thank you in advance.
6863523 is offline   Reply With Quote

Old   December 10, 2016, 07:12
Smile to use the dimensionless N-S equations in OpenFOAM
  #4
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
The files have been updated after some modifications.
Thanks.

Last edited by 6863523; December 11, 2016 at 06:52.
6863523 is offline   Reply With Quote

Old   December 10, 2016, 08:41
Smile to use the dimensionless N-S equations in OpenFOAM
  #5
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
Just updated.
If I deleted the 'nu' in my new solver and the 'nu' definition in the transportProperties.
It will show that the nu is undefined in the transportProperties.
The problem seems to be clear, the newly defined solver may not be used at all by the program.
I will define the new solver again to test if it can be used or not.
6863523 is offline   Reply With Quote

Old   December 10, 2016, 11:57
Smile
  #6
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
Newly updated.
I have tried the new solver for other cases, the new solver can be used, showing that the undefined Re, which I have used in the new solver.
The problem seems that the dimension checking was not turned off at all by copying out the controlDict and switching the dimensionSet from 1 to 0.
Since I do not have the privilege to modify this file, I wonder if there exists any other solution, except copying the controlDict out.
Still, any suggestion would be appreicated.
6863523 is offline   Reply With Quote

Old   December 10, 2016, 14:11
Smile to use the dimensionless N-S equations in OpenFOAM
  #7
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
Any information missed, please let me know.
Could any one try to turn off the dimension checking and let me know if it is possible?
Thank you in advance.

Last edited by 6863523; December 11, 2016 at 09:42.
6863523 is offline   Reply With Quote

Old   December 10, 2016, 14:28
Smile to use the dimensionless N-S equations in OpenFOAM
  #8
Senior Member
 
Bill Wang
Join Date: Aug 2014
Posts: 109
Rep Power: 8
6863523 is on a distinguished road
The dimensionless Navier-Stokes equations that I want to solve in the new created solver is:
fvVectorMatrix UEqn
(
Re*(fvm::ddt(U))
+ Re*(fvm::div(phi, U))
- fvm::laplacian(nu,U)
+ F
);

solve(UEqn == -fvc::grad(P));

The '' F '' is a unit vector, standing the gradient of the nondimensionalized mean pressure, as can be seen in the attached equations. I have it like an input (I will study my flow in function of given values of Re and F).

I have decomposed the traditional p into F+ grad(p), as in the attachment.
Thank you. Any information needed, please let me know.
Attached Images
File Type: jpg equations.jpg (13.9 KB, 53 views)
File Type: jpg decomposed p.jpg (1.5 KB, 20 views)

Last edited by 6863523; December 13, 2016 at 10:07.
6863523 is offline   Reply With Quote

Reply

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
Getting Started with OpenFOAM wyldckat OpenFOAM 24 October 2, 2019 22:35
OpenFOAM for Solving Full Transonic Potential Equations jotter OpenFOAM Programming & Development 1 November 25, 2017 19:47
OpenFOAM v3.0+ ?? SBusch OpenFOAM 22 December 26, 2016 15:24
OpenFOAM Training, London, Chicago, Munich, Houston 2016-2017 cfd.direct OpenFOAM Announcements from Other Sources 0 September 14, 2016 04:19
Actual drag force from dimensionless equations slaxmi CFX 10 September 14, 2007 20:20


All times are GMT -4. The time now is 17:29.