## Abstract

We perform a numerical study of the longevity of a super-hydrophobic surface (SHS) in under-saturated, stationary liquid. We numerically solve the spatial-temporal evolution of the gas concentration in the liquid, the time-variation of mass flux of gas out of the plastron, as well as the time required for the gas in the plastron to be fully dissolved (i.e., the plastron lifetime). We find that the profiles of gas concentration at different times are self-similar, and the mass flux reduces with time (t) at a rate of 1/t^{0.5}. In addition, we examine the impact of texture parameters, including pitch, gas fraction, texture height, and advancing contact angle, on the diffusion process. Our results show that both plastron lifetime and diffusion length increase with increasing the gas fraction or increasing the texture height and are independent of the advancing contact angle and pitch. We propose simple analytical models for plastron lifetime and diffusion length. We show that the model has a fair agreement with the experimental data reported in the literature, and can predict the longevity for SHS with various texture geometries, texture sizes, and under different degrees of under-saturations. Our models could guide the design of long-life SHS for underwater applications such as reducing skin-friction drag and preventing biofouling.

## 1 Introduction

The super-hydrophobic surface (SHS) has received growing attention in recent years due to its potential for wide engineering applications. When submerging in liquid, the SHS can entrap a layer of gas bubbles (or plastrons) between the solid surface and the liquid, forming the so-called Cassie-Baxter state [1]. The gas layer could protect underwater surfaces from corrosion and biofouling [2,3]. Moreover, the gas layer could support a shear-free boundary, resulting in a reduction of skin frictions in both laminar [4–8] and turbulent flows [9–18]. However, implementing the SHS in real engineering systems is still a challenge [19]. One main reason is that the entrapped gas on the SHS could sometimes be replaced by the liquid, leading to a transition from the Cassie-Baxter state to the Wenzel state [20] (i.e., a wetting transition). A wetted SHS loses the most benefits [2,21]. In addition to pressure [22–29] and flow-induced shear stress [30,31], gas dissolution into under-saturated liquid [32,33] is one of the main reasons causing the wetting transition. This study aims to better understand this gas dissolution process, and develop an analytical model to predict the SHS longevity in under-saturated liquid.

Over the past decade, the dissolution of gas from an SHS into the under-saturated liquid has been extensively studied through a combination of numerical, experimental, and analytical approaches. In early numerical studies [34–38], the rate of mass transfer from the plastron to surrounding liquid (i.e., the speed of plastron loss) was modeled by assuming an *invasion coefficient* of gas [39]. As a result, the mass transfer rate was found to be constant during the depinning stage [23] where the interface slides down along the sidewalls with a constant shape. Based on this model, the impacts of liquid pressure [34,36,38] and texture geometry [34,35,38] on the plastron longevity were investigated. Later, Kadoko et al. [40] calculated the rate of mass transfer analytically by approximating the gas dissolution problem in one dimension and solving the dissolved gas concentration in the liquid. Different from the previous numerical studies, Kadoko et al. found that the mass transfer rate reduces during the depinning stage. They examined the plastron longevity as a function of texture height and initial gas concentration in the liquid.

In the experimental studies, the wetting transition induced by the gas dissolution has been observed through various optical techniques, such as Total-Internal-Reflection [41–44], Confocal Microscopy [45–47], X-Ray Diffraction [48], and direct imaging [49,50]. The plastron longevity was measured either directly by tracking the movement of the gas–liquid interface [45] or indirectly based on the intensity of the light reflected from the surface [41]. Longevity for SHSs with both regular and irregular texture geometries has been investigated. Most experimental data showed that the plastron longevity decreases with increasing the hydrostatic pressure [43–45] or reducing the dissolved gas in the liquid [50] (in both cases, the degree of liquid under-saturation increases). For low hydrostatic pressures and small texture size, an equilibrium state with zero mass flux can be achieved [49,51].

*t*for an SHS with regular pores in under-saturated liquid

_{D}where *h* is the height of the texture, *L _{D}* is the diffusion length,

*D*is the diffusion coefficient,

_{G}*c*

_{0}and

*c*

_{∞}are the initial gas concentration in the plastron and in the liquid, respectively. By fitting the experimental data, Lv et al. Reported

*L*= 450

_{D}*μ*m. They showed that when normalizing the time by

*t*, the movements of the gas–liquid interface under different pressures collapse. Recently, Kim and Park [52] reported a much smaller diffusion length

_{D}*L*= 4 to 6

_{D}*μ*m for a different SHS, suggesting that

*L*may vary with texture geometry and liquid pressure. By assuming a known invasion coefficient, Hemeda and Tafreshi [38] derived equations of the plastron longevity in two regimes, defined by pressure smaller or larger than the critical pressure for interface to depin from the edge of the texture.

_{D}The acceleration of the gas dissolution by external flows has also been studied intensively. For example, by using alternate shear-free and no-slip boundaries [53,54] or a wall with a prescribed slip coefficient [55], several researchers numerically solved the coupled mass and momentum transport equations in laminar flows over SHS. By using an approximate integral method, Barth et al. [56,57] developed first-principles models of the convective gas dissolution in laminar and turbulent flows over SHS. In experimental studies, by measuring the total-internal-reflection light from the SHS, several researchers [58–61] reported a reduction of plastron longevity with increasing *Reynolds number*. By using fluorescence microscopy [53], the spatial distributions of gas concentration were measured in laminar flows over SHS. Recently, by tracking the movement of the gas–liquid interface, scaling laws [52,62,63] of *Sherwood number* were proposed to quantify the convective enhancement of gas dissolution.

Here, we focus on gas dissolution in stationary liquid. We aim to propose a general formula for the plastron longevity, which does not require a per-knowledge of the diffusion length or invasion coefficient. We also aim to propose a model for the diffusion length, which can explain the different values of *L _{D}* reported in the literature. To achieve these goals, we numerically solve the spatial and temporal evolutions of the gas concentration in the liquid. The mass transfer rate and plastron lifetime are calculated based on the resolved gas concentration field near the interface. We will show that as gas being slowly dissolved into the liquid, the gas concentration in the liquid increases and the mass transfer rate reduces. In addition, we examine the impacts of texture size, texture shape, and interface shape on plastron longevity. Last, we propose analytical models for the plastron lifetime and diffusion length, which can work for SHS with various texture geometries, texture sizes, and under different degrees of under-saturations.

## 2 Numerical Methods

### 2.1 Problem Description.

In this study, we focus on the gas transfer from a SHS with microgrooves to the under-saturated liquid, as shown in Fig. 2(a). The geometry of the groove is characterized by three parameters, groove pitch *λ*, groove width *w*, and groove height *h*. The SHS is installed on the bottom wall of a closed chamber, which has a height of *H = *100 mm and is filled with water. The gas trapped in the plastron is air, with a density of ρ_{gas,0}* *= 1.2 kg/m^{3} (evaluated at atmosphere pressure). The diffusion coefficient of the air in the water is *D _{G} *= 2.0 × 10

^{−9}m

^{2}/s. The water is initially saturated with air at atmosphere pressure so that

*c*

_{∞}

*= 0.023 kg/m*

^{3}. The pressure in the water above the SHS is denoted as

*P*.

_{L}When *P _{L}* is equal to the atmosphere pressure, the gas–liquid interface is flat and there is no gas flux. As increasing

*P*, the interface becomes curved,

_{L}*c*

_{0}increases, and the gas dissolution takes place when

*c*

_{0}is larger than

*c*

_{∞}. Depending on whether

*P*is larger than the critical pressure for the interface to depin, the wetting transition can experience either a two-stage or a one-stage process [34]. When

_{L}*P*is smaller than the critical pressure, the wetting process consists of two stages: an initial sagging and a depining impalement. During the first stage of initial sagging, the interface increases its curvature while pinning to the tip of the texture. During the second stage of depining impalement, the interface detaches from the tip of the texture, slides down along the sidewalls of the groove, and maintains a constant shape. The interface depin only after the local contact angle reaches the advancing contact angle

_{L}*θ*

_{adv}. For SHS with single-length scale texture,

*θ*

_{adv}is typically smaller than 120 deg [23,24,45]. For SHS with double-length scale hierarchical texture,

*θ*

_{adv}can exceed 120 deg [63–65]. On the other hand, when

*P*is larger than the critical pressure, the wetting process only includes the second stage of depining impalement.

_{L}In this study, we consider only the depining impalement for simplicity. We set *c*_{0}* *= 1.2*c*_{∞}* *= 0.028 kg/m^{3}, resembling a case with *P _{L}* ≈ 1.2 atm, much larger than the critical pressure for the interface to depin. Although ignoring the initial sagging, our simulations provide an approximation of the plastron lifetime for SHS with larger ratios of groove height to width. As listed in Table 1, we vary the pitch, ratio of groove width to pitch, and ratio of groove height to pitch to study the impact of texture geometry on the gas dissolution. In addition, we vary

*θ*

_{adv}from 90 deg to 150 deg to study the impact of surface hydrophobicity and hierarchical texture on gas dissolution.

### 2.2 Numerical Model.

To simulate the gas dissolution process described above, we use a simplified model shown in Fig. 2(b). First, we assume that the rate of gas dissolution is identical for different grooves on the SHS. Thus, instead of solving for the entire surface, the model only simulates the gas transfer for a single groove and uses a symmetric boundary condition on the two sides of the domain. This assumption is valid when the surface is much larger than the diffusion length. Second, considering that the displacement of the gas–liquid interface is negligible compared to the height of simulation domain, we assume that the interface is fixed in space during the wetting process. To validate this assumption, we performed two simulations where the interface is fixed either at the tip of the groove or at the bottom of the groove. We find that the plastron longevities calculated from the two simulations differ by less than 1%.

*x*and

*y*as the Cartesian coordinates (

*y*is the direction normal to the surface), and

*t*is the time. We set

*y*=

*0 at the tip of the groove, and*

*x*=

*0 at the middle of the groove. The evolution of dissolved gas concentration in the liquid*

*c*(

*x, y, t*) is governed by the following equation:

It should be noted that replacing this Neumann boundary condition to Dirichlet boundary condition, *c*(*y *=* H*) = *c*_{∞}, would not change the simulation results. The reason is that the height of the simulation domain is sufficiently large such that at the end of the simulations (i.e., when all the gas in the plastron has been transferred into the liquid), the gas concentration at *y *=* H* remains unchanged.

*c*

_{0}at the gas–liquid interface

The reason a fixed gas concentration is used at the interface is because the gas pressure in the plastron remains as a constant during the depining impalement.

*J*

_{out}, can be calculated by based on the Fick's first law

*A*=−

*w*(

*θ*

_{adv}−π/2)/cos

*θ*

_{adv}is the interface area per unit depth into the page, ∂

*c*/∂

*n*is the gradient of gas concentration normal to the interface, and the integral is evaluated along the interface

*S*. The mass of gas remaining in the plastron (per unit depth into the page), denoted as

*m*, is calculated as

*m*

_{0}= ρ

_{gas,0}

*wh*is initial mass of gas in the plastron. We define and calculate the plastron lifetime

*t*as the time when

_{f}*m*is equal to zero

The above-mentioned governing equation with the prescribed boundary and initial conditions is numerically solved in COMSOL Multiphysics 5.4. For simulation domain near the surface (0<*y *<* *5 mm), a small square element of size 2 *μ*m × 2 *μ*m is used. For the domain near the curved interface (*y *<* *0), a triangle element of size 2 *μ*m is used. Far away from the surface (*y *>* *5 mm), a larger element is used. A time interval ranging from 4 to 2500 s is used depending on the plastron lifetime (a larger time interval is selected for cases with a larger plastron lifetime). The number of time-step is about 700. Further reducing element size or time interval does not change the plastron lifetime by more than 3%.

*J*

_{ave}as

*L*is the diffusion length which indicates how far the gas in the plastron propagates into the liquid at the end of wetting process. Solving this equation, the diffusion length can be expressed as

_{D}Note, Eq. (12) is very similar to the analytical model of characterize time scale proposed by Lv et al. [45] shown in Eq. (1). When the diffusion length is known, the plastron lifetime can be calculated from Eq. (12). We attempt to understand what parameters impact *L _{D}*, propose an analytical model for

*L*as well as for

_{D}*t*.

_{f}## 3 Results and Discussion

### 3.1 A Typical Case.

First, we present the simulation results of a typical gas dissolution process for a grooved SHS with *λ *= 100 *μ*m, *w/λ *= 0.5, *h/λ *= 1.0, and *θ*_{adv}* *= 120 deg. Figure 3 shows the spatial distributions of gas concentration near the SHS at the early stage of wetting transition. At the very beginning *D _{G}t*/

*λ*

^{2}< 10, only the gas concentration near the gas–liquid interface is impacted, while the gas concentration far away from the plastron is unchanged. As time goes on, the gas from the plastron propagates further into the liquid in both

*x*and

*y*directions, and the gas concentration in the liquid further away from the surface increases. For

*D*/

_{G}t*λ*

^{2}> 10, the gas concentration near the surface is nearly uniform along the

*x*-direction. The gas propagation can be approximated in one dimension.

Figure 4(a) shows *x*-averaged gas concentration as a function of *y* at three different times, *t *=* *0.2, 0.5, and 0.8*t _{f}*. As expected, the gas concentration changes from

*c*

_{0}to

*c*

_{∞}as moving from the surface to the far-field of the liquid. Moreover, as increasing time, the

*y*-gradient of gas concentration reduces, which indicates a reduction of mass transfer rate during the wetting transition. Figure 4(b) shows the same profiles in Fig. 4(a) plotting as a function of

*y*/(

*D*)

_{G}t^{0.5}. Interestingly, the three profiles collapse and agree well with the analytical solution 1−erf[

*y*/2(

*D*)

_{G}t^{0.5}]. This result implies that the concentration profiles at different times are self-similar and ∂

*c*/∂

*y*reduces at a rate of 1/

*t*

^{0.5}.

Figure 4(c) shows the time-variation of the mass of gas remaining in the plastron. For comparison, we plot the solution by Kadoko et al. [40], who approximated the gas dissolution in one dimension and analytically solved the dissolved gas concentration in the liquid. Kadoko et al. also considered the initial sagging and the change of gas pressure in the plastron. Although our simulation ignored the initial sagging, the trend of mass variation shows a good agreement with their result. As increasing time, the mass in the plastron reduces and reaches to 0 at *t *= *t*_{f}. Moreover, as increasing time, the rate of mass loss decreases. As shown in Fig. 4(d), we find that *J*_{out} reduces at a rate of 1/*t*^{0.5}, which is consistent to the self-similar behavior of the gas concentration profiles. This trend is expected since as increasing time, the gas concentration in the liquid increases, and the gradient of gas concentration near the SHS reduces.

It should be emphasized that the current result is different from previous numerical studies [34–38] which assumed a constant invasion coefficient and found a constant mass flux during the depining wetting stage. And the current result only applies to cases where the liquid is stationary. When the SHS is exposed to external flows, a concentration boundary layer could be developed over the SHS [53], and the gas concentration in the boundary layer may not vary during the wetting transition.

*J*

_{out}:

*G*is a nondimensional parameter depending on four texture parameters:

*λ/H*,

*w/λ*,

*h/λ*, and

*w*/

*A*. Here,

*λ/H*denotes the ratio of texture size to domain height,

*w/λ*and

*h/λ*are related to the texture shape, and

*w*/

*A*= (cos

*θ*

_{adv})/(π/2-

*θ*

_{adv}) is the ratio of groove width to interface area depending only on

*θ*

_{adv}. For flat interface,

*θ*

_{adv}

*= 90 deg and*

*w*/

*A*=

*1. For curved interface,*

*θ*

_{adv}

*>*90 deg and

*w*/

*A*<

*1. The larger the*

*θ*

_{adv}, the smaller the

*w*/

*A*. Substituting Eq. (13) into Eqs. (9) and (12), the following two scaling relationships for the plastron lifetime and the diffusion length can be found

### 3.2 Effects of Texture Parameters on Gas Diffusion.

First, we run a series of simulations where *λ/H* varies, and *w*/*λ*, *h*/*λ*, and *θ*_{adv} are fixed. This resembles a situation where SHSs of different texture sizes are submerged in water at fixed depth. Figure 5(a) shows the time variations of *J*_{out} for different *λ/H*. Interestingly, regardless of *λ/H*, the profiles of *J*_{out} for all cases overlap. This result indicates that *G* is independent of *λ/H*. Therefore, according to Eqs. (14) and (15), we expect scaling relations: *t _{f}* ∼(

*λ/H*)

^{2}, and

*L*∼(

_{D}*λ/H*), which are confirmed by Figs. 5(b)–5(c).

Second, we modify *w*/*λ* from 0.1 to 0.9, and fix the values of *λ/H*, *h*/*λ*, and *θ*_{adv}. This means that the fraction of surface area covered by gas (or the gas fraction *ϕ _{g}*) varies. For most SHS,

*ϕ*is larger than 0.5. A larger

_{g}*ϕ*is preferred for drag reduction [27]. Figure 6(a) shows the time variations of

_{g}*J*

_{out}for different

*w*/

*λ*. It is found that with increasing

*w*/

*λ*,

*J*

_{out}reduces. This trend is expected since as increasing

*w*/

*λ*, more liquid is exposed to the plastron, the gas concentration in the liquid can adjust more quickly to

*c*

_{0}, and leading to a faster reduction of the normal-gradient of gas concentration. Moreover, as shown in Fig. 6(b), the profiles of

*J*

_{out}×(

*w*/

*λ*) collapse. It suggests that the

*J*and

_{out}*G*are inversely proportional to the gas fraction, i.e.,

*J*

_{out}∼(

*w/λ*)

^{−1},

*G*∼(

*w/λ*)

^{−1}. Thus, according to Eqs. (14) and (15), we expect scaling relations:

*t*∼(

_{f}*w*/

*λ*)

^{2}and

*L*∼

_{D}*G*

^{−2}∼(

*w*/

*λ*)

^{2}, which are confirmed in Figs. 6(c)–6(d).

Third, we run a number of simulations where *h*/*λ* varies from 0.2 to 4, and *λ/H*, *w*/*λ*, and *θ*_{adv} remain as constants. This means that the texture height varies, while all other parameters are kept same. Figure 7(a) shows the time variations of *J*_{out} for these cases. As shown, the profiles of *J*_{out} for all cases overlap, indicating that *J*_{out} and *G* are independent of *h*/*λ*. This result is expected since *h* is not included in the boundary conditions for solving the dissolved gas concentration in the liquid. As a result, according to Eqs. (14) and (15), we expect scaling relations: *t _{f}* ∼(

*h*/

*λ*)

^{2}and

*L*∼

_{D}*h*/

*λ*, which are confirmed by Figs. 7(b)–7(c).

Last, we perform simulations where *θ*_{adv} varies from 90 deg to 150 deg, and *λ/H*, *w*/*λ*, and *h*/*λ* are fixed. Larger *θ*_{adv} can be obtained by either coating the surface with more hydrophobic materials or adding nanoscale roughness on the side walls of the texture. With increasing *θ*_{adv}, the interface is more curved, and *A*/*w* increases. Figure 8(a) shows the time variations of *J*_{out} for these cases. As shown, *J*_{out} reduces slightly with increasing *θ*_{adv}. Similar to increasing *w*, increasing *θ*_{adv} effectively increases the contact area between the liquid and plastron. Thus, the liquid can adjust more quickly to *c*_{0}, leading to smaller *J*_{out}. Moreover, as shown in Fig. 8(b), when plotting *J*_{out}×(*A*/*w*), the profiles collapse, indicating that *J*_{out} and *G* are linearly proportional to *w*/*A*. As a result, the mass transfer rate (*AJ*_{out}) is same for different *θ*_{adv}. And according to Eqs. (14) and (15), we expect that *t _{f}* and

*L*are independent of

_{D}*θ*

_{adv}, which are confirmed in Figs. 8(c)–8(d). Although increasing

*θ*

_{adv}allows the SHS to sustain higher pressures [22], it may not impact the SHS longevity against gas diffusion.

### 3.3 Analytical Models of Plastron Lifetime and Diffusion Length.

*G*only depends on

*w/λ*and

*w*/

*A*, and is independent of

*λ/H*and

*h*/

*λ*

where the constant 0.81 is obtained by fitting the simulation results. We replace *w/λ* by *ϕ _{g}*, since gas fraction is a more general term for characterizing the texture geometry of SHS. As shown in Figs. 9(a) and 9(b), the plastron lifetimes and diffusion lengths obtained from all simulations in this study can be predicted by Eqs. (17) and (18), although these SHSs have different magnitudes of

*λ*,

*w*,

*h*, and

*θ*

_{adv}.

*ϕ*→ 1 and considering a flat interface, the gas diffusion problem can be approximated as one-dimensional (1D), where the gas concentration

_{g}*c*and the mass flux of gas out of plastron

^{1D}*J*

_{out}

*can be expressed as*

^{1D}Thus, the models of plastron lifetime and diffusion length provided in Eqs. (17) and (18) are very similar to those obtained for 1D diffusion. The only difference is that the models for SHS have an additional term of *ϕ _{g}*

^{2}compared to the models for 1D diffusion where

*ϕ*= 1. The two constants 0.81 and π/4 differ by less than 5%.

_{g}The model provided in Eq. (17) can be readily used to predict the longevity of SHS in under-saturated liquid if the gas fraction, roughness height, and initial gas concentrations in the liquid and plastron are known. A comparison between the experimental measured *t _{f}* and model predicted

*t*is shown in Fig. 10. The experimental data are taken from five different studies. Since Lv et al. [45] measured the plastron lifetime based on the time when the interface reaches to 0.72

_{f}*h*(

*h*=

*40*

*μ*m for their SHS), we thus use

*h*=

*28.8*

*μ*m in Eq. (17) to estimate

*t*. For the randomly roughed SHS used by Hokmabad and Ghaemi [60], we use the reported mean peak to trough height as

_{f}*h*, and a gas fraction of 0.6 to 0.9 to estimate the range of

*t*. For the randomly roughed SHS studied in Poetes et al. [43], we use the reported plastron thickness of a fresh immersed SHS as

_{f}*h*, and a gas fraction of 0.6–0.9 to estimate the range of

*t*. As shown in Fig. 10, regardless of the texture geometry, there is a fair agreement between the predictions and measurements for most cases. The mismatch between prediction and measurement could attribute to the measurement uncertainty, irregular textures, and interface instabilities in experiments [45]. In particular, the measurement based on light reflected from the SHS [41,43,60] may underestimate

_{f}*t*, since light reflected from a curved or depined interface may not be detectable by the sensor. Interface instability caused by surface energy minimization could trigger a sudden collapse of the plastron [45].

_{f}The diffusion length model in Eq. (18) explains how *L _{D}* varies with the degree of under-saturation and the texture parameters. First, the model predicts that

*L*reduces when increasing the degree of under-saturation. This is expected since as increasing

_{D}*c*

_{0}−

*c*

_{∞}, per amount of liquid can absorb more gas, less liquid will be impacted by the plastron, so that the gas propagates for a shorter distance. Second, our model predicts that

*L*increases with increasing the texture height and the gas fraction. This is also expected since as increasing texture height or gas fraction, the amount of gas in the plastron increases, more liquid is required to absorb all the gas, so that the gas propagates for a larger distance. For the experimental conditions in Lv et al. [45], the model predicts an

_{D}*L*varying from 400 to 1800

_{D}*μ*m depending on the pressure, which is close to their reported

*L*= 450

_{D}*μ*m. For the experimental conditions in Kim and Park [52],

*L*predicted by our model is two orders of magnitude larger than the measured 4–6

_{D}*μ*m. This disagreement is probably because the wetting transition observed in their experiments was impacted not only by the gas diffusion but also by the interface instability induced by surface energy minimization. In fact, for the SHS used by Kim & Park, the Cassie-Baxter state has higher interfacial energy than the Wenzel state and is thus not thermodynamically stable [66].

Note, in the current study, *H* is sufficiently large so that all gas in the plastron is dissolved in the liquid. However, when reducing *H*, the liquid could be fully saturated before the gas is completely dissolved, and the plastron lifetime can be infinite [40]. We expect that the proposed model of plastron lifetime holds true as long as *H* is sufficiently larger than *L _{D}* so that the gas concentration in the .far-field is not affected by the plastron. To examine the limit of

*H*, we performed simulations where

*H*varies from 10 to 100 mm, and the groove dimensions are fixed (λ

*= 0.1 mm,*

*w*/λ

*= 0.5,*

*h*/λ

*= 1). Figure 11 shows*

*t*as a function of

_{f}*H/L*. As expected, for

_{D}*H*>

*3*

*L*, the model predicts the plastron lifetime very well, and the plastron lifetime does not change with

_{D}*H*. But when

*H*<

*3*

*L*, the model under-estimates the plastron lifetime. The limit where the proposed model is valid is

_{D}*H*>

*3*

*L*(

_{D}*L*can be predicted by Eq. (18)).

_{D}## 4 Conclusions

In this work, we performed a numerical study of the longevity of super-hydrophobic surfaces fully submerged in under-saturated liquid. We numerically solved the spatial and temporal variations of the dissolved gas concentration in the liquid as it slowly increases during the wetting process. Based on the resolved gas concentrations and Fick's first law, we calculated the rate of mass transfer, the time-variation of the mass remaining in the plastron, and the time required for all gas to be dissolved (i.e., the plastron lifetime). We found that the gas concentration profiles at different times display a self-similarity, and the mass flux reduces with time (*t*) at a rate of 1/*t*^{0.5}. These observations agree with the one-dimensional gas diffusion and could be attributed to the small texture size compared to the diffusion length.

Furthermore, we performed a series of numerical simulations to examine the impacts of texture parameters, including pitch, gas fraction, texture height, and advancing contact angle, on the mass flux, plastron lifetime, and diffusion length. Based on the simulation results, we proposed simple analytical model for the plastron lifetime, which does not require a per-knowledge of diffusion length or invasion coefficient. Our model showed that the SHS longevity increases with increasing the gas fraction or increasing the texture height and is independent of the advancing contact angle and texture pitch. We found that the proposed model agrees fairly well with the experimental data reported in the literature, and can predict the longevity of SHS with various texture sizes, texture geometries, and under different degrees of under-saturations. We also proposed an analytical model for the diffusion length, which showed that the diffusion length increases with decreasing the degree of under-saturation, increasing gas fraction, or increasing the texture height. We believe our work can provide a general guideline for designing SHS with extended lifetimes when submerged in under-saturated liquid.

It should be noted that the proposed models and conclusions in this study only apply to these wetting processes which mainly consist of the depining impalement, not the sagging. These wetting processes are often observed when the pressure above the surface is larger than the critical pressure for interface to depin from the tip of texture. If the sagging stage is mainly involved, the effect of surface tension needs to be accounted. As interface curvature increases during the sagging, the gas pressure in the plastron reduces due to the surface tension, which could result in a lower mass flux. Although ignoring the sagging, the proposed model could provide a good approximation for the longevity of these SHSs which have a large ratio of texture height to width.

## Acknowledgment

We thank the support by UMass Dartmouth's Marine and Undersea Technology (MUST) Research Program funded by the Office of Naval Research (ONR) under Grant No. N00014-20-1-2170. We also thank the support by National Science Foundation (NSF) under Grant No. 2041479.

## Funding Data

Office of Naval Research (ONR) (Grant No. N00014-20-1-2170; Funder ID: 10.13039/100000006).

National Science Foundation (NSF) (Grant No. 2041479; Funder ID: 10.13039/100000001).

## Nomenclature

*A*=area of gas-liquid interface per unit depth into the page

*c*=gas concentration in the liquid

*c*_{∞}=initial gas concentration in the liquid

*c*_{0}=initial gas concentration in the plastron

*c*=^{1D}gas concentration in the liquid for one-dimensional gas diffusion

*D*=_{G}diffusion coefficient

*G*=A non-dimensional parameter depending on SHS texture parameters

*h*=groove height of the SHS texture

*H*=height of the water column above the SHS

*I*=all variables should appear in italics

*J*_{ave}=time-average mass flux

*J*_{out}=mass flux averaged over gas-liquid interface

*J*_{out}=^{1D}mass flux averaged over gas-liquid interface for 1D gas diffusion

*L*=_{D}diffusion length scale

*L*=_{D}^{1D}diffusion length scale for 1D gas diffusion

*m*=mass of gas in the plastron

*P*=_{L}pressure in the water above the SHS

*S*=gas-liquid interface

*t*=time

*t*=_{D}characterize time for a SHS in under-saturated liquid

*t*=_{f}time when the mass in the plastron reduces to zero (plastron lifetime)

*t*=_{f}^{1D}plastron lifetime for 1D gas diffusion

*w*=groove width of the SHS texture

*x*=Cartesian coordinate, direction parallel to the SHS

*y*=Cartesian coordinate, direction normal to the SHS

*λ*=groove pitch of the SHS texture

- ρ
_{gas,0}= density of gas evaluated at atmosphere pressure

*θ*_{adv}=advancing contact angle of the SHS

*ϕ*=_{g}fraction of surface area covered by gas on the SHS (gas fraction)