Stable water isotopes are natural tracers quantifying the contribution of moisture recycling to local precipitation,i.e.,the moisture recycling ratio,but various isotope-based models usually lead to different results,...Stable water isotopes are natural tracers quantifying the contribution of moisture recycling to local precipitation,i.e.,the moisture recycling ratio,but various isotope-based models usually lead to different results,which affects the accuracy of local moisture recycling.In this study,a total of 18 stations from four typical areas in China were selected to compare the performance of isotope-based linear and Bayesian mixing models and to determine local moisture recycling ratio.Among the three vapor sources including advection,transpiration,and surface evaporation,the advection vapor usually played a dominant role,and the contribution of surface evaporation was less than that of transpiration.When the abnormal values were ignored,the arithmetic averages of differences between isotope-based linear and the Bayesian mixing models were 0.9%for transpiration,0.2%for surface evaporation,and–1.1%for advection,respectively,and the medians were 0.5%,0.2%,and–0.8%,respectively.The importance of transpiration was slightly less for most cases when the Bayesian mixing model was applied,and the contribution of advection was relatively larger.The Bayesian mixing model was found to perform better in determining an efficient solution since linear model sometimes resulted in negative contribution ratios.Sensitivity test with two isotope scenarios indicated that the Bayesian model had a relatively low sensitivity to the changes in isotope input,and it was important to accurately estimate the isotopes in precipitation vapor.Generally,the Bayesian mixing model should be recommended instead of a linear model.The findings are useful for understanding the performance of isotope-based linear and Bayesian mixing models under various climate backgrounds.展开更多
This paper is concerned with the aging and dependence properties in the additive hazard mixing models including some stochastic comparisons. Further, some useful bounds of reliability functions in additive hazard mixi...This paper is concerned with the aging and dependence properties in the additive hazard mixing models including some stochastic comparisons. Further, some useful bounds of reliability functions in additive hazard mixing models are obtained.展开更多
The quasi-steady methods based on mixing models have been widely applied to flow computations of turbomachinery multi- stages in aerospace engineering. Meanwhile, the unsteady numerical simulation has also been used d...The quasi-steady methods based on mixing models have been widely applied to flow computations of turbomachinery multi- stages in aerospace engineering. Meanwhile, the unsteady numerical simulation has also been used due to its ability in obtaining time-dependent flow solutions. In the paper, two different mixing treatments and the corresponding flux balanced ones are presented to exchange the flow solutions on the interfaces between adjacent blade rows. The four mixing treatments are then used for flow computations of a subsonic 1.5-stage axial turbine and a quasi-l.5-stage transonic compressor rotor. The results are compared with those by unsteady numerical method, which is implemented by using the sliding mesh technique. The effects of the quasi-steady and unsteady computation methods on the conservation of flow solutions across the interfaces are presented and addressed. Furthermore, the influence of mixing treatments on shock wave and flow separation of the transonic compressor rotor is presented in detail. All the results demonstrate that the flux balanced mixing treatments can be used for multi-stage flow computations with improved performance on interface conservation, even in the complex flows.展开更多
Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochast...Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.展开更多
The mixed distribution model is often used to extract information from heteroge-neous data and perform modeling analysis.When the density function of mixed distribution is complicated or the variable dimension is high...The mixed distribution model is often used to extract information from heteroge-neous data and perform modeling analysis.When the density function of mixed distribution is complicated or the variable dimension is high,it usually brings challenges to the parameter es-timation of the mixed distribution model.The application of MM algorithm can avoid complex expectation calculations,and can also solve the problem of high-dimensional optimization by decomposing the objective function.In this paper,MM algorithm is applied to the parameter estimation problem of mixed distribution model.The method of assembly and decomposition is used to construct the substitute function with separable parameters,which avoids the problems of complex expectation calculations and the inversion of high-dimensional matrices.展开更多
Tree height (H) in a natural stand or forest plantation is a fundamental variable in management, and the use of mathematical expressions that estimate H as a function of diameter at breast height (d) or variables at t...Tree height (H) in a natural stand or forest plantation is a fundamental variable in management, and the use of mathematical expressions that estimate H as a function of diameter at breast height (d) or variables at the stand level is a valuable support tool in forest inventories. The objective was to fit and propose a generalized H-d model for Pinus montezumae and Pinus pseudostrobus established in forest plantations of Nuevo San Juan Parangaricutiro, Michoacan, Mexico. Using nonlinear least squares (NLS), 10 generalized H-d models were fitted to 883 and 1226 pairs of H-d data from Pinus montezumae and Pinus pseudostrobus, respectively. The best model was refitted with the maximum likelihood mixed effects model (MEM) approach by including the site as a classification variable and a known variance structure. The Wang and Tang equation was selected as the best model with NLS;the MEM with an additive effect on two of its parameters and an exponential variance function improved the fit statistics for Pinus montezumae and Pinus pseudostrobus, respectively. The model validation showed equality of means among the estimates for both species and an independent subsample. The calibration of the MEM at the plot level was efficient and might increase the applicability of these results. The inclusion of dominant height in the MEM approach helped to reduce bias in the estimates and also to better explain the variability among plots.展开更多
In this paper, the measurement of liquid mixing in a downcomer of segmental type of distillation column is presented. The extent of liquid mixing is calculated by means of a mixing pool model. The results indicate tha...In this paper, the measurement of liquid mixing in a downcomer of segmental type of distillation column is presented. The extent of liquid mixing is calculated by means of a mixing pool model. The results indicate that liquid mixing in a downcomer is actually incomplete. It is a significant correction to the assumption of complete downcomer mixing or no downcomer mixing which is generally adopted in many distillation calculations. Besides, the present results are used in a two dimensional eddy diffusion model to calculate the distillation tray efficiency. It is shown that the assumption of complete downcomer mixing is closer to the actual situation than that of no downcomer mixing.展开更多
The quantitative trait loci(QTL)-by-environment(Q × E) interaction effect is hard to detect because there are no effective ways to control the genomic background. In this study, we propose a linear mixed model th...The quantitative trait loci(QTL)-by-environment(Q × E) interaction effect is hard to detect because there are no effective ways to control the genomic background. In this study, we propose a linear mixed model that simultaneously analyzes data from multiple environments to detect Q × E interactions. This model incorporates two different kinship matrices derived from the genome-wide markers to control both main and interaction polygenic background effects. Simulation studies demonstrate that our approach is more powerful than the meta-analysis and inclusive composite interval mapping methods. We further analyze four agronomic traits of rice across four environments. A main effect QTL is identified for 1000-grain weight(KGW), while no QTL are found for tiller number. Additionally, a large QTL with a significant Q × E interaction is detected on chromosome 7 affecting grain number, yield, and KGW. This region harbors two important genes, PROG1 and Ghd7. Furthermore, we apply our mixed model to analyze lodging in barley across six environments. The six regions exhibiting Q × E interaction effects identified by our approach overlap with the SNPs previously identified using EM and MCMC-based Bayesian methods, further validating the robustness of our approach. Both simulation studies and empirical data analyses show that our method outperforms all other methods compared.展开更多
Over the past decades,the expansion of natu-ral secondary forests has played a crucial role in offsetting the loss of primary forests and combating climate change.Despite this,there is a gap in our understanding of ho...Over the past decades,the expansion of natu-ral secondary forests has played a crucial role in offsetting the loss of primary forests and combating climate change.Despite this,there is a gap in our understanding of how tree species’growth and mortality patterns vary with eleva-tion in these secondary forests.In this study,we analyzed data from two censuses(spanning a five-year interval)conducted in both evergreen broadleaved forests(EBF)and temperate coniferous forests(TCF),which have been recovering for half a century,across elevation gradients in a subtropical mountain region,Mount Wuyi,China.The results indicated that the relative growth rate(RGR)of EBF(0.028±0.001 cm·cm^(-1)·a^(-1))and the mortality rate(MR)(20.03%±1.70%)were 27.3%and 16.4%higher,respec-tively,than those of TCF.Interestingly,the trade-off between RGR and MR in EBF weakened as elevation increased,a trend not observed in TCF.Conversely,TCF consistently showed a stronger trade-off between RGR and MR compared to EBF.Generalized linear mixed models revealed that ele-vation influences RGR both directly and indirectly through its interactions with slope,crown competition index(CCI),and tree canopy height(CH).However,tree mortality did not show a significant correlation with elevation.Additionally,DBH significantly influenced both tree growth and mortal-ity,whereas and CH and CCI had opposite effects on tree growth between EBF and TCF.Our study underscores the importance of elevation in shaping the population dynamics and the biomass carbon sink balance of mountain forests.These insights enhance our understanding of tree species’life strategies,enabling more accurate predictions of forest dynamics and their response to environmental changes.展开更多
Industrial ebullated-bed is an important device for promoting the cleaning and upgrading of oil products. The lumped kinetic model is a powerful tool for predicting the product yield of the ebullated-bed residue hydro...Industrial ebullated-bed is an important device for promoting the cleaning and upgrading of oil products. The lumped kinetic model is a powerful tool for predicting the product yield of the ebullated-bed residue hydrogenation (EBRH) unit, However, during the long-term operation of the device, there are phenomena such as low frequency of material property analysis leading to limited operating data and diverse operating modes at the same time scale, which poses a huge challenge to building an accurate product yield prediction model. To address these challenges, a data augmentation-based eleven lumped reaction kinetics mechanism model was constructed. This model combines generative adversarial networks, outlier elimination, and L2 norm data filtering to expand the dataset and utilizes kernel principal component analysis-fuzzy C-means for operating condition partitioning. Based on the hydrogenation reaction mechanism, a single and sub operating condition eleven lumped reaction kinetics model of an ebullated-bed residue hydrogenation unit, comprising 55 reaction paths and 110 parameters, was constructed before and after data augmentation. Compared to the single model before data enhancement, the average absolute error of the sub-models under data enhancement division was reduced by 23%. Thus, these findings can help guide the operation and optimization of the production process.展开更多
The genus Beta encompasses economically important root crops such as sugar and table beet.A Beta diversity set including the wild relative B.vulgaris ssp.maritima was grown in the field,and a large phenotypic diversit...The genus Beta encompasses economically important root crops such as sugar and table beet.A Beta diversity set including the wild relative B.vulgaris ssp.maritima was grown in the field,and a large phenotypic diversity was observed.The genomes of 290 accessions were sequenced,and more than 10 million high-quality SNPs were employed to study genetic diversity.A genome-wide association study was performed,and marker-trait associations were found for nine phenotypic traits.The candidate gene within the M locus controlling monogermity on chromosome 4 was previously unknown.The most significant association for monogermity was identified at the end of chromosome 4.Within this region,a non-synonymous mutation within the zinc-finger domain of the WIP2 gene co-segregated with monogermity.This gene plays a regulatory role in AGL8/FUL in Arabidopsis.Intriguingly,commercial hybrids are in a heterozygous state at this position.Thus,the long-sought gene for monogermity was identified in this study.Red and yellow pigmentation due to betalain accumulation in shoots and roots is an important characteristic of table and leaf beets.The strongest associations were found upstream or downstream of two genes encoding Cytochrome P450 and anthocyanin MYB-like transcription factor proteins involved in betalain biosynthesis.Significant associations for Cercospora leaf spot resistance were identified on chromosomes 1,2,7,and 9.The associated regions harbor genes encoding proteins with leucinerich repeats and nucleotide binding sites whose homologs are major constituents of plant-pathogen defense.展开更多
DNA microarray technology is an extremely effective technique for studying gene expression patterns in cells, and the main challenge currently faced by this technology is how to analyze the large amount of gene expres...DNA microarray technology is an extremely effective technique for studying gene expression patterns in cells, and the main challenge currently faced by this technology is how to analyze the large amount of gene expression data generated. To address this, this paper employs a mixed-effects model to analyze gene expression data. In terms of data selection, 1176 genes from the white mouse gene expression dataset under two experimental conditions were chosen, setting up two conditions: pneumococcal infection and no infection, and constructing a mixed-effects model. After preprocessing the gene chip information, the data were imported into the model, preliminary results were calculated, and permutation tests were performed to biologically validate the preliminary results using GSEA. The final dataset consists of 20 groups of gene expression data from pneumococcal infection, which categorizes functionally related genes based on the similarity of their expression profiles, facilitating the study of genes with unknown functions.展开更多
Demand-responsive transportation has been introduced in many cities around the world.However,whether it is applicable in the railway is still questionable,an exploration of passenger choice behavior between demandresp...Demand-responsive transportation has been introduced in many cities around the world.However,whether it is applicable in the railway is still questionable,an exploration of passenger choice behavior between demandresponsive trains and pre-scheduled trains is pivotal in addressing this issue.To delve into passengers’choice preferences when facing demand-responsive trains and to dissect the feasibility of implementing demandresponsive service in high-speed railways,the stated preference survey method is employed to investigate travel intention of passengers.Based on the survey data obtained in China,the heterogeneity of passengers is analyzed from three aspects:personal socio-economic characteristics,travel characteristics,and travel mode choice.Considering the situation that demand-responsive train cannot operate,the risk attributes are considered.To bolster the appeal of demand-responsive trains,personalized service product attributes are added.Mixed Logit mode,which takes into account the heterogeneous travel choice behavior of passengers,is developed,and Maximum Likelihood Estimation and the Monte Carlo method are used to calibrate model parameters.The willingness to pay in terms of different factors of passengers is determined.The results indicate that early arrival deviation time,late arrival deviation time,demand response time,and success rate of ticket purchase remarkable influence passengers’decision regarding demand-responsive train,with only the success rate of ticket purchase positively impacting train choice.Moreover,the significant difference in train ticket price is observed solely in the self-funded long distance scenario,while demand-responsive trains are found to be particularly appealing in self-funded short distance scenario.Through the Willingness To Pay(WTP)analysis,it is discovered that by shortening demand response time,enhancing the success rate of ticket purchase,and minimizing the deviation times of early arrival and late arrival of trains,the attractiveness of the demand-responsive train to passengers under three travel scenarios can be augmented.This study provides profound insights into the possibility of railway enterprises operating demand-responsive trains.展开更多
Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to en...Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.展开更多
The rapid evolution of unmanned aerial vehicle(UAV)technology and autonomous capabilities has positioned UAV as promising last-mile delivery means.Vehicle and onboard UAV collaborative delivery is introduced as a nove...The rapid evolution of unmanned aerial vehicle(UAV)technology and autonomous capabilities has positioned UAV as promising last-mile delivery means.Vehicle and onboard UAV collaborative delivery is introduced as a novel delivery mode.Spatiotemporal collaboration,along with energy consumption with payload and wind conditions play important roles in delivery route planning.This paper introduces the traveling salesman problem with time window and onboard UAV(TSPTWOUAV)and emphasizes the consideration of real-world scenarios,focusing on time collaboration and energy consumption with wind and payload.To address this,a mixed integer linear programming(MILP)model is formulated to minimize the energy consumption costs of vehicle and UAV.Furthermore,an adaptive large neighborhood search(ALNS)algorithm is applied to identify high-quality solutions efficiently.The effectiveness of the proposed model and algorithm is validated through numerical tests on real geographic instances and sensitivity analysis of key parameters is conducted.展开更多
The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow fi...The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry （PIV） and laser Doppler velocimetry （LDV）. The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.展开更多
The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transp...The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transported Probability Density Function(TPDF) turbulence combustion model is promising in engineering applications. In flame region, the impact of chemical reaction should be considered in TPDF molecular mixing model. Based on pioneer research, three new TPDF turbulence-chemistry dual time scale molecular mixing models were proposed tentatively by adding the chemistry time scale in molecular mixing model for nonpremixed flame. The Aero-Engine Combustor Simulation Code(AECSC) which is based on LES-TPDF method was combined with the three new models. Then the Sandia laboratory's methane-air jet flames: Flame D and Flame E were simulated. Transient simulation results show that all the three new models can predict the instantaneous combustion flow pattern of the jet flames. Furthermore,the average scalar statistical results were compared with the experimental data. The simulation result of the new TPDF arithmetic mean modification model is the closest to the experimental data:the average error in Flame D is 7.6% and 6.6% in Flame E. The extinction and re-ignition phenomena of the jet flames especially Flame E were captured. The turbulence time scale and the chemistry time scale are in different order in the whole flow field. The dual time scale TPDF combustion model has ability to deal with both the turbulence effect and the chemistry reaction effect, as well as their interaction more accurately for nonpremixed flames.展开更多
Background:We used mixed models with random components to develop height-diameter(h-d) functions for mixed,uneven-aged stands in northwestern Durango(Mexico),considering the breast height diameter(d) and stand variabl...Background:We used mixed models with random components to develop height-diameter(h-d) functions for mixed,uneven-aged stands in northwestern Durango(Mexico),considering the breast height diameter(d) and stand variables as predictors.Methods:The data were obtained from 44 permanent plots used to monitor stand growth under forest management in the study area.Results:The generalized Bertalanffy-Richards model performed better than the other generalized models in predicting the total height of the species under study.For the genera Pinus and Quercus,the models were successfully calibrated by measuring the height of a subsample of three randomly selected trees close to the mean d,whereas for species of the genera Cupressus,Arbutus and Alnus,three trees were also selected,but they are specifically the maximum,minimum and mean d trees.Conclusions:The presented equations represent a new tool for the evaluation and management of natural forest in the region.展开更多
Fluid flow and mixing of molten steel in a twin-slab-strand continuous casting tundish were investigated using a mixing model under non-isothermal conditions.This model led to a set of ordinary differential equations ...Fluid flow and mixing of molten steel in a twin-slab-strand continuous casting tundish were investigated using a mixing model under non-isothermal conditions.This model led to a set of ordinary differential equations that were solved with a Runge-Kutta algorithm.Steady state water modeling was carried out under non-isothermal conditions.Experimental data obtained from the water model were used to calibrate the mixing model.Owing to the presence of a mixed convection in the non-isothermal conditions,a channelizing flow would be created in the fluid inside the tundish.A mixing model was designed that was capable of predicting RTD（residence time distribution）curves for different cases in non-isothermal conditions.The relationship between RTD parameters and the Tu（tundish Richardson number）was obtained for various cases under non-isothermal conditions.The results show that the RTD parameters were completely different under isothermal and non-isothermal conditions.The comparison of the RTD curves between the isothermal and non-isothermal conditions presents that the extent of mixing in the tundish in non-isothermal conditions is lower than the mixing extent in isothermal conditions.展开更多
In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic ...In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic mean, respectively, are proposed. It is shown that these new estimators dominate the unbiased estimator under the squared error loss function. Finally, some simulation results to compare the performance of the proposed estimators with that of the unbiased estimator are reported. The simulation results indicate that these new shrinkage estimators provide a substantial improvement in risk under most situations.展开更多
基金This study was supported by the National Natural Science Foundation of China(42261008,41971034)the Natural Science Foundation of Gansu Province,China(22JR5RA074).
文摘Stable water isotopes are natural tracers quantifying the contribution of moisture recycling to local precipitation,i.e.,the moisture recycling ratio,but various isotope-based models usually lead to different results,which affects the accuracy of local moisture recycling.In this study,a total of 18 stations from four typical areas in China were selected to compare the performance of isotope-based linear and Bayesian mixing models and to determine local moisture recycling ratio.Among the three vapor sources including advection,transpiration,and surface evaporation,the advection vapor usually played a dominant role,and the contribution of surface evaporation was less than that of transpiration.When the abnormal values were ignored,the arithmetic averages of differences between isotope-based linear and the Bayesian mixing models were 0.9%for transpiration,0.2%for surface evaporation,and–1.1%for advection,respectively,and the medians were 0.5%,0.2%,and–0.8%,respectively.The importance of transpiration was slightly less for most cases when the Bayesian mixing model was applied,and the contribution of advection was relatively larger.The Bayesian mixing model was found to perform better in determining an efficient solution since linear model sometimes resulted in negative contribution ratios.Sensitivity test with two isotope scenarios indicated that the Bayesian model had a relatively low sensitivity to the changes in isotope input,and it was important to accurately estimate the isotopes in precipitation vapor.Generally,the Bayesian mixing model should be recommended instead of a linear model.The findings are useful for understanding the performance of isotope-based linear and Bayesian mixing models under various climate backgrounds.
基金Supported by the Scientific Research Foundation of Hebei University of Science and Technology
文摘This paper is concerned with the aging and dependence properties in the additive hazard mixing models including some stochastic comparisons. Further, some useful bounds of reliability functions in additive hazard mixing models are obtained.
基金supported by the National Natural Science Foundation of China(Grant Nos.51376009&51676003)
文摘The quasi-steady methods based on mixing models have been widely applied to flow computations of turbomachinery multi- stages in aerospace engineering. Meanwhile, the unsteady numerical simulation has also been used due to its ability in obtaining time-dependent flow solutions. In the paper, two different mixing treatments and the corresponding flux balanced ones are presented to exchange the flow solutions on the interfaces between adjacent blade rows. The four mixing treatments are then used for flow computations of a subsonic 1.5-stage axial turbine and a quasi-l.5-stage transonic compressor rotor. The results are compared with those by unsteady numerical method, which is implemented by using the sliding mesh technique. The effects of the quasi-steady and unsteady computation methods on the conservation of flow solutions across the interfaces are presented and addressed. Furthermore, the influence of mixing treatments on shock wave and flow separation of the transonic compressor rotor is presented in detail. All the results demonstrate that the flux balanced mixing treatments can be used for multi-stage flow computations with improved performance on interface conservation, even in the complex flows.
基金supported by the National Natural Science Foundation of China(Grant Nos.82173620 to Yang Zhao and 82041024 to Feng Chen)partially supported by the Bill&Melinda Gates Foundation(Grant No.INV-006371 to Feng Chen)Priority Academic Program Development of Jiangsu Higher Education Institutions.
文摘Deterministic compartment models(CMs)and stochastic models,including stochastic CMs and agent-based models,are widely utilized in epidemic modeling.However,the relationship between CMs and their corresponding stochastic models is not well understood.The present study aimed to address this gap by conducting a comparative study using the susceptible,exposed,infectious,and recovered(SEIR)model and its extended CMs from the coronavirus disease 2019 modeling literature.We demonstrated the equivalence of the numerical solution of CMs using the Euler scheme and their stochastic counterparts through theoretical analysis and simulations.Based on this equivalence,we proposed an efficient model calibration method that could replicate the exact solution of CMs in the corresponding stochastic models through parameter adjustment.The advancement in calibration techniques enhanced the accuracy of stochastic modeling in capturing the dynamics of epidemics.However,it should be noted that discrete-time stochastic models cannot perfectly reproduce the exact solution of continuous-time CMs.Additionally,we proposed a new stochastic compartment and agent mixed model as an alternative to agent-based models for large-scale population simulations with a limited number of agents.This model offered a balance between computational efficiency and accuracy.The results of this research contributed to the comparison and unification of deterministic CMs and stochastic models in epidemic modeling.Furthermore,the results had implications for the development of hybrid models that integrated the strengths of both frameworks.Overall,the present study has provided valuable epidemic modeling techniques and their practical applications for understanding and controlling the spread of infectious diseases.
基金Supported by the National Natural Science Foundation of China(12261108)the General Program of Basic Research Programs of Yunnan Province(202401AT070126)+1 种基金the Yunnan Key Laboratory of Modern Analytical Mathematics and Applications(202302AN360007)the Cross-integration Innovation team of modern Applied Mathematics and Life Sciences in Yunnan Province,China(202405AS350003).
文摘The mixed distribution model is often used to extract information from heteroge-neous data and perform modeling analysis.When the density function of mixed distribution is complicated or the variable dimension is high,it usually brings challenges to the parameter es-timation of the mixed distribution model.The application of MM algorithm can avoid complex expectation calculations,and can also solve the problem of high-dimensional optimization by decomposing the objective function.In this paper,MM algorithm is applied to the parameter estimation problem of mixed distribution model.The method of assembly and decomposition is used to construct the substitute function with separable parameters,which avoids the problems of complex expectation calculations and the inversion of high-dimensional matrices.
文摘Tree height (H) in a natural stand or forest plantation is a fundamental variable in management, and the use of mathematical expressions that estimate H as a function of diameter at breast height (d) or variables at the stand level is a valuable support tool in forest inventories. The objective was to fit and propose a generalized H-d model for Pinus montezumae and Pinus pseudostrobus established in forest plantations of Nuevo San Juan Parangaricutiro, Michoacan, Mexico. Using nonlinear least squares (NLS), 10 generalized H-d models were fitted to 883 and 1226 pairs of H-d data from Pinus montezumae and Pinus pseudostrobus, respectively. The best model was refitted with the maximum likelihood mixed effects model (MEM) approach by including the site as a classification variable and a known variance structure. The Wang and Tang equation was selected as the best model with NLS;the MEM with an additive effect on two of its parameters and an exponential variance function improved the fit statistics for Pinus montezumae and Pinus pseudostrobus, respectively. The model validation showed equality of means among the estimates for both species and an independent subsample. The calibration of the MEM at the plot level was efficient and might increase the applicability of these results. The inclusion of dominant height in the MEM approach helped to reduce bias in the estimates and also to better explain the variability among plots.
文摘In this paper, the measurement of liquid mixing in a downcomer of segmental type of distillation column is presented. The extent of liquid mixing is calculated by means of a mixing pool model. The results indicate that liquid mixing in a downcomer is actually incomplete. It is a significant correction to the assumption of complete downcomer mixing or no downcomer mixing which is generally adopted in many distillation calculations. Besides, the present results are used in a two dimensional eddy diffusion model to calculate the distillation tray efficiency. It is shown that the assumption of complete downcomer mixing is closer to the actual situation than that of no downcomer mixing.
基金supported by the National Key Research and Development Programs of China(2024YFF1000100 and 2021YFD1301102)the National Natural Science Foundation of China (32172702)+3 种基金the State Key Laboratory of Animal Biotech Breeding (XQSWYZQZ-KFYX-4)Zaozhuang Elite Industrial Innovation ProgramAgricultural Science and Technology Innovation Program (ASTIP-IAS-TS-6)supported by the United States National Science Foundation (NSF) Collaborative Research Grant (DBI-1458515)
文摘The quantitative trait loci(QTL)-by-environment(Q × E) interaction effect is hard to detect because there are no effective ways to control the genomic background. In this study, we propose a linear mixed model that simultaneously analyzes data from multiple environments to detect Q × E interactions. This model incorporates two different kinship matrices derived from the genome-wide markers to control both main and interaction polygenic background effects. Simulation studies demonstrate that our approach is more powerful than the meta-analysis and inclusive composite interval mapping methods. We further analyze four agronomic traits of rice across four environments. A main effect QTL is identified for 1000-grain weight(KGW), while no QTL are found for tiller number. Additionally, a large QTL with a significant Q × E interaction is detected on chromosome 7 affecting grain number, yield, and KGW. This region harbors two important genes, PROG1 and Ghd7. Furthermore, we apply our mixed model to analyze lodging in barley across six environments. The six regions exhibiting Q × E interaction effects identified by our approach overlap with the SNPs previously identified using EM and MCMC-based Bayesian methods, further validating the robustness of our approach. Both simulation studies and empirical data analyses show that our method outperforms all other methods compared.
基金funded by the National Natural Science Foundation of China(Grant No.32271872).
文摘Over the past decades,the expansion of natu-ral secondary forests has played a crucial role in offsetting the loss of primary forests and combating climate change.Despite this,there is a gap in our understanding of how tree species’growth and mortality patterns vary with eleva-tion in these secondary forests.In this study,we analyzed data from two censuses(spanning a five-year interval)conducted in both evergreen broadleaved forests(EBF)and temperate coniferous forests(TCF),which have been recovering for half a century,across elevation gradients in a subtropical mountain region,Mount Wuyi,China.The results indicated that the relative growth rate(RGR)of EBF(0.028±0.001 cm·cm^(-1)·a^(-1))and the mortality rate(MR)(20.03%±1.70%)were 27.3%and 16.4%higher,respec-tively,than those of TCF.Interestingly,the trade-off between RGR and MR in EBF weakened as elevation increased,a trend not observed in TCF.Conversely,TCF consistently showed a stronger trade-off between RGR and MR compared to EBF.Generalized linear mixed models revealed that ele-vation influences RGR both directly and indirectly through its interactions with slope,crown competition index(CCI),and tree canopy height(CH).However,tree mortality did not show a significant correlation with elevation.Additionally,DBH significantly influenced both tree growth and mortal-ity,whereas and CH and CCI had opposite effects on tree growth between EBF and TCF.Our study underscores the importance of elevation in shaping the population dynamics and the biomass carbon sink balance of mountain forests.These insights enhance our understanding of tree species’life strategies,enabling more accurate predictions of forest dynamics and their response to environmental changes.
基金supported by National Natural Science Foundation of China(Basic Science Center Program:61988101)National Natural Science Foundation of China(62394345,62373155,62173147)the Major Science and Technology Project of Xinjiang(No.2022A01006-4).
文摘Industrial ebullated-bed is an important device for promoting the cleaning and upgrading of oil products. The lumped kinetic model is a powerful tool for predicting the product yield of the ebullated-bed residue hydrogenation (EBRH) unit, However, during the long-term operation of the device, there are phenomena such as low frequency of material property analysis leading to limited operating data and diverse operating modes at the same time scale, which poses a huge challenge to building an accurate product yield prediction model. To address these challenges, a data augmentation-based eleven lumped reaction kinetics mechanism model was constructed. This model combines generative adversarial networks, outlier elimination, and L2 norm data filtering to expand the dataset and utilizes kernel principal component analysis-fuzzy C-means for operating condition partitioning. Based on the hydrogenation reaction mechanism, a single and sub operating condition eleven lumped reaction kinetics model of an ebullated-bed residue hydrogenation unit, comprising 55 reaction paths and 110 parameters, was constructed before and after data augmentation. Compared to the single model before data enhancement, the average absolute error of the sub-models under data enhancement division was reduced by 23%. Thus, these findings can help guide the operation and optimization of the production process.
基金the financial support from the German Research Foundation(Deutsche Forschungsgemeinschaft,DFG)—Project Number 400993799(Project 2 within the Research Training Group 2501 Translational Evolutionary Research,https://gepris.dfg.de/gepris/projekt/400993799)supported by the BMBF-funded de.NBI Cloud,part of the German Network for Bioinformatics(de.NBI)。
文摘The genus Beta encompasses economically important root crops such as sugar and table beet.A Beta diversity set including the wild relative B.vulgaris ssp.maritima was grown in the field,and a large phenotypic diversity was observed.The genomes of 290 accessions were sequenced,and more than 10 million high-quality SNPs were employed to study genetic diversity.A genome-wide association study was performed,and marker-trait associations were found for nine phenotypic traits.The candidate gene within the M locus controlling monogermity on chromosome 4 was previously unknown.The most significant association for monogermity was identified at the end of chromosome 4.Within this region,a non-synonymous mutation within the zinc-finger domain of the WIP2 gene co-segregated with monogermity.This gene plays a regulatory role in AGL8/FUL in Arabidopsis.Intriguingly,commercial hybrids are in a heterozygous state at this position.Thus,the long-sought gene for monogermity was identified in this study.Red and yellow pigmentation due to betalain accumulation in shoots and roots is an important characteristic of table and leaf beets.The strongest associations were found upstream or downstream of two genes encoding Cytochrome P450 and anthocyanin MYB-like transcription factor proteins involved in betalain biosynthesis.Significant associations for Cercospora leaf spot resistance were identified on chromosomes 1,2,7,and 9.The associated regions harbor genes encoding proteins with leucinerich repeats and nucleotide binding sites whose homologs are major constituents of plant-pathogen defense.
文摘DNA microarray technology is an extremely effective technique for studying gene expression patterns in cells, and the main challenge currently faced by this technology is how to analyze the large amount of gene expression data generated. To address this, this paper employs a mixed-effects model to analyze gene expression data. In terms of data selection, 1176 genes from the white mouse gene expression dataset under two experimental conditions were chosen, setting up two conditions: pneumococcal infection and no infection, and constructing a mixed-effects model. After preprocessing the gene chip information, the data were imported into the model, preliminary results were calculated, and permutation tests were performed to biologically validate the preliminary results using GSEA. The final dataset consists of 20 groups of gene expression data from pneumococcal infection, which categorizes functionally related genes based on the similarity of their expression profiles, facilitating the study of genes with unknown functions.
基金supported by the National Natural Science Foundation of China(No.72471023,71971019)the Fundamental Research Funds for the Central Universities(No.2024QYBS025).
文摘Demand-responsive transportation has been introduced in many cities around the world.However,whether it is applicable in the railway is still questionable,an exploration of passenger choice behavior between demandresponsive trains and pre-scheduled trains is pivotal in addressing this issue.To delve into passengers’choice preferences when facing demand-responsive trains and to dissect the feasibility of implementing demandresponsive service in high-speed railways,the stated preference survey method is employed to investigate travel intention of passengers.Based on the survey data obtained in China,the heterogeneity of passengers is analyzed from three aspects:personal socio-economic characteristics,travel characteristics,and travel mode choice.Considering the situation that demand-responsive train cannot operate,the risk attributes are considered.To bolster the appeal of demand-responsive trains,personalized service product attributes are added.Mixed Logit mode,which takes into account the heterogeneous travel choice behavior of passengers,is developed,and Maximum Likelihood Estimation and the Monte Carlo method are used to calibrate model parameters.The willingness to pay in terms of different factors of passengers is determined.The results indicate that early arrival deviation time,late arrival deviation time,demand response time,and success rate of ticket purchase remarkable influence passengers’decision regarding demand-responsive train,with only the success rate of ticket purchase positively impacting train choice.Moreover,the significant difference in train ticket price is observed solely in the self-funded long distance scenario,while demand-responsive trains are found to be particularly appealing in self-funded short distance scenario.Through the Willingness To Pay(WTP)analysis,it is discovered that by shortening demand response time,enhancing the success rate of ticket purchase,and minimizing the deviation times of early arrival and late arrival of trains,the attractiveness of the demand-responsive train to passengers under three travel scenarios can be augmented.This study provides profound insights into the possibility of railway enterprises operating demand-responsive trains.
基金supported by the National Natural Science Foundation of China(No.92371206)the Postgraduate Scientific Research Innovation Project of Hunan Province,China(No.CX2023063).
文摘Satellite Component Layout Optimization(SCLO) is crucial in satellite system design.This paper proposes a novel Satellite Three-Dimensional Component Assignment and Layout Optimization(3D-SCALO) problem tailored to engineering requirements, aiming to optimize satellite heat dissipation while considering constraints on static stability, 3D geometric relationships between components, and special component positions. The 3D-SCALO problem is a challenging bilevel combinatorial optimization task, involving the optimization of discrete component assignment variables in the outer layer and continuous component position variables in the inner layer,with both influencing each other. To address this issue, first, a Mixed Integer Programming(MIP) model is proposed, which reformulates the original bilevel problem into a single-level optimization problem, enabling the exploration of a more comprehensive optimization space while avoiding iterative nested optimization. Then, to model the 3D geometric relationships between components within the MIP framework, a linearized 3D Phi-function method is proposed, which handles non-overlapping and safety distance constraints between cuboid components in an explicit and effective way. Subsequently, the Finite-Rectangle Method(FRM) is proposed to manage 3D geometric constraints for complex-shaped components by approximating them with a finite set of cuboids, extending the applicability of the geometric modeling approach. Finally, the feasibility and effectiveness of the proposed MIP model are demonstrated through two numerical examples"and a real-world engineering case, which confirms its suitability for complex-shaped components and real engineering applications.
基金Fundamental Research Funds for the Central Universities(2024JBZX038)National Natural Science F oundation of China(62076023)。
文摘The rapid evolution of unmanned aerial vehicle(UAV)technology and autonomous capabilities has positioned UAV as promising last-mile delivery means.Vehicle and onboard UAV collaborative delivery is introduced as a novel delivery mode.Spatiotemporal collaboration,along with energy consumption with payload and wind conditions play important roles in delivery route planning.This paper introduces the traveling salesman problem with time window and onboard UAV(TSPTWOUAV)and emphasizes the consideration of real-world scenarios,focusing on time collaboration and energy consumption with wind and payload.To address this,a mixed integer linear programming(MILP)model is formulated to minimize the energy consumption costs of vehicle and UAV.Furthermore,an adaptive large neighborhood search(ALNS)algorithm is applied to identify high-quality solutions efficiently.The effectiveness of the proposed model and algorithm is validated through numerical tests on real geographic instances and sensitivity analysis of key parameters is conducted.
基金supported by National Natural Science Foundation of China (Grant Nos. 51139007, 51079151, 51079152)Research Fundfor the Doctoral Program of Higher Education of China (Grant No. 0100008110012)
文摘The current research of large eddy simulation （LES） of turbulent flow in pumps mainly concentrates in applying conventional subgrid-scale （SGS） model to simulate turbulent flow, which aims at obtaining the flow field in pump. The selection of SGS model is usually not considered seriously, so the accuracy and efficiency of the simulation cannot be ensured. Three SGS models including Smagorinsky-Lilly model, dynamic Smagorinsky model and dynamic mixed model are comparably studied by using the commercial CFD code Fluent combined with its user define function. The simulations are performed for the turbulent flow in a centrifugal pump impeller. The simulation results indicate that the mean flows predicted by the three SGS models agree well with the experimental data obtained from the test that detailed measurements of the flow inside the rotating passages of a six-bladed shrouded centrifugal pump impeller performed using particle image velocimetry （PIV） and laser Doppler velocimetry （LDV）. The comparable results show that dynamic mixed model gives the most accurate results for mean flow in the centrifugal pump impeller. The SGS stress of dynamic mixed model is decompose into the scale similar part and the eddy viscous part. The scale similar part of SGS stress plays a significant role in high curvature regions, such as the leading edge and training edge of pump blade. It is also found that the dynamic mixed model is more adaptive to compute turbulence in the pump impeller. The research results presented is useful to improve the computational accuracy and efficiency of LES for centrifugal pumps, and provide important reference for carrying out simulation in similar fluid machineries.
基金co-supported by the National Key R&D Program of China(Nos.2017YFB0202400 and 2017YFB0202402)the National Natural Science Foundation of China(No.91741125)the Project of Newton International Fellowship Alumnus from Royal Society(No.AL120003)
文摘The numerical simulation of modern aero-engine combustion chamber needs accurate description of the interaction between turbulence and chemical reaction mechanism. The Large Eddy Simulation(LES) method with the Transported Probability Density Function(TPDF) turbulence combustion model is promising in engineering applications. In flame region, the impact of chemical reaction should be considered in TPDF molecular mixing model. Based on pioneer research, three new TPDF turbulence-chemistry dual time scale molecular mixing models were proposed tentatively by adding the chemistry time scale in molecular mixing model for nonpremixed flame. The Aero-Engine Combustor Simulation Code(AECSC) which is based on LES-TPDF method was combined with the three new models. Then the Sandia laboratory's methane-air jet flames: Flame D and Flame E were simulated. Transient simulation results show that all the three new models can predict the instantaneous combustion flow pattern of the jet flames. Furthermore,the average scalar statistical results were compared with the experimental data. The simulation result of the new TPDF arithmetic mean modification model is the closest to the experimental data:the average error in Flame D is 7.6% and 6.6% in Flame E. The extinction and re-ignition phenomena of the jet flames especially Flame E were captured. The turbulence time scale and the chemistry time scale are in different order in the whole flow field. The dual time scale TPDF combustion model has ability to deal with both the turbulence effect and the chemistry reaction effect, as well as their interaction more accurately for nonpremixed flames.
基金financially supported by the"Programa de Mejoramiento del Profesorado"(project:Seguimiento y Evaluacion de Sitios Permanentes de Investigación Forestal y el Impacto Socioeconómico delManejo Forestal en Norte de México)supported by"Programa Banco Santander-USC"(becas para estancias predoctorales destinadas a docentes e investigadores de America Latina)
文摘Background:We used mixed models with random components to develop height-diameter(h-d) functions for mixed,uneven-aged stands in northwestern Durango(Mexico),considering the breast height diameter(d) and stand variables as predictors.Methods:The data were obtained from 44 permanent plots used to monitor stand growth under forest management in the study area.Results:The generalized Bertalanffy-Richards model performed better than the other generalized models in predicting the total height of the species under study.For the genera Pinus and Quercus,the models were successfully calibrated by measuring the height of a subsample of three randomly selected trees close to the mean d,whereas for species of the genera Cupressus,Arbutus and Alnus,three trees were also selected,but they are specifically the maximum,minimum and mean d trees.Conclusions:The presented equations represent a new tool for the evaluation and management of natural forest in the region.
文摘Fluid flow and mixing of molten steel in a twin-slab-strand continuous casting tundish were investigated using a mixing model under non-isothermal conditions.This model led to a set of ordinary differential equations that were solved with a Runge-Kutta algorithm.Steady state water modeling was carried out under non-isothermal conditions.Experimental data obtained from the water model were used to calibrate the mixing model.Owing to the presence of a mixed convection in the non-isothermal conditions,a channelizing flow would be created in the fluid inside the tundish.A mixing model was designed that was capable of predicting RTD（residence time distribution）curves for different cases in non-isothermal conditions.The relationship between RTD parameters and the Tu（tundish Richardson number）was obtained for various cases under non-isothermal conditions.The results show that the RTD parameters were completely different under isothermal and non-isothermal conditions.The comparison of the RTD curves between the isothermal and non-isothermal conditions presents that the extent of mixing in the tundish in non-isothermal conditions is lower than the mixing extent in isothermal conditions.
基金supported by the Funding Project for Academic Human Resources Development in Institutions of Higher Learning Under the Jurisdiction of Beijing Municipality (0506011200702)National Natural Science Foundation of China+2 种基金Tian Yuan Special Foundation (10926059)Foundation of Zhejiang Educational Committee (Y200803920)Scientific Research Foundation of Hangzhou Dianzi University(KYS025608094)
文摘In this article, the problem of estimating the covariance matrix in general linear mixed models is considered. Two new classes of estimators obtained by shrinking the eigenvalues towards the origin and the arithmetic mean, respectively, are proposed. It is shown that these new estimators dominate the unbiased estimator under the squared error loss function. Finally, some simulation results to compare the performance of the proposed estimators with that of the unbiased estimator are reported. The simulation results indicate that these new shrinkage estimators provide a substantial improvement in risk under most situations.
