Performance Evaluation of Wire Cloth Micro Heat Exchangers

The purpose of this study is to validate a thermal-hydraulic simulation model for a new type of heat exchanger for mass, volume, and coolant/refrigerant charge reduction. The new heat exchanger consists of tubes with diameters in the range of 1 m m and wires in the range of 100 m , woven together to form a 200 × 200 × 80 m m 3 wire cloth heat exchanger. Performance of the heat exchanger has been experimentally evaluated using water as inner and air as outer heat transfer medium. A computational thermal and fluid dynamic model has been implemented in OpenFOAM®. The model allows variation of geometry and operating conditions. The validation of the model is based on one single geometry with an opaque fabric and air-side velocities between 1 and 7 m / s . The simulated and measured pressure drops are found to be in good agreement with a relative difference of less than 16%. For the investigated cases, the effective heat transfer coefficients are in very good agreement (less than 5%) when adapting the contact resistance between tubes and wires. The numerical model describes the fluid flow and heat transfer of the tested heat exchanger with adequate precision and can be used for future wire cloth heat exchanger dimensioning for a variety of applications.


Introduction
In recent decades, heat exchanger design and development focused on various aspects in order to improve heat transfer processes. In general, the improvements are related to (i) pumping or fan power reduction, (ii) volume reduction, and (iii) mass reduction, with an increase in (iv) heat transfer surface area and (v) in heat transfer coefficients.
Extended surface plate-fin and tube-fin heat exchangers experienced enhanced miniaturization to achieve the mentioned improvements. A reduction in pressure drop can be achieved by increasing the frontal face area and decreasing the flow length. Volume reduction is presently essential for charge reduction of new and natural refrigerants. Mass reduction is realized by decreasing tube and fin thicknesses down to 100 µm and less. A heat transfer coefficient and surface area increase is achieved by decreasing the fin pitches to dimensions around or below 1 mm and by geometrical modification of the fins (e.g., louvered, wavy, or offset strip fins).
In order to significantly increase heat transfer coefficients, wire structures are being investigated as heat transfer enhancements. Metallic woven-wire mesh structures [1][2][3][4][5], pin fin structures [6,7] and screen-fin structures [8] were designed for flat primary surfaces. One air-to-air heat exchanger is manufactured and marketed by Vision4Energy [9], numerical performance evaluations of metallic woven-wire mesh structures were presented in [10,11]. Kumra et al. and Lee et al. [12,13] analyzed wire-on-tube type heat exchangers with wires being welded on top of both sides of a heat transfer tube.
Previous studies concentrated either on woven-wire meshes or single wires with diameter in the order of 1 mm. The first design lacks in good thermal conductivity within the wire mesh due to insufficient thermal contact between the wires; the second design does not show high heat transfer coefficients on the wires due to the relatively large wire diameter. Attempts to decrease the wire diameter significantly for the second design are effortful [6,7]. New designs are needed, which make use of the fast and well developed weaving process, without having wires in the mesh which do not contribute to the heat transfer process.
Such possible design concepts are detailed in [14,15]. The wire diameters in these (patented) heat exchangers primarily range from 100 µm to 1 mm. The heat exchanger design, analyzed in this paper, follows these designs and uses a wire diameter of 200 µm. In the manufacturing process of the wire cloth the micro tubes (hollow wires) and the wire structure (solid wires) are woven together on an adopted loom, similar to the weaving process known from textile manufacturing [16]: the tubes are processed as weft threads, the wires as warp threads in a standard linen weave. The large-scale wire cloth is thereafter cut to size. The micro tubes can be manufactured with diameters well below 2 mm [17] and wall thicknesses below 100 µm. The resulting geometry provides small volume specific masses and high volume specific surface areas on both fluid sides. In addition, the pressure stability of the micro tubes is ensured at high pressures. A variety of applications are addressed with the wire cloth micro heat exchanger. These include applications such as heating and cooling of automotive as well as air conditioning of buildings. The heat exchanger might be beneficial when volume in air flow direction is strongly limited (e.g., in engine compartments or building walls).
In this paper, we present thermal and flow measurements of a sample heat exchanger and computational results from a fluid dynamic simulation of the air-side flow through the heat exchanger coupled with a heat transfer simulation.

Methods
Performance of the heat exchanger has been determined by simulation and by experimental analysis. Therefore several samples of wire fabrics have been manufactured. A schematic of a section of the wire cloth heat exchanger is shown in Figure 1.

Test Hardware
In this study only one of the fabrics is analyzed. The sample, as implemented into the heat exchanger, is shown in Figure 2. The heat exchanger properties are listed in Table 1. The main difficulty in manufacturing the wire cloth was wire breakage due to the large size differences of tubes and wires and the resulting mechanical stress in the wire; especially for aluminium wires. The second challenge was to roll up the wire cloth continuously after the weaving process without displacement of the wires on the tubes. Again due to the size differences between tubes and wires, the rolling process for the wire cloth differed from standard fabrics. The third difficulty was to provide sufficient mechanical contact between the tubes and the wires through the contacting process. The contacting was realized by diffusion bonding for the stainless steel variant. Lastly, the insertion of the tubes into the header/collector plates was very time consuming. In particular, burs on the tube endings and non-uniform tube distances had to be dealt with.  the standard deviation is 0.07 mm. Thus, a strong influence on the fluid flow due to geometrical irregularities is not expected.

Analysis Model
Within the analysis model, three different regions are distinguished; (i) the fluid flow within the tubes, (ii) the structure itself, including tubes and wires, and (iii) the gas flow outside of the structure. Flow simulations are carried out, using a finite-volume, multi domain solver for conjugated heat  Table 1. the standard deviation is 0.07 mm. Thus, a strong influence on the fluid flow due to geometrical irregularities is not expected.

Analysis Model
Within the analysis model, three different regions are distinguished; (i) the fluid flow within the tubes, (ii) the structure itself, including tubes and wires, and (iii) the gas flow outside of the structure. Flow simulations are carried out, using a finite-volume, multi domain solver for conjugated heat  The standard deviation of the wire pitch t 1 of the final wire fabric is 5 µm; analyzed via microscope images of five different areas on the wire fabric, each area representing one tube and more than 30 wires. Maximum deviations from the target value of t 1 = 200 µm is 10 µm. For the tube pitch t 2 = 3.5 mm the standard deviation is 0.07 mm. Thus, a strong influence on the fluid flow due to geometrical irregularities is not expected.

Analysis Model
Within the analysis model, three different regions are distinguished; (i) the fluid flow within the tubes, (ii) the structure itself, including tubes and wires, and (iii) the gas flow outside of the structure. Flow simulations are carried out, using a finite-volume, multi domain solver for conjugated heat transfer processes available in OpenFOAM R . For efficiency and simplification purpose, the fluid flow in the pipes region is modeled using correlations for pressure loss (Euler number) and heat transfer (Nusselt number) available from literature [18].

Governing Equations and Assumptions
Heat transfer within the solid phase of the structure is described by a simple heat diffusion equation, i.e., The solution to this equation and the relevant boundary conditions at the material surfaces is the time-dependent temperature distribution in the fabric, which depends on the thermal conductivity k and heat capacity c p of the solid. The governing equations for the gas phase are the three-dimensional Navier-Stokes equations, which weresimplified by assuming laminar, incompressible flow of a Newtonian fluid with constant fluid properties and neglecting volumetric forces and radiation heat transfer. Accordingly, the equations for conservation of mass, momentum and energy are given by and where u, τ, p, ρ and h denote the velocity vector, shear stress tensor, static pressure, gas density and gas enthalpy, respectively. Due to the rather low flow velocity, which result in specific Reynolds numbers lower than 1000, the fluid flow is suspected to be purely laminar inside of the fabric zone. Random fluctuations in the upstream flow would be reduced by the rapid acceleration of the flow when entering the fabric and also by the high shear stress in the fabric. For the case of a turbulent flow entering the wire-cloth heat exchanger, air-side heat transfer coefficients can be expected to be larger.

Simulation Domain and Boundary Conditions
To make the numerical analysis manageable, only a small geometrically periodic section of the wire cloth was analyzed, as shown in Figure 3. Effects due to flow and geometric non-uniformity near the distribution plates and the heat exchanger frame are not captured by the analysis. The analyzed domain consists of 2 half-wires and three (half-)tubes with the direction of the gas flow being perpendicular to the plane in which the tube axes lie, i.e. in positive x-direction according to Figure 3. As shown in Figure 3, the gas flow enters the domain at X 0 with fixed velocity and temperature, at X 1 the outlet pressure is fixed, while temperature and velocity gradients are set to zero. In y-direction, cyclic boundary conditions (BC) are applied at Y 0 and Y 1 , while in z-direction (Z 0 , Z 1 ) wall slip/symmetry BCs are applied for velocity and zero gradient for temperature. In the solid domain the temperature at the inner walls of the tubes (I) is fixed. In other words, the heat transfer phenomenon on the liquid side within the tube has been excluded from the numerical analysis. At the interface (W) between solid and gas phase, a special coupling BC is applied by solving Equation (5), ensuring energy conservation across the interface.
Lastly, a solid-solid interface (K) is introduced between the tubes and wires in order to account for the contact resistance between them. At the position of the red contact, shown in Figure 3, a heat flow dependent temperature difference results, as a consequence of assuming a wall of defined thermal conductivity and defined width and depth. the prescribed contact resistance is a result of the manufacturing process, where the wire cloth structure is mechanically and thermally enhanced by diffusion welding or soldering the tubes and the wires together. Ideally the effect of such measures is a very low thermal resistance; but due to deficiencies in the new manufacturing process of the thermal fabric, the resistance is quite substantial and therefore has to be taken into account in the simulation model.  near the distribution plates and the heat exchanger frame are not captured by the analysis. The analyzed domain consists of 2 half-wires and three (half-)tubes with the direction of the gas flow being perpendicular to the plane in which the tube axes lie, i.e. in positive x-direction according to Figure 3. As shown in Figure 3, the gas flow enters the domain at X 0 with fixed velocity and temperature, at X 1 the outlet pressure is fixed, while temperature and velocity gradients are set to zero. In y-direction, cyclic boundary conditions (BC) are applied at Y 0 and Y 1 , while in z-direction (Z 0 , Z 1 ) wall slip / symmetry BCs are applied for velocity and zero gradient for temperature. In the solid domain the temperature at the inner walls of the tubes (I) is fixed. In other words, the heat transfer phenomenon on the liquid side within the tube has been excluded from the numerical analysis. At the interface (W) between solid and gas phase, a special coupling BC is applied by solving equation (5), ensuring energy conservation across the interface.
Lastly, a solid-solid interface (K) is introduced between the tubes and wires in order to account for the contact resistance between them. At the position of the red contact, shown in Figure 3, a heat flow dependent temperature difference results, as a consequence of assuming a wall of defined thermal conductivity and defined width and depth. the prescribed contact resistance is a result of the manufacturing process, where the wire cloth structure is mechanically and thermally enhanced by For the simulations the Grid Convergence Method based on the Richardson extrapolation was used to determine the numerical uncertainty. Two key variables were defined: (i) the air-side pressure drop and (ii) the air-side outlet temperature. For all simulations, including additional geometries, the grid was refined such that the numerical uncertainty (fine-grid convergence index) for the key variables was less than 5%.Thus, uncertainties related to the choice of grid are expected to be small. More information on the meshing, refinement and uncertainties can be found in Martens ([19], p. 37-43).

One-Way Coupling between Momentum-and Energy-Equation
The solution of the differential equations in the two domains is performed in two consecutive steps. In a first step, only the unsteady mass-and momentum conservation equations for the gas phase (Equations (2) and (3)) are solved assuming laminar incompressible flow. From this solution, the average pressure loss ∆p across the wire cloth is calculated, and the temporal mean value of the velocity field u is stored. While velocity fluctuations affect the average pressure loss, their effect on heat transfer is found to be small. This is partly due to the large difference in heat capacities between the solid and the gas phase, which causes the temperature in the solid to remain almost constant over time, even at pronounced transient operating conditions, for example with vortex shedding in the gas phase behind the tubes and wires. The predicted separation behavior of the mean flow was found to result in an adequate prediction of the mean heat transfer. In a second step, the energy equations for the fluid (Equation (4)) and solid (Equation (1)) domains are solved on the basis of the time-averaged velocity field calculated in the first step. By decoupling the momentum and energy equations, it is possible to change properties of the solid phase and solve the heat transfer problem to good accuracy without having to recalculate the velocity field solution again. Figure 4 shows an example of the predicted gas phase velocity field. In the left picture the distribution of the velocity u(t) at a specific time t during the simulation. Upstream and within the wire cloth structure, the velocity distribution appears symmetrical. Downstream, irregular flow structures are found due to vortex shedding. The calculated time-averaged velocity u is displayed in the middle picture of Figure 4 . The velocity distribution upstream and inside the wire cloth is virtually identical to the instantaneous velocity distribution. Downstream, behind the tubes, large symmetric flow separation zones are found. The standard deviation of the velocityσ(u) (right image in Figure 4) shows that significant transient structures are only present in the wake and thus have little effect on the gas-to-solid heat transfer. small. More information on the meshing, refinement and uncertainties can be found in Martens [19, p. 37-43].

One-Way Coupling between Momentum-and Energy-Equation
The solution of the differential equations in the two domains is performed in two consecutive steps. In a first step, only the unsteady mass-and momentum conservation equations for the gas phase ( (2) and (3)) are solved assuming laminar incompressible flow. From this solution, the average pressure loss ∆p across the wire cloth is calculated, and the temporal mean value of the velocity field u is stored. While velocity fluctuations affect the average pressure loss, their effect on heat transfer is found to be small. This is partly due to the large difference in heat capacities between the solid and the gas phase, which causes the temperature in the solid to remain almost constant over time, even at pronounced transient operating conditions, for example with vortex shedding in the gas phase behind the tubes and wires. The predicted separation behavior of the mean flow was found to result in an adequate prediction of the mean heat transfer. In a second step, the energy equations for the fluid (4) and solid In order to validate (a) the symmetry and periodicity assumptions imposed onto the selected computational domain and (b) the uncoupling of momentum and energy equations, several simulation runs on a larger domain were carried out with and without coupled solutions of Equations (1)-(5). The larger domain consisted of 5 times 10 symmetry units in y-and z-direction, respectively, in order to enable vortex shedding in a larger unconstrained scale. By using large diameters for tubes and wires in addition to wide spacing and a high velocity we achieve large Reynolds numbers, which results in strongly transient flow behaviour. By using large diameter tubes and wires in addition to wider tube and wire spacings, high air-side inflow velocities occur and result in larger Reynolds numbers and potentially locally-turbulent flow through the cloth. Therefore the investigate scenario together with the employed assumptions represents a worst-case scenario with regards to the heat transfer capabilities of the system. Table 2 summarizes the results of the described validation runs. With regards to the reference case, with large domain and coupled solution, the uncoupled solution differs by 5.3% in the air-side heat transfer coefficient h air , while the pressure loss is the same for both cases, because fluid temperature change due to the coupling is small for the investigated case. On the smaller domain, the pressure loss is found to be 3.6% higher than on the larger domain, while the heat transfer coefficient is lower in the coupled and uncoupled solution (2.2% and 2.1%, respectively). Overall, the differences in the solutions are small, thus justifying the reduction of the computational domain (to one symmetry unit), the imposed BCs and the decoupling of the momentum and energy equations for the purpose of estimating the thermal performance of the wire-cloth.

Experimental Set-Up
The thermal-hydraulic performance of the wire-cloth was evaluated at a test rig at the Fraunhofer ISE. Air and water are used as working fluids. The test rig is described in detail in [3]. The specific measurement procedure for the wire cloth micro heat exchanger shall be explained briefly based on the schematic in Figure 5  In the test section, fast-response thin-film Pt100 temperature sensors are used. Four sensors are installed at the inlet and 12 at the outlet channel, respectively. The pressure drop across the sample is measured with two differential pressure transmitters with a measuring range of 0-250 Pa and 0-1000 Pa, respectively. The entire testing section was leakage tested and thermally insulated. On the water-side, fluid temperature and volume flow rate is controlled by a heater unit and set to 85 • C and 1.7 m 3 /h. Here, Pt100 resistance temperature sensors and an electromagnetic flow sensor are employed. Pressure drop is measured via a differential pressure transmitter. For each measurement point, water and air side temperatures and velocities are kept constant for approximately 40 min prior to data collection, in order to assure steady conditions. After this time interval all quantities are averaged over a shorter time period to obtain a single performance value for each measurement point. More details on the major sensor technology is given in the Appendix B.

Performance Parameters
The heat transfer rateQ HX of the air to water heat exchanger required to calculate air-side performance characteristics can be expressed using three different equations, i.e., Q HX =ṁ air c p,air ∆T air , Q HX =ṁ water c p,water ∆T water , Q HX = U HX A HTS,air ∆T m .
Here, the mass flow rates on the air-sideṁ air and water-sideṁ water and the temperature changes within the fluids across the heat exchanger ∆T air and ∆T water are determined experimentally. Thus Equation (6) or (7) can be used to calculate the heat transfer rateQ HX from the measured data. For the purpose of calculating heat transfer from simulation data, Equation (6) is used. The overall heat transfer coefficient U HX in Equation (8) relates to the total heat transfer surface area A HTS,air on the air side and the true (or effective) mean temperature difference ∆T m [20], i.e., The correction factor F in the above equation can be determined based on the flow configuration and the operating conditions and the well-known logarithmic mean temperature difference (LMTD) ∆T log which is given by Therefore a combination of either Equation (6) or (7) with Equation (8) yields U HX . The air-side effective heat transfer coefficient U air,eff can be obtained from The second term on the right-hand-side of Equation (11) indicates the water-side thermal resistance, with h water known from previous experimental work for the experimental data (cf. Appendix B); The second term is set to zero for simulation data with T water = T water,in = T water,out used in Equation (10). The third term on the right-hand-side of Equation (11) indicates the tube wall thermal resistance for thin walled tubes. All terms on the right hand side of Equation (11) can be determined from experimental or simulation data. The effective heat transfer coefficient U air,eff includes the conductive heat transfer through the wire and the convective heat transfer from the tube and the wire to the air. In general, U air,eff can be expressed by the product of extended surface efficiency η 0 and convection heat transfer coefficient h air . The efficiency η 0 is given by where η f represents the fin efficiency of the wire structure and A HTS,air = A wall,air + A wire,air .
A separation into the two factors η 0 and h air requires a valid calculation method for the efficienc η f . Such a method is only available for well defined fin types [20]. For a wire cloth, several simplifying assumptions would be necessary to make use of the given methods. Thus, the experimental thermal performance data will be compared to simulation data by using the dimensional effective heat transfer coefficient U air,eff . To match the tested hardware a diffusion bonding zone between wires and tubes has been introduced in the simulation model. As an estimation of the bonding quality was not readily available, a parameter study on the thermal resistances R of the bonding zone has been performed. The thermal behavior of the diffusion-bonding zone has been modelled like a wall as shown in Figure 6. The relevant properties of the wall are the geometry (height h b , length l b and thickness w b ) and its thermal conductivity k b . Equation (14) defines the effective conductivity of the bonding zone: The length of the wall l b (cf. Figure A1) is derived from the geometry of the fabric and is equal to the bent part of the wire: The height of the diffusion-bonding zone h b has been discussed in [21][22][23]. Following the literature a value of 10 µm has been used here. Height and width are coupled via geometric relations: The resistance value R used in the numerical simulations has been calibrated to match heat transfer rates measured for one operating point, see Figure 7b. Once R has been determined, the effective mean thermal conductivity of the bonding zone k b in the fabric can be calculated from Equation (14). See also Yovanovitch et al. [24] in this context. They describe the phenomenon of higher thermal resistance at soldered joints. Aside from benchmarking thermal performance, as a consequence of (decoupled) fluid flow behavior, accuracy of the simulated velocity and pressure fields has been assessed by comparing air-side mean pressure drop ∆p air with experimental data for a range of inflow condition at fixed outflow temperature and outflow static pressure.

Results
The thermal-hydraulic performance of the wire cloth heat exchanger is addressed in this section. Effective heat transfer coefficient and pressure drop are used to show the performance. The uncertainty of performance data is expressed first as a combined standard uncertainty based on measurement data uncertainty in Appendix B and propagation of uncertainty. Second, an expanded uncertainty is obtained by multiplying the combined standard uncertainty by a coverage factor of 2. Assuming that the measurement data are normally distributed, 95% of the data lies in this interval of expanded uncertainty. The expanded uncertainty for U air,eff is maximum 8 W/(m 2 K); for ∆p air it is maximum 1.2 Pa. The energy balance for these measurements is checked by comparing Equation (6) with Equation (7). The values ofQ HX have a relative difference of 1% to maximum 8%.
Thermal performance as a function of contact resistance R between wires and tubes has been analyzed. In the simulations the contact resistance R has been varied between 10 and 10 5 K/W. The best agreement between measured and simulated effective heat transfer coefficient U air,eff was obtained for a contact resistance of R = 625 K/W. Figure 7a shows U air,eff as a function of inlet air velocity u for different contact resistances R. U air,eff increases monotonously with increasing velocity. Figure 7b depicts U air,eff as a function of reciprocal contact resistances for a fixed air velocity u. Minimum and maximum values for U air,eff are up to 40% lower and 35% higher than the measured values, within the considered range or R.
The contact resistance R was selected so that experimental and numerical values of U air,eff match for the analyzed case, i.e., R = 625 K/W (equivalent to R −1 = 1.6 × 10 −3 W/K). This yields k b = 0.26 W/(m K) based on Equation (14) and Equations (A1) and (A3) of the appendix. This value for the thermal conductivity is higher than the maximum value expected for a soldered joint [24], but it is still one to two orders of magnitude lower in comparison to a homogeneous solid body of the same material.

Discussion
It can be shown that wire cloth micro heat exchangers can be manufactured with a well-developed weaving process and that thermal contact between wires and tubes is acceptable, albeit not ideal. Thus the disadvantages of a wire mesh with wires not contributing to the heat transfer process (e.g., [1-5]) or wires with low heat transfer coefficients due to large wire diameter (e.g., [12,13]) can be overcome.
Predicted pressure drop values are in good agreement with the measurements. For an upstream air velocity of 6 m/s the relative difference between calculated and measured values is 16%. This accuracy might be sufficient for preliminary design of sizing considerations. The size of the pressure drop is in a reasonable range for an opaque wire cloth; however, for a heat exchanger of this specific size the values are too high. The pressure drop at 2 m/s per length of heat exchanger in air flow direction for very compact louvered fin heat exchangers (β > 2000 m 2 /m 3 ) can be in the range of 2000 to 3500 Pa/m [6]. For measured wire cloth the value was 23,000 Pa/m. However, the produced sample heat exchanger was not designed to provide good thermal-hydraulic performance, but to allow evaluation of the wire cloth manufacturing process and validation of the analysis model. Less dense wire cloths promise a much better hydraulic performance.
Good agreement between predicted and measured thermal performance can be achieved by properly calculating the value of thermal contact resistance R between wires and tubes due to the employed diffusion-bonding process. The value of the contact resistance R is in a reasonable range according to [24]. However, we were not able to verify its value by means of another measurement method. It remains to be shown whether the assumptions employed in the simulation model (such as laminar flow through the wire structure, uniform inflow velocity, uniform wire and tube distances) are adequate or their influence on the result is balanced by the fitted value for the thermal resistance R. Additional heat exchanger samples with different geometric parameters and a different material (aluminium) will be available for analysis in the near future. Thus, a matching of thermal resistance values can be done. The main challenge for the manufacturing process will be to produce wire cloths with consistent properties and with very good thermal contact between tubes and wires.
The heat transfer on the air side can be separated in terms of heat transfer on the primary surface (tubes) and on the secondary surface (wires). For the flow around the 2 mm tubes, Reynolds numbers (based on the tube diameter) from 100 to 1000 occur, with eddies shedding from each side of the tube, forming rows of vortices in its wake. This phenomenon allows a mixing of air passing the wires in a main flow area and air in the wake of the tube. Thus, heat transfer on the primary surface can contribute significantly to the overall heat transfer. Apart from that, the secondary surface contributes significantly to the heat transfer. Reynolds numbers (based on wire diameter) range from 10 to 100. Thus, a steady closed wake flow regime and a laminar vortex shedding regime occur for low and high velocities, respectively. The wire diameter is by a factor of 10 smaller than the tube diameter. This has the effect of a smaller thermal boundary layer on the wire surface and thus a higher local heat transfer coefficient on the wires.
Values for the effective heat transfer coefficient U air,eff and volumetric surface area density β of the measured wire cloth heat exchanger are at the upper range of compact heat exchangers. In comparison to louvered fins the product of U air,eff and β is more than 15% to 30% higher for the wire cloth analyzed here [6]. It is assumed that this value will decrease when considering more open structures for the purpose of reducing pressure drop. However, overall dimensions could be preserved by means of different wire arrangements (e.g., twill weave), meandering arrangement of the wire cloth and/or further reduction of tube and wire sizes.

Conclusions
The main conclusions of this experimental and numerical work on wire cloth micro heat exchangers are as follows: 1. Simulation results showed reasonable agreement with measured data for pressure drop across the wire cloth with relative differences below 16%. 2. The thermal contact between tubes and wires could not be determined independently, thus the thermal contact resistance was fitted with measured data. The fitted value is in good agreement with the values reported in literature. 3. Simulated and measured values for the efficient heat transfer coefficient data are in very good agreement based on the fitted thermal contact resistance. 4. Values for the efficient heat transfer coefficient are strongly dependent on the thermal contact resistance between tubes and wires. A very good contact could yield up to 40% higher values for the efficient heat transfer coefficient compared to the investigated sample (cf. Figure 7b). 5. The produced and tested sample heat exchanger was not intended to yield a good thermal-hydraulic performance, but did allow evaluation of the employed manufacturing wire cloth process and benchmarking of the simulation model. Therefore, the realised performance is only to a limited extent comparable with other types of heat exchangers. Specifically: • The achieved effective heat transfer coefficient of 180 W K/m 2 at u = 2 m/s is in the range of very compact louvered fin heat exchangers, but not superior. • The pressure drop of the wire cloth heat exchanger sample is very high (by a factor of 7 to 11) compared to louvered fin heat exchangers and with respect to the provided heat transfer surface area.
6. Wire cloth micro heat exchangers provide flexibility in heat exchanger construction. They allow, for example, curved shapes as well as different spacings of the tubes within one heat exchanger.
Correlations for pressure drop and heat transfer characteristics have been developed for a variety of geometrical parameters based on the present simulation model in [19]. These correlations facilitate the design of wire cloth micro heat exchangers for specific applications. Future work on the heat exchanger should focus on the thermal contact between wires and tubes and on more open wire cloth to decrease the pressure drop significantly. Moreover, the possibility to bend the planar wire cloth easily allows considering different shapes of heat exchangers that could be interesting from a design and construction point of view.