{"id":4865,"date":"2012-11-22T16:56:43","date_gmt":"2012-11-22T05:56:43","guid":{"rendered":"https:\/\/www.computationalfluiddynamics.com.au\/?p=687"},"modified":"2012-11-22T16:56:43","modified_gmt":"2012-11-22T05:56:43","slug":"interpreting-results-for-superficial-and-true-velocity","status":"publish","type":"post","link":"https:\/\/www.leapaust.com.au\/blog\/cfd\/interpreting-results-for-superficial-and-true-velocity\/","title":{"rendered":"Tips &#038; Tricks: How to interpret results for multiphase &#038; porous domains using true velocity and superficial velocity"},"content":{"rendered":"<div id=\"bsf_rt_marker\"><\/div><p>A commonly asked question is: <em>What is Superficial Velocity and when do I use it<\/em>?<\/p>\n<p>&nbsp;<\/p>\n<p>If a fluid flows through a region that is occupied by either fixed structures (a porous region, pipe rack, catalyst bed, etc&#8230;),\u00a0or shares the channel with other fluids (e.g. gas-liquid flow), there are two ways to describe the fluid velocity.<\/p>\n<p>&nbsp;<\/p>\n<p>The first is to use the &#8220;True Velocity&#8221; which is the actual velocity of the fluid particles. This velocity will vary with location in the porous matrix. This is the velocity you would measure experimentally if you\u00a0focused on a small region of fluid.<\/p>\n<p>&nbsp;<\/p>\n<p>The second is the &#8220;Superficial Velocity&#8221;, which is the velocity the fluid would have if it were flowing through the same domain without any obstructions. This is a very useful quantity as in incompressible flow it is conserved regardless of the variation of the porosity. Therefore, even at a boundary between a porous region and continuous fluid the superficial velocity is unchanged, whereas the true velocity must increase in the porous region so that mass is conserved.<\/p>\n<p>&nbsp;<\/p>\n<p>These two velocities are easily related via <a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-eq0a.png\"><img decoding=\"async\" class=\"alignnone size-full wp-image-732\" title=\"porous-eq0a\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-eq0a.png\" alt=\"\" width=\"77\" height=\"23\" \/><\/a>, where <em>V<sub>s<\/sub><\/em> is the superficial velocity, <em>V<sub>true<\/sub><\/em> is the true velocity and\u00a0 <em>\u03b5<\/em> is the porosity (the local fraction of the volume occupied by the fluid).<\/p>\n<p>&nbsp;<\/p>\n<p>If you study <a title=\"10 questions to ask yourself when tackling your first (or a new) Multiphase CFD project\" href=\"https:\/\/www.computationalfluiddynamics.com.au\/10-questions-when-tackling-your-first-multiphase-cfd-project\/\">multiphase flows<\/a> you will certainly encounter superficial velocity as it is used to characterise a flow system. Using superficial velocity has the benefit that it is conserved (for an incompressible flow with no phase change) regardless of the complexity of the flow regime, e.g. if the flow regime changes from bubbly to slug flow, the superficial velocity stays constant even though the local velocity varies. Maps plotting the gas superficial velocity on one axis and liquid superficial velocity on the other\u00a0 are known as regime maps and are used to define the boundary between different regimes.\u00a0\u00a0When visualising the results of\u00a0multiphase flow simulations, ANSYS CFD-Post will automatically give the user the choice of plotting either Superficial Velocity or (True) Velocity variables.<\/p>\n<p>&nbsp;<\/p>\n<p>The use of a superficial velocity is also often encountered when dealing with pressure drop correlations for porous regions, be it a true porosity or a porosity used to represent flow obstructions.\u00a0 An experimentalist can use either superficial or true velocities to characterise their system, but when reviewing experimental data it is worth knowing that it is more common to use superficial velocity as this can be measured outside of the porous region.<\/p>\n<p>&nbsp;<\/p>\n<p>Consider Darcy&#8217;s law for slow flow (negligible inertia) in a porous medium, which relates the volumetric flow through a given face to the pressure drop via\u00a0<a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-eq1a.png\"><img decoding=\"async\" class=\"alignnone size-full wp-image-733\" title=\"porous-eq1a\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-eq1a.png\" alt=\"\" width=\"77\" height=\"36\" \/><\/a> where q is the volumetric flowrate,\u00a0 <em>\u03ba<\/em> is the permeability, <em>\u03bc<\/em> is the dynamic viscosity and <em>p<\/em> is the pressure. Here <em>q<\/em> can be replaced by the flow area multiplied by the superficial velocity, so that the above\u00a0 equation becomes: <a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-eq2a.png\"><img decoding=\"async\" class=\"alignnone size-full wp-image-734\" title=\"porous-eq2a\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-eq2a.png\" alt=\"\" width=\"82\" height=\"31\" \/><\/a>\u00a0where now<em> v<\/em> must be the superficial velocity.<\/p>\n<p>&nbsp;<\/p>\n<p>To help demonstrate the difference between superficial velocity and true velocity in a CFD model with a porous domain, let\u2019s look at the flow through a catalytic converter. A porous domain is used to represent the ceramic monolith substrate.\u00a0Exhaust gas\u00a0flows through the inlet, passes through an expansion which leads to the porous domain, and then flows through a contraction and the outlet. The substrate is impermeable in the X and Y direction, and the loss in the Z-direction is represented by specifying the volume porosity and higher order loss coefficients.<\/p>\n<p><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image11.png\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-medium wp-image-726\" title=\"porous-image1\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image11-300x177.png\" alt=\"\" width=\"300\" height=\"177\" srcset=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image11-300x177.png 300w, https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image11.png 391w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>When post-processing the model in CFD-Post, we can visualize the pressure drop due to the substrate by inserting a two-dimensional plane, or similarly, using a chart.<\/p>\n<p><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image22.png\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-medium wp-image-697\" title=\"porous-image2\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image22-300x116.png\" alt=\"\" width=\"300\" height=\"116\" \/><\/a><\/p>\n<p>The pressure drop due to the porous domain is quite evident using both visualization techniques. We can use the same technique to plot the flow velocity:<\/p>\n<p><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image3.png\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-medium wp-image-698\" title=\"porous-image3\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image3-300x116.png\" alt=\"\" width=\"300\" height=\"116\" \/><\/a><\/p>\n<p>&nbsp;<\/p>\n<p>It is\u00a0evident that a jump in the plotted velocity occurs within the porous domain. This is because the velocity visualised by CFD-Post is the actual or true velocity, and due to the porous obstructions the flow will speed up within the porous region. \u00a0To plot the representative or superficial velocity of the flow through the substrate, we can create a User Defined Variable in CFD-Post which references an expression of velocity multiplied by area porosity (or equivalently the volume porosity in an isotropic material), as follows:<\/p>\n<p>Insert an expression by clicking the <em>Expressions<\/em> Tab in CFD-Post and naming the expression appropriately. Set the value as Velocity*Volume Porosity for an isotropic model. Then insert a new variable by clicking the <em>Variables <\/em>Tab, <em>Insert New Variable<\/em> and reference your previous expression.<\/p>\n<p><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image4.png\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-medium wp-image-699\" title=\"porous-image4\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image4-300x80.png\" alt=\"\" width=\"300\" height=\"80\" \/><\/a><\/p>\n<p>&nbsp;<\/p>\n<p>We can now plot Superficial Velocity as a flow variable on our central plane and in a line chart (as below) and, as you would expect, we no longer observe any discontinuity at the porous domain boundaries.<\/p>\n<p><a href=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image5.png\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-medium wp-image-700\" title=\"porous-image5\" src=\"https:\/\/www.leapaust.com.au\/blog\/wp-content\/uploads\/2012\/11\/porous-image5-300x115.png\" alt=\"\" width=\"300\" height=\"115\" \/><\/a><\/p>\n<p>We hope that this\u00a0post has helped\u00a0clear up some common misconceptions about the use of superficial and true velocities in CFD simulations.\u00a0 Please let us know in the comments below if this has helped you with your CFD model, or if you have any further questions.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>A commonly asked question is: What is Superficial Velocity and when do I use it? &nbsp; If a fluid flows through a region that is occupied by either fixed structures (a porous region, pipe rack, catalyst bed, etc&#8230;),\u00a0or shares the channel with other fluids (e.g. gas-liquid flow), there are two ways to describe the fluid&hellip;&nbsp;<a href=\"https:\/\/www.leapaust.com.au\/blog\/cfd\/interpreting-results-for-superficial-and-true-velocity\/\" rel=\"bookmark\">Read More &raquo;<span class=\"screen-reader-text\">Tips &#038; Tricks: How to interpret results for multiphase &#038; porous domains using true velocity and superficial velocity<\/span><\/a><\/p>\n","protected":false},"author":3,"featured_media":724,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"neve_meta_sidebar":"","neve_meta_container":"","neve_meta_enable_content_width":"","neve_meta_content_width":0,"neve_meta_title_alignment":"","neve_meta_author_avatar":"","neve_post_elements_order":"","neve_meta_disable_header":"","neve_meta_disable_footer":"","neve_meta_disable_title":"","neve_meta_reading_time":"","footnotes":""},"categories":[323],"tags":[171,174,396,399,456,468],"class_list":["post-4865","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-cfd","tag-ansys","tag-ansys-cfd","tag-cfd-simulation-accuracy","tag-common-misconceptions","tag-multiphase-flow","tag-porous-domains"],"_links":{"self":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4865","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/comments?post=4865"}],"version-history":[{"count":0,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/posts\/4865\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/media\/724"}],"wp:attachment":[{"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/media?parent=4865"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/categories?post=4865"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.leapaust.com.au\/blog\/wp-json\/wp\/v2\/tags?post=4865"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}