Abstract
This paper deals with a new inverse problem of estimating textile fabrics thickness based on a one-dimensional steady-state heat and moisture transfer model. First, we give an existence and uniqueness theorem for the solution to the heat and moisture transfer problem in the Banach space. Then we present an inverse problem for the determination of textile fabrics thickness and formulate the inverse problem as an optimization problem with maximum probability, according to the requirement of heat-moisture comfort of clothing. Moreover, considering that the objective function of the proposed optimization problem may not be continuous and the points achieving maximum value may not be unique, we proceed to recast this inverse problem as a minimum norm with a maximum probability constraint problem. Finally, we use a static penalty method to convert the constrained problem into an equivalent unconstrained minimization problem and obtain the solution for the optimization problem by a stochastic search method known as particle swarm optimization algorithm. Numerical experiments show that our new model and its equivalent transformation is quite acceptable, and the proposed numerical method is valid and robust.
The notation used in the paper is defined as follows:
Nomenclature | ||
= | distance from the inner covering fabric (m) | |
= | thickness of the fabric batting (m) | |
= | temperature (K) | |
= | temperature of the surface next to the human body (K) | |
= | temperature of the surrounding air (K) | |
= | water vapour pressure (Pa) | |
= | water vapour pressure of the surrounding air (Pa) | |
= | mass flux of water vapour | |
= | mass flux of water vapour of the surface next to the human body | |
= | radius of cylindrical pore (m) | |
= | heat transfer resistance of textile material | |
= | the convective heat transfer coefficient | |
= | constants which are related with molecular weight and gas constant | |
= | relative humidity | |
= | saturated water vapour pressure at temperature (Pa) |
Greek symbols | ||
= | latent heat of sorption and condensation of water vapour (J ) | |
= | effective thermal conductivity of the porous batting | |
= | effective tortuosity of the textile | |
= | porosity of textile surface (%) | |
= | rate of condensation |
1 Introduction
Clothing has served a variety of social purposes. The primary function of clothing is to improve the comfort of the wearer. In hot climates, clothing provides protection from sunburn or wind damage, while in cold climates its thermal insulation properties are generally more important. On this account, happiness and health of humans largely depends on clothing. The main effect of clothing will be the influence on the heat and moisture exchange between the skin and the environment. An understanding of the heat and moisture transfer in clothing system has developed over the past few decades. Henry [Citation1] developed one of the first theories of coupled heat and moisture transfer through clothing considering accumulation effects. Shitzer and Chato [Citation2] studied the heat and mass transfer of the air-clothing-skin system, where the heat and mass transfer in the system was considered as a steady-state problem in a one-dimensional model composed of five layers. Jones [Citation3] and Li and Holcombe [Citation4] introduced dynamic models of heat and moisture characteristics in interaction with two node human thermoregulation model of Gagge et al. [Citation5], respectively. Fan et al. [Citation6, Citation7], Fan and Wen [Citation8] developed more general dynamic models in textile materials. Min et al. [Citation9] set up a mathematical model to simulate the heat and moisture transfer from skin to environment through fabrics by including heat radiation between surfaces and the surface diffusion along fibres. Wu and Fan [Citation10] and Wissler and Havenith [Citation11] presented mathematical models for heat and moisture transport through multi-layer clothing, respectively.
Heat-moisture comfort refers to the condition of the mind which expresses satisfaction with the heat and moisture environment. The heat-moisture comfort indices have been developed with the help of coupled heat and moisture transfer equations designed to predict heat-moisture comfort. The heat-moisture comfort feeling of the human body is influenced by various parameters. The three main factors are: environmental climate, clothing and physical activity. Clothing factors such as thickness of material, clothing insulation value, type of material, volume of air enclosed in the fibres, resistance to diffusion of water vapour, layers of clothing, body fit, air permeability, etc influence the thermal and moisture comfort. A considerable amount of research has been done in order to ascertain the influence of various factors of clothing on comfort, the interested readers may refer to some articles [Citation12–Citation16] for further understanding.
The direct heat and moisture transfer problems consist of determining the temperature and moisture distribution within the fabrics when the boundary and initial conditions, physical properties and geometric parameters are known. In contrast, the inverse heat and moisture transfer problems consider the identification of boundary and initial conditions, physical properties and geometric parameters by given measurements. These inverse problems are beneficial for the development of advanced textile materials to protect against harmful environmental stresses including the climatic conditions.
Only few works have been done on the inverse problems in textile material design. Xu et al. [Citation17, Citation18] studied the inverse heat and moisture transfer problems based on a one-dimensional steady-state model of heat and moisture transfer through parallel cylinder pore textiles for estimation of material thickness and type. Xu and Ge [Citation19] also presented an inverse problem of textile thickness determination based on a heat and moisture transfer model with phase change in porous textile materials. Numerical methods and simulations for the heat and moisture transport in porous textile materials have been studied in their papers. However, the authors haven’t explored the existence and uniqueness analysis for the heat and moisture transfer problem, and the according inverse problems are formulated as unconstrained minimization problems by choosing the ideal relative humidity of the surface close to skin as a fixed percentage, for example 50%, and the numerical treatment for that inverse problem is mainly applied by the direct search methods: Hooke–Jeeves (H–J) algorithm and Golden-Section (G-S) method.
In this paper, based on some assumptions, we establish an existence and uniqueness result for the solution to the heat and moisture transfer problem. Considering that, the ideal relative humility in the inverse heat and moisture transfer problems can vary from 40 to 60%, we design a new maximum probability optimization model for the inverse heat and moisture transfer problem, and then we recast this inverse problem as a minimum norm with a maximum probability constraint problem. We transform the proposed constrained problem into an equivalent unconstrained problem by the static penalty method, and then adopt the particle swarm optimization algorithm to obtain the approximate solution for the new optimization problem. Three representative fabrics are selected for our model and algorithm testification. The numerical results show that our new optimization model is quite reasonable and the proposed algorithm is quite robust.
This paper is organized as follows. In Section 2, the heat and moisture transfer problem based on a steady-state model is described and solved. The Banach fixed-point theorem is applied to prove the existence and uniqueness of a solution of the direct problem considered. In Section 3, an inverse problem of thickness determination for textile fabric and its transformation into a constrained optimization problem are presented. An overview of the particle swarm optimization (PSO) algorithm is also illustrated in this part. The numerical examples including the noise perturbation for the optimization problem are discussed in Section 4. Section 5 ends the paper by making some conclusions.
2 The direct heat and moisture transfer problem
We first introduce the direct problem for the heat and moisture transfer problem. And then we give an existence and uniqueness theorem for the solution to the above problem. Finally, we illustrate the finite difference method for solving it.
2.1 The direct problem
Let us consider a system that consists of human skin, fabric and environment as shown in Figure . Considering that the fabric may be distributed uniformly and the mass and heat influx may uniformly fall on the fabric as well. Furthermore, assuming that the distribution of temperature, water vapour pressure and mass flux of water vapour proliferation within the fabric are all steady as the time changes.
Mass flux of water vapour proliferation within the parallel cylindrical pore can be represented as1 1 Combining the Equation (Equation11 1 ) with the mass conservation law, the energy conservation law and the rate of condensation equation, one obtains a one-dimensional model of the coupled heat and mass transfer problem through parallel cylinder pore textiles as follows [Citation20] :2 2 which subjected to the following boundary conditions:3 3 where the saturation vapour pressure is related to temperature by:4 4 Our model is based on the assumptions listed below.
(1) | The inner pore structure of the fabric is assumed to be parallel cylinder, and there are no voids in the pores. | ||||
(2) | The fibrous porous batting is isotropic and homogeneous in fibre arrangement and material properties. | ||||
(3) | The air and vapour are continuous and treated as ideal gases. | ||||
(4) | Volume changes of the fibres due to changing moisture and water content are neglected. |
2.2 Theoretical result for the direct problem
Let and . We decouple the ordinary differential Equation (Equation22 2 ) to obtain the following two-point boundary value problem of nonlinear integral-differential equation:5 5 where
The two-point boundary value problem (Equation55 5 ) can be rewritten aswhere .
Before we present our main theorem, some basic assumptions are given as follows:
(1) | |||||
(2) | belongs to , i.e. there exists a constant , such that | ||||
(3) | There exist constants such that | ||||
(4) | There exists a constant such that |
Theorem 3.1
Assume (1)–(4) hold and the constant satisfieswhereThen there exists a unique solution in the space such that6 6
Proof
For any , we have By the assumption (Equation11 1 ) and mean value theorem of differentials, we getwhere lies between and
In terms of the above basic assumptions (Equation22 2 ) and (Equation33 3 ) and some classic inequalities, we obtain the following estimatesandTherefore, we haveBy the Banach fixed-point theorem, there exists a unique solution in the space such thatwhich completes the proof.
Under the same assumptions as that in Theorem 3.1, we conclude that there exists a unique solution to the direct problem (Equation55 5 ).
2.3 The finite difference method to solve the direct problem
The direct problem can be solved by discretizing the above equations with the finite difference method. Let be the step length in space coordinates, and denote partitions of . The finite difference approximation for the differential equation of (Equation55 5 ) may be written in the formand7 7 where is the approximate value of .
Using and , we approximateUsing , we find8 8 and the relative humidity at the boundary next to the human body9 9
3 The inverse heat and moisture transfer problem of determining textile fabrics thickness
In this section, we introduce a new inverse heat and moisture transfer problem and the corresponding PSO algorithm.
3.1 Formulation of the inverse problem of thickness determination of textile fabrics
If the parameters such as thickness, radius of cylindrical pore, conductivity and porosity are given in the Equation (Equation22 2 ), combining with the boundary conditions (Equation33 3 ), one obtains the temperature and moisture distribution within the fabrics, which is the direct problem. But if any one or more parameters are unknown, it turns into an inverse problem. In what follows, we assume that the other parameters are known except for the textiles thickness .
The heat and moisture comfort indices can be used as an additional information to determine the unknown parameter. Huang [Citation21] showed that when the temperature of microclimate, which is located between the skin surface and the inner surface of fabric, varies from 31 to C, the relative humidity from 40 to 60% and the wind speed from 0.1 to 0.4 m s, the human body feels comfortable. By using these comfort indices, the inverse problem can be recast as an optimization problem.
3.2 Minimum norm with maximum probability constraint problem
Suppose and are minimum and maximum average temperature at a particular place during a particular period, respectively. and are minimum and maximum average relative humidity, respectively. We discretize the intervals and into a series of , cells, namelythus we obtain combinations about the environmental temperature and relative humidityIf is given, we can obtain the according relative humility of the surface close to skin calculated from direct problem, denoted by which vary from different and . Suppose be the probability of an event that where the thickness of fabric is . We want to find the optimal thickness, which has the maximum probability to make the human body feels comfortable, i.e. the maximum probability of value lies in [40%, 60%]. We obtain the following maximum optimization problem:10 10 The objective function may not be continuous and its points achieve maximum value may not be unique. Let be the maximum value of . Considering that the thinner fabric the better the human body feels. Therefore, the original optimization problem (Equation1010 10 ) can be modified as the following constrained minimization problem11 11 12 12 This objective function is linear and combined with the nonlinear constraint. Suppose be the solution to this constrained optimization problem, then we call it as a minimum norm with maximum probability constraint problem in textile fabric.
The static penalty method [Citation22] can be used to transform a constrained problem into an unconstrained problem. For this method, the constraints are placed into the objective function via a penalty parameter in a way that penalizes any violation of the constraints. Using static penalty method, the constrained problem can be replaced by the following unconstrained problem13 13 where is a large positive constant. We can intuitively see that an optimal solution to the above problem must have very close to , because otherwise a large penalty will be incurred.
Due to the existence of probability function , the objective function of the optimization problem (Equation1313 13 ) may not be continuous. Conventional optimization technique such as gradient based algorithm are not good enough to solve the problems involving the probability function because it depends on first and second derivatives of the objective function. Thus, it’s important to find an appropriate method for solving such a complicate problem. In what follows, we use a stochastic method known as particle swarm optimization algorithm to solve the above optimization problem. The advantages of the PSO algorithm is its capability in searching for the global optimum and no computation of the complicated gradients.
3.3 Particle swarm optimization
The PSO method [Citation23–Citation26] is an efficient technique for solving many complex optimization problems which are nonlinear, nondifferentiable and multi-modal. It is becoming very popular because of its simplicity of implementation as well as ability to swiftly converge to a good solution. It does not require any gradient information of the function to be optimized and uses only primitive mathematical operators.
In PSO, the set of particles is called a swarm. A swarm consists of N particles moving around in a D-dimensional search space. The position of the ith particle can be represented byThe velocity for the ith particle can be written asThe positions and velocities of the particles are confined within and , respectively. Each particle coexists and evolves simultaneously based on knowledge shared with neighbouring particles, and it makes use of its own memory and knowledge gained by the swarm as a whole to find the best solution. The best previously encountered position of the th particle is denoted its individual best positiona value called . The best value of the all individual values is denoted the global best positionand called . The PSO process is initialized with a population of random particles, and the algorithm then executes a search for optimal solutions by continuously updating generations. At each generation, the position and velocity of the ith particle are updated by and in the swarm. The update equations can be formulated as:14 14 15 15 where and are random numbers between , and and are acceleration constants, which control how far a particle will move in a single generation. Velocities and denote the velocities of the new and old particle, respectively. is the current particle position, and is the new, updated particle position. The inertia weight w controls the impact of the previous velocity of a particle on its current one.
The typical pseudo-code of the PSO process is shown below.
The algorithm only requires fitness function for each of the particle, without assumption such as continuity and differentiability, which makes it very useful for a discontinuous function.
4 Numerical examples and discussion
We select three representative fabrics: wool, polyester and polypropylene fibre for the proposed optimization problem (Equation1313 13 ). In the model, the parameters and for these three kinds of fabrics are listed in Table . The parameters used in the study are taken asThe boundary conditions next to skin areThe numbers of grid points are all taken as 30.
The parameters in PSO are set as and the penalty parameter
We chose two different combination of the temperature and humidity in the environment for the numerical calculation.
Case 1 The temperature and humidity in the environment are assumed to and
Case 2 The temperature and humidity in the environment are assumed to and
It can be seen from Figures and that thickness of textile fabric have significant implication for heat-moisture comfort, as it means, too thin or too thick of a textile fabric will make human feeling uncomfortable.
In order to compare the results involving random measurement noise, we add a uniform distribution uncorrelated errors. The simulated inexact measurement data can be expressed as where are the exact data, and is the noise level. In Figures and , we plot the objective function L with the noisy measurements and for wool, polyester and polypropylene fibre in Cases 1 and 2, respectively.
Except for the number of particles, the computational time involved in the PSO method mainly depends on the structure of proposed optimization function. To obtain a global minima for the optimization problem (Equation1313 13 ), in each run of the PSO searching, one needs to solve direct problems, each direct problem is composed by the nonlinear integral-differential Equation (Equation55 5 ) with different boundary conditions. Thus, the computational time normally needs about 1 h in our PC computer. Tables and show the estimated results with different noise levels for three kinds of fabrics by the PSO algorithm in Cases 1 and 2, respectively. It is noted that the thickness of textile fabric changes with the environment conditions.
5 Conclusions
We present a new inverse problem that can determine the unknown thickness based on a minimum norm with the maximum probability constraint optimization problem. Comparing with the existing inverse problem to determine textile fabrics thickness, the new constrained optimization model is more reasonable and practical, which provides an effective way to study the inverse problems of heat and moisture transport in porous textile materials.
We adopt PSO algorithm to solve the proposed optimization problem. The advantages of this stochastic search method is that it does not require the gradient computation and the choice of initial guess. The numerical results is in relatively good agreement with the values as shown in Figures –. We can verify that the results in Tables and are the minimum norm solutions with maximum probability by direct computation. Therefore, one can conclude that the PSO is stable to deal with this inverse problem even with small noisy measurements and the proposed new optimization is theoretically well-founded.
Further works in the future are still needed. One is to establish the stability result on the solution to the strong nonlinearity boundary value problem (Equation22 2 ) and (Equation33 3 ), and the other is to adopt our maximum probability constraint technique to other inverse heat and moisture transfer problems, especially for the inverse problem considering fabric composed of multiple layers.
Acknowledgements
The authors kindly acknowledge the help of the anonymous referees in improving the readability of the paper.
Additional information
Funding
References
- Henry P. Diffusion in absorbing media. Proc. Royal Soc. London. Ser. A Math. Phys. Sci. 1939;171:215–241.
- Shitzer A, Chato JC. Thermal interaction with garments, heat transfer in medicine and biology: analysis and applications. New York, NY: Plenum Press; 1985.
- Jones B. Transient interaction between the human and the thermal environment. ASHRAE Trans. 1992;98:189–195.
- Li Y, Holcombe B. Mathematical simulation of heat and moisture transfer in a human-clothing-environment system. Text. Res. J. 1998;68:389–397.
- Gagge A, Stolwijk J, Nishi Y. An effective temperature scale based on a simple model of human physiological regulatiry response. Vol. 13; Sapporo: Memoirs of the Faculty of Engineering, Hokkaido University; 1972. p. 21–36.
- Fan J, Cheng X, Wen X, Sun W. An improved model of heat and moisture transfer with phase change and mobile condensates in fibrous insulation and comparison with experimental results. Int. J. Heat Mass Transfer. 2004;47:2343–2352.
- Fan J, Luo Z, Li Y. Heat and moisture transfer with sorption and condensation in porous clothing assemblies and numerical simulation. Int. J. Heat Mass Transfer. 2000;43:2989–3000.
- Fan J, Wen X. Modeling heat and moisture transfer through fibrous insulation with phase change and mobile condensates. Int. J. Heat Mass Transfer. 2002;45:4045–4055.
- Min K, Son Y, Kim C, Lee Y, Hong K. Heat and moisture transfer from skin to environment through fabrics: a mathematical model. Int. J. Heat Mass Transfer. 2007;50:5292–5304.
- Wu H, Fan J. Study of heat and moisture transfer within multi-layer clothing assemblies consisting of different types of battings. Int. J. Therm. Sci. 2008;47:641–647.
- Wissler E, Havenith G. A simple theoretical model of heat and moisture transport in multi-layer garments in cool ambient air. Eur. J. Appl. Physiol. 2009;105:797–808.
- Amrit U. Bedding textiles and their influence on thermal comfort and sleep. AUTEX Res. J. 2007;8:252–254.
- Li Y. Perceptions of temperature, moisture and comfort in clothing during environmental transients. Ergonomics. 2005;48:234–248.
- Scheurell D, Spivak S, Hollies N. Dynamic surface wetness of fabrics in relation to clothing comfort. Text. Res. J. 1985;55:394–399.
- Wang Z, Li Y, Kowk Y, Yeung C. Mathematical simulation of the perception of fabric thermal and moisture sensations. Text. Res. J. 2002;72:327–334.
- Zhang P, Gong R, Yanai Y, Tokura H. Influence of clothing material properties on rectal temperature in different environments. Int. J. Cloth. Sci. Tech. 2002;14:299–306.
- Xu D, Chen Y, Zhou X. Type design for bilayer textile materials under low temperature: modeling, numerical algorithm and simulation. Int. J. Heat Mass Transfer. 2013;60:582–590.
- Xu D, Cheng J, Zhou X. An inverse problem of thickness design for single layer textile material under low temperature. J. Math Ind. 2010;2:139–146.
- Xu D, Ge M. Thickness determination in textile material design: dynamic modeling and numerical algorithms. Inverse Prob. 2012;28:35011–35032.
- Xu D, Cheng J, Zhou X. A model of heat and moisture transfer through the parallel pore textiles. J. Fiber Bioeng. Inf. 2011;4:250–255.
- Huang J. Clothing comfort. Beijing: Science Press; 2008.
- Bazaraa M, Sherali H, Shetty C. Nonlinear programming. New York, NY: Wiley Online Library;1993.
- Kennedy J, Eberhart R. Particle swarm optimization. In: 1995 Proceedings of the IEEE International Conference on Neural Networks. Vol. 4; Perth; 1995. p. 1942–1948.
- Clerc M. Particle swarm optimization. Vol. 67. London: Recherche; 2006.
- Lazinica A. Particle swarm optimization. Kirchengasse: InTech; 2009.
- Parsopoulos K, Vrahatis M. Particle swarm optimization and intelligence: advances and applications. Chicago, IL: Information Science Reference; 2010.