Accurately predicting bare‐soil evaporation requires the proper characterization of the near‐surface atmospheric conditions. These conditions, dependent on factors such as surface microtopography and wind velocity, vary greatly and therefore require high‐resolution datasets to be fully incorporated into evaporation models. These factors are oftentimes parameterized in models through the aerodynamic resistance (
Evaporation from undulating soil surfaces is rarely studied due to limited modeling theory and inadequate experimental data linking dynamic soil and atmospheric interactions. The goal of this paper is to provide exploratory insights into evaporation behavior from undulating soil surfaces under turbulent conditions through numerical and experimental approaches. A previously developed and verified coupled free flow and porous media flow model was extended by incorporating turbulent airflow through Reynolds‐averaged Navier–Stokes equations. The model explicitly describes the relevant physical processes and the key properties in the free flow, porous media, and at the interface, allowing for the analysis of coupled exchange fluxes. An experiment aiming to simultaneously collect the data of the boundary layer and soil evaporation in the area around the soil–atmosphere interface was conducted using a wind tunnel integrated with a soil tank. The turbulent boundary layer above the undulating soil surface was captured using high‐resolution hot‐wire anemometry, confirming the presence of recirculation zones in the undulating valleys and locally low evaporative flux. Experimental data were used to validate the extended model, and modeling results demonstrate that turbulent airflow enhances evaporation and shortens the duration of Stage I. The surface geometry significantly affects the local evaporative flux by influencing the vapor distribution, concentration gradient, and water availability at the soil surface, especially when recirculation zones form in the valleys. As a joint result of turbulence and surface undulations, the influence of wind speed on both the local and system‐level evaporation rate is restricted.
more » « less- NSF-PAR ID:
- 10448353
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Vadose Zone Journal
- Volume:
- 19
- Issue:
- 1
- ISSN:
- 1539-1663
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract r a), in which the vapor roughness length (z 0v) and the momentum roughness length (z 0m) are two crucial parameters that describe the transport near the soil‐atmosphere interface. Typically, when evaluating bare‐soil evaporation, these two characteristic lengths are assumed equal, although differences are likely to occur especially in turbulent flows over undulating surfaces. Thus, this study aims to investigate the relationship betweenz 0vandz 0mabove undulating surfaces to ultimately improve accuracy in estimating evaporation rate. To achieve this goal, four uniquely designed wind tunnel—soil tank experiments were conducted considering different wind speeds and undulation spacings. Particle image velocimetry (PIV) was used to measure the velocity field above the undulating surface in high resolution. Using the high‐fidelity data set, the logarithmic ratio ofz 0vtoz 0mis determined and used to estimater a. Results confirm that these lengths differ significantly, with the logarithmic ratio roughly ranging from −15 to −5 under the conditions tested. PIV‐measured results demonstrate this ratio is closely tied to the mass and momentum transport behaviors influenced by surface undulations. Using the data‐integrated formulation ofr a, predictions of evaporation rate were prepared for both the laboratory and lysimeter experiments, demonstrating the efficacy of the proposed approach in this study. -
more » « less
Abstract Drying front propagation and coupled heat and mass transfer analysis from porous media is critical for soil–water dynamics, electronics cooling, and evaporative drying. In this study, de-ionized water was evaporated from three 3D printed porous structures (with 0.41 mm, 0.41 mm, and 0.16 mm effective radii, respectively) created out of acrylonitrile butadiene styrene (ABS) plastic using stereolithography technology. The structures were immersed in water until all the pores were invaded and then placed on the top of a sensitive scale to record evaporative mass loss. A 1000 W/m2 heat flux was applied with a solar simulator to the top of each structure to accelerate evaporation. The evaporative mass losses were recorded at 15 min time intervals and plotted against time to compare evaporation rates from the three structures. The evaporation phenomena were captured with a high-speed camera from the side of the structures to observe the drying front propagation during evaporation, and a high-resolution thermal camera was used to capture images to visualize the thermal gradients during evaporation. The 3D-structure with the smallest effective pore radius (i.e., 0.16 mm) experienced the sharpest decrease in the mass loss as the water evaporated from 0.8 g to 0.1 g within 180 min. The designed pore structures influenced hydraulic linkages, and therefore, evaporation processes. A coupled heat-and-mass-transfer model modeled constant rate evaporation, and the falling rate period was modeled through the normalized evaporation rate.
-
more » « less
Abstract Fundamental process understanding and description of heat, mass, and momentum exchanges across the land‐atmosphere interface in model boundary forcing parameterizations is critical to the simulation of near‐surface soil moisture dynamics (e.g., bare‐soil evaporation). This study explores the sensitivity of a continuum‐scale porous media heat and mass transfer model to the spatial‐discretization length‐scales (i.e., spatial‐resolution) of near‐surface atmospheric data; the goal is to determine how much data are needed to force the model and adequately capture evaporative water losses and subsurface state variable distributions. The requisite atmospheric forcing data were taken from the high‐resolution, precision bare‐soil evaporation experiments of Trautz et al. (2018,
https://doi.org/10.1029/2018WR023102 ). Simulation results demonstrated that shallow subsurface mass and heat transfer dynamics can be adequately captured with forcing data averaged over large length‐scales, or a minimal number of measurements, provided that soil conditions are properly described. The soil moisture spatial distributions were found to be insensitive to horizontal variations in the forcing data. The model failed to capture small‐scale trends observed experimentally; this did not impact the accuracy of total evaporative water loss estimates however. These results indicate that in future physical experimental efforts conducted at 1–10‐m length‐scales, there is no need to focus on the generation of high‐spatial resolution atmospheric measurements—time and effort would be better spent in characterizing soil conditions and properties. Even though a theoretical foundation was not provided to directly extrapolate this work to the field scale, these findings have practical value in designing field data collection strategies. -
Flow separation caused by steep topography remains a significant obstacle in accurately predicting turbulent boundary-layer flows over complex terrain, despite the utilization of sophisticated numerical models. The addition of atmospheric thermal stability, in conjunction with steep topography, further complicates the determination of disrupted turbulent wind patterns. The turbulent separated flows over a two-dimensional (2D) steep hill under thermal stratification has not been extensively addressed in previous experimental studies. Such measurements are crucial for enhancing our comprehension of flow physics and validating numerical models. We measured the turbulent wind flows over a 2D steep hill immersed in a stable boundary layer (of the bulk Richardson Number = 0.256) in a thermally-stratified boundary-layer wind tunnel. The flow separation, re-circulation zone and flow reattachment were characterized by the planar particle image velocimetry technique. Vertical profiles of mean air temperature and its fluctuations are also quantified at representative locations above the 2D steep hill and in the near wake region. Results indicate that the separated shear layer, initiated near the crest of the 2D steep hill, dominates the physical process leading to high turbulence levels and the turbulent kinetic energy production in the wake region for both stable and neutral thermal stability. Although the stable boundary layer does not dramatically change the turbulent flow pattern around the hill, the mean separation bubble is elongated by 13%, and its vertical extent is decreased by approximately 20%. Furthermore, the reduced turbulence intensities and turbulent kinetic energy of the near wake flow are attributed to the relatively low turbulence intensity and low momentum of the stable boundary layer due to buoyancy damping, compared to the neutral boundary layer. Additionally, a distinct low-temperature region—a cold pool—is extended beyond the separation bubble, reflecting the significant sheltering effect of the 2D steep hill on the downwind flow and temperature field.more » « less
-
more » « less
Abstract An experimental apparatus was designed to study the impacts of wettability on evaporation of water from Ottawa sand. Evaporation rates were measured for: (1) a 5.7-cm-thick layer of hydrophilic Ottawa sand; (2) a 5.7-cm-thick layer with 12% hydrophobic content, consisting of a 0.7-cm-layer of n-Octyltriethoxysilane-coated hydrophobic sand buried 1.8 cm below the surface of hydrophilic sand; and (3) a 5.7-cm-thick layer with mixed wettabilities, consisting of 12% n-Octyltriethoxysilane-coated hydrophobic sand mixed into hydrophilic sand. The sand–water mixtures experienced forced convection above and through the sand layer, while a simulated solar flux (i.e., 112±20 W/m2) was applied. Evaporation from homogeneous porous media is classified into the constant-rate, falling-rate, and slow-rate periods. Wettability affected the observed evaporation mechanisms, including the transition from constant-rate to falling-rate periods. Evaporation entered the falling-rate period at 12%, 20%, and 24% saturations for the all hydrophilic sand, hydrophobic layer, and hydrophobic mixture, respectively. Wettability affected the duration of the experiments, as the all hydrophilic sand, hydrophobic layer, and hydrophobic mixture lasted 17, 20, and 26 trials, respectively. Both experiments with hydrophobic particles lasted longer than the all hydrophilic experiment and had shorter constant-rate evaporation periods, suggesting hydrophobic material interrupts capillary action of water to the soil surface and reduces evaporation. Sand temperatures suggest more evaporation occurred near the test section inlet for higher saturations and the hydrophobic layer experienced more evaporation occur near the outlet. Evaporation fluxes were up to 12× higher than the vapor diffusion flux due to enhanced vapor diffusion and forced convection.
