Showing 1–24 of 24 results for author: Carmeliet, J

Highly inhomogeneous interactions between background climate and urban warming across typical local climate zones in heatwave and non-heatwave days

Authors: Jing Kong, Yongling Zhao, Kai Gao, Dominik Strebel, Jan Carmeliet, Chengwang Lei

Abstract: Urban heat island (UHI) in conjunction with heatwave (HW) leads to exacerbation of thermal stress in urban areas. Prior research on UHI and HW has predominantly concentrated on examining the thermal conditions at the surface and near-surface, with few investigations extending to the radiative and dynamical interactions of UHI and HW, particularly with a focus on the inhomogeneities across local cl… ▽ More Urban heat island (UHI) in conjunction with heatwave (HW) leads to exacerbation of thermal stress in urban areas. Prior research on UHI and HW has predominantly concentrated on examining the thermal conditions at the surface and near-surface, with few investigations extending to the radiative and dynamical interactions of UHI and HW, particularly with a focus on the inhomogeneities across local climate zones (LCZs). Here, we analyse the temperature disparity between HW and non-HW conditions across LCZs in the Sydney area by quantifying the contributions of individual radiative and dynamical processes using the coupled surface-atmosphere climate feedback-response analysis method (CFRAM). Three HW events in 2017, 2019, and 2020 are simulated using the Weather Research and Forecasting (WRF) model coupled with the Single-Layer Urban Canopy Model (SLUCM). The maximum temperature difference between HW and non-HW days may reach up to 10 K, with the increased net solar radiation during HWs being comparable to the typical level of anthropogenic heat flux in urban areas. It is also found that the reduction of clouds, the presence of vapor, and the increase of sensible heat contribute to the warming effect at different levels, with the contribution of clouds being the most dominant. Conversely, the generation of dry convection and the increase of latent heat flux lead to mitigating effects, with the latter being more dominant and capable of causing up to 10 K surface temperature difference between LCZ1 (compact high-rise) and LCZ9 (sparsely built). The differences in the contributions of climate feedback processes across different LCZs become more evident during more severe and humid HWs. These findings underscore the necessity of implementing local climate zone-tailored heat mitigation strategies. △ Less

Submitted 27 May, 2024; originally announced May 2024.

A WRF-UCM-SOLWEIG framework of 10m resolution to quantify the intra-day impact of urban features on thermal comfort

Authors: Xiaotian Ding, Yongling Zhao, Yifan Fan, Jian Ge, Jan Carmeliet

Abstract: City-scale outdoor thermal comfort diagnostics are essential for understanding actual heat stress. However, previous research primarily focused on the street scale. Here, we present the WRF-UCM-SOLWEIG framework to achieve fine-grained thermal comfort mapping at the city scale. The background climate condition affecting thermal comfort is simulated by the Weather Research and Forecasting (WRF) mod… ▽ More City-scale outdoor thermal comfort diagnostics are essential for understanding actual heat stress. However, previous research primarily focused on the street scale. Here, we present the WRF-UCM-SOLWEIG framework to achieve fine-grained thermal comfort mapping at the city scale. The background climate condition affecting thermal comfort is simulated by the Weather Research and Forecasting (WRF) model coupled with the urban canopy model (UCM) at a local-scale (500m). The most dominant factor, mean radiant temperature, is simulated using the Solar and Longwave Environmental Irradiance Geometry (SOLWEIG) model at the micro-scale (10m). The Universal Thermal Climate Index (UTCI) is calculated based on the mean radiant temperature and local climate parameters. The influence of different ground surface materials, buildings, and tree canopies is simulated in the SOLWEIG model using integrated urban morphological data. We applied this proposed framework to the city of Guangzhou, China, and investigated the intra-day variation in the impact of urban morphology during a heat wave period. Through statistical analysis, we found that the elevation in UTCI is primarily attributed to the increase in the fraction of impervious surface (ISF) during daytime, with a maximum correlation coefficient of 0.80. Tree canopy cover has a persistent cooling effect during the day. Implementing 40% of tree cover can reduce the daytime UTCI by 1.5 to 2.0 K. At nighttime, all urban features have a negligible contribution to outdoor thermal comfort. Overall, the established framework provides essential input data and references for studies and urban planners in the practice of urban (micro)climate diagnostics and planning. △ Less

Submitted 27 October, 2023; originally announced October 2023.

An experimental investigation of the heat and flow features in street canyons: Impacts of the approaching turbulent boundary layer flow

Authors: Yunpeng Xue, Yongling Zhao, Shuo-Jun Mei, Yuan Chao, Jan Carmeliet

Abstract: The study of turbulent boundary layer flow holds significant importance in urban climate research, particularly concerning numerical simulation studies where it serves as a crucial inflow boundary condition. However, understanding the turbulent boundary layer's influence on flow and heat features within canyon and canopy flow remains incomplete. To address this knowledge gap, our current work empl… ▽ More The study of turbulent boundary layer flow holds significant importance in urban climate research, particularly concerning numerical simulation studies where it serves as a crucial inflow boundary condition. However, understanding the turbulent boundary layer's influence on flow and heat features within canyon and canopy flow remains incomplete. To address this knowledge gap, our current work employs simultaneous Particle Image Velocimetry and Laser-Induced Fluorescence (PIV-LIF) measurements within a large closed-circuit water tunnel. Through this approach, we obtain valuable flow information under various flow and thermal conditions, allowing us to explore the impacts of three distinct turbulent boundary layer flows. The three chosen turbulent boundary layer flows display distinct influences on flow characteristics and heat removal capacity. The ventilation rate exhibits a maximum difference of 80% among the tested boundary layer flows. Additionally, the most significant variation in heat removal capacity is approximately 45%. Moreover, the different turbulence inlet profiles result in diverse fluctuating features at the canyon opening, while the deeper region of the canyon remains less affected. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Using Clustering to Understand Intra-city Warming in Heatwaves: Insights into Paris, Montreal, and Zurich

Authors: Yongling Zhao, Dominik Strebel, Dominique Derome, Igor Esau, Qi Li, Jan Carmeliet

Abstract: We introduce a novel methodological advancement by clustering paired near-surface air temperature with the planetary boundary layer height (PBLH) to characterize intra-city clusters for analytics. To illustrate this approach, we analyze three heatwaves (HW): the 2019 HW in Paris, the 2018 HW in Montreal, and the 2017 HW in Zurich. We assess cluster-based characteristics before, during, and after h… ▽ More We introduce a novel methodological advancement by clustering paired near-surface air temperature with the planetary boundary layer height (PBLH) to characterize intra-city clusters for analytics. To illustrate this approach, we analyze three heatwaves (HW): the 2019 HW in Paris, the 2018 HW in Montreal, and the 2017 HW in Zurich. We assess cluster-based characteristics before, during, and after heatwave events. Using the objective hysteresis model, we determine the overall strength coefficient of the hysteresis loop between ground storage flux and all-wave downward radiative flux, ranging from 0.414 to 0.457 for urban clusters and from 0.126 to 0.157 for rural clusters during the heatwave periods. Across all cities, we observe a consistent refueling-restoration mode in the cumulative ground heat flux as the heatwaves progress. Future developments of this proposed two-component clustering approach, with the integration of more influential physics, will offer a more comprehensive characterization of cities for urban climate analytics. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Critical transitions on route to chaos of natural convection on a heated horizontal circular surface

Authors: Yuhan Jiang, Yongling Zhao, Jan Carmeliet, Bingchuan Nie, Feng Xu

Abstract: The transition route and bifurcations of the buoyant flow developing on a heated circular horizontal surface are elaborated using direct numerical simulations and direct stability analysis. A series of bifurcations, as a function of Rayleigh numbers (Ra) ranging from $10^1$ to $6\times10^7$, are found on the route to the chaos of the flow at $Pr=7$. When $Ra<1.0\times10^3$, the buoyant flow above… ▽ More The transition route and bifurcations of the buoyant flow developing on a heated circular horizontal surface are elaborated using direct numerical simulations and direct stability analysis. A series of bifurcations, as a function of Rayleigh numbers (Ra) ranging from $10^1$ to $6\times10^7$, are found on the route to the chaos of the flow at $Pr=7$. When $Ra<1.0\times10^3$, the buoyant flow above the heated horizontal surface is dominated by conduction, because of which distinct thermal boundary layer and plume are not present. At $Ra=1.1\times10^6$, a Hopf bifurcation occurs, resulting in the flow transition from a steady state to a periodic puffing state. As Ra increases further, the flow enters a periodic rotating state at $Ra=1.9\times10^6$, which is a unique state that was rarely discussed in the literature. These critical transitions, leaving from a steady state and subsequently entering a series of periodic states (puffing, rotating, flapping and doubling) and finally leading to chaos, are diagnosed using spectral analysis and two-dimensional Fourier Transform (2DFT). Moreover, direct stability analysis is conducted by introducing random numerical perturbations into the boundary condition of the surface heating. We find that when the state of a flow is in the vicinity of bifurcation points (e.g., $Ra=2.0\times10^6$), the flow is conditionally unstable to perturbations, and it can bifurcate from the rotating state to the flapping state in advance. However, for relatively stable flow states, such as at $Ra=1.5\times10^6$, the flow remains its periodic puffing state even though it is being perturbed. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Impact of street canyon morphology on heat and fluid flow-an experimental water tunnel study using simultaneous PIV-LIF technique

Authors: Yunpeng Xue, Yongling Zhao, Shuo-Jun Mei, Yuan Chao, Jan Carmeliet

Abstract: Urban areas are known for their complex atmospheric environments, with the building morphology having a significant impact on local climate patterns, air quality, and overall urban microclimate. Understanding the heat transport and fluid flow in complex urban environments is crucial for improving urban climate resilience, which remains an open frontier in the field of urban studies. To gain a more… ▽ More Urban areas are known for their complex atmospheric environments, with the building morphology having a significant impact on local climate patterns, air quality, and overall urban microclimate. Understanding the heat transport and fluid flow in complex urban environments is crucial for improving urban climate resilience, which remains an open frontier in the field of urban studies. To gain a more profound insight into the physical processes occurring in urban areas, particularly within street canyons, we conducted an experimental investigation in a large-scale water tunnel. This study involved the simultaneous examination of heat and flow fields, carried out at high spatial and temporal resolutions, utilizing Laser-induced Fluorescence (LIF) for heat analysis and Particle Image Velocimetry (PIV) for flow analysis. Our results of heat and flow in different street canyons indicate that the flow is significantly influenced by a combination of factors, including canyon configuration, the presence of buoyant force, and the magnitude of the approaching flow. The ventilation rate and heat flux from the street canyon, which are key factors shaping the urban microclimate, are found dominated significantly by the street canyon morphology. For instance, changing the aspect ratio of a street canyon results in a significant change of air ventilation rate, ranging from as low as 0.02 to as high as 1.5 under the same flow conditions. Additionally, canyons with high air ventilation rates exhibit significant heat flux removal at the canyon roof level, which is accurately described by the local Richardson number. △ Less

Submitted 19 October, 2023; v1 submitted 8 June, 2023; originally announced June 2023.

Journal ref: Experimental Thermal and Fluid Science. 2024, 150, 111066

Understanding the Impact of Heatwave on Urban Heat Island in Greater Sydney: Temporal Surface Energy Budget Change with Land Types

Authors: Jing Kong, Yongling Zhao, Dominik Strebel, Kai Gao, Jan Carmeliet, Chengwang Lei

Abstract: The impact of heatwaves (HWs) on urban heat island (UHI) is a contentious topic with contradictory research findings. A comprehensive understanding of the response of urban and rural areas to HWs, considering the underlying cause of surface energy budget changes, remains elusive. This study attempts to address this gap by investigating a 2020 HW event in the Greater Sydney Area using the Advanced… ▽ More The impact of heatwaves (HWs) on urban heat island (UHI) is a contentious topic with contradictory research findings. A comprehensive understanding of the response of urban and rural areas to HWs, considering the underlying cause of surface energy budget changes, remains elusive. This study attempts to address this gap by investigating a 2020 HW event in the Greater Sydney Area using the Advanced Weather Research and Forecasting (WRF) model. Findings indicate that the HW intensifies the nighttime surface UHI by approximately 4°C. An analysis of surface energy budgets reveals that urban areas store more heat during the HW due to receiving more solar radiation and less evapotranspiration compared to rural areas. The maximum heat storage flux in urban during the HW can be around 200 W/m2 higher than that during post-HW. The stored heat is released at nightime, raising the air temperature in the urban areas. Forests and savannas have relatively lower storage heat fluxes due to high transpiration and albedo, and the maximum heat storage flux is only around 50 W/m2 higher than that during post-HW. In contrast, a negative synergistic effect is detected between the 2-m UHI and HW. This may be because other meteorological conditions including wind have substantial impacts on the air temperature pattern. The strong hot and dry winds coming from the west and the proximity of tall buildings to the coast diminish the sea breeze coming from the east, resulting in a higher air temperature in the western urban district. Meanwhile, the western forest area also experiences higher temperatures due to the westward winds. In addition, changes in wind direction alter the temperature distribution in the northern rural region. Based on the present study, urban climate simulation data and associated findings can be used to develop urban heat mitigation strategies for UHI during HW. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Aligning Three-Decade Surge in Urban Cooling with Global Warming

Authors: Haiwei Li, Yongling Zhao, Ronita Bardhan, Pak-Wai Chan, Dominique Derome, Zhiwen Luo, Diana Urge-Vorsatz, Jan Carmeliet

Abstract: Rising demand for space cooling has been placing enormous strain on various technological, environmental, and societal dimensions, resulting in issues related to energy consumption, environmental sustainability, health and well-being, affordability, and equity. Holistic approaches that combine energy efficiency optimization, policy-making, and societal adaptation must be rapidly promoted as viable… ▽ More Rising demand for space cooling has been placing enormous strain on various technological, environmental, and societal dimensions, resulting in issues related to energy consumption, environmental sustainability, health and well-being, affordability, and equity. Holistic approaches that combine energy efficiency optimization, policy-making, and societal adaptation must be rapidly promoted as viable, timely solutions. We interpret the 30-year upward trend and spikes in urban cooling demand from the perspective of climate change, urbanization, and background climates, focusing on five representative cities: Hong Kong, Sydney, Montreal, Zurich, and London. An unequivocal, worrying upward trend in cooling demand is observed in meteorological data from 1990 to 2021, using cooling degree hours (CDH) as a city-scale metric. The surge in cooling energy demand can be largely attributed to global warming, urban heat islands, and extreme heat events. Further, our quantification of the impact of the base temperature, in relation to the historical CDH, reveals that a 20% energy saving could be achieved instantly within a rather broad range of temperature and humidity by increasing the setpoint temperature by one degree, while characteristic sensational and physiological levels can be maintained at 'acceptable' and 'physiological thermal neutrality' respectively. However, the potential of reducing cooling demand can be nonlinearly and significantly lowered due to the presence of compound high relative humidity and high air temperature. To reduce cooling energy demand rapidly in a warming climate, we highlight the necessity of promoting hard and soft behavioral adaptation along with regulatory intervention for the operation of space cooling systems. △ Less

Submitted 16 June, 2023; v1 submitted 27 January, 2023; originally announced January 2023.

Comments: 9 pages, 2 figures

Universal scaling of nanoparticle deposition by colloidal droplet drying

Authors: Feifei Qin, Linlin Fei, Jianlin Zhao, Qinjun Kang, Sauro Succi, Dominique Derome, Jan Carmeliet

Abstract: We present a comprehensive study of nanoparticle deposition from drying of colloidal droplets. By means of lattice Boltzmann modeling and theoretical analysis, various deposition patterns, including mountain-like, uniform and coffee ring, as well as un-/symmetrical multiring/mountain-like patterns are achieved. The ratio of nanoparticles deposited at droplet peripheries and center is proposed to q… ▽ More We present a comprehensive study of nanoparticle deposition from drying of colloidal droplets. By means of lattice Boltzmann modeling and theoretical analysis, various deposition patterns, including mountain-like, uniform and coffee ring, as well as un-/symmetrical multiring/mountain-like patterns are achieved. The ratio of nanoparticles deposited at droplet peripheries and center is proposed to quantify different patterns. Its value is controlled by the competition between the capillary flow and nanoparticle diffusion, leading to a linear dependence on an effective Péclet number, across over three orders of magnitude. Remarkably, the final deposition pattern can be predicted based on material properties only. △ Less

Submitted 25 January, 2023; originally announced January 2023.

Time-evolving Impact of Trees on Street Canyon Microclimate

Authors: Haiwei Li, Yongling Zhao, Ronita Bardhan, Aytac Kubilay, Dominique Derome, Jan Carmeliet

Abstract: Nowadays, cities are frequently exposed to heatwaves, worsening the outdoor thermal comfort and increasing cooling energy demand in summer. Urban forestry is seen as one of the viable and preferable solutions to combating extreme heat events and urban heat island (UHI) in times of climate change. While many cities have initiated tree-planting programmes in recent years, the evolving impact of tree… ▽ More Nowadays, cities are frequently exposed to heatwaves, worsening the outdoor thermal comfort and increasing cooling energy demand in summer. Urban forestry is seen as one of the viable and preferable solutions to combating extreme heat events and urban heat island (UHI) in times of climate change. While many cities have initiated tree-planting programmes in recent years, the evolving impact of trees on street microclimate, in a time span of up to several decades, remains unclear. We investigate the cooling effects of linden trees in five groups, i.e., 10-20, 20-30, 30-40, 40-60, and 60-100 years old. The leaf area index (LAI) and leaf area density (LAD) vary nonlinearly as the trees grow, peaking at different ages. Computational fluid dynamics (CFD) simulations solving microclimate are performed for an idealized street canyon with trees of varied age groups. Turbulent airflow, heat and moisture transport, shortwave and longwave radiation, shading and transpiration are fully coupled and solved in OpenFOAM. The meteorological data, including air temperature, wind speed, moisture, and shortwave radiation of the heatwave in Zurich (June 2019), are applied as boundary conditions. The results show that young trees in the age group of 10-20 years old provide little heat mitigation at the pedestrian level in an extreme heat event. Optimal heat mitigation by trees is observed for the group of 30-60 years old trees. Finally, the potential impact of growing trees as a heat mitigation measure on air ventilation is evaluated. △ Less

Submitted 24 January, 2023; originally announced January 2023.

Comments: 4 figures, 8 pages

The time-evolving impact of tree size on nighttime street canyon microclimate: Wind tunnel modeling of aerodynamic effects and heat removal

Authors: Yongling Zhao, Haiwei Li, Ronita Bardhan, Aytac Kubilay, Qi Li, Jan Carmeliet

Abstract: Urban trees play a crucial role in urban climate in many aspects. However, existing research has not adequately explored the impact from a time-evolving perspective, that is, tree growth over time. To bridge this research gap, this study investigates in a wind tunnel the effects of tree-to-canyon foliage cover and relative height (0.32-1.1 times canyon height), mimicking growth of trees, on condit… ▽ More Urban trees play a crucial role in urban climate in many aspects. However, existing research has not adequately explored the impact from a time-evolving perspective, that is, tree growth over time. To bridge this research gap, this study investigates in a wind tunnel the effects of tree-to-canyon foliage cover and relative height (0.32-1.1 times canyon height), mimicking growth of trees, on conditions in street canyons during moderate and extreme heat. The results reveal that trees may affect canyon-wide ventilation and heat removal in two different scenarios. First, when canyons are in isothermal conditions, medium and large trees, that fill half the canyon height or reach slightly above the canyon, decelerate the shear layer and weaken the vortical flow, as a result reducing the canyon-wide ventilation. Second, in extreme heat conditions, medium and large trees trap heat at the pedestrian level due to the blockage of air entrainment and the suppression of upward buoyancy-driven flow from the ground surface. An air temperature rise that corresponds to 1.5 degree Celsius in a full-scale urban setting is observed in measurements. These observations suggest that urban trees' foliage cover must be managed for a canyon's optimal ventilation and heat removal during nighttime. △ Less

Submitted 13 January, 2023; originally announced January 2023.

Comments: 29 pages, 10 figures

Fluid Tunnel Research for Challenges of Urban Climate

Authors: Yongling Zhao, Lup Wai Chew, Yifan Fan, Christof Gromke, Jian Hang, Yichen Yu, Alessio Ricci, Yan Zhang, Yunpeng Xue, Sofia Fellini, Parham A. Mirzaei, Naiping Gao, Matteo Carpentieri, Pietro Salizzoni, Jianlei Niu, Jan Carmeliet

Abstract: Experimental investigations using wind and water tunnels have long been a staple of fluid mechanics research for a large number of applications. These experiments often single out a specific physical process to be investigated, while studies involving multiscale and multi-physics processes are rare due to the difficulty and complexity in the experimental setup. In the era of climate change, there… ▽ More Experimental investigations using wind and water tunnels have long been a staple of fluid mechanics research for a large number of applications. These experiments often single out a specific physical process to be investigated, while studies involving multiscale and multi-physics processes are rare due to the difficulty and complexity in the experimental setup. In the era of climate change, there is an increasing interest in innovative experimental studies in which fluid (wind and water) tunnels are employed for modelling multiscale, multi-physics phenomena of the urban climate. High-quality fluid tunnel measurements of urban-physics related phenomena are also much needed to facilitate the development and validation of advanced multi-physics numerical models. As a repository of knowledge in modelling these urban processes, we cover fundamentals, recommendations and guidelines for experimental design, recent advances and outlook on eight selected research areas, including (i) thermal buoyancy effects of urban airflows, (ii) aerodynamic and thermal effects of vegetation, (iii) radiative and convective heat fluxes over urban materials, (iv) influence of thermal stratification on land-atmosphere interactions, (v) pollutant dispersion, (vi) indoor and outdoor natural ventilation, (vii) wind thermal comfort, and (viii) urban winds over complex urban sites. Further, three main challenges, i.e., modelling of multi-physics, modelling of anthropogenic processes, and combined use of fluid tunnels, scaled outdoor and field measurements for urban climate studies, are discussed. △ Less

Submitted 8 January, 2023; originally announced January 2023.

Comments: 36 pages,9 figures, 1 table

Betweenness centrality illuminates intermittent frictional dynamics

Authors: Omid Dorostkar, Karen E. Daniels, Dominik Strebel, Jan Carmeliet

Abstract: Dense granular systems subjected to an imposed shear stress undergo stick-slip dynamics with systematic patterns of dilation-compaction. During each stick phase, as the frictional strength builds up, the granular system dilates to accommodate shear strain, developing stronger force networks. During each slip event, when the stored energy is released, particles experience large rearrangements and t… ▽ More Dense granular systems subjected to an imposed shear stress undergo stick-slip dynamics with systematic patterns of dilation-compaction. During each stick phase, as the frictional strength builds up, the granular system dilates to accommodate shear strain, developing stronger force networks. During each slip event, when the stored energy is released, particles experience large rearrangements and the granular network can significantly change. Here, we use numerical simulations of 3D, sheared frictional packings to show that the mean betweenness centrality -- a property of network of interparticle connections -- follows consistent patterns during the stick-slip dynamics, showing sharp spikes at each slip event. We identify the source of this behavior as arising from the connectivity and contact arrangements of granular network during dilation-compaction cycles, and find that a lower potential for connection between particles leads to an increase of mean betweenness centrality in the system. Furthermore, we show that at high confinements, few particles lose contact during slip events, leading to a smaller change in granular connectivity and betweenness centrality. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Droplet impact of Newtonian fluids and blood on simple fabrics: effect of fabric pore size and underlying substrate

Authors: T. C. de Goede, A. M. Moqaddam, K. C. M. Limpens, S. A. Kooij, D. Derome, J. Carmeliet, N. Shahidzadeh, D. Bonn

Abstract: When a droplet impacts a fabric mesh at a sufficiently high impact velocity, it not only spreads over the fabric but also penetrate its pores. To determine the influence of this liquid penetration of the fabric on droplet spreading on thin fabric meshes, we measured the droplet spreading ratio on fabric with and without an underlying substrate using a high-speed camera. For fabrics without a subst… ▽ More When a droplet impacts a fabric mesh at a sufficiently high impact velocity, it not only spreads over the fabric but also penetrate its pores. To determine the influence of this liquid penetration of the fabric on droplet spreading on thin fabric meshes, we measured the droplet spreading ratio on fabric with and without an underlying substrate using a high-speed camera. For fabrics without a substrate, the droplet spreading ratio is reduced as the fabric penetration by the liquid reduces the droplet volume spreading on top of the fabric. Using entropic lattice Boltzmann simulations, we find that the lower droplet spreading ratio on fabrics, both with and without a substrate, is due to an increase of viscous losses inside the droplet during spreading. Comparing droplet impact of blood with its Newtonian counterpart, we show that for spreading on fabrics, just like on smooth surfaces, blood can be approximated as a Newtonian fluid. △ Less

Submitted 2 January, 2020; originally announced January 2020.

Comments: 14 pages, 9 figures, 2 Supplemental Movies (not included)

Machine Learning Reveals the State of Intermittent Frictional Dynamics in a Sheared Granular Fault

Authors: C. X. Ren, O. Dorostkar, B. Rouet-Leduc, C. Hulbert, D. Strebel, R. A. Guyer, P. A. Johnson, J. Carmeliet

Abstract: Seismogenic plate boundaries are presumed to behave in a similar manner to a densely packed granular medium, where fault and blocks systems rapidly rearrange the distribution of forces within themselves, as particles do in slowly sheared granular systems. We use machine learning and show that statistical features of velocity signals from individual particles in a simulated sheared granular fault c… ▽ More Seismogenic plate boundaries are presumed to behave in a similar manner to a densely packed granular medium, where fault and blocks systems rapidly rearrange the distribution of forces within themselves, as particles do in slowly sheared granular systems. We use machine learning and show that statistical features of velocity signals from individual particles in a simulated sheared granular fault contain information regarding the instantaneous global state of intermittent frictional stick-slip dynamics. We demonstrate that combining features built from the signals of more particles can improve the accuracy of the global model, and discuss the physical basis behind decrease in error. We show that the statistical features such as median and higher moments of the signals that represent the particle displacement in the direction of shearing are among the best predictive features. Our work provides novel insights into the applications of machine learning in studying frictional processes that take place in geophysical systems. △ Less

Submitted 26 March, 2019; originally announced March 2019.

Unveiling dynamic changes in the diurnal microclimate of a Buxus sempervirens with non-intrusive imaging of flow field, leaf temperature, and plant microstructure

Authors: Lento Manickathan, Thijs Defraeye, Stephan Carl, Henning Richter, Jonas Allegrini, Dominique Derome, Jan Carmeliet

Abstract: Plant response is not only dependent on the atmospheric evaporative demand due to the combined effects of wind speed, air temperature, humidity, and solar radiation, but is also dependent on the water transport within the leaf-xylem-root system. Therefore, a detailed understanding of such dynamics is key to the development of appropriate mitigation strategies and numerical models. In this study, w… ▽ More Plant response is not only dependent on the atmospheric evaporative demand due to the combined effects of wind speed, air temperature, humidity, and solar radiation, but is also dependent on the water transport within the leaf-xylem-root system. Therefore, a detailed understanding of such dynamics is key to the development of appropriate mitigation strategies and numerical models. In this study, we unveil the diurnal dynamics of the microclimate of a Buxus sempervirens plant using multiple high-resolution non-intrusive imaging techniques. The wake flow field is measured using stereoscopic particle image velocimetry, the spatiotemporal leaf temperature history is obtained using infrared thermography, and additionally, the plant porosity is obtained using X-ray tomography. We find that the wake velocity statistics is not directly linked with the distribution of the porosity but depends mainly on the geometry of the plant foliage which generates the shear flow. The interaction between the shear regions and the upstream boundary layer profile is seen to have a dominant effect on the wake turbulent kinetic energy distribution. Furthermore, the leaf area density distribution has a direct impact on the short-wave radiative heat flux absorption inside the foliage where 50% of the radiation is absorbed in the top 20% of the foliage. This localized radiation absorption results in a high local leaf and air temperature. Furthermore, a comparison of the diurnal variation of leaf temperature and the net plant transpiration rate enabled us to quantify the diurnal hysteresis resulting from the stomatal response lag. The day of this plant is seen to comprise of four stages of climatic conditions: no-cooling, high-cooling, equilibrium, and decaying-cooling stages. △ Less

Submitted 6 March, 2019; originally announced March 2019.

Comments: 44 pages, 19 figures, preprint

Sprays from droplets impacting a mesh

Authors: Stefan Kooij, Ali Mazloomi Moqaddam, Thijs C. de Goede, Dominique Derome, Jan Carmeliet, Noushine Shahidzadeh, Daniel Bonn

Abstract: In liquid spray applications, the sprays are often created by the formation and destabilization of a liquid sheet or jet. The disadvantage of such atomization processes is that the breakup is often highly irregular, causing a broad distribution of droplet sizes. As these sizes are controlled by the ligament corrugation and size, a monodisperse spray should consist of ligaments that are both smooth… ▽ More In liquid spray applications, the sprays are often created by the formation and destabilization of a liquid sheet or jet. The disadvantage of such atomization processes is that the breakup is often highly irregular, causing a broad distribution of droplet sizes. As these sizes are controlled by the ligament corrugation and size, a monodisperse spray should consist of ligaments that are both smooth and of equal size. A straightforward way of creating smooth and equally sized ligaments is by droplet impact on a mesh. In this work we show that this approach does however not produce monodisperse droplets, but instead the droplet size distribution is very broad, with a large number of small satellite drops. We demonstrate that the fragmentation is controlled by a jet instability, where initial perturbations caused by the injection process result in long-wavelength disturbances that determine the final ligament breakup. During destabilization the crests of these disturbances are connected by thin ligaments which are the leading cause of the large number of small droplets. A secondary coalescence process, due to small relative velocities between droplets, partly masks this effect by reducing the amount of small droplets. Of the many parameters in this system, we describe the effect of varying the mesh size, mesh rigidity, impact velocity, wetting properties, keeping the liquid properties the same by focusing on water droplets only. We further perform Lattice Boltzmann modeling of the impact process that reproduces key features seen in the experimental data. △ Less

Submitted 6 December, 2018; originally announced December 2018.

Comments: 11 pages, 12 figures, 6 supplemental movies (not included)

Modeling of Stick-Slip Behavior in Sheared Granular Fault Gouge Using the Combined Finite-Discrete Element Method

Authors: Ke Gao, Bryan J. Euser, Esteban Rougier, Robert A. Guyer, Zhou Lei, Earl E. Knight, Jan Carmeliet, Paul A. Johnson

Abstract: Sheared granular layers undergoing stick slip behavior are broadly employed to study the physics and dynamics of earthquakes. Here, a two dimensional implementation of the combined finite discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular fault system including both gouge and plate, a… ▽ More Sheared granular layers undergoing stick slip behavior are broadly employed to study the physics and dynamics of earthquakes. Here, a two dimensional implementation of the combined finite discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular fault system including both gouge and plate, and to investigate the influence of different normal loads on seismic moment, macroscopic friction coefficient, kinetic energy, gouge layer thickness, and recurrence time between slips. In the FDEM model, the deformation of plates and particles is simulated using the FEM formulation while particle particle and particle plate interactions are modeled using DEM derived techniques. The simulated seismic moment distributions are generally consistent with those obtained from the laboratory experiments. In addition, the simulation results demonstrate that with increasing normal load, (i) the kinetic energy of the granular fault system increases; (ii) the gouge layer thickness shows a decreasing trend; and (iii) the macroscopic friction coefficient does not experience much change. Analyses of the slip events reveal that, as the normal load increases, more slip events with large kinetic energy release and longer recurrence time occur, and the magnitude of gouge layer thickness decrease also tends to be larger; while the macroscopic friction coefficient drop decreases. The simulations not only reveal the influence of normal loads on the dynamics of sheared granular fault gouge, but also demonstrate the capabilities of FDEM for studying stick slip dynamic behavior of granular fault systems. △ Less

Submitted 15 May, 2018; v1 submitted 4 April, 2018; originally announced April 2018.

Dynamic induced softening in frictional granular material investigated by DEM simulation

Authors: Laure Lemrich, Paul A. Johnson, Robert Guyer, Xiaoping Jia, Jan Carmeliet

Abstract: A granular system composed of frictional glass beads is simulated using the Discrete Element Method. The inter-grain forces are based on the Hertz contact law in the normal direction with frictional tangential force. The damping due to collision is also accounted for. Systems are loaded at various stresses and their quasi-static elastic moduli are characterized. Each system is subjected to an exte… ▽ More A granular system composed of frictional glass beads is simulated using the Discrete Element Method. The inter-grain forces are based on the Hertz contact law in the normal direction with frictional tangential force. The damping due to collision is also accounted for. Systems are loaded at various stresses and their quasi-static elastic moduli are characterized. Each system is subjected to an extensive dynamic testing protocol by measuring the resonant response to a broad range of AC drive amplitudes and frequencies via a set of diagnostic strains. The system, linear at small AC drive amplitudes has resonance frequencies that shift downward (i.e., modulus softening) with increased AC drive amplitude. Detailed testing shows that the slipping contact ratio does not contribute significantly to this dynamic modulus softening, but the coordination number is strongly correlated to this reduction. This suggests that the softening arises from the extended structural change via break and remake of contacts during the rearrangement of bead positions driven by the AC amplitude. △ Less

Submitted 19 June, 2017; originally announced June 2017.

Comments: acoustics, nonlinearity, granular media

MSC Class: 74E20 (Primary) 70K30 (Secondary) ACM Class: I.6.6

Journal ref: Phys. Rev. E 96, 062901 (2017)

The role of water in the behavior of wood

Authors: Dominique Derome, Ahmad Rafsanjani, Stefan Hering, Martin Dressler, Alessandra Patera, Christian Lanvermann, Marjan Sedighi-Gilani, Falk K. Wittel, Peter Niemz, Jan Carmeliet

Abstract: Wood, due to its biological origin, has the capacity to interact with water. Sorption/desorption of moisture is accompanied with swelling/shrinkage and softening/hardening of its stiffness. The correct prediction of the behavior of wood components undergoing environmental loading requires that the moisture behavior and mechanical behavior of wood are considered in a coupled manner. We propose a co… ▽ More Wood, due to its biological origin, has the capacity to interact with water. Sorption/desorption of moisture is accompanied with swelling/shrinkage and softening/hardening of its stiffness. The correct prediction of the behavior of wood components undergoing environmental loading requires that the moisture behavior and mechanical behavior of wood are considered in a coupled manner. We propose a comprehensive framework using a fully coupled poromechanical approach, where its multiscale implementation provides the capacity to take into account, directly, the exact geometry of the wood cellular structure, using computational homogenization. A hierarchical model is used to take into account the subcellular composite-like organization of the material. Such advanced modeling requires high resolution experimental data for the appropriate determination of inputs and for its validation. High-resolution x-ray tomography, digital image correlation, and neutron imaging are presented as valuable methods to provide the required information. △ Less

Submitted 17 September, 2015; originally announced September 2015.

Journal ref: Journal of Building Physics 36 (4/2013) 398-421

Computational up-scaling of anisotropic swelling and mechanical behavior of hierarchical cellular material

Authors: Ahmad Rafsanjani, Dominique Derome, Falk K. Wittel, Jan Carmeliet

Abstract: The hygro-mechanical behavior of a hierarchical cellular material, i.e. growth rings of softwood is investigated using a two-scale micro-mechanics model based on a computational homogenization technique. The lower scale considers the individual wood cells of varying geometry and dimensions. Honeycomb unit cells with periodic boundary conditions are utilized to calculate the mechanical properties a… ▽ More The hygro-mechanical behavior of a hierarchical cellular material, i.e. growth rings of softwood is investigated using a two-scale micro-mechanics model based on a computational homogenization technique. The lower scale considers the individual wood cells of varying geometry and dimensions. Honeycomb unit cells with periodic boundary conditions are utilized to calculate the mechanical properties and swelling coefficients of wood cells. Using the cellular scale results, the anisotropy in mechanical and swelling behavior of a growth ring in transverse directions is investigated. Predicted results are found to be comparable to experimental data. It is found that the orthotropic swelling properties of the cell wall in thin-walled earlywood cells produce anisotropic swelling behavior while, in thick latewood cells, this anisotropy vanishes. The proposed approach provides the ability to consider the complex microstructure when predicting the effective mechanical and swelling properties of softwood. △ Less

Submitted 4 September, 2015; originally announced September 2015.

Journal ref: Composites Science and Technology 72 (6/2012) 744-751

Acoustically-induced slip in sheared granular layers: application to dynamic earthquake triggering

Authors: Behrooz Ferdowsi, Michele Griffa, Robert A. Guyer, Paul A. Johnson, Chris Marone, Jan Carmeliet

Abstract: A fundamental mystery in earthquake physics is ``how can an earthquake be triggered by distant seismic sources?'' Here, we use discrete element method simulations of a granular layer, during stick-slip, that is subject to transient vibrational excitation to gain further insight into the physics of dynamic earthquake triggering. Using Coulomb friction law for grains interaction, we observe delayed… ▽ More A fundamental mystery in earthquake physics is ``how can an earthquake be triggered by distant seismic sources?'' Here, we use discrete element method simulations of a granular layer, during stick-slip, that is subject to transient vibrational excitation to gain further insight into the physics of dynamic earthquake triggering. Using Coulomb friction law for grains interaction, we observe delayed triggering of slip in the granular gouge. We find that at a critical vibrational amplitude (strain) there is an abrupt transition from negligible time-advanced slip (clock advance) to full clock advance, {\it i.e.}, transient vibration and triggered slip are simultaneous. The critical strain is order of $10^{-6}$, similar to observations in the laboratory and in Earth. The transition is related to frictional weakening of the granular layer due to a dramatic decrease in coordination number and the weakening of the contact force network. Associated with this frictional weakening is a pronounced decrease in the elastic modulus of the layer. The study has important implications for mechanisms of triggered earthquakes and induced seismic events and points out the underlying processes in response of the fault gouge to dynamic transient stresses. △ Less

Submitted 3 November, 2015; v1 submitted 16 July, 2015; originally announced July 2015.

Metastability limit for the nucleation of NaCl crystals in confinement

Authors: Julie Desarnaud, Hannelore Derluyn, Jan Carmeliet, Daniel Bonn, Noushine Shahidzadeh

Abstract: We study the spontaneous nucleation and growth of sodium chloride crystals induced by controlled evaporation in confined geometries (microcapillaries) spanning several orders of magnitude in volume. In all experiments, the nucleation happens reproducibly at a very high supersaturation S~1.6 and is independent of the size, shape and surface properties of the microcapillary. We show from classical n… ▽ More We study the spontaneous nucleation and growth of sodium chloride crystals induced by controlled evaporation in confined geometries (microcapillaries) spanning several orders of magnitude in volume. In all experiments, the nucleation happens reproducibly at a very high supersaturation S~1.6 and is independent of the size, shape and surface properties of the microcapillary. We show from classical nucleation theory that this is expected: S~1.6 corresponds to the point where nucleation first becomes observable on experimental time scales. A consequence of the high supersaturations reached at the onset of nucleation is the very rapid growth of a single skeletal (Hopper) crystal. Experiments on porous media reveal also the formation of Hopper crystals in the entrapped liquid pockets in the porous network and consequently underline the fact that sodium chloride can easily reach high supersaturations, in spite of what is commonly assumed for this salt. △ Less

Submitted 30 January, 2014; originally announced January 2014.

Comments: 16 pages, 6 figures

Effect of boundary vibration on the frictional behavior of a dense sheared granular layer

Authors: B. Ferdowsi, M. Griffa, R. A. Guyer, P. A. Johnson, J. Carmeliet

Abstract: We report results of 3D Discrete Element Method (DEM) simulations aiming at investigating the role of the boundary vibration in inducing frictional weakening in sheared granular layers. We study the role of different vibration amplitudes applied at various shear stress levels, for a granular layer in the stick-slip regime and in the steady-sliding regime. Results are reported in terms of friction… ▽ More We report results of 3D Discrete Element Method (DEM) simulations aiming at investigating the role of the boundary vibration in inducing frictional weakening in sheared granular layers. We study the role of different vibration amplitudes applied at various shear stress levels, for a granular layer in the stick-slip regime and in the steady-sliding regime. Results are reported in terms of friction drops and kinetic energy release associated with frictional weakening events. We find that larger vibration amplitude induces larger frictional weakening events. The results show evidence of a threshold below which no induced frictional weakening takes place. Friction drop size is found to be dependent on the shear stress at the time of vibration. A significant increase in the ratio between the number of slipping contacts to the number of sticking contacts in the granular layer is observed for large vibration amplitudes. These vibration-induced contact rearrangements enhance particle mobilization and induces a friction drop and kinetic energy release. This observation provides some insight into the grain-scale mechanisms of frictional weakening by boundary vibration in a dense sheared granular layer. In addition to characterizing the basic physics of vibration induced shear weakening, we are attempting to understand how a fault fails in the earth under seismic wave forcing. This is the well know phenomenon of dynamic earthquake triggering. We believe that the granular physics are key to this understanding. △ Less

Submitted 24 January, 2014; originally announced January 2014.