WUFI® 2D computes the time-dependent temperature and moisture fields in a two-dimensional cross-section of a building component. Such a two-dimensional computation automatically includes so-called geometrical and structural thermal bridge effects. These are the effects which the component shape (e.g. a corner) and variations of the thermal properties within the component (e.g. reinforcement bars) have on the resulting temperature field. The modifications of the temperature field and the associated heat flows by thermal bridges may have important consequences for energy loss, mold growth in damp corners etc.
WUFI® 2D is not intended to compete with dedicated thermal bridge programs which usually offer more flexible modelling interfaces and provide specific thermal bridge properties. But it can investigate the effect of thermal bridges on energy losses and, in particular, on the hygric conditions in and on building components (mold growth, damage due to condensation etc.), which purely thermal programs can not.
The international standard ISO 10211 provides a series of two- and three-dimensional test cases for validating thermal bridge software. Of course, WUFI® 2D should be able to reproduce the two-dimensional test cases.
Test case 1
Case 1 considers one half of a symmetrical square column with known constant surface temperatures. The steady-state temperature distribution over the cross-section can be computed analytically. 28 temperatures on an equidistant grid are given by the standard as the reference solution; the software to be validated must reproduce these temperatures within 0.1° C.
The given boundary conditions (20° C along the top edge, 0° C along the left and bottom edges, and adiabatic along the right edge) result in a temperature field with strong variation close to the upper left corner and little variation towards the bottom of the cross-section.
Usually, one would therefore create a more efficient variable-size computational grid with fine grid elements close to the upper left corner and progressively larger grid elements towards the bottom. However, since for the present purpose the temperatures must be evaluated at precisely given coordinates and WUFI® computes the temperatures for the centers of the grid elements, a grid has to be created which ensures that a small element is centered at each of the requested positions.
To this end, the monolithic component has been built up from 28 separate blocks which fill up the space between the reference points and are subdivided by relatively coarse grids. The blocks are separated by 4 mm wide gaps which are subdivided by very fine grids in such a way that one small grid element is precisely centered on each reference point (one at each intersection of the gaps).
Since in the present case of a monolithic component with given surface temperatures the steady-state solution for the temperature field does not depend on the thermal properties, any arbitrary material data may be used. The data of concrete were chosen for this exercise.
The prescribed surface temperatures are applied to the component by setting the ambient air to the desired temperatures and the heat transfer coefficients for the surfaces to very large values. For the adiabatic right-hand surface (the symmetry plane which allows to limit the calculation to one half of the original square column), the heat transfer coefficient has been set to zero in order to suppress any heat exchange.
WUFI® 2D has no mode for steady-state solutions, but such a solution is approached to arbitrary precision if a transient computation with constant boundary conditions is performed for a sufficient number of time steps. Here, 10 steps of 48 hours each were found to be sufficient. Moisture transport was switched off for this purely thermal computation.
A graphical postprocessor for the calculation results allows to extract the final temperatures from any grid element. The comparison with the reference temperatures shows that WUFI’s temperatures deviate by 0.05° C or less and are thus well within the allowed deviation of 0.1° C.
You can download a WUFI® 2D project file (65 KB) for the benchmark calculation. Use the “Import…” function of WUFI® 2D to read this compressed archive file.
Test case 2
Case 2 considers the heat flow through a building component which contains materials with widely differing thermal conductivities.
| Dimensions (mm): |
| AB= 500 | | CD = 15 | | EM = 40 | | IM = 1.5 |
| AC = 6 | | CF = 5 | | GJ = 1.5 | | FG-KJ = 1.5 |
| Thermal conductivities (W/mK): |
| 1 (concrete): 1.15 | | 2 (wood): 0.12 | | 3 (insulation): 0.029 | | 4 (aluminum): 230 |
| Boundary conditions: |
| AB: 0° C with Rse = 0.06 m²K/W | | HI: 20° C with Rsi = 0.11 m²K/W | |
The prescribed building component can easily be assembled from rectangles, using WUFI’s graphical component editor. The automatically generated grid with the fineness setting “coarse” is sufficient for this computation.
The thermal conductivities for the four involved materials are specified by the standard and have been entered in WUFI® accordingly. Since WUFI® needs a full set of thermal porperties for each material (including heat capacity etc.) the missing data have been taken from similar materials in WUFI’s material database. The steady-state result only depends on the prescribed thermal conductivities, not on the added properties.
The ambient air temperatures and the heat transfer coefficients for the top and bottom surfaces are entered in WUFI® 2D as specified by the standard. The left and right surfaces are treated as adiabatic by setting the respective heat transfer coefficients to zero.
In this case, too, the steady-state solution must be approximated by a transient calculation with constant boundary conditions. 30 steps of one hour each were found sufficient.
The resulting temperatures at the specified locations can again be extracted with the graphical postprocessor. However, the standard asks for temperatures on material boundaries whereas WUFI® computes the temperatures for the centers of the grid elements and any material boundaries must always coincide with boundaries between grid elements. So in this case it is not possible to center grid elements on the requested locations (such grid elements would contain two or more different materials which is not allowed).
Close to the four corners (points A, B, H and I), the temperature variation is so small that the center of the outermost grid element instead of the true geometric corner can be taken as sufficiently representative.
Where the requested location lies between two materials (points C, E and F), the temperature at this location must be computed from the temperatures at the centers of the two grid elements straddling the location. The temperature ϑm for a location m between locations 1 and 2 can be computed by ϑm = ((λ1/s1) ϑ1 + (λ2/s2) ϑ2) / ((λ1/s1)+(λ2/s2)), where ϑi is the temperature at location i, λi is the thermal conductivity between locations i and m, and si is the distance between locations i and m. Since the grid elements straddling a material boundary all have the same size, the si cancel and the expression reduces to ϑm = (λ1 ϑ1 + λ2 ϑ2) / (λ1+λ2).
For locations where three materials meet (points D and G), the temperature has been computed from the temperatures in the four adjacent grid elements by the following generalisation of the above formula: ϑm = (λ1 ϑ1 + λ2 ϑ2 + λ3 ϑ3 + λ4 ϑ4) / (λ1 + λ2 + λ3 + λ4).
The comparison with the reference temperatures shows that WUFI’s temperatures deviate by 0.1° C or less and are thus within the allowed deviation of 0.1° C. The heat flow through the component is 9.5 W/m and thus within the required (9.5 ± 0.1) W/m.
You can download a WUFI 2D project file (30 KB) for the benchmark calculation. Use the “Import…” function of WUFI® 2D to read this compressed archive file.
Last Update: December 20, 2024 at 17:05
The new version WUFI® Pro 7 is now available for testing in a new look and with new features. We are working hard on the release, which will be available soon. Meanwhile we are offering the opportunity to test our release candidate in advance and try out the new features. 
Version 6.3 of WUFI
There are still places available for the WUFI
The given boundary conditions (
Under unfavorable ambient conditions, microbial growth may occur on the surfaces of building components. The most important factors are temperature, relative humidity, a substrate with sufficient nutrients and the daily duration of the period where all growth conditions prevail simultaneously (coincidence time). While bacteria need relative humidities of at least 90% for being able to grow, certain xerophilic fungi can thrive at humidities as low as 65%, and most fungi can cope with humidities as low as 80%. Mold fungi also tolerate a larger temperature bandwidth than other organisms, they may grow between 0°C and 50°C. Therefore the whole range of humidities and temperatures mentioned above may be considered to pose a potential mold growth risk.
Converting attic space into living space becomes increasingly popular. Since the unvented type of insulation is energetically more favorable and since this type is easier to install if the insulation has to be added to an existing roof, it should be preferred over a vented variant, unless possible moisture problems are a concern. Since older pitched roofs usually have a relatively vapor-tight exterior lining (e.g. bituminous roofing felt on wood sheathing) they require an analysis of possible condensation problems. German standard DIN 4108-3 dispenses with calculational approval if the room-side vapor barrier has a very high vapor resistance (sd-value > 100 m). However, because such constructions which are then vapor-tight on both sides run a considerable risk of severe moisture damage if small defects or leaks occur, it is advisable not to follow the standard here as was already correctly pointed out in [1]. It has also been recommended there to use vapor retarders instead whose diffusion resistance has been selected to limit condensation in winter to an uncritical level while some drying of moisture that has infiltrated in the component is allowed in summer.
The drying of a cellular concrete flat roof is a classic example for the application of modern computational simulation methods. Since traditional methods considering only vapor diffusion were not sufficient to explain the drying behavior of these flat roofs, a new calculation method [1] was introduced in the early 80s which allowed for capillary transport, too, and was thus able to quantitatively reproduce experimental results obtained in the 60s [2]. Since this type of construction has recently attracted the attention of experts because of moisture problems during the drying phase in more southern climates, it will be used here as an example for investigating the hygrothermal behaviour of massive flat roofs.