Fill out the form to download

Required field
Required field
Not a valid email address
Required field

Solar Load

The electromagnetic radiation emitted by the sun is also referred to as solar radiation. In Computational Fluid Dynamics (CFD), solar radiation acts as an additional heat source. For some applications, such as thermal comfort studies, ignoring the solar radiation factor may impact the accuracy of the results.

In such cases, SimScale has a solar radiation module, available in the Conjugate Heat Transfer v2.0 (CHT v2.0) analysis type, under the global settings, as shown in Figure 1:

toggling on solar radiation in simscale cht v2.0
Figure 1: The Solar Load toggle is under the global simulation settings for the CHT v2.0 analysis type.

The solar radiation settings are found under Solar calculator in the simulation tree. Two main settings are required in the setup: the Sun direction, and the Solar load definitions. We will go through the details below.

Important

The current version of the solar load model has three limitations:

1. It is not coupled with ordinary radiation. Therefore, from Figure 1, it’s not possible to have both Solar Load and Radiation toggles enabled at the same time
2. Once a ray hits an opaque boundary, no reflection is considered
3. Solid volumes are not supported when using Solar Load. For this reason, even though we are using the CHT v2.0 analysis type, the CAD requirements are the same as in a Convective Heat Transfer simulation – that is, there should only be a flow region and no solid parts in the final CAD model.

solar load cad requirements
Figure 2: The CAD model should contain only the flow volume, and no solid parts (like walls). Please note that there should not be any material assignment under “Solids”.

Sun Direction

Under sun direction, the user can determine the position of the sun for their simulation. This is important because the direction of the solar rays changes based on the position of the sun.

In the Workbench, a set of yellow arrows will indicate the current direction of the sun rays, based on the Sun direction configuration. There are two options available for this setting: Time and place and Custom.

Time and Place

This configuration determines the position of the sun relative to the geometry for a given date, time, and location. It is worth mentioning that the sky is assumed to be in the positive z-direction, and the north direction is assumed to be in the positive y-direction. In case your model has a different orientation for the north direction, you can adjust the North angle in the settings.

Figure 3 shows the Time and place configuration steps:

time and place solar calculator settings
Figure 3: In the Workbench, the current north direction is shown by the black arrow, and the sun rays are represented by the yellow arrows.
  1. If the north direction is not positive y in your model, you can adjust the North angle.
  2. In the interactive map, type the location where the model that you are simulating is located.
  3. Choose a Date and time of interest.

SimScale will calculate the position of the sun relative to your model, given the date, time, and location.

Custom

The second option is a custom definition of the Sun direction vector. Based on the orientation cube at the bottom-right of the Workbench, the user can define the direction of the sun’s rays:

custom definition of sun direction solar load
Figure 4: With a custom sun direction, the user can define any direction of interest for the solar rays.

Solar Load

After setting up the sun direction, we can define the Direct and the Diffuse solar loads that will be considered. Please check the next section in this documentation to see how different boundary conditions react to the solar load.

When an external opaque boundary of the model is exposed to either direct or diffuse radiation, the rays will act as a heat source to the boundary. With that being said, both direct and diffuse solar loads are considered external when they are absorbed by the external boundaries of the computational domain.

Furthermore, a direct solar load can also be absorbed by internal boundaries – This happens when a ray enters the domain via a transparent or semi-transparent boundary and hits an opaque boundary before exiting through another transparent or semi-transparent boundary. In this case, the solar load is internal.

For example, let’s consider a meeting room with 3 transparent/semi-transparent boundaries, shown in the left-hand side room from Figure 5. In this case, the yellow arrows on the right-hand side room represent the internal solar load, while the blue arrows represent the external solar load.

internal and external solar load direct and diffuse
Figure 5: A direct solar load may enter the model through a transparent/semitransparent boundary, such as a window, and hit an internal opaque boundary, such as occupants or a table.

Find below the two supported configurations for a solar load.

Custom

With a Custom definition for the solar load, the user can explicitly define the direct and diffuse solar load:

solar load custom definition direct and diffuse
Figure 6: With a custom definition, the user explicitly defines the solar load values.

This option is useful in case you are interested in a specific value for the solar load.

Fair Weather Conditions

The Fair Weather Conditions method is an alternative way to define the diffuse and direct solar loads. For the underlying equations, please visit Table 14 from Chapter 29 in the 2001 ASHRAE Fundamentals Handbook\(^{[1]}\).

This approach assumes that the sky is in the positive z-direction, and the positive y-direction points to the north. For this reason, if you are using a Custom sun direction, the Sun direction vector must be negative in z. Otherwise, the sun would be below the horizon.

For the fair weather method, the user needs to provide a series of parameters, as shown below:

fair weather conditions method input
Figure 7: Based on the input settings, the platform will calculate the direct and diffuse solar loads.
  • The Sky cloud cover fraction represents the fraction of the sky which is taken by clouds.
  • The Ground reflectivity depends on the color and texture of the ground around the computational domain.
  • The remaining terms are related to extraterrestrial solar irradiance. Find below an adaptation of table 7, Chapter 30, from the 2001 ASHRAE Fundamentals Handbook\(^{[1]}\):
Apparent Solar Irradiation \([\frac{W}{m²}]\)Atmospheric ExtinctionDiffuse Radiation Factor
January12300.1420.058
February12150.1440.060
March11860.1560.071
April11360.1800.097
May11040.1960.121
June10880.2050.134
July10850.2070.136
August11070.2010.122
September11510.1770.092
October11920.1600.073
November12210.1490.063
December12330.1420.057
Table 1: Values for Apparent solar irradiation, Atmospheric extinction, and Diffuse radiation factor, obtained on the 21st of each month, during the year of 1964. Adapted from 2001 ASHRAE Fundamentals Handbook\(^{[1]}\).

Interactions with Boundary Conditions

The coupling of the solar load with the temperature field of the domain is done via boundary conditions. Once a boundary is hit by a ray, it will either absorb or transmit a fraction, or all of the incoming solar load. All boundary conditions transmit solar load in the same way, however, depending on their Temperature type definition, boundary conditions react differently to solar loads absorbed internally and externally.

temperature type definition cht v2
Figure 8: The Temperature type definition determines how internal and external solar loads impact the boundary.

Furthermore, under Solar radiative behavior, we define if a boundary is Transparent, Semi-transparent, or Opaque. Typically, inlets and outlets are transparent, whereas boundaries that do not allow light to travel through them are opaque. Certain materials, such as glass, are semi-transparent since they allow some light through.

Below, we are going to expand on the different Temperature type settings for a wall boundary condition.

Fixed Temperature

Since the boundary has a fixed temperature, the solar load does not have any effect.

Adiabatic

For internally absorbed solar load, the adiabatic boundary condition becomes a fixed gradient, forcing the wall to transfer heat to the adjacent fluid cells.

In the event of external solar load, the adiabatic configuration of the boundary prevents any external heat flux from entering the computational domain. Therefore, the external rays are completely blocked.

External Wall Heat Flux

When selecting External wall heat flux the user needs to define the Heat flux type, as seen in Figure 8. There are three options available:

  1. Fixed heat flux: Since we define the fixed heat flux, we do not take any external solar load into account. However, the internal solar load is added to the fixed heat flux, resulting in more heat being dissipated to the adjacent flow region cells.
  2. Fixed power: The behavior is the same as fixed heat flux. The external solar load does not affect this boundary condition, whereas internal solar loads will be added as a heat flux based on the surface area of the boundary.
  3. Derived heat flux: Both internal and external solar loads have an impact when using this configuration. The heat flux through this boundary is defined by the Heat transfer coefficient, the temperature differences, and the solar loads.

Note

When using the Derived heat flux configuration, if the boundary is not exposed to the external environment (for example, an internal wall of a building), the Outer face absorptivity should be set to zero.

outer surface absorptivity definition
Figure 9: Outer surface absorptivity for faces in the shadow should be zero.

Post-processing

When simulating solar loads, it is possible to obtain the internal and external loads in the post-processor. For example, using the solar load settings from Figure 4, we have the following results:

external solar load results simscale
Figure 10: For the external solar load, notice how the shadow created by the outlet section above the ceiling affects the external solar load.

Now, inspecting the internal solar load, only the rays that go through the window (semi-transparent material) and the inlet/outlet (transparent) are taken into account.

internal solar radiation results simscale
Figure 11: The internal opaque boundaries hit with incoming rays through the window will dissipate additional heat.

It’s also possible to evaluate the thermal comfort parameters, such as PMV and PPD. For more details, please visit this documentation page.

Last updated: July 12th, 2021

Contents
Data Privacy