A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a si...A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.展开更多
Effective temperature level of stream, namely stream pseudo temperature, is determined by its actual temperature and heat transfer temperature difference contribution value. Heat transfer temperature difference con-tr...Effective temperature level of stream, namely stream pseudo temperature, is determined by its actual temperature and heat transfer temperature difference contribution value. Heat transfer temperature difference con-tribution value of a stream depends on its heat transfer film coefficient, cost per unit heat transfer area, actual tem-perature, and so on. In the determination of the suitable heat transfer temperature difference contribution values of the stream, the total annual cost of multistream heat exchanger network (MSHEN) is regarded as an objective func-tion, and genetic/simulated annealing algorithm (GA/SA) is adopted for optimizing the heat transfer temperature difference contribution values of the stream. The stream pseudo temperatures are subsequently obtained. On the ba-sis of stream pseudo temperature, optimized MSHEN can be attained by the temperature-enthalpy (T-H) diagram method. This approach is characterized with fewer decision variables and higher feasibility of solutions. The calcu-lation efficiency of GA/SA can be remarkably enhanced by this approach and more probability is shown in search-ing the global optimum solution. Hence this approach is presented for solving industrial-sized MSHEN which is difficult to deal by traditional algorithm. Moreover, in the optimization of stream heat transfer temperature differ-ence contribution values, the effects of the stream temperature, the heat transfer film coefficient, and the construc-tion material of heat exchangers are considered, therefore this approach can be used to optimize and design heat exchanger network (HEN) with unequal heat transfer film coefficients and different of construction materials. The performance of the proposed approach has been demonstrated with three examples and the obtained solutions are compared with those available in literatures. The results show that the large-scale MSHEN synthesis problems can be solved to obtain good solutions with the modest computational effort.展开更多
Propagation properties of spatially pseudo-partially coherent Gaussian Schell-model beams through the atmo- spheric turbulence over a long-distance uplink path are studied by numerical simulation. A linear coordinatio...Propagation properties of spatially pseudo-partially coherent Gaussian Schell-model beams through the atmo- spheric turbulence over a long-distance uplink path are studied by numerical simulation. A linear coordination trans- formation is introduced to overcome the window effect and the loss-of-resolution problem. The beam spreading, beam wandering, and intensity scintillation as functions of turbulence strength, source correlation length, and change fre- quency of random phase that models the partial coherence of the source are analyzed. It is found that the beam spreading and the intensity scintillation of the partially coherent beam are less affected by the turbulence than those of the coherent one, but it suffers from a more severe diffractive effect, and the change frequency of random phase has no evident influence on it. The beam wandering is insensitive to the variation of source correlation length, and decreases firstly then goes to a fixed value as the change frequency increases.展开更多
概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generali...概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generalized probability density evolution equation,GDEE).对于GDEE的求解通常采用有限差分法,然而,由于GDEE是初始条件间断的变系数一阶双曲偏微分方程,通过有限差分法求解GDEE可能会面临网格敏感性问题、数值色散和数值耗散现象.文章从全局逼近的角度出发,基于Chebyshev拟谱法为GDEE构造了全局插值格式,解决了数值色散、数值耗散以及网格敏感性问题.考虑GDEE的系数在每个时间步长均为常数,推导了GDEE在每一个时间步长内时域上的序列矩阵指数解.由于序列矩阵指数解形式上是解析的,从而很好地克服了数值稳定性问题.两个数值算例表明,通过Chebyshev拟谱法结合时域的序列矩阵指数解求解GDEE得到的结果与精确解以及Monte Carlo模拟的结果非常吻合,且数值耗散和数值色散现象几乎可以忽略.此外,拟谱法具有高效的收敛性且序列矩阵指数解不受CFL (Courant-Friedrichs-Lewy)条件的限制,因此该方法具有良好的数值稳定性和计算效率.展开更多
The method of studying crystals using divergent X-ray beams emitted by a point source is known as the divergent beam method. If the X-ray source is inside a crystal and the characteristic radiation is emitted by the c...The method of studying crystals using divergent X-ray beams emitted by a point source is known as the divergent beam method. If the X-ray source is inside a crystal and the characteristic radiation is emitted by the crystal itself, it is called the Kossel method. If the source is outside the crystal and the divergent展开更多
文摘A supersonic turbulent flowfield involving the pseudo-shock waves in an isolator of a supersonic combustion ramjet is computed using two different LES codes which are a high-order upwind finite volume scheme, and a sixth order compact differencing scheme utilizing the localized artificial diffusivity method for stabilizing shock waves and employing a wall model to enable the use of coarse mesh. In the validation study where a supersonic turbulent boundary layer flow over a flat plate is examined, both LES codes are well validated using velocity profile in the boundary layer given by the hot-wire anemometry and normal stress given by the laser Doppler anemometry. In particular, the sixth order compact differencing scheme gives closer agreements with these experimental data. Then, the validated LES codes are applied to solve the Mach number 2.5 supersonic turbulent flowfield involving the pseudo-shock waves. It is shown that typical features of unsteady flowfield of the pseudo-shock waves are well obtained by both schemes. Again, it is indicated that the sixth order compact differencing scheme gives closer agreements with the existing velocity data obtained by particle image velocimetry and pressure fluctuation data on the wall surface. Besides, the computational cost of the compact differencing scheme is found to be 1/7 of that for the upwind finite volume scheme, even though a wall model is solved at each grid point on the wall surface. Therefore, the obtained results in the present study allow recommending the sixth order compact differencing scheme with a wall model for simulating supersonic turbulent flowfield in an isolator involving the pseudo-shock waves.
基金Acknowledgment： The work is supported by the National Natural Science Foundation of China #60675006 and the National Key Technology R＆D Program during the l lth Five-Year Plan Period #2006BAI03A09.
基金Supported by the Deutsche Forschungsgemeinschaft (DFG No.RO 294/9).
文摘Effective temperature level of stream, namely stream pseudo temperature, is determined by its actual temperature and heat transfer temperature difference contribution value. Heat transfer temperature difference con-tribution value of a stream depends on its heat transfer film coefficient, cost per unit heat transfer area, actual tem-perature, and so on. In the determination of the suitable heat transfer temperature difference contribution values of the stream, the total annual cost of multistream heat exchanger network (MSHEN) is regarded as an objective func-tion, and genetic/simulated annealing algorithm (GA/SA) is adopted for optimizing the heat transfer temperature difference contribution values of the stream. The stream pseudo temperatures are subsequently obtained. On the ba-sis of stream pseudo temperature, optimized MSHEN can be attained by the temperature-enthalpy (T-H) diagram method. This approach is characterized with fewer decision variables and higher feasibility of solutions. The calcu-lation efficiency of GA/SA can be remarkably enhanced by this approach and more probability is shown in search-ing the global optimum solution. Hence this approach is presented for solving industrial-sized MSHEN which is difficult to deal by traditional algorithm. Moreover, in the optimization of stream heat transfer temperature differ-ence contribution values, the effects of the stream temperature, the heat transfer film coefficient, and the construc-tion material of heat exchangers are considered, therefore this approach can be used to optimize and design heat exchanger network (HEN) with unequal heat transfer film coefficients and different of construction materials. The performance of the proposed approach has been demonstrated with three examples and the obtained solutions are compared with those available in literatures. The results show that the large-scale MSHEN synthesis problems can be solved to obtain good solutions with the modest computational effort.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 61107066 and 40805006)
文摘Propagation properties of spatially pseudo-partially coherent Gaussian Schell-model beams through the atmo- spheric turbulence over a long-distance uplink path are studied by numerical simulation. A linear coordination trans- formation is introduced to overcome the window effect and the loss-of-resolution problem. The beam spreading, beam wandering, and intensity scintillation as functions of turbulence strength, source correlation length, and change fre- quency of random phase that models the partial coherence of the source are analyzed. It is found that the beam spreading and the intensity scintillation of the partially coherent beam are less affected by the turbulence than those of the coherent one, but it suffers from a more severe diffractive effect, and the change frequency of random phase has no evident influence on it. The beam wandering is insensitive to the variation of source correlation length, and decreases firstly then goes to a fixed value as the change frequency increases.
文摘概率密度演化方法(probability density evolution equation,PDEM)为非线性随机结构的动力响应分析提供了新的途径.通过PDEM获得结构响应概率密度函数(probability density function,PDF)的关键步骤是求解广义概率密度演化方程(generalized probability density evolution equation,GDEE).对于GDEE的求解通常采用有限差分法,然而,由于GDEE是初始条件间断的变系数一阶双曲偏微分方程,通过有限差分法求解GDEE可能会面临网格敏感性问题、数值色散和数值耗散现象.文章从全局逼近的角度出发,基于Chebyshev拟谱法为GDEE构造了全局插值格式,解决了数值色散、数值耗散以及网格敏感性问题.考虑GDEE的系数在每个时间步长均为常数,推导了GDEE在每一个时间步长内时域上的序列矩阵指数解.由于序列矩阵指数解形式上是解析的,从而很好地克服了数值稳定性问题.两个数值算例表明,通过Chebyshev拟谱法结合时域的序列矩阵指数解求解GDEE得到的结果与精确解以及Monte Carlo模拟的结果非常吻合,且数值耗散和数值色散现象几乎可以忽略.此外,拟谱法具有高效的收敛性且序列矩阵指数解不受CFL (Courant-Friedrichs-Lewy)条件的限制,因此该方法具有良好的数值稳定性和计算效率.
文摘The method of studying crystals using divergent X-ray beams emitted by a point source is known as the divergent beam method. If the X-ray source is inside a crystal and the characteristic radiation is emitted by the crystal itself, it is called the Kossel method. If the source is outside the crystal and the divergent
