volume mesh generation in pHYP #86
Replies: 13 comments 6 replies
-
|
Using |
Beta Was this translation helpful? Give feedback.
-
|
Thank you for the response .If I'm reducing nConstantStart,I'm getting negative volume cells and then segmentation fault error. |
Beta Was this translation helpful? Give feedback.
-
|
I recommend visualizing the mesh to see where the cells are going negative. Then change the options dictionary accordingly. |
Beta Was this translation helpful? Give feedback.
-
|
When I'm getting segmentation fault error ,I'm unable to visualize it. |
Beta Was this translation helpful? Give feedback.
-
|
For future reference, instead of an image, please post the full terminal output as text. |
Beta Was this translation helpful? Give feedback.
-
|
Yes,I'm sure that all surface normals are pointing out This is the terminal output for reference ,I want to know why I'm getting NaN ,what it means? Even,I'm unable to visualize the output file |
Beta Was this translation helpful? Give feedback.
-
|
You can use the |
Beta Was this translation helpful? Give feedback.
-
|
yes,I'm able to visualize,Now how can I eliminate negative volumes and NaN
…On Thu, 18 Apr, 2024, 03:22 Sabet Seraj, ***@***.***> wrote:
You can use the nTruncate option to stop the extrusion before the NaNs
occur but after the negative volumes occur (for example, nTruncate: 83).
This will allow you to visualize the negative volumes.
—
Reply to this email directly, view it on GitHub
<#85 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/BCIUJYVUSWJUMI2WNFVVV23Y53VIPAVCNFSM6AAAAABGJ46TKOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANRSGQ3TQOJZGY>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
-
|
Please attach a picture of the negative volumes |
Beta Was this translation helpful? Give feedback.
-
|
Usually I check negative volumes in pointwise,But I'm unable to open the output file in pointwise ,I'm attaching paraview results And I updated few options in pyhyp script and I added both script and terminal output below rst Importsfrom pyhyp import pyHyp rst generaloptions = { } rst run pyHyphyp = pyHyp(options=options) terminal output Grid | CPU | Sub | KSP | nAvg | Sl | Sensor | Sensor | Min | Min | deltaS | March | cMax | Ratio |Lvl | Time | Its | Its | | | Max | Min | Quality | Volume | | Distance | | kMax |#-------------------------------------------------------------------------------------------------------------------# Grid | CPU | Sub | KSP | nAvg | Sl | Sensor | Sensor | Min | Min | deltaS | March | cMax | Ratio |Lvl | Time | Its | Its | | | Max | Min | Quality | Volume | | Distance | | kMax |#-------------------------------------------------------------------------------------------------------------------# Grid | CPU | Sub | KSP | nAvg | Sl | Sensor | Sensor | Min | Min | deltaS | March | cMax | Ratio |Lvl | Time | Its | Its | | | Max | Min | Quality | Volume | | Distance | | kMax |#-------------------------------------------------------------------------------------------------------------------# |
Beta Was this translation helpful? Give feedback.
-
|
Can anyone help me with this?Why I'm getting Negative Volumes and NaN's
…On Fri, 19 Apr, 2024, 02:02 Galen Ng, ***@***.***> wrote:
Send a picture of the negative volumes
—
Reply to this email directly, view it on GitHub
<#85 (comment)>, or
unsubscribe
<https://github.com/notifications/unsubscribe-auth/BCIUJYQKIE7UQTC2AOV5U7DY6AUWNAVCNFSM6AAAAABGJ46TKOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANRVGI3DOMJYGE>
.
You are receiving this because you authored the thread.Message ID:
***@***.***>
|
Beta Was this translation helpful? Give feedback.
-
|
Your mesh is probably folding in on itself because the wing body junction is concave. Check the options in the documentation for concave geometries. |
Beta Was this translation helpful? Give feedback.
-
|
We're sorry for your frustration. Please report how you modified the smoothing parameters in a systematic way to help us help you debug this problem. |
Beta Was this translation helpful? Give feedback.
-
|
Also wing-body junctions are difficult to do with a single surface mesh extrusion because of the topology. Look into overset meshing. |
Beta Was this translation helpful? Give feedback.
-
|
I modified smoothing parameters as mentioned in the pyhyp options documentation. |
Beta Was this translation helpful? Give feedback.
-
|
If I do overset meshing ,should i keep overset boundary condition at the junction of body and wing? |
Beta Was this translation helpful? Give feedback.
-
|
Thank you so much for the help,I'm able to generate Volume mesh by overset approach |
Beta Was this translation helpful? Give feedback.
-
|
I generated Volume mesh of body+wing geometry by overset approach with the following code ,but while solving the Optimization problem in MACH-Aero ,I'm getting this error,""Trying to create the zipper mesh,but "wall" is not a family in the cgns file or has not been added a combination of families"". Code to generate Volume mesh in pyhyp from collections import OrderedDict parser = argparse.ArgumentParser() Near-Fieldreference first off wall spacing for L2 level meshes#s0 = 1.4e-7 number of Levels in the near-FieldnNearfield = {"L3": 31, "L2": 61, "L1": 121}[args.level] Farfieldbackground mesh spacingdhStar = {"L3": 0.178, "L2": 0.09, "L1": 0.045}[args.level] nFarfield = {"L3": 13, "L2": 25, "L1": 49}[args.level] Generalfactor for spacingsfact = {"L3": 1.0, "L2": 2.0, "L1": 4.0}[args.level] levels of coarsening for the surface meshescoarsen = {"L1": 1, "L2": 2, "L3": 3}[args.level] commonOptions = { body optionsbody_dict = { wing optionswing_dict = { figure out what grids we will generate againoptions = OrderedDict() Run pyHypMultihyp = pyHypMulti(options=options, commonOptions=commonOptions) read the gridsbody = "%s/near_body_vol_%s.cgns" % (args.output_dir, args.level) bodyGrid = readGrid(body) gridList = [bodyGrid, wingGrid] combine gridscombinedGrid = combineGrids(gridList) move to z=0combinedGrid.symmZero("z") Write nearfield meshnearfield = "%s/near_%s.cgns" % (args.output_dir, args.level) MPI.COMM_WORLD.barrier() farfield = "%s/far_%s.cgns" % (args.output_dir, args.level) we can do the stuff in one proc after this pointif rank == 0: what should i change in the volume mesh script to avoid the errror?? Thanks in advance |
Beta Was this translation helpful? Give feedback.
-
I'm trying to generate volume mesh for wing+body using pyHYP .I'm able to run pyhyp script without errors,but the grid is only generating for some distance instead of given march distance. this is the script
what changes can avoid this.Thanks in advance
Beta Was this translation helpful? Give feedback.
All reactions