Solving a surfaceScalarfield - how to implement boundary?

stockzahn · March 9, 2020, 13:23

Dear FOAMers,

I would like to solve an energy equation knowing all the respective values on the cells faces, like examplarily shown is this small piece of code:

Code:


fvm::(rho, h)
==
fvc::div(kappaS * snGradT * faceArea)

kappaS ... surfaceScalarField containing the thermal conductivity values
snGradT ... surfaceScalarField containing the surface normal temperature gradients
snGradT ... surfaceScalarField containing face Areas

This code works fine for the internal field, e.g. with adiabatic boundaries and an internal heat source. Unfortunately, I'm not able to integrate the boundary conditions in an elegant way; until now I calculate and add an additional source term, whose calculation is different depending on the type of BC. Does anybody know if there is another way to integrate the values of the patch faces into the LES?

Have a nice evening/morning/night/lunch,

stockzahn

stockzahn · March 11, 2020, 11:01

Solved. The boundary fields have to be updated seperately. Problem was my lack of knowledge of the mesh structure in OF.

March 9, 2020, 13:23	Solving a surfaceScalarfield - how to implement boundary?	#1
stockzahn Member Join Date: Oct 2018 Location: France Posts: 34 Rep Power: 8	Dear FOAMers, I would like to solve an energy equation knowing all the respective values on the cells faces, like examplarily shown is this small piece of code: Code: fvm::(rho, h) == fvc::div(kappaS * snGradT * faceArea) kappaS ... surfaceScalarField containing the thermal conductivity values snGradT ... surfaceScalarField containing the surface normal temperature gradients snGradT ... surfaceScalarField containing face Areas This code works fine for the internal field, e.g. with adiabatic boundaries and an internal heat source. Unfortunately, I'm not able to integrate the boundary conditions in an elegant way; until now I calculate and add an additional source term, whose calculation is different depending on the type of BC. Does anybody know if there is another way to integrate the values of the patch faces into the LES? Have a nice evening/morning/night/lunch, stockzahn

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
problem with Min/max rho	tH3f0rC3	OpenFOAM	8	July 31, 2019 10:48
conjugate heat transfer in OpenFOAM	skuznet	OpenFOAM Running, Solving & CFD	99	March 16, 2017 06:07
Moving mesh	Niklas Wikstrom (Wikstrom)	OpenFOAM Running, Solving & CFD	122	June 15, 2014 07:20
pisoFoam with k-epsilon turb blows up - Some questions	Heroic	OpenFOAM Running, Solving & CFD	26	December 17, 2012 04:34
Error while running rhoPisoFoam..	nileshjrane	OpenFOAM Running, Solving & CFD	8	August 26, 2010 13:50

March 11, 2020, 11:01		#2
stockzahn Member Join Date: Oct 2018 Location: France Posts: 34 Rep Power: 8	Solved. The boundary fields have to be updated seperately. Problem was my lack of knowledge of the mesh structure in OF.