Comments (7)
Hello -
Based on your post, I see you have two regions, named fluid
and solid
.
Therefore, please try
postProcess -func wallHeatFlux -region fluid
The polyMesh error is due to your not specifying the region.
Your use of the function in controlDict should also require you specify the region. If memory serves, it is as simple as adding the following to the script wallHeatFlux:
region fluid;
from hybridcentralsolvers.
Thanks for quick reply. Though, this doesn't work. I did try this in the variations I mentioned about.
Edit: If it helps, here's the corresponding error
postProcess -func wallHeatFlux -region fluid
--> FOAM FATAL ERROR:
Unable to find compressible turbulence model in the databaseFrom function virtual bool Foam::functionObjects::wallHeatFlux::execute() in file wallHeatFlux/wallHeatFlux.C at line 237.
FOAM exiting
from hybridcentralsolvers.
Hmm. Well, you can do some programming, or you can just trick OpenFOAM. If you don't mind not being elegant...
cp the contents of fluid
into constant
. That includes the polyMesh
directory plus thermoPhyiscalProperties
and turbulenceProperties
.
cp system/fvSchemes
and system/fvSolution
to system
For the time step of interest, cp the contents of the fluid directory (e.g., p, phi, rho, T, and U) up one "directory level" so that, for example, under time 0.001
you find p, phi, rho, etc. etc.
use rhoPimpleFoam to get the wall heat flux.
e.g.,
rhoPimpleFoam -postProcess -func wallHeatFlux -latestTime
I didn't run the laminar superSonicCorner case to completion... just ran to 0.0066 seconds. Here's the result of the rhoPimpleFoam postProcess:
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : rhoPimpleFoam -postProcess -func wallHeatFlux -latestTime
Date : May 17 2019
Time : 11:55:55
Host : "aris"
PID : 4550
I/O : uncollated
Case : /home/eadaymo/OpenFOAM/eadaymo-6/applications/hybridCentralSolvers/Tutorials/superSonicCone/laminar-redo
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create mesh for time = 0.0066
PIMPLE: No convergence criteria found
PIMPLE: No corrector convergence criteria found
Calclations will do 2 corrections
wallHeatFlux wallHeatFlux:
processing all wall patches
Time = 0.0066
Reading thermophysical properties
Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}
Reading field U
Reading/calculating face flux field phi
Creating turbulence model
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating field dpdt
Creating field kinetic energy K
No MRF models present
No finite volume options present
wallHeatFlux wallHeatFlux write:
writing object wallHeatFlux
min/max/integ(cone_fluid) = -718.076, -6.62934, -2.197
End
Hope this helps... Like I said it's not elegant, but may just get you the result you need.
from hybridcentralsolvers.
Here's the fancier approach.
I added postProcess.H to the chtMultiRegionSolver... I had to do a few more edits to get it to compile. Use at your own risk - I did not test thoroughly.
When you compile the attached you'll get chtMultiRegionCentralFoam-postProcess.
chtMultiRegionCentralFoam-postProcess.tar.gz
Following a tutorial, you can run:
chtMultiRegionCentralFoam-postProcess -postProcess -func wallHeatFlux -region
And get heat fluxes like the following:
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
Build : 6-e29811f5dff8
Exec : chtMultiRegionCentralFoam-postProcess -postProcess -func wallHeatFlux -region fluid
Date : May 17 2019
Time : 12:02:23
Host : "aris"
PID : 4786
I/O : uncollated
Case : /home/eadaymo/OpenFOAM/eadaymo-6/applications/hybridCentralSolvers/Tutorials/superSonicCone/laminar-redo
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create fluid mesh for region fluid for time = 0
Create solid mesh for region solid for time = 0
wallHeatFlux wallHeatFlux:
processing all wall patches
Time = 0
Reading thermophysical properties
Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}
Reading field U
Creating field kinetic energy K
PIMPLE: No convergence criteria found
PIMPLE: No corrector convergence criteria found
Calclations will do 2 corrections
Selecting finite volume kappaFunction type standardMachToAcCourantRatio
Creating turbulence model
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No finite volume options present
No MRF models present
Marking cells with bad quality to limit pressure gradient
max/min FaceAcCourant: 0.72/1.0903e-23
max/min Maf: 5.00866/0
max/min kappa: 1/0
*** Reading solid mesh thermophysical properties for region solid
Adding to thermos
Selecting thermodynamics package
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}
Adding to radiations
Selecting radiationModel opaqueSolid
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel none
Selecting sootModel none
Adding fvOptions
wallHeatFlux wallHeatFlux write:
writing object wallHeatFlux
min/max/integ(cone_fluid) = 0, 0, 0
Time = 0.0002
Reading thermophysical properties
Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleEnthalpy;
}
Reading field U
Creating field kinetic energy K
PIMPLE: No convergence criteria found
PIMPLE: No corrector convergence criteria found
Calclations will do 2 corrections
Selecting finite volume kappaFunction type standardMachToAcCourantRatio
Creating turbulence model
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No finite volume options present
No MRF models present
Marking cells with bad quality to limit pressure gradient
max/min FaceAcCourant: 0.946666/1.43302e-23
max/min Maf: 5.22189/0
max/min kappa: 1/0
*** Reading solid mesh thermophysical properties for region solid
Adding to thermos
Selecting thermodynamics package
{
type heSolidThermo;
mixture pureMixture;
transport constIso;
thermo hConst;
equationOfState rhoConst;
specie specie;
energy sensibleEnthalpy;
}
Adding to radiations
Selecting radiationModel opaqueSolid
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting scatterModel none
Selecting sootModel none
Adding fvOptions
wallHeatFlux wallHeatFlux write:
writing object wallHeatFlux
min/max/integ(cone_fluid) = -2.89064e-12, 3.1529e-12, -3.54824e-15
.
.
.
etc.
from hybridcentralsolvers.
Thanks a lot for this! This works for now.
There is some issue with parallel running, but that's another topic.
from hybridcentralsolvers.
Glad it works. Not sure what to say on the parallel issue, though. Hope this keeps you moving forward. Cheers
from hybridcentralsolvers.
Dear colleagues, i would like to stress, that only OpenFOAM+ fork of OpenFOAM technology is now supported, so please use appropriate branches of repository
from hybridcentralsolvers.
Related Issues (15)
- Is it correct to initialize this here and not before? HOT 1
- Termination in parallel, but not serial calculation HOT 27
- OpenFOAM v4>pisoCentral>tools>correctCentralACMIInterpolation Class HOT 1
- Openfoam-v1706 HOT 2
- pisoCentral OF6 HOT 2
- Some error when running parallel on OpenFOAM-v2212 HOT 9
- Max Mach number? HOT 1
- How to implement QGDFoam's readScNumbers.H to reactingPimpleCentralFoam? HOT 3
- compilation errors for OF2306 HOT 2
- stabilization problem HOT 3
- Issue with wedge domains HOT 5
- Compiling pimpleCentralFoam problem on openfoam v1912 HOT 8
- coupling with adjointShapeOptimizationFoam HOT 7
- Trying with other geometries HOT 4
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from hybridcentralsolvers.