{"id":42891,"date":"2020-04-17T17:43:27","date_gmt":"2020-04-17T17:43:27","guid":{"rendered":"\/forum\/forums\/topic\/diverging-solution-from-pressure-based-to-density-based\/"},"modified":"2020-04-17T17:43:27","modified_gmt":"2020-04-17T17:43:27","slug":"diverging-solution-from-pressure-based-to-density-based","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/diverging-solution-from-pressure-based-to-density-based\/","title":{"rendered":"Diverging solution from Pressure-based to Density-based"},"content":{"rendered":"

Hello everyone,<\/p>\n

<\/p>\n

I try to simulate a convergent-divergent nozzle flow based on this subject https:\/\/ntrs.nasa.gov\/search.jsp?R=19820006179<\/a>
Following the method used by these people: https:\/\/tfaws.nasa.gov\/TFAWS11\/Proceedings\/TFAWS2011-AE-001.pdf?fbclid=IwAR28cKI5Ghk_hSVQ5ckgZiXs3ibXEyAWBguCsy7cCUrd8aXHURRLa3gRX1E<\/a><\/p>\n

<\/p>\n

The methods used are:<\/span><\/p>\n

<\/p>\n

– Pressure Based Coupled Solver (PBCS) with 1\/ 2nd order for all equations OR 2\/ PRESTO for Pressure and QUICK for other equations<\/p>\n

<\/p>\n

– Density Based Solver (DBNS) with 2nd order for all equations<\/p>\n

<\/p>\n

Input parameters:<\/span><\/p>\n

<\/p>\n

Material:<\/strong> Air (ideal-gas)<\/p>\n

<\/p>\n

Model:<\/strong> SST k-omega (2 equ.)<\/p>\n

<\/p>\n

Boundary conditions:<\/strong> 1 pressure inlet (2.5 atm), 1 pressure outlet (1 atm)<\/p>\n

<\/p>\n

Solution method:<\/strong> ROE-FDS – Least square – 2nd order<\/p>\n

<\/p>\n

Initialization method:<\/strong> Hybrid or Standard<\/p>\n

<\/p>\n

I successfully represented the model with the Pressure-based solver both for 2nd order and PRESTO\/QUICK as you can see below:<\/p>\n

<\/p>\n

\"\" \"\"<\/p>\n

<\/p>\n

But as soon as I move to the Density based solver, I get a diverging solution, even if I change the Hybrid initialization to the Standard initialization. It always starts pretty well but diverges after around 250 iterations, when the mach disk appears. As you can see on these pictures (I stopped the simulation right before the divergence):<\/p>\n

<\/p>\n

\"\" \"\"<\/p>\n

<\/p>\n

I tried to change settings such as turbulent model (going from k-omega to k-epsilon etc.) and also inlet pressure, pseudo transient on\/off… but it always diverges. I have the same problem on other projects but here, I don’t understand why I have good results with Pressure based solver and not with Density based solver.<\/p>\n

<\/p>\n

At first, I thought it was related to Mesh quality but I tried on an other project (with good results on density based solver) to move from a good mesh (with a converged solution) to coarse mesh, and the solution still converges. I think the problem is related to the shocks, but I can’t verify it.<\/p>\n

<\/p>\n

Have you any thoughts about from where does the problem comes from?<\/p>\n

<\/p>\n

Best regards,<\/p>\n

<\/p>\n

Killian<\/p>\n","protected":false},"template":"","class_list":["post-42891","topic","type-topic","status-closed","hentry"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_old_topic_id":["15072"],"_bbp_old_topic_author_name_id":["Anonymous"],"_bbp_old_is_topic_anonymous_id":["false"],"_bbp_old_closed_status_id":["publish"],"_bbp_author_ip":[null],"_bbp_old_sticky_status_id":["normal"],"_bbp_likes_count":["0"],"_btv_view_count":["1384"],"_bbp_subscription":["240687"],"_bbp_topic_status":["unanswered"],"_bbp_status":["publish"],"_bbp_topic_id":["42891"],"_bbp_forum_id":["27792"],"_bbp_engagement":["118","166009"],"_bbp_voice_count":["2"],"_bbp_reply_count":["3"],"_bbp_last_reply_id":["112012"],"_bbp_last_active_id":["112012"],"_bbp_last_active_time":["2020-04-21 21:19:07"]},"test":"killian153"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/42891","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/42891\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=42891"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}