A novel method for automated grid generation of ice shapes for local‐flow analysis

Modelling a complex geometry, such as ice roughness, plays a key role for the computational flow analysis over rough surfaces. This paper presents two enhancement ideas in modelling roughness geometry for local flow analysis over an aerodynamic surface. The first enhancement is use of the leading‐edge region of an airfoil as a perturbation to the parabola surface. The reasons for using a parabola as the base geometry are: it resembles the airfoil leading edge in the vicinity of its apex and it allows the use of a lower apparent Reynolds number. The second enhancement makes use of the Fourier analysis for modelling complex ice roughness on the leading edge of airfoils. This method of modelling provides an analytical expression, which describes the roughness geometry and the corresponding derivatives. The factors affecting the performance of the Fourier analysis were also investigated. It was shown that the number of sine–cosine terms and the number of control points are of importance. Finally, these enhancements are incorporated into an automated grid generation method over the airfoil ice accretion surface. The validations for both enhancements demonstrate that they can improve the current capability of grid generation and computational flow field analysis around airfoils with ice roughness. Copyright © 2004 John Wiley & Sons, Ltd.     

Least‐squares finite element method on adaptive grid for PDEs with shocks

We apply the least‐squares finite element method with adaptive grid to nonlinear time‐dependent PDEs with shocks. The least‐squares finite element method is also used in applying the deformation method to generate the adaptive moving grids. The effectiveness of this method is demonstrated by solving a Burgers' equation with shocks. Computational results on uniform grids and adaptive grids are compared for the purpose of evaluation. The results show that the adaptive grids can capture the shock more sharply with significantly less computational time. For moving shock, the adaptive grid moves correctly with the shock. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Meshless simulation of equilibrium swelling/deswelling of pH‐sensitive hydrogels

Hydrogels have been widely used in microelectromechanical systems (MEMS) and Bio‐MEMS devices. In this article, the equilibrium swelling/deswelling of the pH‐stimulus cylindrical hydrogel in the microchannel is studied and simulated by the meshless method. The multi‐field coupling model, called multi‐effect‐coupling pH‐stimulus (MECpH) model, is presented and used to describe the chemical field, electric field, and the mechanical field involved in the problem. The partial differential equations (PDEs) describing these three fields are either nonlinear or coupled together. This multi‐field coupling and high nonlinear characteristics produce difficulties for the conventional numerical methods (e.g., the finite element method or the finite difference method), so an alternative—meshless method is developed to discretize the PDEs, and the efficient iteration technique is adopted to solve the nonlinear problem. The computational results for the swelling/deswelling diameter of the hydrogel under the different pH values are firstly compared with experimental results, and they have a good agreement. The influences of other parameters on the mechanical properties of the hydrogel are also investigated in detail. It is shown that the multi‐field coupling model and the developed meshless method are efficient, stable, and accurate for simulation of the properties of the stimuli‐sensitive hydrogel. © 2005 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 326–337, 2006     

Electrophoresis simulations using Chebyshev pseudo-spectral method on a moving mesh

We present the implementation and demonstration of the Chebyshev pseudo-spectral method coupled with an adaptive mesh method for performing fast and highly accurate electrophoresis simulations. The Chebyshev pseudo-spectral method offers higher numerical accuracy than all other finite difference methods and is applicable for simulating all electrophoresis techniques in channels with open or closed boundaries. To improve the computational efficiency, we use a novel moving mesh scheme that clusters the grid points in the regions with poor numerical resolution. We demonstrate the application of the Chebyshev pseudo-spectral method on a moving mesh for simulating nonlinear electrophoretic processes through examples of isotachophoresis (ITP), isoelectric focusing (IEF), and electromigration-dispersion in capillary zone electrophoresis (CZE) at current densities as high as 1000 A/m. We also show the efficacy of our moving mesh method over existing methods that cluster the grid points in the regions with large concentration gradients. We have integrated the adaptive Chebyshev pseudo-spectral method in the open-source SPYCE simulator and verified its implementation with other electrophoresis simulators.     

改进型单晶炉提拉系统偏心平衡研究

随着光伏行业的快速发展, 对硅单晶的品质和长晶装备的稳定性的要求也不断提高。直拉法是生产硅单晶的主要方法,通过提高单晶炉副室的高度以扩大单晶硅的生产规模。由于副室高度的大幅增加,且单晶炉提拉头质心相对于旋转轴心有一定距离,对单晶炉整体稳定性有较大影响,从而降低了单晶硅的生产质量。针对此问题,对单晶炉建立可靠的力学分析模型,采用数值仿真方法,对单晶炉整体进行动力学响应分析,计算得到副室高度增加后的单晶炉工作时中钨丝绳下端晶棒的运动规律以及最大摆动幅度,为改进设计提供依据。数值仿真分析表明提高单晶炉副室高度后,提拉头较大的质心偏心是单晶炉提拉系统发生摆动的主要原因。在此基础上提出在提拉头上添加质心调节装置,通过控制系统调节可保证提拉头质心位置在旋转轴线上以降低提拉系统的摆动。     

On Physically Unacceptable Numerical Solutions Yielding Strong Chaotic Signals

Physically unacceptable chaotic numerical solutions of nonlinear circuits and systems are discussed in this paper. First, as an introduction, a simple example of a wrong choice of a numerical solver to deal with a second-order linear ordinary differential equation is presented. Then, the main result follows with the analysis of an ill-designed numerical approach to solve and analyze a particular nonlinear memristive circuit. The obtained trajectory of the numerical solution is unphysical (not acceptable), as it violates the presence of an invariant plane in the continuous systems. Such a poor outcome is then turned around, as we look at the unphysical numerical solution as a source of strong chaotic sequences. The 0–1 test for chaos and bifurcation diagrams are applied to prove that the unacceptable (from the continuous system point of view) numerical solutions are, in fact, useful chaotic sequences with possible applications in cryptography and the secure transmission of data.     

Orbital-Free Density Functional Theory for Molecular Structure Calculations

We give here an overview of the orbital-free density functional theory that is used for modeling atoms and molecules.We review typical approximations to the kinetic energy,exchange-correlation corrections to the kinetic and Hartree energies, and constructions of the pseudopotentials.We discuss numerical discretizations for the orbital-free methods and include several numerical results for illustrations.     

Numerical Study on Entropy Generation of the Multi-Stage Centrifugal Pump

The energy loss of the multi-stage centrifugal pump was investigated by numerical analysis using the entropy generation method with the RNG k-ε turbulence model. Entropy generation due to time-averaged motion and velocity fluctuation was mainly considered. It was found that the entropy generation of guide vanes and impellers account for 71.2% and 23.3% of the total entropy generation under the designed flow condition. The guide vanes are the main hydraulic loss domains and their entropy generation is about 9 W/K, followed by impellers. There are vortices at the tongue of the guide vane inlet as well as flow separations in the impellers, which lead to entropy generation. The fluid impacts the outer surface of the guide vanes, resulting in the increase in entropy generation. There are refluxes near the guide vane tongues which also increase the entropy generation of this part. The entropy generation distribution of the guide vanes and impellers was investigated, which found that the positive guide vane has more entropy generation compared with the reverse guide. The entropy generation of the blade suction surface is higher compared with the pressure surface. This study indicated that the entropy generation method has distinct advantages in the assessment of hydraulic loss.     

耦合FE/WB法在声分析中的应用

简要描叙FE法(finite element method)和WB法(wave based method)的理论背景以及耦合FE/WB法的数学基础.耦合FE/WB法利用两者的优势——FE法的广泛应用和WB法的高收敛特性,将FE模型中较大且几何简单的部分采用WB法代替.耦合模型具有相对较少的自由度.对于较高的频率还可以进行细分得到更高的计算精度,并利用模态缩减法进一步减少自由度数.数值算例结果表明,该耦合方法有潜力覆盖中频段的声分析.     

The stability analysis of rolling motion of hypersonic vehicles and its validations

The stability of the rolling motion of near space hypersonic vehicles with rudder control is studied using method of qualitative analysis of nonlinear differential equations, and the stability criteria of the deflected rolling motions are improved. The outcomes can serve as the basis for further study regarding the influence of pitching and lateral motion on the stability of rolling motion. To validate the theoretical results, numerical simulations were done for the rolling motion of two hypersonic vehicles with typical configurations. Also, wind tunnel experiments for four aircraft models with typical configurations have been done. The results show that: 1) there exist two dynamic patterns of the rolling motion under statically stable condition. The first one is point attractor, for which the motion of aircraft returns to the original state. The second is periodic attractor, for which the aircraft rolls periodically. 2) Under statically unstable condition, there exist three dynamic patterns of rolling motion, namely, the point attractor, periodic attractor around deflected state of rolling motion, and double periodic attractors or chaotic attractors.