ACU-T: 3201 Solar Radiation and Thermal Shell Tutorial

This tutorial introduces you to setting up a CFD simulation involving solar radiation and thermal shells using AcuSolve and HyperWorks CFD. Prior to starting this tutorial, you should have already run through the introductory tutorial, ACU-T: 1000 Basic Flow Set Up, and have a basic understanding of HyperWorks CFD and AcuSolve. To run this simulation, you will need access to a licensed version of HyperWorks CFD and AcuSolve.

Prior to running through this tutorial, click here to download the tutorial models. Extract ACU-T3201_Atrium.x_t and SolarLoad.dat from HyperWorksCFD_tutorial_inputs.zip.

Problem Description

The problem to be addressed is shown schematically in Figure 1. The model consists of an atrium with a couch and chairs in the center. Air flows into the atrium through the inlet vent and exits through the outlet. The front portion of the atrium consists of glass walls supported by an aluminum frame. This aluminum frame will be modeled as a thermal shell; hence, this tutorial introduces you to the process of setting up a transient solar radiation simulation and thermal shells in HyperWorks CFD.


Figure 1.

Solar Radiation Parameters

AcuSolve uses an ideal gray surface solar radiation model to calculate the solar heat flux. The fluxes are computed using a ray trace algorithm and five optical properties of the surface, specular transmissivity ( τ s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiaadohaaeqaaaaa@38DF@ ), diffuse transmissivity ( τ d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiaadsgaaeqaaaaa@38D0@ ), specular reflectivity ( ρ s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aaS baaSqaaiaadohaaeqaaaaa@38DA@ ), diffuse reflectivity ( ρ d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aaS baaSqaaiaadsgaaeqaaaaa@38CB@ ) and absorptivity ( α MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3795@ ).


Figure 2.
A specular transmission occurs when a photon passes straight through a surface with no change of direction. In a diffuse transmission the photon penetrates the surface, but its outgoing energy is uniformly distributed in solid angle over the hemisphere, weighted by projected surface area. For a specular reflection, the angle of reflection is equal to the angle of incidence. Diffuse reflections are similar to diffuse transmissions, except the hemisphere over which the outgoing energy is distributed is on the same side of the surface as the incident photon. Finally, the photon may be absorbed by the surface. These five interactions are associated with five surface properties that together must obey the following constraint:(1) τ s ( θ ) + τ d ( θ ) + ρ s ( θ ) + ρ d ( θ ) + α ( θ ) = 1 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiaadohaaeqaaOGaaiikaiabeI7aXjaacMcacqGHRaWkcqaH epaDdaWgaaWcbaGaamizaaqabaGccaGGOaGaeqiUdeNaaiykaiabgU caRiabeg8aYnaaBaaaleaacaWGZbaabeaakiaacIcacqaH4oqCcaGG PaGaey4kaSIaeqyWdi3aaSbaaSqaaiaadsgaaeqaaOGaaiikaiabeI 7aXjaacMcacqGHRaWkcqaHXoqycaGGOaGaeqiUdeNaaiykaiabg2da 9iaaigdaaaa@57CD@
Where,
τ s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiaadohaaeqaaaaa@38DF@
Specular transmissivity
τ d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiXdq3aaS baaSqaaiaadsgaaeqaaaaa@38D0@
Diffuse transmissivity
ρ s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aaS baaSqaaiaadohaaeqaaaaa@38DA@
Specular reflectivity
ρ d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdi3aaS baaSqaaiaadsgaaeqaaaaa@38CB@
Diffuse reflectivity
α MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3795@
Absorptivity
θ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqiUdehaaa@37AC@
Angle of incidence

For the solar radiative heat fluxes to be computed, a solar radiation surface needs to be defined on that given surface.

In this tutorial, the solar flux loading is given in the form of a data file which was generated using the acuSflux script available in AcuSolve. The script can be used to generate a data file with a four-column array of solar flux vector data values. The piecewise linear type is used in this tutorial to emulate the pattern of sunrise to sunset over the atrium.

For example, to generate the solar load data file for a location with known geological coordinates, enter the following command in the AcuSolve Command Prompt: acuSflux -time "dec-3-2019 11:00:00" -tinc 1800 -nts 25 -lat 42.6064 -lon -83.1498 -ndir "1,0,0" -udir "0,0,1"

Here,
time
The start time in GMT (ex: “dec-3-2019 21:00:00”)
tinc
The time increment in seconds
nts
Number of discrete time steps
lat
Latitude coordinates of the location in degrees North (ex: 45.112 or -37.56 (equal to 37.56 S))
lon
Longitude coordinates of the location in degrees East (ex: 86.26 or -54.84 (equal to 54.84 W))
ndir
The north direction unit vector in model coordinates (should be enclosed in double quotes) (ex: “0,1,0”)
udir
The upward direction unit vector in model coordinates (should be enclosed in double quotes) (ex: “0,1,0”)

Thermal Shell Modeling

The thermal shell in AcuSolve is a feature that creates zero physical thickness volumetric shell elements from surface elements. This is useful when the thickness of the component is too small to be modeled as a solid medium. The thermal shell can have multiple layers, each with different thicknesses and material models. A schematic of the thermal shell is shown below.


Figure 3.

When defining a thermal shell on a surface, two sets of boundary conditions are needed. One for the Primary Wall surface i.e. Shell Inner and one for the Shell Outer Wall surface. In this tutorial, a solar radiation surface will be defined on the outer shell surface so that it receives solar heat flux, whereas the inner shell surface will be modeled as a default wall.

Start HyperWorks CFD and Create the HyperMesh Model Database

  1. Start HyperWorks CFD from the Windows Start menu by clicking Start > Altair <version> > HyperWorks CFD.
    When HyperWorks CFD is loaded, the Geometry ribbon is open by default.
  2. Create a new .hm database in one of the following ways:
    • From the menu bar, click File > Save.
    • From the Home tools, Files tool group, click the Save As tool.


      Figure 4.
  3. In the Save File As dialog, navigate to the directory where you would like to save the database.
  4. Enter Atrium_Solar as the name for the database then click Save.
    This will be your problem directory and all the files related to the simulation will be stored in this location.

Import and Validate the Geometry

Import the Geometry

  1. From the menu bar, click File > Import > Geometry Model.
  2. In the Import File dialog, browse to your working directory then select ACU-T3201_Atrium.x_t and click Open.
  3. In the Geometry Import Options dialog, leave all the default options unchanged then click Import.


    Figure 5.


    Figure 6.

    The model contains an atrium with glass panes supported by an aluminum frame in the front. Air enters from the opening on the roof in the front and exits through the outlet in the rear.

Validate the Geometry

  1. From the Geometry ribbon, click the Validate tool.


    Figure 7.
    The Validate tool scans through the entire model, performs checks on the surfaces and solids, and flags any defects in the geometry, such as free edges, closed shells, intersections, duplicates, and slivers.

    The current model doesn’t have any of the issues mentioned above. Alternatively, if any issues are found, they are indicated by the number in the brackets adjacent to the tool name.

    Observe that a blue check mark appears on the top-left corner of the Validate icon. This indicates that the tool found no issues with the geometry model.


    Figure 8.
  2. Press Esc or right-click in the modeling window and swipe the cursor over the green check mark from right to left.
  3. Save the database.

Set Up Flow

Set Up the Simulation Parameters and Solver Settings

  1. From the Flow ribbon, click the Physics tool.


    Figure 9.
    The Setup dialog opens.
  2. Under the Physics models setting:
    1. Set Time marching to Transient.
    2. Set the Time step size to 1 and the Final time to 30.
    3. Set the Turbulence model to Laminar.
    4. Activate the Heat transfer checkbox.


    Figure 10.
  3. Click the Solver controls setting and set the Maximum stagger iterations to 3.


    Figure 11.
  4. Close the dialog and save the model.

Assign Material Properties

  1. From the Flow ribbon, click the Material tool.


    Figure 12.
  2. Verify that Air is assigned as the material for the fluid domain.
    The legend in the top-left corner of the modeling window lists all the material models assigned to the current model.
    Since this model has a single volume, by default air is assigned as the material for the fluid domain.


    Figure 13.
  3. On the guide bar, click to execute the command and exit the tool.

Define Thin Solids

In this simulation, you will model the aluminum frames as a thin solid.

  1. From the Flow ribbon, click the Thin tool.


    Figure 14.
  2. In the modeling window, select the four surfaces highlighted in the image below.


    Figure 15.
  3. In the microdialog, set the Layer thickness to 0.025 and the Material to Aluminum.
  4. On the guide bar, verify that the number of Parent Surfaces selected is 4 then click to execute the command.
    Once the command is executed, the legend should be updated accordingly to reflect the changes.
  5. Save the model.

Define Flow Boundary Conditions

  1. From the Flow ribbon, click the Constant tool.


    Figure 16.
  2. In the modeling window, click the inlet surface highlighted in the figure below.


    Figure 17.
  3. In the microdialog, enter the values shown below.


    Figure 18.
  4. Click on the guide bar to execute the changes.
  5. Click the Outlet tool.


    Figure 19.
  6. Select the surface highlighted in the figure below then click on the guide bar.


    Figure 20.
  7. Click the No Slip tool.


    Figure 21.
  8. Select all three wall surfaces of the atrium, the roof, and the front glass walls.


    Figure 22. Front-Side of the Model


    Figure 23. Back-Side of the Model

    In total, 21 surfaces should be selected.

  9. In the microdialog, enter the values shown in the figure below.


    Figure 24.
  10. Click . In the new microdialog that appears, set the Nodal Output frequency to 1.


    Figure 25.
  11. On the guide bar, click to execute the command and remain in the tool.
  12. On the guide bar, click the drop-down menu next to Surfaces and change the selection entity to Thin Solids.
    The visibility of all the surfaces except the thin solids is made transparent.
  13. Select all the thin solid surfaces using the window selection method.


    Figure 26.
  14. In the microdialog, enter the values shown in the figure below.


    Figure 27.
  15. Click . In the new microdialog that appears, set the Nodal Output frequency to 1.
  16. On the guide bar, verify that the number of Thin Solids selected is 4 and the Direction is set to Away from Parent Surface, then click .
  17. Change the selection entity on the guide bar back to Surfaces.
  18. From the Boundaries legend, right-click on Default Wall and select Isolate.


    Figure 28.
  19. Select all the surfaces except the aluminum frame, as highlighted in the figure below.


    Figure 29.
  20. In the microdialog, enter the values shown in the figure below.


    Figure 30.
  21. Click . In the new microdialog that appears, set the Nodal Output frequency to 1.
  22. From the Boundaries legend, rename Wall 1 to Floor by double-clicking on it.
  23. Click on the guide bar.
    The updated Boundaries legend should look similar to the one shown below.


    Figure 31.
  24. Turn on the visibility of all the surfaces by right-clicking in the modeling window and selecting Show All or by simply pressing the A key.
  25. Save the model.

Set Up Solar Radiation

Set Up the Solar Radiation Parameters

  1. From the Radiation ribbon, Solar Radiation tools, click the Physics tool.


    Figure 32.
  2. In the Solar Radiation Settings dialog, activate the Solar radiation equation.
  3. Click to load the solar flux input from a file.
  4. In the Open file dialog, set the filter to Dat file (.dat) and select the SolarLoad.dat file provided with the input file for this tutorial.
  5. Click Open.
    The plot in the dialog should look like the one shown in the figure below.


    Figure 33.
  6. Close the dialog and save the model.

Define the Solar Radiation Models

  1. From the Radiation ribbon, Solar Radiation tools, click the Model tool.


    Figure 34.
  2. In the Solar radiation model library, click to add a new solar radiation model.
  3. In the Name column, enter BB out and set the Side to Outward.
  4. Similarly, create the other models and enter the values as shown in the figure below.


    Figure 35.
  5. Close the dialog and save the model.

Assign the Solar Radiation Models

  1. From the Radiation ribbon, click the Surface tool.


    Figure 36.
  2. Select the three wall surfaces, the inlet, and the roof surface, as highlighted in the figures below.


    Figure 37. Front-Side of the Model


    Figure 38. Back-Side of the Model
  3. In the microdialog, set the Solar radiation model to BB out then click on the guide bar.
  4. Select all the glass surfaces shown in the figure below, assign the Glass model to them, then click on the guide bar.


    Figure 39.
  5. Rotate the model and select the floor surface. In the microdialog, assign the BB in model then click on the guide bar.


    Figure 40.
  6. From the Solar Radiation Model legend, right-click on Unassigned and select Isolate.
  7. Select the surfaces of the couch, table, and chairs, assign the BB def model to them, then click on the guide bar.


    Figure 41.
  8. On the guide bar, change the selection entity to Thin Solids.
  9. Using the window selection method, select the four thin solid surfaces and assign the BB out model to them.


    Figure 42.
  10. On the guide bar, verify that the Direction is set to Away from Parent Surface then click to execute the changes.
  11. Turn on the display of all the surfaces and save the model.

Generate the Mesh

In this step, you will define the mesh controls and then generate the mesh.

Define the Surface Mesh Controls

  1. From the Mesh ribbon, click the Surface tool.


    Figure 43.
  2. Using the window selection method, select all the surfaces in the model.
  3. In the microdialog, set the Average element size to 0.15 and the Mesh growth rate to 1.0.


    Figure 44.
  4. On the guide bar, click to execute the command and exit the tool.
  5. Save the model.

Generate the Mesh

  1. From the Mesh ribbon, click the Volume tool.


    Figure 45.
  2. In the Meshing Operations dialog, set the Mesh growth rate to 1.1 then click Mesh to start the meshing process.
    The Run Status dialog opens and the status of the meshing process is shown.
  3. Once the mesh is generated, close the Run Status dialog and save the model.
    Note: Considering the run time of the simulation, a very coarse mesh with no boundary layers is used for this tutorial. Otherwise, a relatively fine mesh with boundary layers to adequately resolve the gradients in the flow and temperature fields should be used.

Compute the Solution

Define the Nodal Output Frequency

  1. From the Solution ribbon, click the Field tool.


    Figure 46.
  2. In the Field Output dialog, activate the check box for Write initial conditions.
  3. Set the Time step interval to 1.


    Figure 47.

Define the Nodal Initial Conditions and Compute the Solution

  1. From the Solution ribbon, click the Run tool.


    Figure 48.
  2. In the Launch AcuSolve dialog, set the Parallel processing option to Intel MPI.
  3. Optional: Set the number of processors to 4 or 8 based on availability.
  4. Deactivate the Automatically define pressure reference option.
  5. Expand the Default initial conditions menu.
  6. Deactivate the Pre-compute flow option.
  7. Set the Temperature to 288.15.
  8. Verify that all the values are set as shown in the figure below.


    Figure 49.
  9. Click Run.
    Once the solution process is started, the Run Status dialog appears.
  10. In the dialog, right-click on the AcuSolve run and select View log file.
    Once the run is complete, a summary of the solution process is shown in the log file.


    Figure 50.

Post-Process the Results with HW-CFD Post

Create Surface Groups

  1. Once the solution is completed, navigate to the Post ribbon.
  2. From the menu bar, click File > Open > Results.
  3. Select the AcuSolve log file in your problem directory to load the results for post-processing.
    The solid and all the surfaces are loaded in the Post Browser.
  4. Click the Boundary Groups tool.


    Figure 51.
  5. Click on the guide bar to open Advanced Selection.
  6. Set the drop-down to By Boundaries then select Thin Solid Wall and Floor.


    Figure 52.
  7. Exit the dialog then click on the guide bar.
    The two boundaries are grouped together under the User Defined surface group in the Post Browser.
  8. Isolate the new surface group and rename it to ThinSolid_Floor.


    Figure 53.

Plot the Temperature Contour

  1. In the Post Browser, right-click on the ThinSolid_Floor surface group and select Edit.
  2. In the display properties microdialog, set the display to temperature and activate the Legend toggle.
  3. Set the legend bounds to 287 and 292.
  4. Click and set the Colormap Name to Rainbow Uniform.


    Figure 54.
  5. Click on the guide bar.
  6. Drag the animation slider at the bottom of the modeling window to the 31st frame.


    Figure 55.

Summary

In this tutorial, you learned how to set up and solve a CFD analysis involving solar radiation. You started by importing a geometry model into HyperWorks CFD and setting up the simulation parameters and boundary conditions. Once you computed the solution, you post-processed the results using the Post ribbon.