Mastering the influence laws of parameters on the solution structure of nonlinear systems is the basis of carrying out vibration isolation and control.Many researches on solution structure and bifurcation phenomenon i...Mastering the influence laws of parameters on the solution structure of nonlinear systems is the basis of carrying out vibration isolation and control.Many researches on solution structure and bifurcation phenomenon in parameter spaces are carried out broadly in many fields,and the research on nonlinear gear systems has attracted the attention of many scholars.But there is little study on the solution domain boundary of nonlinear gear systems.For a periodic non-autonomous nonlinear dynamic system with several control parameters,a solution domain boundary analysis method of nonlinear systems in parameter spaces is proposed,which combines the cell mapping method based on Poincaré point mapping in phase spaces with the domain decomposition technique of parameter spaces.The cell mapping is known as a global analysis method to analyze the global behavior of a nonlinear dynamic system with finite dimensions,and the basic idea of domain decomposition techniques is to divide and rule.The method is applied to analyze the solution domain boundaries in parameter spaces of a nonlinear gear system.The distribution of different period domains,chaos domain and the domain boundaries between different period domains and chaotic domain are obtained in control parameter spaces constituted by meshing damping ratio with excitation frequency,fluctuation coefficient of meshing stiffness and average exciting force respectively by calculation.The calculation results show that as the meshing damping increases,the responses of the system change towards a single motion,while the variations of the excitation frequency,meshing stiffness and exciting force make the solution domain presenting diversity.The proposed research contribution provides evidence for vibration control and parameter design of the gear system,and confirms the validity of the solution domain boundary analysis method.展开更多
The magnetic prospection is one of the most useful methods to determine buried geological structures such as shallow fracture zones. The investigation of vertical and horizontal gradient and total magnetic field varia...The magnetic prospection is one of the most useful methods to determine buried geological structures such as shallow fracture zones. The investigation of vertical and horizontal gradient and total magnetic field variations over geological structures, which have been used for many years, may reveal their locations, geometries and physical characteristics. In this study, a proposed iterative 3-D rectangular prismatic model inversion algorithm was modified to interpret vertical magnetic gradient data defining the boundaries and the physical parameters of the anomalous structure. Vertical magnetic gradient measurements were carried out at the Tuzla fault, an active fault system located along NE-SW direction in Izmir (Turkey). Boundary analysis studies were applied to data in order to obtain boundaries of the structures, afterwards the inversion process was carried out considering these geometries. As a result, location, direction and other physical and geometrical features of the fault are achieved.展开更多
The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic diff...The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.展开更多
The effects of axisymmetric flow of a Powell-Eyring fluid over an impermeable radially stretching surface are presented. Characteristics of the heat transfer process are analyzed with a more realistic condition named ...The effects of axisymmetric flow of a Powell-Eyring fluid over an impermeable radially stretching surface are presented. Characteristics of the heat transfer process are analyzed with a more realistic condition named the convective boundary condition. Governing equations for the flow problem are derived by the boundary layer approximations. The modeled highly coupled partial differential system is converted into a system of ordinary differential equations with acceptable similarity transformations. The convergent series solutions for the resulting system are constructed and analyzed. Optimal values are obtained and presented in a numerical form using an optimal homotopy analysis method (OHAM). The rheological characteristics of different parameters of the velocity and temperature profiles are presented graphically. Tabular variations of the skin friction coefficient and the Nusselt number are also calculated. It is observed that the temperature distribution shows opposite behavior for Prandtl and Biot numbers. Furthermore, the rate of heating/cooling is higher for both the Prandtl and Biot numbers.展开更多
In this paper, the boundary element method is applied to investigate the internal state of stress of autofretted tube with notch and the calculated results are important in the practical design.
This paper discusses three kinds of BVPs of elliptic system (E1) on the complementary domain of disk. The results, based on Fourier analysis method. not only enrich the results of reference [1]. but also are of both t...This paper discusses three kinds of BVPs of elliptic system (E1) on the complementary domain of disk. The results, based on Fourier analysis method. not only enrich the results of reference [1]. but also are of both theoretical and practical interest.展开更多
After volume fracturing of horizontal wells in shale gas reservoir, an extremely complex fracture system is formed. The space area of the fracture system is the reservoir reconstruction volume of shale gas reservoir. ...After volume fracturing of horizontal wells in shale gas reservoir, an extremely complex fracture system is formed. The space area of the fracture system is the reservoir reconstruction volume of shale gas reservoir. The geometric parameters such as crack length, crack width, crack height, and characteristic parameters such as crack permeability and fracture conductivity proposed for a single crack in conventional fracturing are insufficient to describe and characterize the complex network fracture system after volume fracturing. In this paper, the discrete fracture modeling method is used to establish the volume fracturing network fracture model of horizontal wells in shale gas reservoir by using the random modeling method within the determined reservoir space. The model is random and selective, and can fully provide different forms of volume fracturing fracture expansion, such as conventional fracture morphology, line network model and arbitrarily distributed network fractures. The research results provide a theoretical basis for the development plan and stimulation plan of shale gas reservoir, and have important reference value and significance for other unconventional gas reservoir fracturing.展开更多
A new numerical method,scaled boundary isogeometric analysis(SBIGA)combining the concept of the scaled boundary finite element method(SBFEM)and the isogeometric analysis(IGA),is proposed in this study for 2D elastosta...A new numerical method,scaled boundary isogeometric analysis(SBIGA)combining the concept of the scaled boundary finite element method(SBFEM)and the isogeometric analysis(IGA),is proposed in this study for 2D elastostatic problems with both homogenous and inhomogeneous essential boundary conditions.Scaled boundary isogeometric transformation is established at a specified scaling center with boundary isogeometric representation identical to the design model imported from CAD system,which can be automatically refined without communication with the original system and keeping geometry invariability.The field variable,that is,displacement,is constructed by the same basis as boundary isogeometric description keeping analytical features in radial direction.A Lagrange multiplier scheme is suggested to impose the inhomogeneous essential boundary conditions.The new proposed method holds the semi-analytical feature inherited from SBFEM,that is,discretization only on boundaries rather than the entire domain,and isogeometric boundary geometry from IGA,which further increases the accuracy of the solution.Numerical examples,including circular cavity in full plane,Timoshenko beam with inhomogenous boundary conditions and infinite plate with circular hole subjected to remotely tension,demonstrate that SBIGA can be applied efficiently to elastostatic problems with various boundary conditions,and powerful in accuracy of solution and less degrees of freedom(DOF)can be achieved in SBIGA than other methods.展开更多
For the multi-frequency acoustic analysis, a series expansion method has been introduced to reduce the computation time of the frequency-independent parts, but the Runge phenomenon will arise when this method is emplo...For the multi-frequency acoustic analysis, a series expansion method has been introduced to reduce the computation time of the frequency-independent parts, but the Runge phenomenon will arise when this method is employed in high frequency band. Therefore, this method is improved by analyzing the application condition and proposing the selection principle of the series truncation number. The argument interval can be adjusted with the wavenumber factor. Therefore, the problem of unstable numeration and poor precision can be solved, and the application scope of this method is expanded. The numerical example of acoustic radiation shows that the improved method is correct for acoustic analysis in wider frequency band with less series truncation number and computation amount.展开更多
The near crack line field analysis method has been used io investigate into theexact elastic-plastic solutions of a mode II crack under plane stress condilion in anelastic-perfectly plastic solid. The assumptions of t...The near crack line field analysis method has been used io investigate into theexact elastic-plastic solutions of a mode II crack under plane stress condilion in anelastic-perfectly plastic solid. The assumptions of the usual small scale yielding theory.hare been completely. dbandoned and the correct formulations of matching conditionsat the elaslic-plastic boundary. have been given. By, matching the general solution ofthe plastic slress field (bul not the special solution used to be adopted) with the exactelastic stress field (but not the crack tip K-dominant field) at the elastic-plasticboundary, near the crack line, the plastic stresses. the length of the plastic =one and theunit normal vector of the elastic-plastic boundary.which are sufficiently precise near the crack line region ,have been given.展开更多
The near crack line field analysis method has been used to investigate into the exact elastic-plastic solutions of a mode Ⅱ crack under plane stress condition in anelastic-perfectly plastic solid. The assumptions of ...The near crack line field analysis method has been used to investigate into the exact elastic-plastic solutions of a mode Ⅱ crack under plane stress condition in anelastic-perfectly plastic solid. The assumptions of the usual small scale yielding theoryhave been completely. dbandoned and the correct .formulations of matching conditionsat the elastic-plastic boundary have been given. By matching the general solution of the plastic stress field (but not the special solution used to be adopted) will the exactelastic stress field (but not the crack tip K-dominant field) at the elastic-plasticboundary. near the crack line, the plastic .stresses, the length of the plastic zone and theunit normal vector of the elaslic-plastic boundary. which sufficiently precise nearthe crack line region, hare been given.展开更多
An integration processing system of three-dimensional laser scanning information visualization in goaf was developed. It is provided with multiple functions, such as laser scanning information management for goaf, clo...An integration processing system of three-dimensional laser scanning information visualization in goaf was developed. It is provided with multiple functions, such as laser scanning information management for goaf, cloud data de-noising optimization, construction, display and operation of three-dimensional model, model editing, profile generation, calculation of goaf volume and roof area, Boolean calculation among models and interaction with the third party soft ware. Concerning this system with a concise interface, plentiful data input/output interfaces, it is featured with high integration, simple and convenient operations of applications. According to practice, in addition to being well-adapted, this system is favorably reliable and stable.展开更多
Accurate edge localization of bilevel images is of primary importance in barcode decoding.In the sub-pixel edge location algorithm for bilevel images,the bilevel image(barcode) imaging process is modeled as a square...Accurate edge localization of bilevel images is of primary importance in barcode decoding.In the sub-pixel edge location algorithm for bilevel images,the bilevel image(barcode) imaging process is modeled as a square wave convoluted with a Gaussian kernel and then sampled discretely by pixel arrays.Based on the gray levels of the pixels,assumed sub-pixel edge locations are set and adjusted so that the discrepancy of the theoretical gray level of pixels and the actual gray level of pixels reaches the minimum and then the best approximation of the actual sub-pixel edges of the bilevel image is obtained.Examples are presented to illustrate the techniques of the algorithm which can solve the problems of edge location or signal recovery of bilevel images in the case of the two features:one is that the support of the Gaussian kernel is comparable to the distance of the adjacent edges;the other is that the distance between the adjacent edges is comparable to the distance of the adjacent pixels.展开更多
In order to analyze the stress and strain fields in the fibers and the matrix in composite materials,a fiber-scale unit cell model is established and the corresponding periodical boundary conditions are introduced.Ass...In order to analyze the stress and strain fields in the fibers and the matrix in composite materials,a fiber-scale unit cell model is established and the corresponding periodical boundary conditions are introduced.Assuming matrix cracking as the failure mode of composite materials,an energy-based fatigue damage parameter and a multiaxial fatigue life prediction method are established.This method only needs the material properties of the fibers and the matrix to be known.After the relationship between the fatigue damage parameter and the fatigue life under any arbitrary test condition is established,the multiaxial fatigue life under any other load condition can be predicted.The proposed method has been verified using two different kinds of load forms.One is unidirectional laminates subjected to cyclic off-axis loading,and the other is filament wound composites subjected to cyclic tension-torsion loading.The fatigue lives predicted using the proposed model are in good agreements with the experimental results for both kinds of load forms.展开更多
We present a mean field study of a propagation-tumover lattice model, which was proposed by Hodges and Crabtree [Proc. Nat. Acad. Sci. 109, 13296 (2012)] for understanding how posttranslational histone marks modulat...We present a mean field study of a propagation-tumover lattice model, which was proposed by Hodges and Crabtree [Proc. Nat. Acad. Sci. 109, 13296 (2012)] for understanding how posttranslational histone marks modulate gene expression in mammalian ceils. The kinetics of the lattice model consists of nucleation, propagation and turnover mechanisms, and exhibits second-order phase transition for the histone marking domain. We showed rigorously that the dynamics essentially depends on a non-dimensional parameter k = k+/k-, the ratio between the propagation and turnover rates, which has been observed in the simulations. We then studied the lowest order mean field approximation, and observed the phase transition with an analytically obtained critical parameter. The boundary layer analysis was utilized to investigate the structure of the decay profile of the mark density. We also studied the higher order mean field approximation to achieve sharper estimate of the critical transition parameter and more detailed features. The comparison between the simulation and theoretical results shows the validity of our theory.展开更多
Agriculture faces the dual challenges of food security and environmental sustainability.Here,we investigate current maize production at the field scale,analyze the yield gaps and impacting factors,and recommend measur...Agriculture faces the dual challenges of food security and environmental sustainability.Here,we investigate current maize production at the field scale,analyze the yield gaps and impacting factors,and recommend measures for sustainably closing yield gaps.An experiment was conducted on a 3.9-ha maize seed production field in arid north-west China,managed with border and drip irrigation,respectively,in 2015 and 2016.The relative yield reached 70%in both years.However,drip irrigation saved 227 mm irrigation water during a drier growing season compared with traditional border irrigation,accounting for 44%of the maize evapotranspiration(ET).Yield variability under drip irrigation was12.1%,lower than the 18.8%under border irrigation.Boundary line analysis indicates that a relative yield increase of 8%to 10%might be obtained by optimizing the yield-limiting factors.Plant density and soil available water content and available nitrogen were the three major factors involved.In conclusion,closing yield gaps with agricultural sustainability may be realized by optimizing agronomic,irrigation and fertilizer management,using water-saving irrigation methods and using site-specific management.展开更多
基金supported by National Hi-tech Research and Development Program of China (863 Program,Grant No.2009AA04Z404)
文摘Mastering the influence laws of parameters on the solution structure of nonlinear systems is the basis of carrying out vibration isolation and control.Many researches on solution structure and bifurcation phenomenon in parameter spaces are carried out broadly in many fields,and the research on nonlinear gear systems has attracted the attention of many scholars.But there is little study on the solution domain boundary of nonlinear gear systems.For a periodic non-autonomous nonlinear dynamic system with several control parameters,a solution domain boundary analysis method of nonlinear systems in parameter spaces is proposed,which combines the cell mapping method based on Poincaré point mapping in phase spaces with the domain decomposition technique of parameter spaces.The cell mapping is known as a global analysis method to analyze the global behavior of a nonlinear dynamic system with finite dimensions,and the basic idea of domain decomposition techniques is to divide and rule.The method is applied to analyze the solution domain boundaries in parameter spaces of a nonlinear gear system.The distribution of different period domains,chaos domain and the domain boundaries between different period domains and chaotic domain are obtained in control parameter spaces constituted by meshing damping ratio with excitation frequency,fluctuation coefficient of meshing stiffness and average exciting force respectively by calculation.The calculation results show that as the meshing damping increases,the responses of the system change towards a single motion,while the variations of the excitation frequency,meshing stiffness and exciting force make the solution domain presenting diversity.The proposed research contribution provides evidence for vibration control and parameter design of the gear system,and confirms the validity of the solution domain boundary analysis method.
文摘The magnetic prospection is one of the most useful methods to determine buried geological structures such as shallow fracture zones. The investigation of vertical and horizontal gradient and total magnetic field variations over geological structures, which have been used for many years, may reveal their locations, geometries and physical characteristics. In this study, a proposed iterative 3-D rectangular prismatic model inversion algorithm was modified to interpret vertical magnetic gradient data defining the boundaries and the physical parameters of the anomalous structure. Vertical magnetic gradient measurements were carried out at the Tuzla fault, an active fault system located along NE-SW direction in Izmir (Turkey). Boundary analysis studies were applied to data in order to obtain boundaries of the structures, afterwards the inversion process was carried out considering these geometries. As a result, location, direction and other physical and geometrical features of the fault are achieved.
基金supported by the National Natural Science Foundation of China
文摘The calculations of unsteady flow to a multiple well system with the application of boundary elementmethod (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem ofparabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time varia-ble. The image function of water head H can be solved by BEM. We derived the boundary integral equation ofthe transformed variable H and the discretization form of it, so that there is no need to discretize the bounda-ries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.
文摘The effects of axisymmetric flow of a Powell-Eyring fluid over an impermeable radially stretching surface are presented. Characteristics of the heat transfer process are analyzed with a more realistic condition named the convective boundary condition. Governing equations for the flow problem are derived by the boundary layer approximations. The modeled highly coupled partial differential system is converted into a system of ordinary differential equations with acceptable similarity transformations. The convergent series solutions for the resulting system are constructed and analyzed. Optimal values are obtained and presented in a numerical form using an optimal homotopy analysis method (OHAM). The rheological characteristics of different parameters of the velocity and temperature profiles are presented graphically. Tabular variations of the skin friction coefficient and the Nusselt number are also calculated. It is observed that the temperature distribution shows opposite behavior for Prandtl and Biot numbers. Furthermore, the rate of heating/cooling is higher for both the Prandtl and Biot numbers.
文摘In this paper, the boundary element method is applied to investigate the internal state of stress of autofretted tube with notch and the calculated results are important in the practical design.
文摘This paper discusses three kinds of BVPs of elliptic system (E1) on the complementary domain of disk. The results, based on Fourier analysis method. not only enrich the results of reference [1]. but also are of both theoretical and practical interest.
文摘After volume fracturing of horizontal wells in shale gas reservoir, an extremely complex fracture system is formed. The space area of the fracture system is the reservoir reconstruction volume of shale gas reservoir. The geometric parameters such as crack length, crack width, crack height, and characteristic parameters such as crack permeability and fracture conductivity proposed for a single crack in conventional fracturing are insufficient to describe and characterize the complex network fracture system after volume fracturing. In this paper, the discrete fracture modeling method is used to establish the volume fracturing network fracture model of horizontal wells in shale gas reservoir by using the random modeling method within the determined reservoir space. The model is random and selective, and can fully provide different forms of volume fracturing fracture expansion, such as conventional fracture morphology, line network model and arbitrarily distributed network fractures. The research results provide a theoretical basis for the development plan and stimulation plan of shale gas reservoir, and have important reference value and significance for other unconventional gas reservoir fracturing.
基金supported by the National Natural Science Foundation of China(Grant Nos.51138001,51009019,51109134)
文摘A new numerical method,scaled boundary isogeometric analysis(SBIGA)combining the concept of the scaled boundary finite element method(SBFEM)and the isogeometric analysis(IGA),is proposed in this study for 2D elastostatic problems with both homogenous and inhomogeneous essential boundary conditions.Scaled boundary isogeometric transformation is established at a specified scaling center with boundary isogeometric representation identical to the design model imported from CAD system,which can be automatically refined without communication with the original system and keeping geometry invariability.The field variable,that is,displacement,is constructed by the same basis as boundary isogeometric description keeping analytical features in radial direction.A Lagrange multiplier scheme is suggested to impose the inhomogeneous essential boundary conditions.The new proposed method holds the semi-analytical feature inherited from SBFEM,that is,discretization only on boundaries rather than the entire domain,and isogeometric boundary geometry from IGA,which further increases the accuracy of the solution.Numerical examples,including circular cavity in full plane,Timoshenko beam with inhomogenous boundary conditions and infinite plate with circular hole subjected to remotely tension,demonstrate that SBIGA can be applied efficiently to elastostatic problems with various boundary conditions,and powerful in accuracy of solution and less degrees of freedom(DOF)can be achieved in SBIGA than other methods.
基金supported by the National Natural Science Foundation of China(51379083,51479079,51579109)the Specialized Research Fund for the Doctoral Program of Higher Education of China(20120142110051)
文摘For the multi-frequency acoustic analysis, a series expansion method has been introduced to reduce the computation time of the frequency-independent parts, but the Runge phenomenon will arise when this method is employed in high frequency band. Therefore, this method is improved by analyzing the application condition and proposing the selection principle of the series truncation number. The argument interval can be adjusted with the wavenumber factor. Therefore, the problem of unstable numeration and poor precision can be solved, and the application scope of this method is expanded. The numerical example of acoustic radiation shows that the improved method is correct for acoustic analysis in wider frequency band with less series truncation number and computation amount.
文摘The near crack line field analysis method has been used io investigate into theexact elastic-plastic solutions of a mode II crack under plane stress condilion in anelastic-perfectly plastic solid. The assumptions of the usual small scale yielding theory.hare been completely. dbandoned and the correct formulations of matching conditionsat the elaslic-plastic boundary. have been given. By, matching the general solution ofthe plastic slress field (bul not the special solution used to be adopted) with the exactelastic stress field (but not the crack tip K-dominant field) at the elastic-plasticboundary, near the crack line, the plastic stresses. the length of the plastic =one and theunit normal vector of the elastic-plastic boundary.which are sufficiently precise near the crack line region ,have been given.
文摘The near crack line field analysis method has been used to investigate into the exact elastic-plastic solutions of a mode Ⅱ crack under plane stress condition in anelastic-perfectly plastic solid. The assumptions of the usual small scale yielding theoryhave been completely. dbandoned and the correct .formulations of matching conditionsat the elastic-plastic boundary have been given. By matching the general solution of the plastic stress field (but not the special solution used to be adopted) will the exactelastic stress field (but not the crack tip K-dominant field) at the elastic-plasticboundary. near the crack line, the plastic .stresses, the length of the plastic zone and theunit normal vector of the elaslic-plastic boundary. which sufficiently precise nearthe crack line region, hare been given.
基金Project(51274250)supported by the National Natural Science Foundation of ChinaProject(2012BAK09B02-05)supported by the National Key Technology R&D Program during the 12th Five-year Plan of China
文摘An integration processing system of three-dimensional laser scanning information visualization in goaf was developed. It is provided with multiple functions, such as laser scanning information management for goaf, cloud data de-noising optimization, construction, display and operation of three-dimensional model, model editing, profile generation, calculation of goaf volume and roof area, Boolean calculation among models and interaction with the third party soft ware. Concerning this system with a concise interface, plentiful data input/output interfaces, it is featured with high integration, simple and convenient operations of applications. According to practice, in addition to being well-adapted, this system is favorably reliable and stable.
基金Supported by the Postdoctoral Science Fund of China (20070410940)the Open Fund of Liaoning Key Laboratory of Intelligent Information Processing, Dalian University (2005-8)
文摘Accurate edge localization of bilevel images is of primary importance in barcode decoding.In the sub-pixel edge location algorithm for bilevel images,the bilevel image(barcode) imaging process is modeled as a square wave convoluted with a Gaussian kernel and then sampled discretely by pixel arrays.Based on the gray levels of the pixels,assumed sub-pixel edge locations are set and adjusted so that the discrepancy of the theoretical gray level of pixels and the actual gray level of pixels reaches the minimum and then the best approximation of the actual sub-pixel edges of the bilevel image is obtained.Examples are presented to illustrate the techniques of the algorithm which can solve the problems of edge location or signal recovery of bilevel images in the case of the two features:one is that the support of the Gaussian kernel is comparable to the distance of the adjacent edges;the other is that the distance between the adjacent edges is comparable to the distance of the adjacent pixels.
基金the supports from the Jiangsu Province Key Laboratory of Aerospace Power System of China(No.NJ20140019)the National Natural Science Foundation of China(No.51205190)
文摘In order to analyze the stress and strain fields in the fibers and the matrix in composite materials,a fiber-scale unit cell model is established and the corresponding periodical boundary conditions are introduced.Assuming matrix cracking as the failure mode of composite materials,an energy-based fatigue damage parameter and a multiaxial fatigue life prediction method are established.This method only needs the material properties of the fibers and the matrix to be known.After the relationship between the fatigue damage parameter and the fatigue life under any arbitrary test condition is established,the multiaxial fatigue life under any other load condition can be predicted.The proposed method has been verified using two different kinds of load forms.One is unidirectional laminates subjected to cyclic off-axis loading,and the other is filament wound composites subjected to cyclic tension-torsion loading.The fatigue lives predicted using the proposed model are in good agreements with the experimental results for both kinds of load forms.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11174011, 11021463, 11421101, and 91530322)
文摘We present a mean field study of a propagation-tumover lattice model, which was proposed by Hodges and Crabtree [Proc. Nat. Acad. Sci. 109, 13296 (2012)] for understanding how posttranslational histone marks modulate gene expression in mammalian ceils. The kinetics of the lattice model consists of nucleation, propagation and turnover mechanisms, and exhibits second-order phase transition for the histone marking domain. We showed rigorously that the dynamics essentially depends on a non-dimensional parameter k = k+/k-, the ratio between the propagation and turnover rates, which has been observed in the simulations. We then studied the lowest order mean field approximation, and observed the phase transition with an analytically obtained critical parameter. The boundary layer analysis was utilized to investigate the structure of the decay profile of the mark density. We also studied the higher order mean field approximation to achieve sharper estimate of the critical transition parameter and more detailed features. The comparison between the simulation and theoretical results shows the validity of our theory.
基金funded by the National Natural Science Foundation of China(51725904,51621061,51861125103)the Discipline Innovative Engineering Plan(111 Program,B14002)。
文摘Agriculture faces the dual challenges of food security and environmental sustainability.Here,we investigate current maize production at the field scale,analyze the yield gaps and impacting factors,and recommend measures for sustainably closing yield gaps.An experiment was conducted on a 3.9-ha maize seed production field in arid north-west China,managed with border and drip irrigation,respectively,in 2015 and 2016.The relative yield reached 70%in both years.However,drip irrigation saved 227 mm irrigation water during a drier growing season compared with traditional border irrigation,accounting for 44%of the maize evapotranspiration(ET).Yield variability under drip irrigation was12.1%,lower than the 18.8%under border irrigation.Boundary line analysis indicates that a relative yield increase of 8%to 10%might be obtained by optimizing the yield-limiting factors.Plant density and soil available water content and available nitrogen were the three major factors involved.In conclusion,closing yield gaps with agricultural sustainability may be realized by optimizing agronomic,irrigation and fertilizer management,using water-saving irrigation methods and using site-specific management.