Cell-dynamical simulation of magnetic hysteresis in the two-dimensional Ising system

We present results from numerical simulations using a ‘‘cell-dynamical system’’ to obtain solutions to the time-dependent Ginzburg-Landau equation for a scalar, two-dimensional (2D), (Φ2)2 model in the presence of a sinusoidal external magnetic field. Our results confirm a recent scaling law proposed by Rao, Krishnamurthy, and Pandit [Phys. Rev. B 42, 856 (1990)], and are also in excellent agreement with recent Monte Carlo simulations of hysteretic behavior of 2D Ising spins by Lo and Pelcovits [Phys. Rev. A 42, 7471 (1990)].

Conjectures about phase turbulence in the complex Ginzburg-Landau equation

In the complex Ginzburg-Landau equation, we consider possible ''phase turbulent'' regimes, where asymptotic correlations are controlled by phase fluctuations rather than by topological defects. Conjecturing that the decay of such correlations is governed by the Kardar-Parisi-Zhang (KPZ) model of growing interfaces, we derive the following results: (1) A scaling ansatz implies that equal-time spatial correlations in 1d, 2d, and 3d decay like e(-Ax2 zeta), where A is a nonuniversal constant, and zeta=1/2 in 1d. (2) Temporal correlations decay as exp(-t(2 beta)h(t/L(z))), with the scaling law <(beta)over bar> = <(zeta)over bar>/z, where z = 3/2, 1.58..., and 1.66..., for d = 1,2, and 3 respectively. The scaling function h(y) approaches a constant as y --> 0, and behaves like y(2(beta-<(beta)over bar>)), for large y. If in 3d the associated KPZ model turns out to be in its weak-coupling (''smooth'') phase, then, instead of the above behavior, the CGLE exhibits rotating long-range order whose connected correlations decay like 1/x in space or 1/t(1/2) in time. (3) For system sizes, L, and times t respectively less than a crossover length, L(c), and time, t(c), correlations are governed by the free-field or Edwards-Wilkinson (EW) equation, rather than the KPZ model. In 1d, we find that L(c) is large: L(c) similar to 35,000; for L < L(c) we show numerical evidence for stretched exponential decay of temporal correlations with an exponent consistent with the EW value beta(EW)= 1/4.

Critical properties of the double-exchange ferromagnet Nd(0.6)Pb(0.4)MnO(3)

Results of a study of dc magnetization M(T,H), performed on a Nd(0.6)Pb(0.4)MnO(3) single crystal in the temperature range around T(C) (Curie temperature) which embraces the supposed critical region \epsilon\=\T-T(C)\/T(C)less than or equal to0.05 are reported. The magnetic data analyzed in the critical region using the Kouvel-Fisher method give the values for the T(C)=156.47+/-0.06 K and the critical exponents beta=0.374+/-0.006 (from the temperature dependence of magnetization) and gamma=1.329+/-0.003 (from the temperature dependence of initial susceptibility). The critical isotherm M(T(C),H) gives delta=4.54+/-0.10. Thus the scaling law gamma+beta=deltabeta is fulfilled. The critical exponents obey the single scaling equation of state M(H,epsilon)=epsilon(beta)f(+/-)(H/epsilon(beta+gamma)), where f(+) for T>T(C) and f(-) for T

Linear stability, transient growth and the role of viscosity stratification in compressible Couette flow

Linear stability and the nonmodal transient energy growth in compressible plane Couette flow are investigated for two prototype mean flows: (a) the uniform shear flow with constant viscosity, and (b) the nonuniform shear flow with stratified viscosity. Both mean flows are linearly unstable for a range of supersonic Mach numbers (M). For a given M, the critical Reynolds number (Re) is significantly smaller for the uniform shear flow than its nonuniform shear counterpart; for a given Re, the dominant instability (over all streamwise wave numbers, α) of each mean flow belongs to different modes for a range of supersonic M. An analysis of perturbation energy reveals that the instability is primarily caused by an excess transfer of energy from mean flow to perturbations. It is shown that the energy transfer from mean flow occurs close to the moving top wall for “mode I” instability, whereas it occurs in the bulk of the flow domain for “mode II.” For the nonmodal transient growth analysis, it is shown that the maximum temporal amplification of perturbation energy, Gmax, and the corresponding time scale are significantly larger for the uniform shear case compared to those for its nonuniform counterpart. For α=0, the linear stability operator can be partitioned into L∼L̅ +Re2 Lp, and the Re-dependent operator Lp is shown to have a negligibly small contribution to perturbation energy which is responsible for the validity of the well-known quadratic-scaling law in uniform shear flow: G(t∕Re)∼Re2. In contrast, the dominance of Lp is responsible for the invalidity of this scaling law in nonuniform shear flow. An inviscid reduced model, based on Ellingsen-Palm-type solution, has been shown to capture all salient features of transient energy growth of full viscous problem. For both modal and nonmodal instability, it is shown that the viscosity stratification of the underlying mean flow would lead to a delayed transition in compressible Couette flow.

Investigations on micro-blast wave assisted metal foil forming for biomedical applications

The deformation dynamics of metal foils (<0.25 mm thick) subjected to micro-blast wave are presented in this paper. The energy of micro-blast wave emanating from the open end of a polymer tube is used to deliver micro-particles for bio-medical applications. In these experiments metal foils are used to transfer the energy of the micro-blast wave to the micro-particles. Using cubic root scaling law the over pressure of the blast wave at the open end of the polymer tube is estimated and using this peak plate over pressure is estimated. The finite element analysis is used to estimate the velocity profile of the deforming metal foils. The finite element analysis results are compared with experimental results for the maximum deformation and deformed shape. Based on the deformation velocity, metal foil to be used for experiments is selected. Among the materials investigated 0.1 mm thick brass foil has the maximum velocity of 205 m/s and is used in the experiments. It is found from finite element analysis that the particles deposited within a radius of 0.5 mm will leave the foil with nearly equal velocity (error < 5%). The spray cone angle which is the angle of deviation of the path of particles from the axis of the polymer tube is also estimated and found to be less than 7 degrees up to a radius of 0.75 mm. Illustrative experiments are carried out to deliver micro particles (0.7 mu m diameter tungsten) into plant tissues. Particle penetration depth up to 460 mu m was achieved in ground tissue of potato tuber. (C) 2012 Elsevier Ltd. All rights reserved.

Evolution of ferromagnetism from a frustrated state in LixNi(2-x)O2 (0.67 < x < 0.98)

Two-dimensional triangular-lattice antiferromagnetic systems continue to be an interesting area in condensed matter physics and LiNiO2 is one such among them. Here we present a detailed experimental magnetic study of the quasi-stoichiometric LixNi2-xO2 system (0.67scaling law: delta = 1 + gamma beta(-1). Furthermore, the universality class of the scaling relations is verified, where the scaled m and scaled h collapse into two branches. Finally, based on our observations, a phase diagram is constructed.

Evidence for 3D isotropic long range spin-spin interaction near the ferromagnetic transition in bulk and thin film SrRuO3

In the case of metallic ferromagnets there has always been a controversy, i.e. whether the magnetic interaction is itinerant or localized. For example SrRuO3 is known to be an itinerant ferromagnet where the spin-spin interaction is expected to be mean field in nature. However, it is reported to behave like Ising, Heisenberg or mean field by different groups. Despite several theoretical and experimental studies and the importance of strongly correlated systems, the experimental conclusion regarding the type of spin-spin interaction in SrRuO3 is lacking. To resolve this issue, we have investigated the critical behaviour in the vicinity of the paramagnetic-ferromagnetic phase transition using various techniques on polycrystalline as well as (001) oriented SrRuO3 films. Our analysis reveals that the application of a scaling law in the field-cooled magnetization data extracts the value of the critical exponent only when it is measured at H -> 0. To substantiate the actual nature without any ambiguity, the critical behavior is studied across the phase transition using the modified Arrott plot, Kouvel-Fisher plot and M-H isotherms. The critical analysis yields self-consistent beta, gamma and delta values and the spin interaction follows the long-range mean field model. Further the directional dependence of the critical exponent is studied in thin films and it reveals the isotropic nature. It is elucidated that the different experimental protocols followed by different groups are the reason for the ambiguity in determining the critical exponents in SrRuO3.

Evidence for 3D isotropic long range spin-spin interaction near the ferromagnetic transition in bulk and thin film SrRuO3

In the case of metallic ferromagnets there has always been a controversy, i.e. whether the magnetic interaction is itinerant or localized. For example SrRuO3 is known to be an itinerant ferromagnet where the spin-spin interaction is expected to be mean field in nature. However, it is reported to behave like Ising, Heisenberg or mean field by different groups. Despite several theoretical and experimental studies and the importance of strongly correlated systems, the experimental conclusion regarding the type of spin-spin interaction in SrRuO3 is lacking. To resolve this issue, we have investigated the critical behaviour in the vicinity of the paramagnetic-ferromagnetic phase transition using various techniques on polycrystalline as well as (001) oriented SrRuO3 films. Our analysis reveals that the application of a scaling law in the field-cooled magnetization data extracts the value of the critical exponent only when it is measured at H -> 0. To substantiate the actual nature without any ambiguity, the critical behavior is studied across the phase transition using the modified Arrott plot, Kouvel-Fisher plot and M-H isotherms. The critical analysis yields self-consistent beta, gamma and delta values and the spin interaction follows the long-range mean field model. Further the directional dependence of the critical exponent is studied in thin films and it reveals the isotropic nature. It is elucidated that the different experimental protocols followed by different groups are the reason for the ambiguity in determining the critical exponents in SrRuO3.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Threshold diversity and trans-scales sensitivity in a finite nonlinear evolution model of materials failure

We present a slice-sampling method and study the ensemble evolution of a large finite nonlinear system in order to model materials failure. There is a transitional region of failure probability. Its size effect is expressed by a slowly decaying scaling law. In a meso-macroscopic range (similar to 10(5)) in realistic failure, the diversity cannot be ignored. Sensitivity to mesoscopic details governs the phenomena. (C) 1997 Published by Elsevier Science B.V.

Effects of physical parameter range on dimensionless variable sensitivity in water flooding reservoirs

The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.

Effects of physical parameter range on dimensionless variable sensitivity in water flooding reservoirs

The similarity criterion for water flooding reservoir flows is concerned with in the present paper. When finding out all the dimensionless variables governing this kind of flow, their physical meanings are subsequently elucidated. Then, a numerical approach of sensitivity analysis is adopted to quantify their corresponding dominance degree among the similarity parameters. In this way, we may finally identify major scaling law in different parameter range and demonstrate the respective effects of viscosity, permeability and injection rate.

Ablation of hydrogen isotopic spherical pellet due to the impact of energetic ions of the respective isotopes

The ablation rate of a hydrogen isotopic spherical pellet G(is) due to the impact of energetic ions of the respective isotopes and its scaling law are obtained using the transsonic neutral-shielding model, where subscript s might refer to either hydrogen or deuterium. Numerical results show that if E0s/E0e2 greater-than-or-equal-to 1.5, G(is)/G(es) greater-than-or-equal-to 20%, where E0s and E0e are the energy of undisturbed ion and electron, respectively, and G(es) is the ablation rate of a pellet due to the impact of electrons. Hence, under the NBI heating, the effect of the impact of energetic ions on the pellet ablation should be taken into consideration. This result also gives an explanation of the observed enhancement of pellet ablation during NBIH.

二维槽道湍流的标度律分析

对不同Reynolds数的二维不可压槽道湍流（即二维槽道流动的扰动饱和态）进行了标度律分析。指出二维槽道湍流中存在明显的标度律及扩展自相似性；分析了标度数随着流动Reynolds数的变化，指出随着Reynolds数的增长，标度指数的奇异性并没有减小的趋势；由此推测充分发展的二维槽道湍流的标度指数也是奇异的。将所得计算值和佘振苏等人的SL标度律预测值相比较，认为佘等人的SL标度律公式能很好描述二维湍流。

槽道湍流的直接数值模拟

通过直接数值模拟(DNS)研究槽道湍流的性质和机理。包含五个部分：1）湍流直接数值模拟的差分方法研究。2）求解不可压N－S方程的高效算法和不可压槽道湍流的直接数值模拟。3）可压缩槽道湍流的直接数值模拟和压缩性机理分析。4）“二维湍流”的机理分析。5）槽道湍流的标度律分析。1．针对壁湍流计算网格变化剧烈的特点，构造了基于非等距网格的的迎风紧致格式。该方法直接针对计算网格构造格式中的系数，克服了传统方法采用 Jacobian 变换因网格变化剧烈而带来的误差。针对湍流场的多尺度特性分析了差分格式的精度、网格尺度与数值模拟能分辨的最小尺度的关系，给出不同差分格式对计算网格步长的限制。同时分析了计算中混淆误差的来源和控制方法，指出了迎风型紧致格式能很好地控制混淆误差。2．将上述格式与三阶精度的Adams半隐格式相结合，构造了不可压槽道湍流直接数值模拟的高效算法。该算法利用基于交错网格的离散形式的压力Poisson方程求解压力项，避免了压力边界条件处理的困难。利用FFT对方程中的隐式部分进行解耦，解耦后的方程采用追赶法（LU分解法）求解，大大减少了计算量。为了检验该方法，进行了三维不可压槽道湍流的直接数值模拟，得到了Re＝2800的充分发展不可压槽道湍流，并对该湍流场进行了统计分析。包括脉动速度偏斜因子在内的各阶统计量与实验结果及Kim等人的计算结果吻合十分理想，说明本方法是行之有效的。3．进行了三维充分发展的可压缩槽道湍流的直接数值模拟。得到了 Re=3300，Ma=0.8的充分发展可压槽道湍流的数据库。流场的统计特征（如等效平均速度分布，“半局部”尺度无量纲化的脉动速度均方根）和他人的数值计算结果吻合。得到了可压槽道湍流的各阶统计量，其中脉动速度的偏斜因子和平坦因子等高阶统计量尚未见其他文献报道。同时还分析了压缩性效应对壁湍流影响的机理，指出近壁处的压力－膨胀项将部分湍流脉动的动能转换成内能，使得可压湍流近壁速度条带结构更加平整。4．模拟了二维不可压槽道流动的饱和态（所谓“二维湍流”），分析了“二维槽道湍流”的非线性行为特征。分析了流场中的上抛－下扫和间歇现象，研究了“二维湍流”与三维湍流的区别。指出“二维湍流”反映了三维湍流的部分特征，同时指出了展向扰动对于湍流核心区发展的重要性。5．首次对可压缩槽道湍流及“二维槽道湍流”标度律进行了分析，得出了以下结论：a）槽道湍流中，在槽道中心线附近较宽的区域，存在标度律。b）该区域流场存在扩展自相似性（ESS）。c）在Mach数不是很高时，压缩性对标度指数影响不大。本文结果同SL标度律的理论值吻合较好，有效支持了该理论。对“二维槽道湍流”也有相似的结论，但与三维湍流不同的是，“二维槽道湍流”存在标度律的区域更宽，近壁处的标度指数比中心处有所升高。