Hydraulic Diameter Research Articles

An experimental investigation of boundary layer over permeable interfaces in Hele-Shaw micromodels

This study experimentally investigates boundary layer development over permeable interfaces using Hele-Shaw micromodels and high-resolution micro-particle image velocimetry (micro-PIV). Velocity vectors, captured at a 5 μm scale, reveal the flow behavior at the interface between free-flow and porous media with ordered structures and porosities ranging from 50% to 85%. The results show that the boundary layer streamline alignment decreases with increasing porosity, while lower permeability fosters more uniform and parallel flow near the interface. Flow channeling occurs along paths of the least resistance, with more flow directed through the Hele-Shaw free-flow region as the solid fraction of the porous material increases. The Reynolds number (0.14–0.94), based on the Hele-Shaw hydraulic diameter, has a minimal effect on the normalized velocity distribution. Furthermore, an analytical solution for the external boundary layer thickness exhibited good agreement with experimental data, confirming a thickness of 2–4 times the square root of the free-flow Hele-Shaw permeability. Additionally, a Q-criterion analysis identified, for the first time, distinct zones within the external boundary layer, capturing the balance between rotational and deformation components as a function of permeability. These findings offer insight into flow dynamics in porous media systems, with implications for both natural and industrial applications, and contribute to the improved modeling of fluid dynamics and momentum transport in coupled free-flow and porous media environments.

Boiling heat transfer of water on smooth and square pin fins surfaces in minichannel

In this investigation, flow boiling heat transfer and hydrodynamic instability using deionized water in a minichannel with smooth and square pin fin surfaces was experimentally examined. The study is aimed to assess local and overall heat transfer coefficients, heat transfer coefficients, pressure drop and Nusselt number across various liquid mass flux ranges (50–93 kg m−2 s−1) and a constant heat flux of 46.91 kW/m2. The minichannel, with a hydraulic diameter of 2.85 mm, features a rectangular shape. Experimental testing was conducted on aluminum surface within a minichannel of dimensions 270 mm in length, 30 mm in width, and 1.5 mm height. The square pin fins are integrated at the base of the boiling surfaces, arranged in staggered form. The results revealed that the boiling heat transfer coefficient on square pin fin surfaces was approximately 88 % higher than that on smooth surfaces. Additionally, a notable increase in local heat transfer coefficient was observed for square pin fin surfaces (i.e., 34.64 kW/m2.K) compared to smooth surfaces (8.09 kW/m2.K) at the same mass flux of 50 kg/m2.s. Besides, the square pin fins surface exhibited a lower and stable pressure drop, particularly near annular flow regimes, compared to the smooth surface. Notably, the highest pressure drop observed was 4.8 kPa for the square pin fins surface, while the smooth surface experienced a higher pressure drop of 6.36 kPa, accompanied by more pronounced pressure fluctuations during annular flow.

Unravelling Different Water Management Strategies in Three Olive Cultivars: The Role of Osmoprotectants, Proteins, and Wood Properties

Understanding the responses of olive trees to drought stress is crucial for improving cultivation and developing drought-tolerant varieties. Water transport and storage within the plant is a key factor in drought-tolerance strategies. Water management can be based on a variety of factors such as stomatal control, osmoprotectant molecules, proteins and wood properties. The aim of the study was to evaluate the water management strategy under drought stress from an anatomical and biochemical point of view in three young Italian olive cultivars (Giarraffa, Leccino and Maurino) previously distinguished for their physiological and metabolomic responses. For each cultivar, 15 individuals in pots were exposed or not to 28 days of water withholding. Every 7 days, the content of sugars (including mannitol), proline, aquaporins, osmotins, and dehydrins, in leaves and stems, as well as the chemical and anatomical characteristics of the wood of the three cultivars, were analyzed. ‘Giarraffa’ reduced glucose levels and increased mannitol production, while ‘Leccino’ accumulated more proline. Both ‘Leccino’ and ‘Maurino’ increased sucrose and aquaporin levels, possibly due to their ability to remove embolisms. ‘Maurino’ and ‘Leccino’ accumulated more dehydrins and osmotins. While neither genotype nor stress affected wood chemistry, ‘Maurino’ had a higher vessel-to-xylem area ratio and a larger hydraulic diameter, which allows it to maintain a high transpiration rate but may make it more susceptible to cavitation. The results emphasized the need for an integrated approach, highlighting the importance of the relative timing and sequence of each parameter analyzed, allowing, overall, to define a “strategy” rather than a “response” to drought of each cultivar.

A platinum-coated staggered reactor to intensify lean hydrogen/air combustion: A large eddy simulation study

Catalytic-aided combustion has been proven effective for premixed hydrogen/air mixtures, particularly under lean to ultra-lean conditions. However, minimising the required catalyst sets a significant challenge because noble metals with high catalytic activity are rare and expensive. Therefore, this study aims to intensify the catalytic combustion process by investigating a non-planar reactor comprising an array of platinum-coated half- and full-cylinders through large eddy simulation. A premixed mixture with a fuel-lean equivalence ratio of 0.15 and an incoming Reynolds number of 3500 based on hydraulic diameter is used. For comparison, a planar reactor without cylinders is also studied under the same operating conditions and with the same amount of platinum-coated surface area. The simulation employs the turbulent kinetic energy sub-grid model and the eddy dissipation concept to model the turbulent catalytic reacting flow. The discrete ordinate model is used to account for radiation heat transfer in the catalytic process. Numerical simulations are validated against experimental results prior to analysis. The findings indicate that the placement of cylinders along the reactor length enhances convective mass transfer and intensifies catalytic combustion, resulting in effective combustion over a smaller catalytic surface. Compared to planar models, non-planar reactors demonstrate a much better H2 conversion efficiency throughout the reactor length, saving nearly 62.5 % of the catalyst.

Comparative Analysis of Modular Pin-Fin Heat Sink Performance: Influence of Geometric Variation on Heat Transfer Characteristics

The advancement of digital technology in the age of Industrial Revolution 4.0 has driven other engineering sectors to keep up with the progressive demand for high-performance computing and electronics. Thermal management plays a critical role in maintaining the performance of electronic devices by keeping the system temperature at an optimal level and preventing overheating. A heat dissipation device widely employed in computing and other electronic systems’ heat management is a fin-type heat sink, owing to its robust and simple design. In this work, the authors evaluated heat sinks of different pin-fin cross-sectional geometry, arranged in a staggered configuration, in terms of heat transfer characteristics under forced convection. Three basic geometries were chosen on the grounds of manufacturing easiness and market availability: circular, square, and hexagonal. All the tested geometries had approximately the same hydraulic diameter. Therefore, the results could be compared to each other. The investigation revealed that the square cross-section induced an excellent convection coefficient, hence higher heat dissipation, compared to the counterparts. The tradeoff between heat transfer performance and size-related parameters, such as surface area and material volume, is also discussed in this paper.

Impacts of Material and Machine on the Variation of Additively Manufactured Cooling Channels

Abstract While additive manufacturing (AM) can reduce component development time and create unique internal cooling designs, the AM process also introduces several sources of variability, such as the selection of machine, material, and print parameters. Because of these sources, wide variations in a part's geometrical accuracy and surface roughness levels can occur, especially for small internal cooling features that are difficult to post-process. This study investigates how the selection of machine and material in the AM process influences variations in surface quality and deviations from the design intent. Two microscale cooling geometries were tested: wavy channels and diamond-shaped pin fins. Test coupons were fabricated with five different additive machines and four materials using process parameters recommended by the manufacturers. The as-built geometry was measured non-destructively with computed tomography scans. To evaluate surface roughness, the coupons were cut open and examined using a laser microscope. Three distinct roughness profiles on the coupon surfaces were captured including upskin, downskin, and channel walls built at 90 deg to the build plate. Results indicated that both material and machine contribute to producing different roughness levels and very different surface morphologies. The roughness levels on the downskin surfaces are significantly greater than on the upskin or sidewall surfaces. Geometric analysis revealed that while the hydraulic diameter of all coupons was well captured, the pin cross section varied considerably. Along with characterizing the coupon surfaces, cooling performance was investigated by experimentally measuring friction factor and heat transfer. The variations in surface morphology as a function of material and machine resulted in heat transfer fluctuating by up to 50% between coupons featuring wavy channels and 26% for coupons with pin fin arrays. Increased arithmetic mean surface roughness led to increased heat transfer and pressure drop; however, a secondary driver in the performance of the wavy channels was found to be the roughness morphology, which could be described using the surface skewness and kurtosis.

Thermal performance of a novel hybrid heat pipe with composite heat transfer characteristics

Considering the capillary threshold of oscillating heat pipes (OHP), a novel hybrid heat pipe (HHP) was designed with channels with hydraulic diameters that alternately exceed and fall below this threshold, parallelly arrayed within a unified copper substrate. This technology attains a comprehensive heat transfer performance within the HHP, seamlessly transitioning from pure phase change heat transfer to hybrid pulsating flow heat transfer to effectively accommodate the escalating demands of heat load. Performance and visualization are used to analyze the operational characteristics of HHP. The results indicate that the HHP mainly transfers heat through the phase change in the large channels in lower heat input, while the working fluid in the small channels remains calm. Once the heat input is sufficient, the geyser boiling in the large channels and the liquid slug pulsating in the small channels jointly carry out heat transfer. Furthermore, the liquid content in the large and small channels can be adaptively adjusted to each other. The high-speed liquid in the small channels can be squeezed into the large channels to ensure the sufficient pulsating space and achieve efficient heat transfer under high heat input. In addition, the operation of the HHP can be significantly impacted by the filling ratio (FR). In the phase change stage, the FR of 47.6 % can be considered as the turning point. After transitioning to the hybrid pulsating flow stage, the HHP with low FR can get a better performance. This new type of heat pipe with gradually changing operation characteristics may be critical to solve thermal problems in variable heat load domain.

On the prediction and optimization of the flow boiling heat transfer in mini and micro channel heat sinks

The most common methods for predicting flow boiling heat transfer in mini/micro channels-based heat sinks rely on semi-empirical correlations derived from experimental data. However, these correlations are often limited to specific testing conditions. This study proposes a novel approach using deep learning and genetic algorithms (GA) to predict and optimize refrigerants' flow boiling heat transfer coefficients (FBHTC) in mini/microchannels-based heat sinks. The dataset used in this study includes FBHTC observations from the literature for seven refrigerants (R1234yf, R1234ze, R134A, R513A, R410A, R22, and R32). The optimal input parameters identified include hydraulic diameters ranging from 1 to 7 mm, saturation temperature from 0 to 20 °C, flow qualities from 0.006 to 0.972, heat flux from 3 to 78.8 kW/m2, and mass fluxes between 100 and 1200 kg/m2s. Gradient-boost regression trees were employed to develop the deep learning and GA models for accurate estimation and optimization. Correlation analysis and feature engineering selected the most influential parameters to construct a precise and simple model. The results demonstrate that the models could estimate refrigerants' FBHTC with high accuracy, achieving an R2 of 0.988 and a mean squared error (MSE) of 0.05%. The GA-based method effectively optimized the FBHTC for each refrigerant by determining the appropriate input parameters, including the saturation temperature, heat and mass fluxes, quality, and hydraulic diameter. Additionally, a parametric analysis using explainable artificial intelligence was conducted to interpret the impact of each input parameter on the FBHTC.

Non-uniform maldistribution of boiling flow in two parallel microchannels with entry effects

Abstract Current theoretical studies of parallel channel instability are restricted to the full development assumption, and the entry effect of the existing developmental stage in miniaturized heat dissipation scenarios is not taken into consideration, leading to some deviations in the predicted results of the stability in practical applications. Using a kinetic model of phase transition and pressure drop, we study microchannels with an L/D ratio of 50 and a hydraulic diameter of 200 μm. The significant influence of the entry effect during both the development and fully developed stages is revealed by our results, which demonstrates a narrowing of the overall flow discrepancy and, consequently, an enhancement in the system’s stability. Particularly, the mitigating effect of the entry effect on the flow instability and maldistribution within the parallel channels was gradually enhanced when the qin increased from 20 W/cm2 to 100 W/cm2.

Two-phase flow evolution and interfacial area transport downstream of the mixing-vane spacer grid in rod bundle channels

This study investigated the axial two-phase flow evolution and interfacial area transport characteristics downstream of the mixing vane spacer grid (MVSG) in a tight lattice rod bundle channel with a subchannel hydraulic diameter of Dh = 17.27 mm. The experiment was conducted under the air-water two-phase flow at room temperature and atmospheric pressure. The phase distributions of 6 positions downstream of MVSG (Z/Dh = 9.15, 14.94, 20.73, 26.52, 32.31, and 61.26) were measured utilizing a self-developed double-layer wire-mesh sensor (WMS). 42 cases were obtained, involving the bubbly flow, cap-bubbly flow, and slug flow. Void fraction, bubble velocity, interfacial area concentration (IAC), and bubble size distribution (BSD) databases were built. Under the group-1 (G-1) flow, due to the swirling flow generated by MVSG, the G-1 bubbles were drawn from the bundle surface and the gap region into the subchannel center to form a spindle core-peak distribution. BSD results demonstrate that the swirling flow promotes bubbles’ coalescence. The one-dimensional (1-D) void fraction and IAC downstream of the MVSG decrease first and then recover, which contrasts with the non-mixing vane spacer grid (N-MVSG) effect. It was attributed to the increment of the bubbles’ velocity after MVSG due to the bubbles’ migration to the channel center and coalescence. Under the group-2 (G-2) flow, MVSG intensifies the core peak distribution of void fraction, and the void fraction profile is also spindle-shaped. The obvious bubbles’ break-up occurs at Z/Dh = 14.94 attributed to the swirling decay. The 1-D void fraction and IAC transport indicate, that with the liquid-velocity increment, the MVSG effect is different from the N-MVSG effect, which is mainly achieved by affecting the G-1 bubbles’ dynamic behaviors. The model evaluation utilizing the interfacial area transport equation (IATE) closing bubble interaction source/sink terms indicates that the advection effect dominates the IAC evolution and the contribution of the pressure effect is weak. The remarkable bubbles’ coalescence occurs under the high liquid velocity. In future studies, the additional bubble break-up and coalescence source/sink term model considering the MVSG effect should be developed based on these cases to improve the IATE prediction ability.

Effect of pin fins on heat transfer during condensation in minichannel heat exchanger

In this study, condensation heat transfer of water vapor in a minichannel heat exchanger on smooth and pin fins surfaces was investigated experimentally. The experimental study was carried out to evaluate heat transfer coefficient, overall heat transfer coefficient, and Nusselt number over different ranges of vapor mass flux from 0.0064 to 0.0368 kg m−2 s−1 and cold-water flow rate range between 0.0013 kg s−1 to 0.0057 kg s−1. The minichannel has a rectangular shape with a hydraulic diameter of 1.3 mm. The experimental testing is carried out on aluminum surface with a channel that has a length of 270, width of 30 mm, and height of 1.3 mm. The pin fins surface is on the bottom of the condensing channel and the fins are circular and have diameter, height, and spacing of 1 mm, and are in-inline arrangement. It is found that condensation heat transfer coefficient on pin fins surface is about 15% higher than that on smooth surface. It is found that the condensation heat transfer coefficient increases significantly with increasing the vapor mass flux. In addition, the effect of increasing the cold fluid flow rate is lower than that of the hot vapor flow rate, but it becomes more significant for higher vapor flux.

Variation on hydraulic diameter for the different PETG condenser tube thickness

Variation on hydraulic diameter for the different PETG condenser tube thickness

A fast evaluation method for surface area, volume fraction, and hydraulic diameter of TPMS with different geometric characteristics

A fast evaluation method for surface area, volume fraction, and hydraulic diameter of TPMS with different geometric characteristics

Assessment of Additive Manufactured Micro-Channel Characteristics: Impact of Hydraulic Diameter Evaluation

Abstract Among internal cooling techniques for gas turbines components, skin cooling is recognized as one of the most promising, especially thanks to the spread of additive manufacturing techniques. In this regard, several studies have tried to characterize heat transfer and pressure loss performance of additive manufactured micro-channels, as well as the impact of AM characteristic parameters, mainly in the form of building angle, and resulting surface roughness and channel shape. The open literature offers several correlations in terms of friction and heat transfer coefficient using a representative hydraulic diameter. Despite that, little attention is given on the importance of such characteristic length definition which may lead to under/over estimation of the cooling rates when correlations are used in the design In this work, coupons featuring circular micro-channelswith diameters of 0.5 and 1 mm and different building angles have been tested, to retrieve pressure losses and an average heat transfer coefficient through a lumped approach. The coupons were additive manufactured using the Laser Powder Bed Fusion (L-PBF) technique. Exploiting the experimental survey, the impact of the considered hydraulic diameter was investigated, by comparing different approaches based on a direct measurement. An additional and new approach, aimed at the identification of an “effective” diameter, based on fluid-dynamic considerations, was also considered. The data correlation and the comparison with available correlations from different authors showed the newly proposed approach to provide superior scaling capabilities and, in turn, allowing for the development of more accurate prediction methods

Heat transfer performance of a multiport mini-channel loop thermosyphon for cooling miniaturized electronic devices

Multiport mini-channel thermosyphon with hydraulic diameters less than 2 mm undergoes a severe entrainment problem both in vapor and liquid regions, impeding the flow and affecting the performance characteristics. To address this issue, a novel multiport mini-channel loop thermosyphon (MPMC-LT) design is proposed. This design separates the vapor and liquid flow paths, reducing entrainment and enhancing the performance characteristics. Tests are conducted with heat loads ranging from 5 to 45W at four inclinations (0°, 30°, 60°, and 90°) relative to the horizontal plane. Results revealed that the MPMC-LT design extends the entrainment limit by 10.6 times compared to conventional designs, reducing the interaction between the vapor and condensed working fluid. Additionally, the total entropy and thermal resistance drops by 20.7 % and 19.9 %, respectively at an optimum inclination angle of 30°. It is also noted that the effective thermal conductivity at the optimal inclination is 5.8 times higher compared to the conventional multiport thermosyphon designs. The design effectivity uses the buoyancy, inertial and surface tension forces to reduce resistance in the fluid flow, thereby enhancing the heat transfer characteristics. These results highlight that MPMC-LT is a suitable thermal management device for power electronic systems.

Revised Friction Groups for Evaluating Hydraulic Parameters: Pressure Drop, Flow, and Diameter Estimation

Suitable friction groups are provided for solving three typical hydraulic problems. While the friction group based on viscous forces is used for calculating the pressure drop or head loss in pipes and open channels, commonly referred to as the Type 1 problem in hydraulic engineering, additional friction groups with similar behaviors are introduced for calculating steady flow discharge as the Type 2 problem and, for estimating hydraulic diameter as the Type 3 problem. Contrary to the viscous friction group, the traditional Darcy–Weisbach friction factor demonstrates a negative correlation with the Reynolds number. This results in curves that slope downward from small to large Reynolds numbers on the well-known Moody chart. In contrast, the friction group used here, based on viscous forces, establishes a more appropriate relationship. In this case, the friction and Reynolds number are positively correlated, meaning that both increase or decrease simultaneously. Here, rearranged diagrams for all three mentioned problems show similar behaviors. This paper compares the Moody diagram with the diagram for the viscous force friction group. The turbulent parts of both diagrams are based on the Colebrook equation, with the newly reformulated version using the viscous force friction group. As the Colebrook equation is implicit with respect to friction, requiring an iterative solution, an explicit solution using the Lambert W-function for the reformulated version is offered. Examples are provided for both pipes and open channel flow.

Experimental investigation on pressure drop characteristics of adiabatic two-phase flow in a Gyroid-structured channel

Recent advancements in additive manufacturing techniques have enabled the fabrication of intricate structures. Among these structures, triply periodic minimal surfaces (TPMSs) such as the gyroid, are particularly promising for heat and mass transfer applications owing to their higher surface area to volume ratios compared to conventional structures such as heat exchangers. This study experimentally investigates the pressure drop characteristics of single- and two-phase flows in a gyroid-structured channel under adiabatic conditions. In particular, using an additively manufactured test section with a gyroid-structured channel, the pressure drop characteristics of both single- and two-phase flows are analyzed. The results ofsingle-phase flow experiments reveal that the friction factor depends on the hydraulic diameter, which is defined by the internal volume and surface area of the channel. This suggests that in addition to the hydraulic diameter, other parameters such as porosity and wall thickness must also be considered. Subsequently, the two-phase pressure drop predictions of homogeneous and separated models are compared with the pressure drop data obtained from two-phase flow experiments. The results reveal that gas–liquid separation must be considered to accurately predict the pressure drop in regions influenced by gravitational effects. Furthermore, correlations for predicting the pressure drops both single- and two-phase flows within the operating constraints are proposed.

Microscale heat transfer and phase change mechanisms of carbon dioxide at near-critical conditions

Recent studies highlight the potential benefits of two-phase carbon dioxide (CO2) systems operating near their thermodynamic critical point. Significant heat transfer coefficients and resistance to critical heat flux from improved vapor phase properties can enable heat transfer and reliability enhancements. However, the physical mechanisms of phase change that govern heat transfer behavior in the near-critical regime of reduced pressures (0.85<PR<1, where PR=P/Pcrit) are not well understood. In the present study, subcooled microscale (hydraulic diameter, Dh=500μm) flow boiling experiments are conducted to elucidate the heat transfer behavior near the critical point and the underlying mechanisms that influence it. Experimental data reveals a deviation from the trends predicted by the Cheng correlation (Cheng, L et al., 2008, “New prediction methods for CO2 evaporation inside tubes: Part II—An updated general flow boiling heat transfer model based on flow patterns.” Int. J. Heat Mass Transfer), in particular, displaying a decreasing peak heat transfer coefficient as pressure increases. Heat transfer in the near-critical regime can be further partitioned into two regions: near-critical boiling (0.85<PR<0.96) that exhibits heat transfer highly dependent on surface temperature, and supercritical-like behavior (0.96<PR<1) exhibiting a strong heat transfer dependence on mass flux. High-speed side-view flow visualization reveals the dominant nucleation mechanism (heterogeneous or homogeneous) governs the heat transfer behavior and is adequately predicted using flow inlet conditions. The present study provides a method for predicting nucleation behavior, which, in turn, facilitates the prediction of heat transfer near the critical point. The continuous alteration of heat transfer in the near-critical regime as a function of PR is underscored by this mechanistic insight, bridging the gap between typical flow boiling (PR<0.85) and supercritical heat transfer (PR>1).

Potential of Wake Scaling Techniques for Vertical-Axis Wind Turbine Wake Prediction

Analytical wake models are widely used to predict wind turbine wakes. While these models are well-established for horizontal-axis wind turbines (HAWTs), the analytical wake models for vertical-axis wind turbines (VAWTs) remain under-explored in the wind energy community. In this study, the accuracy of two wake scaling techniques is evaluated to predict the change in the normalized maximum wake velocity deficit behind VAWTs by re-scaling the maximum wake velocity deficit behind an actuator disk with the same thrust coefficient. The wake scaling is defined in terms of equivalent diameter, considering the geometrical properties of the wake-generating object. Two different equivalent diameters are compared, namely the momentum diameter and hydraulic diameter. Different approaches are used to calculate the change in the normalized wake velocity deficit behind a disk with the same thrust coefficient as the VAWT. The streamwise distance is scaled with the equivalent diameter to predict the normalized maximum wake velocity deficit behind the desired VAWT. The performance of the proposed framework is assessed using large-eddy simulation data of VAWTs operating in a turbulent boundary layer with varying operating conditions and aspect ratios. For all of the cases, the momentum diameter scaling provides reasonable predictions of the VAWT normalized maximum wake velocity deficit.

Responses to drought of two Mediterranean ring-porous, deciduous species: Searching for climate smart trees and shrubs

Drought-tolerant tree species with high growth rates and a good capacity for carbon storage in woody tissues (dense wood) are searched for due to aridification. Deciduous, ring-porous tree and shrub species could show such drought tolerance and growth traits, thus representing good candidates for climate-smart rewilding. However, we still do not know the long-term growth rates of these species and how they respond to drought, particularly in climate change hotspots such as the Mediterranean Basin. We analysed these issues at the site and individual levels in two ring-porous, deciduous species (Pistacia terebinthus, Celtis australis) using dendroecology and wood anatomy. The ring width, earlywood vessel diameter, vessel density (VD) and area (%) were measured in two focal sites, one per species, and then growth data were compared with two secondary sites to test if site-to-site synchrony changed through time. Ring-width indices (RWI) and the hydraulic diameter (Dh) of earlywood were calculated. Growth rates (ring width), Dh and vessel area were higher in C. australis (1.03−2.26 mm, 269 μm, 33.9 %) than in P. terebinthus (0.57−0.72 mm, 146 μm, 21.5 %). Consequently, VD was higher in P. terebinthus than in C. australis (104 vs. 61 vessels mm−2). The ring width and Dh were more coupled in P. terebinthus (r = 0.43) than in C. australis (r = 0.32). RWI series of the focal and secondary sites have been synchronized since the 1990s as temperatures rose. Precipitation during the growing season (May, June) enhanced growth and VD of both species. P. terebinthus was more responsive to a drought index than C. australis. The two study species show high growth rates and tolerate drought being thus suitable candidates for climate-smart rewilding.