A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using...A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using the well-known EBU-Arrhenius model and the original second-order moment model. The comparison shows the advantage of the new model that it requires almost the same computational storage and time as that of the original second-order moment model, but its modeling results are in better agreement with experiments than those using other models. Hence, the new second-order moment model is promising in modeling turbulent combustion with NOx formation with finite reaction rate for engineering application.展开更多
The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collisi...The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θ model in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θ model has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θ model was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.展开更多
A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the b...A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.展开更多
A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale flu...A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the sin- gle-scale second-order moment two-phase turbulence model superposed with a particle collision model (USM-θ model) in most regions.展开更多
A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluc...A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.展开更多
Turbulent dispersed multiphase flows,including gas-particle,gas-droplet and bubble-liquid flows,are widely encountered in various engineering facilities.Modeling of two-phase turbulence,in particular the dispersed pha...Turbulent dispersed multiphase flows,including gas-particle,gas-droplet and bubble-liquid flows,are widely encountered in various engineering facilities.Modeling of two-phase turbulence,in particular the dispersed phase turbulence,is the key problem in the Eulerian-Eulerian simulation of practical dispersed multiphase flows.Although different models were developed and used,the experimental validation shows that they cannot always give satisfactory prediction results.In this paper the present author give a detailed review of the unified second-order moment (USM),k-k p and nonlinear k-k p two-phase turbulence models,proposed by him.The derivation and closure of these models are described in detail and the experimental validation and application of these models are extensively discussed.展开更多
This work presents a comprehensive second-order predictive modeling (PM) methodology based on the maximum entropy (MaxEnt) principle for obtaining best-estimate mean values and correlations for model responses and par...This work presents a comprehensive second-order predictive modeling (PM) methodology based on the maximum entropy (MaxEnt) principle for obtaining best-estimate mean values and correlations for model responses and parameters. This methodology is designated by the acronym 2<sup>nd</sup>-BERRU-PMP, where the attribute “2<sup>nd</sup>” indicates that this methodology incorporates second- order uncertainties (means and covariances) and second (and higher) order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best-Estimate Results with Reduced Uncertainties” and the last letter (“P”) in the acronym indicates “probabilistic,” referring to the MaxEnt probabilistic inclusion of the computational model responses. This is in contradistinction to the 2<sup>nd</sup>-BERRU-PMD methodology, which deterministically combines the computed model responses with the experimental information, as presented in the accompanying work (Part I). Although both the 2<sup>nd</sup>-BERRU-PMP and the 2<sup>nd</sup>-BERRU-PMD methodologies yield expressions that include second (and higher) order sensitivities of responses to model parameters, the respective expressions for the predicted responses, for the calibrated predicted parameters and for their predicted uncertainties (covariances), are not identical to each other. Nevertheless, the results predicted by both the 2<sup>nd</sup>-BERRU-PMP and the 2<sup>nd</sup>-BERRU-PMD methodologies encompass, as particular cases, the results produced by the extant data assimilation and data adjustment procedures, which rely on the minimization, in a least-square sense, of a user-defined functional meant to represent the discrepancies between measured and computed model responses.展开更多
This work presents a comprehensive second-order predictive modeling (PM) methodology designated by the acronym 2<sup>nd</sup>-BERRU-PMD. The attribute “2<sup>nd</sup>” indicates that this met...This work presents a comprehensive second-order predictive modeling (PM) methodology designated by the acronym 2<sup>nd</sup>-BERRU-PMD. The attribute “2<sup>nd</sup>” indicates that this methodology incorporates second-order uncertainties (means and covariances) and second-order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best- Estimate Results with Reduced Uncertainties” and the last letter (“D”) in the acronym indicates “deterministic,” referring to the deterministic inclusion of the computational model responses. The 2<sup>nd</sup>-BERRU-PMD methodology is fundamentally based on the maximum entropy (MaxEnt) principle. This principle is in contradistinction to the fundamental principle that underlies the extant data assimilation and/or adjustment procedures which minimize in a least-square sense a subjective user-defined functional which is meant to represent the discrepancies between measured and computed model responses. It is shown that the 2<sup>nd</sup>-BERRU-PMD methodology generalizes and extends current data assimilation and/or data adjustment procedures while overcoming the fundamental limitations of these procedures. In the accompanying work (Part II), the alternative framework for developing the “second- order MaxEnt predictive modelling methodology” is presented by incorporating probabilistically (as opposed to “deterministically”) the computed model responses.展开更多
Dense gas-particle flows are frequently encountered in fluidized beds,riser and downer reactors,pneumatic transport and the near-wall zone of dilute gas-particle flows.Particle-particle collision plays an important ro...Dense gas-particle flows are frequently encountered in fluidized beds,riser and downer reactors,pneumatic transport and the near-wall zone of dilute gas-particle flows.Particle-particle collision plays an important role in the behavior of two-phase flows.In this paper a USM-Q two-phase turbulence model for dense gas-particle flows is proposed to account for both two-phase turbulence and inter-particle collision.For two-fluid large-eddy simulation of gas-particle flows,the author proposed a unified second-order moment(USM) two-phase SGS stress model and a two-phase k-kp SGS energy-equation stress model.The proposed models can fully account for the interaction between the gas and particle SGS stresses.展开更多
This work illustrates the innovative results obtained by applying the recently developed the 2<sup>nd</sup>-order predictive modeling methodology called “2<sup>nd</sup>- BERRU-PM”, where the ...This work illustrates the innovative results obtained by applying the recently developed the 2<sup>nd</sup>-order predictive modeling methodology called “2<sup>nd</sup>- BERRU-PM”, where the acronym BERRU denotes “best-estimate results with reduced uncertainties” and “PM” denotes “predictive modeling.” The physical system selected for this illustrative application is a polyethylene-reflected plutonium (acronym: PERP) OECD/NEA reactor physics benchmark. This benchmark is modeled using the neutron transport Boltzmann equation (involving 21,976 uncertain parameters), the solution of which is representative of “large-scale computations.” The results obtained in this work confirm the fact that the 2<sup>nd</sup>-BERRU-PM methodology predicts best-estimate results that fall in between the corresponding computed and measured values, while reducing the predicted standard deviations of the predicted results to values smaller than either the experimentally measured or the computed values of the respective standard deviations. The obtained results also indicate that 2<sup>nd</sup>-order response sensitivities must always be included to quantify the need for including (or not) the 3<sup>rd</sup>- and/or 4<sup>th</sup>-order sensitivities. When the parameters are known with high precision, the contributions of the higher-order sensitivities diminish with increasing order, so that the inclusion of the 1<sup>st</sup>- and 2<sup>nd</sup>-order sensitivities may suffice for obtaining accurate predicted best- estimate response values and best-estimate standard deviations. On the other hand, when the parameters’ standard deviations are sufficiently large to approach (or be outside of) the radius of convergence of the multivariate Taylor-series which represents the response in the phase-space of model parameters, the contributions stemming from the 3<sup>rd</sup>- and even 4<sup>th</sup>-order sensitivities are necessary to ensure consistency between the computed and measured response. In such cases, the use of only the 1<sup>st</sup>-order sensitivities erroneously indicates that the computed results are inconsistent with the respective measured response. Ongoing research aims at extending the 2<sup>nd</sup>-BERRU-PM methodology to fourth-order, thus enabling the computation of third-order response correlations (skewness) and fourth-order response correlations (kurtosis).展开更多
Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for ...Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.展开更多
The selection of wave force models will significantly impact the structural responses of floating wind turbines.In this study,comparisons of wave force model effects on the structural responses and fatigue loads of a ...The selection of wave force models will significantly impact the structural responses of floating wind turbines.In this study,comparisons of wave force model effects on the structural responses and fatigue loads of a semi-submersible floating wind turbine(SFWT)were conducted.Simulations were performed by employing the Morison equation(ME)with linear or second-order wave kinematics and potential flow theory(PFT)with first-or second-order wave forces.A comparison of regular waves,irregular waves,and coupled wind/waves analyses with the experimental data showed that many of the simulation results and experimental data are relatively consistent.However,notable discrepancies are found in the response amplitude operators for platform heave,tower base bending moment,and tension in mooring lines.PFT models give more satisfactory results of heave but more significant discrepan-cies in tower base bending moment than the ME models.In irregular wave analyses,low-frequency resonances were captured by PFT models with second-order difference-frequency terms,and high-frequency resonances were captured by the ME models or PFT models with second-order sum-frequency terms.These force models capture the response frequencies but do not reasonably predict the response amplitudes.The coupled wind/waves analyses showed more satisfactory results than the wave-only analyses.However,an important detail to note is that this satisfactory result is based on the overprediction of wind-induced responses.展开更多
With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Ya...With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples:turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder (GL) model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.展开更多
We propose that the QCD vacuum pion tetrahedron condensate density vary in space and drops to extremely low values in the Kennan, Barger and Cowie (KBC) void in analogy to earth’s atmospheric density drop with elevat...We propose that the QCD vacuum pion tetrahedron condensate density vary in space and drops to extremely low values in the Kennan, Barger and Cowie (KBC) void in analogy to earth’s atmospheric density drop with elevation from earth. We propose a formula for the gravitation acceleration based on the non-uniform pion tetrahedron condensate. Gravity may be due to the underlying microscopic attraction between quarks and antiquarks, which are part of the vacuum pion tetrahedron condensate. We propose an electron tetrahedron model, where electrons are comprised of tetraquark tetrahedrons, and . The quarks determine the negative electron charge and the or quarks determine the electron two spin states. The electron tetrahedron may perform a high frequency quark exchange reactions with the pion tetrahedron condensate by tunneling through the condensation gap creating a delocalized electron cloud with a fixed spin. The pion tetrahedron may act as a QCD glue bonding electron pairs in atoms and molecules and protons to neutrons in the nuclei. Conservation of valence quarks and antiquarks is proposed.展开更多
In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is deriv...In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.展开更多
We apply the second-order Born-Oppenheimer (BO) approximation to investigate the dynamics of the Rabi model, which describes the interaction between a two-level system and a single bosonic mode beyond the rotating w...We apply the second-order Born-Oppenheimer (BO) approximation to investigate the dynamics of the Rabi model, which describes the interaction between a two-level system and a single bosonic mode beyond the rotating wave approxi- mation. By comparing with the numerical results, we find that our approach works well when the frequency of the two-level system is much smaller than that of the bosonic mode.展开更多
In this paper, an attitude maneuver control problem is investigated for a rigid spacecraft using an array of two variable speed control moment gyroscopes (VSCMGs) with gimbal axes skewed to each other. A mathematica...In this paper, an attitude maneuver control problem is investigated for a rigid spacecraft using an array of two variable speed control moment gyroscopes (VSCMGs) with gimbal axes skewed to each other. A mathematical model is constructed by taking the spacecraft and the gyroscopes together as an integrated system, with the coupling interaction between them considered. To overcome the singular issues of the VSCMGs due to the conventional torque-based method, the first-order derivative of gimbal rates and the second-order derivative of the rotor spinning velocity, instead of the gyroscope torques, are taken as input variables. Moreover, taking external disturbances into account, a feedback control law is designed for the system based on a method of nonlinear model predictive control (NMPC). The attitude maneuver can be realized fast and smoothly by using the proposed controller in this paper.展开更多
In this paper, we investigate the solution moment stability for a Harrison-type predator-prey model with parametric dichotomous noises. Using the Shapiro-Loginov formula, the equations for the first-order and second-o...In this paper, we investigate the solution moment stability for a Harrison-type predator-prey model with parametric dichotomous noises. Using the Shapiro-Loginov formula, the equations for the first-order and second-order moments are obtained and the corresponding stable conditions are given. It is found that the solution moment stability depends on the noise intensity and correlation time of noise. The first-order and second-order moments become unstable with the decrease of correlation time. That is, the dichotomous noise can improve the solution moment stability with respect to Gaussian white noise. Finally, some numerical results are presented to verify the theoretical analyses.展开更多
This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>...This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.展开更多
A two-moment bulk stratiform microphysics scheme, including recently developed physically-based droplet activation/ice nucleation parameterizations has been implemented into the Grid-point Atmospheric Model of IAP LA...A two-moment bulk stratiform microphysics scheme, including recently developed physically-based droplet activation/ice nucleation parameterizations has been implemented into the Grid-point Atmospheric Model of IAP LASG (GAMIL) as an effort to enhance the model's capability to simulate aerosol indirect effects. Unlike the previous one-moment cloud microphysics scheme, the new scheme produces a reasonable rep- resentation of cloud particle size and number concentration. This scheme captures the observed spatial variations in cloud droplet number concentrations. Simulated ice crystal number concentrations in cirrus clouds qualitatively agree with in situ observations. The longwave and shortwave cloud forcings are in better agreement with observations. Sensitivity tests show that the column cloud droplet number concentrations calculated from two different droplet activation parameterizations are similar. However, ice crystal number concentration in mixed-phased clouds is sensitive to different heterogeneous ice nucleation formulations. The simulation with high ice crystal number concentration in mixed-phase clouds has less liquid water path and weaker cloud forcing. ~rthermore, ice crystal number concentration in cirrus clouds is sensitive to different ice nucleation parameterizations. Sensitivity tests also suggest that the impact of pre-existing ice crystals on homogeneous freezing in old clouds should be taken into account.展开更多
基金The project sponsored by the Foundation for Doctorate Thesis of Tsinghua Universitythe National Key Project in 1999-2004 sponsored by the Ministry of Science and Technology of China
文摘A new second-order moment model for turbulent combustion is applied in the simulation of methane-air turbulent jet flame. The predicted results are compared with the experimental results and with those predicted using the well-known EBU-Arrhenius model and the original second-order moment model. The comparison shows the advantage of the new model that it requires almost the same computational storage and time as that of the original second-order moment model, but its modeling results are in better agreement with experiments than those using other models. Hence, the new second-order moment model is promising in modeling turbulent combustion with NOx formation with finite reaction rate for engineering application.
基金Project supported by the National Key Basic Research and Development Program of China(No.G1999-0222-08)
文摘The USM-θ model of Bingham fluid for dense two-phase turbulent flow was developed, which combines the second-order moment model for two-phase turbulence with the particle kinetic theory for the inter-particle collision. In this model, phases interaction and the extra term of Bingham fluid yield stress are taken into account. An algorithm for USM-θ model in dense two-phase flow was proposed, in which the influence of particle volume fraction is accounted for. This model was used to simulate turbulent flow of Bingham fluid single-phase and dense liquid-particle two-phase in pipe. It is shown USM-θ model has better prediction result than the five-equation model, in which the particle-particle collision is modeled by the particle kinetic theory, while the turbulence of both phase is simulated by the two-equation turbulence model. The USM-θ model was then used to simulate the dense two-phase turbulent up flow of Bingham fluid with particles. With the increasing of the yield stress, the velocities of Bingham and particle decrease near the pipe centre. Comparing the two-phase flow of Bingham-particle with that of liquid-particle, it is found the source term of yield stress has significant effect on flow.
基金Supported by the Special Funds for Major State Basic Research Projects, PRC(G1999-0222-08) and the National Natural Science Foundation of China(No. 19872039).
文摘A full second-order moment (FSM) model and an algebraic stress (ASM) two-phase turbulence modelare proposed and applied to predict turbulent bubble-liquid flows in a 2D rectangular bubble column. Predictiongives the bubble and liquid velocities, bubble volume fraction, bubble and liquid Reynolds stresses and bubble-liquidvelocity correlation. For predicted two-phase velocities and bubble volume fraction there is only slight differencebetween these two models, and the simulation results using both two models are in good agreement with the particleimage velocimetry (PIV) measurements. Although the predicted two-phase Reynolds stresses using the FSM are insomewhat better agreement with the PIV measurements than those predicted using the ASM, the Reynolds stressespredicted using both two models are in general agreement with the experiments. Therefore, it is suggested to usethe ASM two-phase turbulence model in engineering application for saving the computation time.
基金The project supported by the Special Funds for Major State Basic Research,China(G-1999-0222-08)the Postdoctoral Science Foundation(2004036239)
文摘A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concepts of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision and through a unified treatment of these two kinds of fluctuations. The proposed model is used to simulate gas-particle flows in a channel and in a downer. Simulation results are in agreement with the experimental results reported in references and are near the results obtained using the sin- gle-scale second-order moment two-phase turbulence model superposed with a particle collision model (USM-θ model) in most regions.
基金Project supported by China Post-Doctoral Science Foundation(No.2004036239)
文摘A two-scale second-order moment two-phase turbulence model accounting for inter-particle collision is developed, based on the concept of particle large-scale fluctuation due to turbulence and particle small-scale fluctuation due to collision. The proposed model is used to simulate gas-particle downer reactor flows. The computational results of both particle volume fraction and mean velocity are in agreement with the experimental results. After analyzing effects of empirical coefficient on prediction results, we can come to a conclusion that, inside the limit range of empirical coefficient, the predictions do not reveal a large sensitivity to the empirical coefficient in the downer reactor, but a relatively great change of the constants has important effect on the prediction.
基金supported by the National Key Project of Fundamental Research of China (Grant No.G1999-0222-07-08)the National Natural Science Foundation of China (Grant Nos.50736006 and 50606026)the Foundation of State Key Laboratory of Engines,Tianjin University (Grant No.K2010-07)
文摘Turbulent dispersed multiphase flows,including gas-particle,gas-droplet and bubble-liquid flows,are widely encountered in various engineering facilities.Modeling of two-phase turbulence,in particular the dispersed phase turbulence,is the key problem in the Eulerian-Eulerian simulation of practical dispersed multiphase flows.Although different models were developed and used,the experimental validation shows that they cannot always give satisfactory prediction results.In this paper the present author give a detailed review of the unified second-order moment (USM),k-k p and nonlinear k-k p two-phase turbulence models,proposed by him.The derivation and closure of these models are described in detail and the experimental validation and application of these models are extensively discussed.
文摘This work presents a comprehensive second-order predictive modeling (PM) methodology based on the maximum entropy (MaxEnt) principle for obtaining best-estimate mean values and correlations for model responses and parameters. This methodology is designated by the acronym 2<sup>nd</sup>-BERRU-PMP, where the attribute “2<sup>nd</sup>” indicates that this methodology incorporates second- order uncertainties (means and covariances) and second (and higher) order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best-Estimate Results with Reduced Uncertainties” and the last letter (“P”) in the acronym indicates “probabilistic,” referring to the MaxEnt probabilistic inclusion of the computational model responses. This is in contradistinction to the 2<sup>nd</sup>-BERRU-PMD methodology, which deterministically combines the computed model responses with the experimental information, as presented in the accompanying work (Part I). Although both the 2<sup>nd</sup>-BERRU-PMP and the 2<sup>nd</sup>-BERRU-PMD methodologies yield expressions that include second (and higher) order sensitivities of responses to model parameters, the respective expressions for the predicted responses, for the calibrated predicted parameters and for their predicted uncertainties (covariances), are not identical to each other. Nevertheless, the results predicted by both the 2<sup>nd</sup>-BERRU-PMP and the 2<sup>nd</sup>-BERRU-PMD methodologies encompass, as particular cases, the results produced by the extant data assimilation and data adjustment procedures, which rely on the minimization, in a least-square sense, of a user-defined functional meant to represent the discrepancies between measured and computed model responses.
文摘This work presents a comprehensive second-order predictive modeling (PM) methodology designated by the acronym 2<sup>nd</sup>-BERRU-PMD. The attribute “2<sup>nd</sup>” indicates that this methodology incorporates second-order uncertainties (means and covariances) and second-order sensitivities of computed model responses to model parameters. The acronym BERRU stands for “Best- Estimate Results with Reduced Uncertainties” and the last letter (“D”) in the acronym indicates “deterministic,” referring to the deterministic inclusion of the computational model responses. The 2<sup>nd</sup>-BERRU-PMD methodology is fundamentally based on the maximum entropy (MaxEnt) principle. This principle is in contradistinction to the fundamental principle that underlies the extant data assimilation and/or adjustment procedures which minimize in a least-square sense a subjective user-defined functional which is meant to represent the discrepancies between measured and computed model responses. It is shown that the 2<sup>nd</sup>-BERRU-PMD methodology generalizes and extends current data assimilation and/or data adjustment procedures while overcoming the fundamental limitations of these procedures. In the accompanying work (Part II), the alternative framework for developing the “second- order MaxEnt predictive modelling methodology” is presented by incorporating probabilistically (as opposed to “deterministically”) the computed model responses.
基金supported by the National Key Project of Fundamental Research of China (Grant No. G1999-0222-07-08)the Projects of the National Natural Science Foundation of China (Grant Nos. 50736006 and 50606026)the Foundation of the State Key Laboratory of Engines, Tianjin University (Grant No. K-2010-07)
文摘Dense gas-particle flows are frequently encountered in fluidized beds,riser and downer reactors,pneumatic transport and the near-wall zone of dilute gas-particle flows.Particle-particle collision plays an important role in the behavior of two-phase flows.In this paper a USM-Q two-phase turbulence model for dense gas-particle flows is proposed to account for both two-phase turbulence and inter-particle collision.For two-fluid large-eddy simulation of gas-particle flows,the author proposed a unified second-order moment(USM) two-phase SGS stress model and a two-phase k-kp SGS energy-equation stress model.The proposed models can fully account for the interaction between the gas and particle SGS stresses.
文摘This work illustrates the innovative results obtained by applying the recently developed the 2<sup>nd</sup>-order predictive modeling methodology called “2<sup>nd</sup>- BERRU-PM”, where the acronym BERRU denotes “best-estimate results with reduced uncertainties” and “PM” denotes “predictive modeling.” The physical system selected for this illustrative application is a polyethylene-reflected plutonium (acronym: PERP) OECD/NEA reactor physics benchmark. This benchmark is modeled using the neutron transport Boltzmann equation (involving 21,976 uncertain parameters), the solution of which is representative of “large-scale computations.” The results obtained in this work confirm the fact that the 2<sup>nd</sup>-BERRU-PM methodology predicts best-estimate results that fall in between the corresponding computed and measured values, while reducing the predicted standard deviations of the predicted results to values smaller than either the experimentally measured or the computed values of the respective standard deviations. The obtained results also indicate that 2<sup>nd</sup>-order response sensitivities must always be included to quantify the need for including (or not) the 3<sup>rd</sup>- and/or 4<sup>th</sup>-order sensitivities. When the parameters are known with high precision, the contributions of the higher-order sensitivities diminish with increasing order, so that the inclusion of the 1<sup>st</sup>- and 2<sup>nd</sup>-order sensitivities may suffice for obtaining accurate predicted best- estimate response values and best-estimate standard deviations. On the other hand, when the parameters’ standard deviations are sufficiently large to approach (or be outside of) the radius of convergence of the multivariate Taylor-series which represents the response in the phase-space of model parameters, the contributions stemming from the 3<sup>rd</sup>- and even 4<sup>th</sup>-order sensitivities are necessary to ensure consistency between the computed and measured response. In such cases, the use of only the 1<sup>st</sup>-order sensitivities erroneously indicates that the computed results are inconsistent with the respective measured response. Ongoing research aims at extending the 2<sup>nd</sup>-BERRU-PM methodology to fourth-order, thus enabling the computation of third-order response correlations (skewness) and fourth-order response correlations (kurtosis).
基金Project supported by the National Natural Science Foundation of China (No. 12272323)。
文摘Second-order axially moving systems are common models in the field of dynamics, such as axially moving strings, cables, and belts. In the traditional research work, it is difficult to obtain closed-form solutions for the forced vibration when the damping effect and the coupling effect of multiple second-order models are considered.In this paper, Green's function method based on the Laplace transform is used to obtain closed-form solutions for the forced vibration of second-order axially moving systems. By taking the axially moving damping string system and multi-string system connected by springs as examples, the detailed solution methods and the analytical Green's functions of these second-order systems are given. The mode functions and frequency equations are also obtained by the obtained Green's functions. The reliability and convenience of the results are verified by several examples. This paper provides a systematic analytical method for the dynamic analysis of second-order axially moving systems, and the obtained Green's functions are applicable to different second-order systems rather than just string systems. In addition, the work of this paper also has positive significance for the study on the forced vibration of high-order systems.
基金funded by the National Natural Science Foundation of China(No.51809135)the Shandong Provincial Natural Science Foundation(No.ZR2018BEE 047)+1 种基金the National Natural Science Foundation of China–Shandong Joint Fund(No.U2006229)the SKL of HESS(No.HESS-1808).
文摘The selection of wave force models will significantly impact the structural responses of floating wind turbines.In this study,comparisons of wave force model effects on the structural responses and fatigue loads of a semi-submersible floating wind turbine(SFWT)were conducted.Simulations were performed by employing the Morison equation(ME)with linear or second-order wave kinematics and potential flow theory(PFT)with first-or second-order wave forces.A comparison of regular waves,irregular waves,and coupled wind/waves analyses with the experimental data showed that many of the simulation results and experimental data are relatively consistent.However,notable discrepancies are found in the response amplitude operators for platform heave,tower base bending moment,and tension in mooring lines.PFT models give more satisfactory results of heave but more significant discrepan-cies in tower base bending moment than the ME models.In irregular wave analyses,low-frequency resonances were captured by PFT models with second-order difference-frequency terms,and high-frequency resonances were captured by the ME models or PFT models with second-order sum-frequency terms.These force models capture the response frequencies but do not reasonably predict the response amplitudes.The coupled wind/waves analyses showed more satisfactory results than the wave-only analyses.However,an important detail to note is that this satisfactory result is based on the overprediction of wind-induced responses.
基金supported by the National Natural Science Foundation of China (10872192)
文摘With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples:turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder (GL) model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.
文摘We propose that the QCD vacuum pion tetrahedron condensate density vary in space and drops to extremely low values in the Kennan, Barger and Cowie (KBC) void in analogy to earth’s atmospheric density drop with elevation from earth. We propose a formula for the gravitation acceleration based on the non-uniform pion tetrahedron condensate. Gravity may be due to the underlying microscopic attraction between quarks and antiquarks, which are part of the vacuum pion tetrahedron condensate. We propose an electron tetrahedron model, where electrons are comprised of tetraquark tetrahedrons, and . The quarks determine the negative electron charge and the or quarks determine the electron two spin states. The electron tetrahedron may perform a high frequency quark exchange reactions with the pion tetrahedron condensate by tunneling through the condensation gap creating a delocalized electron cloud with a fixed spin. The pion tetrahedron may act as a QCD glue bonding electron pairs in atoms and molecules and protons to neutrons in the nuclei. Conservation of valence quarks and antiquarks is proposed.
基金Project supported by the National Natural Science Foundation of China(Nos.11971303 and 11871330)。
文摘In this paper,the static output feedback stabilization for large-scale unstable second-order singular systems is investigated.First,the upper bound of all unstable eigenvalues of second-order singular systems is derived.Then,by using the argument principle,a computable stability criterion is proposed to check the stability of secondorder singular systems.Furthermore,by applying model reduction methods to original systems,a static output feedback design algorithm for stabilizing second-order singular systems is presented.A simulation example is provided to illustrate the effectiveness of the design algorithm.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10975181 and 11175247)the National Basic Research Program of China(Grant No.2012CB922104)
文摘We apply the second-order Born-Oppenheimer (BO) approximation to investigate the dynamics of the Rabi model, which describes the interaction between a two-level system and a single bosonic mode beyond the rotating wave approxi- mation. By comparing with the numerical results, we find that our approach works well when the frequency of the two-level system is much smaller than that of the bosonic mode.
基金supported by the National Natural Science Foundation of China(Nos.11372130,11290153,and 11290154)
文摘In this paper, an attitude maneuver control problem is investigated for a rigid spacecraft using an array of two variable speed control moment gyroscopes (VSCMGs) with gimbal axes skewed to each other. A mathematical model is constructed by taking the spacecraft and the gyroscopes together as an integrated system, with the coupling interaction between them considered. To overcome the singular issues of the VSCMGs due to the conventional torque-based method, the first-order derivative of gimbal rates and the second-order derivative of the rotor spinning velocity, instead of the gyroscope torques, are taken as input variables. Moreover, taking external disturbances into account, a feedback control law is designed for the system based on a method of nonlinear model predictive control (NMPC). The attitude maneuver can be realized fast and smoothly by using the proposed controller in this paper.
基金Project supported by the National Natural Science Foundation of China(Grant No.11272051)
文摘In this paper, we investigate the solution moment stability for a Harrison-type predator-prey model with parametric dichotomous noises. Using the Shapiro-Loginov formula, the equations for the first-order and second-order moments are obtained and the corresponding stable conditions are given. It is found that the solution moment stability depends on the noise intensity and correlation time of noise. The first-order and second-order moments become unstable with the decrease of correlation time. That is, the dichotomous noise can improve the solution moment stability with respect to Gaussian white noise. Finally, some numerical results are presented to verify the theoretical analyses.
文摘This work presents the “Second-Order Comprehensive Adjoint Sensitivity Analysis Methodology (2<sup>nd</sup>-CASAM)” for the efficient and exact computation of 1<sup>st</sup>- and 2<sup>nd</sup>-order response sensitivities to uncertain parameters and domain boundaries of linear systems. The model’s response (<em>i.e.</em>, model result of interest) is a generic nonlinear function of the model’s forward and adjoint state functions, and also depends on the imprecisely known boundaries and model parameters. In the practically important particular case when the response is a scalar-valued functional of the forward and adjoint state functions characterizing a model comprising N parameters, the 2<sup>nd</sup>-CASAM requires a single large-scale computation using the First-Level Adjoint Sensitivity System (1<sup>st</sup>-LASS) for obtaining all of the first-order response sensitivities, and at most N large-scale computations using the Second-Level Adjoint Sensitivity System (2<sup>nd</sup>-LASS) for obtaining exactly all of the second-order response sensitivities. In contradistinction, forward other methods would require (<em>N</em>2/2 + 3 <em>N</em>/2) large-scale computations for obtaining all of the first- and second-order sensitivities. This work also shows that constructing and solving the 2<sup>nd</sup>-LASS requires very little additional effort beyond the construction of the 1<sup>st</sup>-LASS needed for computing the first-order sensitivities. Solving the equations underlying the 1<sup>st</sup>-LASS and 2<sup>nd</sup>-LASS requires the same computational solvers as needed for solving (<em>i.e.</em>, “inverting”) either the forward or the adjoint linear operators underlying the initial model. Therefore, the same computer software and “solvers” used for solving the original system of equations can also be used for solving the 1<sup>st</sup>-LASS and the 2<sup>nd</sup>-LASS. Since neither the 1<sup>st</sup>-LASS nor the 2<sup>nd</sup>-LASS involves any differentials of the operators underlying the original system, the 1<sup>st</sup>-LASS is designated as a “<u>first-level</u>” (as opposed to a “first-order”) adjoint sensitivity system, while the 2<sup>nd</sup>-LASS is designated as a “<u>second-level</u>” (rather than a “second-order”) adjoint sensitivity system. Mixed second-order response sensitivities involving boundary parameters may arise from all source terms of the 2<sup>nd</sup>-LASS that involve the imprecisely known boundary parameters. Notably, the 2<sup>nd</sup>-LASS encompasses an automatic, inherent, and independent “solution verification” mechanism of the correctness and accuracy of the 2nd-level adjoint functions needed for the efficient and exact computation of the second-order sensitivities.
基金supported by the National Natural Science Funds of China(Grant No.41205071)the Ministry of Science and Technology of China for the National Basic Research Program of China(973 Program:Grant No.2011CB309704)the funding support from the U.S.Department of Energy(DOE),Office of Science,Earth System Modeling Program
文摘A two-moment bulk stratiform microphysics scheme, including recently developed physically-based droplet activation/ice nucleation parameterizations has been implemented into the Grid-point Atmospheric Model of IAP LASG (GAMIL) as an effort to enhance the model's capability to simulate aerosol indirect effects. Unlike the previous one-moment cloud microphysics scheme, the new scheme produces a reasonable rep- resentation of cloud particle size and number concentration. This scheme captures the observed spatial variations in cloud droplet number concentrations. Simulated ice crystal number concentrations in cirrus clouds qualitatively agree with in situ observations. The longwave and shortwave cloud forcings are in better agreement with observations. Sensitivity tests show that the column cloud droplet number concentrations calculated from two different droplet activation parameterizations are similar. However, ice crystal number concentration in mixed-phased clouds is sensitive to different heterogeneous ice nucleation formulations. The simulation with high ice crystal number concentration in mixed-phase clouds has less liquid water path and weaker cloud forcing. ~rthermore, ice crystal number concentration in cirrus clouds is sensitive to different ice nucleation parameterizations. Sensitivity tests also suggest that the impact of pre-existing ice crystals on homogeneous freezing in old clouds should be taken into account.