Comparative study on phase transition behaviors of fractional molecular field theory and random-site Ising model

Fractional molecular field theory(FMFT)is a phenomenological theory that describes phase transitions in crystals with randomly distributed components,such as the relaxor-ferroelectrics and spin glasses.In order to verify the feasibility of this theory,this paper fits it to the Monte Carlo simulations of specific heat and susceptibility versus temperature of two-dimensional(2D)random-site Ising model(2D-RSIM).The results indicate that the FMFT deviates from the 2D-RSIM significantly.The main reason for the deviation is that the 2D-RSIM is a typical system of component random distribution,where the real order parameter is spatially heterogeneous and has no symmetry of space translation,but the basic assumption of FMFT means that the parameter is spatially uniform and has symmetry of space translation.

Simulation of two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Global Two-Dimensional Chemistry Model and Simulation of Atmospheric Chemical Composition

A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.

Universality Class of the Nonequilibrium Phase Transition in Two-Dimensional Ising Ferromagnet Driven by Propagating Magnetic Field Wave

The purpose of this work is to identify the universality class of the nonequilibrium phase transition in the two-dimensional kinetic Ising ferromagnet driven by propagating magnetic field wave. To address this issue, the finite size analysis of the nonequilibrium phase transition, in two-dimensional Ising ferromagnet driven by plane propagating magnetic wave, is studied by Monte Carlo simulation. It is observed that the system undergoes a nonequilibrium dynamic phase transition from a high temperature dynamically symmetric (propagating) phase to a low temperature dynamically symmetry-broken (pinned) phase as the system is cooled below the transition temperature. This transition temperature is determined precisely by studying the fourth-order Binder Cumulant of the dynamic order parameter as a function of temperature for different system sizes (L). From the finite size analysis of dynamic order parameter ?and the dynamic susceptibility , we have estimated the critical exponents and ?(measured from the data read at the critical temperature obtained from Binder cumulant), and (measured from the peak positions of dynamic susceptibility). Our results indicate that such driven Ising ferromagnet belongs to the same universality class of the two-dimensional equilibrium Ising ferromagnet (where and ), within the limits of statistical errors.

Impact of Large Fresh Water Discharges into Hakata Bay Due to Torrential Rain Using Nested Two-Dimensional Convective-Dispersive Model

In this study, numerical simulation of a two-dimensional convective-dispersive model in Hakata Bay, Japan, is performed to analyze the impact of major river discharges due to torrential rain in Fukuoka City. Tank models are applied to calculate river discharges, which are taken into consideration as river inflow in the hydrodynamic model of Hakata Bay. A two-way nesting “edge” technique is developed and applied in the model in order to consider the influence of narrow and complex geographical features. The area around “Island City” and Imazu Bay are calculated in high resolution. The resulting model has high reproducibility since the calculated river discharges, tidal current, and salinity show good agreement with observed data. To analyze the impact of large river discharges, the calculation period is set from 11 September 2002 to 21 September 2002 since there was torrential rain on September 16 in the given year in Fukuoka City (163.5 mm/d). The results show that low-salinity water covered the whole of the inner part of Hakata Bay, and water of lower salinity than outer sea water (<34.0 psu) spread out to the bay’s mouth two days after the torrential rain event. Fresh water covered the entire area of Imazu Bay and flowed out from the mouth of the Bay after the torrential rain event. The behavior of fresh water after a few days of torrential rain was remarkably different from normal discharge river flow. These results indicate that the environment in Imazu Bay can be degraded severely by torrential rain. Therefore, countermeasures to protect ecosystems in Hakata Bay must be examined immediately.

Fast forward modeling of gravity anomalies for two-dimensional bodies of arbitrary shape and density distribution

A fast and high precision spatial domain algorithm is presented for forward modeling of two-dimensional(2D)body gravity anomalies of arbitrary shape and density distribution.The new algorithm takes advantage of the convolution properties of the expression for 2D gravity anomalies,uses a rectangular cell as a grid subdivision unit,and then 2D bodies with irregular cross-sections are approximated by a combination of 2D bodies with a rectangular cross section.The closed-form expression is used to calculate the gravitational anomalies of the combination of 2D bodies with a rectangular cross section.To improve computing effi ciency,the new algorithm uses a fast algorithm for the implementation of the Toeplitz matrix and vector multiplication.The synthetic 2D models with rectangular and circular cross-sections and constant and variable densities are designed to evaluate the computational accuracy and speed of the new algorithm.The experiment results show that the computation costs less than 6 s for a grid subdivision with 10000×10000 elements.Compared to the traditional forward modeling methods,the proposed method significantly improved computational effi ciency while guaranteeing computational accuracy.

The Model for Linear Magnetoresistance of Two-Dimensional Metal-Semiconductor Composites with Interfacial Shells

A metal-semiconductor composite with the interracial shells is investigated theoretically for the large linear mag- netoresistance effect of high doping Ag2＋δ Se and Ag2＋δ te materials. The magnetoresistance （MR） of composites is a function of the magnetic field, temperature, the conductivities of two phases without magnetic field, and the thickness and conductivity of the interracial shells. The MR increases with the increase of the magnetic field and with the decrease of temperature, and no saturation is found even under the high magnetic field. Moreover, it is interestingly found that the interracial shell is an important factor for the MR of the composites. The MR increases with the thickness and the conductivity of the interfacial shells. Lastly, the theoretical results on the MR are compared with the experimental data. It is found that the value of the MR of the composite with the interfacial shell is larger than that without the interfacial shell.

Nucleation of Kinetic Ising Model Under Oscillating Field

We have studied the nucleation process of a two-dimensional kinetic Ising model subject to a bias oscillating external field, focusing on how the nucleation time depends on the oscillation frequency. It is found that the nucleation time shows a clear-cut minimum with the variation of oscillation frequency, wherein the average size of the critical nuclei is the smallest, indicating that an oscillating external field with an optimal frequency can be much more favorable to the nucleation process than a constant field. We have also investigated the effect of the initial phase of the external field, which helps to illustrate the occurrence of such an interesting finding.

Ising model on evolution networks and its application on opinion formation

Many phenomena show that in a favorable circumstance an agent still has an updating possibility, and in an unfavor- able circumstance an agent also has a possibility of holding its own state and reselecting its neighbors. To describe this kind of phenomena an Ising model on evolution networks was presented and used for consensus formation and separation of opinion groups in human population. In this model the state-holding probability p and selection-rewiring probability q were introduced. The influence of this mixed dynamics of spin flips and network rewiring on the ordering behavior of the model was investigated, p hinders ordering of opinion networks and q accelerates the dynamical process of networks. Influence of q on the ordering and separating stems from its effect on average path length of networks.

Development and application of a two-dimensional water quality model for the Daqinghe River Mouth of the Dianchi Lake

Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.

Performance assessment of two-dimensional hydraulic models for generation of flood inundation maps in mountain river basins

Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.

A Solvable Decorated Ising Lattice Model

A decorated lattice is suggested and the Ising model on it with three kinds of interactions K1, K2, and K3 is studied. Using an equivalent transformation, the square decorated Ising lattice is transformed into a regular square Ising lattice with nearest-neighbor, next-nearest-nelghbor, and four-spin interactions, and the critical fixed point is found at K1 = 0.5769, K2= -0.0671, and K3 = 0.3428, which determines the critical temperature of the system. It is also found that this system and the regular square Ising lattice, and the eight-vertex model belong to the same universality class.

Mathematical structure of the three-dimensional(3D) Ising model

An overview of the mathematical structure of the three-dimensional（3D） Ising model is given from the points of view of topology,algebra,and geometry.By analyzing the relationships among transfer matrices of the 3D Ising model,Reidemeister moves in the knot theory,Yang-Baxter and tetrahedron equations,the following facts are illustrated for the 3D Ising model.1） The complex quaternion basis constructed for the 3D Ising model naturally represents the rotation in a（3＋1）-dimensional space-time as a relativistic quantum statistical mechanics model,which is consistent with the 4-fold integrand of the partition function obtained by taking the time average.2） A unitary transformation with a matrix that is a spin representation in 2 n·l·o-space corresponds to a rotation in 2n·l·o-space,which serves to smooth all the crossings in the transfer matrices and contributes the non-trivial topological part of the partition function of the 3D Ising model.3） A tetrahedron relationship would ensure the commutativity of the transfer matrices and the integrability of the 3D Ising model,and its existence is guaranteed by the Jordan algebra and the Jordan-von Neumann-Wigner procedures.4） The unitary transformation for smoothing the crossings in the transfer matrices changes the wave functions by complex phases φx,φy,and φz.The relationship with quantum field and gauge theories and the physical significance of the weight factors are discussed in detail.The conjectured exact solution is compared with numerical results,and the singularities at/near infinite temperature are inspected.The analyticity in β=1/（kBT） of both the hard-core and the Ising models has been proved only for β〉0,not for β=0.Thus the high-temperature series cannot serve as a standard for judging a putative exact solution of the 3D Ising model.

Ground-State Phase Diagram of Transverse Spin-2 Ising Model with Longitudinal Crystal-Field

The transverse spin-2 Ising ferromagnetic model with a longitudinal crystal-field is studied within the mean-field theory based on Bogoliubov inequality for the Gibbs free energy. The ground-state phase diagram and the tricritical point are obtained in the transverse field Ω/ zJ-longitudinal crystal D / zJ field plane. We find that there are the first order-order phase transitions in a very small range of D /zJ besides the usual first order-disorder phase transitions and the second order-disorder phase transitions,

Comment on 'Mathematical structure of the three-dimensional (3D) Ising model'

The review paper by Zhang Zhi-Dong （Zhang Z D 2013 Chin. Phys. B 22 030513, arXiv：1305.2956） contains many errors and is based on several earlier works that are equally wrong.

Design and Analysis of Integrated Predictive Iterative Learning Control for Batch Process Based on Two-dimensional System Theory

Based on the two-dimensional （2D） system theory, an integrated predictive iterative learning control （2D-IPILC） strategy for batch processes is presented. First, the output response and the error transition model predictions along the batch index can be calculated analytically due to the 2D Roesser model of the batch process. Then, an integrated framework of combining iterative learning control （ILC） and model predictive control （MPC） is formed reasonably. The output of feedforward ILC is estimated on the basis of the predefined process 2D model. By min- imizing a quadratic objective function, the feedback MPC is introduced to obtain better control performance for tracking problem of batch processes. Simulations on a typical batch reactor demonstrate that the satisfactory tracking performance as well as faster convergence speed can be achieved than traditional proportion type （P- t-we） ILC despite the model error and disturbances.

Critical behaviour of the ferromagnetic Ising model on a triangular lattice

We use a new updated algorithm scheme to investigate the critical behaviour of the two-dimensional ferromagnetic Ising model on a triangular lattice with the nearest neighbour interactions. The transition is examined by generating accurate data for lattices with L=8, 10, 12, 15, 20, 25, 30, 40 and 50. The updated spin algorithm we employ has the advantages of both a Metropolis algorithm and a single-update method. Our study indicates that the transition is continuous at Тc=3.6403（2）. A convincing finite-size scaling analysis of the model yields ν=0.9995（21）, β/ν=0.12400（17）, γ/v=1.75223（22）, γ^1/ν=1.7555（22）, α/ν=0.00077（420） （scaling） and α/ν=0.0010（42） （hyperscaling）. The present scheme yields more accurate estimates for all the critical exponents than the Monte Carlo method, and our estimates are shown to be in excellent agreement with their predicted values.

Study of Depolarization Field Influence on Ferroelectric Films Within Transverse Ising Model

An improved transverse Ising model is proposed by taking the depolarization field effect into account. Within the framework of mean-held theory we investigate the behavior of the ferroelectric thin film. Our results show that the influence of the depolarization field is to flatten the spontaneous polarization profile and make the films more homogeneous, which is consistent with Ginzburg Landau theory. This fact shows that this model can be taken as an effective model to deal with the ferroelectric film and can be further extended to refer to quantum effect. The competition between quantum effect and depolarization field induces some interesting phenomena on ferroelectric thin films.

Kinetic Ising model in a time-dependent oscillating external magnetic field:effective-field theory

Recently, Shiet al. [2008 Phys. Left. A 372 5922] have studied the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field and presented the dynamic phase diagrams by using an effective-field theory （EFT） and a mean-field theory （MFT）. The MFT results are in conflict with those of the earlier work of Tome and de Oliveira, [1990 Phys. Rev. A 41 4251]. We calculate the dynamic phase diagrams and find that our results are similar to those of the earlier work of Tome and de Oliveira; hence the dynamic phase diagrams calculated by Shiet al. are incomplete within both theories, except the low values of frequencies for the MFT calculation. We also investigate the influence of external field frequency （w） and static external field amplitude （h0） for both MFT and EFT calculations. We find that the behaviour of the system strongly depends on the values of w and h0.

HIGH RESOLUTION PARAMETRIC MODELLING FOR TWO-DIMENSIONAL RADAR TARGET USING PRONY ALGORITHM

On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high resolution estimation algorithm to extract the model parameters is also developed by building the relation of the scattering model and Prony model. The analysis of Cramer-Rao bound and simulation show that the method here has better statistical performance. The simulated analysis also indicates that the accurate extraction of the diffraction coefficient of scattering center is restricted by signal to noise ratio, radar center frequency and radar bandwidth.