Estimate diffuse and direct component from global irradiance












0














I am looking to separate the diffuse and direct component of global irradiance and found the Erbs model to do this in pvlib (see pvlib.irradiance.erbs) however, I am getting very strange results. I would expect the Direct Normal Irradiance (DNI) to be lower than the Global Horizontal Irradiance (GHI); or am I missing something? Values of GHI are not above 800 W m^2 for these days.



enter image description here





EDIT: As per Cliff H advice, I have limited the solar zenith to less than 85 arc degrees; the results have improved however, there are large spikes in DNI values that do not appear very reasonable, e.g. start of 07-16.



enter image description here










share|improve this question





























    0














    I am looking to separate the diffuse and direct component of global irradiance and found the Erbs model to do this in pvlib (see pvlib.irradiance.erbs) however, I am getting very strange results. I would expect the Direct Normal Irradiance (DNI) to be lower than the Global Horizontal Irradiance (GHI); or am I missing something? Values of GHI are not above 800 W m^2 for these days.



    enter image description here





    EDIT: As per Cliff H advice, I have limited the solar zenith to less than 85 arc degrees; the results have improved however, there are large spikes in DNI values that do not appear very reasonable, e.g. start of 07-16.



    enter image description here










    share|improve this question



























      0












      0








      0







      I am looking to separate the diffuse and direct component of global irradiance and found the Erbs model to do this in pvlib (see pvlib.irradiance.erbs) however, I am getting very strange results. I would expect the Direct Normal Irradiance (DNI) to be lower than the Global Horizontal Irradiance (GHI); or am I missing something? Values of GHI are not above 800 W m^2 for these days.



      enter image description here





      EDIT: As per Cliff H advice, I have limited the solar zenith to less than 85 arc degrees; the results have improved however, there are large spikes in DNI values that do not appear very reasonable, e.g. start of 07-16.



      enter image description here










      share|improve this question















      I am looking to separate the diffuse and direct component of global irradiance and found the Erbs model to do this in pvlib (see pvlib.irradiance.erbs) however, I am getting very strange results. I would expect the Direct Normal Irradiance (DNI) to be lower than the Global Horizontal Irradiance (GHI); or am I missing something? Values of GHI are not above 800 W m^2 for these days.



      enter image description here





      EDIT: As per Cliff H advice, I have limited the solar zenith to less than 85 arc degrees; the results have improved however, there are large spikes in DNI values that do not appear very reasonable, e.g. start of 07-16.



      enter image description here







      pvlib solar






      share|improve this question















      share|improve this question













      share|improve this question




      share|improve this question








      edited Nov 19 at 12:36

























      asked Nov 12 at 20:23









      Kievit

      374




      374
























          1 Answer
          1






          active

          oldest

          votes


















          2














          DNI > GHI is common at low solar elevation. GHI decreases much faster than DNI as solar elevation decreases. For example, think of a clear day with the sun right near the horizon. DNI will large because it's measured on a plane normal to the sun vector, but GHI will be near zero.



          The values of DNI that are much greater than 1000 W/m2 are likely at very high zenith, since the Erbs model basically divides by cos(zenith). In practice, I limit using decomposition models like Erbs to zenith<85 degrees, to avoid the non-physical results.






          share|improve this answer























          • Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
            – Kievit
            Nov 12 at 22:53












          • The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
            – Cliff H
            Nov 14 at 16:13










          • Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
            – Kievit
            Nov 19 at 12:31











          Your Answer






          StackExchange.ifUsing("editor", function () {
          StackExchange.using("externalEditor", function () {
          StackExchange.using("snippets", function () {
          StackExchange.snippets.init();
          });
          });
          }, "code-snippets");

          StackExchange.ready(function() {
          var channelOptions = {
          tags: "".split(" "),
          id: "1"
          };
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function() {
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled) {
          StackExchange.using("snippets", function() {
          createEditor();
          });
          }
          else {
          createEditor();
          }
          });

          function createEditor() {
          StackExchange.prepareEditor({
          heartbeatType: 'answer',
          autoActivateHeartbeat: false,
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          });


          }
          });














          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53269553%2festimate-diffuse-and-direct-component-from-global-irradiance%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes









          2














          DNI > GHI is common at low solar elevation. GHI decreases much faster than DNI as solar elevation decreases. For example, think of a clear day with the sun right near the horizon. DNI will large because it's measured on a plane normal to the sun vector, but GHI will be near zero.



          The values of DNI that are much greater than 1000 W/m2 are likely at very high zenith, since the Erbs model basically divides by cos(zenith). In practice, I limit using decomposition models like Erbs to zenith<85 degrees, to avoid the non-physical results.






          share|improve this answer























          • Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
            – Kievit
            Nov 12 at 22:53












          • The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
            – Cliff H
            Nov 14 at 16:13










          • Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
            – Kievit
            Nov 19 at 12:31
















          2














          DNI > GHI is common at low solar elevation. GHI decreases much faster than DNI as solar elevation decreases. For example, think of a clear day with the sun right near the horizon. DNI will large because it's measured on a plane normal to the sun vector, but GHI will be near zero.



          The values of DNI that are much greater than 1000 W/m2 are likely at very high zenith, since the Erbs model basically divides by cos(zenith). In practice, I limit using decomposition models like Erbs to zenith<85 degrees, to avoid the non-physical results.






          share|improve this answer























          • Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
            – Kievit
            Nov 12 at 22:53












          • The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
            – Cliff H
            Nov 14 at 16:13










          • Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
            – Kievit
            Nov 19 at 12:31














          2












          2








          2






          DNI > GHI is common at low solar elevation. GHI decreases much faster than DNI as solar elevation decreases. For example, think of a clear day with the sun right near the horizon. DNI will large because it's measured on a plane normal to the sun vector, but GHI will be near zero.



          The values of DNI that are much greater than 1000 W/m2 are likely at very high zenith, since the Erbs model basically divides by cos(zenith). In practice, I limit using decomposition models like Erbs to zenith<85 degrees, to avoid the non-physical results.






          share|improve this answer














          DNI > GHI is common at low solar elevation. GHI decreases much faster than DNI as solar elevation decreases. For example, think of a clear day with the sun right near the horizon. DNI will large because it's measured on a plane normal to the sun vector, but GHI will be near zero.



          The values of DNI that are much greater than 1000 W/m2 are likely at very high zenith, since the Erbs model basically divides by cos(zenith). In practice, I limit using decomposition models like Erbs to zenith<85 degrees, to avoid the non-physical results.







          share|improve this answer














          share|improve this answer



          share|improve this answer








          edited Nov 19 at 16:02









          Kievit

          374




          374










          answered Nov 12 at 21:44









          Cliff H

          711




          711












          • Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
            – Kievit
            Nov 12 at 22:53












          • The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
            – Cliff H
            Nov 14 at 16:13










          • Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
            – Kievit
            Nov 19 at 12:31


















          • Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
            – Kievit
            Nov 12 at 22:53












          • The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
            – Cliff H
            Nov 14 at 16:13










          • Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
            – Kievit
            Nov 19 at 12:31
















          Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
          – Kievit
          Nov 12 at 22:53






          Thanks for this. I have limited the zenith to < 85 arc degrees. This has definitely improved the results however, I still get funny artefacts.
          – Kievit
          Nov 12 at 22:53














          The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
          – Cliff H
          Nov 14 at 16:13




          The results appear credible to me. The Erbs model is estimating the diffuse component of GHI, and returns DNI by a closure calculation, i.e., DNI = (GHI - DHI) / cos(Z). The empirical equation that obtains DHI from GHI is a fit through a fairly broad scatter of data, so it should be regarded as one of many possible values of DHI that correspond to the observed GHI value.
          – Cliff H
          Nov 14 at 16:13












          Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
          – Kievit
          Nov 19 at 12:31




          Thanks. I did not consider that DNI assumes the surface perpendicular to the rays. As this is slightly different, I have asked a new question here on how to estimate DNI and DHI for a horizontal surface from GHI.
          – Kievit
          Nov 19 at 12:31


















          draft saved

          draft discarded




















































          Thanks for contributing an answer to Stack Overflow!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          To learn more, see our tips on writing great answers.





          Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


          Please pay close attention to the following guidance:


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53269553%2festimate-diffuse-and-direct-component-from-global-irradiance%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          Xamarin.iOS Cant Deploy on Iphone

          Glorious Revolution

          Dulmage-Mendelsohn matrix decomposition in Python