Water electrolyzers play a crucial role in green hydrogen production.However,their efficiency and scalability are often compromised by bubble dynamics across various scales,from nanoscale to macroscale components.This...Water electrolyzers play a crucial role in green hydrogen production.However,their efficiency and scalability are often compromised by bubble dynamics across various scales,from nanoscale to macroscale components.This review explores multi-scale modeling as a tool to visualize multi-phase flow and improve mass transport in water electrolyzers.At the nanoscale,molecular dynamics(MD)simulations reveal how electrode surface features and wettability influence nanobubble nucleation and stability.Moving to the mesoscale,models such as volume of fluid(VOF)and lattice Boltzmann method(LBM)shed light on bubble transport in porous transport layers(PTLs).These insights inform innovative designs,including gradient porosity and hydrophilic-hydrophobic patterning,aimed at minimizing gas saturation.At the macroscale,VOF simulations elucidate two-phase flow regimes within channels,showing how flow field geometry and wettability affect bubble discharging.Moreover,artificial intelligence(AI)-driven surrogate models expedite the optimization process,allowing for rapid exploration of structural parameters in channel-rib flow fields and porous flow field designs.By integrating these approaches,we can bridge theoretical insights with experimental validation,ultimately enhancing water electrolyzer performance,reducing costs,and advancing affordable,high-efficiency hydrogen production.展开更多
Based on the theory of superimposed deformation and the regional tectonic background,the multi-phase non-coaxial superimposed structures in Junggar Basin were systematically analyzed using seismic interpretation,field...Based on the theory of superimposed deformation and the regional tectonic background,the multi-phase non-coaxial superimposed structures in Junggar Basin were systematically analyzed using seismic interpretation,field outcrop observation,and paleo-stress field recovery methods according to the characteristics of the current tectonic framework.Moreover,the tectonic evolution process of the basin was reconstructed using sandbox analogue modelling technology.The results showed that the study area has experienced five phases of non-coaxial deformation with superimposition:The first phase of deformation(D1)is characterized by NNE-SSW extension during late Carboniferous to early Permian,which formed large graben,half graben and other extensional structure style around the basin.The second phase of deformation(D2)is represented by NE-SW compression during the middle to late Permian,and it comprised numerous contraction structures that developed based on D1.The basic form of the entire basin is alternating uplift and depression.The third phase of deformation(D3)is the NW-SE transpressional strike-slip in the Triassic-Jurassic,which produced numerous strike-slip structural styles in the middle part of the basin.The fourth phase of deformation(D4)is the uniform sedimentation during Cretaceous,and the fifth phase(D5)is the compression along NNE-SSW due to the North Tianshan northward thrust,which produced three rows of fold thrust belts and tear faults in the front of the mountain in the southern margin of the basin.The newly established three-dimensional tectonic evolution model shows that,based on the large number of NW-trending grabens and half grabens in the Carboniferous basement of Junggar Basin,multiple level NE trending uplifts have formed with the joint superposition of the late structural inversion and multiple stress fields.This has resulted in the current tectonic units of alternating uplifts and depressions in different directions in the study area.展开更多
A thermodynamically complete multi-phase equation of state(EOS)applicable to both dense and porous metals at wide ranges of temperature and pressure is constructed.A standard three-term decomposition of the Helmholtz ...A thermodynamically complete multi-phase equation of state(EOS)applicable to both dense and porous metals at wide ranges of temperature and pressure is constructed.A standard three-term decomposition of the Helmholtz free energy as a function of specific volume and temperature is presented,where the cold component models both compression and expansion states,the thermal ion component introduces the Debye approximation and melting entropy,and the thermal electron component employs the Thomas-Fermi-Kirzhnits(TFK)model.The porosity of materials is considered by introducing the dynamic porosity coefficientαand the constitutive P-αrelation,connecting the thermodynamic properties between dense and porous systems,allowing for an accurate description of the volume decrease caused by void collapse while maintaining the quasi-static thermodynamic properties of porous systems identical to the dense ones.These models enable the EOS applicable and robust at wide ranges of temperature,pressure and porosity.A systematic evaluation of the new EOS is conducted with aluminum(Al)as an example.300 K isotherm,shock Hugoniot,as well as melting curves of both dense and porous Al are calculated,which shows great agreements with experimental data and validates the effectiveness of the models and the accuracy of parameterizations.Notably,it is for the first time Hugoniot P-σcurves up to 10~6 GPa and shock melting behaviors of porous Al are derived from analytical EOS models,which predict much lower compression limit and shock melting temperatures than those of dense Al.展开更多
Based on the fundamental equations of the mechanics of solid continuum, the paper employs an analytical model for determination of elastic thermal stresses in isotropic continuum represented by periodically distribute...Based on the fundamental equations of the mechanics of solid continuum, the paper employs an analytical model for determination of elastic thermal stresses in isotropic continuum represented by periodically distributed spherical particles with different distributions in an infinite matrix, imaginarily divided into identical cells with dimensions equal to inter-particle distances, containing a central spherical particle with or without a spherical envelope on the particle surface. Consequently, the multi-particle-（envelope）- matrix system, as a model system regarding the analytical modelling, is applicable to four types of multi-phase materials. As functions of the particle volume fraction v, the inter-particle distances dl, d2, d3 along three mutually per- pendicular axes, and the particle and envelope radii, R1 and R2, respectively, the thermal stresses within the cell, are originated during a cooling process as a consequence of the difference in thermal expansion coefficients of phases rep- resented by the matrix, envelope and particle. Analytical-（experimental）-computational lifetime prediction methods for multi-phase materials are proposed, which can be used in engineering with appropriate values of parameters of real multi-phase materials.展开更多
A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection condit...A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection conditions by coupling a multi-phase field model and a lattice Boltzmann model.The elimination mechanism in the evolution process is analyzed for the three conformational schemes constituting converging bi-crystals under pure diffusion and forced convection conditions,respectively,expanding the research of the competitive growth of columnar dendrites under melt convection conditions.The results show that the elimination mechanism for the competitive growth of converging bi-crystals of all three configurations under pure diffusion conditions follows the conventional Walton-Chalmers model.When there is forced convection with lateral flow in the liquid phase,the anomalous elimination phenomenon of unfavorable dendrites eliminating favorable dendrites occurs in the grain boundaries.In particular,the anomalous elimination phenomenon is relatively strong in conformation 1 and conformation 2 when the orientation angle of unfavorable dendrites is small,and relatively weak in conformation 3.Moreover,the presence of convection increases the tip growth rate of both favorable and unfavorable dendrites in the grain boundary.In addition,the parallelization of the multi-phase-field-lattice Boltzmann model is achieved by designing the parallel computation of the model on the GPU platform concerning the computerunified-device-architecture parallel technique,and the results show that the parallel computation of this model based on the GPU has absolute advantages,and the parallel acceleration is more obvious as the computation area increases.展开更多
To determine the time-independent constitutive modeling for porous and multi- phase nanocrystalline materials and understand the effects of grain size and porosity on their mechanical behavior, each phase was treated ...To determine the time-independent constitutive modeling for porous and multi- phase nanocrystalline materials and understand the effects of grain size and porosity on their mechanical behavior, each phase was treated as a mixture of grain interior and grain bound- ary, and pores were taken as a single phase, then Budiansky＇s self-consistent method was used to calculate the Young＇s modulus of porous, possible multi-phase, nanocrystalline materials, the prediction being in good agreement with the results in the literature. Further, the established method is extended to simulate the constitutive relations of porous and possible multi-phase nanocrystalline materials with small plastic deformation in conjunction with the secant-moduli approach and iso-strain assumption. Comparisons between the experimental grain size and porosity dependent mechanical data and the corresponding predictions using the established model show that it appears to be capable of describing the time-independent mechanical behaviors for porous and multi-phase nanocrystalline materials in a small plastic strain range. Further discussion on the modification factor, the advantages and limitations of the model developed were present.展开更多
The multi-phase field model of grain competitive growth during directional solidification of alloy is established.Solving multi-phase field models for thin interface layer thickness conditions,the grain boundary evolu...The multi-phase field model of grain competitive growth during directional solidification of alloy is established.Solving multi-phase field models for thin interface layer thickness conditions,the grain boundary evolution and grain elimination during the competitive growth of SCN-0.24-wt%camphor model alloy bi-crystals are investigated.The effects of different crystal orientations and pulling velocities on grain boundary microstructure evolution are quantitatively analyzed.The obtained results are shown below.In the competitive growth of convergent bi-crystals,when favorably oriented dendrites are in the same direction as the heat flow and the pulling speed is too large,the orientation angle of the bi-crystal from small to large size is the normal elimination phenomenon of the favorably oriented dendrite,blocking the unfavorably oriented dendrite,and the grain boundary is along the growth direction of the favorably oriented dendrite.When the pulling speed becomes small,the grain boundary shows the anomalous elimination phenomenon of the unfavorably oriented dendrite,eliminating the favorably oriented dendrite.In the process of competitive growth of divergent bi-crystal,when the growth direction of favorably oriented dendrites is the same as the heat flow direction and the orientation angle of unfavorably oriented grains is small,the frequency of new spindles of favorably oriented grains is significantly higher than that of unfavorably oriented grains,and as the orientation angle of unfavorably oriented dendrites becomes larger,the unfavorably oriented grains are more likely to have stable secondary dendritic arms,which in turn develop new primary dendritic arms to occupy the liquid phase grain boundary space,but the grain boundary direction is still parallel to favorably oriented dendrites.In addition,the tertiary dendritic arms on the developed secondary dendritic arms may also be blocked by the surrounding lateral branches from further developing into nascent main axes,this blocking of the tertiary dendritic arms has a random nature,which can have aninfluence on the generation of nascent primary main axes in the grain boundaries.展开更多
Smoothed particle hydrodynamics(SPH)is a Lagrangian,mesh-free numerical method renowned for its ability to handle fluid dynamics problems with large interface deformations and multiphase flow coupling.This study intro...Smoothed particle hydrodynamics(SPH)is a Lagrangian,mesh-free numerical method renowned for its ability to handle fluid dynamics problems with large interface deformations and multiphase flow coupling.This study introduces an SPH-based multiphase flow model for simulating bubbly flows involving various immiscible fluids.The model uses a volume-based density discretization equation,ensuring numerical accuracy near interfaces,independent of density ratios,thus maintaining accuracy and stability even at high density ratios.By integrating the Continuous Surface Force method for surface tension into the multiphase SPH model,the study simulates interfacial behavior between phases.The model accurately predicts Laplace pressure differences across interfaces,maintaining interface stability at density ratios up to 100.0.Simulations of single and double bubble ascents elucidate the influence of the Bond number on bubble shape,rising distance,and velocity.As the Bond number increases,bubbles flatten and develop tails,affecting their integrity.The study also simulates multiple bubbles ascents in water,showcasing the model's ability to capture complex interfacial behaviors in bubbly flows,including deformation,adsorption,coalescence,and tearing.展开更多
In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as...In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.展开更多
A mechanical model of liquid crystals （LCs） is applied to study the polymorphism of homologous series of terphenyl compounds. With a senti-experimental molecular orbit method, we calculate the moment of inertia whic...A mechanical model of liquid crystals （LCs） is applied to study the polymorphism of homologous series of terphenyl compounds. With a senti-experimental molecular orbit method, we calculate the moment of inertia which represents the rotation state to describe the phase transition temperature obtained from experimental data. We propose a novel explanation of the phase sequence or polymorphism of LC materials using the two key parameters, the moment of inertia and critical rotational velocity. The effect of molecular polarity on the appearance of liquid crystalline is also discussed.展开更多
In this paper, we present an improved multi-order parameter model for multi-component model of polycrystalline solidification. We introduce an interpolation function in the phase field dynamical equation to obtain con...In this paper, we present an improved multi-order parameter model for multi-component model of polycrystalline solidification. We introduce an interpolation function in the phase field dynamical equation to obtain controllable grain boundary energy at large undercooling. The same interpolation function is also employed in the kinetics coefficient to allow for better control of grain boundary migration. Temperature dependent phase field parameters and noise terms are consistently coupled into the dynamics of a binary system in a manner that allows for quantitative simulations in the thin interface limit. The model is applied to multi-phase solidification in Al-Cu alloy, where a parabolic fitting method is employed to model the free energy of Al-Cu phases and two-phase nucleation is demonstrated in directional solidification.展开更多
The perturbation in the magnetic field generated by the rotation or oscillation of magnetic domains in magnetic materials can emit low-frequency electromagnetic waves,which are expected to be used in low-frequency com...The perturbation in the magnetic field generated by the rotation or oscillation of magnetic domains in magnetic materials can emit low-frequency electromagnetic waves,which are expected to be used in low-frequency communications.However,the magnetic emission intensity,defined by the perturbation ability,of current commercially applied amorphous alloys,such as Metglas,cannot meet the application requirements for low-frequency antennas due to the domain motion energy loss.Herein,a multi-phase Metglas amorphous alloy was constructed by incorporatingα-Fe nanocrystals using rapid annealing to manipulate the domain movement.It was found that 3.89 times higher magnetic emission intensity is obtained compared to the pristine due to the synergism of the deformation and displacement mechanisms.Moreover,the low-frequency magnetic emission performance verification was carried out by preparing magnetoelectric composites as the antenna vibrator by assembling the alloy and macro piezoelectric fiber composites(MFC).Enhancements of magnetic emission intensity are found at 93.3%and 49.2%at the first and second harmonic frequencies compared with the unmodified alloy vibrator.Therefore,the approach leads to the development of high-performance communication with a novel standard for evaluation.展开更多
Wax deposition during crude oil production,transportation,and processing has been a headache since the early days of oil utilization.It may lead to low mobility ratios,blockage of production tubing/pipelines as well a...Wax deposition during crude oil production,transportation,and processing has been a headache since the early days of oil utilization.It may lead to low mobility ratios,blockage of production tubing/pipelines as well as fouling of surface and processing facilities,among others.These snags cause massive financial constraints increasing projects’turnover.Decades of meticulous research have been dedicated to this problem that is worth a review.Thus,this paper reviews the mechanisms,experimentation,thermodynamic and kinetic modeling,prediction,and remediation techniques of wax deposition.An overall assessment suggests that available models are more accurate for single than multi-phase flows while the kind of remediation and deployment depend on the environment and severity level.In severe cases,both chemical and mechanical are synergistically deployed.Moreover,future prospective research areas that require attention are proposed.Generally,this review could be a valuable tool for novice researchers as well as a foundation for further research on this topic.展开更多
Aiming at the simulation of multi-phase flow in the wellbore during the processes of gas kick and well killing of complex-structure wells(e.g.,directional wells,extended reach wells,etc.),a database including 3561 gro...Aiming at the simulation of multi-phase flow in the wellbore during the processes of gas kick and well killing of complex-structure wells(e.g.,directional wells,extended reach wells,etc.),a database including 3561 groups of experimental data from 32 different data sources is established.Considering the effects of fluid viscosity,pipe size,interfacial tension,fluid density,pipe inclination and other factors on multi-phase flow parameters,a new gas-liquid two-phase drift flow relation suitable for the full flow pattern and full dip range is established.The distribution coefficient and gas drift velocity models with a pipe inclination range of-90°–90°are established by means of theoretical analysis and data-driven.Compared with three existing models,the proposed models have the highest prediction accuracy and most stable performance.Using a well killing case with the backpressure method in the field,the applicability of the proposed model under the flow conditions with a pipe inclination range of-90°–80°is verified.The errors of the calculated shut in casing pressure,initial back casing pressure and casing pressure when adjusting the displacement are 2.58%,3.43%,5.35%,respectively.The calculated results of the model are in good agreement with the field backpressure data.展开更多
The Cobourg limestone is a very low porosity rock consisting of lighter nodular regions that are predominantly calcite and darker regions consisting of calcite,quartz,dolomite,and an appreciable clay fraction.This pap...The Cobourg limestone is a very low porosity rock consisting of lighter nodular regions that are predominantly calcite and darker regions consisting of calcite,quartz,dolomite,and an appreciable clay fraction.This paper presents the application of the theory of multi-phasic composites to estimate the possible maximum effective thermal conductivity of the heterogeneous rocks.The thermal conductivity estimates are expected to be representative of the intact rock,without fractures or fissures that can influence the heat conduction process.The estimates are therefore indicative of the thermal properties of the rock in undisturbed regions unaffected by the influences of stress relief and excavation damage during construction of deep ground repositories for the disposal of heat-emitting nuclear waste.展开更多
Interfacial energy anisotropy plays an important role in tilted growth of eutectics. However, previous studies mainly focused on the solid-solid interface energy anisotropy, and whether the solid-liquid interface ener...Interfacial energy anisotropy plays an important role in tilted growth of eutectics. However, previous studies mainly focused on the solid-solid interface energy anisotropy, and whether the solid-liquid interface energy anisotropy can significantly affect the tilted growth of eutectics still remains unclear. In this study, a multi-phase field model is employed to investigate both the effect of solid-liquid interfacial energy anisotropy and the effect of solid-solid interfacial energy anisotropy on tilted growth of eutectics. The findings reveal that both the solid-liquid interfacial energy anisotropy and the solid-solid interfacial energy anisotropy can induce the tilted growth of eutectics. The results also demonstrate that when the rotation angle is within a range of 30°-60°, the growth of tilted eutectics is governed jointly by the solid-solid interfacial energy anisotropy and the solid-liquid interfacial energy anisotropy;otherwise, it is mainly controlled by the solid-solid interfacial energy anisotropy. Further analysis shows that the unequal pinning angle at triple point caused by the adjustment of the force balance results in different solute-diffusion rates on both sides of triple point. This will further induce an asymmetrical concentration distribution along the pulling direction near the solid-liquid interface and the tilted growth of eutectics. Our findings not only shed light on the formation mechanism of tilted eutectics but also provide theoretical guidance for controlling the microstructure evolution during eutectic solidification.展开更多
We present a multi-phase image segmentation method based on the histogram of the Gabor feature space,which consists of a set of Gabor-filter responses with various orientations,scales and frequencies.Our model replace...We present a multi-phase image segmentation method based on the histogram of the Gabor feature space,which consists of a set of Gabor-filter responses with various orientations,scales and frequencies.Our model replaces the error function term in the original fuzzy region competition model with squared 2-Wasserstein distance function,which is a metric to measure the distance of two histograms.The energy functional is minimized by alternative minimization method and the existence of closed-form solutions is guaranteed when the exponent of the fuzzy membership term being 1 or 2.We test our model on both simple synthetic texture images and complex natural images with two or more phases.Experimental results are shown and compared to other recent results.展开更多
The phase-field method has emerged as the method of choice for the description of microstructure evolution and phase transitions in metallic materials.Following general thermodynamic laws a set of evolution equations ...The phase-field method has emerged as the method of choice for the description of microstructure evolution and phase transitions in metallic materials.Following general thermodynamic laws a set of evolution equations for the structural variables of the system,the so called phase-fields,are derived.The paper reviews shortly the theoretical background of the multi-phase-field.Different examples demonstrating the applicability of the method to technical steels will be presented ranging from deformation of the dendritic strand shell during peritectic transformation,grain growth in Austenite to stress driven growth of Pearlite.展开更多
In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the phys...In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the physical domain to resolvethe relevant scales in multiscale physical systems while minimizing computationalcosts. The algorithm is a generalization of the moving mesh methods basedon harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To make 3D moving mesh simulations possible,the key is to develop an efficient mesh redistribution procedure so that this part willcost as little as possible comparing with the solution evolution part. Since the meshredistribution procedure normally requires to solve large size matrix equations, wewill describe a procedure to decouple the matrix equation to a much simpler blocktridiagonaltype which can be efficiently solved by a particularly designed multi-gridmethod. To demonstrate the performance of the proposed 3D moving mesh strategy,the algorithm is implemented in finite element simulations of fluid-fluid interface interactionsin multiphase flows. To demonstrate the main ideas, we consider the formationof drops by using an energetic variational phase field model which describesthe motion of mixtures of two incompressible fluids. Numerical results on two- andthree-dimensional simulations will be presented.展开更多
基金supported by a grant from the Research Grants Council of the Hong Kong Special Administrative Region,China(Project No.15308024)a grant from Research Centre for Carbon-Strategic Catalysis,The Hong Kong Polytechnic University(CE2X).
文摘Water electrolyzers play a crucial role in green hydrogen production.However,their efficiency and scalability are often compromised by bubble dynamics across various scales,from nanoscale to macroscale components.This review explores multi-scale modeling as a tool to visualize multi-phase flow and improve mass transport in water electrolyzers.At the nanoscale,molecular dynamics(MD)simulations reveal how electrode surface features and wettability influence nanobubble nucleation and stability.Moving to the mesoscale,models such as volume of fluid(VOF)and lattice Boltzmann method(LBM)shed light on bubble transport in porous transport layers(PTLs).These insights inform innovative designs,including gradient porosity and hydrophilic-hydrophobic patterning,aimed at minimizing gas saturation.At the macroscale,VOF simulations elucidate two-phase flow regimes within channels,showing how flow field geometry and wettability affect bubble discharging.Moreover,artificial intelligence(AI)-driven surrogate models expedite the optimization process,allowing for rapid exploration of structural parameters in channel-rib flow fields and porous flow field designs.By integrating these approaches,we can bridge theoretical insights with experimental validation,ultimately enhancing water electrolyzer performance,reducing costs,and advancing affordable,high-efficiency hydrogen production.
基金supported by the National Natural Science Foundation of China,(Grant No.42072144)Shengli Oilfield,SINOPEC,China(Nos.30200018-21-ZC0613-0030 and 30200018-20-ZC0613-0116)。
文摘Based on the theory of superimposed deformation and the regional tectonic background,the multi-phase non-coaxial superimposed structures in Junggar Basin were systematically analyzed using seismic interpretation,field outcrop observation,and paleo-stress field recovery methods according to the characteristics of the current tectonic framework.Moreover,the tectonic evolution process of the basin was reconstructed using sandbox analogue modelling technology.The results showed that the study area has experienced five phases of non-coaxial deformation with superimposition:The first phase of deformation(D1)is characterized by NNE-SSW extension during late Carboniferous to early Permian,which formed large graben,half graben and other extensional structure style around the basin.The second phase of deformation(D2)is represented by NE-SW compression during the middle to late Permian,and it comprised numerous contraction structures that developed based on D1.The basic form of the entire basin is alternating uplift and depression.The third phase of deformation(D3)is the NW-SE transpressional strike-slip in the Triassic-Jurassic,which produced numerous strike-slip structural styles in the middle part of the basin.The fourth phase of deformation(D4)is the uniform sedimentation during Cretaceous,and the fifth phase(D5)is the compression along NNE-SSW due to the North Tianshan northward thrust,which produced three rows of fold thrust belts and tear faults in the front of the mountain in the southern margin of the basin.The newly established three-dimensional tectonic evolution model shows that,based on the large number of NW-trending grabens and half grabens in the Carboniferous basement of Junggar Basin,multiple level NE trending uplifts have formed with the joint superposition of the late structural inversion and multiple stress fields.This has resulted in the current tectonic units of alternating uplifts and depressions in different directions in the study area.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12205023,U2230401,12374056,U23A20537,11904027)。
文摘A thermodynamically complete multi-phase equation of state(EOS)applicable to both dense and porous metals at wide ranges of temperature and pressure is constructed.A standard three-term decomposition of the Helmholtz free energy as a function of specific volume and temperature is presented,where the cold component models both compression and expansion states,the thermal ion component introduces the Debye approximation and melting entropy,and the thermal electron component employs the Thomas-Fermi-Kirzhnits(TFK)model.The porosity of materials is considered by introducing the dynamic porosity coefficientαand the constitutive P-αrelation,connecting the thermodynamic properties between dense and porous systems,allowing for an accurate description of the volume decrease caused by void collapse while maintaining the quasi-static thermodynamic properties of porous systems identical to the dense ones.These models enable the EOS applicable and robust at wide ranges of temperature,pressure and porosity.A systematic evaluation of the new EOS is conducted with aluminum(Al)as an example.300 K isotherm,shock Hugoniot,as well as melting curves of both dense and porous Al are calculated,which shows great agreements with experimental data and validates the effectiveness of the models and the accuracy of parameterizations.Notably,it is for the first time Hugoniot P-σcurves up to 10~6 GPa and shock melting behaviors of porous Al are derived from analytical EOS models,which predict much lower compression limit and shock melting temperatures than those of dense Al.
基金the Slovak Research and Development Agency under the contract No.COST-0022-06,APVV-51-061505the 6th FP EU NESPA+5 种基金the Slovak Grant Agency VEGA (2/7197/27,2/7194/27,2/7195/27)NANOSMART,Centre of Excellence (1/1/2007-31/12/2010)Slovak Academy of Sciences,by KMM-NoE 502243-2 (10/2004-9/2008)NENAMAT INCO-CT-2003-510363COST Action 536 and COST Action 538János Bolyai Research Grant NSF-MTA-OTKA grant-MTA:96/OTKA:049953,OTKA 63609
文摘Based on the fundamental equations of the mechanics of solid continuum, the paper employs an analytical model for determination of elastic thermal stresses in isotropic continuum represented by periodically distributed spherical particles with different distributions in an infinite matrix, imaginarily divided into identical cells with dimensions equal to inter-particle distances, containing a central spherical particle with or without a spherical envelope on the particle surface. Consequently, the multi-particle-（envelope）- matrix system, as a model system regarding the analytical modelling, is applicable to four types of multi-phase materials. As functions of the particle volume fraction v, the inter-particle distances dl, d2, d3 along three mutually per- pendicular axes, and the particle and envelope radii, R1 and R2, respectively, the thermal stresses within the cell, are originated during a cooling process as a consequence of the difference in thermal expansion coefficients of phases rep- resented by the matrix, envelope and particle. Analytical-（experimental）-computational lifetime prediction methods for multi-phase materials are proposed, which can be used in engineering with appropriate values of parameters of real multi-phase materials.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.52161002,51661020,and 11364024)the Postdoctoral Science Foundation of China(Grant No.2014M560371)the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China(Grant No.J201304).
文摘A graphics-processing-unit(GPU)-parallel-based computational scheme is developed to realize the competitive growth process of converging bi-crystal in two-dimensional states in the presence of forced convection conditions by coupling a multi-phase field model and a lattice Boltzmann model.The elimination mechanism in the evolution process is analyzed for the three conformational schemes constituting converging bi-crystals under pure diffusion and forced convection conditions,respectively,expanding the research of the competitive growth of columnar dendrites under melt convection conditions.The results show that the elimination mechanism for the competitive growth of converging bi-crystals of all three configurations under pure diffusion conditions follows the conventional Walton-Chalmers model.When there is forced convection with lateral flow in the liquid phase,the anomalous elimination phenomenon of unfavorable dendrites eliminating favorable dendrites occurs in the grain boundaries.In particular,the anomalous elimination phenomenon is relatively strong in conformation 1 and conformation 2 when the orientation angle of unfavorable dendrites is small,and relatively weak in conformation 3.Moreover,the presence of convection increases the tip growth rate of both favorable and unfavorable dendrites in the grain boundary.In addition,the parallelization of the multi-phase-field-lattice Boltzmann model is achieved by designing the parallel computation of the model on the GPU platform concerning the computerunified-device-architecture parallel technique,and the results show that the parallel computation of this model based on the GPU has absolute advantages,and the parallel acceleration is more obvious as the computation area increases.
基金Project supported by the National Natural Science Foundation of China (No. 10502025)Fok Ying Tong Education Foundation (No.101005)University Foundation of Jiangsu Province (No.05KJB1300421)
文摘To determine the time-independent constitutive modeling for porous and multi- phase nanocrystalline materials and understand the effects of grain size and porosity on their mechanical behavior, each phase was treated as a mixture of grain interior and grain bound- ary, and pores were taken as a single phase, then Budiansky＇s self-consistent method was used to calculate the Young＇s modulus of porous, possible multi-phase, nanocrystalline materials, the prediction being in good agreement with the results in the literature. Further, the established method is extended to simulate the constitutive relations of porous and possible multi-phase nanocrystalline materials with small plastic deformation in conjunction with the secant-moduli approach and iso-strain assumption. Comparisons between the experimental grain size and porosity dependent mechanical data and the corresponding predictions using the established model show that it appears to be capable of describing the time-independent mechanical behaviors for porous and multi-phase nanocrystalline materials in a small plastic strain range. Further discussion on the modification factor, the advantages and limitations of the model developed were present.
基金supported by the National Natural Science Foundation of China(Grant Nos.52161002,51661020,and 11504149)the Postdoctoral Science Foundation of China(Grant No.2014M560371)the Funds for Distinguished Young Scientists of Lanzhou University of Technology,China(Grant No.J201304)。
文摘The multi-phase field model of grain competitive growth during directional solidification of alloy is established.Solving multi-phase field models for thin interface layer thickness conditions,the grain boundary evolution and grain elimination during the competitive growth of SCN-0.24-wt%camphor model alloy bi-crystals are investigated.The effects of different crystal orientations and pulling velocities on grain boundary microstructure evolution are quantitatively analyzed.The obtained results are shown below.In the competitive growth of convergent bi-crystals,when favorably oriented dendrites are in the same direction as the heat flow and the pulling speed is too large,the orientation angle of the bi-crystal from small to large size is the normal elimination phenomenon of the favorably oriented dendrite,blocking the unfavorably oriented dendrite,and the grain boundary is along the growth direction of the favorably oriented dendrite.When the pulling speed becomes small,the grain boundary shows the anomalous elimination phenomenon of the unfavorably oriented dendrite,eliminating the favorably oriented dendrite.In the process of competitive growth of divergent bi-crystal,when the growth direction of favorably oriented dendrites is the same as the heat flow direction and the orientation angle of unfavorably oriented grains is small,the frequency of new spindles of favorably oriented grains is significantly higher than that of unfavorably oriented grains,and as the orientation angle of unfavorably oriented dendrites becomes larger,the unfavorably oriented grains are more likely to have stable secondary dendritic arms,which in turn develop new primary dendritic arms to occupy the liquid phase grain boundary space,but the grain boundary direction is still parallel to favorably oriented dendrites.In addition,the tertiary dendritic arms on the developed secondary dendritic arms may also be blocked by the surrounding lateral branches from further developing into nascent main axes,this blocking of the tertiary dendritic arms has a random nature,which can have aninfluence on the generation of nascent primary main axes in the grain boundaries.
基金the State Key Laboratory of Intelligent Optimized Manufacturing in Mining&Metallurgy Process Open Research Fund(No.JTKY202404622)and(No.BGRIMM-KZSKL-2023-12).
文摘Smoothed particle hydrodynamics(SPH)is a Lagrangian,mesh-free numerical method renowned for its ability to handle fluid dynamics problems with large interface deformations and multiphase flow coupling.This study introduces an SPH-based multiphase flow model for simulating bubbly flows involving various immiscible fluids.The model uses a volume-based density discretization equation,ensuring numerical accuracy near interfaces,independent of density ratios,thus maintaining accuracy and stability even at high density ratios.By integrating the Continuous Surface Force method for surface tension into the multiphase SPH model,the study simulates interfacial behavior between phases.The model accurately predicts Laplace pressure differences across interfaces,maintaining interface stability at density ratios up to 100.0.Simulations of single and double bubble ascents elucidate the influence of the Bond number on bubble shape,rising distance,and velocity.As the Bond number increases,bubbles flatten and develop tails,affecting their integrity.The study also simulates multiple bubbles ascents in water,showcasing the model's ability to capture complex interfacial behaviors in bubbly flows,including deformation,adsorption,coalescence,and tearing.
文摘In this paper, a finite element method （FEM）-based multi-phase problem based on a newly proposed thermal elastoplastic constitutive model for saturated/unsaturated geomaterial is discussed. A program of FEM named as SOFT, adopting unified field equations for thermo-hydro-mechanical-air （THMA） behavior of geomaterial and using finite element-finite difference （FE-FD） scheme for so/l-water-air three-phase coupling problem, is used in the numerical simulation. As an application of the newly proposed numerical method, two engineering problems, one for slope failure in unsaturated model ground and another for in situ heating test related to deep geological repository of high-level radioactive waste （HLRW）, are simulated. The model tests on slope failure in unsaturated Shirasu ground, carried out by Kitamura et al. （2007）, is simulated in the framework of soil-water-air three-phase coupling under the condition of constant temperature. While the in situ heating test reported by Munoz （2006） is simulated in the same framework under the conditions of variable temperature hut constant air pressure.
文摘A mechanical model of liquid crystals （LCs） is applied to study the polymorphism of homologous series of terphenyl compounds. With a senti-experimental molecular orbit method, we calculate the moment of inertia which represents the rotation state to describe the phase transition temperature obtained from experimental data. We propose a novel explanation of the phase sequence or polymorphism of LC materials using the two key parameters, the moment of inertia and critical rotational velocity. The effect of molecular polarity on the appearance of liquid crystalline is also discussed.
基金financial support by the National Foundation of China under Grant No. 51875131Nikolas Provatas acknowledges financial support by the Canada Research Chairs (CRC) Programsupported by the China Scholarship Council as a Graduate Trainee at McGill University。
文摘In this paper, we present an improved multi-order parameter model for multi-component model of polycrystalline solidification. We introduce an interpolation function in the phase field dynamical equation to obtain controllable grain boundary energy at large undercooling. The same interpolation function is also employed in the kinetics coefficient to allow for better control of grain boundary migration. Temperature dependent phase field parameters and noise terms are consistently coupled into the dynamics of a binary system in a manner that allows for quantitative simulations in the thin interface limit. The model is applied to multi-phase solidification in Al-Cu alloy, where a parabolic fitting method is employed to model the free energy of Al-Cu phases and two-phase nucleation is demonstrated in directional solidification.
基金supported by the Key Research and Development Program of Hubei Province(No.2021BAA214)the Open Fund of Sanya Science and Education Innovation Park of Wuhan University of Technology(Nos.2021KF0022,2021KF0013,and 2020KF0026)+2 种基金Independent Innovation Projects of the Hubei Longzhong Laboratory(Nos.2022ZZ-34 and 2022ZZ-35)the National Science Fund for Distinguished Young Scholars of Hubei Province(No.201CFA067)the National innovation and entrepreneurship training program for college students(Nos.202310497010 and S202310497026).
文摘The perturbation in the magnetic field generated by the rotation or oscillation of magnetic domains in magnetic materials can emit low-frequency electromagnetic waves,which are expected to be used in low-frequency communications.However,the magnetic emission intensity,defined by the perturbation ability,of current commercially applied amorphous alloys,such as Metglas,cannot meet the application requirements for low-frequency antennas due to the domain motion energy loss.Herein,a multi-phase Metglas amorphous alloy was constructed by incorporatingα-Fe nanocrystals using rapid annealing to manipulate the domain movement.It was found that 3.89 times higher magnetic emission intensity is obtained compared to the pristine due to the synergism of the deformation and displacement mechanisms.Moreover,the low-frequency magnetic emission performance verification was carried out by preparing magnetoelectric composites as the antenna vibrator by assembling the alloy and macro piezoelectric fiber composites(MFC).Enhancements of magnetic emission intensity are found at 93.3%and 49.2%at the first and second harmonic frequencies compared with the unmodified alloy vibrator.Therefore,the approach leads to the development of high-performance communication with a novel standard for evaluation.
基金contributions from colleagues and support from Sinopec Company limited(Project P19018-2)the National Natural Science Foundation of China(52174047)。
文摘Wax deposition during crude oil production,transportation,and processing has been a headache since the early days of oil utilization.It may lead to low mobility ratios,blockage of production tubing/pipelines as well as fouling of surface and processing facilities,among others.These snags cause massive financial constraints increasing projects’turnover.Decades of meticulous research have been dedicated to this problem that is worth a review.Thus,this paper reviews the mechanisms,experimentation,thermodynamic and kinetic modeling,prediction,and remediation techniques of wax deposition.An overall assessment suggests that available models are more accurate for single than multi-phase flows while the kind of remediation and deployment depend on the environment and severity level.In severe cases,both chemical and mechanical are synergistically deployed.Moreover,future prospective research areas that require attention are proposed.Generally,this review could be a valuable tool for novice researchers as well as a foundation for further research on this topic.
基金Supported by the Project of National Natural Science Foundation of China(51991363,51974350)Young Changjiang Scholars Award Program(Q2016135)Ministry of Education Innovation Team Project(IRT_14R58)。
文摘Aiming at the simulation of multi-phase flow in the wellbore during the processes of gas kick and well killing of complex-structure wells(e.g.,directional wells,extended reach wells,etc.),a database including 3561 groups of experimental data from 32 different data sources is established.Considering the effects of fluid viscosity,pipe size,interfacial tension,fluid density,pipe inclination and other factors on multi-phase flow parameters,a new gas-liquid two-phase drift flow relation suitable for the full flow pattern and full dip range is established.The distribution coefficient and gas drift velocity models with a pipe inclination range of-90°–90°are established by means of theoretical analysis and data-driven.Compared with three existing models,the proposed models have the highest prediction accuracy and most stable performance.Using a well killing case with the backpressure method in the field,the applicability of the proposed model under the flow conditions with a pipe inclination range of-90°–80°is verified.The errors of the calculated shut in casing pressure,initial back casing pressure and casing pressure when adjusting the displacement are 2.58%,3.43%,5.35%,respectively.The calculated results of the model are in good agreement with the field backpressure data.
文摘The Cobourg limestone is a very low porosity rock consisting of lighter nodular regions that are predominantly calcite and darker regions consisting of calcite,quartz,dolomite,and an appreciable clay fraction.This paper presents the application of the theory of multi-phasic composites to estimate the possible maximum effective thermal conductivity of the heterogeneous rocks.The thermal conductivity estimates are expected to be representative of the intact rock,without fractures or fissures that can influence the heat conduction process.The estimates are therefore indicative of the thermal properties of the rock in undisturbed regions unaffected by the influences of stress relief and excavation damage during construction of deep ground repositories for the disposal of heat-emitting nuclear waste.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 51871183 and 51571165)。
文摘Interfacial energy anisotropy plays an important role in tilted growth of eutectics. However, previous studies mainly focused on the solid-solid interface energy anisotropy, and whether the solid-liquid interface energy anisotropy can significantly affect the tilted growth of eutectics still remains unclear. In this study, a multi-phase field model is employed to investigate both the effect of solid-liquid interfacial energy anisotropy and the effect of solid-solid interfacial energy anisotropy on tilted growth of eutectics. The findings reveal that both the solid-liquid interfacial energy anisotropy and the solid-solid interfacial energy anisotropy can induce the tilted growth of eutectics. The results also demonstrate that when the rotation angle is within a range of 30°-60°, the growth of tilted eutectics is governed jointly by the solid-solid interfacial energy anisotropy and the solid-liquid interfacial energy anisotropy;otherwise, it is mainly controlled by the solid-solid interfacial energy anisotropy. Further analysis shows that the unequal pinning angle at triple point caused by the adjustment of the force balance results in different solute-diffusion rates on both sides of triple point. This will further induce an asymmetrical concentration distribution along the pulling direction near the solid-liquid interface and the tilted growth of eutectics. Our findings not only shed light on the formation mechanism of tilted eutectics but also provide theoretical guidance for controlling the microstructure evolution during eutectic solidification.
文摘We present a multi-phase image segmentation method based on the histogram of the Gabor feature space,which consists of a set of Gabor-filter responses with various orientations,scales and frequencies.Our model replaces the error function term in the original fuzzy region competition model with squared 2-Wasserstein distance function,which is a metric to measure the distance of two histograms.The energy functional is minimized by alternative minimization method and the existence of closed-form solutions is guaranteed when the exponent of the fuzzy membership term being 1 or 2.We test our model on both simple synthetic texture images and complex natural images with two or more phases.Experimental results are shown and compared to other recent results.
文摘The phase-field method has emerged as the method of choice for the description of microstructure evolution and phase transitions in metallic materials.Following general thermodynamic laws a set of evolution equations for the structural variables of the system,the so called phase-fields,are derived.The paper reviews shortly the theoretical background of the multi-phase-field.Different examples demonstrating the applicability of the method to technical steels will be presented ranging from deformation of the dendritic strand shell during peritectic transformation,grain growth in Austenite to stress driven growth of Pearlite.
基金the Joint Applied Mathematics Research Institute of Peking University and Hong Kong Baptist University.Li was also partially supported by the National Basic Research Program of China under the grant 2005CB321701The research of Tang was supported by CERG Grants of Hong Kong Research Grant Council,FRG grants of Hong Kong Baptist University,and NSAF Grant#10476032 of National Science Foundation of China.He was supported in part by the Chinese Academy of Sciences while visiting its Institute of Computational Mathematics.
文摘In this paper, we present an adaptive moving mesh algorithm for meshesof unstructured polyhedra in three space dimensions. The algorithm automaticallyadjusts the size of the elements with time and position in the physical domain to resolvethe relevant scales in multiscale physical systems while minimizing computationalcosts. The algorithm is a generalization of the moving mesh methods basedon harmonic mappings developed by Li et al. [J. Comput. Phys., 170 (2001), pp. 562-588, and 177 (2002), pp. 365-393]. To make 3D moving mesh simulations possible,the key is to develop an efficient mesh redistribution procedure so that this part willcost as little as possible comparing with the solution evolution part. Since the meshredistribution procedure normally requires to solve large size matrix equations, wewill describe a procedure to decouple the matrix equation to a much simpler blocktridiagonaltype which can be efficiently solved by a particularly designed multi-gridmethod. To demonstrate the performance of the proposed 3D moving mesh strategy,the algorithm is implemented in finite element simulations of fluid-fluid interface interactionsin multiphase flows. To demonstrate the main ideas, we consider the formationof drops by using an energetic variational phase field model which describesthe motion of mixtures of two incompressible fluids. Numerical results on two- andthree-dimensional simulations will be presented.
