The physical and qualitative analysis of fluctuations in air and vapour concentrations in a porous medium

This work presents the development and physical analysis of a sweat transport model that couples the fluctuations in air and vapour concentrations, and temperature, in a one-dimensional porous clothing assembly. The clothing is exposed to inherent time-varying conditions due to variations in the body temperature and ambient conditions. These fluctuations are governed by a coupled system of nonlinear relaxation–transport–diffusion PDEs of Petrovskii parabolic type. A condition for the well-posedness of the resulting system of equations is derived. It is shown that the energy of the diffusion part of the system is exponentially decreasing. The boundedness and stability of the system of equations is thus confirmed. The variational formulation of the system is derived, and the existence and uniqueness of a weak solution is demonstrated analytically. This system is shown to conserve positivity. The difficulty of obtaining an analytical solution due to the complexity of the problem, urges for a numerical approach. A comparison of three cases is made using the Crank–Nicolson finite difference method (FDM). Numerical experiments show the existence of singular coefficient matrices at the site of phase change. Furthermore, the steady-state profiles of temperature, air and vapour concentrations influence the attenuation of fluctuations. Numerical results verify the analytical findings of this work.


Introduction
Porous media has been a vital subject of research in recent years with applications spanning from scientific to engineering fields. One such area of interest is the movement of sweat in the form of vapour and air through textile fabrics. The behaviour of clothing systems under various environmental conditions is of utter importance. Clothing does not solely cover the human body, but also protects and provides thermal comfort to the body. It is necessary to choose the appropriate textile to suit the external environment. Various structural properties including the choice between hygroscopic and nonhygroscopic fabric plays a major role in the dissemination of sweat through the clothing system. Studies of the transport structure in fibrous porous textiles have been published from the 1930s. Over the years, many numerical models have been presented with improved modelling and more realistic features.
Henry [1,2] studied the diffusion of moisture in cotton bales for the first time. A linear coupled system of equations with respect to vapour and temperature was formulated. To keep the equations simple, some unrealistic assumptions were made. The diffusion constant was kept independent of the vapour concentration and temperature. It was not until the 1980s that significant progress was made in this field in terms of theoretical modelling and numerical analysis. The analysis of the convection-diffusion mechanism together with phase change in a porous slab was conducted by Ogniewicz & Tien [3]. Limited to a one-dimensional steady-state formulation, the humidity conditions were varied on both sides of the slab. A quasi-steady state model introduced by Motakef & El-Masri [4] was later extended by Shapiro & Motakef [5] to incorporate time dependence on heat and moisture transport. However, the application of this model was restricted to the choice of a larger advection time scale compared to that of diffusion. Further development was noticed in the works of Farnworth [6], Vafai & Tien [7] and Fan et al. [8].
A sweating guarded hot plate was used to conduct experiments on the distribution of temperature and water content in a three-layer clothing assembly comprising an inner fibrous batting sandwiched between an inner and outer thin covering fabric [9]. The experimental values obtained with a polyester batting have been used in this paper. This paper focuses on the work Huang et al. [10]. A more realistic model was presented where vapour and air motion are treated as separate components. An analysis of the time scale was conducted to understand the occurrence of convection, conduction, phase change and diffusion. Owing to a much shorter time scale compared to that observed during experiments, a quasi-steady system of equations for the vapour and air concentrations were chosen.
However, in reality air and vapour do not have such a smooth profile as shown in the numerical experiments of Huang et al. [10]. There always exists some fluctuation in the body temperature and external environment which, in turn, influences the concentration of vapour and air in the clothing. The main aim of this paper is to add an unsteadiness in terms of fluctuations to the air and vapour concentrations in the batting. Such a study provides an insight on the mechanism of relaxation, advection and diffusion with respect to air and vapour concentrations. The variables at steady state can be linearized to simplify the existing system of equations. A qualitative analysis is made on the nature of the solution to this problem. Theoretical analysis on the system of equations for the transfer of heat and moisture in porous fabrics can be found in Wang & Sun [11] and Kelly [12]. To the best of our knowledge, to date the system has not been treated with fluctuations. In this work, a qualitative analysis of the resulting coupled nonlinear equations with fluctuations is carried out. Each of the processes of relaxation, transport and diffusion, respectively, is influenced by fluctuations in both water vapour concentration (ĉ 1 ) and air concentration (ĉ 2 ). It is observed that the cross effect of the fluctuating variables on each other plays a significant role in the model. The coefficients of the matrices in the resulting relaxationadvection-diffusion equation have been found to influence the solution qualitatively, but at the same time add to the difficulty of the qualitative analysis. The existence of a positive and unique solution eliminates the need of a numerical approach to validate the new system of equations. However, the numerics are considered owing to the difficulty of finding an analytic solution. Varying situations at the inner and external environments are taken into consideration to understand the importance of the linearized profile taken for vapour and air concentration at steady state. The behaviour of the system of equations under different environmental conditions is tested using a semi-implicit finite difference scheme. Figure 1 shows the clothing assembly used in the experiments of Fan et al. [9]. Here, the batting of this assembly is studied under different internal and external conditions. Heat moves by convection in gas and by conduction in all phases. The vapour/air mixture is transferred through the porous textile media  by convection and diffusion. Finally, phase change is induced by absorption of the fibre, evaporation and condensation. The direction of the above-mentioned processes depends on the internal and external environments. These processes are governed by the equations (2.1)-(2.4), which model heat and mass transfer through porous media. They are acquired from the balance of mass, balance of energy and Fick's law. These equations proposed in Huang et al. [10] are employed as an initial scenario in this work.

Mathematical description
where C, C v and C a are the total, vapour and air concentrations (mol m −3 ), respectively, T is the temperature (K) andW is the percentage liquid content on the fibre surface. D g and D l are the gas and liquid diffusion coefficients (m 2 s −1 ), and u g and u w represent the vapour-air mixture and water velocity (m s −1 ) given by Darcy's Law, respectively. C vt and C vg are the effective volumetric heat capacity of the batting and volumetric heat capacity of the gas mixture (J (m −3 K −1 )). Γ = Γ s + Γ ce (kg (s −1 m −3 )) is the phase change due to absorption, condensation and evaporation, and λ T (J kg −1 ) is the latent heat of phase change. ρ and ρ w are the densities for fibre and water (kg m −3 ). κ (W (m −1 K −1 )) is the effective thermal conductivity, τ c is the effective tortuosity of the batting and M (kg mol −1 ) is the molecular weight of water. The relationship between the porosity with water content, , and without water content, , is given by = − (ρ/ρ w )W(1 − ). A few assumptions are made based on the nature of the problem under investigation and these are described in the next subsection.

Assumptions
As the textile fabric being used in the batting is polyester, which is non-hygroscopic, the liquid absorption by the fibre can be neglected. The resulting non-dimensional equations are as follows: where θ, c, c 1 and c 2 are the non-dimensionalized terms for temperature, total concentration, vapour and air concentration, respectively, for 0 ≤ ξ ≤ 1. ξ and τ are, respectively, the non-dimensionalized spatial and temporal components. In addition, c = c 1 + φc 2 , where the constant φ comes from the scalings used to non-dimensionalize C v and C a . The non-dimensionalization technique has been detailed in Huang et al. [10], and the values of the parameters can be found in table 1.
Note that and p s (θ ) = 0.675p s0 (exp(1.37α(θ − 1)) − 0.0538) for θ < θ c , Looking at the advection coefficient of temperature v T /c vt which is approximately of order 10 −6 whereas that of both concentrations (vapour and air) v G = 15.8, the dynamical evolution of θ is completely different from that of c 1 and c 2 . The time scale for the evolution of air and vapour concentration is not comparable to the temperature time scale. Thus, it can be deduced that θ does not vary about the steady state. The fluctuations in temperature are very small compared to those of air and vapour concentration. As a result, a stationary temperature is maintained within the polyester batting.

Fluctuations
The human body and environment are subject to changes in temperature, air movement and humidity among others throughout the day. Even in a healthy person, the body temperature oscillates according to the time of the day and depending on their exposure [13]. The core temperature fluctuates by approximately 0.6 • C and is lowest at approximately 03.00, and highest at approximately 06.00 [14]. These fluctuations in temperature, in turn, cause fluctuations in air and vapour concentrations.

Methodology
Incorporating (2.11) and (2.12) in (2.5) and (2.6) gives the type of phenomena to be studied. The resulting system of equations is a nonlinear relaxation-transport-diffusion partial differential equation. The name Petrovskii parabolic has been given to such types of PDEs [15]. This equation governs the processes involved in the evolution of air and vapour in the porous batting. To study and understand the relaxation, transport and diffusion phenomena, the resulting equation is linearized about the steady-state values. The interactions between the different phenomena and between air and vapour are thus studied in equation (3.1). The presence of the total concentration at steady state,c, in most of the coefficient terms of (3.1) indicates a strong coupling between the evolution of the fluctuations in concentration of air and vapour in the batting. Equation (3.1) indicates that the steady state plays a vital role in the evolution of c 1 andĉ 2 .
are the relaxation, transport and diffusion matrices, respectively.
To study the evolution of the fluctuated concentrations of each of the components, the following quantities are introduced:ĉ whereĉ is the fluctuation in the total concentration,ĉ 1 is the fluctuation in vapour concentration and c 2 is the fluctuation in air concentration. After some algebra and simplification, the equation for total fluctuated concentration is given by

Boundary conditions
The boundary conditions depend on the setting of the problem. A profile at steady state is accompanied by its own set of boundary conditions depending on the conditions of the surroundings. As far as the fluctuations are concerned, Neumann boundary conditions are used. Here, care is taken to prevent the accumulation of fluctuations in the system. Hence, at the inner boundary (ξ = 0), ∂ĉ 1 ∂ξ (0, τ ) = 0 and ∂ĉ 2 ∂ξ (0, τ ) = 0, (3.3) and at the outer boundary (ξ = 1), are taken.

Theoretical analysis
It can be observed that the evolution involves relaxation, transport and diffusion. According to definition 3.1, the system of equations (3.1) is well-posed provided the diffusion matrix, d, is positive definite.
represent the Petrovskii parabolic equation (3.1). The initial value problem for the system (3.5) is wellposed if for any time Boundary conditions are said to be well-posed if the solution of the partial differential equation depends continuously on the boundary data [15]. Neumann boundary conditions are known to satisfy these criteria.
Let D be the diffusion matrix in (3.1). For the problem to be well-posed, the matrix where D sym = 1 2 (D + D T ) should be positive definite. Therefore, the well-posedness of the problem heavily depends on the values ofc,c 1 ,c 2 andθ .
Proof. Firstly, the trace of (3.7) which is given by the sum of its diagonal elements is calculated.
Hence, it is important to prove that its determinant is positive to ensure the positive definiteness of (3.7) and thus the well-posedness of (3.1). The determinant of D sym is given by Upon simplifying (3.9), the following is obtained: where v G > 0, δ G > 0 and φ > 0. The only varying terms will be the concentrations and temperature at steady state.
which proves the Theorem.
By definition, a well-posed problem is one where a solution exists, the solution is unique and the solution depends continuously on the data. However, this is rarely achieved in applied mathematical problems [16]. The problem is influenced by experimental values which are only measured at a discrete number of points. Therefore, the problem cannot depend continuously on the data, and hence becomes ill-posed. A suitable choice of initial and boundary conditions can remedy the ill-posed nature of the problem.
Next, the three main phenomena involved in the evolution of air and vapour concentration, that is, relaxation, transport and diffusion, are discussed. These processes are examined using the principle of superposition. If the sum of these three functions is a solution of a linear PDE, then by the principle of superposition, each of the functions is individually also a solution of the PDE (http://www.owlnet.rice.edu/∼ceng501/Chap7.pdf (accessed 11 May 2017)).

Relaxation, transport and diffusion
The relaxation effect ofĉ 1 onĉ 1τ andĉ 2 onĉ 2τ is governed by the equations below: and ∂ĉ 2 ∂τ = r 22ĉ2 , (3.12) and the solutions areĉ where a 1 and a 2 are constants, respectively. Irrespective of the signs of the constants a 1 and a 2 , if r 11 < 0 and r 22 < 0, (3.11) and (3.12) bring about natural relaxation.ĉ 1 andĉ 2 tend towards zero, as required. When r 11 > 0 and r 22 > 0, the solutionĉ 1 can increase infinitely and thus become unbounded. However, diffusion present in the system given by (3.1) damps the increasing relaxation term, as shown later.
The cross effect ofĉ 2 on the evolution ofĉ 1 , and vice versa, is now studied. Consider the following equations: and (3.14) With r 12 < 0 and r 21 < 0, it is indicated that an increase inĉ 2 results in a decrease inĉ 1τ , whereas decreasingĉ 2 increasesĉ 1τ . A similar observation is true for (3.14). On the other hand, if r 12 > 0 and r 21 > 0, an increase or decrease on the r.h.s. of (3.13) and (3.14) directly increases or decreases, respectively, c 1 andĉ 2 . A relaxing system causes the fluctuations to die out with time. In this situation, the wearer experiences better comfort due to the presence of a steady amount of vapour and air in the batting. Increasing amplitude of fluctuations mean the presence of more air/vapour in the fabric, which results in discomfort.
The following equations: and are pure advection equations. They do not cause any amplification or decay of the solution. These equations solely transport the fluctuations towards or away from the skin, depending on the signs of the coefficients. and Equations (3.17) and (3.18) give Equation (3.19) indicates that, for t 12 t 21 > 0, diffusion is added to the advection ofĉ 1 ; otherwise antidiffusion is included. The same findings hold forĉ 2 .
The pure diffusion equations are as follows: and Positive diffusion coefficients dissipate the fluctuations whereas negative coefficients act as antidiffusion terms. It can be seen in the following equations thatĉ 1 is affected by the diffusive nature of c 2 and vice versa. The concerned equations are and Differentiating (3.22) with respect to τ , the following is obtained: Equation (3.24) can be written as The last term in equation (3.25) adds or reduces the dissipation property of the system depending on the terms d 12 and d 21 . The same holds for the influence ofĉ 1ξξ onĉ 2τ . As pointed out earlier, this dissipative nature of the system has a major role in the boundedness of the system. Positive diffusion coefficients for the problem in 0 ≤ ξ ≤ 1 can ensure that the system (3.1) remains bounded and stable. This can be seen through the energy estimate of the diffusion process. Consider the following coupled diffusion equations: The square of the L 2 -norm onĉ 1 is defined as From (3.26) and (3.28), Similarly, it is easily seen that The energy, E(τ ) is defined as follows: where λ > 0 is a coupling parameter chosen so as to obtain the most appropriate stability result [17]. The energy equation for the coupled equations (3.26) and (3.27) is thus given by Using the Poincaré inequality [17], (3.32) can be written as The following result is thus obtained: where ω = e 1 + e 2 > 0. As a result, the diffusion equations are monotone decreasing. Clearly, the role of these phenomena depend on the system (3.1) and the steady-state scenario.

Qualitative properties of the relaxation-diffusion system
Each phenomenon has a distinct role in the evolution of the fluctuations. However, the transport phenomenon is the least dominant one. This can be seen through the Peclet number [18], which is very small, as given in the following: where l is the mean diameter of the pores. The same is confirmed by the numerical experiments in §4. As a result, the transport term will not affect the solution of (3.1) significantly. This permits the omission of the advection term in the analysis conducted under this section. Consider the relaxation-diffusion equation a.e in ∂Ω, (3.37) where v = (ˆc 1 Let Ω denote the domain such that 0 ≤ ξ ≤ 1 and τ final be such that 0 < τ < τ final . Now, the variational formulation of the problem given by (3.37) is derived. The idea of a variational formulation is to express (3.37) as a second-order ODE problem amenable to an existence and uniqueness theory.
Multiplying (3.37) by a test function ν(ξ ), which does not depend on the time τ , and integrating with respect to ξ considering the boundary conditions, the following equation is obtained: As neither Ω nor ν(ξ ) depend on τ , (3.38) can be rewritten as Owing to the fact that the variables ξ and τ have very distinct roles, they can be separated considering the solution v(ξ , τ ) as a function of time that is evaluated on a functional space defined on Ω.
and v(ξ , τ ) still takes the value v(τ )(ξ ). Next, the scalar product L 2 (Ω) and bilinear a(ω, ν) are defined as The variational formulation of (3.38), for ν(ξ ) ∈ H 1 (Ω), is given as follows: In the section that follows, it is established that the solution of (3.42) exists and is unique.
It remains to show that the series , and that its sum, v(τ ), is a solution of (3.42).
Consider the partial sum of order h of the series (3.47) We have ϕ h ∈ C([0, τ final ]; L 2 (Ω)) because each α i (τ ) is continuous. For m > h, using the orthonormal property of the eigenfunction v i , the following results: provided d −λ i − r + ≥ 0 and because the eigenvalues (λ i ) have a strictly positive increasing sequence.   influence ofĉ 1 on the influence ofĉ 2 on the influence ofĉ 1 on influence ofĉ 2 on process evolution ofĉ 1 evolution ofĉ 1 the evolution ofĉ 2 the evolution ofĉ 2 relaxation

Next, for m > h,
Hence, (3.55) The fact that ϕ h (0) converges to v 0 in L 2 (Ω), the desired initial condition v(0) = v 0 can be deduced. Clearly, v(τ ) being the sum of the series (3.47) implies that it also satisfies the variational formulation (3.42) for every test function ν = v h . As a result, v(τ ) is a solution of (3.42).
The positivity of the solution is given by proposition 3.5, which can be stated as follows.    Integrating (3.59) with respect to time, As v − (0) = (v 0 ) − = 0, for r + < 0

Verification of variational formulation
The result from proposition 3.4 is applied to the original PDE (3.37). The aim is to prove that the variational formulation can be solved for (3.37). This is done in proposition 3.8 below. Theorem 3.6 and proposition 3.7, respectively, which are quoted below, will be useful in order to establish the proof for proposition 3.8.
If the solution v is sufficiently regular, which is true by proposition 3.7, then by integration by parts the variational formulation (3.42) is equivalent to In the following section, the properties of (3.1) are verified numerically under three distinct situations, namely: -Case 1. The environmental temperature and relative humidity are considered to be lower than that of the human body. -Case 2. The relative humidity is taken to be higher in the environment compared to that of the body, while the same temperature profile as in Case 1 is maintained. -Case 3. A constant temperature and relative humidity profile is taken from the internal to the external environments.

Numerical discussion and results
The numerical solution of (3.1) is sought through the semi-implicit Crank-Nicolson finite difference scheme. This numerical scheme gives second-order accuracy both in time and space, and is unconditionally stable. The results presented in this section have been generated in MATLAB. Here, the phenomena discussed under three distinct situations are illustrated. The non-dimensional τ is taken such that 0 ≤ τ ≤ 1. A Gaussian profilê is chosen as an initial configuration. This choice is simply made to confirm the behaviour of the fluctuations with respect to the different processes discussed in §3.3 and to compare the processes under three distinct cases. The initial condition is given by the solid black line in this section.

Case 1
The steady-state profiles by Huang et al. [10] are interpolated by a least squares approximation and are given in figure 2. In this interpolation, 11 points are used from ξ = 0 to ξ = 1, with ξ = 0.1. The inner layer is exposed to the human body in contact with sweat in the form of vapour at approximately 303 K and higher humidity. The outer layer is in contact with cold moving air at approximately 253 K and lower humidity. In addition, there is no air supply at the inner cover. The equations are thus given bȳ  The behaviour of each process under these conditions is summarized as follows:

Relaxation process
From table 2, r 11 shifts from a positive to a negative value somewhere in the range 0.4 < ξ < 0.5. An overshoot is also seen in the condition number, and the relaxation matrix is singular at that location as shown in figure 3. This transition can be considered as a site of exchange between air and vapour. At 0.4 < ξ < 0.5, r 11 shifts from a positive to a negative value, that is, vapour concentration starts relaxing. On the other hand, as r 12 > 0 for ξ ≥ 0.2, it reduces the relaxing capacity ofĉ 1 . The cross effect ofĉ 1 andĉ 2 on each other is illustrated in figure 6. A decreasingĉ 2 results in an increasingĉ 1 and that is due to the values of r 12 and r 21 , as explained by (3.13) and (3.14).

Transport process
The large condition number of the transport matrix can pose a problem while inverting the matrix. However, the determinant suggests that the matrix is invertible at every point in the domain 0 ≤ ξ ≤ 1. The condition number and determinant of this transport matrix can be seen in Figure 7. t 12 has a considerable effect on the advection of vapour. This is true because air motion provides extra resistance to vapour as explained by Huang et al. [10]. However, |t 21 | |t 22 | does not contribute

Diffusion process
As d 11 ≈ d 12 ,ĉ 2 has an almost equal contribution in the diffusion of vapour. The bigger positive values of d 22 mean that air damps at a faster rate compared to vapour, as seen in figures 10 and 11. Some over damping at 0.4 < ξ < 0.6 in figure 11 is a result of extra diffusion added by (3.25).

Case 2
A different scenario is taken, where the outer clothing layer is exposed to a higher level of humidity compared to the inner layer. As a result, an inverted situation is taken for the vapour concentration at   Figure 12 shows the profiles at steady state for this case. Such a situation can exist in reality during a cold rainy day. Table 3 summarizes the signs of the coefficient matrices for Case 2.

Relaxation process
Referring to (3.11) and (3.12), r 11 < 0 and r 22 < 0 mean that both air and vapour concentrations undergo natural relaxation. Nevertheless, r 22 < r 11 < 0 results in air concentration relaxing faster than that of vapour. Over relaxation may also occur due to the large relaxation coefficients ofĉ 2 . The positive values of r 12 in the range 0.1 ≤ ξ ≤ 0.7 reduce the relaxation ofĉ 1 . The above is observed in figures 13 and 14. Clearly,ĉ 1 does not relax to zero in this case as seen in figure 13. This is in fact logical because there will always exist vapour in the batting. The high level of humidity outside and the production of sweat by the body maintains a certain concentration of vapour in the batting, and hence fluctuations will always exist.
In this case, the relaxation matrix is invertible at every point in the batting. However, at around ξ = 0.7 where the transition in r 12 takes place, a peak is observed in the condition number and the determinant reaches a minimum point as shown in figure 15.

Transport process
The transport process in Case 2 is very distinct from that of Case 1. Almost all elements are positive, which implies advection in the opposite direction. Figures 16 and 17 show the evolution ofĉ 1 andĉ 2 within the batting. Here, the fluctuationsĉ 1 andĉ 2 travel to the l.h.s., that is, towards the body. Table 3 points out the negative values of t 11 at ξ ≤ 0.2. Thus,ĉ 1 has a tendency to move away from the body.  However, as t 12 > 0,ĉ 2 has a dominant influence on the advection ofĉ 1 which prevents this movement. As the elements are almost of the same order, air and vapour fluctuations move at almost the same speed.
In this case, the transport process hints at an exchange site between air and vapour, which is shown in figure 18. In the range 0.2 < ξ < 0.3, the determinant becomes zero, and this is where t 11 shifts from negative to positive. In the previous case, this phenomenon was noted in the relaxation process.

Diffusion process
The diffusion process is the least affected process by the change made at steady state. The trend of the elements is similar to those in Case 1, d 21   This case is also dissipative in nature. The diffusion and relaxation offer enough damping to avoid the growth of the fluctuations over time.

Case 3
The final situation considered is when the vapour and air concentrations, and temperature are constant at steady state. That is, the internal and external environments are exposed to the same conditions. The profiles of such a situation is given in figure 21. All processes are constant throughout the batting. The coefficient matrices are given in table 4.

Relaxation process
The only non-zero element is r 11 . This is caused by the last term in the first element of the relaxation coefficient, in (3.1). The value of −(β G / )A √θ givesĉ 1 a slight relaxing power. In spite of this, very little relaxation is observed in the fluctuations of vapour, as seen in figure 22. Evidently, no change happens toĉ 2 .

Transport process
In this case, the fluctuations do not travel to either side of the batting. The zero coefficient elements suggest thatĉ 1 andĉ 2 remain stationary.

Diffusion process
The only visible process in this situation is diffusion. A uniform damping takes place.ĉ 2 has an equal contribution in the diffusion ofĉ 1τ because d 11 ≈ d 12 . As observed in the earlier cases, air damps more rapidly than vapour. This is seen in figures 23 and 24.
Next, the combined behaviour of all three processes in each of the three different cases is investigated.

Relaxation-transport-diffusion process
The three physical processes in some way compensate for the numerical error produced by the finite difference method, as no oscillation is noticed. The numerical results with respect to to vapour concentration are obtained at τ = 0.  Nevertheless, lim τ →0ĉ1 = 0 in Case 2. Figures 34-36 show the fluctuations damping with time, but a small amount of fluctuation will always exist in the batting.
The evolution ofĉ 2 is similar in Case 1 and 2. Air concentration diffuses and relaxes very quickly while being transported out of the polyester batting.
As observed in Cases 1-3, the steady state has an important influence on the relaxation, transport and diffusion processes. The slightest change in profile of any one parameter at steady state influences the evolution ofĉ 1 andĉ 2 significantly. The relaxation, transport and diffusion processes provide useful insight on the choice of fabric with respect to the environmental conditions. In this study, the three-layered clothing assembly is ideal for the situation in Case 1. The 10-layer polyester batting damps fluctuations in both air and vapour concentrations completely, hence giving the wearer a uniform feel. The clothing assembly also constantly eliminates wetness in the fabric, as seen in figures 31-33, providing better comfort. In Case 2, the relaxation and transport phenomena do not permit complete evacuation of sweat. In places of high humidity, a fabric with higher vapour resistance at the cover layers should be chosen. An appropriate clothing assembly should ensure a negative transport matrix in (3.1). A three-layered clothing assembly is not appropriate in Case 3. With a constant temperature of 303K between the body and environment, a thinner and hygroscopic garment will be preferable.

Conclusion
In this work, the model given in Huang et al. [ concentrations at steady state are investigated. A complex nonlinear system of Petrovskii parabolic PDE resulted from this formulation. To analyse the resulting PDE, the profiles at steady state were linearized. The importance of these profiles was depicted in three cases where it was clearly noticed that a small change made at steady state influenced the three phenomena significantly. It is shown that the system has a positive and unique solution, and would always remain bounded due to its dissipative nature. The semi-implicit Crank-Nicolson was chosen to solve the system of PDEs numerically. Numerical investigation supported the arguments made for the three distinct cases, highlighting the importance of the steady-state values in the evolution of the fluctuation of air and vapour concentration, respectively. Ideally, a clothing assembly should filter out the fluctuations and advect the vapour away from the skin. The fabric properties, which are incorporated in the values at steady state, determine its effectiveness when exposed to various conditions. It is seen that the clothing assembly with negative relaxation, negative transport and positive diffusion matrices is more capable of eliminating the fluctuations.
The fluctuations that are inherent to both the human and environmental conditions were the main focus of this study, hence the one-dimensional setting used. However, in reality the fabric in a clothing assembly may not be distributed uniformly and hence the heat and moisture transfer may not be isotropic in the clothing. As a result, a three-dimensional model is more realistic. This limitation will be addressed in future studies. In addition, the model is currently being extended in order to accommodate actual fluctuations in the ambient temperature.
Data accessibility. All the experimental parameters used in this paper have been cited from Fan et al. [9] and Huang et al. [10]. The data available in these two cited articles are sufficient for a reader to replicate our work.