The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female populat...The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female population to replace the traditional standard that utilizes simplified breast models,for the accurate estimation of the mean glandular dose of a patient undergoing digital breast tomosynthesis(DBT).The first set of detailed Chinese female breast models and representative breast parameters was constructed.Considering backscatter radiation and computational efficiency,we improved the combination of these models and the Chinese reference adult female whole-body voxel phantom.Image acquisition for four commercial DBT systems that are widely employed in China was simulated using the Monte Carlo method to obtain the normalized glandular dose coefficients of DBT(D_(gN)^(DBT))and the glandular depth dose(D_(g)^(dep)(z))for different breast characteristics and X-ray spectra.We calculated a series of D_(gN)^(DBT) values for breasts with different percentage mass glandularities(5%,25%,50%,75%,and 100%)and compressed breast thicknesses(2,3,4,5,6,and 7 cm)at various tube potentials(25,28,30,32,35,and 49 kV)and target/filter combinations(W/Rh,W/Al,Mo/Mo,Rh/Rh,and Rh/Ag).The parameter dependence of the breast characteristics and beam conditions on D_(gN)^(DBT) in detailed breast models was investigated.The D_(gN)^(DBT) results were 14.6-51.0%lower than those of the traditional dosimetry standard in China.The difference in D_(gN)^(DBT) was mainly due to a decrease in the depth of the main energy deposition area caused by the glandular distribution along the depth direction.The results obtained in this study may be used to improve breast dosimetry in China and provide more detailed information on risk assessment during DBT.展开更多
This study assesses the capability of a coarse-resolution ocean model to replicate the response of the Southern Ocean Meridional Overturning Circulation(MOC) to intensified westerlies,focusing on the role of the eddy ...This study assesses the capability of a coarse-resolution ocean model to replicate the response of the Southern Ocean Meridional Overturning Circulation(MOC) to intensified westerlies,focusing on the role of the eddy transfer coefficient(κ).κ is a parameter commonly used to represent the velocities induced by unresolved eddies.Our findings reveal that a stratification-dependent κ,incorporating spatiotemporal variability,leads to the most robust eddy-induced MOC response,capturing 82% of the reference eddy-resolving simulation.Decomposing the eddy-induced velocity into its vertical variation(VV) and spatial structure(SS) components unveils that the enhanced eddy compensation response primarily stems from an augmented SS term,while the introduced VV term weakens the response.Furthermore,the temporal variability of the stratification-dependent κ emerges as a key factor in enhancing the eddy compensation response to intensified westerlies.The experiment with stratification-dependent κ exhibits a more potent eddy compensation response compared to the constant κ,attributed to the structure of κ and the vertical variation of the density slope.These results underscore the critical role of accurately representing κ in capturing the response of the Southern Ocean MOC and emphasize the significance of the isopycnal slope in modulating the eddy compensation mechanism.展开更多
Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using general...Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using generalized linear models in fixed effects/coefficients. Correlations are modeled using random effects/coefficients. Nonlinearity is addressed using power transforms of primary (untransformed) predictors. Parameter estimation is based on extended linear mixed modeling generalizing both generalized estimating equations and linear mixed modeling. Models are evaluated using likelihood cross-validation (LCV) scores and are generated adaptively using a heuristic search controlled by LCV scores. Cases covered include linear, Poisson, logistic, exponential, and discrete regression of correlated continuous, count/rate, dichotomous, positive continuous, and discrete numeric outcomes treated as normally, Poisson, Bernoulli, exponentially, and discrete numerically distributed, respectively. Example analyses are also generated for these five cases to compare adaptive random effects/coefficients modeling of correlated outcomes to previously developed adaptive modeling based on directly specified covariance structures. Adaptive random effects/coefficients modeling substantially outperforms direct covariance modeling in the linear, exponential, and discrete regression example analyses. It generates equivalent results in the logistic regression example analyses and it is substantially outperformed in the Poisson regression case. Random effects/coefficients modeling of correlated outcomes can provide substantial improvements in model selection compared to directly specified covariance modeling. However, directly specified covariance modeling can generate competitive or substantially better results in some cases while usually requiring less computation time.展开更多
This work is a simulation model with the LAMMPS calculation code of an electrode based on alkali metal oxides (lithium, sodium and potassium) using the Lennard Jones potential. For a multiplicity of 8*8*8, we studied ...This work is a simulation model with the LAMMPS calculation code of an electrode based on alkali metal oxides (lithium, sodium and potassium) using the Lennard Jones potential. For a multiplicity of 8*8*8, we studied a gap-free model using molecular dynamics. Physical quantities such as volume and pressure of the Na-O and Li-O systems exhibit similar behaviors around the thermodynamic ensembles NPT and NVE. However, for the Na2O system, at a minimum temperature value, we observe a range of total energy values;in contrast, for the Li2O system, a minimum energy corresponds to a range of temperatures. Finally, for physicochemical properties, we studied the diffusion coefficient and activation energy of lithium and potassium oxides around their melting temperatures. The order of magnitude of the diffusion coefficients is given by the relation Dli-O >DNa-O for the multiplicity 8*8*8, while for the activation energy, the order is well reversed EaNa-O > EaLi-O.展开更多
The populations of urban centers in Congo-Brazzaville have decided to develop various methods of water storage (concrete or masonry underground tanks) for domestic use, due to shortages in the distribution of water th...The populations of urban centers in Congo-Brazzaville have decided to develop various methods of water storage (concrete or masonry underground tanks) for domestic use, due to shortages in the distribution of water through the public network. However, questions remain as to the physico-chemical quality of the water stored in these tanks, when these structures are built in wet and relatively polluted areas. This paper presents a model of pollutant diffusion through the cementitious matrix (concrete) of tank walls simulated at a buried reservoir. The results of the experimental and numerical simulations show that certain concrete parameters, such as porosity, permeability and diffusivity, have a significant influence on the transfer of pollutants through the concrete walls, thus altering the physico-chemical quality of the stored water. The numerical models (1D) used to predict pollutant transfer and the quality of the stored water are consistent with those of the optimal control for identifying the diffusion coefficient. Major ion concentrations appear to be correlated with system porosity and diffusion coefficient. Nevertheless, the identification of the diffusion coefficient from the optimal control method, based on an explicit numerical resolution of a finite volume PDE for the approximation of the experiment, is not consistent with that of the optimal control method.展开更多
Remaining useful life(RUL) prediction is one of the most crucial elements in prognostics and health management(PHM). Aiming at the imperfect prior information, this paper proposes an RUL prediction method based on a n...Remaining useful life(RUL) prediction is one of the most crucial elements in prognostics and health management(PHM). Aiming at the imperfect prior information, this paper proposes an RUL prediction method based on a nonlinear random coefficient regression(RCR) model with fusing failure time data.Firstly, some interesting natures of parameters estimation based on the nonlinear RCR model are given. Based on these natures,the failure time data can be fused as the prior information reasonably. Specifically, the fixed parameters are calculated by the field degradation data of the evaluated equipment and the prior information of random coefficient is estimated with fusing the failure time data of congeneric equipment. Then, the prior information of the random coefficient is updated online under the Bayesian framework, the probability density function(PDF) of the RUL with considering the limitation of the failure threshold is performed. Finally, two case studies are used for experimental verification. Compared with the traditional Bayesian method, the proposed method can effectively reduce the influence of imperfect prior information and improve the accuracy of RUL prediction.展开更多
This study is an extension of the previous work done with ARS-680 Environmental Chamber. Drying is a complex operation that demands much energy and time. Drying is essentially important for preservation of ginger rhiz...This study is an extension of the previous work done with ARS-680 Environmental Chamber. Drying is a complex operation that demands much energy and time. Drying is essentially important for preservation of ginger rhizome. Drying of ginger was modeled, and then the effective diffusion coefficient and activation energy were determined. For this purpose, the experiments were done at six levels of varied temperatures: 10°C, 20°C, 30°C, 40°C, 50°C and 60°C. The values of effective diffusion coefficients obtained in this work for the variously treated ginger rhizomes closely agreed with the average effective diffusion coefficients of other notable authors who determined the drying kinetics and convective heat transfer coefficients of ginger slices.展开更多
In this study,the effects of surface exchange coefficients on simulations of Super Typhoon Megi(2010)are investigated using a fully coupled ocean-atmosphere-wave model.Several experiments are conducted using different...In this study,the effects of surface exchange coefficients on simulations of Super Typhoon Megi(2010)are investigated using a fully coupled ocean-atmosphere-wave model.Several experiments are conducted using different parameterization schemes for the drag(C_(D))and enthalpy exchange(C_(K))coefficients.For the selected case,considering only the leveling-off of C_(D)at high wind speeds does not effectively improve the simulated typhoon track,intensity,or size.We found that increasing C_(K)monotonically with wind speed(Komori et al.,2018)yields stronger winds and deeper pressures by enhancing latent and sensible heat fluxes,but typhoon intensity remains underestimated.We propose a new higher C_(K)than that from Komori et al.(2018)based on the theory of Emanuel(1995).This approach produces a greater modeled typhoon intensity that is in good agreement with the best track data and effectively improves the track error for the simulation.Improved accuracy for modeled typhoon intensity is achieved with the new coefficient because C_(K)/C_(D)reaches the threshold of about 0.75 predicted by Emanuel(1995).The new proposed C_(K)also results in a reasonably accurate modeled sea surface temperature.However,typhoon size and surface wave height are overestimated.This finding implies that more numerical tests for tropical cyclones of different nature(such as strong,weak,dissipating,rapidly intensifying,or weakening tropical cyclones)should be studied,and more physical processes should be explored in future coupled models.展开更多
It is common for datasets to contain both categorical and continuous variables. However, many feature screening methods designed for high-dimensional classification assume that the variables are continuous. This limit...It is common for datasets to contain both categorical and continuous variables. However, many feature screening methods designed for high-dimensional classification assume that the variables are continuous. This limits the applicability of existing methods in handling this complex scenario. To address this issue, we propose a model-free feature screening approach for ultra-high-dimensional multi-classification that can handle both categorical and continuous variables. Our proposed feature screening method utilizes the Maximal Information Coefficient to assess the predictive power of the variables. By satisfying certain regularity conditions, we have proven that our screening procedure possesses the sure screening property and ranking consistency properties. To validate the effectiveness of our approach, we conduct simulation studies and provide real data analysis examples to demonstrate its performance in finite samples. In summary, our proposed method offers a solution for effectively screening features in ultra-high-dimensional datasets with a mixture of categorical and continuous covariates.展开更多
The traditional model selection criterions try to make a balance between fitted error and model complexity. Assumptions on the distribution of the response or the noise, which may be misspecified, should be made befor...The traditional model selection criterions try to make a balance between fitted error and model complexity. Assumptions on the distribution of the response or the noise, which may be misspecified, should be made before using the traditional ones. In this ar- ticle, we give a new model selection criterion, based on the assumption that noise term in the model is independent with explanatory variables, of minimizing the association strength between regression residuals and the response, with fewer assumptions. Maximal Information Coe^cient (MIC), a recently proposed dependence measure, captures a wide range of associ- ations, and gives almost the same score to different type of relationships with equal noise, so MIC is used to measure the association strength. Furthermore, partial maximal information coefficient (PMIC) is introduced to capture the association between two variables removing a third controlling random variable. In addition, the definition of general partial relationship is given.展开更多
The research of groundwater vulnerability is the basic work to protect the groundwater. For utilizing groundwater resource continuably, groundwater vulnerability evaluation is necessary. Useful reference to protect, e...The research of groundwater vulnerability is the basic work to protect the groundwater. For utilizing groundwater resource continuably, groundwater vulnerability evaluation is necessary. Useful reference to protect, exploit and utilize on groundwater resource are provided rationally. According to the real condition of Sanjiang Plain, the indexes system is established based on the traditional DRASTIC model. The new system includes the following seven indexes: Depth of Water, Net Recharge, Aquifer Media, Soil Media, Conductivity of the Aquifer, Land Utilizing Ratio and Populace Density. The related analysis appears that the system is rather reasonable. Because traditional methods, such as analytic hierarchy process and fuzzy mathematics theory, can't be avoided human interference in selection of weights, they could lead to an imprecise result. In order to evaluate the groundwater vulnerability reasonably, entropy weight coefficient method is applied for the first time, which provides a new way to groundwater vulnerability evaluation. The method is a model whose weights are insured by the calculation process, so the artificial disturb can be avoided. It has been used to evaluate the groundwater vulnerability in Sanjiang Plain. The satisfied result is acquired. Comparably, the same result is acquired by the other method named projection pursuit evaluation based on real-coded accelerating genetic algorithm. It shows that entropy weight coefficient method is applicable on groundwater vulnerability evaluation. The evaluation result can provide reference on the decision-making departments.展开更多
According to the experimental data of the orifice discharge coefficient for the flow through a vertical sharp-edged orifice obtained in the previous study of this work,a theoretical criterion for flow mechanisms of sm...According to the experimental data of the orifice discharge coefficient for the flow through a vertical sharp-edged orifice obtained in the previous study of this work,a theoretical criterion for flow mechanisms of small orifice(viz.thick-walled orifice and nozzle) and large orifice(viz.thin-walled orifice) was proposed based on the ratio of orifice diameter to plate thickness.It can help explain the dissipation of the mechanical energy loss in the flow process for the two flow mechanisms under different operating regimes.The main parameters such as orifice diameter,plate thickness and liquid head were correlated,and a semi-empirical model for orifice coefficient and an empirical model with high precision at the stable region were developed.展开更多
Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient (Ca) and oceanic drag coefficient (Cw). Parameterizations of drag coefficients rather than constant values provide us a ...Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient (Ca) and oceanic drag coefficient (Cw). Parameterizations of drag coefficients rather than constant values provide us a way to look insight into the dependence of these characteristics on sea ice conditions. In the present study, the parameterized ice drag coefficients are included into a free-drift sea ice dynamic model, and the wind factor a and the deflection angle θ between sea ice drift and wind velocity as well as the ratio of Ca to Cw are studied to investigate their dependence on the impact factors such as local drag coefficients, floe and ridge geometry. The results reveal that in an idealized steady ocean, Ca/Cw increases obviously with the increasing ice concentration for small ice floes in the marginal ice zone, while it remains at a steady level (0.2-0.25) for large floes in the central ice zone. The wind factor a increases rapidly at first and approaches a steady level of 0.018 when A is greater than 20%. And the deflection angle ~ drops rapidly from an initial value of approximate 80° and decreases slowly as A is greater than 20% without a steady level like a. The values of these parameters agree well with the previously reported observations in Arctic. The ridging intensity is an important parameter to determine the dominant contribution of the ratio of skin friction drag coefficient (Cs'/Cs) and the ratio of ridge form drag coefficient (Cr'/Cr) to the value of Ca/Cw, a, and 8, because of the dominance of ridge form drag for large ridging intensity and skin friction for small ridging intensity among the total drag forces. Parameterization of sea ice drag coefficients has the potential to be embedded into ice dynamic models to better account for the variability of sea ice in the transient Arctic Ocean.展开更多
This letter presents a theoretical model of the normal (head-on) collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. Aft...This letter presents a theoretical model of the normal (head-on) collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. After the contact force law between the contacted spheres during a collision is phenomenologically formulated in terms of the compression or overlap displacement under considera- tion of an elastic-plastic loading and a plastic unloading subprocesses, the coefficient of restitution is gained by the dynamic equation of the contact process once an initial impact velocity is input. It is found that the theoretical predictions of the coefficient of restitution varying with the impact velocity are well in agreement with the existing experimental characteristics which are fitted by the explicit formula.展开更多
Based on medium plate runout table ultra-fast cooling( UFC)-accelerated cooling equipment( ACC) system,a heat transfer coefficient model was constructed. Firstly,according to the measured data,heat transfer coefficien...Based on medium plate runout table ultra-fast cooling( UFC)-accelerated cooling equipment( ACC) system,a heat transfer coefficient model was constructed. Firstly,according to the measured data,heat transfer coefficients under different roll speed and water volume were calculated by using an inverse heat conduction method. Secondly,a monofactorial heat transfer coefficient calculation formula was obtained. Finally,the heat transfer coefficient model based on medium plate runout table UFC-ACC system was constructed by intercept function,slope function,interaction influence function and linear or nonlinear influencing factors. The precision of these models was validated by comparing model prediction value with measured data,and the results were in good agreement with practical needs,and the average deviation was less than 5%.展开更多
The modified sub regular solution model was used for a calculation of the activity coefficient of immiscible binary alloy systems. The parameters needed for the calculation are the interaction parameters, λ 1 a...The modified sub regular solution model was used for a calculation of the activity coefficient of immiscible binary alloy systems. The parameters needed for the calculation are the interaction parameters, λ 1 and λ 2, which are represented as a linear function of temperature, T . The molar excess Gibbs free energy, G m E, can be written in the form G m E= x A x B[( λ 11 + λ 12 T )+( λ 21 + λ 22 T ) x B ] The calculation is carried out numerically for three immiscible binary alloy systems, Al Pb, Cu Tl and In V. The agreement between the calculated and experimentally determined values of activity coefficient is excellent.展开更多
Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivi...Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivity index, relying on the meso-scale structure of porous material and the property of liquid, denotes the connectivity of pores in Representative Element Area (REA). If the conductivity of the porous material is anisotropic, the equivalent connectivity index is a second order tensor. Based on the basic theories of continuous mechanics and tensor analysis, relationship between area porosity and volumetric porosity of porous materials is deduced. Then a generalized expression, describing the relation between effective stress coefficient tensor and equivalent connectivity tensor of pores, is proposed, and the expression can be applied to isotropic media and also to anisotropic materials. Furthermore, evolution of porosity and equivalent connectivity index of the pore are studied in the strain space, and the method to determine the corresponding functions in expressions above is proposed using genetic algorithm and genetic programming. Two applications show that the results obtained by the method in this paper perfectly agree with the test data. This paper provides an important theoretical support to the coupled hydro-mechanical research.展开更多
A real-time channel flood forecast model was developed to simulate channel flow in plain rivers based on the dynamic wave theory. Taking into consideration channel shape differences along the channel, a roughness upda...A real-time channel flood forecast model was developed to simulate channel flow in plain rivers based on the dynamic wave theory. Taking into consideration channel shape differences along the channel, a roughness updating technique was developed using the Kalman filter method to update Manning's roughness coefficient at each time step of the calculation processes. Channel shapes were simplified as rectangles, triangles, and parabolas, and the relationships between hydraulic radius and water depth were developed for plain rivers. Based on the relationship between the Froude number and the inertia terms of the momentum equation in the Saint-Venant equations, the relationship between Manning's roughness coefficient and water depth was obtained. Using the channel of the Huaihe River from Wangjiaba to Lutaizi stations as a case, to test the performance and rationality of the present flood routing model, the original hydraulic model was compared with the developed model. Results show that the stage hydrographs calculated by the developed flood routing model with the updated Manning's roughness coefficient have a good agreement with the observed stage hydrographs. This model performs better than the original hydraulic model.展开更多
The drag coefficient is important in meteorological studies of the boundary layer because it describes the air-sea momentum flux. Eight drag coefficient schemes were assessed. These parametrizations were compared taki...The drag coefficient is important in meteorological studies of the boundary layer because it describes the air-sea momentum flux. Eight drag coefficient schemes were assessed. These parametrizations were compared taking into account data from in situ and laboratory observations.The drag coefficients determined using three schemes were consistent with the level-off phenomenon, supported by the results of laboratory studies. The drag coefficient determined using one scheme decreased at wind speeds higher than approximately 30 m s-1, in agreement with indirect measurements under typhoon conditions. In contrast, the drag coefficients determined using the other four schemes increased with wind speed, even under high wind regimes. Sensitivity tests were performed using simulations of two super typhoons in the Weather Research and Forecasting model. While the typhoon tracks were negligibly sensitive to the parametrization used, the typhoon intensities (the maximum lO-m wind speed and the minimum sea level pressure), sizes, and structure, were very sensitive to it.展开更多
基金supported by the National Natural Science Foundation of China(Nos.U2167209 and 12175114)the National Key R&D Program of China(No.2021YFF0603600).
文摘The rise in breast cancer diagnoses among Chinese women has necessitated the use of X-ray breast screening,which carries a radiation risk.This study aimed to provide a dosimetry protocol for the Chinese female population to replace the traditional standard that utilizes simplified breast models,for the accurate estimation of the mean glandular dose of a patient undergoing digital breast tomosynthesis(DBT).The first set of detailed Chinese female breast models and representative breast parameters was constructed.Considering backscatter radiation and computational efficiency,we improved the combination of these models and the Chinese reference adult female whole-body voxel phantom.Image acquisition for four commercial DBT systems that are widely employed in China was simulated using the Monte Carlo method to obtain the normalized glandular dose coefficients of DBT(D_(gN)^(DBT))and the glandular depth dose(D_(g)^(dep)(z))for different breast characteristics and X-ray spectra.We calculated a series of D_(gN)^(DBT) values for breasts with different percentage mass glandularities(5%,25%,50%,75%,and 100%)and compressed breast thicknesses(2,3,4,5,6,and 7 cm)at various tube potentials(25,28,30,32,35,and 49 kV)and target/filter combinations(W/Rh,W/Al,Mo/Mo,Rh/Rh,and Rh/Ag).The parameter dependence of the breast characteristics and beam conditions on D_(gN)^(DBT) in detailed breast models was investigated.The D_(gN)^(DBT) results were 14.6-51.0%lower than those of the traditional dosimetry standard in China.The difference in D_(gN)^(DBT) was mainly due to a decrease in the depth of the main energy deposition area caused by the glandular distribution along the depth direction.The results obtained in this study may be used to improve breast dosimetry in China and provide more detailed information on risk assessment during DBT.
基金supported by the National Key R&D Program for Developing Basic Sciences(2022YFC3104802)the National Natural Science Foundation of China(Nos.42306219 and 42106020)+3 种基金the Tai Shan Scholar Pro-gram(Grant No.tstp20231237)Part of computing resources are financially supported by Laoshan Laboratory(No.LSKJ202300301)Dr.Eric P.CHASSIGNET is supported by the CAS President’s International Fellowship Initiative(PIFI)NOAA Climate Program Office MAPP Program(Award NA15OAR4310088).
文摘This study assesses the capability of a coarse-resolution ocean model to replicate the response of the Southern Ocean Meridional Overturning Circulation(MOC) to intensified westerlies,focusing on the role of the eddy transfer coefficient(κ).κ is a parameter commonly used to represent the velocities induced by unresolved eddies.Our findings reveal that a stratification-dependent κ,incorporating spatiotemporal variability,leads to the most robust eddy-induced MOC response,capturing 82% of the reference eddy-resolving simulation.Decomposing the eddy-induced velocity into its vertical variation(VV) and spatial structure(SS) components unveils that the enhanced eddy compensation response primarily stems from an augmented SS term,while the introduced VV term weakens the response.Furthermore,the temporal variability of the stratification-dependent κ emerges as a key factor in enhancing the eddy compensation response to intensified westerlies.The experiment with stratification-dependent κ exhibits a more potent eddy compensation response compared to the constant κ,attributed to the structure of κ and the vertical variation of the density slope.These results underscore the critical role of accurately representing κ in capturing the response of the Southern Ocean MOC and emphasize the significance of the isopycnal slope in modulating the eddy compensation mechanism.
文摘Adaptive fractional polynomial modeling of general correlated outcomes is formulated to address nonlinearity in means, variances/dispersions, and correlations. Means and variances/dispersions are modeled using generalized linear models in fixed effects/coefficients. Correlations are modeled using random effects/coefficients. Nonlinearity is addressed using power transforms of primary (untransformed) predictors. Parameter estimation is based on extended linear mixed modeling generalizing both generalized estimating equations and linear mixed modeling. Models are evaluated using likelihood cross-validation (LCV) scores and are generated adaptively using a heuristic search controlled by LCV scores. Cases covered include linear, Poisson, logistic, exponential, and discrete regression of correlated continuous, count/rate, dichotomous, positive continuous, and discrete numeric outcomes treated as normally, Poisson, Bernoulli, exponentially, and discrete numerically distributed, respectively. Example analyses are also generated for these five cases to compare adaptive random effects/coefficients modeling of correlated outcomes to previously developed adaptive modeling based on directly specified covariance structures. Adaptive random effects/coefficients modeling substantially outperforms direct covariance modeling in the linear, exponential, and discrete regression example analyses. It generates equivalent results in the logistic regression example analyses and it is substantially outperformed in the Poisson regression case. Random effects/coefficients modeling of correlated outcomes can provide substantial improvements in model selection compared to directly specified covariance modeling. However, directly specified covariance modeling can generate competitive or substantially better results in some cases while usually requiring less computation time.
文摘This work is a simulation model with the LAMMPS calculation code of an electrode based on alkali metal oxides (lithium, sodium and potassium) using the Lennard Jones potential. For a multiplicity of 8*8*8, we studied a gap-free model using molecular dynamics. Physical quantities such as volume and pressure of the Na-O and Li-O systems exhibit similar behaviors around the thermodynamic ensembles NPT and NVE. However, for the Na2O system, at a minimum temperature value, we observe a range of total energy values;in contrast, for the Li2O system, a minimum energy corresponds to a range of temperatures. Finally, for physicochemical properties, we studied the diffusion coefficient and activation energy of lithium and potassium oxides around their melting temperatures. The order of magnitude of the diffusion coefficients is given by the relation Dli-O >DNa-O for the multiplicity 8*8*8, while for the activation energy, the order is well reversed EaNa-O > EaLi-O.
文摘The populations of urban centers in Congo-Brazzaville have decided to develop various methods of water storage (concrete or masonry underground tanks) for domestic use, due to shortages in the distribution of water through the public network. However, questions remain as to the physico-chemical quality of the water stored in these tanks, when these structures are built in wet and relatively polluted areas. This paper presents a model of pollutant diffusion through the cementitious matrix (concrete) of tank walls simulated at a buried reservoir. The results of the experimental and numerical simulations show that certain concrete parameters, such as porosity, permeability and diffusivity, have a significant influence on the transfer of pollutants through the concrete walls, thus altering the physico-chemical quality of the stored water. The numerical models (1D) used to predict pollutant transfer and the quality of the stored water are consistent with those of the optimal control for identifying the diffusion coefficient. Major ion concentrations appear to be correlated with system porosity and diffusion coefficient. Nevertheless, the identification of the diffusion coefficient from the optimal control method, based on an explicit numerical resolution of a finite volume PDE for the approximation of the experiment, is not consistent with that of the optimal control method.
基金supported by National Natural Science Foundation of China (61703410,61873175,62073336,61873273,61773386,61922089)。
文摘Remaining useful life(RUL) prediction is one of the most crucial elements in prognostics and health management(PHM). Aiming at the imperfect prior information, this paper proposes an RUL prediction method based on a nonlinear random coefficient regression(RCR) model with fusing failure time data.Firstly, some interesting natures of parameters estimation based on the nonlinear RCR model are given. Based on these natures,the failure time data can be fused as the prior information reasonably. Specifically, the fixed parameters are calculated by the field degradation data of the evaluated equipment and the prior information of random coefficient is estimated with fusing the failure time data of congeneric equipment. Then, the prior information of the random coefficient is updated online under the Bayesian framework, the probability density function(PDF) of the RUL with considering the limitation of the failure threshold is performed. Finally, two case studies are used for experimental verification. Compared with the traditional Bayesian method, the proposed method can effectively reduce the influence of imperfect prior information and improve the accuracy of RUL prediction.
文摘This study is an extension of the previous work done with ARS-680 Environmental Chamber. Drying is a complex operation that demands much energy and time. Drying is essentially important for preservation of ginger rhizome. Drying of ginger was modeled, and then the effective diffusion coefficient and activation energy were determined. For this purpose, the experiments were done at six levels of varied temperatures: 10°C, 20°C, 30°C, 40°C, 50°C and 60°C. The values of effective diffusion coefficients obtained in this work for the variously treated ginger rhizomes closely agreed with the average effective diffusion coefficients of other notable authors who determined the drying kinetics and convective heat transfer coefficients of ginger slices.
基金supported by the National Natural Science Foundation of China(Nos.41906014,U20A2099 and 41976017)。
文摘In this study,the effects of surface exchange coefficients on simulations of Super Typhoon Megi(2010)are investigated using a fully coupled ocean-atmosphere-wave model.Several experiments are conducted using different parameterization schemes for the drag(C_(D))and enthalpy exchange(C_(K))coefficients.For the selected case,considering only the leveling-off of C_(D)at high wind speeds does not effectively improve the simulated typhoon track,intensity,or size.We found that increasing C_(K)monotonically with wind speed(Komori et al.,2018)yields stronger winds and deeper pressures by enhancing latent and sensible heat fluxes,but typhoon intensity remains underestimated.We propose a new higher C_(K)than that from Komori et al.(2018)based on the theory of Emanuel(1995).This approach produces a greater modeled typhoon intensity that is in good agreement with the best track data and effectively improves the track error for the simulation.Improved accuracy for modeled typhoon intensity is achieved with the new coefficient because C_(K)/C_(D)reaches the threshold of about 0.75 predicted by Emanuel(1995).The new proposed C_(K)also results in a reasonably accurate modeled sea surface temperature.However,typhoon size and surface wave height are overestimated.This finding implies that more numerical tests for tropical cyclones of different nature(such as strong,weak,dissipating,rapidly intensifying,or weakening tropical cyclones)should be studied,and more physical processes should be explored in future coupled models.
文摘It is common for datasets to contain both categorical and continuous variables. However, many feature screening methods designed for high-dimensional classification assume that the variables are continuous. This limits the applicability of existing methods in handling this complex scenario. To address this issue, we propose a model-free feature screening approach for ultra-high-dimensional multi-classification that can handle both categorical and continuous variables. Our proposed feature screening method utilizes the Maximal Information Coefficient to assess the predictive power of the variables. By satisfying certain regularity conditions, we have proven that our screening procedure possesses the sure screening property and ranking consistency properties. To validate the effectiveness of our approach, we conduct simulation studies and provide real data analysis examples to demonstrate its performance in finite samples. In summary, our proposed method offers a solution for effectively screening features in ultra-high-dimensional datasets with a mixture of categorical and continuous covariates.
基金partly supported by National Basic Research Program of China(973 Program,2011CB707802,2013CB910200)National Science Foundation of China(11201466)
文摘The traditional model selection criterions try to make a balance between fitted error and model complexity. Assumptions on the distribution of the response or the noise, which may be misspecified, should be made before using the traditional ones. In this ar- ticle, we give a new model selection criterion, based on the assumption that noise term in the model is independent with explanatory variables, of minimizing the association strength between regression residuals and the response, with fewer assumptions. Maximal Information Coe^cient (MIC), a recently proposed dependence measure, captures a wide range of associ- ations, and gives almost the same score to different type of relationships with equal noise, so MIC is used to measure the association strength. Furthermore, partial maximal information coefficient (PMIC) is introduced to capture the association between two variables removing a third controlling random variable. In addition, the definition of general partial relationship is given.
基金Supported by the National Natural Science Foundation of China(30400275)the Tackle Key Problems of Heilongjiang Province(the Hobbledehoy Science Fund of Heilongjiang Province)(QC04C28)
文摘The research of groundwater vulnerability is the basic work to protect the groundwater. For utilizing groundwater resource continuably, groundwater vulnerability evaluation is necessary. Useful reference to protect, exploit and utilize on groundwater resource are provided rationally. According to the real condition of Sanjiang Plain, the indexes system is established based on the traditional DRASTIC model. The new system includes the following seven indexes: Depth of Water, Net Recharge, Aquifer Media, Soil Media, Conductivity of the Aquifer, Land Utilizing Ratio and Populace Density. The related analysis appears that the system is rather reasonable. Because traditional methods, such as analytic hierarchy process and fuzzy mathematics theory, can't be avoided human interference in selection of weights, they could lead to an imprecise result. In order to evaluate the groundwater vulnerability reasonably, entropy weight coefficient method is applied for the first time, which provides a new way to groundwater vulnerability evaluation. The method is a model whose weights are insured by the calculation process, so the artificial disturb can be avoided. It has been used to evaluate the groundwater vulnerability in Sanjiang Plain. The satisfied result is acquired. Comparably, the same result is acquired by the other method named projection pursuit evaluation based on real-coded accelerating genetic algorithm. It shows that entropy weight coefficient method is applicable on groundwater vulnerability evaluation. The evaluation result can provide reference on the decision-making departments.
基金supported by the National Natural Science Foundation of China(20806090)
文摘According to the experimental data of the orifice discharge coefficient for the flow through a vertical sharp-edged orifice obtained in the previous study of this work,a theoretical criterion for flow mechanisms of small orifice(viz.thick-walled orifice and nozzle) and large orifice(viz.thin-walled orifice) was proposed based on the ratio of orifice diameter to plate thickness.It can help explain the dissipation of the mechanical energy loss in the flow process for the two flow mechanisms under different operating regimes.The main parameters such as orifice diameter,plate thickness and liquid head were correlated,and a semi-empirical model for orifice coefficient and an empirical model with high precision at the stable region were developed.
基金The National Natural Science Foundation of China under contracts Nos 41276191 and 41306207the Public Science and Technology Research Funds Projects of Ocean under contract No.201205007-05the Global Change Research Program of China under contract No.2015CB953901
文摘Many interesting characteristics of sea ice drift depend on the atmospheric drag coefficient (Ca) and oceanic drag coefficient (Cw). Parameterizations of drag coefficients rather than constant values provide us a way to look insight into the dependence of these characteristics on sea ice conditions. In the present study, the parameterized ice drag coefficients are included into a free-drift sea ice dynamic model, and the wind factor a and the deflection angle θ between sea ice drift and wind velocity as well as the ratio of Ca to Cw are studied to investigate their dependence on the impact factors such as local drag coefficients, floe and ridge geometry. The results reveal that in an idealized steady ocean, Ca/Cw increases obviously with the increasing ice concentration for small ice floes in the marginal ice zone, while it remains at a steady level (0.2-0.25) for large floes in the central ice zone. The wind factor a increases rapidly at first and approaches a steady level of 0.018 when A is greater than 20%. And the deflection angle ~ drops rapidly from an initial value of approximate 80° and decreases slowly as A is greater than 20% without a steady level like a. The values of these parameters agree well with the previously reported observations in Arctic. The ridging intensity is an important parameter to determine the dominant contribution of the ratio of skin friction drag coefficient (Cs'/Cs) and the ratio of ridge form drag coefficient (Cr'/Cr) to the value of Ca/Cw, a, and 8, because of the dominance of ridge form drag for large ridging intensity and skin friction for small ridging intensity among the total drag forces. Parameterization of sea ice drag coefficients has the potential to be embedded into ice dynamic models to better account for the variability of sea ice in the transient Arctic Ocean.
基金supported by the Innovation Team Fund of the National Natural Science Foundation of China(11121202)
文摘This letter presents a theoretical model of the normal (head-on) collisions between two soft spheres for predicting the experimental characteristic of the coefficient of restitution dependent on impact velocity. After the contact force law between the contacted spheres during a collision is phenomenologically formulated in terms of the compression or overlap displacement under considera- tion of an elastic-plastic loading and a plastic unloading subprocesses, the coefficient of restitution is gained by the dynamic equation of the contact process once an initial impact velocity is input. It is found that the theoretical predictions of the coefficient of restitution varying with the impact velocity are well in agreement with the existing experimental characteristics which are fitted by the explicit formula.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51104045)
文摘Based on medium plate runout table ultra-fast cooling( UFC)-accelerated cooling equipment( ACC) system,a heat transfer coefficient model was constructed. Firstly,according to the measured data,heat transfer coefficients under different roll speed and water volume were calculated by using an inverse heat conduction method. Secondly,a monofactorial heat transfer coefficient calculation formula was obtained. Finally,the heat transfer coefficient model based on medium plate runout table UFC-ACC system was constructed by intercept function,slope function,interaction influence function and linear or nonlinear influencing factors. The precision of these models was validated by comparing model prediction value with measured data,and the results were in good agreement with practical needs,and the average deviation was less than 5%.
文摘The modified sub regular solution model was used for a calculation of the activity coefficient of immiscible binary alloy systems. The parameters needed for the calculation are the interaction parameters, λ 1 and λ 2, which are represented as a linear function of temperature, T . The molar excess Gibbs free energy, G m E, can be written in the form G m E= x A x B[( λ 11 + λ 12 T )+( λ 21 + λ 22 T ) x B ] The calculation is carried out numerically for three immiscible binary alloy systems, Al Pb, Cu Tl and In V. The agreement between the calculated and experimentally determined values of activity coefficient is excellent.
基金supported by the Yalongjiang River Joint Fund by the National Natural Science Foundation of China(NSFC)Ertan Hydropower Development Company,LTD(Nos.50579091 and 50539090)+1 种基金NSFC(No.10772190)Major State Basic Research Project of China(No.2002CB412708)
文摘Physical mechanisms and influencing factors on the effective stress coefficient for rock/soil-like porous materials are investigated, based on which equivalent connectivity index is proposed. The equivalent connectivity index, relying on the meso-scale structure of porous material and the property of liquid, denotes the connectivity of pores in Representative Element Area (REA). If the conductivity of the porous material is anisotropic, the equivalent connectivity index is a second order tensor. Based on the basic theories of continuous mechanics and tensor analysis, relationship between area porosity and volumetric porosity of porous materials is deduced. Then a generalized expression, describing the relation between effective stress coefficient tensor and equivalent connectivity tensor of pores, is proposed, and the expression can be applied to isotropic media and also to anisotropic materials. Furthermore, evolution of porosity and equivalent connectivity index of the pore are studied in the strain space, and the method to determine the corresponding functions in expressions above is proposed using genetic algorithm and genetic programming. Two applications show that the results obtained by the method in this paper perfectly agree with the test data. This paper provides an important theoretical support to the coupled hydro-mechanical research.
基金supported by the Special Fund for Public Welfare (Meteorology) of China (Grants No. GYHY201006037 and GYHY200906007)
文摘A real-time channel flood forecast model was developed to simulate channel flow in plain rivers based on the dynamic wave theory. Taking into consideration channel shape differences along the channel, a roughness updating technique was developed using the Kalman filter method to update Manning's roughness coefficient at each time step of the calculation processes. Channel shapes were simplified as rectangles, triangles, and parabolas, and the relationships between hydraulic radius and water depth were developed for plain rivers. Based on the relationship between the Froude number and the inertia terms of the momentum equation in the Saint-Venant equations, the relationship between Manning's roughness coefficient and water depth was obtained. Using the channel of the Huaihe River from Wangjiaba to Lutaizi stations as a case, to test the performance and rationality of the present flood routing model, the original hydraulic model was compared with the developed model. Results show that the stage hydrographs calculated by the developed flood routing model with the updated Manning's roughness coefficient have a good agreement with the observed stage hydrographs. This model performs better than the original hydraulic model.
基金supported by the National Key Basic Research Program of China(973 Program)[grant number 2012CB417402]the Strategic Priority Research Program of the Chinese Academy of Sciences[grant number XDA11010104]+1 种基金the National Natural Science Foundation of China[grant numbers 41576013,41476021,41506023]the National High Technology Research and Development Program of China(863 Program)[grant number2013AA122803]
文摘The drag coefficient is important in meteorological studies of the boundary layer because it describes the air-sea momentum flux. Eight drag coefficient schemes were assessed. These parametrizations were compared taking into account data from in situ and laboratory observations.The drag coefficients determined using three schemes were consistent with the level-off phenomenon, supported by the results of laboratory studies. The drag coefficient determined using one scheme decreased at wind speeds higher than approximately 30 m s-1, in agreement with indirect measurements under typhoon conditions. In contrast, the drag coefficients determined using the other four schemes increased with wind speed, even under high wind regimes. Sensitivity tests were performed using simulations of two super typhoons in the Weather Research and Forecasting model. While the typhoon tracks were negligibly sensitive to the parametrization used, the typhoon intensities (the maximum lO-m wind speed and the minimum sea level pressure), sizes, and structure, were very sensitive to it.