Study on the critical stress threshold of weakly cemented sandstone damage based on the renormalization group method

During the microstructural analysis of weakly cemented sandstone,the granule components and ductile structural parts of the sandstone are typically generalized.Considering the contact between granules in the microstructure of weakly cemented sandstone,three basic units can be determined:regular tetrahedra,regular hexahedra,and regular octahedra.Renormalization group models with granule-and pore-centered weakly cemented sandstone were established,and,according to the renormalization group transformation rule,the critical stress threshold of damage was calculated.The results show that the renormalization model using regular octahedra as the basic units has the highest critical stress threshold.The threshold obtained by iterative calculations of the granule-centered model is smaller than that obtained by the pore-centered model.The granule-centered calculation provides the lower limit(18.12%),and the pore-centered model provides the upper limit(36.36%).Within this range,the weakly cemented sandstone is in a phase-like critical state.That is,the state of granule aggregation transforms from continuous to discrete.In the relative stress range of 18.12%-36.36%,the weakly cemented sandstone exhibits an increased proportion of high-frequency signals(by 83.3%)and a decreased proportion of low-frequency signals(by 23.6%).The renormalization calculation results for weakly cemented sandstone explain the high-low frequency conversion of acoustic emission signals during loading.The research reported in this paper has important significance for elucidating the damage mechanism of weakly cemented sandstone.

Renormalization Group Theory and Its Application to Thermally-Induced Turbulence

Renormalization group theory applied to turbulence will be reviewed in this article.Techniques associated are used for analyzing thermally-induced turbulence.Transport properties such as effective viscosity and thermal diffusivity are derived.

Derivation of a second-order model for Reynolds stress using renormalization group analysis and the two-scale expansion technique

With the two-scale expansion technique proposed by Yoshizawa,the turbulent fluctuating field is expanded around the isotropic field.At a low-order two-scale expansion,applying the mode coupling approximation in the Yakhot-Orszag renormalization group method to analyze the fluctuating field,the Reynolds-average terms in the Reynolds stress transport equation,such as the convective term,the pressure-gradient-velocity correlation term and the dissipation term,are modeled.Two numerical examples：turbulent flow past a backward-facing step and the fully developed flow in a rotating channel,are presented for testing the efficiency of the proposed second-order model.For these two numerical examples,the proposed model performs as well as the Gibson-Launder （GL） model,giving better prediction than the standard k-ε model,especially in the abilities to calculate the secondary flow in the backward-facing step flow and to capture the asymmetric turbulent structure caused by frame rotation.

Specific Heat and Magnetic Susceptibility of Graphene:A Renormalization Group Study

In the present paper, we study effect of the long-range Coulomb interaction on the thermodynamic propertiesof graphene by renormalization group methods.Our calculations show that both the specific heat and the magneticsusceptibility of the material behave differently from the Landau Fermi liquid.More precisely, we find that thesequantities are logarithmically suppressed with respect to its noninteracting counterpart when temperature is low.

Time-Dependent Density Matrix Renormalization Group Coupled with n-Mode Representation Potentials for the Excited State Radiationless Decay Rate:Formalism and Application to Azulene

We propose a method for calculating the nonradiative decay rates for polyatomic molecules including anharmonic effects of the potential energy surface(PES)in the Franck-Condon region.The method combines the n-mode repre-sentation method to construct the ab initio PES and the nearly exact time-dependent density matrix renormalization group method(TD-DMRG)to simulate quantum dynamics.In addition,in the framework of TD-DMRG,we further develop an algorithm to calculate the final-state-resolved rate coefficient which is very useful to analyze the contribution from each vibrational mode to the transition process.We use this method to study the internal conversion(IC)process of azulene after taking into account the anharmonicity of the ground state PES.The results show that even for this semi-rigid molecule,the intramode anharmonicity enhances the IC rate significantly,and after considering the two-mode coupling effect,the rate increases even further.The reason is that the anharmonicity enables the C-H vibrations to receive electronic energy while C-H vibrations do not contribute on the harmonic PES as the Huang-Rhys factor is close to 0.

Real-space parallel density matrix renormalization group with adaptive boundaries

We propose an improved real-space parallel strategy for the density matrix renormalization group(DMRG)method,where boundaries of separate regions are adaptively distributed during DMRG sweeps.Our scheme greatly improves the parallel efficiency with shorter waiting time between two adjacent tasks,compared with the original real-space parallel DMRG with fixed boundaries.We implement our new strategy based on the message passing interface(MPI),and dynamically control the number of kept states according to the truncation error in each DMRG step.We study the performance of the new parallel strategy by calculating the ground state of a spin-cluster chain and a quantum chemical Hamiltonian of the water molecule.The maximum parallel efficiencies for these two models are 91%and 76%in 4 nodes,which are much higher than the real-space parallel DMRG with fixed boundaries.

Renormalization group methods for a Mathieu equation with delayed feedback

This paper presents the application of the renormalization group （RG） methods to the delayed differential equation. By analyzing the Mathieu equation with time delay feedback, we get the amplitude and phase equations, and then obtain the approximate solutions by solving the corresponding RG equations. It shows that the approximate solutions obtained from the RG method are superior to those from the conventionally perturbation methods.

Improved hybrid parallel strategy for density matrix renormalization group method

We propose a new heterogeneous parallel strategy for the density matrix renormalization group(DMRG)method in the hybrid architecture with both central processing unit(CPU)and graphics processing unit(GPU).Focusing on the two most time-consuming sections in the finite DMRG sweeps,i.e.,the diagonalization of superblock and the truncation of subblock,we optimize our previous hybrid algorithm to achieve better performance.For the former,we adopt OpenMP application programming interface on CPU and use our own subroutines with higher bandwidth on GPU.For the later,we use GPU to accelerate matrix and vector operations involving the reduced density matrix.Applying the parallel scheme to the Hubbard model with next-nearest hopping on the 4-leg ladder,we compute the ground state of the system and obtain the charge stripe pattern which is usually observed in high temperature superconductors.Based on simulations with different numbers of DMRG kept states,we show significant performance improvement and computational time reduction with the optimized parallel algorithm.Our hybrid parallel strategy with superiority in solving the ground state of quasi-two dimensional lattices is also expected to be useful for other DMRG applications with large numbers of kept states,e.g.,the time dependent DMRG algorithms.

Phase Diagrams of the Semi-Infinite Blume-Capel Model with Mixed Spins (SA = 1 and SB = 3/2) by Migdal Kadanoff Renormalization Group

We study the mixed spin-1 and spin-3/2 Blume-Capel model under crystal field in the tridimensional semi-infinite case. This has been done by using the real-space renormalization group approximation and specifically the Migdal-Kadanoff technique. As a function of the ratio R of bulk and surface interactions and the ratios R1 and R2 of bulk and surface crystals fields on the spin-1 and spin-3/2 respectively, we have determined various types of phase diagrams. Besides second- order transition lines, first-order phase transition lines terminating at tricritical points are obtained. We found that there existed nine main types of phase diagram showing a variety of phase transitions associated with the surface, including ordinary, extraordinary, surface and special phase transitions.

Renormalization group theory for temperature-driven first-order phase transitions in scalar models

We study the scaling and universal behavior of temperature-driven first-order phase transitions in scalar models. These transitions are found to exhibit rich phenomena, though they are controlled by a single complex-conjugate pair of imaginary fixed points of φ3 theory. Scaling theories and renormalization group theories are developed to account for the phenomena, and three universality classes with their own hysteresis exponents are found： a field-like thermal class, a partly thermal class, and a purely thermal class, designated, respectively, as Thermal Classes I, II, and III. The first two classes arise from the opposite limits of the scaling forms proposed and may cross over to each other depending on the temperature sweep rate. They are both described by a massless model and a purely massive model, both of which are equivalent and are derived from φ3 theory via symmetry. Thermal Class III characterizes the cooling transitions in the absence of applied external fields and is described by purely thermal models, which include cases in which the order parameters possess different symmetries and thus exhibit different universality classes. For the purely thermal models whose free energies contain odd-symmetry terms, Thermal Class III emerges only at the mean-field level and is identical to Thermal Class II. Fluctuations change the model into the other two models. Using the extant three- and two- loop results for the static and dynamic exponents for the Yang-Lee edge singularity, respectively, which falls into the same universality class as φ3 theory, we estimate the thermal hysteresis exponents of the various classes to the same precision. Comparisons with numerical results and experiments are briefly discussed.

Natural orbitals renormalization group approach to a Kondo singlet

A magnetic impurity embedded in a Fermi sea is collectively screened by a cloud of conduction electrons to form a Kondo singlet below a characteristic energy scale TK,the Kondo temperature,through the mechanism of the Kondo effect.We have reinvestigated the Kondo singlet by means of the newly developed natural orbitals renormalization group(NORG)method.We find that,in the framework of natural orbitals formalism,the Kondo screening mechanism becomes transparent and simple,while the intrinsic structure of a Kondo singlet is clearly resolved.For a single impurity Kondo system in whichever case of either finite size or thermodynamic limit,there exists a single active natural orbital that screens the magnetic impurity dominantly.In the perspective of entanglement,the magnetic impurity is entangled dominantly with the active natural orbital,i.e.,the subsystem formed by the active natural orbital and the magnetic impurity basically disentangles from the remaining system.We have also studied the structures of the active natural orbital respectively projected into real space and momentum space.Moreover,the dynamical properties,represented by one-particle Green’s functions defined at the active natural orbital,are obtained by the correction vector method.Meanwhile,the well-known Kondo resonance is clearly observed in the spectral function at the active natural orbital.To realize the thermodynamic limit,the Wilson chains with the numerical renormalization group approach are employed.

No-core Monte Carlo shell model calculations with unitary correlation operator method and similarity renormalization group

The unitary correlation operator method （UCOM） and the similarity renormalization group theory （SRG） are compared and discussed in the framework of the no-core Monte Carlo shell model （MCSM） calculations for ^3H and ^4He. The treatment of spurious center-of-mass motion by Lawson＇s prescription is performed in the MCSM calculations. These results with both transformed interactions show good suppression of spurious center-of-mass motion with proper Lawson＇s prescription parameter βc.m. values. The UCOM potentials obtain faster convergence of total energy for the ground state than that of SRG potentials in the MCSM calculations, which differs from the cases in the no-core shell model calculations （NCSM）. These differences are discussed and analyzed in terms of the truncation scheme in the MCSM and NCSM, as well as the properties of the potentials of SRG and UCOM.

Singularly Perturbed Renormalization Group Method and Its Significance in Dynamical Systems Theory

In this paper,we mainly investigate three topics on the renormalization group(RG)method to singularly perturbed problems:1)We will present an explicit strategy of RG procedure to get the approximate solution up to any order.2)We will refer that the RG procedure can,in fact,be used to get the normal form of differential dynamical systems.3)We will also present the approximating center manifolds of the perturbed systems,and investigate the long time asymptotic behavior by means of RG formula.

Wilsonian Renormalization Group and the Lippmann-Schwinger Equation with a Multitude of Cutoff Parameters

The Wilsonian renormalization group approach to the Lippmann-Schwinger equation with a multitude of cutoff parameters is introduced.A system of integro-differential equations for the cutoff-dependent potential is obtained.As an illustration,a perturbative solution of these equations with two cutoff parameters for a simple case of an S-wave low-energy potential in the form of a Taylor series in momenta is obtained.The relevance of the obtained results for the effective field theory approach to nucleon-nucleon scattering is discussed.

Study on critical conditions for rock failure by means of group renormalization

A study of the characteristics of the accumulative rock failure and its evolution byapplication of the group renormalization method were presented. In addition, the interactionand long-range correlated effects between the immediate neighboring units was studied.The concept of mechanical transference for model OFC, employed in the study ofself-organized criticality, and the coefficient a were introduced into the calculation model forgroup renormalization. With the introduction, mechanisms for the drastic increase and decrease of failure intensity of rocks were investigated under similar macro-conditions.

Stability of the topological quantum critical point between multi-Weyl semimetal and band insulator

One could tune a topological double-Weyl semimetal or a topological triple-Weyl semimetal to become a topologically trivial insulator by opening a band gap.This kind of quantum phase transition is characterized by the change of certain topological invariant.A new gapless semimetallic state emerges at each topological quantum critical point.Here we perform a renormalization group analysis to investigate the stability of such critical points against perturbations induced by random scalar potential and random vector potential.We find that the quantum critical point between double-Weyl semimetal and band insulator is unstable and can be easily turned into a compressible diffusive metal by any type of weak disorder.The quantum critical point between triple-Weyl semimetal and band insulator flows to a stable strong-coupling fixed point if the system contains a random vector potential merely along the z-axis,but becomes a compressible diffusive metal when other types of disorders exist.

Effects of quantum quench on entanglement dynamics in antiferromagnetic Ising model

We study the relationship between quench dynamics of entanglement and quantum phase transition in the antiferromagnetic Ising model with the Dzyaloshinskii–Moriya(DM)interaction by using the quantum renormalization-group method and the definition of negativity.Two types of quench protocols(i)adding the DM interaction suddenly and(ii)rotating the spins around x axis are considered to drive the dynamics of the system,respectively.By comparing the behaviors of entanglement in both types of quench protocols,the effects of quench on dynamics of entanglement are studied.It is found that there is the same characteristic time at which the negativity firstly reaches its maximum although the system shows different dynamical behaviors.Especially,the characteristic time can accurately reflect the quantum phase transition from antiferromagnetic to saturated chiral phases in the system.In addition,the correlation length exponent can be obtained by exploring the nonanalytic and scaling behaviors of the derivative of the characteristic time.

Renormalization-group theory of first-order phase transition dynamics in field-driven scalar model

Through a detailed study of the mean-field approximation, the Gaussian approximation, the perturbation expansion, and the field-theoretic renormalization-group analysis of a φ^3 theory, we show that the instability fixed points of the theory, together with their associated instability exponents, are quite probably relevant to the scaling and universality behavior exhibited by the first-order phase transitions in a field-driven scalar Ca model, below its critical temperature and near the instability points. Finite- time scaling and leading corrections to the scaling are considered. We also show that the instability exponents of the first-order phase transitions are equivalent to those of the Yang-Lee edge singularity, and employ the latter to improve our estimates of the former. The outcomes agree well with existing numerical results.

Mechanism of water-inrush from fault induced by mining near the working face

Adopted the fractal tree-like failure model, and established the renormalization group transform function of fractured fault, and investigated the mechanism of water-inrush from fault, and found out the critical probability of water-inrush from fault caused by fault fracture. The results indicate： when the failure rate P is less than the critical failure rate Pc=0.206 3, the failure of the system is just partial. When P is more than the critical failure rate Pc=0.206 3, the random distributed crannies concentrate to certain domain of attraction （such as the maximum shear stress face in the fault） gradually. The process will continue until the crannies run-through, forming conductivity channel, and cause water-inrush from fault.

Numerical Simulation and Experimental Study of the Effects of Disposal Space on the Flow Field Around the Combined Three-Tube Reefs

The artificial reefs placed on the seabed with different layouts and disposal spaces will produce variational flow field. The intensity and scale of the combined three-tube artificial reefs with different layouts at five Reynolds numbers（Re） are numerically investigated by use of the RNG k-ε turbulent model and SIMPLEC algorithm. A stationary no-slip boundary condition is used on the models and the bottoms, and the free surface is treated as a ＂moving wall＂ with zero shear force and the same velocity with inflow. In order to validate the simulation results, a particle image velocimetry（PIV） experiment is carried out to analyze the flow field. The numerical simulation results are consistent with the data obtained from experiment. The corresponding errors are all below 20%. Based on the validation, the effects of disposal space on flow field are simulated and analyzed. According to the simulation, in a parallel combination, a better artificial reef effect is obtained when the disposal space between two parallel reefs is 1.0L（L is the length of the combined three-tube reef model）. In a vertical combination, when the disposal space between two vertical reefs is 1.0L to 2.0L, the artificial reef effect is better.