|
[Sponsors] |
December 7, 2013, 10:23 |
Cell minimum side length
|
#1 |
New Member
M. Sabouri
Join Date: Nov 2011
Posts: 24
Rep Power: 15 |
Dear Foamers,
Does anybody know a simple way to determine the minimum edge length for a single cell an for the whole mesh? Thanks. |
|
December 7, 2013, 11:53 |
|
#2 |
Senior Member
Niels Gjoel Jacobsen
Join Date: Mar 2009
Location: Copenhagen, Denmark
Posts: 1,903
Rep Power: 37 |
Hello,
Yes, you can do something along these lines: Code:
// Get needed reference to the mesh const edgeList& edges = mesh.edges(); const pointField& pp = mesh.points(); const labelListList& cEdges = mesh.cellEdges(); // Make a list with all the edge lenghts scalarField eLengths(edges.size(), 0.0); forAll (edges, edgei) { eLengths[edgei] = edges[edgei].mag(pp); } // Find the minimum edge length in the whole mesh scalar minLength = Foam::min(eLengths); // Find the minimum per cell scalarField minLengthCells(cEdges.size(), GREAT); forAll (cEdges, celli) { const labelList& ce = cEdges[celli]; forAll (ce, edgei) { minLengthCells[celli] = Foam::min(minLengthCells[celli], eLengths[ce[edgei]]); } } Kind regards, Niels P.S. Please note that I have not tried to compile it, so there might be minor blops, but the overall idea gives the requested result.
__________________
Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request. |
|
December 7, 2013, 12:31 |
|
#3 |
New Member
M. Sabouri
Join Date: Nov 2011
Posts: 24
Rep Power: 15 |
Thank you Niels!
I will try it. Moslem |
|
February 23, 2015, 06:30 |
|
#4 |
Member
Naresh Yathuru
Join Date: Feb 2015
Posts: 66
Rep Power: 11 |
hi excuse me for restarting the thread after a long time. i m new to openfoam can some one tell me where to type this above code? i mean in which directory? or in controldict?
Thanks |
|
February 23, 2015, 14:16 |
|
#5 |
New Member
M. Sabouri
Join Date: Nov 2011
Posts: 24
Rep Power: 15 |
It must be added to a source code.
|
|
February 23, 2015, 14:41 |
|
#6 |
Senior Member
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 740
Rep Power: 28 |
This code would need to be implemented in a utility or solver, then compiled and run on the desired case. You will need to add some kind of output capability for your results, either as a logFile or some kind of field.
|
|
Tags |
mesh and grid |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Cell Reynolds Number | laliong | Main CFD Forum | 12 | September 17, 2019 04:18 |
Turbulence length scale and integral length scale | rizhang | CFX | 2 | April 22, 2016 08:22 |
Error finding variable "THERMX" | sunilpatil | CFX | 8 | April 26, 2013 08:00 |
code for designing a minimum length shock free nozzle | staleyjj1981 | Main CFD Forum | 0 | December 4, 2009 16:53 |
How to find the length of cell | zhou1 | FLUENT | 0 | October 26, 2003 13:42 |