• PVPMC Blog
  • PV_LIB Toolbox
  • Events and Workshops
  • Contact Us
Header Menu
PV Performance Modeling Collaborative
PV Performance Modeling Collaborative
An Industry and National Laboratory collaborative to improve Photovoltaic Performance Modeling
An Industry and National Laboratory collaborative to Improve Photovoltaic Performance Modeling
  • Home
  • Modeling Steps
    • 1. Weather and Design
      • Sun Position
      • Irradiance & Insolation
      • Weather Observations
      • Array Orientation
      • Plane of Array (POA) Irradiance
      • Shading, Soiling, and Reflection Losses
    • 2. DC Module IV Characteristics
      • Module Temperature
      • Cell Temperature
      • Effective Irradiance
      • Single Diode Equivalent Circuit Models
      • Point-value models
    • 3. DC Array IV
      • Mismatch Losses
      • DC Component Health
      • DC Wiring Losses
      • Array Utilization
    • 4. DC to AC Conversion
      • CEC Inverter Test Protocol
      • Operating Temperature
      • Sandia Inverter Model
      • Driesse Inverter Model
      • Inverter Saturation or “Clipping”
      • Loss of Grid
      • Advanced Inverter Features
    • 5. AC System Output
      • AC Wiring Losses
      • PV Performance Metrics
  • PV Research Projects
    • Bifacial PV Project
      • Outdoor Bifacial PV Performance Data
        • Field Example of Bifacial Gain at Sandia
      • Bifacial PV Characterization and Rating Standards
      • Bifacial PV Performance Models
        • Ray Tracing Models for Backside Irradiance
          • RADIANCE models
          • COMSOL models
        • View Factor Models
          • Sandia View Factor Model Implementation
          • NREL View Factor Model Implementation
    • PV Lifetime Project
      • PV Lifetime Modules
      • PV Lifetime Field Deployments
      • Indoor Module Flash Testing
      • PV Degradation Modeling
  • Applications & Tools
    • PV_LIB Toolbox
    • GridPV Toolbox
    • Wavelet Variability Model
      • WVM Square Plant Example
      • Polygon Vertices to Define Plant Footprint Example
      • WVM Discrete Point Example
    • Solar Variability Zones
    • Spectral irradiance dataset from Albuquerque
  • Resources and Events
    • Events and Workshops
    • PVPMC Blog
    • Sign up for Email Newsletter
    • IEA PVPS Task 13 Reports
    • 2021 Blind PV Modeling Comparison
MENU CLOSE back  

COMSOL models

You are here:
  1. Home
  2. PV Research Projects
  3. Bifacial PV Project
  4. Bifacial PV Performance Models
  5. Ray Tracing Models for Backside Irradiance
  6. COMSOL models

Ray Tracing Work at the University of Iowa

COMSOL is a finite element method (FEM) solver and the Ray Optics Module can be used to model electromagnetic wave propagation in systems in which the wavelength is much smaller than the smallest geometric detail in the model. The electromagnetic waves are treated as rays that can propagate through homogeneous or graded media. Because it is not necessary to resolve the wavelength with a finite element mesh, ray trajectories can be computed over long distances at a low computational cost. Rays can also undergo reflection and refraction at boundaries between different media.[1] COMSOL’s Ray Optics package uses a forward ray tracing model that can run on a high performance computing (HPC) platform to model large domains and many separate rays.

Preliminary simulations of backside irradiance have been made to explore the effect of cell spacing within a bifacial module.  For example, Prism Solar offers two versions of its bifacial solar modules.  One fills the module area with 60 PV cells and the other uses 48 cells in the same area but includes transparent spaces between cells to allow more light to pass through the module.

Comsol_Prism_gridView of the 60 cell (Left) and 48 cell (Right) modules

The modeling was performed assuming a single module was installed in landscape orientation at 35° tilt (Albuquerque’s latitude) angle and 1 m above the ground surface (measured from the bottom edge of the module).  Three different types of roof surfaces were considered: soiled beige roof (37% albedo), clean beige roof (55% albedo) and white beige roof (83% albedo), based on measurements made by NREL.

In this preliminary example, rays are released with an incidence angle of 0° to the module surface (this represents the solar position at solar noon on the equinox).  At this stage, only direct beam radiation is simulated (diffuse radiation is ignored).

Released rays get scattered at the module and ground surfaces and reach the back of the module later in time (in this case, because rays are normal to module, there is no refraction by module). Released rays find their way to ground except when they are absorbed by the cells in the module.

Comsol_Prism_modelPosition of the module at 1m above of the ground (left) and cross section of the module and released rays with incidence angle of 35° (normal to module) with initial intensity of 1000 W/ (right)

Comsol_Prism_model2Due to different reflectivity of ground, intensity of scattered rays from the ground are different. These figures shows the simulation results at 13.6 ns for soiled beige roof (left), clean beige roof (middle) and white roof (right), when the rays are scattered by ground and travel upward.

Initial simulations were run on a workstation and thus used a coarse grid and a sparse set of rays.  This results in non-uniform irradiance patterns caused by the coarse grid.  Future runs will use the Univ. of Iowa’s HPC, which will result in more accurate and smoother results.  The initial coarse results are shown below.

The initial results are shown as maps of cumulative light passing through an imaginary plane coplanar with the back of the bifacial module.  Simulated sun rays hitting the module get scattered by the module and ground. Scattered rays from the ground reach back-side of the panel later in time and deposit additional power on cells. Therefore, we determined the power intensity on the back of the module at a time (8 ns) when the sun rays pass through the panel through the spacing between the cells (left figures at 8 ns) and after (20 ns) the rays are reflected back from the ground to the back of the panel (right figures). Due to high energy absorption of the silicon cells, rays only passes through the gap between the cells, which is filled with EVA (ethyl vinyl acetate). Difference of power intensity in the left and right figures shows how much additional power reaches the back of the cells in the case of a soiled beige roof, clean beige roof and white roof.

48cell_result_soiled_beige48 cell module: Deposited ray power intensity (/2) after front-side irradiance and before backside irradiance (left) and after backside irradiance (right) for soiled beige roof (reflectivity of ~37%)

60cell_result_soiled_beige60 cell module: Deposited ray power intensity (/2) after front-side irradiance and before backside irradiance (left) and after backside irradiance (right) for soiled beige roof (reflectivity of ~37%)

48cell_result_white

48 cell module: Deposited ray power intensity (W/) after front-side irradiance and before backside irradiance (left) and after backside irradiance (right) for white roof (reflectivity of ~83%)

60cell_result_white

60 cell module: Deposited ray power intensity (W/) after front-side irradiance and before backside irradiance (left) and after backside irradiance (right) for white roof (reflectivity of ~83%)

Below are backside irradiance maps (difference of top figures) for two types of roofs.

Backside_Irradiance_soiled_beigeBackside irradiance maps for the souled beige roof (37% albedo).  Left is for the 48 cell module with average intensity of 57.11 W/m^2 and right is for 60 cell with average intensity of 52.49 W/m^2.

Backside_Irradiance_whiteBackside irradiance maps for the white roof (83% albedo).  Left is for the 48 cell module with average intensity of 289.70 W/ m^2 and right is for 60 cell module with average intensity of 257.66 W/m^2

Including gaps between cells in bifacial modules increases the bifacial efficiency but is unlikely to offer significant economic advantages in energy costs over tightly spaced cell versions.  However, in certain applications, such as carports, patio covers and skylights, the added light transparency may offer benefits beyond the pure cost of energy.

Ray tracing methods, such as COMSOL’s Ray Optics Module are a promising method for evaluating back side irradiance and thus bifacial PV performance.  The University of Iowa will continue to develop more complex and sophisticated models of back side irradiance.  Stay tuned for updates on results.

[1] https://www.comsol.com/ray-optics-module

 

This modeling work was performed by Amir Asgharzadeh Shishavan, a PhD. graduate student at the University of Iowa under the supervision of Dr. Fatima Toor.

PVPMC
  • PVPMC Home
  • Modeling Steps
  • PV Research Projects
  • Applications & Tools
Recent Posts
  • Call for Abstracts for the 2023 PVPMC Workshop – Due March 6, 2023
    January 26, 2023
  • Save-the-Date for the 2023 PV Performance Modeling Collaborative Workshop – May 9-10 in Salt Lake City, UT.
    December 15, 2022
  • Registration Open for 2022 PV Performance Modeling Workshop
    May 27, 2022
Word Cloud
Irradiance Sensors Photon SAND2014-17810 W SAND2014-19827 W SAND2015-9872E SAND2016-0054 W SAND2016-0054W SAND2016-2333 W SAND2016-10263 W SAND2019-15366 W SAND2021-7933 S Sandia
Sign up for our email newsletter
Click to subscribe
PV Performance Modeling Collaborative
  • PVPMC Home
  • Modeling Steps
  • PV Research Projects
  • Applications & Tools
Footer Links

© 2018 National Technology and Engineering Solutions of Sandia, LLC.
Questions & Comments | Privacy & Security

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC., a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.