{"id":41,"date":"2021-10-02T23:21:34","date_gmt":"2021-10-02T23:21:34","guid":{"rendered":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/chapter\/viscous-flux\/"},"modified":"2022-01-10T17:56:34","modified_gmt":"2022-01-10T17:56:34","slug":"viscous-flux","status":"publish","type":"chapter","link":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/chapter\/viscous-flux\/","title":{"raw":"2.1  Viscous Flux","rendered":"2.1  Viscous Flux"},"content":{"raw":"<div class=\"viscous-flux\">\r\n<p class=\"import-Normal\">The viscous flux appearing throughout this book is the volume flux calculated by Darcy\u2019s law. This calculation is most familiar in the context of water flow and appears in various forms throughout the GW-Project books. Here, we express Darcy\u2019s law in a form suitable for the calculation of the viscous contribution to the flow of the gas phase. As is the case for water and other fluids, the driving forces are the pressure gradient and the body force due to gravity. The latter may be important when the gas column of interest is very thick (Thorstenson and Pollock, 1989), but we include only the pressure gradient in the following developments and write Darcy\u2019s law as Equation\u00a01.<\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle v=-\\frac{k_{g}}{\\mu }\\frac{dp}{dl}[\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(1)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nwhere:\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\" border=\"0\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>v<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">macroscopic volume flux (specific discharge, L<sup>3<\/sup> \/ L<sup>2<\/sup> T = L\/T)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>p<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">pressure of the gas (F\/L<sup>2<\/sup>)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>l<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">coordinate along which the motion occurs (L)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03bc<\/em><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">dynamic viscosity (FT\/L<sup>2<\/sup>) (which we treat as a constant in all subsequent developments)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>k<\/em><sub><em>g<\/em><\/sub><\/td>\r\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\r\n<td style=\"width: 83%; vertical-align: top;\">permeability to gas (L<sup>2<\/sup>) (associated with the resistance to gas motion that arises solely from viscous shear at the pore scale)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">The subscript, <em class=\"import-Eqinline\">g<\/em>, distinguishes the gas permeability (<em class=\"import-Eqinline\">k<\/em><sub class=\"import-Eqinline\"><em>g<\/em><\/sub>, also known as effective permeability) from the intrinsic permeability, <em class=\"import-Eqinline\">k<\/em>, of the porous medium. This distinction is required because the co-existence of liquids in the pore space causes the gas permeability to be smaller than the intrinsic permeability, sometimes dramatically so (e.g., Brooks and Corey, 1966). The gas permeability <em class=\"import-Eqinline\">k<\/em><sub class=\"import-Eqinline\"><em>g<\/em><\/sub> is equal to the intrinsic permeability when the porous medium is dry. The reduction in gas permeability due to the presence of liquid in the pore space is related to the concepts of effective and relative permeability as described by Corey (1994).<\/p>\r\n<p class=\"import-Normal\">Multiplication of <em class=\"import-Eqinline\">v<\/em> by the total molar concentration of the gas <em class=\"import-Eqinline\">C<\/em> (total moles of all components per unit volume), followed by use of the ideal gas law, <em class=\"import-Eqinline\">p<\/em><em class=\"import-Eqinline\">\u00a0<\/em><em class=\"import-Eqinline\">=<\/em><em class=\"import-Eqinline\">\u00a0<\/em><em class=\"import-Eqinline\">RTC<\/em>, wherein <em class=\"import-Eqinline\">R<\/em> is the gas constant (FL\/moles T) and <em class=\"import-Eqinline\">T<\/em> is absolute temperature, gives Equation\u00a02 for the viscous mole flux of the gas phase.<a id=\"equation_2\"><\/a><\/p>\r\n\r\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\r\n<tbody>\r\n<tr>\r\n<td style=\"width: 10%;\"><\/td>\r\n<td style=\"width: 80%; text-align: center;\">[latex]\\displaystyle N^{v}=-\\frac{k_{g}p}{\\mu RT}\\frac{dp}{dl}=-\\frac{k_{g}p}{\\mu }\\frac{dC}{dl}[\/latex]<\/td>\r\n<td style=\"width: 10%; text-align: right;\">(2)<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<p class=\"import-Normal\">The flux calculated by Equation\u00a02 is a contributor to the motion of the mixture as a whole. Being a component of the mixture, an individual species is carried along (advected) by the viscous contribution to the phase motion. The advection mole flux of an individual species due to the viscous flux is calculated according to [latex]N_{i}^{v}=x_{i}N^{v}[\/latex], <em>i<\/em><em> = A, B<\/em> where <em>x<\/em><sub><em>i<\/em><\/sub> is the mole fraction of component <em>i<\/em> in the mixture. The mole fraction is defined by <em>x<\/em><sub><em>i<\/em><\/sub> <em>=<\/em> <em>C<\/em><sub><em>i<\/em><\/sub><em>\/C<\/em>, where <em>C<\/em><sub><em>i<\/em><\/sub> is the molar concentration of constituent <em>i<\/em>. The sum of mole fractions over all constituents is always unity.<\/p>\r\n\r\n<\/div>","rendered":"<div class=\"viscous-flux\">\n<p class=\"import-Normal\">The viscous flux appearing throughout this book is the volume flux calculated by Darcy\u2019s law. This calculation is most familiar in the context of water flow and appears in various forms throughout the GW-Project books. Here, we express Darcy\u2019s law in a form suitable for the calculation of the viscous contribution to the flow of the gas phase. As is the case for water and other fluids, the driving forces are the pressure gradient and the body force due to gravity. The latter may be important when the gas column of interest is very thick (Thorstenson and Pollock, 1989), but we include only the pressure gradient in the following developments and write Darcy\u2019s law as Equation\u00a01.<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-content\/ql-cache\/quicklatex.com-764583f2e9c5f371b75cae9f1a8d717d_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#118;&#61;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#107;&#95;&#123;&#103;&#125;&#125;&#123;&#92;&#109;&#117;&#32;&#125;&#92;&#102;&#114;&#97;&#99;&#123;&#100;&#112;&#125;&#123;&#100;&#108;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"38\" width=\"83\" style=\"vertical-align: -15px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(1)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>where:<\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>v<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">macroscopic volume flux (specific discharge, L<sup>3<\/sup> \/ L<sup>2<\/sup> T = L\/T)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>p<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">pressure of the gas (F\/L<sup>2<\/sup>)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>l<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">coordinate along which the motion occurs (L)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>\u03bc<\/em><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">dynamic viscosity (FT\/L<sup>2<\/sup>) (which we treat as a constant in all subsequent developments)<\/td>\n<\/tr>\n<tr>\n<td style=\"width: 15%; text-align: right; vertical-align: top;\"><em>k<\/em><sub><em>g<\/em><\/sub><\/td>\n<td style=\"width: 2%; text-align: center; vertical-align: top;\">=<\/td>\n<td style=\"width: 83%; vertical-align: top;\">permeability to gas (L<sup>2<\/sup>) (associated with the resistance to gas motion that arises solely from viscous shear at the pore scale)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">The subscript, <em class=\"import-Eqinline\">g<\/em>, distinguishes the gas permeability (<em class=\"import-Eqinline\">k<\/em><sub class=\"import-Eqinline\"><em>g<\/em><\/sub>, also known as effective permeability) from the intrinsic permeability, <em class=\"import-Eqinline\">k<\/em>, of the porous medium. This distinction is required because the co-existence of liquids in the pore space causes the gas permeability to be smaller than the intrinsic permeability, sometimes dramatically so (e.g., Brooks and Corey, 1966). The gas permeability <em class=\"import-Eqinline\">k<\/em><sub class=\"import-Eqinline\"><em>g<\/em><\/sub> is equal to the intrinsic permeability when the porous medium is dry. The reduction in gas permeability due to the presence of liquid in the pore space is related to the concepts of effective and relative permeability as described by Corey (1994).<\/p>\n<p class=\"import-Normal\">Multiplication of <em class=\"import-Eqinline\">v<\/em> by the total molar concentration of the gas <em class=\"import-Eqinline\">C<\/em> (total moles of all components per unit volume), followed by use of the ideal gas law, <em class=\"import-Eqinline\">p<\/em><em class=\"import-Eqinline\">\u00a0<\/em><em class=\"import-Eqinline\">=<\/em><em class=\"import-Eqinline\">\u00a0<\/em><em class=\"import-Eqinline\">RTC<\/em>, wherein <em class=\"import-Eqinline\">R<\/em> is the gas constant (FL\/moles T) and <em class=\"import-Eqinline\">T<\/em> is absolute temperature, gives Equation\u00a02 for the viscous mole flux of the gas phase.<a id=\"equation_2\"><\/a><\/p>\n<table style=\"border: none; border-collapse: collapse; width: 100%;\">\n<tbody>\n<tr>\n<td style=\"width: 10%;\"><\/td>\n<td style=\"width: 80%; text-align: center;\"><img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-content\/ql-cache\/quicklatex.com-2fae2e353828c96a73653f9c7897b8ce_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#92;&#100;&#105;&#115;&#112;&#108;&#97;&#121;&#115;&#116;&#121;&#108;&#101;&#32;&#78;&#94;&#123;&#118;&#125;&#61;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#107;&#95;&#123;&#103;&#125;&#112;&#125;&#123;&#92;&#109;&#117;&#32;&#82;&#84;&#125;&#92;&#102;&#114;&#97;&#99;&#123;&#100;&#112;&#125;&#123;&#100;&#108;&#125;&#61;&#45;&#92;&#102;&#114;&#97;&#99;&#123;&#107;&#95;&#123;&#103;&#125;&#112;&#125;&#123;&#92;&#109;&#117;&#32;&#125;&#92;&#102;&#114;&#97;&#99;&#123;&#100;&#67;&#125;&#123;&#100;&#108;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"38\" width=\"204\" style=\"vertical-align: -15px;\" \/><\/td>\n<td style=\"width: 10%; text-align: right;\">(2)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p class=\"import-Normal\">The flux calculated by Equation\u00a02 is a contributor to the motion of the mixture as a whole. Being a component of the mixture, an individual species is carried along (advected) by the viscous contribution to the phase motion. The advection mole flux of an individual species due to the viscous flux is calculated according to <img loading=\"lazy\" decoding=\"async\" src=\"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-content\/ql-cache\/quicklatex.com-7c7697f2f0eac6a024a14b286c60d87c_l3.png\" class=\"ql-img-inline-formula quicklatex-auto-format\" alt=\"&#78;&#95;&#123;&#105;&#125;&#94;&#123;&#118;&#125;&#61;&#120;&#95;&#123;&#105;&#125;&#78;&#94;&#123;&#118;&#125;\" title=\"Rendered by QuickLaTeX.com\" height=\"17\" width=\"82\" style=\"vertical-align: -5px;\" \/>, <em>i<\/em><em> = A, B<\/em> where <em>x<\/em><sub><em>i<\/em><\/sub> is the mole fraction of component <em>i<\/em> in the mixture. The mole fraction is defined by <em>x<\/em><sub><em>i<\/em><\/sub> <em>=<\/em> <em>C<\/em><sub><em>i<\/em><\/sub><em>\/C<\/em>, where <em>C<\/em><sub><em>i<\/em><\/sub> is the molar concentration of constituent <em>i<\/em>. The sum of mole fractions over all constituents is always unity.<\/p>\n<\/div>\n","protected":false},"author":1,"menu_order":6,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":[],"pb_section_license":""},"chapter-type":[],"contributor":[],"license":[],"class_list":["post-41","chapter","type-chapter","status-publish","hentry"],"part":86,"_links":{"self":[{"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/chapters\/41","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/wp\/v2\/users\/1"}],"version-history":[{"count":8,"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/chapters\/41\/revisions"}],"predecessor-version":[{"id":397,"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/chapters\/41\/revisions\/397"}],"part":[{"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/parts\/86"}],"metadata":[{"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/chapters\/41\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/wp\/v2\/media?parent=41"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/pressbooks\/v2\/chapter-type?post=41"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/wp\/v2\/contributor?post=41"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/books.gw-project.org\/flux-equations-for-gas-diffusion-in-porous-media\/wp-json\/wp\/v2\/license?post=41"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}