New element lops time off CFD simulations

March 9, 2006
"CFD simulations are finite-volume based, so the solver calculates mass, momentum and energy fluxes through the surface of these control volumes," says CD-adapco Vice President Dennis Nagy.

The test case has a laminar incompressible flow in a 100-m pipe with a bend. The mesh-dependence study focuses on a static pressure drop from end to end. Inlet flow is at 10 m/sec and outlet at 0 pressure. The computation on any mesh is considered converged when the average inlet-pressure change is less then 0.01 Pa for 10 iterations.


The same tube is meshed with tets of similar surface dimensions.


Results for tets and polyhedrons show similar results despite a widely different number of elements.


The last time finite-element meshing was in the headlines, developers were debating the advantages of tetrahedral (pyramid) versus hexahedral (brick) elements. Simulation experts today are talking about polyhedral elements, those with 12 to 14 faces and shaped like soccer balls. The unusual element comes from CD-adapco Inc., Melville, N.Y., and London, England, and is said to work well on large CFD problems, but also shows promise for other simulations such as heat transfer/conduction and stresses in solids.

"And the fewer element surfaces in the grid, the shorter the solution time."

Tets and hexes were used because meshers could build grids of them and compute shape functions in the volumes. In the 1990s, researchers recognized it would be useful to extend meshing to more than six faces, the number on a brick or hex element. A polyhedral mesh is a good candidate for the next preferred element (or in the case of CFD, control volume) because it can be formed by combining several tets into a single polyhedron and is thus more flexible for filling the space with well-shaped elements.

Low-order elements, such as tets, are often deformed during meshing to look more like wedges or slivers. "Some of these point upstream, letting fluid flows hit their very oblique surfaces. The disparity between the small inflow area and large outflow areas leads to excessive numerical diffusion. Ideally, flow should be perpendicular to the element surfaces. Hence, tetrahedrons are the least ideal computational cell. And they are terribly inefficient at filling space. More tets are needed to fill a void than higher order elements with similar dimensions. Hexahedrons are an improvement, but they are difficult and expensive to generate automatically," says Nagy.

Polyhedral meshes take up to 30% more time to generate than tet meshes, but they fill space more efficiently, with fewer surfaces of a given size than other elements would require. "A 14-sided polyhedral element, for example, always has a couple faces close to the flow normal. Hence, the solver finds accurate answers and sooner," says Nagy. "Tests show they produce the same answers as traditional meshes but in only 20% of the solution time and with half the memory required."

So-called face-addressing CFD solvers benefit from their ability to run polyhedral meshes, which contain fewer faces in a computational grid, for a given mesh resolution. Minimizing the number of faces boosts solver speed without sacrificing accuracy. Because each polyhedral cell has more faces than a tet or hex element, it also has more neighbors than traditional cell

types. The result, says Nagy, is that data propagates more quickly through polyhedral meshes, which leads to higher rates of convergence.

A Tale Of Tets And Polyhedrons

Surface size (m)
Tets
Polys
4.0
95,050
19,288
2.00
94,542
19,187
1.00
191,873
36,586
0.50
609,565
112,471
0.25
3,351,276
597,342
Surface size refers to the height of a tet or diameter of a polyhedron.

Of bubbles and polyhedrons

The graph shows that even a relatively coarse polyhedron mesh which requires 1.6 hr to converge approaches a grid-independent solution. A finer mesh that needs 6 hr to converge provides a solution with an acceptably small error. With tetrahedral meshes, results overshoot the grid-independent solution. Only a very fine mesh that needs over 33 hr to converge provides a solution with the same acceptably small error.


In 1887, Lord Kelvin found that a 14-sided polyhedron (tetrakaidecahedron) would most efficiently pack bubbles into foam when all bubbles are of equal size. In 1994, physicists Denis Weaire and Robert Phelan found that a mixture of 12 and 14-sided polyhedrons partition space 3% more efficiently than Kelvin's foam.

The recent polyhedral meshes actually use an unrestricted number of faces, so they fill space in most efficiently. In fact, for a given resolution level, a mesh of polyhedral cells has fewer faces than a mesh of any other cell type.

CD-adapco Inc., cd-adepco.com

Sponsored Recommendations

Altech's Liquid Tight Strain Relifs Catalog

March 13, 2024
With experienced Product Engineers and Customer Service personnel, Altech provides solutions to your most pressing application challenges. All with one thought in mind - to ensure...

Industrial Straight-Through Cable Gland

March 13, 2024
Learn more about Altech's cable glands and all they have to offer for your needs!

Smooth Sorting with SEW-EURODRIVE!

Feb. 22, 2024
Sorting systems are essential when it comes to warehouse automation, material handling, and distribution. SEW-EURODRIVE’s automated sorting solutions increase capacity, reliability...

Setting the Agricultural Industry in Motion

Feb. 22, 2024
Agricultural products go through specialized processes of harvesting, storing, conveying, cleaning, drying, sorting, and handling to ensure that seasonal products reach the processing...

Voice your opinion!

To join the conversation, and become an exclusive member of Machine Design, create an account today!