Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

Smoothing LDA Model for Text Categorization

Abstract. Latent Dirichlet Allocation (LDA) is a document level language model. In general, LDA employ the symmetry Dirichlet distribution as prior of the topic-words’ distributions to implement model smoothing. In this paper, we propose a data-driven smoothing strategy in which probability mass is allocated from smoothing-data to latent variables by the intrinsic inference procedure of LDA. In such a way, the arbitrariness of choosing latent variables'priors for the multi-level graphical model is overcome. Following this data-driven strategy,two concrete methods, Laplacian smoothing and Jelinek-Mercer smoothing, are employed to LDA model. Evaluations on different text categorization collections show data-driven smoothing can significantly improve the performance in balanced and unbalanced corpora.

Parallel Fully-Implicit Computation of Magnetohydrodynamics Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETS,: package. The results of convergence tests arc plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

中美大学生归因偏差的跨文化研究

This research systematically compared Chinese undergraduates with American undergraduates on four kinds of attributional bias: correspondent bias, siuational overattribution, intergroup attributional bias and self-serving attributional bias, and examined the effect of the implicit theories reflecting the cultures on attributional bias. First is analyzed three pairs of opposite implicit theories: dispositionalist theory and situationalist theory, generality and particularity, stressing the positive evaluation of self and despising the positive evaluation of self. It developed the Modern Implicit Theories Inventory and Traditional Implicit Theories Inventory to measure these implicit theories, and the results showed these inventories had satisfactory validity and reliability, and they were suitable for the group comparison of Chinese implicit theories with European-American. At the same time through the test it found Chinese undergraduates agreed all these opposite implicit theories more than American undergraduates. Second, it studied Chinese and American undergraduates' attributional accuracy on locus of causality. The results showed: Chinese and American undergraduates both had the correspondent bias under the different salient situational constraints, and the degree of Chinese and American undergraduates' correspondent bias under the different salient situational constraints had no significant difference' Chinese and American undergraduates both showed the situational overattribution; Chinese undergraduates had more the correspondent bias and situational overattribution than American undergraduates. Third, on the research of Chinese and American undergraduates' intergroup attributional bias, it found Chinese and American undergraduates both had no intergroup attributional bias among kin, friends and strangers, while they both show some favorable outcome effects for these three group actors. The favorable outcome effects were significant on the attributional dimensions of locus of causality and controllability for strangers' behavior, and stability for kin and friends' behavior rating by Chinese undergraduates, and stability for friends' behavior rating by American undergraduates. Fourth, it explored Chinese and American undergraduates' self-serving attributional bias, and the result indicated that Chinese and American undergraduates both showed significant self-serving attributional bias: for outcome effects, Chinese undergraduates' self-serving attributional bias were reflected on the attributional dimensions of locus of causality, stability, controllability and globality, and American undergraduates were reflected on the attributional dimensions of locus of causality, stability and globality; for categorization effects, both Chinese and American undergraduates' self-serving attributional bias were reflected on attributional difference between self's negative behavior and others', but Chinese undergraduates were embodied on the attributional dimensions of locus of causality, stability and globality while American undergraduates were reflected on the attributional dimensions of stability and globality. It also found Chinese undergraduates had more self-serving attributional bias than American undergraduates. This was reflected on the attributional dimensions of locus of causality, stability and controllability for outcome effects, and for categorization effects, locus of causality, stability and globality rating for self and others' negative behavior. All studies indicated that Chinese and American undergraduates' implicit theories had no significant effects on all their four attributional bias. These findings' potentially important implications were discussed and the further research was suggested.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

A Multimoment Finite-Volume Shallow-Water Model On The Yin-Yang Overset Spherical Grid

A numerical model for shallow-water equations has been built and tested on the Yin-Yang overset spherical grid. A high-order multimoment finite-volume method is used for the spatial discretization in which two kinds of so-called moments of the physical field [i.e., the volume integrated average ( VIA) and the point value (PV)] are treated as the model variables and updated separately in time. In the present model, the PV is computed by the semi-implicit semi-Lagrangian formulation, whereas the VIA is predicted in time via a flux-based finite-volume method and is numerically conserved on each component grid. The concept of including an extra moment (i.e., the volume-integrated value) to enforce the numerical conservativeness provides a general methodology and applies to the existing semi-implicit semi-Lagrangian formulations. Based on both VIA and PV, the high-order interpolation reconstruction can only be done over a single grid cell, which then minimizes the overlapping zone between the Yin and Yang components and effectively reduces the numerical errors introduced in the interpolation required to communicate the data between the two components. The present model completely gets around the singularity and grid convergence in the polar regions of the conventional longitude-latitude grid. Being an issue demanding further investigation, the high-order interpolation across the overlapping region of the Yin-Yang grid in the current model does not rigorously guarantee the numerical conservativeness. Nevertheless, these numerical tests show that the global conservation error in the present model is negligibly small. The model has competitive accuracy and efficiency.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Wavelike patterns in one-dimensional coupled map lattices

We investigate the existence of wavelike solution for the logistic coupled map lattices for which the spatiotemporal periodic patterns can be predicted by a simple two-dimensional mapping. The existence of such wavelike solutions is proved by the implicit function theorem with constraints. We also examine the stabilities of these wave solutions under perturbations of uniform small deformation type. We show that in some specific cases these perturbations are completely general. The technique used in this paper is also applicable to investigate other space-time regular patterns.

Compressible laminar boundary-layer flows of a dusty gas over a semi-infinite flat-plate

The compressible laminar boundary-layer flows of a dilute gas-particle mixture over a semi-infinite flat plate are investigated analytically. The governing equations are presented in a general form where more reasonable relations for the two-phase interaction and the gas viscosity are included. The detailed flow structures of the gas and particle phases are given in three distinct regions : the large-slip region near the leading edge, the moderate-slip region and the small-slip region far downstream. The asymptotic solutions for the two limiting regions are obtained by using a seriesexpansion method. The finite-difference solutions along the whole length of the plate are obtained by using implicit four-point and six-point schemes. The results from these two methods are compared and very good agreement is achieved. The characteristic quantities of the boundary layer are calculated and the effects on the flow produced by the particles are discussed. It is found that in the case of laminar boundary-layer flows, the skin friction and wall heat-transfer are higher and the displacement thickness is lower than in the pure-gas case alone. The results indicate that the Stokes-interaction relation is reasonable qualitatively but not correct quantitatively and a relevant non-Stokes relation of the interaction between the two phases should be specified when the particle Reynolds number is higher than unity.

Speckle photography applied to the density field of a flame

are obtained by using implicit four-point and six-point schemes. The results from

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.

高寒草甸生态系统陆地生物圈模式研究及应用

改进了适合于高寒草甸生态系统的陆地生态圈模式，分析了模式中温度变化与水分运动分层的物理原因，说明了气候状况对地表面能量交换的影响，给出了净辐射和蒸散量新的计算方法，提出了有很差分计算中具有二阶精度的Euler隐式格式，介绍了中国科学院海北高寒草甸生态站的气候概况和野外观测情况。最后利用本模式对高寒草甸生态站地区的土壤-植被-大气间水热交换过程进行了数值模拟，模拟值与实测值吻合较好。

任意阶隐式指数时程差分多步法及其在非线性系统中的应用

对非线性系统提出了任意阶隐式指数时程差分多步法 ,实现了任意阶次指数时程差分预测 校正算法 .发展完善了指数时程差分法 .将新算法应用于非线性系统 ,取得了较好的效果 .数值结果表明隐式指数时程差分多步法很好地修正了显式指数时程差分多步法 ,隐式指数时程差分多步法是一种高精度、高效率的方法

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.