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

picture "cracked" on paraFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 29, 2018, 09:06
Default picture "cracked" on paraFoam
  #1
New Member
 
naufal
Join Date: May 2018
Posts: 11
Rep Power: 7
mnaufalazwar is on a distinguished road
Hi guys, I am working on heat transfer in porous media, and I have a problem with my paraview when I run paraFoam on my terminal.

On my case, the flow is driven by pressure gradient, when I put low pressure on my case, it seems fine, but when I put higher pressure, the picture in paraview seems broken like this :


I think my problem is because the temperature is very high and it makes my flow is turbulent, how can I insert the turbulent characteristic on my solver?

I make my own solver named darcyTemperatureFoam, here is my createFields.H :

Info<< "Reading field p\n" << endl;

volScalarField p
(
IOobject
(
"p",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Info<< "Reading field U\n" << endl;

volVectorField U
(
IOobject
(
"U",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedVector("U",dimensionSet(0,1,-1,0,0,0,0),vector::zero)
);

surfaceScalarField phi ("phi", linearInterpolate(U) & mesh.Sf());

Info<< "Reading field T\n" << endl;

volScalarField T
(
IOobject
(
"T",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Info<< "Reading transportProperties\n" << endl;

IOdictionary transportProperties
(
IOobject
(
"transportProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ_IF_MODIFIED,
IOobject::NO_WRITE
)
);

Info<< "Reading diffusivity DT\n" << endl;

dimensionedScalar DT
(
transportProperties.lookup("DT")
);

Info<< "Reading porosity eps\n" << endl;

dimensionedScalar eps
(
transportProperties.lookup("eps")
);

Info<< "Reading solidHeatCapacities rhoCps\n" << endl;

dimensionedScalar rhoCps
(
transportProperties.lookup("rhoCps")
);

Info<< "Reading fluidHeatCapacities rhoCpf\n" << endl;

dimensionedScalar rhoCpf
(
transportProperties.lookup("rhoCpf")
);

Info<< "Reading permeability k\n" << endl;

dimensionedScalar k
(
transportProperties.lookup("k")
);

Info<< "Reading fluid viscosity mu\n" << endl;

dimensionedScalar mu
(
transportProperties.lookup("mu")
);

and here is my darcyTemperatureFoam.C

#include "fvCFD.H"
#include "simpleControl.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])
{
#include "setRootCase.H"

#include "createTime.H"
#include "createMesh.H"

simpleControl simple(mesh);

#include "createFields.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nCalculating pressure distribution\n" << endl;

while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
solve
(
fvm::laplacian(k/mu, p)
);
}

U = -k/mu*fvc::grad(p);

phi = linearInterpolate(U) & mesh.Sf();

solve
(
(eps*rhoCpf+(1.-eps)*rhoCps)*fvm::ddt(T)
+ rhoCpf*fvm::div(phi, T)
==
fvm::laplacian(DT, T)
);

runTime.write();

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}

Info<< "End\n" << endl;

return 0;
}
Attached Images
File Type: jpg broken picture.JPG (11.8 KB, 38 views)
mnaufalazwar 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
[OpenFOAM] Unable to launch parafoam & boonhow ParaView 5 March 23, 2022 09:09
[OpenFOAM.com] paraFoam cannot open due to Qt issues - [Solved/Information] u2berggeist OpenFOAM Installation 2 July 2, 2018 18:03
[OpenFOAM] Visualization problem on ParaFoam Rider ParaView 10 June 27, 2016 11:19
[OpenFOAM] How to use paraFoam on a cluster andreas ParaView 1 March 6, 2013 18:11
[OpenFOAM] OpenFoam (Ubuntu): paraFoam via Xming+PuTTY raketenmaid ParaView 4 February 5, 2013 06:20


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