CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (http://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   dsmcFieldsCalc not calculating (http://www.cfd-online.com/Forums/openfoam-pre-processing/109314-dsmcfieldscalc-not-calculating.html)

yellowCamel November 14, 2012 17:23

dsmcFieldsCalc not calculating
 
Hello, I am using dsmcFoam to simulate a cold gas thruster that vents to a near vaccum. I used the supersonic wedge as my basis and adjusted the initial conditions for my problem. The geometry is just three concentric cylinders of growing diameter. The first cylinder having a inlet on its face and the final cylinder having an outlet on all faces except the area that interests with the middle cylinder. I am able to intialize and run the program, but when it comes to running the dsmcFieldsCalc, it breaks down. This is the error message that I recieve when I try running the command: Small value (min(mag(rhoNMean)) [0 -3 0 0 0 0 0] 9.99e-301) found in rhoNMean field. Not calculating dsmcFields to avoid division by zero. Has anyone come across this problem before and can explain how they fixed this? Thanks, yc

vinz January 15, 2013 04:42

I can't help you much since I have the same problem.
Did you find a solution?

T-1 February 14, 2013 07:14

Yellow / vinz this happens in the very low density regions. (In paraview you can identify these cells using threshold whilst viewing the density fields.)
I had a similar problem and used some scripting as a quick fix. If you havent resolved the issue since your last message, post your case and i'll check if i can help.

vinz February 14, 2013 09:30

Hello T-1,

The problem looks to be related to the fact that we are simulating wedges which apparently not handled by dsmcFoam so far.

murali May 1, 2013 15:47

dsmcfield not calculating
 
hello T-1

i am also facing same problem, Small value (min(mag(rhoNMean)) [0 -3 0 0 0 0 0] 9.99e-301) found in rhoNMean field. Not calculating dsmcFields to avoid division by zero. calculating averages


can you help me to solve this error

CFDnewbie147 February 11, 2014 05:20

same problem :/
 
Hello, does anybod fix this problem and have a solution?
Please answer if you know it.

Best regards
CFDNewbie147

wyldckat February 16, 2014 11:40

Greetings CFDNewbie147,

I'm sorry to say that everyone who asked about this on this thread, did not follow the basic instructions on http://www.cfd-online.com/Forums/ope...-get-help.html

The very most basic question is: what were the exact steps that were taken to reach that message?

Beyond this:
  1. Which OpenFOAM version is being used?
  2. Which tutorial are you using?
  3. If not a tutorial, on which tutorial did you base yourself on?
  4. Can you share a simple test case in which this can be reproduced?
edit: The following thread seems to be a long discussion on a similar topic: http://www.cfd-online.com/Forums/ope...smcfoam-5.html

Best regards,
Bruno

CFDnewbie147 February 17, 2014 02:49

Hello Bruno,

thank you for answering.

1. I use the OpenFOAM version 2.2.1
2./3. I use the the wedge15Ma5 dsmc tutorial as base for my simulations
I took the simulation case and adapted it to an extern built mesh which I've imported to OpenFOAM without any failures in checkMesh. It looks like a windtunnel with inlet, outlet, walls and the geometry inside the windtunnel. As in the wedge15Ma5 tutorial i made the sides of the windtunnel to the "flow" region with "patch" and the geometry itself as the "obstacle".

With this i do the dsmcInitialise where particles are inserted and then I do dsmcFoam. But when doing the dsmcFieldsCalc there's always the mentioned failure message.

4. It's "only" a windtunner with an object inside it. So yes I can or anybody can built it quickly. (e.g from motorbike tutorial)

I hope this will be enough information.

Best regards
CFDNewbie147

wyldckat February 22, 2014 11:30

Hi CFDnewbie147,

I'm not very familiar with these solvers. But I've found the file "discreteMethods/dsmcFoam/README" in the "tutorials" folder, which has some crucial information on what fields are used by this solver.

Beyond this, please provide a simple case that reproduces this problem. That way it'll take me considerably little time to try and diagnose the problem. Otherwise, don't expect me to give any answers on this before July 2014 :D

Best regards,
Bruno

CFDnewbie147 February 27, 2014 03:11

Hello again,

I think I solved the problem. When you look at the file "discreteMethods/dsmcFoam/README" as you've said: [...]"setting zeroGradient boundary conditions on walls" [...] for iDof, internalE, linearKe, momentum, rhoM and rhoN.
I changed the originally "calculated" for my wall to zeroGradient and the dsmcFields-calculation works and i get overallT, UMean and so on...

I don't know exactly why there was calculated because these parameters are only relevant for cell data but now it works. If anybody does have some hints, please answer :)

Best regards
CFDNewbie148

wyldckat March 1, 2014 06:26

Hi CFDnewbie147,

Quote:

Originally Posted by CFDnewbie147 (Post 477034)
I changed the originally "calculated" for my wall to zeroGradient and the dsmcFields-calculation works and i get overallT, UMean and so on...

I don't know exactly why there was calculated because these parameters are only relevant for cell data but now it works. If anybody does have some hints, please answer :)

It's somewhat simple to explain:
  • "zeroGradient" means that the patch imposes that the derivative must be zero, therefore there is a constraint in the boundary that will make the field values in cells adjacent to the patch to bend the flow in order to make the first derivative equal to zero. Note: this is usually necessary in situations where the equations need a constraint on a boundary, in order to be able to find a solution.
  • "calculated" means that anything that goes through the patch is inferred from the adjacent cells, without further restrictions. As you've mentioned, this is usually used in situations where either only the cell data is relevant or if the fields are fully calculated based on other fields.
Best regards,
Bruno


All times are GMT -4. The time now is 23:11.