Noise-induced phase transition in the Vicsek model through eigen microstate methodology

This paper presents a comprehensive framework for analyzing phase transitions in collective models such as theVicsek model under various noise types. The Vicsek model, focusing on understanding the collective behaviors of socialanimals, is known due to its discontinuous phase transitions under vector noise. However, its behavior under scalar noiseremains less conclusive. Renowned for its efficacy in the analysis of complex systems under both equilibrium and nonequilibriumstates, the eigen microstate method is employed here for a quantitative examination of the phase transitions inthe Vicsek model under both vector and scalar noises. The study finds that the Vicsek model exhibits discontinuous phasetransitions regardless of noise type. Furthermore, the dichotomy method is utilized to identify the critical points for thesephase transitions. A significant finding is the observed increase in the critical point for discontinuous phase transitions withescalation of population density.

A feedback control method for phase signal demodulation in fber-optic hydrophones

In the realm of acoustic signal detection,the identification of weak signals,particularly in the presence of negative signal-to-noise ratios,poses a significant challenge.This challenge is further heightened when signals are acquired through fiber-optic hydrophones,as these signals often lack physical significance and resist clear systematic modeling.Conventional processing methods,e.g.,low-pass filter(LPF),require a thorough understanding of the effective signal bandwidth for noise reduction,and may introduce undesirable time lags.This paper introduces an innovative feedback control method with dual Kalman filters for the demodulation of phase signals with noises in fiber-optic hydrophones.A mathematical model of the closed-loop system is established to guide the design of the feedback control,aiming to achieve a balance with the input phase signal.The dual Kalman filters are instrumental in mitigating the effects of signal noise,observation noise,and control execution noise,thereby enabling precise estimation for the input phase signals.The effectiveness of this feedback control method is demonstrated through examples,showcasing the restoration of low-noise signals,negative signal-to-noise ratio signals,and multi-frequency signals.This research contributes to the technical advancement of high-performance devices,including fiber-optic hydrophones and phase-locked amplifiers.

Emergent topological ordered phase for the Ising-XY model revealed by cluster-updating Monte Carlo method

The two-component cold atom systems with anisotropic hopping amplitudes can be phenomenologically described by a two-dimensional Ising-XY coupled model with spatial anisotropy.At low temperatures,theoretical predictions[Phys.Rev.A 72053604(2005)]and[arXiv:0706.1609]indicate the existence of a topological ordered phase characterized by Ising and XY disorder but with 2XY ordering.However,due to ergodic difficulties faced by Monte Carlo methods at low temperatures,this topological phase has not been numerically explored.We propose a linear cluster updating Monte Carlo method,which flips spins without rejection in the anisotropy limit but does not change the energy.Using this scheme and conventional Monte Carlo methods,we succeed in revealing the nature of topological phases with half-vortices and domain walls.In the constructed global phase diagram,Ising and XY-type transitions are very close to each other and differ significantly from the schematic phase diagram reported earlier.We also propose and explore a wide range of quantities,including magnetism,superfluidity,specific heat,susceptibility,and even percolation susceptibility,and obtain consistent and reliable results.Furthermore,we observed first-order transitions characterized by common intersection points in magnetizations for different system sizes,as opposed to the conventional phase transition where Binder cumulants of various sizes share common intersections.The critical exponents of different types of phase transitions are reasonably fitted.The results are useful to help cold atom experiments explore the half-vortex topological phase.

Data-driven optimization study of the multi-relaxation-time lattice Boltzmann method for solid-liquid phase change

Sharp phase interfaces and accurate temperature distributions are important criteria in the simulation of solid-liquid phase changes.The multi-relaxation-time lattice Boltzmann method(MRT-LBM)shows great numerical performance during simulation;however,the value method of the relaxation parameters needs to be specified.Therefore,in this study,a random forest(RF)model is used to discriminate the importance of different relaxation parameters to the convergence,and a support vector machine(SVM)is used to explore the decision boundary of the convergent samples in each dimensional model.The results show that the convergence of the samples is consistent with the sign of the decision number,and two types of the numerical deviations appear,i.e.,the phase mushy zone and the non-physical heat transfer.The relaxation parameters chosen on the decision boundary can further suppress the numerical bias and improve numerical accuracy.

Volumetric lattice Boltzmann method for pore-scale mass diffusionadvection process in geopolymer porous structures

Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.

Phase composition,conductivity,and sensor properties of cerium-doped indium oxide

The hydrothermal synthesis of In_(2)O_(3)and CeO_(2)–In_(2)O_(3)is investigated as well as the properties of sensor layers based on these compounds.During the synthesis of In_(2)O_(3),intermediate products In(OH)_(3)and InOOH are formed,which are the precursors of stable cubic(c-In_(2)O_(3))and metastable rhombohedral(rh-In_(2)O_(3))phases,respectively.A transition from c-In_(2)O_(3)to rh-In_(2)O_(3)is observed with the addition of CeO_(2).The introduction of cerium into rh-In_(2)O_(3)results in a decrease in the sensor response to hydrogen,while it increases in composites based on c-In_(2)O_(3).The data on the sensor activity of the composites correlate with XPS results in which CeO_(2)causes a decrease in the concentrations of chemisorbed oxygen and oxygen vacancies in rh-In_(2)O_(3).The reverse situation is observed in composites based on c-In_(2)O_(3).Compared to In_(2)O_(3)and CeO_(2)–In_(2)O_(3)obtained by other methods,the synthesized composites demonstrate maximum response to H_(2)at low temperatures by 70–100℃,and have short response time(0.2–0.5 s),short recovery time(6–7 s),and long-term stability.A model is proposed for the dependence of sensitivity on the direction of electron transfer between In_(2)O_(3)and CeO_(2).

Phase equilibrium data prediction and process optimizationin butadiene extraction process

In response to the lack of reliable physical parameters in the process simulation of the butadiene extraction,a large amount of phase equilibrium data were collected in the context of the actual process of butadiene production by acetonitrile.The accuracy of five prediction methods,UNIFAC(UNIQUAC Functional-group Activity Coefficients),UNIFAC-LL,UNIFAC-LBY,UNIFAC-DMD and COSMO-RS,applied to the butadiene extraction process was verified using partial phase equilibrium data.The results showed that the UNIFAC-DMD method had the highest accuracy in predicting phase equilibrium data for the missing system.COSMO-RS-predicted multiple systems showed good accuracy,and a large number of missing phase equilibrium data were estimated using the UNIFAC-DMD method and COSMO-RS method.The predicted phase equilibrium data were checked for consistency.The NRTL-RK(non-Random Two Liquid-Redlich-Kwong Equation of State)and UNIQUAC thermodynamic models were used to correlate the phase equilibrium data.Industrial device simulations were used to verify the accuracy of the thermodynamic model applied to the butadiene extraction process.The simulation results showed that the average deviations of the simulated results using the correlated thermodynamic model from the actual values were less than 2%compared to that using the commercial simulation software,Aspen Plus and its database.The average deviation was much smaller than that of the simulations using the Aspen Plus database(>10%),indicating that the obtained phase equilibrium data are highly accurate and reliable.The best phase equilibrium data and thermodynamic model parameters for butadiene extraction are provided.This improves the accuracy and reliability of the design,optimization and control of the process,and provides a basis and guarantee for developing a more environmentally friendly and economical butadiene extraction process.

Dual-phase coexistence enables to alleviate resistance drift in phase-change films

The amorphous phase-change materials with spontaneous structural relaxation leads to the resistance drift with the time for phase-change neuron synaptic devices. Here, we modify the phase change properties of the conventional Ge_2Sb_2Te_5(GST) material by introducing an SnS phase. It is found that the resistance drift coefficient of SnS-doped GST was decreased from 0.06 to 0.01. It can be proposed that the origin originates from the precipitation of GST nanocrystals accompanied by the precipitation of SnS crystals compared to single-phase GST compound systems. We also found that the decrease in resistance drift can be attributed to the narrowed bandgap from 0.65 to 0.43 eV after SnS-doping. Thus, this study reveals the quantitative relationship between the resistance drift and the band gap and proposes a new idea for alleviating the resistance drift by composition optimization, which is of great significance for finding a promising phase change material.

Stress-corrosion coupled damage localization induced by secondary phases in bio-degradable Mg alloys:phase-field modeling

In this study,a phase-field scheme that rigorously obeys conservation laws and irreversible thermodynamics is developed for modeling stress-corrosion coupled damage(SCCD).The coupling constitutive relationships of the deformation,phase-field damage,mass transfer,and electrostatic field are derived from the entropy inequality.The SCCD localization induced by secondary phases in Mg is numerically simulated using the implicit iterative algorithm of the self-defined finite elements.The quantitative evaluation of the SCCD of a C-ring is in good agreement with the experimental results.To capture the damage localization,a micro-galvanic corrosion domain is defined,and the buffering effect on charge migration is explored.Three cases are investigated to reveal the effect of localization on corrosion acceleration and provide guidance for the design for resistance to SCCD at the crystal scale.

Determining beam transverse absolute position by triangulation of multi-electrode signal phase differences

Accurate measurement of the transverse position of a beam is crucial in particle accelerators because it plays a key role in determining the beam parameters.Existing methods for beam-position measurement rely on the detection of image currents induced on electrodes or narrow-band wake field induced by a beam passing through a cavity-type structure.However,these methods have limitations.The indirect measurement of multiple parameters is computationally complex,requiring external calibration to determine the system parameters in advance.Furthermore,the utilization of the beam signal information is incomplete.Hence,this study proposes a novel method for measuring the absolute electron beam transverse position.By utilizing the geometric relationship between the center position of the measured electron beam and multiple detection electrodes and by analyzing the differences in the arrival times of the beam signals detected by these electrodes,the absolute transverse position of the electron beam crossing the electrode plane can be calculated.This method features absolute position measurement,a position sensitivity coefficient independent of vacuum chamber apertures,and no requirement for a symmetrical detector electrode layout.The feasibility of this method is validated through numerical simulations and beam experiments.

Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

Ageing process of pre-precipitation phase in Ni_(0.75)Al_(0.05)Fe_(0.2) alloy based on phase field method

By utilizing phase field method combined with analysis on free energy and interatomic potentials, pre-precipitation phase formation and transformation process of Ni0.75Al0.05Fe0.2 alloy in early precipitation stage during the ageing process under 1 000 K were studied. And free energy, microstructures, compositions and volume fractions of pre-precipitation phase and equilibrium phase were analyzed. The simulation results indicate that nonstoichiometric Llo pre-precipitation phase formed first, and then would gradually transform into L12 equilibrium phase. It is discovered that the phase transformation process was closely related to free energy and interatomic potentials. Additionally, it is revealed that free energy of Llo pre-precipitation phase was higher and interatomic potential was smaller than that of L12 equilibrium phase. Therefore, it is concluded that Llo phase was unstable, and phase transformation would occur to L12 which was more stable.

Effect of forced lamina flow on microsegregation simulated by phase field method quantitatively

The influence of supercooled melt forced lamina flow on microsegregation was investigated. The concentration distribution at solid-liquid boundary of binary alloy Ni-Cu was simulated using phase field model coupled with flow field. The microsegregation, concentration maximum value, boundary thickness of concentration near upstream dendrite and normal to flow dendrite, and downstream dendrite were studied quantitatively in the case of forced lamia flow. The simulation results show that solute field and flow field interact complexly. Compared with melt without flow, in front of upstream dendrite tip, the concentration boundary thickness is the lowest and the concentration maximum value is the smallest for melt with flow. However, in front of downstream dendrite tip, the results are just the opposite. The zone of poor Cu in upstream dendrite where is the most severely microsegregation and shrinkage cavity is wider and the concentration is lower for melt with flow than that without flow.

Phase field method simulation of faceted dendrite growth with arbitrary symmetries

A numerical simulation based on a regularized phase field model is developed to describe faceted dendrite growth morphology. The effects of mesh grid, anisotropy, supersaturation and fold symmetry on dendrite growth morphology were investigated, respectively. These results indicate that the nucleus grows into a hexagonal symmetry faceted dendrite. When the mesh grid is above 640×640, the size has no much effect on the shape. With the increase in the anisotropy value, the tip velocities of faceted dendrite increase and reach a balance value, and then decrease gradually. With the increase in the supersaturation value, crystal evolves from circle to the developed faceted dendrite morphology. Based on the Wulff theory and faceted symmetry morphology diagram, the proposed model was proved to be effective, and it can be generalized to arbitrary crystal symmetries.

RESEARCH ON METHOD TO CALCULATE VELOCITIES OF SOLID PHASE AND LIQUID PHASE IN DEBRIS FLOW

Velocities of solid phase and liquid phase in debris flow are one key problem to research on impact and abrasion mechanism of banks and control structures under action of debris flow. Debris flow was simplified as two-phase liquid composed of solid phase with the same diameter particles and liquid phase with the same mechanical features. Assume debris flow was one-dimension two-phase liquid moving to one direction, then general equations of velocities of solid phase and liquid phase were founded in two-phase theory. Methods to calculate average pressures, volume forces and surface forces of debris flow control volume were established. Specially, surface forces were ascertained using Bingham＇s rheology equation of liquid phase and Bagnold＇s testing results about interaction between particles of solid phase. Proportional coefficient of velocities between liquid phase and solid phase was put forward, meanwhile, divergent coefficient between theoretical velocity and real velocity of solid phase was provided too. To state succinctly before, method to calculate velocities of solid phase and liquid phase was obtained through solution to general equations. The method is suitable for both viscous debris flow and thin debris flow. Additionally, velocities every phase can be identified through analyzing deposits in-situ after occurring of debris flow. It is obvious from engineering case the result in the method is consistent to that in real-time field observation.

Low side lobe pattern synthesis using projection method with genetic algorithm for truncated cone conformal phased arrays

A hybrid method for synthesizing antenna＇s three dimensional （3D） pattern is proposed to obtain the low sidelobe feature of truncated cone conformal phased arrays. In this method, the elements of truncated cone conformal phased arrays are projected to the tangent plane in one generatrix of the truncated cone. Then two dimensional （2D） Chebyshev amplitude distribution optimization is respectively used in two mutual vertical directions of the tangent plane. According to the location of the elements, the excitation current amplitude distribution of each element on the conformal structure is derived reversely, then the excitation current amplitude is further optimized by using the genetic algorithm （GA）. A truncated cone problem with 8x8 elements on it, and a 3D pattern desired side lobe level （SLL） up to 35 dB, is studied. By using the hybrid method, the optimal goal is accomplished with acceptable CPU time, which indicates that this hybrid method for the low sidelobe synthesis is feasible.

Method for array gain and phase uncertainties calibration based on ISM and ESPRIT

A new method for array calibration of array gain and phase uncertainties, which severely degrade the performance of spatial spectrum estimation, is presented. The method is based on the idea of the instrumental sensors method （ISM）, two well-calibrated sensors are added into the original array. By applying the principle of estimation of signal parameters via rotational invariance techniques （ESPRIT）, the direction-of-arrivals （DOAs） and uncertainties can be estimated simultaneously through eigen-decomposition. Compared with the conventional ones, this new method has less computational complexity while has higher estimation precision, what＇s more, it can overcome the problem of ambiguity. Both theoretical analysis and computer simulations show the effectiveness of the proposed method.

Crack propagation simulation in brittle elastic materials by a phase field method

To overcome the difficulties of re-meshing and tracking the crack-tip in other computational methods for crack propagation simulations,the phase field method based on the minimum energy principle is introduced by defining a continuous phase field variable(x)∈[0,1]to characterize discontinuous cracks in brittle materials.This method can well describe the crack initiation and propagation without assuming the shape,size and orientation of the initial crack in advance.In this paper,a phase field method based on Miehe's approach[Miehe et al.,Comp.Meth.App.Mech.Eng.(2010)]is applied to simulate different crack propagation problems in twodimensional(2D),isotropic and linear elastic materials.The numerical implementation of the phase field method is realized within the framework of the finite element method(FEM).The validity,accuracy and efficiency of the present method are verified by comparing the numerical results with other reference results in literature.Several numerical examples are presented to show the effects of the loading type(tension and shear),boundary conditions,and initial crack location and orientation on the crack propagation path and force-displacement curve.Furthermore,for a single edge-cracked bi-material specimen,the influences of the loading type and the crack location on the crack propagation trajectory and force-displacement curve are also investigated and discussed.It is demonstrated that the phase field method is an efficient tool for the numerical simulation of the crack propagation problems in brittle elastic materials,and the corresponding results may have an important relevance for predicting and preventing possible crack propagations in engineering applications.

Preparation and Characterization of Sodium Sulfate/Silica Composite as a Shape-stabilized Phase Change Material by Sol-gel Method

A sodium sulfate （NaeSO4）/silica （SiO2） composite was prepared as a shape-stabilized solid-liquid phase change material by a sol-gel procedure using Na2SiO3 as the silica source. Na2SO4 in the composite acts as a latent heat storage substance for solid-liquid phase change, while SiO2 acts as a support material to provide structural strength and prevent leakage of melted NazSO4. The microstructure and composition of the prepared composite were characterized by the N2 adsorption, transmission electron microscope （TEM）, scanning electron microscope （SEM）, Fourier transform infrared spectroscopy （FTIR） and X-ray diffraction. The results show that the prepared Na2SOJSiO2 composite is a nanostructured hybrid of NazSO4 and SiO2 without new substances produced during the phase change. The macroscopic shape of the NazSO4/SiO2 composite after the melting and freezing cycles does not change and there is no leakage of Na2SO4. Determined by differential scanning calorimeter （DSC） analysis, the values of phase change latent heat of melting and freezing of the prepared NazSO4/SiO2 （50%, by mass） composite are 82.3 kJ.kg i and 83.7 kJ.kg-1, and temperatures of melting and freezing are 886.0 ℃ and 880.6 ℃, respectively. Furthermore, the Na2SOJSiO2 composite maintains good thermal energy storage and release ability even after 100 cycles of melting and freezing. The satisfactory thermal storage performance renders this composite a versatile tool for high-temperature thermal energy storage.

Calculation of Phase Equilibria Based on the Levenberg-Marquardt Method

The Levenberg-Marquardt method, the best algorithm to obtain the least-square solution of nonlinear equations, is applied to calculate the stable phase equilibria. It can get the best combination between robustness and speed of the calculations. Its application to ternary AI-Si-Mg system is executed in detail. The calculated phase equilibria agree well with the experimental results. Furthermore, the Levenberg-Marquardt method is not sensitive to the initial values.