This article lists all possible conversion equations of influence coefficients(ICs)with different formats when some of these are known to handle two-plane balancing tasks with no or less trial weight runs.There are tw...This article lists all possible conversion equations of influence coefficients(ICs)with different formats when some of these are known to handle two-plane balancing tasks with no or less trial weight runs.There are two approaches for two-plane balancing.One can treat it as a multi-plane balance problem involving a 2×2 matrix of complex ICs where two direct ICs along with two cross-effect ICs are generated so that correction weights at one or two balance planes can be determined.One can also apply a static pare(in-phase)and/or couple pair(180 degrees out-of-phase)weights for balancing.The latter approach has been used quite often in the field,especially on steam turbine and generator rotors.Dependent on vibration mode shapes and combinations as well as balance plane accessibility,sometimes applying static or couple pair weights can be a wise choice;other times weights at one or two end planes are needed.There are totally 4 possible sets of IC data due to weights at plane 1,plane 2,static pair,and couple pair.IC data would typically be obtained by applying trial weights followed by trial weight runs.It is found,however,that all these IC data can be converted easily without trial weight runs once any two of 4 sets are known.The above findings and conversion equations have been obtained analytically and verified by experimental results.Real cases are given to demonstratetheirapplications.展开更多
Measurement error of unbalance's vibration response plays a crucial role in calibration and on-line updating of influence coefficient(IC). Focusing on the two problems that the moment estimator of data used in cali...Measurement error of unbalance's vibration response plays a crucial role in calibration and on-line updating of influence coefficient(IC). Focusing on the two problems that the moment estimator of data used in calibration process cannot fulfill the accuracy requirement under small sample and the disturbance of measurement error cannot be effectively suppressed in updating process, an IC calibration and on-line updating method based on hierarchical Bayesian method for automatic dynamic balancing machine was proposed. During calibration process, for the repeatedly-measured data obtained from experiments with different trial weights, according to the fact that measurement error of each sensor had the same statistical characteristics, the joint posterior distribution model for the true values of the vibration response under all trial weights and measurement error was established. During the updating process, information obtained from calibration was regarded as prior information, which was utilized to update the posterior distribution of IC combined with the real-time reference information to implement online updating. Moreover, Gibbs sampling method of Markov Chain Monte Carlo(MCMC) was adopted to obtain the maximum posterior estimation of parameters to be estimated. On the independent developed dynamic balancing testbed, prediction was carried out for multiple groups of data through the proposed method and the traditional method respectively, the result indicated that estimator of influence coefficient obtained through the proposed method had higher accuracy; the proposed updating method more effectively guaranteed the measurement accuracy during the whole producing process, and meantime more reasonably compromised between the sensitivity of IC change and suppression of randomness of vibration response.展开更多
In order to develop a general calculating rotor’s torsional stiffness based on stiffness influence coefficient for different rotor assembling, the calculation method of the torsional stiffness influence coefficient o...In order to develop a general calculating rotor’s torsional stiffness based on stiffness influence coefficient for different rotor assembling, the calculation method of the torsional stiffness influence coefficient of equal thickness disc is researched in this paper at first. Then the torsional stiffness influence coefficient λ of equal thickness disc is fit to a binary curved face and a calculation equation is obtained based on a large quantity of calculating data, which lays the foundation for research on a general calculating method of rotor torsional stiffness. Thirdly a simplified calculation method for equivalent stiffness diameter of stepped equal thickness disc and cone disc in the steam turbine generators is suggested. Finally a general calculating program for calculating rotor’s torsional vibration features is developed, and the torsional vibration features of a verity of steam turbine rotors are calculated for verification. The calculating results show that stiffness influence coefficient λ of equal-thickness disc depends on parameters of B and H, as well as the stiffness influence coefficient λ;and discs with complex structure can be simplified to equal-thickness discs with little error by using the method suggested in this paper;error can be controlled within 1% when equivalent diameter of stiffness is calculated by this method.展开更多
In this paper, we study local influence analysis for Zhang's generalized correlation coefficients and Hotelling's generalized correlation coefficient by using approach. of local influence analysis suggested by...In this paper, we study local influence analysis for Zhang's generalized correlation coefficients and Hotelling's generalized correlation coefficient by using approach. of local influence analysis suggested by Shi (1991), i.e., generalized influence function (GIF) and generalized Cook distance (GCD). An example is given to illustrate our results.展开更多
In order to execute geometric analysis for planar deployable mechanism of scissors unit,the dynamic analysis model of scissor planar deployable structure is created based on the Cartesian coordinate system,the influen...In order to execute geometric analysis for planar deployable mechanism of scissors unit,the dynamic analysis model of scissor planar deployable structure is created based on the Cartesian coordinate system,the influence coefficient is acquired by means of the coordinate transformation,combining the D 'Alembert 's principle with Dynamic-Static method,the dynamic characteristic analysis is completed finally. Moreover,specific calculating examples are adopted to verify the effectiveness of proposed method,and the result shows that the movement of each component of scissors unit mechanism is more smooth during initial deployment stage,however,when the configuration angle θ of unit mechanism is approaching π,some comparative large variations would appear on movement parameters and hinge constraint force.展开更多
Based on the model structure of the influence coefficient method analyzed in depth by matrix theory ,it is explained the reason why the unreasonable and instable correction masses with bigger MSE are obtained by LS in...Based on the model structure of the influence coefficient method analyzed in depth by matrix theory ,it is explained the reason why the unreasonable and instable correction masses with bigger MSE are obtained by LS influence coefficient method when there are correlation planes in the dynamic balancing. It also presencd the new ridge regression method for solving correction masses according to the Tikhonov regularization theory, and described the reason why the ridge regression can eliminate the disadvantage of the LS method. Applying this new method to dynamic balancing of gas turbine, it is found that this method is superior to the LS method when influence coefficient matrix is ill-conditioned,the minimal correction masses and residual vibration are obtained in the dynamic balancing of rotors.展开更多
According to the concept of virtual bending force,a rational explanation for SHOHET's model was presented. Considering the deformation characters of the work rolls in four-high mill,the deformation model of the wo...According to the concept of virtual bending force,a rational explanation for SHOHET's model was presented. Considering the deformation characters of the work rolls in four-high mill,the deformation model of the work roll was regarded as a cantilever beam and new influence coefficients were deduced.The effect of the bending force was taken into account independently. Therefore,the contribution to work roll deflection caused by rolling load,rolling pressure between rolls and bending force can be got from the new formulas.To validate the accuracy of the formulas,the results obtained from the new formulas were compared with those from SHOHET's formulas.It is found that they highly coincide,which illustrates that the formulas are reliable.展开更多
The influence of the pressure transmission medium(PTM)on the excitonic interband transitions in monolayer tungsten diselenide(WSe2)is investigated using photoluminescence(PL)spectra under hydrostatic pressure up...The influence of the pressure transmission medium(PTM)on the excitonic interband transitions in monolayer tungsten diselenide(WSe2)is investigated using photoluminescence(PL)spectra under hydrostatic pressure up to 5GPa.Three kinds of PTMs,condensed argon(Ar),1:1 n-pentane and isopentane mixture(PM),and4:1 methanol and ethanol mixture(MEM,a PTM with polarity),are used.It is found that when either Ar or PM is used as the PTM,the PL peak of exciton related to the direct K-K interband transition shows a pressure-induced blue-shift at a rate of 32±4 or 32±1 meV/GPa,while it turns to be 50±9meV/GPa when MEM is used as the PTM.The indirect A-K interband transition presents almost no shift with increasing pressure up to approximatel.y 5 GPa when Ar and PM are used as the PTM,while it shows a red-shift at the rate of-17±7meV/GPa by using MEM as the PTM.These results reveal that the optical interband transitions of monolayer WSe2 are very sensitive to the polarity of the PTM.The anomalous pressure coefficient obtained using the polar PTM of MEM is ascribed to the existence of hydrogen-like bonds between hydroxyl in MEM and Se atoms under hydrostatic pressure.展开更多
Typical rotors such as those on steam turbine or generator are often supported by two bearings with two balance planes at both ends.Vibrations are monitored by a pair of proximity probes at each bearing.There are two ...Typical rotors such as those on steam turbine or generator are often supported by two bearings with two balance planes at both ends.Vibrations are monitored by a pair of proximity probes at each bearing.There are two approaches to reduce 1X vibration due to unbalance at both ends via balancing with influence coefficient method.The first approach is to treat it as a multiple-plane balancing problem involving 2x2 matrix of complex influence coefficients.The second approach is to treat it as two single-plane balance problems using static(in-phase)and couple(180 degree out-of-phase)components,respectively.Conversion equations of influence coefficients between these two approaches have been found previously by the author.The corresponding spreadsheets that convert influence coefficients between these two formats are presented in the current paper.The paper shows effectiveness of these conversion equations in dealing with real balancing problems in the field.A detailed balance case is presented to demonstrate how the conversion equations are used to reduce vibration effectively.展开更多
The steel slag (SS) permeable concrete was prepared by SS. The influences of the aggregate-cement rate, the aggregate particle size, the water-cement rate, the admixture dosage and other factors on the permeability ...The steel slag (SS) permeable concrete was prepared by SS. The influences of the aggregate-cement rate, the aggregate particle size, the water-cement rate, the admixture dosage and other factors on the permeability coefficient of SS permeable concrete were analyzed. The law of influence was also investigated. The study serves as a technological reference for the construction and design of SS permeable concrete.展开更多
Due to the diversity of work requirements and environment,the number of degrees of freedom(DOFs)and the complexity of structure of industrial robots are constantly increasing.It is difficult to establish the accurate ...Due to the diversity of work requirements and environment,the number of degrees of freedom(DOFs)and the complexity of structure of industrial robots are constantly increasing.It is difficult to establish the accurate dynamical model of industrial robots,which greatly hinders the realization of a stable,fast and accurate trajectory tracking control.Therefore,the general expression of the constraint relation in the explicit dynamic equation of the multi-DOF industrial robot is derived on the basis of solving the Jacobian matrix and Hessian matrix by using the kinematic influence coefficients method.Moreover,an explicit dynamic equation with general constraint relation expression is established based on the Udwadia-Kalaba theory.The problem of increasing the time of establishing constraint relationship when the multi-DOF industrial robots complete complex task constraints is solved.With the SCARA robot as the research object,the simulation results show that the proposed method can provide a new idea for industrial robot system modeling with complex constraints.展开更多
Varying-coefficient single-index model( VCSIM) avoids the so-called "curse of dimensionality " and is flexible enough to include several important statistical models. This paper considers statistical diagnos...Varying-coefficient single-index model( VCSIM) avoids the so-called "curse of dimensionality " and is flexible enough to include several important statistical models. This paper considers statistical diagnosis for VCSIM. First,the parametric estimation equation is established based on empirical likelihood. Then,some diagnosis statistics are defined. At last, an example is given to illustrate all the results.展开更多
Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering appl...Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering applications. A three-dimensional anisotropic heterogeneous steady random seepage model was developed on the basis of the finite element method. A statistical analysis of the distribution characteristics of soil parameters sampled from the main embankment of the Yangtze River in the Southern Jingzhou zone of China was conducted. The Kolomogorov-Smirnov test verified the statistical hypothesis that the permeability coefficient tensor has a Gaussian distribution. With the help of numerical analysis of the stochastic seepage field using the developed model, various statistical and random characteristics of the stochastic seepage field of the main embankment of the Yangtze River in the Southern Jingzhou zone of China were investigated. The model was also examined with statistical testing. Through the introduction of random variation of the upstream and downstream water levels into the model, the effects of the boundary randomness due to variation of the downstream and upstream water levels on the variation of simulated results presented with a vector series of the random seepage field were analyzed. Furthermore, the combined influence of the variation of the soil permeability coefficient and such seepage resistance measures as the cut-off wall and relief ditch on the hydraulic head distribution was analyzed and compared with the results obtained by determinate analysis. Meanwhile, sensitivities of the hydraulic gradient and downstream exit height to the variation of boundary water level were studied. The validity of the simulated results was verified by stochastic testing and measured data. The developed model provides more detail and a full stochastic algorithm to characterize and analyze three-dimensional stochastic seepage field problems.展开更多
A three-dimensional, first order turbulence closure, thermal diffusion model is described in this paper. The governing equations consist of an equation of continuity, three components of momentum, conservation equatio...A three-dimensional, first order turbulence closure, thermal diffusion model is described in this paper. The governing equations consist of an equation of continuity, three components of momentum, conservation equations for salt, temperature and subgridscale energy, and an equation of state. In the model, according to the hypothesis of Kolmogorov and Prandtl, the viscosity coefficient of turbulent flow of homogeneous fluid is related to the local turbulent energy, and the horizontal and vertical exchange coefficients of mass, heat and momentum are computed with the introduction of subgridscale turbulence energy. The governing equations are solved by finite difference techniques. This model is applied to the Jiaozhou bay to predict thermal pollution by the Huangdao power plant. An instantaneous tidal current field is computed, then the distribution of temperature increment is predicted, and finally the effect of wind stress on thermal discharge is discussed.展开更多
In this paper, a varying-coefficient density-ratio model for case-control studies is developed. We investigate the local empirical likelihood diagnosis of varying coefficient density-ratio model for case-control data....In this paper, a varying-coefficient density-ratio model for case-control studies is developed. We investigate the local empirical likelihood diagnosis of varying coefficient density-ratio model for case-control data. The local empirical log-likelihood ratios for the nonparametric coefficient functions are introduced. First, the estimation equations based on empirical likelihood method are established. Then, a few of diagnostic statistics are proposed. At last, we also examine the performance of proposed method for finite sample sizes through simulation studies.展开更多
文摘This article lists all possible conversion equations of influence coefficients(ICs)with different formats when some of these are known to handle two-plane balancing tasks with no or less trial weight runs.There are two approaches for two-plane balancing.One can treat it as a multi-plane balance problem involving a 2×2 matrix of complex ICs where two direct ICs along with two cross-effect ICs are generated so that correction weights at one or two balance planes can be determined.One can also apply a static pare(in-phase)and/or couple pair(180 degrees out-of-phase)weights for balancing.The latter approach has been used quite often in the field,especially on steam turbine and generator rotors.Dependent on vibration mode shapes and combinations as well as balance plane accessibility,sometimes applying static or couple pair weights can be a wise choice;other times weights at one or two end planes are needed.There are totally 4 possible sets of IC data due to weights at plane 1,plane 2,static pair,and couple pair.IC data would typically be obtained by applying trial weights followed by trial weight runs.It is found,however,that all these IC data can be converted easily without trial weight runs once any two of 4 sets are known.The above findings and conversion equations have been obtained analytically and verified by experimental results.Real cases are given to demonstratetheirapplications.
基金supported by National Hi-tech Research and Development Program of China (863 Program, Grant No. 2008 AA04Z114)
文摘Measurement error of unbalance's vibration response plays a crucial role in calibration and on-line updating of influence coefficient(IC). Focusing on the two problems that the moment estimator of data used in calibration process cannot fulfill the accuracy requirement under small sample and the disturbance of measurement error cannot be effectively suppressed in updating process, an IC calibration and on-line updating method based on hierarchical Bayesian method for automatic dynamic balancing machine was proposed. During calibration process, for the repeatedly-measured data obtained from experiments with different trial weights, according to the fact that measurement error of each sensor had the same statistical characteristics, the joint posterior distribution model for the true values of the vibration response under all trial weights and measurement error was established. During the updating process, information obtained from calibration was regarded as prior information, which was utilized to update the posterior distribution of IC combined with the real-time reference information to implement online updating. Moreover, Gibbs sampling method of Markov Chain Monte Carlo(MCMC) was adopted to obtain the maximum posterior estimation of parameters to be estimated. On the independent developed dynamic balancing testbed, prediction was carried out for multiple groups of data through the proposed method and the traditional method respectively, the result indicated that estimator of influence coefficient obtained through the proposed method had higher accuracy; the proposed updating method more effectively guaranteed the measurement accuracy during the whole producing process, and meantime more reasonably compromised between the sensitivity of IC change and suppression of randomness of vibration response.
文摘In order to develop a general calculating rotor’s torsional stiffness based on stiffness influence coefficient for different rotor assembling, the calculation method of the torsional stiffness influence coefficient of equal thickness disc is researched in this paper at first. Then the torsional stiffness influence coefficient λ of equal thickness disc is fit to a binary curved face and a calculation equation is obtained based on a large quantity of calculating data, which lays the foundation for research on a general calculating method of rotor torsional stiffness. Thirdly a simplified calculation method for equivalent stiffness diameter of stepped equal thickness disc and cone disc in the steam turbine generators is suggested. Finally a general calculating program for calculating rotor’s torsional vibration features is developed, and the torsional vibration features of a verity of steam turbine rotors are calculated for verification. The calculating results show that stiffness influence coefficient λ of equal-thickness disc depends on parameters of B and H, as well as the stiffness influence coefficient λ;and discs with complex structure can be simplified to equal-thickness discs with little error by using the method suggested in this paper;error can be controlled within 1% when equivalent diameter of stiffness is calculated by this method.
文摘In this paper, we study local influence analysis for Zhang's generalized correlation coefficients and Hotelling's generalized correlation coefficient by using approach. of local influence analysis suggested by Shi (1991), i.e., generalized influence function (GIF) and generalized Cook distance (GCD). An example is given to illustrate our results.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51175422)
文摘In order to execute geometric analysis for planar deployable mechanism of scissors unit,the dynamic analysis model of scissor planar deployable structure is created based on the Cartesian coordinate system,the influence coefficient is acquired by means of the coordinate transformation,combining the D 'Alembert 's principle with Dynamic-Static method,the dynamic characteristic analysis is completed finally. Moreover,specific calculating examples are adopted to verify the effectiveness of proposed method,and the result shows that the movement of each component of scissors unit mechanism is more smooth during initial deployment stage,however,when the configuration angle θ of unit mechanism is approaching π,some comparative large variations would appear on movement parameters and hinge constraint force.
文摘Based on the model structure of the influence coefficient method analyzed in depth by matrix theory ,it is explained the reason why the unreasonable and instable correction masses with bigger MSE are obtained by LS influence coefficient method when there are correlation planes in the dynamic balancing. It also presencd the new ridge regression method for solving correction masses according to the Tikhonov regularization theory, and described the reason why the ridge regression can eliminate the disadvantage of the LS method. Applying this new method to dynamic balancing of gas turbine, it is found that this method is superior to the LS method when influence coefficient matrix is ill-conditioned,the minimal correction masses and residual vibration are obtained in the dynamic balancing of rotors.
基金Project(20050216007) supported by the Specialized Research Fund for the Doctoral Program of Higher Education of China
文摘According to the concept of virtual bending force,a rational explanation for SHOHET's model was presented. Considering the deformation characters of the work rolls in four-high mill,the deformation model of the work roll was regarded as a cantilever beam and new influence coefficients were deduced.The effect of the bending force was taken into account independently. Therefore,the contribution to work roll deflection caused by rolling load,rolling pressure between rolls and bending force can be got from the new formulas.To validate the accuracy of the formulas,the results obtained from the new formulas were compared with those from SHOHET's formulas.It is found that they highly coincide,which illustrates that the formulas are reliable.
基金Supported by the National Key Research and Development Program of China under Grant No 2016YFA0301202the National Natural Science Foundation of China under Grant Nos 11474275,61674135 and 91536101+1 种基金the Strategic Priority Research Program of the Chinese Academy of Sciences under Grant No XDPB0603the China Postdoctoral Science Foundation under Grant No 2017M622400
文摘The influence of the pressure transmission medium(PTM)on the excitonic interband transitions in monolayer tungsten diselenide(WSe2)is investigated using photoluminescence(PL)spectra under hydrostatic pressure up to 5GPa.Three kinds of PTMs,condensed argon(Ar),1:1 n-pentane and isopentane mixture(PM),and4:1 methanol and ethanol mixture(MEM,a PTM with polarity),are used.It is found that when either Ar or PM is used as the PTM,the PL peak of exciton related to the direct K-K interband transition shows a pressure-induced blue-shift at a rate of 32±4 or 32±1 meV/GPa,while it turns to be 50±9meV/GPa when MEM is used as the PTM.The indirect A-K interband transition presents almost no shift with increasing pressure up to approximatel.y 5 GPa when Ar and PM are used as the PTM,while it shows a red-shift at the rate of-17±7meV/GPa by using MEM as the PTM.These results reveal that the optical interband transitions of monolayer WSe2 are very sensitive to the polarity of the PTM.The anomalous pressure coefficient obtained using the polar PTM of MEM is ascribed to the existence of hydrogen-like bonds between hydroxyl in MEM and Se atoms under hydrostatic pressure.
文摘Typical rotors such as those on steam turbine or generator are often supported by two bearings with two balance planes at both ends.Vibrations are monitored by a pair of proximity probes at each bearing.There are two approaches to reduce 1X vibration due to unbalance at both ends via balancing with influence coefficient method.The first approach is to treat it as a multiple-plane balancing problem involving 2x2 matrix of complex influence coefficients.The second approach is to treat it as two single-plane balance problems using static(in-phase)and couple(180 degree out-of-phase)components,respectively.Conversion equations of influence coefficients between these two approaches have been found previously by the author.The corresponding spreadsheets that convert influence coefficients between these two formats are presented in the current paper.The paper shows effectiveness of these conversion equations in dealing with real balancing problems in the field.A detailed balance case is presented to demonstrate how the conversion equations are used to reduce vibration effectively.
文摘The steel slag (SS) permeable concrete was prepared by SS. The influences of the aggregate-cement rate, the aggregate particle size, the water-cement rate, the admixture dosage and other factors on the permeability coefficient of SS permeable concrete were analyzed. The law of influence was also investigated. The study serves as a technological reference for the construction and design of SS permeable concrete.
基金the Beijing Municipal Scienceand Technology Project (No.KM202111417006)the Academic Research Projects of Beijing Union University (Nos.ZK10202305 and ZK80202004)the Beijing Municipal Science and Technology Project (No.KM202111417005)。
文摘Due to the diversity of work requirements and environment,the number of degrees of freedom(DOFs)and the complexity of structure of industrial robots are constantly increasing.It is difficult to establish the accurate dynamical model of industrial robots,which greatly hinders the realization of a stable,fast and accurate trajectory tracking control.Therefore,the general expression of the constraint relation in the explicit dynamic equation of the multi-DOF industrial robot is derived on the basis of solving the Jacobian matrix and Hessian matrix by using the kinematic influence coefficients method.Moreover,an explicit dynamic equation with general constraint relation expression is established based on the Udwadia-Kalaba theory.The problem of increasing the time of establishing constraint relationship when the multi-DOF industrial robots complete complex task constraints is solved.With the SCARA robot as the research object,the simulation results show that the proposed method can provide a new idea for industrial robot system modeling with complex constraints.
文摘Varying-coefficient single-index model( VCSIM) avoids the so-called "curse of dimensionality " and is flexible enough to include several important statistical models. This paper considers statistical diagnosis for VCSIM. First,the parametric estimation equation is established based on empirical likelihood. Then,some diagnosis statistics are defined. At last, an example is given to illustrate all the results.
基金supported by the National Natural Science Foundation of China (Grant No. 50379046)the Doctoral Fund of the Ministry of Education of China (Grant No. A50221)
文摘Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering applications. A three-dimensional anisotropic heterogeneous steady random seepage model was developed on the basis of the finite element method. A statistical analysis of the distribution characteristics of soil parameters sampled from the main embankment of the Yangtze River in the Southern Jingzhou zone of China was conducted. The Kolomogorov-Smirnov test verified the statistical hypothesis that the permeability coefficient tensor has a Gaussian distribution. With the help of numerical analysis of the stochastic seepage field using the developed model, various statistical and random characteristics of the stochastic seepage field of the main embankment of the Yangtze River in the Southern Jingzhou zone of China were investigated. The model was also examined with statistical testing. Through the introduction of random variation of the upstream and downstream water levels into the model, the effects of the boundary randomness due to variation of the downstream and upstream water levels on the variation of simulated results presented with a vector series of the random seepage field were analyzed. Furthermore, the combined influence of the variation of the soil permeability coefficient and such seepage resistance measures as the cut-off wall and relief ditch on the hydraulic head distribution was analyzed and compared with the results obtained by determinate analysis. Meanwhile, sensitivities of the hydraulic gradient and downstream exit height to the variation of boundary water level were studied. The validity of the simulated results was verified by stochastic testing and measured data. The developed model provides more detail and a full stochastic algorithm to characterize and analyze three-dimensional stochastic seepage field problems.
基金This project was financially supported by the National Committee of Science and Technology Grants/903-85-08-05
文摘A three-dimensional, first order turbulence closure, thermal diffusion model is described in this paper. The governing equations consist of an equation of continuity, three components of momentum, conservation equations for salt, temperature and subgridscale energy, and an equation of state. In the model, according to the hypothesis of Kolmogorov and Prandtl, the viscosity coefficient of turbulent flow of homogeneous fluid is related to the local turbulent energy, and the horizontal and vertical exchange coefficients of mass, heat and momentum are computed with the introduction of subgridscale turbulence energy. The governing equations are solved by finite difference techniques. This model is applied to the Jiaozhou bay to predict thermal pollution by the Huangdao power plant. An instantaneous tidal current field is computed, then the distribution of temperature increment is predicted, and finally the effect of wind stress on thermal discharge is discussed.
文摘In this paper, a varying-coefficient density-ratio model for case-control studies is developed. We investigate the local empirical likelihood diagnosis of varying coefficient density-ratio model for case-control data. The local empirical log-likelihood ratios for the nonparametric coefficient functions are introduced. First, the estimation equations based on empirical likelihood method are established. Then, a few of diagnostic statistics are proposed. At last, we also examine the performance of proposed method for finite sample sizes through simulation studies.