Pro/E在测定滑动轴承径向油膜压力中的应用

运用Pro/Engineer软件,作出滑动轴承径向油膜压力分布曲线图,并求出油膜曲面面积,从而得到了滑动轴承径向油膜压力平均值。传统手工作图将滑动轴承径向油膜压力分布曲线画在方格纸上,通过数格法得到近似曲面面积,误差大。通过Pro/Engineer软件作图简单、美观、精确,值得推广。

FLOW FIELD ANALYSES OF PLANE JET AT LOW REYNOLDS NUMBERS USING LATTICE BOLTZMANN METHOD

A two-dimensional（2-D） incompressible plane jet is investigated using the lattice Boltzmann method（LBM） for low Reynolds numbers of 42 and 65 based on the jet-exit-width and the maximum jet-exit-velocity. The results show that the mean centerline velocity decays as x-1/3 and the jet spreads as x2/3 in the self-similar region, which are consistent with the theoretical predictions and the experimental data. The time histories and PSD analyses of the instantaneous centerline velocities indicate the periodic behavior and the interaction between periodic components of velocities should not be neglected in the far field region, although it is invisible in the near field region.

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional （3D） DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method （AGMG）. We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner （AGMG-CG） to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling （3DDCXH）. In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.

Visualization of 3-D Field of Explosion

In order to visualize the 3-D field of explosion and describe the complex physical phenomena of explosion, the 3-D data resulting from numerical simulation by 3-D multi-material in cell (MMIC), and the application of volume visualization is explored, based on the characteristics of explosion and shock. Based on this, a visualization system for 3-D explosion--ViSC3D is designed. Approaches for the visualization of 3-D field of explosion are presented. The algorithm and the functions of ViSC3D are also presented. ViSC3D is thus a useful tool to observe and analyze either the full picture or the details of a 3-D field of explosion, that are difficult to observe and analyze directly. With ViSC3D, the field of explosion between the hill slopes is visualized. The cutaway views and 2-D slices are also given. The full picture and partial details of 3-D field of explosion can be observed clearly. Furthermore, ViSC3D can be used to visualize other similar 3-D data fields.

Singularity-free Green's function for EM sources embedded in a stratified medium

We present a method to unify the calculation of Green＇s functions for an electromagnetic（EM） transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green＇s function are derived by decomposing the partial wave solutions to Helmholtz＇s equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method（MCSEM） demonstrates its simplicity and flexibility.

Finite element numerical simulation of 2.5D direct current method based on mesh refinement and recoarsement

To deal with the problem of low computational precision at the nodes near the source and satisfy the requirements for computational efficiency in inversion imaging and finite-element numerical simulations of the direct current method, we propose a new mesh refinement and recoarsement method for a two-dimensional point source. We introduce the mesh refinement and mesh recoarsement into the traditional structured mesh subdivision. By refining the horizontal grids, the singularity owing to the point source is minimized and the topography is simulated. By recoarsening the horizontal grids, the number of grid cells is reduced significantly and computational efficiency is improved. Model tests show that the proposed method solves the singularity problem and reduces the number of grid cells by 80% compared to the uniform grid refinement.

Research on optimization of valve open time of the launch barge's ballast tanks

Launch barge is an effective tool for transporting ship segments from one place to another in shipyards. During shifting of segments onto a barge, the slideway on the barge＇s deck must be adjusted to maintain the same level as the wharf and also the barge must be kept level by adjusting the water in the ballast tanks. When to open the adjusting valves is an important factor influencing the barge＇s trim during the water-adjustment process. Because these adjustments are complex a mathematical model was formulated,after analyzing the characteristics of the process of moving the segments onto the barges deck, and considering the effects of this movement＇s speed and variations in tidal levels during the move. Then the model was solved by the penalty function method, the grid method, and improved simulated annealing, respectively. The best optimization model and its corresponding solution were then determined. Finally, it was proven that the model and the method adopted are correct and suitable, by calculating and analysing an example.

New Approach to Linear Stackelberg Problems with Multiple Leaders-followers

In this paper, we study linear static Stac kelberg problems with multiple leaders-followers in which each decision maker wi thin his group may or may not cooperate. An exact penalty function method is dev eloped. The duality gaps of the followers’ problems are appended to the leaders’ objective function with a penalty. The structure leads to the decomposition of the composite problem into a series of linear programmings leading to an efficie nt algorithm. We prove that local optimality is reached for an exact penalty fun ction and illustrate the method with three examples. The model in this paper ext ends the stackelberg leader-follower model.

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.

Accelerated Matrix Recovery via Random Projection Based on Inexact Augmented Lagrange Multiplier Method

In this paper, a unified matrix recovery model was proposed for diverse corrupted matrices. Resulting from the separable structure of the proposed model, the convex optimization problem can be solved efficiently by adopting an inexact augmented Lagrange multiplier (IALM) method. Additionally, a random projection accelerated technique (IALM+RP) was adopted to improve the success rate. From the preliminary numerical comparisons, it was indicated that for the standard robust principal component analysis (PCA) problem, IALM+RP was at least two to six times faster than IALM with an insignificant reduction in accuracy; and for the outlier pursuit (OP) problem, IALM+RP was at least 6.9 times faster, even up to 8.3 times faster when the size of matrix was 2 000×2 000.

Gridding cropland data reconstruction over the agricultural region of China in 1820

Recent studies have demonstrated the importance of LUCC change with climate and ecosystem simulation, but the result could only be determined precisely if a high-resolution underlying land cover map is used. While the efforts based satellites have provided a good baseline for present land cover, what the next advancement in the research about LUCC change required is the development of reconstruction of historical LUCC change especially spatially-explicit historical dataset. Being different from other similar studies, this study is based on the analysis of historical land use patterns in the traditional cultivated region of China. Taking no account of the less important factors, altitude, slope and population patterns are selected as the major drivers of reclamation in ancient China, and used to design the HCGM （Historical Cropland Gridding Model, at a 60 km×60 km resolution）, which is an empirical model for allocating the historical cropland inventory data spatially to grid cells in each political unit. Then we use this model to reconstruct cropland distribution of the study area in 1820, and verify the result by prefectural cropland data of 1820, which is from the historical documents. The statistical analyzing result shows that the model can simulate the patterns of the cropland distribution in the historical period in the traditional cultivated region efficiently.

A Numerical Method for Calculating Transmission Coefficients Across Arbitrary Potential Barriers with High Accuracy

We report a new method for calculating transmission coefficients across arbitrary potential barriers based on the Runge-Kutta method. A numerical solution of the Schrodinger equation is calculated using the Runge-Kutta method,and a new model is established to analyze the numerical results to find the transmission coefficient. This technique is applied to various cases, such as parabolic potential barrier and double-barrier structures. Transmission probability with high precision is obtained and discussed. The tunnelling current density through a MOS structure is also explored and the result coincides with the Fowler-Nordheim model,which indicates the applicability of our method.

Numerical modeling of effect of slot on bubble motion in aluminum electrolytic process

A transient three-dimensional（3 D） model was established to understand the bubble motion in an industrial electrolytic process. An anode with a new design was tested. It incorporates two slots that allow an efficient removal of gas bubbles. The electromagnetic fields were described by solving Maxwell＇s equations. The bubble movement was studied with two-way coupling Euler-Lagrange approach. The interplay of current density and bubble nucleation rate was included. The collision and coalescence of bubbles were considered. Random walk module was invoked for involving the chaotic effect of the turbulence. The numerical results were validated by experimental measurements. The results indicate that the current distribution and the bubble nucleation periodically change. Due to the slot, the bubble elimination heavily increases. The contribution of the slot to the bubble removal exceeds 50% in the case of three currents, and the promotion of the slot decays with increasing the current.

Heat Generation by Electric Current in Normal-Metal-Molecular Quantum Dot-Superconductor System

We investigate the heat generation induced by electrical current in a normal-metal-molecular quantum dot-superconductor （NDS） system. By using nonequilibrium Green＇s function method, the heat generation Q is derived and studied in detail. The superconducting lead influences the heat generation significantly. An obvious step appears in Q - eV characteristics and the iocation of this step is related with the phonon frequency ωo. The heat generations exhibit very different behaviour in the condition eV 〈 △ and eV 〉 △ due to different tunneling mechanism. From the study of Q - eVg curves, there is an extra peak as eV 〉 △. The difference in this two cases is also shown in Q - ωo curve, an extra peak emerges as eV 〉 △.

NUMERICAL SIMULATIONS OF EVOLUTIONS FROM WALL PULSE TO TURBULENT COHE RENT STRUCTURES

The high order compact d if ference method is developed for solving the perturbation equations based on Navi er Stokes equations, and is used in studying complex evolution processes from w all negative pulse to the turbulent coherent structure in the channel flow. Th is method contains three dimensional coupling difference scheme with high accur acy and high resolution, and the high order time splitting methods. Compared with the general spectral method, the method can be used to research turbule nt coherent structure under more general boundary conditions and in flow domains . In this paper, the generation and evolution of the turbulent coherent structur es ind uced by wall pulse in the channel flow are simulated, and the basic characterist ics and rules of the turbulent coherent structure are shown. Computational r esults indicate that a wall negative pulse is more convenient than the resonant three wave model.

The Numerical Modeling of 3-D Elastic Wave Equation Using a High-Order,Staggered-Grid, Finite Difference Scheme

This article provides the application of the high-order, staggered-grid, finite-difference scheme to model elastic wave propagation in 3-D isotropic media. Here, we use second-order, tempo- ral- and high-order spatial finite-difference formulations with a staggered grid for discretization of the 3-D elastic wave equations of motion. The set of absorbing boundary conditions based on paraxial approximations of 3-D elastic wave equations are applied to the numerical boundaries. The trial re- sults for the salt model show that the numerical dispersion is decreased to a minimum extent, the accuracy high and diffracted waves abundant. It also shows that this method can be used for modeling wave propagation in complex media with the lateral variation of velocity.

The Interaction of Oblique Flexural Gravity Waves With a Small Bottom Deformation on a Porous Ocean-Bed： Green＇s Function Approach

The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is assumed to be covered by an infinitely extended thin uniform elastic plate, while the lower surface is bounded by a porous bottom surface having a small deformation. By employing a simplified perturbation analysis, involving a small parameter c^（〈〈l ）, which measures the smallness of the deformation, the governing Boundary Value Problem （BVP） is reduced to a simpler BVP for the first-order correction of the potential function. This BVP is solved using a method based on Green＇s integral theorem with the introduction of suitable Green＇s function to obtain the first-order potential, and this potential function is then utilized to calculate the first-order reflection and transmission coefficients in terms of integrals involving the shape function c（x） representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that when the quotient of twice the component of the incident field wave number propagating just below the elastic plate and the ripple wave number approaches one, the theory predicts a resonant interaction between the bed and the surface below the elastic plate. Again, for small angles of incidence, the reflected wave energy is more as compared to the other angles of incidence. It is also observed that the reflected wave energy is somewhat sensitive to the changes in the flexural rigidity of the elastic plate, the porosity of the bed and the ripple wave numbers. The main advantage of the present study is that the results for the values of reflection and transmission coefficients obtained are found to satisfy the energy-balance relation almost accurately.

A Free Surface Frequency Domain Green Function with Viscous Dissipation and Partial Reflections from Side Walls

The free-surface Green function method is widely used in solving the radiation or diffraction problems caused by a ship or ocean structure oscillating on the waves. In the context of inviscid potential flow, hydrodynamic problems such as multi-body interaction and tank side wall effect cannot be properly dealt with based on the traditional free-surface frequency domain Green function method, in which the water viscosity is omitted and the energy dissipation effect is absent. In this paper, an open-sea Green function with viscous dissipation was presented within the theory ofvisco-potential flow. Then the tank Green function with a partial reflection from the side walls in wave tanks was formulated as a formal sum of open-sea Green functions representing the infinite images between two parallel side walls of the source in the tank. The new far-field characteristics of the tank Green function is vitally important fur improving the validity of side-wall effects evaluation, which can be used in supervising the tank model tests.

An Algorithm for Extracting Contour Lines Based on Interval Tree from Grid DEM

This paper proposes a new algorithm for determining the starting points of contour lines. The new algorithm is based on the interval tree. The result improves the algorithm＇s efficiency remarkably. Further, a new strategy is designed to constrain the direction of threading and the resulting contour bears more meaningful information.

Frequency in Middle of Magnon Band Gap in a Layered Ferromagnetic Superlattice

The frequency in middle of magnon energy band in a five-layer ferromagnetic superlattice is studied by using the linear spin-wave approach and Green＇s function technique. It is found that four energy gaps and corresponding four frequencie in middle of energy gaps exist in the magnon band along Kx direction perpendicular to the superlattice plane. The spin quantum numbers and the interlayer exchange couplings all affect the four frequencies in middle of the energy gaps. When all interlayer exchange couplings are same, the effect of spin quantum numbers on the frequency wg1 in middle of the energy gap Δw12 is complicated, and the frequency wg1 depends on the match of spin quantum numbers in each layer. Meanwhile, the frequencies wg2, wg3, and wg4 in middle of other energy gaps increase monotonously with increasing spin quantum numbers. When the spin quantum numbers in each layer are same, the frequencies wg1, wg2, wg3, and wg4 all increase monotonously with increasing interlayer exchange couplings.