Cyclic loads generated by environmental factors,such as winds,waves,and trains,will likely lead to performance degradation in pile foundations,resulting in issues like permanent displacement accumulation and bearing c...Cyclic loads generated by environmental factors,such as winds,waves,and trains,will likely lead to performance degradation in pile foundations,resulting in issues like permanent displacement accumulation and bearing capacity attenuation.This paper presents a semi-analytical solution for predicting the axial cyclic behavior of piles in sands.The solution relies on two enhanced nonlinear load-transfer models considering stress-strain hysteresis and cyclic degradation in the pile-soil interaction.Model parameters are calibrated through cyclic shear tests of the sand-steel interface and laboratory geotechnical testing of sands.A novel aspect involves the meticulous formulation of the shaft loadtransfer function using an interface constitutive model,which inherently inherits the interface model’s advantages,such as capturing hysteresis,hardening,degradation,and particle breakage.The semi-analytical solution is computed numerically using the matrix displacement method,and the calculated values are validated through model tests performed on non-displacement and displacement piles in sands.The results demonstrate that the predicted values show excellent agreement with the measured values for both the static and cyclic responses of piles in sands.The displacement pile response,including factors such as bearing capacity,mobilized shaft resistance,and convergence rate of permanent settlement,exhibit improvements compared to non-displacement piles attributed to the soil squeezing effect.This methodology presents an innovative analytical framework,allowing for integrating cyclic interface models into the theoretical investigation of pile responses.展开更多
We investigate a periodically driven Haldane model subjected to a two-stage driving scheme in the form of a step function.By using the Floquet theory,we obtain the topological phase diagram of the system.We also find ...We investigate a periodically driven Haldane model subjected to a two-stage driving scheme in the form of a step function.By using the Floquet theory,we obtain the topological phase diagram of the system.We also find that anomalous Floquet topological phases exist in the system.Focusing on examining the quench dynamics among topological phases,we analyze the site distribution of the 0-mode and p-mode edge states in long-period evolution after a quench.The results demonstrate that,under certain conditions,the site distribution of the 0-mode can be confined at the edge even in long-period evolution.Additionally,both the 0-mode and p-mode can recover and become confined at the edge in long-period evolution when the post-quench parameters(T,M_(2) /M_(1))in the phase diagram cross away from the phase boundary (M_(2)/ M_(1))=(6√3t2)/ M_(1)−1.Furthermore,we conclude that whether the edge state is confined at the edge in the long-period evolution after a quench depends on the similarity of the edge states before and after the quench.Our findings reveal some new characteristics of quench dynamics in a periodically driven system.展开更多
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 behavior...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.展开更多
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 ver...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.展开更多
At present,the architecture modeling method of fluvial reservoirs are still developing.Traditional methods usually use grids to characterize architecture interbeds within the reservoir.Due to the thin thickness of thi...At present,the architecture modeling method of fluvial reservoirs are still developing.Traditional methods usually use grids to characterize architecture interbeds within the reservoir.Due to the thin thickness of this type of the interlayers,the number of the model grids must be greatly expanded.The number of grids in the tens of millions often makes an expensive computation;however,upscaling the model will generate a misleading model.The above confusion is the major reason that restricts the largescale industrialization of fluvial reservoir architecture models in oilfield development and production.Therefore,this paper explores an intelligent architecture modeling method for multilevel fluvial reservoirs based on architecture interface and element.Based on the superpositional relationship of different architectural elements within the fluvial reservoir,this method uses a combination of multilevel interface constraints and non-uniform grid techniques to build a high-resolution 3D geological model for reservoir architecture.Through the grid upscaling technology of heterogeneous architecture elements,different upscaling densities are given to the lateral-accretion bedding and lateral-accretion bodies to simplify the model gridding.This new method greatly reduces the number of model grids while ensuring the accuracy of lateral-accretion bedding models,laying a foundation for large-scale numerical simulation of the subsequent industrialization of the architecture model.This method has been validated in A layer of X oilfield with meandering fluvial channel sands as reservoirs and B layer of Y oilfield with braided river sands as reservoirs.The simulation results show that it has a higher accuracy of production history matching and remaining oil distribution forecast of the targeted sand body.The numerical simulation results show that in the actual development process of oilfield,the injected water will not displace oil in a uniform diffusive manner as traditionally assumed,but in a more complex pattern with oil in upper part of sand body being left behind as residual oil due to the influences of different levels of architecture interfaces.This investigation is important to guiding reservoir evaluation,remaining oil analysis,profile control and potential tapping and well pattern adjustment.展开更多
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 predictio...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.展开更多
The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems.Here,we develop a new analytical approach based on the bifur...The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems.Here,we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models.Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions.To elucidate the idea,we apply this analytical approach to three representative phase field equations:the Allen-Cahn equation,the Cahn-Hilliard equation,and the Allen-Cahn-Ohta-Kawasaki system.The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.展开更多
We investigated the properties of the phase diagram of high-order susceptibilities,speed of sound,and polytropic index based on an extended Nambu-Jona-Lasinio model with an eight-quark scalar-vector interaction.Non-mo...We investigated the properties of the phase diagram of high-order susceptibilities,speed of sound,and polytropic index based on an extended Nambu-Jona-Lasinio model with an eight-quark scalar-vector interaction.Non-monotonic behavior was observed in all these quantities around the phase transition boundary,which also revealed the properties of the critical point.Further,this study indicated that the chiral phase transition boundary and critical point could vary depending on the scalarvector coupling constant G_(SV).At finite densities and temperatures,the negative G_(SV)term exhibited attractive interactions,which enhanced the critical point temperature and reduced the chemical potential.The G_(SV)term also affected the properties of the high-order susceptibilities,speed of sound,and polytropic index near the critical point.The non-monotonic(peak or dip)structures of these quantities shifted to a low baryon chemical potential(and high temperature)with a negative G_(SV).G_(SV)also changed the amplitude and range of the nonmonotonic regions.Therefore,the scalar-vector interaction was useful for locating the phase boundary and critical point in QCD phase diagram by comparing the experimental data.The study of the non-monotonic behavior of high-order susceptibilities,speed of sound,and polytropic index is of great interest,and further observations related to high-order susceptibilities,speed of sound,and polytropic index being found and applied to the search for critical points in heavy-ion collisions and the study of compact stars are eagerly awaited.展开更多
In contrast to conventional transformers, power electronic transformers, as an integral component of new energy power system, are often subjected to high-frequency and transient electrical stresses, leading to heighte...In contrast to conventional transformers, power electronic transformers, as an integral component of new energy power system, are often subjected to high-frequency and transient electrical stresses, leading to heightened concerns regarding insulation failures. Meanwhile, the underlying mechanism behind discharge breakdown failure and nanofiller enhancement under high-frequency electrical stress remains unclear. An electric-thermal coupled discharge breakdown phase field model was constructed to study the evolution of the breakdown path in polyimide nanocomposite insulation subjected to high-frequency stress. The investigation focused on analyzing the effect of various factors, including frequency, temperature, and nanofiller shape, on the breakdown path of Polyimide(PI) composites. Additionally, it elucidated the enhancement mechanism of nano-modified composite insulation at the mesoscopic scale. The results indicated that with increasing frequency and temperature, the discharge breakdown path demonstrates accelerated development, accompanied by a gradual dominance of Joule heat energy. This enhancement is attributed to the dispersed electric field distribution and the hindering effect of the nanosheets. The research findings offer a theoretical foundation and methodological framework to inform the optimal design and performance management of new insulating materials utilized in high-frequency power equipment.展开更多
An in-depth investigation is made on the problem of an arc-shapedinterface insulating crack in a three-phase concentric circularcylindrical piezoelectric composite constitutive model. An exactsolution in series form i...An in-depth investigation is made on the problem of an arc-shapedinterface insulating crack in a three-phase concentric circularcylindrical piezoelectric composite constitutive model. An exactsolution in series form is derived by employing the complex variablemethod. In addition, the distribution of physical quantities such asstresses, strains, electric displacements and electric fields in thewhole field and along the interface is also presented.展开更多
A novel approach based on the quantitative phase field model was proposed to calculate the interface mobility and applied to the α/β interface of a ternary Ti-6Al-4V alloy.Phase field simulations indicate that the h...A novel approach based on the quantitative phase field model was proposed to calculate the interface mobility and applied to the α/β interface of a ternary Ti-6Al-4V alloy.Phase field simulations indicate that the higher interface mobility leads to the faster transformation rate,but only a unique value of interface mobility matches the diffusion equation under the diffusion-controlled condition.By comparing the transformation kinetics from phase field simulations with that from classical diffusion equation,the interface mobility at different temperatures can be obtained.The results show that the calculated interface mobility increases with increasing temperature and accords with Arrhenius equation very well.展开更多
A thermodynamic consistent phase field model is developed to describe the sintering process with multiphase powders. In this model, the interface region is assumed to be a mixture of different phases with the same che...A thermodynamic consistent phase field model is developed to describe the sintering process with multiphase powders. In this model, the interface region is assumed to be a mixture of different phases with the same chemical potential, but with different compositions. The interface diffusion and boundary diffusion are also considered in the model. As an example, the model is applied to the sintering process with Fe-Cu powders. The free energy of each phase is described by the well-developed thermodynamic models, together with the published optimized parameters. The microstructure and solute distribution during the sintering process can both be obtained quantitively.展开更多
The early precipitation process of Ni(75)Al(14)Mo(11) alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and oc...The early precipitation process of Ni(75)Al(14)Mo(11) alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and occupation probability of the three kinds of atoms were investigated.It is indicated that the non-stoichiometric Ll0(Ⅰ/Ⅱ) phases are found in the precipitation process.With the temperature increasing,the appearance time of Ll0 is brought forward.The Ll0(Ⅱ) structure always precipitates earlier than the Ll0(Ⅰ) structure.Compared with lower temperature,higher temperature brings the formation time of Ll0 phase forward and makes Ll0 phase have a higher order degree.But lower temperature shortens the process time of the Ll0 phase to the Ll2 phase.Al and Mo atoms tend to occupy γ site,Ni atom tends to occupy a and β sites.At the same temperature,Al atom has stronger occupation ability than Mo atom in the same site.Ni,Al and Mo collectively form the composited Ll2 structure.展开更多
Scanning electron microscopy and X-ray energy dispersive spectrum analysis show that the clusters of intermetallic AlFeSi particle are distributed on or near the aluminum foil stock surfaces heterogeneously. 3D finite...Scanning electron microscopy and X-ray energy dispersive spectrum analysis show that the clusters of intermetallic AlFeSi particle are distributed on or near the aluminum foil stock surfaces heterogeneously. 3D finite element modeling shows that these clusters of hard particles induce the fracture of the nano-scale lubricant oil film at first and further lead to severe deformation in the nearby aluminum foil substrate along the rolling direction. Consequently, the optical property in this region differs from that in the surroundings, resulting in surface defects.展开更多
Building a model for a object through model composition is a interesting topics, this paper research the interface composit ion of models in the Zeigler’s modeling methodology DEVS.
The behavior of soil-structure interface plays a major role in the definition of soil-structure interaction. In this paper a bi-potential surface elasto-plastic model for soil-structure interface is proposed in order ...The behavior of soil-structure interface plays a major role in the definition of soil-structure interaction. In this paper a bi-potential surface elasto-plastic model for soil-structure interface is proposed in order to describe the interface deformation behavior,including strain softening and normal dilatancy. The model is formulated in the framework of generalized potential theory,in which the soil-structure interface problem is regard as a two-dimensional mathematical problem in stress field,and plastic state equations are used to replace the traditional field surface. The relation curves of shear stress and tangential strain are fitted by a piecewise function composed by hyperbolic functions and hyperbolic secant functions,while the relation curves of normal strain and tangential strain are fitted by another piecewise function composed by quadratic functions and hyperbolic secant functions. The approach proposed has the advantage of deriving an elastoplastic constitutive matrix without postulating the plastic potential functions and yield surface. Moreover,the mathematical principle is clear,and the entire model parameters can be identified by experimental tests. Finally,the predictions of the model have been compared with experimental results obtained from simple shear tests under normal stresses,and results show the model is reasonable and practical.展开更多
Thixocasting requires manufacturing of billets with non-dendritic microstructure.Aluminum alloy A356 billets were produced by rheocasting in a mould placed inside a linear electromagnetic stirrer.Subsequent heat treat...Thixocasting requires manufacturing of billets with non-dendritic microstructure.Aluminum alloy A356 billets were produced by rheocasting in a mould placed inside a linear electromagnetic stirrer.Subsequent heat treatment was used to produce a transition from rosette to globular microstructure.The current and the duration of stirring were explored as control parameters.Simultaneous induction heating of the billet during stirring was quantified using experimentally determined thermal profiles.The effect of processing parameters on the dendrite fragmentation was discussed.Corresponding computational modeling of the process was performed using phase-field modeling of alloy solidification in order to gain insight into the process of morphological changes of a solid during this process.A non-isothermal alloy solidification model was used for simulations.The morphological evolution under such imposed thermal cycles was simulated and compared with experimentally determined one.Suitable scaling using the thermosolutal diffusion distances was used to overcome computational difficulties in quantitative comparison at system scale.The results were interpreted in the light of existing theories of microstructure refinement and globularisation.展开更多
This article introduces an element diffusion behavior model for a titanium/steel explosive clad plate characterized by a typical curved interface during the heat-treatment process. A series of heat-treatment experimen...This article introduces an element diffusion behavior model for a titanium/steel explosive clad plate characterized by a typical curved interface during the heat-treatment process. A series of heat-treatment experiments were conducted in the temperature range from 750℃ to 950℃, and the effects of heat-treatment parameters on the microstructural evolution and diffusion behavior were investigated by optical microscopy, scanning electron microscopy, X-ray diffraction analysis, and electron-probe microanalysis. Carbon atoms within the steel matrix were observed to diffuse toward the titanium matrix and to aggregate at the bonding interface at 850℃ or lower; in contrast, when the temperature exceeded 850℃, the mutual diffusion of Ti and Fe occurred, along with the diffusion of C atoms, resulting in the for- marion of Ti-Fe intermetallics (Fe2Ti/FeTi). The diffusion distances of C, Ti, and Fe atoms increased with increasing heating temperature and/or holding time. On the basis of this diffusion behavior, a novel diffusion model was proposed. This model considers the effects of various factors, including the curvature radius of the curved interface, the diffusion coefficient, the heating temperature, and the holding rime. The experimental results show good agreement with the calculated values. The proposed model could clearly provide a general prediction of the elements' diffusion at both straight and curved interfaces.展开更多
Several methods representing the evolution of microstructure were introduced, which include the Johnson-Mehl-Avrami-Kolmogorov (JMAK) equation, Internal State Variable (ISV) framework, Koistinen-Marburger (K-M) ...Several methods representing the evolution of microstructure were introduced, which include the Johnson-Mehl-Avrami-Kolmogorov (JMAK) equation, Internal State Variable (ISV) framework, Koistinen-Marburger (K-M) equation, modified Magee's rule and phase field model, etc. By combining calculation of martensite transformation kinetics, considering the selection of parameters with the effect of austenite grain size (AGS), some suitable ways of obtaining better results have been proposed.展开更多
Due to the complexity of soil-structure interaction, simple constitutive models typically used for interface elements in general computer programs cannot satisfy the requirements of discontinuous deformation analysis ...Due to the complexity of soil-structure interaction, simple constitutive models typically used for interface elements in general computer programs cannot satisfy the requirements of discontinuous deformation analysis of structures that contain different interfaces. In order to simulate the strain-softening characteristics of interfaces, a nonlinear strain-softening interface constitutive model was incorporated into fast Lagrange analysis of continua in three dimensions (FLAC3D) through a user-defined program in the FISH environment. A numerical simulation of a direct shear test for geosynthetic interfaces was conducted to verify that the interface model was implemented correctly. Results of the numerical tests show good agreement with the results obtained from theoretical calculations, indicating that the model incorporated into FLAC3D can simulate the nonlinear strain-softening behavior of interfaces involving geosynthetic materials. The results confirmed the validity and reliability of the improved interface model. The procedure and method of implementing an interface constitutive model into a commercial computer program also provide a reference for implementation of a new interface constitutive model in FLAC3D.展开更多
基金the financial support provided by the National Natural Science Foundation of China(Grant No.42272310).
文摘Cyclic loads generated by environmental factors,such as winds,waves,and trains,will likely lead to performance degradation in pile foundations,resulting in issues like permanent displacement accumulation and bearing capacity attenuation.This paper presents a semi-analytical solution for predicting the axial cyclic behavior of piles in sands.The solution relies on two enhanced nonlinear load-transfer models considering stress-strain hysteresis and cyclic degradation in the pile-soil interaction.Model parameters are calibrated through cyclic shear tests of the sand-steel interface and laboratory geotechnical testing of sands.A novel aspect involves the meticulous formulation of the shaft loadtransfer function using an interface constitutive model,which inherently inherits the interface model’s advantages,such as capturing hysteresis,hardening,degradation,and particle breakage.The semi-analytical solution is computed numerically using the matrix displacement method,and the calculated values are validated through model tests performed on non-displacement and displacement piles in sands.The results demonstrate that the predicted values show excellent agreement with the measured values for both the static and cyclic responses of piles in sands.The displacement pile response,including factors such as bearing capacity,mobilized shaft resistance,and convergence rate of permanent settlement,exhibit improvements compared to non-displacement piles attributed to the soil squeezing effect.This methodology presents an innovative analytical framework,allowing for integrating cyclic interface models into the theoretical investigation of pile responses.
基金the National Natural Science Foundation of China(Grant No.12004049).
文摘We investigate a periodically driven Haldane model subjected to a two-stage driving scheme in the form of a step function.By using the Floquet theory,we obtain the topological phase diagram of the system.We also find that anomalous Floquet topological phases exist in the system.Focusing on examining the quench dynamics among topological phases,we analyze the site distribution of the 0-mode and p-mode edge states in long-period evolution after a quench.The results demonstrate that,under certain conditions,the site distribution of the 0-mode can be confined at the edge even in long-period evolution.Additionally,both the 0-mode and p-mode can recover and become confined at the edge in long-period evolution when the post-quench parameters(T,M_(2) /M_(1))in the phase diagram cross away from the phase boundary (M_(2)/ M_(1))=(6√3t2)/ M_(1)−1.Furthermore,we conclude that whether the edge state is confined at the edge in the long-period evolution after a quench depends on the similarity of the edge states before and after the quench.Our findings reveal some new characteristics of quench dynamics in a periodically driven system.
基金the National Natural Science Foundation of China(Grant No.62273033).
文摘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.
基金Project supported by the Open Project of the Key Laboratory of Xinjiang Uygur Autonomous Region,China(Grant No.2021D04015)the Yili Kazakh Autonomous Prefecture Science and Technology Program Project,China(Grant No.YZ2022B021).
文摘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.
文摘At present,the architecture modeling method of fluvial reservoirs are still developing.Traditional methods usually use grids to characterize architecture interbeds within the reservoir.Due to the thin thickness of this type of the interlayers,the number of the model grids must be greatly expanded.The number of grids in the tens of millions often makes an expensive computation;however,upscaling the model will generate a misleading model.The above confusion is the major reason that restricts the largescale industrialization of fluvial reservoir architecture models in oilfield development and production.Therefore,this paper explores an intelligent architecture modeling method for multilevel fluvial reservoirs based on architecture interface and element.Based on the superpositional relationship of different architectural elements within the fluvial reservoir,this method uses a combination of multilevel interface constraints and non-uniform grid techniques to build a high-resolution 3D geological model for reservoir architecture.Through the grid upscaling technology of heterogeneous architecture elements,different upscaling densities are given to the lateral-accretion bedding and lateral-accretion bodies to simplify the model gridding.This new method greatly reduces the number of model grids while ensuring the accuracy of lateral-accretion bedding models,laying a foundation for large-scale numerical simulation of the subsequent industrialization of the architecture model.This method has been validated in A layer of X oilfield with meandering fluvial channel sands as reservoirs and B layer of Y oilfield with braided river sands as reservoirs.The simulation results show that it has a higher accuracy of production history matching and remaining oil distribution forecast of the targeted sand body.The numerical simulation results show that in the actual development process of oilfield,the injected water will not displace oil in a uniform diffusive manner as traditionally assumed,but in a more complex pattern with oil in upper part of sand body being left behind as residual oil due to the influences of different levels of architecture interfaces.This investigation is important to guiding reservoir evaluation,remaining oil analysis,profile control and potential tapping and well pattern adjustment.
基金Project supported by the Hefei National Research Center for Physical Sciences at the Microscale (Grant No.KF2021002)the Natural Science Foundation of Shanxi Province,China (Grant Nos.202303021221029 and 202103021224051)+2 种基金the National Natural Science Foundation of China (Grant Nos.11975024,12047503,and 12275263)the Anhui Provincial Supporting Program for Excellent Young Talents in Colleges and Universities (Grant No.gxyq ZD2019023)the National Key Research and Development Program of China (Grant No.2018YFA0306501)。
文摘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.
基金supported as part of the Computational Materials Sciences Program funded by the U.S.Department of Energy,Office of Science,Basic Energy Sciences,under Award No.DE-SC0020145Y.Z.would like to acknowledge support for his effort by the Simons Foundation through Grant No.357963 and NSF grant DMS-2142500.
文摘The phase field method is playing an increasingly important role in understanding and predicting morphological evolution in materials and biological systems.Here,we develop a new analytical approach based on the bifurcation analysis to explore the mathematical solution structure of phase field models.Revealing such solution structures not only is of great mathematical interest but also may provide guidance to experimentally or computationally uncover new morphological evolution phenomena in materials undergoing electronic and structural phase transitions.To elucidate the idea,we apply this analytical approach to three representative phase field equations:the Allen-Cahn equation,the Cahn-Hilliard equation,and the Allen-Cahn-Ohta-Kawasaki system.The solution structures of these three phase field equations are also verified numerically by the homotopy continuation method.
基金supported by the National Natural Science Foundation of China(Nos.12205158 and 11975132)the Shandong Provincial Natural Science Foundation,China(Nos.ZR2021QA037,ZR2022JQ04 and ZR2019YQ01)。
文摘We investigated the properties of the phase diagram of high-order susceptibilities,speed of sound,and polytropic index based on an extended Nambu-Jona-Lasinio model with an eight-quark scalar-vector interaction.Non-monotonic behavior was observed in all these quantities around the phase transition boundary,which also revealed the properties of the critical point.Further,this study indicated that the chiral phase transition boundary and critical point could vary depending on the scalarvector coupling constant G_(SV).At finite densities and temperatures,the negative G_(SV)term exhibited attractive interactions,which enhanced the critical point temperature and reduced the chemical potential.The G_(SV)term also affected the properties of the high-order susceptibilities,speed of sound,and polytropic index near the critical point.The non-monotonic(peak or dip)structures of these quantities shifted to a low baryon chemical potential(and high temperature)with a negative G_(SV).G_(SV)also changed the amplitude and range of the nonmonotonic regions.Therefore,the scalar-vector interaction was useful for locating the phase boundary and critical point in QCD phase diagram by comparing the experimental data.The study of the non-monotonic behavior of high-order susceptibilities,speed of sound,and polytropic index is of great interest,and further observations related to high-order susceptibilities,speed of sound,and polytropic index being found and applied to the search for critical points in heavy-ion collisions and the study of compact stars are eagerly awaited.
基金supported in part by the National Key R&D Program of China (No.2021YFB2601404)Beijing Natural Science Foundation (No.3232053)National Natural Science Foundation of China (Nos.51929701 and 52127812)。
文摘In contrast to conventional transformers, power electronic transformers, as an integral component of new energy power system, are often subjected to high-frequency and transient electrical stresses, leading to heightened concerns regarding insulation failures. Meanwhile, the underlying mechanism behind discharge breakdown failure and nanofiller enhancement under high-frequency electrical stress remains unclear. An electric-thermal coupled discharge breakdown phase field model was constructed to study the evolution of the breakdown path in polyimide nanocomposite insulation subjected to high-frequency stress. The investigation focused on analyzing the effect of various factors, including frequency, temperature, and nanofiller shape, on the breakdown path of Polyimide(PI) composites. Additionally, it elucidated the enhancement mechanism of nano-modified composite insulation at the mesoscopic scale. The results indicated that with increasing frequency and temperature, the discharge breakdown path demonstrates accelerated development, accompanied by a gradual dominance of Joule heat energy. This enhancement is attributed to the dispersed electric field distribution and the hindering effect of the nanosheets. The research findings offer a theoretical foundation and methodological framework to inform the optimal design and performance management of new insulating materials utilized in high-frequency power equipment.
基金the National Natural Science Foundation of China(No.59635140)the Doctorate Foundation of Xi'an Jiaotong University
文摘An in-depth investigation is made on the problem of an arc-shapedinterface insulating crack in a three-phase concentric circularcylindrical piezoelectric composite constitutive model. An exactsolution in series form is derived by employing the complex variablemethod. In addition, the distribution of physical quantities such asstresses, strains, electric displacements and electric fields in thewhole field and along the interface is also presented.
基金Project (51101059) supported by the National Natural Science Foundation of ChinaProject (20110490874) supported by the China Postdoctoral Science Foundation
文摘A novel approach based on the quantitative phase field model was proposed to calculate the interface mobility and applied to the α/β interface of a ternary Ti-6Al-4V alloy.Phase field simulations indicate that the higher interface mobility leads to the faster transformation rate,but only a unique value of interface mobility matches the diffusion equation under the diffusion-controlled condition.By comparing the transformation kinetics from phase field simulations with that from classical diffusion equation,the interface mobility at different temperatures can be obtained.The results show that the calculated interface mobility increases with increasing temperature and accords with Arrhenius equation very well.
基金Project(2011CB606306)supported by the National Basic Research Program of ChinaProject(51101014)supported by the National Natural Science Foundation of ChinaProject(SKLSP201214)supported by the State Key Laboratory of Solidification Processing in Northwestern Polytechnical University,China
文摘A thermodynamic consistent phase field model is developed to describe the sintering process with multiphase powders. In this model, the interface region is assumed to be a mixture of different phases with the same chemical potential, but with different compositions. The interface diffusion and boundary diffusion are also considered in the model. As an example, the model is applied to the sintering process with Fe-Cu powders. The free energy of each phase is described by the well-developed thermodynamic models, together with the published optimized parameters. The microstructure and solute distribution during the sintering process can both be obtained quantitively.
基金Project(51275486)supported by the National Natural Science Foundation of China
文摘The early precipitation process of Ni(75)Al(14)Mo(11) alloy was simulated by microscopic phase-field model at different temperatures.The microstructure of the alloy,the precipitation time of Llo structure and occupation probability of the three kinds of atoms were investigated.It is indicated that the non-stoichiometric Ll0(Ⅰ/Ⅱ) phases are found in the precipitation process.With the temperature increasing,the appearance time of Ll0 is brought forward.The Ll0(Ⅱ) structure always precipitates earlier than the Ll0(Ⅰ) structure.Compared with lower temperature,higher temperature brings the formation time of Ll0 phase forward and makes Ll0 phase have a higher order degree.But lower temperature shortens the process time of the Ll0 phase to the Ll2 phase.Al and Mo atoms tend to occupy γ site,Ni atom tends to occupy a and β sites.At the same temperature,Al atom has stronger occupation ability than Mo atom in the same site.Ni,Al and Mo collectively form the composited Ll2 structure.
基金Project(51074117)supported by the National Natural Science Foundation of ChinaProject(2009CDA044)supported by the Foundation for Distinguished Young Scientists of Hubei Province,ChinaProjects(201104493,20100471161)supported by the China Postdoctoral Science Foundation
文摘Scanning electron microscopy and X-ray energy dispersive spectrum analysis show that the clusters of intermetallic AlFeSi particle are distributed on or near the aluminum foil stock surfaces heterogeneously. 3D finite element modeling shows that these clusters of hard particles induce the fracture of the nano-scale lubricant oil film at first and further lead to severe deformation in the nearby aluminum foil substrate along the rolling direction. Consequently, the optical property in this region differs from that in the surroundings, resulting in surface defects.
文摘Building a model for a object through model composition is a interesting topics, this paper research the interface composit ion of models in the Zeigler’s modeling methodology DEVS.
基金supported by the National Natural Science Foundation of ChinaYalona River Hydropower Development of Ertan Hydropower Development Company (No.50639050)
文摘The behavior of soil-structure interface plays a major role in the definition of soil-structure interaction. In this paper a bi-potential surface elasto-plastic model for soil-structure interface is proposed in order to describe the interface deformation behavior,including strain softening and normal dilatancy. The model is formulated in the framework of generalized potential theory,in which the soil-structure interface problem is regard as a two-dimensional mathematical problem in stress field,and plastic state equations are used to replace the traditional field surface. The relation curves of shear stress and tangential strain are fitted by a piecewise function composed by hyperbolic functions and hyperbolic secant functions,while the relation curves of normal strain and tangential strain are fitted by another piecewise function composed by quadratic functions and hyperbolic secant functions. The approach proposed has the advantage of deriving an elastoplastic constitutive matrix without postulating the plastic potential functions and yield surface. Moreover,the mathematical principle is clear,and the entire model parameters can be identified by experimental tests. Finally,the predictions of the model have been compared with experimental results obtained from simple shear tests under normal stresses,and results show the model is reasonable and practical.
文摘Thixocasting requires manufacturing of billets with non-dendritic microstructure.Aluminum alloy A356 billets were produced by rheocasting in a mould placed inside a linear electromagnetic stirrer.Subsequent heat treatment was used to produce a transition from rosette to globular microstructure.The current and the duration of stirring were explored as control parameters.Simultaneous induction heating of the billet during stirring was quantified using experimentally determined thermal profiles.The effect of processing parameters on the dendrite fragmentation was discussed.Corresponding computational modeling of the process was performed using phase-field modeling of alloy solidification in order to gain insight into the process of morphological changes of a solid during this process.A non-isothermal alloy solidification model was used for simulations.The morphological evolution under such imposed thermal cycles was simulated and compared with experimentally determined one.Suitable scaling using the thermosolutal diffusion distances was used to overcome computational difficulties in quantitative comparison at system scale.The results were interpreted in the light of existing theories of microstructure refinement and globularisation.
文摘This article introduces an element diffusion behavior model for a titanium/steel explosive clad plate characterized by a typical curved interface during the heat-treatment process. A series of heat-treatment experiments were conducted in the temperature range from 750℃ to 950℃, and the effects of heat-treatment parameters on the microstructural evolution and diffusion behavior were investigated by optical microscopy, scanning electron microscopy, X-ray diffraction analysis, and electron-probe microanalysis. Carbon atoms within the steel matrix were observed to diffuse toward the titanium matrix and to aggregate at the bonding interface at 850℃ or lower; in contrast, when the temperature exceeded 850℃, the mutual diffusion of Ti and Fe occurred, along with the diffusion of C atoms, resulting in the for- marion of Ti-Fe intermetallics (Fe2Ti/FeTi). The diffusion distances of C, Ti, and Fe atoms increased with increasing heating temperature and/or holding time. On the basis of this diffusion behavior, a novel diffusion model was proposed. This model considers the effects of various factors, including the curvature radius of the curved interface, the diffusion coefficient, the heating temperature, and the holding rime. The experimental results show good agreement with the calculated values. The proposed model could clearly provide a general prediction of the elements' diffusion at both straight and curved interfaces.
基金Item Sponsored by State High Technology Research and Development Program(863 Plan) of China (2001AA332020)
文摘Several methods representing the evolution of microstructure were introduced, which include the Johnson-Mehl-Avrami-Kolmogorov (JMAK) equation, Internal State Variable (ISV) framework, Koistinen-Marburger (K-M) equation, modified Magee's rule and phase field model, etc. By combining calculation of martensite transformation kinetics, considering the selection of parameters with the effect of austenite grain size (AGS), some suitable ways of obtaining better results have been proposed.
基金supported by the National Natural Science Foundation of China (Grant No.51079047)
文摘Due to the complexity of soil-structure interaction, simple constitutive models typically used for interface elements in general computer programs cannot satisfy the requirements of discontinuous deformation analysis of structures that contain different interfaces. In order to simulate the strain-softening characteristics of interfaces, a nonlinear strain-softening interface constitutive model was incorporated into fast Lagrange analysis of continua in three dimensions (FLAC3D) through a user-defined program in the FISH environment. A numerical simulation of a direct shear test for geosynthetic interfaces was conducted to verify that the interface model was implemented correctly. Results of the numerical tests show good agreement with the results obtained from theoretical calculations, indicating that the model incorporated into FLAC3D can simulate the nonlinear strain-softening behavior of interfaces involving geosynthetic materials. The results confirmed the validity and reliability of the improved interface model. The procedure and method of implementing an interface constitutive model into a commercial computer program also provide a reference for implementation of a new interface constitutive model in FLAC3D.