Compared with front engine vehicle, the windward side’s flow field in cooling model of rear engine bus is complicated and it can’t be calculated by means of 1D model. For this problem, this paper has used Star-CCM t...Compared with front engine vehicle, the windward side’s flow field in cooling model of rear engine bus is complicated and it can’t be calculated by means of 1D model. For this problem, this paper has used Star-CCM to build a 3D simulation model of cooling system, engine compartment and complete vehicle. Then, it had a 1D/3D coupling calculation on cooling system with Kuli software. It could be helpful in the optimization design of the flow field of rear engine compartment and optimization match of cooling system.展开更多
By means of both a theory for pressure-induced shifts (PS) of energy spectra and a theory for shifts of energy spectra due to electron-phonon interaction (EPI), the 'pure electronic' PS and the PS due to EPI o...By means of both a theory for pressure-induced shifts (PS) of energy spectra and a theory for shifts of energy spectra due to electron-phonon interaction (EPI), the 'pure electronic' PS and the PS due to EPI of R<SUB>1</SUB> line of GSGG:Cr<SUP>3+</SUP> at 70 K have been calculated, respectively. Their physical origins have been revealed. It is found that the admixture of and base-wavefunctions in the wavefunctions of R<SUB>1</SUB> level of GSGG:Cr<SUP>3+</SUP> at 70 K is remarkable under the normal pressure, and the degree of the admixture rapidly decreases with increasing pressure. The change of the degree of the admixture with the pressure plays a key role for not only the pure electronic PS of R<SUB>1</SUB> line but also the PS of R<SUB>1</SUB> line due to EPI. The detailed calculations and analyses show that the pressure-dependent behaviors of the pure electronic PS of R<SUB>1</SUB> line and the PS of R<SUB>1</SUB> line due to EPI are quite different. It is the combined effect of them that gives rise to the total PS of R<SUB>1</SUB> line, which has satisfactorily explained the experimental data (including a reversal of PS of R<SUB>1</SUB> line). In contributions to PS of R<SUB>1</SUB> line due to EPI at 70 K, the temperature-independent contribution is much larger than the temperature-dependent contribution. The former results from the interaction between the zero-point vibration of the lattice and localized electronic state.展开更多
By means of both a theory for pressure-induced shifts (PS) of energy spectra and a theory for shifts of energy spectra due to electron-phonon interaction (EPI), the 'pure electronic' PS and the PS due to EPI o...By means of both a theory for pressure-induced shifts (PS) of energy spectra and a theory for shifts of energy spectra due to electron-phonon interaction (EPI), the 'pure electronic' PS and the PS due to EPI of R<SUB>1</SUB> line, R<SUB>2</SUB> line, and U band of GSGG:Cr<SUP>3+</SUP> at 300 K have been calculated, respectively. The calculated results are in good agreement with all the experimental data. Their physical origins have also been explained. It is found that the mixing-degree of and base-wavefunctions in the wavefunctions of R<SUB>1</SUB> level of GSGG:Cr<SUP>3+</SUP> at 300 K is remarkable under normal pressure, and the mixing-degree rapidly decreases with increasing pressure. The change of the mixing-degree with pressure plays a key role not only for the 'pure electronic' PS of R<SUB>1</SUB> line and R<SUB>2</SUB> line but also the PS of R<SUB>1</SUB> line and R<SUB>2</SUB> line due to EPI. The pressure-dependent behaviors of the 'pure electronic' PS of R<SUB>1</SUB> line (or R<SUB>2</SUB> line) and the PS of R<SUB>1</SUB> line (or R<SUB>2</SUB> line) due to EPI are quite different. It is the combined effect of them that gives rise to the total PS of R<SUB>1</SUB> line (or R<SUB>2</SUB> line). In the range of about 15 kbar ~ 45 kbar, the mergence and/or order-reversal between levels and levels take place, which cause the fluctuation of the rate of PS for with pressure. At 300 K, both the temperature-dependent contribution to R<SUB>1</SUB> line (or R<SUB>2</SUB> line or U band) from EPI and the temperature-independent one are important.展开更多
Although the hydraulic transients in pipe systems are usually simulated by using a one-dimensional (l-D) approach, local three-dimensional (3-D) simulations are necessary because of obvious 3-D flow features in so...Although the hydraulic transients in pipe systems are usually simulated by using a one-dimensional (l-D) approach, local three-dimensional (3-D) simulations are necessary because of obvious 3-D flow features in some local regions of the hydropower systems. This paper combines the 1-D method with a 3-D fluid flow model to simulate the Multi-Dimensional (MD) hydraulic transients in hydropower systems and proposes two methods for modeling the compressible water with the correct wave speed, and two strategies for efficiently coupling the 1-D and 3-D computational domains. The methods are validated by simulating the water hammer waves and the oscillations of the water level in a surge tank, and comparing the results ~with the 1-D solution data. An MD study is conducted for the transient flows in a realistic water conveying system that consists of a draft tube, a tailrace surge tank and a tailrace tunnel. It is shown that the 1-D-3-D coupling approach is an efficient and promising way to simulate the hydraulic transients in the hydropower systems in which the interactions between 1-D hydraulic fluctuations of the pipeline systems and the local 3-D flow patterns should be considered.展开更多
The vent tube is commonly used for the water hammer protection in the hydropower tailrace system. In transient processes, with air entering and exiting the vent tube, one sees complex hydraulic phenomena, which threat...The vent tube is commonly used for the water hammer protection in the hydropower tailrace system. In transient processes, with air entering and exiting the vent tube, one sees complex hydraulic phenomena, which threaten the station's safe operation. It is necessary to investigate the transient mechanisms in the tailrace system with vent tube. In this paper, a 3-D, two-phase numerical model of a vent tube on the connection of the tailrace tunnel and the diversion tunnel, is developed based on the FLUENT with the volume of fluid(VOF) algorithm to investigate the transient air-water flow patterns and the complex hydraulic phenomena in the vent tube of the tailrace system. A 1-D and 3-D unidirectional adjacent coupling(1-D-3-D-UAC) approach with a linear interpolation method is adopted to adjust the timesteps between the 1-D model and the 3-D model on the tunnel inlet and outlet boundaries through the user defined function(UDF), to transmit the data from the 1-D model to the 3-D model. The model is verified by comparing the results obtained by using the 1-D model alone and from the experiments in literature. The transient flow processes under the full load rejection consist of four stages: the water level dropping stage, the air entering stage, the air pocket collapsing stage, and the air exiting stage. Detailed hydraulic phenomena in the air pocket collapsing process are also discussed.展开更多
This paper describes the calculation method for unsteady state conditions in the secondary air systems in gas turbines. The 1D-3D-Structure coupled method was applied. A 1D code was used to model the standard componen...This paper describes the calculation method for unsteady state conditions in the secondary air systems in gas turbines. The 1D-3D-Structure coupled method was applied. A 1D code was used to model the standard components that have typical geometric characteristics. Their flow and heat transfer were described by empirical correlations based on experimental data or CFD calculations. A 3D code was used to model the non-standard components that cannot be described by typical geometric languages, while a finite element analysis was carried out to compute the structural deformation and heat conduction at certain important positions. These codes were coupled through their interfaces. Thus, the changes in heat transfer and structure and their interactions caused by exterior disturbances can be reflected. The results of the coupling method in an unsteady state showed an apparent deviation from the existing data, while the results in the steady state were highly consistent with the existing data. The difference in the results in the unsteady state was caused primarily by structural deformation that cannot be predicted by the 1D method. Thus, in order to obtain the unsteady state performance of a secondary air system more accurately and efficiently, the 1D-3D-Structure coupled method should be used.展开更多
In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC...In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC)method has been proposed.It combines one-dimensional tran-sient calculation of air system,Conventional Sequence Staggered(CSS)method,Time-adaptive Aerothermal Coupling calculation(TAC)method and differential evolution optimization algorithm to obtain an efficient and high-precision aerothermal coupling calculation method of air system.Considering both the heat conduction in the solid domain and the flow in the fluid domain as unsteady states in the OTAC,the interaction of fluid-solid information within a single coupling time step size was implemented based on the CSS method.Furthermore,the coupling time step size was automatically adjusted with the number of iterations by using the Proportional-Integral-Deri vative(PID)controller.Results show that when compared with the traditional loosely coupling method with a fixed time step size,the computational accuracy and efficiency of the OTAC method are improved by 8.9%and 30%,respectively.Compared with the tight coupling calculation,the OTAC method can achieve a speedup of 1 to 2 orders of magnitude,while the calculation error is maintained within 6.1%.展开更多
The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the c...The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the calculation accuracy.In this paper,a Cross-dimensional Data Transmission method(CDT)from 3D to 1D is proposed by introducing flow field uniformity into the data transmission.First,a uniformity index was established to quantify the flow field parameter distribution characteristics,and a uniformity index prediction model based on the locally weighted regression method(Lowess)was established to quickly obtain the flow field information.Then,an information selection criterion in 3D to 1D data transmission was established based on the Spearman rank correlation coefficient between the uniformity index and the accuracy of coupling calculation,and the calculation method was automatically determined according to the established criterion.Finally,a modified function was obtained by fitting the ratio of the 3D mass-average parameters to the analytical solution,which are then used to modify the selected parameters at the 1D-3D interface.Taking a typical disk cavity air system as an example,the results show that the calculation accuracy of the CDT method is greatly improved by a relative 53.88%compared with the traditional 1D-3D coupling method.Furthermore,the CDT method achieves a speedup of 2 to 3 orders of magnitude compared to the 3D calculation.展开更多
文摘Compared with front engine vehicle, the windward side’s flow field in cooling model of rear engine bus is complicated and it can’t be calculated by means of 1D model. For this problem, this paper has used Star-CCM to build a 3D simulation model of cooling system, engine compartment and complete vehicle. Then, it had a 1D/3D coupling calculation on cooling system with Kuli software. It could be helpful in the optimization design of the flow field of rear engine compartment and optimization match of cooling system.
文摘By means of both a theory for pressure-induced shifts (PS) of energy spectra and a theory for shifts of energy spectra due to electron-phonon interaction (EPI), the 'pure electronic' PS and the PS due to EPI of R<SUB>1</SUB> line of GSGG:Cr<SUP>3+</SUP> at 70 K have been calculated, respectively. Their physical origins have been revealed. It is found that the admixture of and base-wavefunctions in the wavefunctions of R<SUB>1</SUB> level of GSGG:Cr<SUP>3+</SUP> at 70 K is remarkable under the normal pressure, and the degree of the admixture rapidly decreases with increasing pressure. The change of the degree of the admixture with the pressure plays a key role for not only the pure electronic PS of R<SUB>1</SUB> line but also the PS of R<SUB>1</SUB> line due to EPI. The detailed calculations and analyses show that the pressure-dependent behaviors of the pure electronic PS of R<SUB>1</SUB> line and the PS of R<SUB>1</SUB> line due to EPI are quite different. It is the combined effect of them that gives rise to the total PS of R<SUB>1</SUB> line, which has satisfactorily explained the experimental data (including a reversal of PS of R<SUB>1</SUB> line). In contributions to PS of R<SUB>1</SUB> line due to EPI at 70 K, the temperature-independent contribution is much larger than the temperature-dependent contribution. The former results from the interaction between the zero-point vibration of the lattice and localized electronic state.
文摘By means of both a theory for pressure-induced shifts (PS) of energy spectra and a theory for shifts of energy spectra due to electron-phonon interaction (EPI), the 'pure electronic' PS and the PS due to EPI of R<SUB>1</SUB> line, R<SUB>2</SUB> line, and U band of GSGG:Cr<SUP>3+</SUP> at 300 K have been calculated, respectively. The calculated results are in good agreement with all the experimental data. Their physical origins have also been explained. It is found that the mixing-degree of and base-wavefunctions in the wavefunctions of R<SUB>1</SUB> level of GSGG:Cr<SUP>3+</SUP> at 300 K is remarkable under normal pressure, and the mixing-degree rapidly decreases with increasing pressure. The change of the mixing-degree with pressure plays a key role not only for the 'pure electronic' PS of R<SUB>1</SUB> line and R<SUB>2</SUB> line but also the PS of R<SUB>1</SUB> line and R<SUB>2</SUB> line due to EPI. The pressure-dependent behaviors of the 'pure electronic' PS of R<SUB>1</SUB> line (or R<SUB>2</SUB> line) and the PS of R<SUB>1</SUB> line (or R<SUB>2</SUB> line) due to EPI are quite different. It is the combined effect of them that gives rise to the total PS of R<SUB>1</SUB> line (or R<SUB>2</SUB> line). In the range of about 15 kbar ~ 45 kbar, the mergence and/or order-reversal between levels and levels take place, which cause the fluctuation of the rate of PS for with pressure. At 300 K, both the temperature-dependent contribution to R<SUB>1</SUB> line (or R<SUB>2</SUB> line or U band) from EPI and the temperature-independent one are important.
基金the National Natural Science Foundation of China (Grant Nos. 51039005, 50909076)
文摘Although the hydraulic transients in pipe systems are usually simulated by using a one-dimensional (l-D) approach, local three-dimensional (3-D) simulations are necessary because of obvious 3-D flow features in some local regions of the hydropower systems. This paper combines the 1-D method with a 3-D fluid flow model to simulate the Multi-Dimensional (MD) hydraulic transients in hydropower systems and proposes two methods for modeling the compressible water with the correct wave speed, and two strategies for efficiently coupling the 1-D and 3-D computational domains. The methods are validated by simulating the water hammer waves and the oscillations of the water level in a surge tank, and comparing the results ~with the 1-D solution data. An MD study is conducted for the transient flows in a realistic water conveying system that consists of a draft tube, a tailrace surge tank and a tailrace tunnel. It is shown that the 1-D-3-D coupling approach is an efficient and promising way to simulate the hydraulic transients in the hydropower systems in which the interactions between 1-D hydraulic fluctuations of the pipeline systems and the local 3-D flow patterns should be considered.
基金Project supported by the National Key R&D Program of China(Grant No.2016YFC0401810)the Research Project of Science and Technology Commission of Shanghai Munici-pality(Grant No.16DZ1202205)and the Fundamental Research Funds for the Central Universities(Grant No.2016B10814)
文摘The vent tube is commonly used for the water hammer protection in the hydropower tailrace system. In transient processes, with air entering and exiting the vent tube, one sees complex hydraulic phenomena, which threaten the station's safe operation. It is necessary to investigate the transient mechanisms in the tailrace system with vent tube. In this paper, a 3-D, two-phase numerical model of a vent tube on the connection of the tailrace tunnel and the diversion tunnel, is developed based on the FLUENT with the volume of fluid(VOF) algorithm to investigate the transient air-water flow patterns and the complex hydraulic phenomena in the vent tube of the tailrace system. A 1-D and 3-D unidirectional adjacent coupling(1-D-3-D-UAC) approach with a linear interpolation method is adopted to adjust the timesteps between the 1-D model and the 3-D model on the tunnel inlet and outlet boundaries through the user defined function(UDF), to transmit the data from the 1-D model to the 3-D model. The model is verified by comparing the results obtained by using the 1-D model alone and from the experiments in literature. The transient flow processes under the full load rejection consist of four stages: the water level dropping stage, the air entering stage, the air pocket collapsing stage, and the air exiting stage. Detailed hydraulic phenomena in the air pocket collapsing process are also discussed.
基金supported by funds from National natural science foundation of China(Grant No.51176004)
文摘This paper describes the calculation method for unsteady state conditions in the secondary air systems in gas turbines. The 1D-3D-Structure coupled method was applied. A 1D code was used to model the standard components that have typical geometric characteristics. Their flow and heat transfer were described by empirical correlations based on experimental data or CFD calculations. A 3D code was used to model the non-standard components that cannot be described by typical geometric languages, while a finite element analysis was carried out to compute the structural deformation and heat conduction at certain important positions. These codes were coupled through their interfaces. Thus, the changes in heat transfer and structure and their interactions caused by exterior disturbances can be reflected. The results of the coupling method in an unsteady state showed an apparent deviation from the existing data, while the results in the steady state were highly consistent with the existing data. The difference in the results in the unsteady state was caused primarily by structural deformation that cannot be predicted by the 1D method. Thus, in order to obtain the unsteady state performance of a secondary air system more accurately and efficiently, the 1D-3D-Structure coupled method should be used.
基金support of the National Natural Science Foundation of China (No.52007002)the Science Center for Gas Turbine Project,China (No.P2022-A-II-007-001)the Fundamental Research Funds for the Central Universities,China (No.NS2023010).
文摘In order to further achieve the balance between the calculation accuracy and efficiency of the transient analysis of the aero-engine disc cavity system,an Optimized Time-adaptive Aerother-mal Coupling calculation(OTAC)method has been proposed.It combines one-dimensional tran-sient calculation of air system,Conventional Sequence Staggered(CSS)method,Time-adaptive Aerothermal Coupling calculation(TAC)method and differential evolution optimization algorithm to obtain an efficient and high-precision aerothermal coupling calculation method of air system.Considering both the heat conduction in the solid domain and the flow in the fluid domain as unsteady states in the OTAC,the interaction of fluid-solid information within a single coupling time step size was implemented based on the CSS method.Furthermore,the coupling time step size was automatically adjusted with the number of iterations by using the Proportional-Integral-Deri vative(PID)controller.Results show that when compared with the traditional loosely coupling method with a fixed time step size,the computational accuracy and efficiency of the OTAC method are improved by 8.9%and 30%,respectively.Compared with the tight coupling calculation,the OTAC method can achieve a speedup of 1 to 2 orders of magnitude,while the calculation error is maintained within 6.1%.
基金supported by the National Science and Technology Major Project,China(No.2017-III-0010-0036).
文摘The Secondary Air System(SAS)plays an important role in the safe operation and performance of aeroengines.The traditional 1D-3D coupling method loses information when used for secondary air systems,which affects the calculation accuracy.In this paper,a Cross-dimensional Data Transmission method(CDT)from 3D to 1D is proposed by introducing flow field uniformity into the data transmission.First,a uniformity index was established to quantify the flow field parameter distribution characteristics,and a uniformity index prediction model based on the locally weighted regression method(Lowess)was established to quickly obtain the flow field information.Then,an information selection criterion in 3D to 1D data transmission was established based on the Spearman rank correlation coefficient between the uniformity index and the accuracy of coupling calculation,and the calculation method was automatically determined according to the established criterion.Finally,a modified function was obtained by fitting the ratio of the 3D mass-average parameters to the analytical solution,which are then used to modify the selected parameters at the 1D-3D interface.Taking a typical disk cavity air system as an example,the results show that the calculation accuracy of the CDT method is greatly improved by a relative 53.88%compared with the traditional 1D-3D coupling method.Furthermore,the CDT method achieves a speedup of 2 to 3 orders of magnitude compared to the 3D calculation.