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

Search Results

Register Blogs Community New Posts Updated Threads Search

Showing results 1 to 11 of 11
Search took 0.00 seconds.
Search: Posts Made By: bigphil
Forum: OpenFOAM Programming & Development April 9, 2018, 05:13
Replies: 47
Views: 6,186
Posted By bigphil
Without seeing the rest of your code, I cannot...

Without seeing the rest of your code, I cannot say whether you need to explicitly perform a parallel sync; however, from looking at your code snippet, the only terms that might need it are those you...
Forum: OpenFOAM Programming & Development April 7, 2018, 10:38
Replies: 47
Views: 6,186
Posted By bigphil
Actually, as "snGrad()" constructs and returns...

Actually, as "snGrad()" constructs and returns the full patch field, it will be more efficient to store the snGrad as follows:

const polyPatchList& patches = mesh.boundaryMesh();...
Forum: OpenFOAM Programming & Development April 7, 2018, 10:33
Replies: 47
Views: 6,186
Posted By bigphil
Sorry, "forAll(patchI, faceI)" should have been...

Sorry, "forAll(patchI, faceI)" should have been "forAll(gradcMinusField, faceI)": I have editted this in my previous post.

Also, for your second error, please be careful, you have written:

else...
Forum: OpenFOAM Programming & Development April 6, 2018, 04:54
Replies: 47
Views: 6,186
Posted By bigphil
Try the following: const polyPatchList&...

Try the following:


const polyPatchList& patches = mesh.boundaryMesh();
forAll(patches,patchI)
{
if (patches[patchI].name() == "midup")
{
fixedGradientFvPatchScalarField&...
Forum: OpenFOAM Programming & Development June 6, 2017, 12:53
Replies: 47
Views: 6,186
Posted By bigphil
Hi Bobi, The block coupled method will help...

Hi Bobi,

The block coupled method will help the convergence of the cPlus and cMinus equations; however, it seems that currently you are limited by the convergence of the Ue equation (because of...
Forum: OpenFOAM Programming & Development May 27, 2017, 11:08
Replies: 47
Views: 6,186
Posted By bigphil
babakflame, You can manually apply...

babakflame,

You can manually apply under-relaxation when setting the boundary conditions e.g.

// No under-relaxation
//gradcMinusField = cMinus.boundaryField()[patchI] *...
Forum: OpenFOAM Programming & Development May 26, 2017, 13:03
Replies: 47
Views: 6,186
Posted By bigphil
Hi babakflame, Can you also post the output...

Hi babakflame,

Can you also post the output from the solver? What exactly is diverging? the linear solver for one of the equations? or the outer pimple loop? or time-steps?

From a quick look...
Forum: OpenFOAM Programming & Development May 22, 2017, 09:27
Replies: 47
Views: 6,186
Posted By bigphil
Hi babakflame, Can you post the code you are...

Hi babakflame,

Can you post the code you are using?
This might help give an idea of how the stability of your method could be improved (also it will help others who might want to try something...
Forum: OpenFOAM Programming & Development March 8, 2017, 03:35
Replies: 47
Views: 6,186
Posted By bigphil
The snGrad at a patch is correctly calculated...

The snGrad at a patch is correctly calculated from:

Ue.boundaryField()[patchI].snGrad()


If you were to use:

fvc::snGrad(Ue)

this would give you the snGrad on every internal face and...
Forum: OpenFOAM Programming & Development March 7, 2017, 15:47
Replies: 47
Views: 6,186
Posted By bigphil
Hi babakflame, A couple of...

Hi babakflame,

A couple of comments/questions:


in the picture of your domain, the normal vectors "n" point inwards (i.e. in to the domain); however, OpenFOAM convention assumes the normals...
Forum: OpenFOAM Programming & Development March 6, 2017, 03:34
Replies: 47
Views: 6,186
Posted By bigphil
Hi babakflame, What equation are you trying...

Hi babakflame,

What equation are you trying to implement as the boundary condition?

I don't think the following makes sense:

cPlus.boundaryField()[patchI].snGrad() =
...
Showing results 1 to 11 of 11

 
Forum Jump

All times are GMT -4. The time now is 05:14.