1. Introduction
Bad weather such as typhoons, blizzards, and freezing rain will bring huge disasters to human beings. It is of great significance to study its genesis and evolution mechanism. In the early warning observation of typhoon disasters, every 1 km of typhoon track forecast error can reduce losses by USD 20 million, and every 1 m/s reduction in intensity forecast error can reduce losses by USD 100 million [
1,
2]. The refined detection of severe weather in the inner core area is an important means to predict its intensity and development trend [
3,
4]. As a key parameter to determine the duration of a typhoon, high-performance measurement of humidity in the inner core region is particularly important [
5,
6].
The internal environment of severe weather is complex, with wind speeds as high as 50 m/s, and temperature and air pressure as low as −90 °C and 10 hPa [
7,
8,
9,
10]. On the one hand, freezing temperatures and sudden changes in air pressure will cause the sensor to condense and repeatedly freeze [
11], adding massive interference to the sensor measurement. We urgently need to improve the sensitivity of the sensor to detect weak changes in the sensor in extreme environments and reduce measurement errors. On the other hand, extremely high wind speeds and drastic changes in the background also put high requirements on the dynamic performance of the sensor [
12,
13,
14]. The state of the environment is constantly changing, and we need sensors that can respond in real-time to sudden changes in the background to obtain all the conditions of the environment and ensure the comprehensiveness of detection. Therefore, it is essential to simultaneously improve sensors’ static sensitivity and dynamic response performance in complex environments.
Humidity sensors mainly include capacitive [
15,
16], resistive [
17], piezoresistive [
18,
19], resonant [
20], and fiber-optic [
21,
22] components. As the main research object and commercial product [
23,
24,
25,
26,
27], the capacitive humidity sensor has the advantages of high sensitivity, fast response speed, good product interchangeability, and easy integration. In the humidity measurement process of capacitive sensors, the adsorption of water molecules by the humidity sensitive medium will cause the dielectric constant to change and thus the capacitance value to change. Capacitive structures of sensitive sensing units include sandwich type [
28] and interdigital type [
29]. The two electrodes of the sandwich capacitor are facing each other, and the moisture-sensitive film is located between the electrodes, which has high sensitivity and slow dynamic response. The electrodes of the interdigital capacitor are in the same plane, and the surface of the electrode is covered with a moisture-sensitive film, which has low sensitivity and a fast dynamic response.
To improve the sensitivity and dynamic response performance of humidity sensors in extreme environments, researchers have made a lot of technical research. In 2014, Luo designed a sandwich capacitive humidity sensor with a heater, which solved the sensor condensation problem by heating. At the same time, the measurement circuit is designed according to the characteristics of the humidity sensor and ground experiments, and the results show that the humidity sensitivity is about 2.3
[
30]. In 2016, Keyvan proposed an interdigitated humidity sensor that uses ZnO nanorods to make a moisture-sensitive film and found that the thinning of the spacing of ZnO nanorods will reduce the density of the moisture-sensitive film and speed up the response speed [
31]. In 2017, Zhou studied a sandwich capacitive humidity sensor with a grid-shaped top plate. The numerical results show that the response time of the microhole upper electrode sensor is 86% better than that of the grid upper electrode sensor. The response time of the humidity sensor can be improved by reducing the hole spacing, increasing the aperture, reducing film thickness, and reasonably enlarging the porosity of the film [
32]. In 2018, Liang realized a capacitive humidity sensor with a moisture-sensitive film thickness of 1.5 μm and 3.0 μm. The experimental comparison found that the sensor’s response time with a thickness of 1.5 μm was shorter than 3.0 μm, shortened from 35 s to 5 s. Still, the sensitivity was reduced from 299.9 pF/% RH to 178.6
[
33]. In 2019, Wu proposed a two-step method to optimize the layout of traditional serpentine film heaters. The substrate area within 96.4% temperature uniformity increases by more than five times compared to the conventional serpentine heater. The optimized serpentine heater realizes a coverage ratio as high as 86.64% on a rectangular membrane [
34].
It can be seen from the above studies that most of the current work starts with heating the sensor and changing the structure of the sensor unit to improve the performance of the sensor in complex environments. However, for capacitive humidity sensors, it is contradictory to improve sensor sensitivity and dynamic response time simultaneously. According to the sensor’s moisture-sensing principle, increasing the sensor’s moisture-sensing area by changing the structure of the sensor unit and the moisture-sensing layer will significantly increase the fundamental capacitance value of the sensor and amplify the sensor’s signal to be measured. Nevertheless, the time for water molecules to travel in the moisture-sensitive material is correspondingly increased. Although the sensor’s sensitivity is improved, its dynamic response time is reduced. Likewise, by lowering the moisture-sensitive area of the sensor, the opposite effect occurs. This motivates research on optimizing the sensor structure design while considering the sensitivity and dynamic response performance. In addition, in the research work on heating the humidity-sensitive structure, although the temperature rise significantly improves the sensor performance, the heating area is small, and the heat transfer is not uniform, which affects the high-performance humidity measurement.
Given the above problems, this paper takes the interdigital capacitive humidity sensor with better dynamic performance as the research object. On the one hand, based on the conformal mapping method [
35] and Fick’s second diffusion law [
36], the sensitivity and dynamic response theoretical model of the humidity sensor is established. The key structural parameters that affect the performance of the sensor are identified, and the accuracy of the theoretical model is verified through simulation modeling. On the other hand, an objective function that takes into account the sensitivity and dynamic response performance of the sensor is constructed, and the optimized design of the sensor structure is realized by using a PSO (Particle Swarm Optimization) algorithm [
37]. At the same time, a new sensor heating-electrode structure is designed, which significantly increases the electrode heating area and heat transfer uniformity, and further improves the sensor sensitivity and dynamic response performance. The main contributions of this paper are as follows:
Construct a sensitivity and dynamic response theoretical mathematical model of the interdigital capacitive humidity sensor, which obtains the interaction relationship between the sensor performance and various structural parameters;
Optimize the key structural parameters of the sensor while taking into account the sensitivity and dynamic response performance;
Design the sensor heating electrode structure and heating strategy to significantly increase the heating area and uniformity while effectively reducing system power consumption.
The rest of this paper is organized as follows. Mathematical Modeling is presented in
Section 2.
Section 3 describes the Optimizing Analysis Methods proposed in this paper. The last two parts,
Section 4 and
Section 5, are evaluation and conclusions, respectively.
3. Optimization Analysis Method
3.1. Structural Optimization Analysis
Based on the sensitivity and dynamic response model of the humidity sensor constructed in the article, the relationship between the sensor sensitivity
S, the steady-state adjustment time
, and the structural parameters is obtained, as shown in
Figure 2 and
Figure 3. Among all the structural parameters of the sensor, the greater the interdigital length
L, width
b, the thickness of the moisture-sensing layer
, and thickness of the electrode layer
, the greater the sensitivity
S of the sensor. The smaller the interdigital distance
a, the smaller the sensor sensitivity
S. In terms of dynamic adjustment time, the change of
L will not affect the adjustment time of the sensor, but smaller
a,
, and
with a larger
b can contribute to a shorter adjustment time and system’s better dynamic performance.
To ensure that the sensor has high dynamic performance with high sensitivity, in the structural design process, the influence of structural parameters on the sensor’s sensitivity and dynamic performance must be considered simultaneously. Through theoretical modeling analysis, the smaller the interdigital spacing a, the larger the interdigital width b, the higher the sensitivity, and the better the dynamic performance. Therefore, when designing and processing the sensor, it is sufficient to ensure that a and b can obtain the optimum extreme values of their performances within the acceptable error range. For the thickness of the moisture-sensitive layer and the electrode layer , the larger the value, the higher the sensitivity and the worse the dynamic performance. The improvement of sensor sensitivity and dynamic performance are contradictory. Therefore, determining the appropriate parameter values of and to make the sensor take into account both high sensitivity and high dynamic performance has become a critical step in the design and processing of the sensor.
Sensitivity
S and conditioning time
were chosen as two key metrics for sensor performance optimization. In order to ensure that the variables of different units are in the same magnitude, the two indicator functions are normalized. The sensitivity
S and adjustment time
are normalized to the interval 0–1, and the normalization function can be expressed as:
where
and
are the normalized sensitivity and steady-state time, respectively.
The objective function taking into account the sensitivity and dynamic response performance of the interdigital capacitive humidity sensor is established:
The value range of the constructed sensor structure optimization objective function is . The larger the S, the higher the sensitivity of the sensor, and the smaller the sensor , the better the dynamic performance. The purpose of the optimization design is to obtain the most suitable sensor structure parameters and to ensure that the value of the objective function is minimized.
3.2. Heat Transfer Analysis
The amount of moisture the air can hold is proportional to the temperature. When moist air encounters a cooler surface, the temperature of the humid air drops accordingly, reducing its ability to retain moisture. If the temperature is lower than a specific limit, moisture will be precipitated on the surface, resulting in condensation and affecting the humidity sensor’s performance [
39]. By adding a heating layer to the humidity sensor to make its surface temperature higher than the ambient temperature, the problem of sensor condensation is solved. In addition, based on the sensor sensitivity and dynamic response model established in the article, there is a strong relationship between temperature and the sensor. Increasing the sensor temperature within a certain range can further improve the sensor’s sensitivity and dynamic response performance, as shown in
Figure 4.
To optimize the design of the sensor heating layer and realize the optimal sensor heating strategy, it is necessary to analyze the heat exchange between the sensor and the external environment. The humidity sensor is abstracted as a uniform cuboid with area
S and height
H. The heat transfer process of the humidity sensor can be divided into three parts: thermal conduction
, thermal convection
, and thermal radiation
[
40], as shown in
Figure 5.
is the surface temperature of the sensor, and
is the ambient temperature.
The total heat dissipation of the entire humidity sensor structure can be expressed as:
Thermal conduction is the process of transferring thermal energy from a high temperature to a low temperature. The heat conduction equation of the humidity sensor can be expressed as:
where
is the thermal conductivity,
S is the area of the heated zone, and
is the temperature gradient. The established sensor model assumes thermal conductivity is a specific value inside the sensor. The sensor heat transfer
is obtained by integrating both sides of the equation:
Thermal radiation refers to the phenomenon that an object radiates electromagnetic waves due to its temperature. The thermal radiation equation of the humidity sensor can be expressed as:
where
is the emissivity of the object. Assuming the sensor is an ideal black body with
.
is the Boltzmann constant,
.
Since the thermal radiation of the system is related to the Boltzmann constant, it is much smaller than thermal convection and thermal conduction. Therefore, in order to simplify the analysis process, the thermal radiation effect of the sensor can be ignored. Thermal convection is the process of transferring heat through a flowing medium. The thermal convection equation of the humidity sensor can be expressed as:
where
is the air convection coefficient, generally
.
Through analysis, the total heat dissipation of the humidity sensor structure is:
It can be seen that when the ambient temperature and heating temperature are constant, the total heat transfer of the sensor is closely related to the structure of the heating layer. By designing the installation of the heater and optimizing the heating strategy, the system can be more suitable for extreme weather detection. While ensuring the elimination of adverse effects such as condensation, the electrodes of the humidity sensor are heated more evenly, the heating area is more expansive, and the power consumption is lower to achieve the optimal heating effect.
4. Evaluation
4.1. Accuracy Evaluation of Theoretical Model
Using the established interdigital capacitive humidity sensor sensitivity and dynamic response theoretical model, the corresponding performance index can be quickly given according to the sensor structure. To better support the optimal design of the sensor structure, the COMSOL finite element analysis software is used to establish a humidity sensor model, calculate various performance indicators of the sensor, and compare it with the theoretical model to verify the accuracy of the model, as shown in
Figure 6.
In the evaluation of the steady-state sensitivity model, the initial parameters are set according to the common humidity sensor structure. The center values of each structure are , , , and . The maximum allowable error variation coefficient is 0.4, so the variation range of each geometric structure parameter is , , , and . We simulate the low temperature and low humidity environment where the sensor is located. The steady-state analysis sets the model ambient temperature as , ambient humidity as . The geometric parameters of the sensor are changed by parametric scanning, and the sensor’s performance under different geometric parameters is obtained.
Figure 7 shows the comparison between the sensitivity calculated by the theoretical model and the results obtained by the finite element analysis at the same scale. The comparison results are verified by using the Pearson correlation coefficient. The correlation coefficient of the two methods in the thickness of the electrode layer
reaches 1, indicating that the established static sensitivity model has the same accuracy as the finite element simulation analysis for this structure. The correlation coefficient of
a,
b,
L, and
are also above 0.95, which can also meet the requirements of the optimization analysis.
In the evaluation of dynamic response models. Although the extension of the interdigital capacitive humidity sensor in the direction of length L can significantly improve the sensitivity of the sensor, it has little effect on the dynamic performance of the sensor. To improve the speed of finite element analysis and calculation and save computing resources, the model of the capacitive humidity sensor is simplified to a two-dimensional model on the y-z-plane. We set the effective diffusion coefficient of water molecules of the moisture-sensitive material as , the relative humidity of the environment changes from 10% to 20% instantaneously at the beginning of the simulation, and the steady-state error of the analysis adjustment time is less than 1%.
Figure 8 shows the comparison between the response time calculated by the theoretical model and the results obtained by the finite element analysis at the same scale. Because the dynamic analysis process of the finite element analysis method is very complex, it is necessary to couple various physical fields and set multi-layer boundary conditions. Some subtle changes will affect the output results. At the same time, the accuracy of meshing will also affect the fluctuation of the dynamic response curve. In this case, the correlation coefficient between theory and simulation for each structure is higher than 0.98.
It can be seen that the static sensitivity and dynamic response theoretical model of the interdigital capacitive humidity sensor established in this paper still has high accuracy in the low temperature and low humidity environment, which can well guide the optimal design of the sensor structure. In addition, establishing the theoretical model effectively solves the problem that traditional finite element simulation analysis takes a long time and occupies a lot of resources. It can flexibly adjust the parameters according to the actual use requirements and quickly obtain optimization results.
4.2. Optimal Design of Sensor Structure
Based on the sensor optimization objective equation established in this paper, the sensor structure was analyzed and designed by using Particle Swarm Optimization (PSO). We set the input sensor interdigital logarithm to
; the variation range of each geometric structure parameter is the same as the accuracy analysis, the ambient temperature to
, and the ambient humidity as
. The optimization problem of the sensor electrode structure can be formulated as:
According to the monotonicity analysis between the sensitivity and adjustment time of the sensor and various structural parameters, it can be determined that when a takes the minimum value, and b takes the maximum value, the utmost sensitivity and the best dynamic performance of the sensor can be obtained. For and , it is contradictory to improve the sensitivity and dynamic performance, and the sensor structure optimization scheme needs to be designed for different situations.
When the structure of the electrode layer of the humidity sensor has been processed, the moisture-sensitive layer has not been uniformly coated and suspended. That is,
has been determined, and
is still unknown. Based on the established sensor optimization objective function, the PSO optimization algorithm is used to realize the optimal design of the thickness
of the moisture-sensitive layer.
Figure 9 shows the iterative optimization process of the algorithm. After 2360 iterations, the optimal objective function
f reaches the lowest value of
, and the optimal parameter result of
is 2.88 μm. The sensor’s sensitivity is
, and the dynamic response time is 2.84 s, which ensures the optimal sensitivity and dynamic response performance under this structure.
When the humidity sensor polyimide thickness has been determined, the electrode structure has not been processed. That is,
has been determined, and
is still unknown.
Figure 10 shows the iterative optimization process of plate thickness
. After 1880 iterations, the optimal objective function
f reaches the lowest value of −0.24, and the optimal parameter result of
is 0.93 μm. At this time, the sensor’s sensitivity is
, and the dynamic response time is 3.01 s.
Through the above discussion and analysis, we have completed the structural optimization design that considers the humidity sensor’s sensitivity and dynamic response performance for different situations. To verify the accuracy and reliability of the scheme, we simulated and calculated the optimal objective function values under other sensor structures, as shown in
Figure 11. It can be seen from the figure that the theoretical calculation result of PSO is the same as that of the simulation analysis and has high accuracy. Under this structural design, the objective function of the sensor has reached the minimum value, and the performance is optimal. Similarly, according to different scenarios and size requirements, the model established in this paper can flexibly change parameters to obtain an optimal sensor design scheme that meets actual needs.
4.3. Optimal Design of Sensor Heating Layer
In the theoretical analysis of the heat transfer of the humidity sensor, the article assumes the uniformity of the heating of the sensor. It analyzes the importance of the heater structure. However, in practice, due to the influence of the heating structure, there are differences in the heating temperature of the sensor, and the difference between high and low temperatures significantly impacts the sensor’s performance. In addition, the energy for high-altitude meteorological detection is very precious, limiting the electrodes’ heating power. To expand the heating area of the sensor when the input voltage is constant, we improve the uniformity of heating, and reduce the system’s power consumption simultaneously. We need to optimize the design of the sensor heating layer structure and its heating strategy to improve sensor performance further.
To design a sensor heating electrode structure that is more in line with the needs of meteorological detection, we use Comsol finite element analysis software to analyze the heat transfer of several common heater structures. We set the ambient temperature to −70 °C, the substrate material to be silicon, the heater to be a nickel-chromium alloy, the electrode material to be silver, and the air convection coefficient to be
. At the same time, we ensure that the total volume of each heater structure is the same, which is
. We assume that the thermal conductivity variation of the sensor material is negligible to simplify the simulation model. The heating power of the four structures is set to be the same, which is 0.02 W.
Figure 12 shows the results of the steady-state finite element analysis of heat transfer for each system. The temperature of the four structures is raised to about −8 °C, much higher than the ambient temperature, which meets the conditions for eliminating sensor condensation. In addition, we drew the temperature distribution curves of each structure’s longitudinal centers to further analyze each system’s heat transfer uniformity, as shown in
Figure 13.
The total area where the surface temperature difference is less than 0.5 °C when the heating reaches a steady state is selected as the optimal heating area of the sensor. The optimum heat transfer areas of the four structures are 5.41, 5.87, 6.42, and , respectively. It can be seen from the figure that compared with the first two structures, the thicker electrode structure has a shorter total circuit length, and the heating area and uniformity are not as good as the thinner heating structure. In addition, the power loss is increased under the same input voltage due to the small resistance value of the thick electrode. The latter two systems increase the complexity of the structural design to increase the heat transfer area and reduce the useless power loss. However, the third structure destroys the uniformity of the longitudinal resistance arrangement, and the temperature is concave in the central area, and the uniformity is poor. The fourth structure is longitudinally symmetrical, and the heating result is more uniform, but the heat transfer area is slightly lower.
In terms of heating strategy, the traditional approach sets the target temperature at 20 °C regardless of the sensor’s environment. In the case of limited energy in the detection system, this solution dramatically increases the power supply pressure of the system. For the humidity sensor in the atmospheric detection system, the design of the heating electrode is mainly to solve the influence of high altitude and low temperature. The condensation phenomenon is caused by the temperature of the sensor being lower than the ambient temperature; so, in the design, it is only necessary to make the temperature of the sensor slightly higher than the ambient temperature while satisfying the regular operation of the sensor. To reduce the useless power loss of the system, we set the temperature of the electrode to be only 10 °C higher than the ambient temperature each time. This real-time closed-loop heating control strategy guarantees the sensor’s performance, saving system energy. Therefore, in the case of the same system input voltage, appropriately increasing the resistance of the heating electrode can meet the heating requirements of the sensor.
Based on the above analysis, we designed a new heater structure with full consideration of each system’s heat transfer and power consumption. It not only ensures the uniformity of temperature changes in all directions but also reduces the useless power loss by reducing the cross-sectional area of the resistance and expanding the heat transfer area. Using the same initial setting conditions, we conduct a heat transfer model analysis on the newly designed heating structure, and the results are shown in
Figure 14. Through calculation, the optimal heating area of this structure is
, and the lateral and longitudinal temperature distribution curves are relatively gentle, which meets the design requirements of heat transfer uniformity and is superior to other traditional electrode structures. In terms of the heating strategy, the closed-loop control reduces the steady-state temperature required for each stage of the system so that the heating layer can reach the preset temperature faster. It not only reduces the warm-up time and power used by the sensor but also dramatically reduces the useless power loss of the system.
5. Conclusions
This paper proposes an optimal design scheme for a highly dynamic and sensitive interdigital capacitive humidity sensor. Aiming at the problem of contradictory constraints between improving the sensitivity of the sensor and the dynamic response performance in complex environments, starting from the sensitive structure of the sensor and the heating structure, the optimal design of the sensor is completed, and the two performances of the sensor are simultaneously improved. Firstly, through the analysis of the physical field of the sensor, a theoretical model of its sensitivity and dynamic response is established, and the influence of various structural parameters and performance is explored. The model is verified by finite element analysis. The fitting coefficient can reach more than 95% under the same scale compared with the simulation results. Moreover, by fully using the established theoretical model, an objective function that satisfies the two performance optimizations of the sensor at the same time is constructed. According to the different usage scenarios of the sensor, the optimal design of the sensor structure is realized. Additionally, compared with the traditional heating structure, the novel sensor heating electrode structure is designed in this paper. It not only increases the heat transfer area of the electrode but also improves its uniformity, further enhancing the sensor’s sensitivity and dynamic response performance. The optimal design scheme proposed in this paper can also be extended to sensors of other materials, structures, and operating environments. Their performance can be optimized by changing the model parameters.