二维磁性系统和准一维磁性纳米管的内能

采用量子统计理论的多体格林函数法计算二维单离子各向异性海森伯铁磁体、反铁磁体以及单壁铁磁纳米管的内能,对比铁磁体和反铁磁体的结果.在相同的参量下,反铁磁能量总是低于相应的铁磁能量(相变点除外).由于反铁磁能量随温度上升的速度较铁磁能量快,当温度升高到居里点TC和奈尔点TN时(TC=TN),铁磁能量和反铁磁能量相等.横向关联效应对系统内能的影响较大,不能忽略.

COMPRESSIBLE FLOW SIMULATION AROUND AIRFOIL BASED ON LATTICE BOLTZMANN METHOD

The flow around airfoil NACA0012 enwrapped by the body-fitted grid is simulated by a coupled doubledistribution-function （DDF） lattice Boltzmann method （LBM） for the compressible Navier-Stokes equations. Firstly, the method is tested by simulating the low Reynolds number flow at Ma =0. 5,a=0. 0, Re=5 000. Then the simulation of flow around the airfoil is carried out at Ma：0. 5, 0. 85, 1.2; a=-0.05, 1.0, 0.0, respectively. And a better result is obtained by using a local refined grid. It reduces the error produced by the grid at Ma=0. 85. Though the inviscid boundary condition is used to avoid the problem of flow transition to turbulence at high Reynolds numbers, the pressure distribution obtained by the simulation agrees well with that of the experimental results. Thus, it proves the reliability of the method and shows its potential for the compressible flow simulation. The suecessful application to the flow around airfoil lays a foundation of the numerical simulation of turbulence.

NOVEL IMMERSED BOUNDARY-LATTICE BOLTZMANN METHOD BASED ON FEEDBACK LAW

The lattice Boltzmann method （LBM） and the immersed boundary method （IBM） are alternative, com- putational techniques for solving complex fluid dynamics systems, and can take the place of the Navier-Stokes（N- S） equation. This paper proposes a novel immersed boundary-lattice Boltzmann method （IB-LBM） based on the feedback law. The method uses the immersed boundary concept in the LBM framework to capture the coupling between a body with complex geometry and a uniform fluid, Then, the flows around a stationary circular cylinder and two circular cylinders in a side by side arrangement are simulated by using the method. Results are agreed well with the benchmark data, so, the capability of the method for complex geometry is demonstrated. Different from the conventional IB-LBM, which uses the Hook＇s law or the direct forcing method to compute the interae- tion force, the method uses the feedback law--the feedback of velocity field and displacement information to calculate the force, thus ensuring the method has advantages of easy implementation and full parallelism.

Robust Delaunay Tetrahedronal Meshing Coupled with Advancing Front Method

A full automatic tetrahedronal mesh generation method for arbitrary 3D domains is described. First, the classic Delaunay method is coupled with simplified advancing front technique (AFT) to obtain the boundary mesh. Then, advancing front high quality point placement is used to generate internal points with optimal positions and a Delaunay method is used to insert them efficiently. Finally, optimization procedures are used for mesh quality improvements. Several application examples are presented to demonstrate the robustness and efficiency of the proposed meshing scheme.

A Parallel Algorithm for Adaptive Local Refinement of Tetrahedral Meshes Using Bisection

Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp：//lsec. cc. ac. cn/phg/）, a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.

TVD SCHEME WITH CHEMICAL REACTION FLOW AND ITS APPLICATION IN COMBUSTION GAS JETTING FLOW

By conjugating features of combustion gas jetting flows of the solid-rocket and using mathematical methods, a numerical scheme is systematically derived based on Harten′s standard TVD scheme, which fits for the flow with high temperature, pressure and velocity. The rational calculation formula of pressure partial derivation is also given out. By using the chemical kinetics knowledge, problems of multi-component and finite rate chemical reaction contained in combustion gas of the rocket flow field are discussed. The method for solving the mass source term of chemical reaction is clarified. Taking 9 reaction equations with 12 components as an example and utilizing the established calculation program, the free jetting flow field of the rocket is simulated. Numerical results show the correctness of the numerical scheme.

Ship Hull Form Optimization by Evolutionary Algorithm in Order to Diminish the Drag

This study presents a numerical method for optimizing hull form in calm water with respect to total drag which contains a viscous drag and a wave drag. The ITTC 1957 model-ship correlation line was used to predict frictional drag and the corrected linearized thin-ship theory was employed to estimate the wave drag The evolution strategy （ES） which is a member of the evolutionary algorithms （EAs） family obtains an optimum hull form by considering some design constraints. Standard Wigley hull is considered as an initial hull in optimization procedures for two test cases and new hull forms were achieved at Froude numbers 0.24, 0.316 and 0.408. In one case the ES technique was ran for the initial hull form, where the main dimensions were fixed and the only variables were the hull offsets. In the other case in addition to hull offsets, the raain dimensions were considered as variables that are optimized simultaneously. The numerical results of optimization procedure demonstrate that the optimized hull forms yield a reduction in total drag.

Statistical Average of Spin Operators for Calculation of Three-Component Magnetization

When one wants to calculate all the three components of magnetization of Heisenberg model under random phase approximation, at least one of the components should be the solution of an ordinary differential equation. In this paper such an equation is established. It is argued that the general expressions of magnetization for any spin quantum number S suggested before are the solution of the ordinary differential equation.

Introduction and evaluation of a novel hybrid brattice for improved dust control in underground mining faces: A computational study

A proper control and management of dust dispersion is essential to ensure safe and productive underground working environment. Brattice installation to direct the flow from main shaft to the mining face was found to be the most effective method to disperse dust particle away from the mining face. However,it limits the movement and disturbs the flexibility of the mining fleets and operators at the tunnel. This study proposes a hybrid brattice system- a combination of a physical brattice together with suitable and flexible directed and located air curtains- to mitigate dust dispersion from the mining face and reduce dust concentration to a safe level for the working operators. A validated three-dimensional computational fluid dynamic model utilizing Eulerian–Lagrangian approach is employed to track the dispersion of dust particle. Several possible hybrid brattice scenarios are evaluated with the objective to improve dust management in underground mine. The results suggest that implementation of hybrid brattice is beneficial for the mining operation: up to three times lower dust concentration is achieved as compared to that of the physical brattice without air curtain.

香港公司是如何管理的

管理一个2万多人的上市公司,实在不是一件容易的事情,然而陆宗霖却自有独特的方法。1990年,陆氏以其“科学管理+精神管理”理论,获美国林肯大学颁授工商管理学荣誉博士;翌年,港陆国际集团亦因其理论而获得亚洲管理院颁发的“人力发展及管理大奖”。

Factors affecting accuracy of radial point interpolation meshfree method for 3-D solid mechanics

Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.

Challenges in Study of Single Particles and Particle Swarms

Numerical simulation of multiphase flows in processing equipment in industry with two-fluid models and Eulerian-Lagrangian approaches requires the constitutive equations describing the interactions between the dispersed phase of high concentration and the continuous phase. The status of research on the forces on dispersed solid and fluid particles is reviewed in this article. As compared with the knowledge on drag of single solid particles study on panicle swarms and on other forces is not sufficient to meet the demand of reliable and efficient numerical simulation of multiphase flows. Thus, thorough study on the panicle swarms becomes the key to accurate multi-scale simulation of multiphase flows. Besides, the development of efficient algorithm dealing with the non-uniformity on both equipment and mesoscopic scales is recognized as an important issue to be resolved. The research topics in the near future are suggested.

Rotational Flow in a Narrow Annular Gap Based on Lattice Boltzmann Method

For comprehensive characteristics of flow in a gas bearing,lattice Boltzmann method(LBM)is applied for study of the two-dimensional flow between two eccentric cylinders with the inner one rotating at a high speed.The flow pattern and circumferential pressure distribution are discussed based on critical issues such as eccentricity ranging from 0.2 to 0.9,clearance ratio varying from 0.005 to 0.01 and rotating speed in the range of 3×104—1.8×105 r/min.The analysis and discussion on the circumferential pressure distribution affirmed the quasilinear relation between the extremum pressure and rotating speed.Furthermore,a high eccentricity and small clearance ratio contributes most to the fluctuation of the circumferential pressure distribution.The flow pattern inside the channel exhibits separation vortex under a large eccentricity.The conclusions drawn in this work give rise to prediction of the flow pattern in the gas bearing which is beneficial for evaluating the performance of as well as instructing the design and development.

Theoretical investigation of micropolar fluid flow between two porous disks

The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparison between OHAM and numerical method shows that OHAM is an exact and high efficient method for solving these kinds of problems. The results are presented to study the velocity and rotation profiles for different physical parameters such as Reynolds number, vortex viscosity parameter, spin gradient viscosity and microinertia density parameter. As an important outcome, the magnitude of the microrotation increases with an increase in the values of injection velocity while it decreases by increasing the values of suction velocity.

SnO_2-based solid solutions for CH_4 deep oxidation: Quantifying the lattice capacity of SnO_2 using an X-ray diffraction extrapolation method

A series of SnO2‐based catalysts modified by Mn, Zr, Ti and Pb oxides with a Sn/M （M=Mn, Zr, Ti and Pb） molar ratio of 9/1 were prepared by a co‐precipitation method and used for CH4 and CO oxidation. The Mn3＋, Zr4＋, Ti4＋and Pb4＋cations are incorporated into the lattice of tetragonal rutile SnO2 to form a solid solution structure. As a consequence, the surface area and thermal stability of the catalysts are improved. Moreover, the oxygen species of the modified catalysts become easier to be reduced. Therefore, the oxidation activity over the catalysts was improved, except for the one modified by Pb oxide. Manganese oxide demonstrates the best promotional effects for SnO2. Using an X‐ray diffraction extrapolation method, the lattice capacity of SnO2 for Mn2O3 was 0.135 g Mn2O3/g SnO2, which indicates that to form stable solid solution, only 21%Sn4＋cations in the lattice can be maximally replaced by Mn3＋. If the amount of Mn3＋cations is over the capacity, Mn2O3 will be formed, which is not favorable for the activity of the catalysts. The Sn rich samples with only Sn‐Mn solid solution phase show higher activity than the ones with excess Mn2O3 species.

Comparison of Improved Meshless Interpolation Schemes for SPH Method and Accuracy Analysis

In the smoothed particle hydrodynamics (SPH) method, a meshless interpolation scheme is needed for the unknown function in order to discretize the governing equation.A particle approximation method has so far been used for this purpose.Traditional particle interpolation (TPI) is simple and easy to do, but its low accuracy has become an obstacle to its wider application.This can be seen in the cases of particle disorder arrangements and derivative calculations.There are many different methods to improve accuracy, with the moving least square (MLS) method one of the most important meshless interpolation methods.Unfortunately, it requires complex matrix computing and so is quite time-consuming.The authors developed a simpler scheme, called higher-order particle interpolation (HPI).This scheme can get more accurate derivatives than the MLS method, and its function value and derivatives can be obtained simultaneously.Although this scheme was developed for the SPH method, it has been found useful for other meshless methods.

Global optimization for ducted coaxial-rotors aircraft based on Kriging model and improved particle swarm optimization algorithm

To improve the operational efficiency of global optimization in engineering, Kriging model was established to simplify the mathematical model for calculations. Ducted coaxial-rotors aircraft was taken as an example and Fluent software was applied to the virtual prototype simulations. Through simulation sample points, the total lift of the ducted coaxial-rotors aircraft was obtained. The Kriging model was then constructed, and the function was fitted. Improved particle swarm optimization(PSO) was also utilized for the global optimization of the Kriging model of the ducted coaxial-rotors aircraft for the determination of optimized global coordinates. Finally, the optimized results were simulated by Fluent. The results show that the Kriging model and the improved PSO algorithm significantly improve the lift performance of ducted coaxial-rotors aircraft and computer operational efficiency.

Numerical Modeling of Local Scour Around Hydraulic Structure in Sandy Beds by Dynamic Mesh Method

Local scour, a non-negligible factor in hydraulic engineering, endangers the safety of hydraulic structures. In this work, a numerical model for simulating local scour was constructed, based on the open source code computational fluid dynamics model Open FOAM. We consider both the bedload and suspended load sediment transport in the scour model and adopt the dynamic mesh method to simulate the evolution of the bed elevation. We use the finite area method to project data between the three-dimensional flow model and the two-dimensional(2D) scour model. We also improved the 2D sand slide method and added it to the scour model to correct the bed bathymetry when the bed slope angle exceeds the angle of repose. Moreover, to validate our scour model, we conducted and compared the results of three experiments with those of the developed model. The validation results show that our developed model can reliably simulate local scour.

On the Plaintiff's Subject Qualification for Environmental Public Interest Litigation Instituted by People's Procuratorates

Although it is not the main duty for the People's Procuratorate to institute environmental public interest litigation proceedings, the existence of "reserve" prosecution right of plaintiff not only is beneficial to timely prevention and cure of environmental damage, but also contributes indispensable judicial supervision force to the environmental legality construction. In order to safeguard environmental public interests, many local prosecutorial organizations have practiced environmental public interests litigation. Since Chinese law has not made explicit stipulation on the function of the People's Procuratorate instituting environmental public interest litigation, the plaintiff's subject qualification for environmental public interest litigation instituted by the People's Procuratorate is still in question. In order to propel the People's Procuratorate to completely and fully exercise their power of supervision and prosecution, there shall be corresponding legislation that provides powerful guarantee for the People's Procuratorate instituting environmental public interest litigation.

Numerical simulation on melt flow with bubble stirring and temperature field in aluminum holding furnace

The numerical model for predicting the flow and temperature fields of the melt in holding furnace with porous brick purging system were set up using Euler-Lagrange approach.In this model,bubbles coalescence and disintegration were ignored based on the dimensionless analysis,and the bubble size was assumed to be obedient to Rosin-Rammler distribution with a mean size of 0.6 mm.The results show that on reference operating condition,during the heating and agitation process,melt mixes well in the furnace,and the melt velocity increases with the increase of gas flux.Holding the melt for 30 min causes the max temperature in the bulk melt to increase to 60 K.After holding the heat,the agitation processing restarts,and it takes 10 min for the stratified melt to retrieve the homogeneous temperature field when the gas flux is 10 L/min,which shows deficient alloying and degassing in the melt.With the increase of gas flux from 10 to 20,30 and 40 L/min,the necessary recovery time decreases from 10 to 6,5 and 4 min gradually,which shows the improvement of the stirring efficiency.Depending on the processing purposes,for both good degassing performance and gas saving,proper operating strategy and parameters (gas flux,primarily) could be adjusted.