[mesh manipulation] refine mesh based on skewness

amuzeshi · July 16, 2019, 07:01

Hi

I'm on simulation of an elastic square prism. Results are obtained to some time, then the solver (solids4Foam) crashes. The reason is that the highlighted cell(s) shown in the figure below becomes concave. (Note: Fluid region mesh is shown only.)
Is there any way for run-time mesh refinement according to skewness in order to figure out this difficulty?

regards
Ali

amuzeshi · July 16, 2019, 07:46

I know that the two classes dynamicMotionSolverFvMesh and dynamicRefineFvMesh can be combined into one class[1]. However, as it is needed to compile a new solver and initialize the skewness field in its *.H file, there remain two main problems:

1- As solids4Foam does not obey solver-oriented programming discipline in regular openfoam, where to declare skewness?

2- How to declare skewness as -to the best of my knowledge- it is not a volScalarField.

regards
Ali

Ref:
[1] Adaptive Mesh Refinement with a Moving Mesh using sprayDyMFoam

amuzeshi · July 17, 2019, 15:04

I know that OpenFOAM defines skewness for faces, not for cells. But I want to refine the mesh based on a user-defined skewness:

cell skewness = max [ ( thetaMax - 90) / 90 , ( 90 - thetaMin) / 90]

ancolli · January 29, 2024, 06:15

Hi, have you been able to solve your problem? Is there anything you'd like to share?
Kind regards

amuzeshi · January 29, 2024, 12:17

Quote:

Originally Posted by ancolli

Hi, have you been able to solve your problem? Is there anything you'd like to share?
Kind regards

Hi
For a similar problem, I use manual re-meshing, so, when the simulation crashes, I start off the latest available time using a new mesh.

July 16, 2019, 07:46		#2
amuzeshi Senior Member Ali Shayegh Join Date: Oct 2015 Posts: 130 Rep Power: 10	I know that the two classes dynamicMotionSolverFvMesh and dynamicRefineFvMesh can be combined into one class[1]. However, as it is needed to compile a new solver and initialize the skewness field in its .H file, there remain two main problems: 1- As solids4Foam does not obey solver-oriented programming discipline in regular openfoam, where to declare skewness? 2- How to declare skewness as -to the best of my knowledge- it is not a volScalarField. regards Ali Ref: [1] Adaptive Mesh Refinement with a Moving Mesh using sprayDyMFoam Last edited by wyldckat; August 3, 2019 at 21:53. Reason: fixed link*

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[snappyHexMesh] Add Mesh Layers doesnt work on the whole surface	Kryo	OpenFOAM Meshing & Mesh Conversion	13	February 17, 2022 07:34
SimpleFoam & Theater	jipai	OpenFOAM Running, Solving & CFD	3	June 18, 2019 10:11
[snappyHexMesh] snappyHexMesh does not create any mesh except one for the reference cell	Arman_N	OpenFOAM Meshing & Mesh Conversion	1	May 20, 2019 17:16
Coarsen and refine mesh in Pointwise	raye_xie	Pointwise & Gridgen	3	June 26, 2017 17:35
Skewness , Minimum Edge Length and Mesh Quality in FLuent	ViLaks	FLUENT	2	March 2, 2017 05:57

July 17, 2019, 15:04		#3
amuzeshi Senior Member Ali Shayegh Join Date: Oct 2015 Posts: 130 Rep Power: 10	I know that OpenFOAM defines skewness for faces, not for cells. But I want to refine the mesh based on a user-defined skewness: cell skewness = max [ ( thetaMax - 90) / 90 , ( 90 - thetaMin) / 90]

January 29, 2024, 06:15		#4
ancolli Senior Member Alejandro Join Date: Jan 2014 Location: Argentina Posts: 128 Rep Power: 12	Hi, have you been able to solve your problem? Is there anything you'd like to share? Kind regards