# Thermal diffusivity identification of distributed parameter systems to sea ice.

1. IntroductionIn recent years, optimal control has been regarded as a powerful method with the ever-increasing demands in science, technology, and engineering [1, 2]. Many examples have an infinite number of degrees of freedom and are often modeled as infinite-dimensional distributed parameter systems (DPSs) such as thermal processes [3], spatially distributed chemical reactions [4], fluid processes [5], and flexible beams or plates [6]. Their states can vary both temporally and spatially. In the view of the complexity of these practical problems, distributed parameter systems, usually described by partial differential equations, have received a special attention and are widely used for modeling in a number of applications, relative to lumped parameter systems.

Consequently, considerable literature on the control theory of distributed parameter systems has been reported since the 1960s [7-9], for example, the solution, optimality conditions, and the analysis of an optimal control problem of DPSs. However, the applications are still limited, because no general system has all specific characteristics of every process. In many cases, it is difficult to obtain an exact partial differential equation of the process only from the physical and chemical laws, and thus, the system identification is often employed to estimate the unknown DPSs from data [10]. In this paper, the method of system identification is employed to the heat transfer problem of sea ice. Well known thermodynamic sea ice processes, especially the freeze and melt cycle in the polar regions, play an important role in the global climate system via their influence on the horizontal extent and the thickness of sea ice. The thermodynamics of sea ice, controlled by the parabolic heat conduction equation, is a typical example of distributed parameter systems.

In fact, the optimal control has been used to deal with ice problems [11-15]. Based on the sea ice temperature distribution, Lv et al. constructed an optimal control model with the thermodynamic parameters (density, specific heat, and thermal conductivity) of snow, ice, and ocean mixed layer as the control variable [11]; Fang et al. obtained the characteristics of the ice thicknesses and the ice temperatures by establishing a distributed parameter coupled system with the ice thickness as a movable inner boundary of the space domain [12]; Tan et al. presented an optimal model with the thicknesses of snow and sea ice as parametric variables [13]. They all emphasized theoretical research on sea ice physical parameters, referring to the thermodynamic coupled system, the existence of the optimal control, and the necessary conditions for optimality, which have provided an academic foundation for the simulation and forecast of sea ice thermodynamics. Furthermore, Bai et al. [14, 15] derived ice-water drag coefficients of an isolated ice floe and the thermal diffusivity of river ice with the method of parameter identification. Their work has realized an application of the optimal control of distributed parameter systems to ice parameters.

The physical parameters are crucial to exactly understand the sea ice thermal behavior. Much previous research in sea ice thermodynamics focused on the thermal conductivity, the specific heat, the latent heat, and so on, through the three determining factors--ice temperature, ice salinity, and ice density [16, 17]. Sea ice porosity which is defined as the ratio of brine and gas volume to the total volume, is also controlled by sea ice temperature, salinity, and density together [18, 19]. With the deepening of research, sea ice mechanical properties became evaluated by using ice porosity, replacing ice temperature, salinity, and density [20-22], but the thermal properties of sea ice are seldom directly studied on the basis of ice porosity. Thus, we intend to examine the thermal diffusivity of sea ice through ice porosity to fully understand the thermal process in sea ice by the method of optimal control of DPSs.

In this paper, we proceed from the one-dimensional heat conduction equation description of sea ice. Based on the field measurements from Antarctica in 2005 and 2006, an optimal control model of distributed parameter systems was put forward for identifying the thermal diffusivity of sea ice by taking thermal characteristics of sea ice and the related research results as specific constraints, and the deviation between the calculated ice temperature and the measurements was used as the performance criterion. The goal is to obtain a general parameterization of the thermal diffusivity and ice porosity for Antarctic sea ice. This identified parameterization of the thermal diffusivity of sea ice and ice porosity shows more direct relation between the thermal diffusivity and ice porosity. Meanwhile, the parameterization for the two parameters not only makes their physical meaning more clear, but also simplifies the complicated calculation by using ice thermal conductivity, specific heat capacity, and ice density.

2. Problem Analysis

Physically, the heat transfer process within ice can be described by the classical one-dimensional heat conduction equation [23]. Usually, a heat source term is added in the heat conduction equation as considering the incoming solar radiation. However, the heat source term could be ignored here, because we only discuss heat transfer in the middle ice layer, where the influence from the solar radiation is limited. The solar radiation is absent especially in polar winter. The distributed parameter systems with the initial condition and the boundary conditions can be given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where [lambda](v) is the thermal diffusivity of sea ice, v is ice porosity, T is sea ice temperature, x is the spatial variable with x [member of] [[x.sub.1], [x.sub.2]], [x.sub.1] and [x.sub.2] are the upper and lower boundaries, t is the time variable with t [member of] (0, [t.sub.f]], and [t.sub.f] is the final time for measurement. Let [OMEGA] = [[x.sub.1], [x.sub.2]], I = (0, [t.sub.f]], Q = {(x, t) [member of] [R.sup.2] :x [member of] [[x.sub.1], [x.sub.2]], and t [member of] (0, [t.sub.f]]} = [OMEGA] x I. The initial and the boundary conditions are all given continuous functions through interpolating and fitting the measured ice temperature. Because [lambda](v), [T.sub.0](x), [T.sub.1](t), and [T.sub.2](t) are all continuous and differentiable, satisfying [T.sub.0]([x.sub.1]) = [T.sub.1](0), [T.sub.0]([x.sub.2]) = [T.sub.2](0), the system (1) exists a unique solution.

The thermal diffusivity of ice is defined as the ratio of thermal conductivity to the product of ice density and specific heat. In a sense, thermal diffusivity is the measure of thermal inertia [24]. It is the most directly observed thermal property, associated with the rate of ice temperature changes. The higher the thermal diffusivity, the faster the propagation of heat into the medium. In other word, substances with high thermal diffusivity rapidly adjust their temperature to that of their surroundings, because they conduct heat quickly in comparison to their volumetric heat capacity or thermal bulk. In the previous system, the thermal diffusivity of sea ice [lambda](v) is an unknown function of ice porosity to be estimated by the system identification [10]. Therefore, it requires an optimal control model to determine the mathematical relationship [lambda](v) between the thermal diffusivity and ice porosity, which not only describes the thermal behavior, but also represents the thermal properties of sea ice. According to the previous research results in sea ice thermodynamics and the limitation of the available field measurements, the following issues should be discussed.

The formulae for the effective specific heat capacity and thermal conductivity of sea ice by Schwerdtfeger [17] are employed to calculate the thermal diffusivity of sea ice by the definition of thermal diffusivity and the in situ available data of measured ice temperature, salinity, and density. Figure 1 shows how the thermal diffusivity of sea ice varies with ice temperature, salinity, and density. According to the calculation, the minimum and the maximum of the thermal diffusivity for the study area in Antarctica are 0.58 x [10.sup.-7] [m.sup.2] x [s.sup.-1] and 10.48 x [10.sup.-7] [m.sup.2] x [s.sup.-1]. And also the statistical relationship of the porosity by Cox and Weeks [25] is used to estimate the minimum and the maximum of ice porosity, which are 11%oand 250[per thousand], respectively. In order to derive the representative results, we restrict the variation of the thermal diffusivity and ice porosity by expanding 10% of their minimum and the maximum calculated through the measured data. Consequently, the lower limits and the upper limits of the thermal diffusivity and the corresponding ice porosity in our model could be represented as

[lambda] [member of] [[[lambda].sub.1], [[lambda].sub.u]] = [0.52 x [10.sup.-7], 11.53 x [10.sup.-7]],

v [member of] [[v.sub.l], [v.sub.u]] = [10, 275], (2)

where the units of the thermal diffusivity [lambda] and ice porosity v are [m.sup.2] x [s.sup.-1] and [per thousand], respectively. The subscripts l and u denote the lower limit and the upper limit.

It is important to note that the thermal diffusivity of sea ice decreases with the increase of ice porosity, at least within the range of ice porosity from 0[per thousand] to 275[per thousand]. Mathematically, the thermal diffusivity of sea ice is a monotone function of ice porosity within its range of 0[per thousand] to 275[per thousand], which can be expressed as

d[lambda](v)/dv < 0, v [member of] (0, 275]. (3)

Another consideration is the extreme case. As ice porosity approaches 0, sea ice becomes to be pure ice, and the thermal diffusivity should be the value of pure ice. This case should not be neglected, because it originates from the concept of ice porosity. This property of sea ice is represented as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)

where the units of the thermal diffusivity [lambda] and ice porosity v are [m.sup.2] o [s.sup.-1] and [per thousand], respectively.

For the form of the parameterization of thermal diffusivity of sea ice and ice porosity is difficult to decide, formulations from literatures which are used to estimate the relations between sea ice properties and ice porosity (Table 1) are tested and modified. With these formulations, the optimal control problem for examining the parameterization of the thermal diffusivity varied with ice porosity is converted into the parameter identification of determining parameter u = (k, a, b). Let

[U.sub.ad] := {u = (k, a,b)| k [member of] [I.sub.7] = {1, 2,...,7} a [member of] [[a.sub.l], [a.sub.u]], b [member of] [[b.sub.l], [b.sub.u]]}. (5)

For all u [member of] [U.sub.ad], in the system (1), [lambda](v) could be expressed as [lambda](v; u). Parameter k [member of] [I.sub.7] = {1, 2,..., 7} denotes the number of seven different forms of the parameterizations in Table 1, which illustrates the trend of the thermal diffusivity with ice porosity. According to the discussion of the extreme previous case, parameter a should be 10.8 x [10.sup.-7] [m.sup.2] x [s.sup.-1], which represents the thermal diffusivity of pure ice, while ice porosity equals 0. And parameter b [member of] [[b.sub.l], [b.sub.u]] decides the rate of change in the thermal diffusivity curve. The related background of sea ice thermodynamics is presented in [30].

3. Thermal Diffusivity Identification of Sea Ice

Optimal control is commonly used for uncertain and complex inverse problems. On the basis of the previous comprehensive analysis, and field measurements in Antarctica, including ice temperature, ice salinity, and ice density, we could propose a succinct optimal control model to determine the parameterization of the thermal diffusivity of sea ice [31]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (6)

where T(x, t; u), is the calculated ice temperature from (1) and [bar.T](x, L) is a fitted function of the measured ice temperature. For the system (1), [lambda](v), [T.sub.0](x), and [T.sub.1](t), [T.sub.2](t) are all continuous and differentiable and satisfy [T.sub.0]([x.sub.1]) = [T.sub.1](0), [T.sub.0]([x.sub.2]) = [T.sub.2](0); thus, the system (1) exists a unique solution. Let T(x, t; u) denotes the solution of the system (1), then for all u [member of] [U.sub.ad], T(x, t; u) is continuous concerning we [U.sub.ad]. Let

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)

Let V = [[a.sub.l], [a.sub.u]] x [[b.sub.l], [b.sub.u]] [subset] [R.sup.2], then V is a compact set. According to the system (1), the mapping [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is continuous; thus, we can get that [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a compact set. Obviously, f(u) is continuous, and the set determined by the constraints in model (6) is compact. We derive [there exist][u.sup.*] [member of] [I.sub.7] x V, such that for all u [member of] [I.sub.7] x V, f([u.sup.*]) [less than or equal to] f(u). In other words, there exists an optimal parameter for the identification problem. With respect to properties and the necessary conditions for optimality and so forth, we will not repeat the detailed course of the proofs because of the limited length, and similar theoretical research has been reported in the literature [11-13].

After overall discussion of the constraints in the model, for any given k [member of] [I.sub.7], parameter a = 10.8 x [10.sup.-7] [m.sup.2] x [s.sup.-1], parameter b could be estimated by the range and the monotonicity of the thermal diffusivity.

In the process of calculation, a simple numerical method--semi-implicit difference scheme--is applied to discretize the parabolic differential equation, for the convenience of calculation, and its good performance. This method is explicit and unconditionally stable. Let [DELTA]t and [DELTA]x be time step and space step, respectively, [n.sub.N] = [t.sub.f]/[DELTA]L and [i.sub.N] = ([x.sub.2] - [x.sub.1])/[DELTA]x denote time mesh points number and space mesh points number. We can get

[partial derivative]T(x, t)/[partial derivative]t[|.sup.n + (1/2).sub.i] = [lambda](v; u)[[partial derivative].sup.]T(x, t)/[partial derivative][x.sup.2][|.sup.n + (1/2).sub.i]. (8)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (9a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (9b)

Therefore, the solution to the system (1) is simplified by

[T.sup.n + 1.sub.i] = [lambda](v; u)[DELTA]t([T.sup.n.sub.i + 1] + [T.sup.n + 1.sub.i - 1]) + ([DELTA][x.sup.2] - [lambda](v; u)[DELTA]t)[T.sup.n.sub.i]/[DELTA][x.sup.2] + [lambda](v; u)[DELTA]t i = 1, 2,..., [i.sub.N] - 1, n = 1,3,5,..., [N.sub.N], (10a)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (10b)

Considering the reality of the problem being solved, the simple genetic algorithm (GA) is employed in this paper. Genetic algorithms belong to the evolutionary algorithms, which generate solutions to optimization problems using techniques inspired by natural evolution and serve as an efficient option in science and engineering for solving practical problems [32]. Identifying the optimal settings for crossover and mutation probabilities is very important, because the variations of crossover and mutation operators influence the convergence and performance of the algorithm significantly [33]. Typical crossover probability and mutation probabilities are in the ranges of 0.5-1.0 and 0.005-0.05, respectively [34]. Therefore, we tried to test the different combinations of crossover probability and mutation probabilities in the typical ranges to determine the optimal crossover and mutation probabilities for seven formulations in Table 1 for the sake of improving the convergence performances of GA. The basic procedures are simply given as follows:

Step 1. Set the generation counter g [left arrow] 0, the maximum generation G for the evolution. Randomly generate an initial population P(0) = {[p.sub.i](u), (i = 1,2,..., N)} including N individuals [p.sub.i](u), (i = 1,2,..., N).

Step 2. Calculate the thermal diffusivity [lambda](v; u) for every individual [p.sub.i](u) in P(k) = {[p.sub.i](u), (i = 1,2,..., N)}, respectively, and judge their values. If [lambda](v; u) [member of] [[[lambda].sub.l], [[lambda].sub.u]] is satisfied, go to Step 4. Or remove the individual [p.sub.i](u) and go to Step 3 to regenerate a new one.

Step 3. Generate a new individual [p.sub.i](u) instead of the old one for the generation and go to Step 2. The goal for this is to make sure the thermal diffusivity of sea ice obtained in our study is in reasonable range.

Step 4. Evaluate the fitness value f(u) of each generated individual by the objective function in the model (6). Apparently, the fitness value is the standard for judging the optimal solution.

Step 5. Produce the next generation P(k + 1) according to the prior generation based on the crossover operation with a prespecified crossover probability and the mutation operation with a prespecified mutation probability. Here, the crossover and mutation probabilities are decided through a series of repeat tests and calculations to improve the performance of the algorithm. And the probabilities are very different for the different field measurements and seven formulations in Table 1.

Step 6. If g [less than or equal to] G, let g [left arrow] and go to Step 2. Or terminate the calculation, and the individual with the minimum of the fitness value is as the optimal solution for the identification problem.

According to the previous procedures and the measurements from seven ice samples collected in Antarctica in 2006, we derived the parameterizations of thermal diffusivity and ice porosity. Because of the various sampling dates, the obtained relations have various degrees of distinction, which are showed in Figure 2. The optimal parameterization of the thermal diffusivity and ice porosity for the Antarctic sea ice is fitted and expressed as

[lambda](v) = 10.8 x [10.sup.-7][(1 + v).sup.-0.302], v [member of] (0, 275], (11)

where the unit of the thermal diffusivity of sea ice [lambda] is [m.sup.2] x [s.sup.-1], and the unit of ice porosity is [per thousand].

In Figure 2, for the seven ice samples, the difference of parameter b causes the slight differences of the identified parameterizations. The reason for this is that the ranges of ice depth selected during identification are very various, although the ice porosity has reflected the predominant factors influencing properties of sea ice: ice temperature, ice salinity, and ice density. Most often, the initial rapid freezing on the upper surface of the ice forms granular crystals, and then, subsequent growth in the lower part of ice develops columnar crystals, which could affect the heat transfer within ice. Thus, a simple analysis was done between the proportion of the length of the granular crystals to the total length [p.sub.gc] for the ice simples used for our identification, and the values of parameter b and their correlation can be seen in Figure 3. It is obvious that the parameter b is highly dependent on ice crystal types, which illustrates that ice crystal has an important influence on the thermal diffusivity of sea ice and should be regarded as a factor to evaluate the thermal diffusivity in later work. [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

4. Numerical Simulations

To ensure the effectiveness of the optimal control model of distributed parameter systems and the validity of the parameterization of the thermal diffusivity of sea ice in this paper, we simulated vertical ice temperature profiles using the general parameterization of the thermal diffusivity varied with ice porosity (11) and compared them with the measured ice temperature from Chinese Antarctic Zhongshan Station on seven different dates in different months in 2005. Figure 4(a) shows the calculated and measured vertical temperature profiles at 13:00 on different dates. The simulated ice temperature profiles seemed to agree well with the measurements. For comparison, we also applied Schwerdtfeger's formulae and the heat conduction equation to simulate the vertical temperature profiles. Figure 4(b) shows the results.

In this paper, we introduce the average error (AE) to quantify the simulated results, which is defined as:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (12)

From Figure 4, we can conclude that the trends of ice temperature deviations are approximately the same. The AE between simulated and measured ice temperature is 1.41% (Figure 4(a)) and 1.71% (Figure 4(b)), respectively. This indicates that the identified parameterization of the thermal diffusivity of sea ice (11) is applicable in sea ice thermodynamics and that the optimal control method we discussed is also effective.

5. Conclusions

Optimal control theory of distributed parameter systems plays an important role in physics, mechanics, engineering, chemical industry, and other fields of practical activity, because of the complexity of the system processes. An optimal control model of distributed parameter systems was proposed in this paper. Instead of words, mathematical expressions were taken as specific constraints of the model by analyzing thermal properties of sea ice and summarizing the previous research results, which linked mathematics with the essence of practical problems. This model can not only depict the heat transfer process within sea ice, but can also fully capture the specific features of sea ice.

The general parameterization of the thermal diffusivity of sea ice in Antarctica was identified by applying the optimal control method: [lambda](v) = 10.8 x [10.sup.-7] [(1 + v).sup.-0,302], v [member of] (0,275]. Using the parameterization, the thermal diffusivity of sea ice could be understood more directly through ice porosity, and simply calculated by ice porosity rather than by thermal conductivity, specific heat capacity, and ice density. Furthermore, this parameterization of the thermal diffusivity has more reasonable physical significance. And the simple analysis of the difference of the parameterizations for the different ice samples shows that ice crystal types have an effect on the thermal diffusivity of sea ice and should be discussed further.

The utility of the optimal control model has been demonstrated by the simulated results. The average errors from the identified parameterization of sea ice thermal diffusivity and Schwerdtfeger's formulae revealed the advantages and growing potentials of the optimal control combined with the background of practical problems.

However, this paper focused on the application of the optimal control method to deal with the thermal diffusivity of sea ice; further efforts should be made to expand the practical application of the method.

http://dx.doi.org/ 10.1155/2013/760378

Acknowledgments

The authors are grateful to the Chinese Arctic and Antarctic Administration, who organized Chinese Antarctic scientific expeditions and provided logistical support to ensure the in situ measurements. The authors were financially supported by different resources of funding. These resources are the National Natural Science Foundation of China (Nos. 51221961, 40930848, and 11171050), the Research Council of Norway and its Norklima programme through the project "Advancing Modelling and Observing solar Radiation of Arctic sea-ice--understanding changes and processes (AMORA)" (No. 193592), and the open fund of State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology (No. LP1217). The authors would like to appreciate Prof. Matti Lepparanta, Dr. Lei Wang, and Dr. Peng Lu for helpful comments. We also thank Dr. Yinke Dou, Dr. Ruibo Lei with their team members, who carried out the fieldwork and collected valuable data during the 21st and 22nd Chinese Antarctic scientific expeditions.

References

[1] Y. Gao, "Demyanov difference of two sets and optimality conditions of Lagrange multiplier type for constrained quasi-differentiable optimization," Journal of Optimization Theory and Applications, vol. 104, no. 2, pp. 377-394, 2000.

[2] J. Chen and M. Gerdts, "Smoothing technique of nonsmooth Newton methods for control-state constrained optimal control problems," SIAM Journal on Numerical Analysis, vol. 50, no. 4, pp. 1982-2011, 2012.

[3] K. A. Hoo and D. Zheng, "Low-order control-relevant models for a class of distributed parameter systems," Chemical Engineering Science, vol. 56, no. 23, pp. 6683-6710, 2001.

[4] N. Beniich, A. El Bouhtouri, and D. Dochain, "Input constrained adaptive tracking for a nonlinear distributed parameter tubular reactor," International Journal of Adaptive Control and Signal Processing, vol. 24, no. 4, pp. 249-260, 2010.

[5] H. Li, J. Liu, C. P. Chen, and H. Deng, "A simple model-based approach for fluid dispensing analysis and control," IEEE/ASME Transactions on Mechatronics, vol. 12, no. 4, pp. 491-503, 2007.

[6] M. A. Demetriou, "Integrated actuator-sensor placement and hybrid controller design of flexible structures under worst case spatiotemporal disturbance variations," Journal of Intelligent Material Systems and Structures, vol. 15, no.12, pp.901-921, 2004.

[7] Y. Sakawa, "Solution of an optimal control problem in a distributed-parameter system," IEEE Transactions on Automatic Control, vol. 9, pp. 420-426, 1964.

[8] G. M. Bahaa, "Optimality conditions for infinite order distributed parabolic systems with multiple time delays given in integral form," Journal of Applied Mathematics, vol. 2012, Article ID 672947, 25 pages, 2012.

[9] G. Wang, "Optimal control of parabolic differential equations with two point boundary state constraints," SIAM Journal on Control and Optimization, vol. 38, no. 5, pp. 1639-1654, 2000.

[10] H. Li and C. Qi, "Modeling of distributed parameter systems for applications-a synthesized review from time-space separation," Journal of Process Control, vol. 20, no. 8, pp. 891-901,2010.

[11] W. Lv, E. Feng, and Z. Li, "A coupled thermodynamic system of sea ice and its parameter identification," Applied Mathematical Modelling, vol. 32, no. 7, pp. 1198-1207, 2008.

[12] H. Fang, J. Wang, E. Feng, and Z. Li, "Parameter identification and application of a distributed parameter coupled system with a movable inner boundary," Computers & Mathematics with Applications, vol. 62, no. 11, pp. 4015-4020, 2011.

[13] B. Tan, P. Lu, Z. Li, E. Feng, and J. Wang, "Simulation and domain identification of sea ice thermodynamic system," Journal of Applied Mathematics, vol. 2012, Article ID 532870, 16 pages, 2012.

[14] Y. Bai, Z. Li, P. Lu, L. Zhang, and E. Feng, "Optimization and identification on current drag coefficients of isolated ice floe based on experimental data," Journal of Dalian University of Technology, vol. 47, no. 6, pp. 885-889, 2007 (Chinese).

[15] Y. Bai, Z. Li, M. Han, and P. Lu, "Approaches to revise ice thermal diffusivity from measured time series of temperature in a river," in Proceedings of the 6th International Symposium on Test and Measurement (ISTM '05), pp. 2980-2986, Dalian, China, 2005.

[16] D. J. Pringle, H. Eicken, H. J. Trodahl, and L. G. E. Backstrom, "Thermal conductivity of landfast Antarctic and Arctic sea ice," Journal of Geophysical Research C, vol. 112, no. 4, Article ID C04017, 2007.

[17] P. Schwerdtfeger, "The thermal properties of sea ice," Journal of Glaciology, vol. 4, no. 36, pp. 789-807, 1963.

[18] M. Lepparanta and T. Manninen, "The brine and gas content of sea ice with attention to low salinities and high temperatures," Internal Report, Finnish Institute of Marine Research, Helsinki, Finland, 1988.

[19] Z. Li, D. S. Sodhi, and P. Lu, "Distribution of ice engineering design criteria of Bohai," Engineering Mechanics, vol. 23, no. 6, pp. 167-172, 2006 (Chinese).

[20] Z. Li, W. Peng, T. Zhang, M. Lu, and S. Gao, "Effect of porosity on the uniaxial compressive strength of ice in Liao dong Gulf," in Proceedings of the 10th International Symposium on Okhotsk Sea, Sea Ice and People, Abstracts, pp. 358-363, Hokkaido, Japan, 1995.

[21] C. Petrich, P. J. Langhorne, and Z. F. Sun, "Modelling the interrelationships between permeability, effective porosity and total porosity in sea ice," Cold Regions Science and Technology, vol. 44, no. 2, pp. 131-144, 2006.

[22] D. M. Cole, "The microstructure of ice and its influence on mechanical properties," Engineering Fracture Mechanics, vol. 68, no. 17-18, pp. 1797-1822, 2001.

[23] A. J. Semtner, "A model for the thermodynamic growth of sea ice in numerical investigations of climate," Journal of Physical Oceanography, vol. 6, no. 3, pp. 379-389, 1976.

[24] B. K. Venkanna, Fundamentals of Heat and Mass Transfer, PHI Learning, New Delhi, ndia, 2010.

[25] G. F. N. Cox and W. F. Weeks, "Equations for determining the gas and brine volumes in sea-ice samples," Journal of Glaciology, vol. 29, no. 102, pp. 306-316, 1983.

[26] D. L. Anderson and W. F. Weeks, "A theoretical analysis of sea ice strength," Transactions, American Geophysical Union, vol. 39, no. 4, pp. 632-640, 1958.

[27] W. F. Weeks and A. Assur, The Mechanical Properties of Sea Ice, Part II, Section C3, V.S. Cold Regions Research and Engineering Laboratory, Cold Regions Science and Engineering, Hanover, Germany, 1967.

[28] H. D. Gans and P. R. Woodmansee, "An improved method of modeling porous materials using finite elements," Computers and Structures, vol. 44, no. 5, pp. 1055-1063, 1992.

[29] P. O. Moslet, "Field testing of uniaxial compression strength of columnar sea ice," Cold Regions Science and Technology, vol.48, no. 1, pp. 1-14, 2007.

[30] L. Shi, Y. Bai, Z. Li, B. Cheng, and M. Lepparanta, "Preliminary results on relationship between thermal diffusivity and porosity of sea ice in the Antarctic," Chinese Journal of Polar Science, vol. 20, no. 1, pp. 72-80, 2009.

[31] J. Nocedal and S. J. Wright, Numerical Optimization, Springer Series in Operations Research and Financial Engineering, Springer, New York, NY, USA, 1999.

[32] M. Melanie, An Introduction to Genetic Algorithms, MIT Press, Cambridge, Mass, USA, 1996.

[33] Y. J. Cao and Q. H. Wu, "Convergence analysis of adaptive genetic algorithms," in Proceedings of the 2nd International Conference on Genetic Algorithms in Engineering Systems: Innovations and Applications (GALESIA '97), No. 446, pp.85-89, September 1997.

[34] M. Srinivas and L. M. Patnaik, "Adaptive probabilities of crossover and mutation in genetic algorithms," IEEE Transactions on Systems, Man and Cybernetics, vol. 24, no. 4, pp. 656-667, 1994.

Liqiong Shi, (1) Zhijun Li, (1) Enmin Feng, (2) Yila Bai, (3) and Yu Yang (1)

(1) State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China

(2) School of Mathematical Sciences, Dalian University of Technology, Dalian 116024, China

(3) Mathematical and Physics Department, Bohai University, Jinzhou 121000, China

Correspondence should be addressed to Liqiong Shi; shiliqiong333@126.com

Received 24 January 2013; Accepted 7 July 2013

Academic Editor: J. Liang

TABLE 1: Formulations of sea ice properties and ice porosity from literatures and modified. Number Formulation Reference and origin 1 [lambda](v) = [26], used for the a[(1 - relation of sea ice b[square root tensile strength and of brine volume v]).sup.-1] 2 [lambda](v) = [27], used for the a[(1 - relation ofthe elastic bv).sup.4] modulus and the porosity 3 [lambda](v) = [27], used for the a[(1 - relation of the effect b[square root stress and the porosity of v]).sup.-2] 4 [lambda](v) = [28], used for the a x elastic modulus and the [e.sup.bv] porosity 5 [lambda](v) = [19], used for the a[(1 + relation of the v).sup.b] mechanical indicator and the porosity 6 [lambda](v) = [29], used for the a[(1 - relation of maximum b[square root intensity and the of v]).sup.2] porosity 7 [lambda](v) = Newly modified form based a[(1 - on references bv).sup.-1]

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Shi, Liqiong; Li, Zhijun; Feng, Enmin; Bai, Yila; Yang, Yu |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2013 |

Words: | 5215 |

Previous Article: | A probabilistic damage identification approach for structures under unknown excitation and with measurement uncertainties. |

Next Article: | Noise suppression in ECG signals through efficient one-step wavelet processing techniques. |

Topics: |