Monte-Carlo simulation of a stochastic differential equation

For solving higher dimensional diffusion equations with an inhomogeneous diffusion coefficient,Monte Carlo（MC） techniques are considered to be more effective than other algorithms, such as finite element method or finite difference method. The inhomogeneity of diffusion coefficient strongly limits the use of different numerical techniques. For better convergence, methods with higher orders have been kept forward to allow MC codes with large step size. The main focus of this work is to look for operators that can produce converging results for large step sizes. As a first step, our comparative analysis has been applied to a general stochastic problem.Subsequently, our formulization is applied to the problem of pitch angle scattering resulting from Coulomb collisions of charge particles in the toroidal devices.

NUMERICAL SIMULATION OF THE FLOW PAST WING-BODY USING EULER EQUATIONS

AbstractThis paper gives H-O type grid generation and applies finite-volume method forwing-body flow.A 3-D Euler code for wing-body has been developed. We improve in thiscode the distribution of the number of grid points over wing and body surfaces,and keepthe training edge of the wing being one of the grid lines even when the edge has sweepbackor sweep forward.The code developed can provide not only the spanwise pressure distribu-tion as usual ou cross-flow planes,but also the chordwise Pressure distributions.The com-putation for NASA TN D-712 wing-body model with the present code shows that thecomputed pressure distributions are in very good agreement with the experiment.

BIFURCATIONS AND NEW EXACT TRAVELLING WAVE SOLUTIONS OF THE COUPLED NONLINEAR SCHRDINGER-KdV EQUATIONS

By using the method of dynamical system, the exact travelling wave solutions of the coupled nonlinear Schrdinger-KdV equations are studied. Based on this method, all phase portraits of the system in the parametric space are given. All possible bounded travelling wave solutions such as solitary wave solutions and periodic travelling wave solutions are obtained. With the aid of Maple software, the numerical simulations are conducted for solitary wave solutions and periodic travelling wave solutions to the coupled nonlinear Schrdinger-KdV equations. The results show that the presented findings improve the related previous conclusions.

Numerical simulation of three-dimensional breaking waves and its interaction with a vertical circular cylinder

Wave breaking plays an important role in wave-structure interaction. A novel control volume finite element method with adaptive unstructured meshes is employed here to study 3-D breaking waves. The numerical framework consists of a ＂volume of fluid＂ type method for the interface capturing and adaptive unstructured meshes to improve computational efficiency. The numerical model is validated against experimental measurements of breaking wave over a sloping beach and is then used to study the breaking wave impact on a vertical circular cylinder on a slope. Detailed complex interfacial structures during wave impact, such as plunging jet formation and splash-up are captured in the simulation, demonstrating the capability of the present method.

Nematic and smectic stripe phases and stripe-Sk X transformations

Based on the findings of skyrmion nature of stripes and the metastability of a state of an arbitrary number of skyrmions,precisely controlled manipulation of stripes of skyrmion number 1 in pre-designed structures and mutual transformation between helical states and skyrmion crystals(Sk Xs)are demonstrated in chiral magnetic films.As a proof of the concept,we show how to use patterned magnetic fields and spin-transfer torques(STTs)to generate nematic and smectic stripe phases,as well as“UST”mosaic from three curved stripes.Cutting one stripe into many pieces and coalescing several skyrmions into one by various external fields are good ways to transform helical states and Sk Xs from each other.

Mechanisms and Kinetics of Reaction CHCIBr＋NO2

The reaction CHClBr＋NO2 was investigated via quantum chemical methods and kinetic calculations. The reaction mechanism on the singlet potential energy surface（PES） was considered by B3LYP method, and the energies were calculated at the CCSD（T） and CASPT2 levels of theory. The rate constants and the ratios of products were obtained by utilizing VTST and RRKM methods over wide temperature and pressure ranges. Our results indicate that carbon-to-nitrogen approach via a barrierless process is preferred in the initial association of CHClBr and NO2. The dominant product is BrNO＋CHCIO（PI）, which agrees well with the experimental observation. P2（ClNO＋CHBrO） and P3（HNO＋CBrClO） may also have minor contributions to the reaction. The calculated overall rate constants are independent of pressure and consistent with the experimental data, which can be fitted with the following equation over the temperature range of 200--1500 K： k（T）=2.31 × 10^-15T^0.99exp（771/T）. Compared with reaction CH2Br＋NO2, reaction CHCIBr＋NO2 has decreased the overall rate constants.

Beam dynamics study of RFQ for CADS with a 3D space-charge-effect

The ADS（accelerator driven subcritical system） project was proposed by the Chinese Academy of Sciences.The initial proton beams delivered from an electron cyclotron resonance ion source can be effectively accelerated by 162.5 MHz 4.2 m long room temperature radio-frequency-quadrupoles（RFQ） operating in CW mode.To test the feasibility of this physical design,a new Fortran code for RFQ beam dynamics study,which is space charge dominated,was developed.This program is based on Particle-In-Cell（PIC） technique in the time domain.Using the RFQ structure designed for the CADS project,the beam dynamics behavior is performed.The well-known simulation code TRACK is used for benchmarks.The results given by these two codes show good agreements.Numerical techniques as well as the results of beam dynamics studies are presented in this paper.