In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with ...In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with prescribed 2-norm has some normalized solutions by introducing variational methods.展开更多
This article studies a class of nonlinear Kirchhoff equations with exponential critical growth,trapping potential,and perturbation.Under appropriate assumptions about f and h,the article obtained the existence of norm...This article studies a class of nonlinear Kirchhoff equations with exponential critical growth,trapping potential,and perturbation.Under appropriate assumptions about f and h,the article obtained the existence of normalized positive solutions for this equation via the Trudinger-Moser inequality and variational methods.Moreover,these solutions are also ground state solutions.Additionally,the article also characterized the asymptotic behavior of solutions.The results of this article expand the research of relevant literature.展开更多
This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solutio...This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.展开更多
The solution diffusion coefficient is a great important intrinsical parameter in crystal growth.On earth,it is impossible to accurately determine the diffusion coefficient since there is nature convection.One of the m...The solution diffusion coefficient is a great important intrinsical parameter in crystal growth.On earth,it is impossible to accurately determine the diffusion coefficient since there is nature convection.One of the marked charateristics of space-crystal growth is to eleminate nature convection,so that purely diffusion-controlled condition of crystal growth could be realized and precise measurement of the diffusion coefficient should be approved.展开更多
Gradiently denitrated gun propellant(GDGP)prepared by a“gradient denitration”strategy is obviously superior in progressive burning performance to the traditional deterred gun propellant.Currently,the preparation of ...Gradiently denitrated gun propellant(GDGP)prepared by a“gradient denitration”strategy is obviously superior in progressive burning performance to the traditional deterred gun propellant.Currently,the preparation of GDGP employed a tedious two-step method involving organic solvents,which hinders the large-scale preparation of GDGP.In this paper,GDGP was successfully prepared via a novelty and environmentally friendly one-step method.The obtained samples were characterized by FT-IR,Raman,SEM and XPS.The results showed that the content of nitrate groups gradiently increased from the surface to the core in the surface layer of GDGP and the surface layer of GDGP exhibited a higher compaction than that of raw gun propellant,with a well-preserved nitrocellulose structure.The denitration process enabled the propellant surface with regressive energy density and good progressive burning performance,as confirmed by oxygen bomb and closed bomb test.At the same time,the effects of different solvents on the component loss of propellant were compared.The result showed that water caused the least component loss.Finally,the stability of GDGP was confirmed by methyl-violet test.This work not only provided environmentally friendly,simple and economic preparation of GDGP,but also confirmed the stability of GDGP prepared by this method.展开更多
Considering three longitudinal displacement functions and uniform axial displacement functions for shear lag effect and uniform axial deformation of thin-walled box girder with varying depths,a simple and efficient me...Considering three longitudinal displacement functions and uniform axial displacement functions for shear lag effect and uniform axial deformation of thin-walled box girder with varying depths,a simple and efficient method with high precision to analyze the shear lag effect of thin-walled box girders was proposed.The governing differential equations and boundary conditions of the box girder under lateral loading were derived based on the energy-variational method,and closed-form solutions to stress and deflection corresponding to lateral loading were obtained.Analysis and calculations were carried out with respect to a trapezoidal box girder under concentrated loading or uniform loading and a rectangular box girder under concentrated loading.The analytical results were compared with numerical solutions derived according to the high order finite strip element method and the experimental results.The investigation shows that the closed-form solution is in good agreement with the numerical solutions derived according to the high order finite strip method and the experimental results,and has good stability.Because of the shear lag effect,the stress in cross-section centroid is no longer zero,thus it is not reasonable enough to assume that the strain in cross-section centroid is zero without considering uniform axial deformation.展开更多
By combining the results of laboratory model tests with relevant flow rules, the failure mode of shallow unsymmetrical loading tunnels and the corresponding velocity field were established. According to the principle ...By combining the results of laboratory model tests with relevant flow rules, the failure mode of shallow unsymmetrical loading tunnels and the corresponding velocity field were established. According to the principle of virtual power, the upper bound solution for surrounding rock pressure of shallow unsymmetrical loading tunnel was derived and verified by an example. The results indicate that the calculated results of the derived upper bound method for surrounding rock pressure of shallow unsymmetrical loading tunnels are relatively close to those of the existing "code method" and test results, which means that the proposed method is feasible. The current code method underestimates the unsymmetrical loading feature of surrounding rock pressure of shallow unsymmetrical loading tunnels, so it is unsafe; when the burial depth is less or greater than two times of the tunnel span and the unsymmetrical loading angle is less than 45°, the upper bound method or the average value of the results calculated by the upper bound method and code method respectively, is comparatively reasonable. When the burial depth is greater than two times of the tunnel span and the unsymmetrical loading angle is greater than 45°, the code method is more suitable.展开更多
The optimization problem is considered in which the objective function is pseudolinear（both pseudoconvex and pseudoconcave） and the constraints are linear. The general expression for the optimal solutions to the pro...The optimization problem is considered in which the objective function is pseudolinear（both pseudoconvex and pseudoconcave） and the constraints are linear. The general expression for the optimal solutions to the problem is derived with the representation theorem of polyhedral sets, and the uniqueness condition of the optimal solution and the computational procedures to determine all optimal solutions （if the uniqueness condition is not satisfied ） are provided. Finally, an illustrative example is also given.展开更多
An exact solution for simply-supported laminated beams with material properties variable with temperature under a combination of uniform thermo-load and mechanical loads was investigated,based on the two-dimensional(2...An exact solution for simply-supported laminated beams with material properties variable with temperature under a combination of uniform thermo-load and mechanical loads was investigated,based on the two-dimensional(2-D)thermo-elasticity theory.Firstly,the beam was divided into a series of layers with uniform material properties along the interfaces of the beam.The uniform thermo-load acted on each layer was transformed into a combination of the normal surface forces acted at the two ends and the transverse thermo-load.Secondly,the state space method was employed to obtain the general solutions of displacements and stresses in an arbitrary layer.Thirdly,based on the interfacial continuity conditions between adjacent layers,the relations of displacement and stress components between the top and bottom layers of the beam were recursively derived by use of the transfer-matrix method.The unknowns in the solutions can be solved by the mechanical loads acted on the top and bottom surfaces.The convergence of the present solutions was checked.The comparative study of the present solutions with the Timoshenko’s solutions and the finite element(FE)solutions was carried out.The effects of material properties variable with temperature on the thermo-elastic behavior of laminated beams were discussed in detail.展开更多
To make heat conduction equation embody the essence of physical phenomenon under study, dimensionless factors were introduced and the transient heat conduction equation and its boundary conditions were transformed to ...To make heat conduction equation embody the essence of physical phenomenon under study, dimensionless factors were introduced and the transient heat conduction equation and its boundary conditions were transformed to dimensionless forms. Then, a theoretical solution model of transient heat conduction problem in one-dimensional double-layer composite medium was built utilizing the natural eigenfunction expansion method. In order to verify the validity of the model, the results of the above theoretical solution were compared with those of finite element method. The results by the two methods are in a good agreement. The maximum errors by the two methods appear when τ（τ is nondimensional time） equals 0.1 near the boundaries of ζ =1 （ζ is nondimensional space coordinate） and ζ =4. As τ increases, the error decreases gradually, and when τ =5 the results of both solutions have almost no change with the variation of coordinate 4.展开更多
Based on the hazard development mechanism, a water solution area is closely related to the supporting effect of pressure-bearing water, the relaxing and collapsing effect of orebody interlayer, the collapsing effect o...Based on the hazard development mechanism, a water solution area is closely related to the supporting effect of pressure-bearing water, the relaxing and collapsing effect of orebody interlayer, the collapsing effect of thawless material in orebody, filling effect caused by cubical expansibility of hydrate crystallization and uplifting effect of hard rock layer over cranny belt. The movement and deformation of ground surface caused by underground water solution mining is believed to be much weaker than that caused by well lane mining, which can be predicted by the stochastic medium theory method. On the basis of analysis on the engineering practice of water solution mining, its corresponding parameters can be obtained from the in-site data of the belt water and sand filling mining in engineering analog approach.展开更多
According to the two-dimensional(2-D) thermo-elasticity theory, the exact elasticity solution of the simply supported laminated beams subjected to thermo-loads was studied. An analytical method was presented to obtain...According to the two-dimensional(2-D) thermo-elasticity theory, the exact elasticity solution of the simply supported laminated beams subjected to thermo-loads was studied. An analytical method was presented to obtain the temperature, displacement and stress fields in the beam. Firstly, the general solutions of temperature, displacements and stresses for a single-layered simply supported beam were obtained by solving the 2-D heat conduction equation and the 2-D elasticity equations, respectively. Then, based on the continuity of temperature, heat flux, displacements and stresses on the interface of two adjacent layers, the formulae of temperature, displacements and stresses between the lowest layer and the top layer of the beam were derived out in a recurrent manner. Finally, the unknown coefficients in the solutions were determined by the use of the upper surface and lower surface conditions of the beam. The distributions of temperature, displacement and stress in the beam were obtained by substituting these coefficients back to the recurrence formulae and the solutions. The excellent convergence of the present method has been demonstrated and the results obtained by the present method agree well with those from the finite element method. The effects of surface temperatures, thickness, layer number and material properties of the plate on the temperature distribution were discussed in detail. Numerical results reveal that the displacements and stresses monotonically increase with the increase of surface temperatures. In particular, the horizontal stresses are discontinuous at the interface.展开更多
This paper carefully analyses the possibility of that a generalized KdV equa- tion possesses a family of solutions of the finite series form,and thus a kind of its explicit complex exact solutions is cxhibited.These c...This paper carefully analyses the possibility of that a generalized KdV equa- tion possesses a family of solutions of the finite series form,and thus a kind of its explicit complex exact solutions is cxhibited.These complex solutions constitute real exact so- lutions to th e system of complex form of the above generalized KdV equation.展开更多
With the expression theorem of convex polyhedron, this paper gives the general expression for the solutions in standard linear programming problems. And the calculation procedures in determining the optimal solutions ...With the expression theorem of convex polyhedron, this paper gives the general expression for the solutions in standard linear programming problems. And the calculation procedures in determining the optimal solutions are also given.展开更多
By using the extended F-expansion method, the exact solutions,including periodic wave solutions expressed by Jacobi elliptic functions, for (2+1)-dimensional nonlinear Schrdinger equation are derived. In the limit c...By using the extended F-expansion method, the exact solutions,including periodic wave solutions expressed by Jacobi elliptic functions, for (2+1)-dimensional nonlinear Schrdinger equation are derived. In the limit cases, the solitary wave solutions and the other type of traveling wave solutions for the system are obtained.展开更多
Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of...Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.展开更多
After diffusion couples were annealed at 260-360 ℃, the concentration profiles of Zn element were measured by EPMA. It has been first quantitatively determined by Matano method that the interdiffusion coefficient in ...After diffusion couples were annealed at 260-360 ℃, the concentration profiles of Zn element were measured by EPMA. It has been first quantitatively determined by Matano method that the interdiffusion coefficient in A1-Zn fee solid solution containing high Zn contents is remarkably decreased due to the small addition of Cu. Also, the interdiffusion coefficient in A1-Zn fee solid solution clearly increases with the increasing of Zn concentration. The interdiffusion activity energy remarkably decreases with the increasing of Zn contents. On the other hand, the interdiffusion activity energy markedly increases due to the small addition of Cu in the A1 Zn fcc solid solution containing high Zn contents.展开更多
基金Supported by the National Natural Science Foundation of China(11671403,11671236,12101192)Henan Provincial General Natural Science Foundation Project(232300420113)。
文摘In this paper,we mainly focus on a type of nonlinear Choquard equations with nonconstant potential.Under appropriate hypotheses on potential function and nonlinear terms,we prove that the above Choquard equation with prescribed 2-norm has some normalized solutions by introducing variational methods.
基金Supported by National Natural Science Foundation of China(11671403,11671236)Henan Provincial General Natural Science Foundation Project(232300420113)。
文摘This article studies a class of nonlinear Kirchhoff equations with exponential critical growth,trapping potential,and perturbation.Under appropriate assumptions about f and h,the article obtained the existence of normalized positive solutions for this equation via the Trudinger-Moser inequality and variational methods.Moreover,these solutions are also ground state solutions.Additionally,the article also characterized the asymptotic behavior of solutions.The results of this article expand the research of relevant literature.
文摘This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.
文摘The solution diffusion coefficient is a great important intrinsical parameter in crystal growth.On earth,it is impossible to accurately determine the diffusion coefficient since there is nature convection.One of the marked charateristics of space-crystal growth is to eleminate nature convection,so that purely diffusion-controlled condition of crystal growth could be realized and precise measurement of the diffusion coefficient should be approved.
文摘Gradiently denitrated gun propellant(GDGP)prepared by a“gradient denitration”strategy is obviously superior in progressive burning performance to the traditional deterred gun propellant.Currently,the preparation of GDGP employed a tedious two-step method involving organic solvents,which hinders the large-scale preparation of GDGP.In this paper,GDGP was successfully prepared via a novelty and environmentally friendly one-step method.The obtained samples were characterized by FT-IR,Raman,SEM and XPS.The results showed that the content of nitrate groups gradiently increased from the surface to the core in the surface layer of GDGP and the surface layer of GDGP exhibited a higher compaction than that of raw gun propellant,with a well-preserved nitrocellulose structure.The denitration process enabled the propellant surface with regressive energy density and good progressive burning performance,as confirmed by oxygen bomb and closed bomb test.At the same time,the effects of different solvents on the component loss of propellant were compared.The result showed that water caused the least component loss.Finally,the stability of GDGP was confirmed by methyl-violet test.This work not only provided environmentally friendly,simple and economic preparation of GDGP,but also confirmed the stability of GDGP prepared by this method.
基金Projects(51078355,50938008) supported by the National Natural Science Foundation of ChinaProject(CX2011B093) supported by the Doctoral Candidate Research Innovation Program of Hunan Province, ChinaProject(20117Q008) supported by the Basic Scientific Research Funds for Central Universities of China
文摘Considering three longitudinal displacement functions and uniform axial displacement functions for shear lag effect and uniform axial deformation of thin-walled box girder with varying depths,a simple and efficient method with high precision to analyze the shear lag effect of thin-walled box girders was proposed.The governing differential equations and boundary conditions of the box girder under lateral loading were derived based on the energy-variational method,and closed-form solutions to stress and deflection corresponding to lateral loading were obtained.Analysis and calculations were carried out with respect to a trapezoidal box girder under concentrated loading or uniform loading and a rectangular box girder under concentrated loading.The analytical results were compared with numerical solutions derived according to the high order finite strip element method and the experimental results.The investigation shows that the closed-form solution is in good agreement with the numerical solutions derived according to the high order finite strip method and the experimental results,and has good stability.Because of the shear lag effect,the stress in cross-section centroid is no longer zero,thus it is not reasonable enough to assume that the strain in cross-section centroid is zero without considering uniform axial deformation.
基金Project(2014M560652)supported by China Postdoctoral Science FoundationProjects(2011CB013802,2013CB036004)supported by the National Basic Research Program of China
文摘By combining the results of laboratory model tests with relevant flow rules, the failure mode of shallow unsymmetrical loading tunnels and the corresponding velocity field were established. According to the principle of virtual power, the upper bound solution for surrounding rock pressure of shallow unsymmetrical loading tunnel was derived and verified by an example. The results indicate that the calculated results of the derived upper bound method for surrounding rock pressure of shallow unsymmetrical loading tunnels are relatively close to those of the existing "code method" and test results, which means that the proposed method is feasible. The current code method underestimates the unsymmetrical loading feature of surrounding rock pressure of shallow unsymmetrical loading tunnels, so it is unsafe; when the burial depth is less or greater than two times of the tunnel span and the unsymmetrical loading angle is less than 45°, the upper bound method or the average value of the results calculated by the upper bound method and code method respectively, is comparatively reasonable. When the burial depth is greater than two times of the tunnel span and the unsymmetrical loading angle is greater than 45°, the code method is more suitable.
文摘The optimization problem is considered in which the objective function is pseudolinear（both pseudoconvex and pseudoconcave） and the constraints are linear. The general expression for the optimal solutions to the problem is derived with the representation theorem of polyhedral sets, and the uniqueness condition of the optimal solution and the computational procedures to determine all optimal solutions （if the uniqueness condition is not satisfied ） are provided. Finally, an illustrative example is also given.
基金Project(2012CB026205)supported by the National Basic Research Program of ChinaProjects(51608264,51778289)supported by the National Natural Science Foundation of ChinaProject(2014Y01)supported by the Transportation Science and Technology Project of Jiangsu Province,China
文摘An exact solution for simply-supported laminated beams with material properties variable with temperature under a combination of uniform thermo-load and mechanical loads was investigated,based on the two-dimensional(2-D)thermo-elasticity theory.Firstly,the beam was divided into a series of layers with uniform material properties along the interfaces of the beam.The uniform thermo-load acted on each layer was transformed into a combination of the normal surface forces acted at the two ends and the transverse thermo-load.Secondly,the state space method was employed to obtain the general solutions of displacements and stresses in an arbitrary layer.Thirdly,based on the interfacial continuity conditions between adjacent layers,the relations of displacement and stress components between the top and bottom layers of the beam were recursively derived by use of the transfer-matrix method.The unknowns in the solutions can be solved by the mechanical loads acted on the top and bottom surfaces.The convergence of the present solutions was checked.The comparative study of the present solutions with the Timoshenko’s solutions and the finite element(FE)solutions was carried out.The effects of material properties variable with temperature on the thermo-elastic behavior of laminated beams were discussed in detail.
基金Projects(50576007,50876016) supported by the National Natural Science Foundation of ChinaProjects(20062180) supported by the National Natural Science Foundation of Liaoning Province,China
文摘To make heat conduction equation embody the essence of physical phenomenon under study, dimensionless factors were introduced and the transient heat conduction equation and its boundary conditions were transformed to dimensionless forms. Then, a theoretical solution model of transient heat conduction problem in one-dimensional double-layer composite medium was built utilizing the natural eigenfunction expansion method. In order to verify the validity of the model, the results of the above theoretical solution were compared with those of finite element method. The results by the two methods are in a good agreement. The maximum errors by the two methods appear when τ（τ is nondimensional time） equals 0.1 near the boundaries of ζ =1 （ζ is nondimensional space coordinate） and ζ =4. As τ increases, the error decreases gradually, and when τ =5 the results of both solutions have almost no change with the variation of coordinate 4.
基金Project(40404001) supported by the National Natural Science Foundation of China
文摘Based on the hazard development mechanism, a water solution area is closely related to the supporting effect of pressure-bearing water, the relaxing and collapsing effect of orebody interlayer, the collapsing effect of thawless material in orebody, filling effect caused by cubical expansibility of hydrate crystallization and uplifting effect of hard rock layer over cranny belt. The movement and deformation of ground surface caused by underground water solution mining is believed to be much weaker than that caused by well lane mining, which can be predicted by the stochastic medium theory method. On the basis of analysis on the engineering practice of water solution mining, its corresponding parameters can be obtained from the in-site data of the belt water and sand filling mining in engineering analog approach.
基金Project(2012CB026205)supported by the National Basic Research Program of ChinaProject(51238003)supported by the National Natural Science Foundation of ChinaProject(2014Y01)supported by the Transportation Department of Jiangsu Province,China
文摘According to the two-dimensional(2-D) thermo-elasticity theory, the exact elasticity solution of the simply supported laminated beams subjected to thermo-loads was studied. An analytical method was presented to obtain the temperature, displacement and stress fields in the beam. Firstly, the general solutions of temperature, displacements and stresses for a single-layered simply supported beam were obtained by solving the 2-D heat conduction equation and the 2-D elasticity equations, respectively. Then, based on the continuity of temperature, heat flux, displacements and stresses on the interface of two adjacent layers, the formulae of temperature, displacements and stresses between the lowest layer and the top layer of the beam were derived out in a recurrent manner. Finally, the unknown coefficients in the solutions were determined by the use of the upper surface and lower surface conditions of the beam. The distributions of temperature, displacement and stress in the beam were obtained by substituting these coefficients back to the recurrence formulae and the solutions. The excellent convergence of the present method has been demonstrated and the results obtained by the present method agree well with those from the finite element method. The effects of surface temperatures, thickness, layer number and material properties of the plate on the temperature distribution were discussed in detail. Numerical results reveal that the displacements and stresses monotonically increase with the increase of surface temperatures. In particular, the horizontal stresses are discontinuous at the interface.
文摘This paper carefully analyses the possibility of that a generalized KdV equa- tion possesses a family of solutions of the finite series form,and thus a kind of its explicit complex exact solutions is cxhibited.These complex solutions constitute real exact so- lutions to th e system of complex form of the above generalized KdV equation.
文摘With the expression theorem of convex polyhedron, this paper gives the general expression for the solutions in standard linear programming problems. And the calculation procedures in determining the optimal solutions are also given.
文摘By using the extended F-expansion method, the exact solutions,including periodic wave solutions expressed by Jacobi elliptic functions, for (2+1)-dimensional nonlinear Schrdinger equation are derived. In the limit cases, the solitary wave solutions and the other type of traveling wave solutions for the system are obtained.
基金Supported by the National Natural Science Foundation of China(12001395)the special fund for Science and Technology Innovation Teams of Shanxi Province(202204051002018)+1 种基金Research Project Supported by Shanxi Scholarship Council of China(2022-169)Graduate Education Innovation Project of Taiyuan Normal University(SYYJSYC-2314)。
文摘Quadratic matrix equations arise in many elds of scienti c computing and engineering applications.In this paper,we consider a class of quadratic matrix equations.Under a certain condition,we rst prove the existence of minimal nonnegative solution for this quadratic matrix equation,and then propose some numerical methods for solving it.Convergence analysis and numerical examples are given to verify the theories and the numerical methods of this paper.
文摘After diffusion couples were annealed at 260-360 ℃, the concentration profiles of Zn element were measured by EPMA. It has been first quantitatively determined by Matano method that the interdiffusion coefficient in A1-Zn fee solid solution containing high Zn contents is remarkably decreased due to the small addition of Cu. Also, the interdiffusion coefficient in A1-Zn fee solid solution clearly increases with the increasing of Zn concentration. The interdiffusion activity energy remarkably decreases with the increasing of Zn contents. On the other hand, the interdiffusion activity energy markedly increases due to the small addition of Cu in the A1 Zn fcc solid solution containing high Zn contents.
