Simulation study of supercritical carbon dioxide jet fracturing for carbonate geothermal reservoir based on fluid-thermo-mechanical coupling model

Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.

Study on inhomogeneous cooling behavior of extruded profile with unequal and large thicknesses during quenching using thermo-mechanical coupling model

The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanical coupling model for simulating the on-line water quenching of extruded profile with unequal and large thicknesses was developed. The temperature field, residual stress field and distortion of profile during quenching were investigated systematically. The results show that heat transfer coefficient increases as water flow rate increases. The peak heat transfer coefficient with higher water flow rates appears at lower interface temperatures. The temperature distribution across the cross-section of profile during quenching is severe nonuniform and the maximum temperature difference is 300 ℃ at quenching time of 3.49 s. The temperature difference through the thickness of different parts of profile first increases sharply to a maximum value, and then gradually decreases. The temperature gradient increases obviously with the increase of thickness of parts. After quenching, there exist large residual stresses on the inner side of joints of profile and the two ends of part with thickness of 10 mm. The profile presents a twisting-type distortion across the cross-section under non-uniform cooling and the maximum twisting angle during quenching is 2.78°.

A thermo-mechanical damage constitutive model for deep rock considering brittleness-ductility transition characteristics

This paper developed a statistical damage constitutive model for deep rock by considering the effects of external load and thermal treatment temperature based on the distortion energy.The model parameters were determined through the extremum features of stress−strain curve.Subsequently,the model predictions were compared with experimental results of marble samples.It is found that when the treatment temperature rises,the coupling damage evolution curve shows an S-shape and the slope of ascending branch gradually decreases during the coupling damage evolution process.At a constant temperature,confining pressure can suppress the expansion of micro-fractures.As the confining pressure increases the rock exhibits ductility characteristics,and the shape of coupling damage curve changes from an S-shape into a quasi-parabolic shape.This model can well characterize the influence of high temperature on the mechanical properties of deep rock and its brittleness-ductility transition characteristics under confining pressure.Also,it is suitable for sandstone and granite,especially in predicting the pre-peak stage and peak stress of stress−strain curve under the coupling action of confining pressure and high temperature.The relevant results can provide a reference for further research on the constitutive relationship of rock-like materials and their engineering applications.

Mechanical behavior and damage constitutive model of sandstone under hydro-mechanical (H-M) coupling

Underground engineering often passes through water-rich fractured rock masses, which are prone to fracture and instability under the long-term coupling of in-situ stress field and pore water(P-W) pressure, ultimately threatening the stability of underground structures. In order to explore the mechanical properties of rocks under H-M coupling, the corresponding damage constitutive(D-C) model has become the focus of attention. Considering the inadequacy of the current research on rock strength parameters,energy evolution characteristics and D-C model under H-M coupling, the mechanical properties of typical sandstone samples are discussed based on laboratory tests. The results show that the variation of characteristic stresses of sandstone under H-M coupling conforms to the normalized attenuation equation and Mohr-Coulomb(M-C) criterion. The P-W pressure mechanism of sandstone exhibits a dynamic change from softening effect to H-M fracturing effect. The closure stress is mainly provided by cohesive strength, while the initiation stress, damage stress, and peak stress are jointly dominated by cohesive strength and friction strength. In addition, residual stress is attributed to the friction strength formed by the bite of the fracture surface. Subsequently, the energy evolution characteristics of sandstone under H-M coupling were studied, and it was found that P-W pressure weakened the energy storage capacity and energy dissipation capacity of sandstone, and H-M fracturing was an important factor in reducing its energy storage efficiency. Finally, combined with energy dissipation theory and statistical damage theory, two types of D-C models considering P-W pressure are proposed accordingly, and the model parameters can be determined by four methods. The application results indicate that the proposed and modified D-C models have high reliability, and can characterize the mechanical behavior of sandstone under H-M coupling, overcome the inconvenience of existing D-C models due to excessive mechanical parameters,and can be applied to the full-range stress–strain process. The results are conducive to revealing the deformation and damage mechanisms of rocks under H-M coupling, and can provide theoretical guidance for related engineering problems.

Degree of Freedom Analysis for Holographic MIMO Based on a Mutual-Coupling-Compliant Channel Model

Degree of freedom(DOF)is a key indicator for spatial multiplexing layers of a wireless channel.Traditionally,the channel of a multiple-input multiple-output(MIMO)half-wavelength dipole array has a DOF that equals the antenna number.However,recent studies suggest that the DOF could be less than the antenna number when strong mutual coupling is considered.We utilize a mutual-coupling-compliant channel model to investigate the DOF of the holographic MIMO(HMIMO)channel and give a upper bound of the DOF with strong mutual coupling.Our numerical simulations demonstrate that a dense array can support more DOF per unit aperture as compared with a half-wavelength MIMO system.

Shear behavior of intact granite under thermo-mechanical coupling and three-dimensional morphology of shear-formed fractures

The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.

Coupled multiphysical model for investigation of influence factors in the application of microbially induced calcite precipitation

The study presents a comprehensive coupled thermo-bio-chemo-hydraulic(T-BCH)modeling framework for stabilizing soils using microbially induced calcite precipitation(MICP).The numerical model considers relevant multiphysics involved in MICP,such as bacterial ureolytic activities,biochemical reactions,multiphase and multicomponent transport,and alteration of the porosity and permeability.The model incorporates multiphysical coupling effects through well-established constitutive relations that connect parameters and variables from different physical fields.It was implemented in the open-source finite element code OpenGeoSys(OGS),and a semi-staggered solution strategy was designed to solve the couplings,allowing for flexible model settings.Therefore,the developed model can be easily adapted to simulate MICP applications in different scenarios.The numerical model was employed to analyze the effect of various factors,including temperature,injection strategies,and application scales.Besides,a TBCH modeling study was conducted on the laboratory-scale domain to analyze the effects of temperature on urease activity and precipitated calcium carbonate.To understand the scale dependency of MICP treatment,a large-scale heterogeneous domain was subjected to variable biochemical injection strategies.The simulations conducted at the field-scale guided the selection of an injection strategy to achieve the desired type and amount of precipitation.Additionally,the study emphasized the potential of numerical models as reliable tools for optimizing future developments in field-scale MICP treatment.The present study demonstrates the potential of this numerical framework for designing and optimizing the MICP applications in laboratory-,prototype-,and field-scale scenarios.

Optimization Control of Multi-Mode Coupling All-Wheel Drive System for Hybrid Vehicle

The all-wheel drive(AWD)hybrid system is a research focus on high-performance new energy vehicles that can meet the demands of dynamic performance and passing ability.Simultaneous optimization of the power and economy of hybrid vehicles becomes an issue.A unique multi-mode coupling(MMC)AWD hybrid system is presented to realize the distributed and centralized driving of the front and rear axles to achieve vectored distribution and full utilization of the system power between the axles of vehicles.Based on the parameters of the benchmarking model of a hybrid vehicle,the best model-predictive control-based energy management strategy is proposed.First,the drive system model was built after the analysis of the MMC-AWD’s drive modes.Next,three fundamental strategies were established to address power distribution adjustment and battery SOC maintenance when the SOC changed,which was followed by the design of a road driving force observer.Then,the energy consumption rate in the average time domain was processed before designing the minimum fuel consumption controller based on the equivalent fuel consumption coefficient.Finally,the advantage of the MMC-AWD was confirmed by comparison with the dynamic performance and economy of the BYD Song PLUS DMI-AWD.The findings indicate that,in comparison to the comparative hybrid system at road adhesion coefficients of 0.8 and 0.6,the MMC-AWD’s capacity to accelerate increases by 5.26%and 7.92%,respectively.When the road adhesion coefficient is 0.8,0.6,and 0.4,the maximum climbing ability increases by 14.22%,12.88%,and 4.55%,respectively.As a result,the dynamic performance is greatly enhanced,and the fuel savings rate per 100 km of mileage reaches 12.06%,which is also very economical.The proposed control strategies for the new hybrid AWD vehicle can optimize the power and economy simultaneously.

Time-dependent behavior and permeability evolution of limestone under hydro-mechanical coupling

Excessively high pore water pressure presents unpredictable risks to the safety of rock tunnels in mountainous regions that are predominantly composed of limestone. Investigating the creep characteristics and permeability evolution of limestone under varying hydrated conditions is crucial for a better understanding of the delayed deformation mechanisms of limestone rock tunnels. To this end, this paper initially conducts a series of multi-stage triaxial creep tests on limestone samples under varying pore water pressures. The experiment examines how pore water pressure affects limestone’s creep strain, strain rate, long-term strength, lifespan, and permeability, all within the context of hydraulicmechanical(HM) coupling. To better describe the creep behavior associated with pore water pressure, this paper proposes a new nonlinear fractional creep constitutive model. This constitutive model depicts the initial, steady-state, and accelerated phases of limestone’s creep behavior. Finally, the proposed model is applied to the numerical realization of deformation in limestone tunnel, validating the effectiveness of the proposed constitutive model in predicting tunnel’s creep deformation. This research enhances our understanding of limestone’s creep characteristics and permeability evolution under HM coupling, laying a foundation for assessing the longterm stability of mountain tunnels.

Contrasts of bimodal tropical instability waves(TIWs)-induced wind stress perturbations in the Pacific Ocean among observations,ocean models,and coupled climate models

The coupling between wind stress perturbations and sea surface temperature(SST)perturbations induced by tropical instability waves(TIWs)in the Pacific Ocean has been revealed previously and proven crucial to both the atmosphere and ocean.However,an overlooked fact by previous studies is that the loosely defined“TIWs”actually consist of two modes,including the Yanai wave-based TIW on the equator(hereafter eTIW)and the Rossby wave-based TIW off the equator(hereafter vTIW).Hence,the individual feedbacks of the wind stress to the bimodal TIWs remain unexplored.In this study,individual coupling relationships are established for both eTIW and v TIW,including the relationship between the TIW-induced SST perturbations and two components of wind stress perturbations,and the relationship between the TIW-induced wind stress perturbation divergence(curl)and the downwind(crosswind)TIW-induced SST gradients.Results show that,due to different distributions of eTIW and vTIW,the coupling strength induced by the eTIW is stronger on the equator,and that by the vTIW is stronger off the equator.The results of any of eTIW and vTIW are higher than those of the loosely defined TIWs.We further investigated how well the coupling relationships remained in several widely recognized oceanic general circulation models and fully coupled climate models.However,the coupling relationships cannot be well represented in most numerical models.Finally,we confirmed that higher resolution usually corresponds to more accurate simulation.Therefore,the coupling models established in this study are complementary to previous research and can be used to refine the oceanic and coupled climate models.

An internal ballistic model of electromagnetic railgun based on PFN coupled with multi-physical field and experimental validation

To accelerate the practicality of electromagnetic railguns,it is necessary to use a combination of threedimensional numerical simulation and experiments to study the mechanism of bore damage.In this paper,a three-dimensional numerical model of the augmented railgun with four parallel unconventional rails is introduced to simulate the internal ballistic process and realize the multi-physics field coupling calculation of the rail gun,and a test experiment of a medium-caliber electromagnetic launcher powered by pulse formation network(PFN)is carried out.Various test methods such as spectrometer,fiber grating and high-speed camera are used to test several parameters such as muzzle initial velocity,transient magnetic field strength and stress-strain of rail.Combining the simulation results and experimental data,the damage condition of the contact surface is analyzed.

Different El Niño Flavors and Associated Atmospheric Teleconnections as Simulated in a Hybrid Coupled Model

A previously developed hybrid coupled model(HCM)is composed of an intermediate tropical Pacific Ocean model and a global atmospheric general circulation model(AGCM),denoted as HCMAGCM.In this study,different El Niño flavors,namely the Eastern-Pacific(EP)and Central-Pacific(CP)types,and the associated global atmospheric teleconnections are examined in a 1000-yr control simulation of the HCMAGCM.The HCMAGCM indicates profoundly different characteristics among EP and CP El Niño events in terms of related oceanic and atmospheric variables in the tropical Pacific,including the amplitude and spatial patterns of sea surface temperature(SST),zonal wind stress,and precipitation anomalies.An SST budget analysis indicates that the thermocline feedback and zonal advective feedback dominantly contribute to the growth of EP and CP El Niño events,respectively.Corresponding to the shifts in the tropical rainfall and deep convection during EP and CP El Niño events,the model also reproduces the differences in the extratropical atmospheric responses during the boreal winter.In particular,the EP El Niño tends to be dominant in exciting a poleward wave train pattern to the Northern Hemisphere,while the CP El Niño tends to preferably produce a wave train similar to the Pacific North American(PNA)pattern.As a result,different climatic impacts exist in North American regions,with a warm-north and cold-south pattern during an EP El Niño and a warm-northeast and cold-southwest pattern during a CP El Niño,respectively.This modeling result highlights the importance of internal natural processes within the tropical Pacific as they relate to the genesis of ENSO diversity because the active ocean–atmosphere coupling is allowed only in the tropical Pacific within the framework of the HCMAGCM.

Multi-scenario Simulation and Spatial-temporal Analysis of LUCC in China's Coastal Zone Based on Coupled SD-FLUS Model

Increased human activities in China's coastal zone have resulted in the depletion of ecological land resources.Thus,conducting current and future multi-scenario simulation research on land use and land cover change(LUCC)is crucial for guiding the healthy and sustainable development of coastal zones.System dynamic(SD)-future land use simulation(FLUS)model,a coupled simulation model,was developed to analyze land use dynamics in China's coastal zone.This model encompasses five scenarios,namely,SSP1-RCP2.6(A),SSP2-RCP4.5(B),SSP3-RCP4.5(C),SSP4-RCP4.5(D),and SSP5-RCP8.5(E).The SD model simulates land use demand on an annual basis up to the year 2100.Subsequently,the FLUS model determines the spatial distribution of land use for the near term(2035),medium term(2050),and long term(2100).Results reveal a slowing trend in land use changes in China's coastal zone from 2000–2020.Among these changes,the expansion rate of construction land was the highest and exhibited an annual decrease.By 2100,land use predictions exhibit high accuracy,and notable differences are observed in trends across scenarios.In summary,the expansion of production,living,and ecological spaces toward the sea remains prominent.Scenario A emphasizes reduced land resource dependence,benefiting ecological land protection.Scenario B witnesses an intensified expansion of artificial wetlands.Scenario C sees substantial land needs for living and production,while Scenario D shows coastal forest and grassland shrinkage.Lastly,in Scenario E,the conflict between humans and land intensifies.This study presents pertinent recommendations for the future development,utilization,and management of coastal areas in China.The research contributes valuable scientific support for informed,long-term strategic decision making within coastal regions.

Hydrologic Response to Future Climate Change in the Dulong-Irra-waddy River Basin Based on Coupled Model Intercomparison Project 6

Within the context of the Belt and Road Initiative(BRI)and the China-Myanmar Economic Corridor(CMEC),the Dulong-Ir-rawaddy(Ayeyarwady)River,an international river among China,India and Myanmar,plays a significant role as both a valuable hydro-power resource and an essential ecological passageway.However,the water resources and security exhibit a high degree of vulnerabil-ity to climate change impacts.This research evaluates climate impacts on the hydrology of the Dulong-Irrawaddy River Basin(DIRB)by using a physical-based hydrologic model.We crafted future climate scenarios using the three latest global climate models(GCMs)from Coupled Model Intercomparison Project 6(CMIP6)under two shared socioeconomic pathways(SSP2-4.5 and SSP5-8.5)for the near(2025-2049),mid(2050-2074),and far future(2075-2099).The regional model using MIKE SHE based on historical hydrologic processes was developed to further project future streamflow,demonstrating reliable performance in streamflow simulations with a val-idation Nash-Sutcliffe Efficiency(NSE)of 0.72.Results showed that climate change projections showed increases in the annual precip-itation and potential evapotranspiration(PET),with precipitation increasing by 11.3%and 26.1%,and PET increasing by 3.2%and 4.9%,respectively,by the end of the century under SSP2-4.5 and SSP5-8.5.These changes are projected to result in increased annual streamflow at all stations,notably at the basin’s outlet(Pyay station)compared to the baseline period(with an increase of 16.1%and 37.0%at the end of the 21st century under SSP2-4.5 and SSP5-8.5,respectively).Seasonal analysis for Pyay station forecasts an in-crease in dry-season streamflow by 31.3%-48.9%and 22.5%-76.3%under SSP2-4.5 and SSP5-8.5,respectively,and an increase in wet-season streamflow by 5.8%-12.6%and 2.8%-33.3%,respectively.Moreover,the magnitude and frequency of flood events are pre-dicted to escalate,potentially impacting hydropower production and food security significantly.This research outlines the hydrological response to future climate change during the 21st century and offers a scientific basis for the water resource management strategies by decision-makers.

Numerical Modelling of Coupled Heat and Mass Transfer in Porous Materials: Application to Cinder Block Bricks

In this work, we present numerical modelling of coupled heat and mass transfer within porous materials. Our study focuses on cinder block bricks generally used in building construction. The material is assumed to be placed in air. Moisture content and temperature have been chosen as the main transfer drivers and the equations governing these transfer drivers are based on the Luikov model. These equations are solved by an implicit finite difference scheme. A Fortran code associated with the Thomas algorithm was used to solve the equations. The results show that heat and mass transfer depend on the temperature of the air in contact with the material. As this air temperature rises, the temperature within the material increases, and more rapidly at the material surface. Also, thermal conductivity plays a very important role in the thermal conduction of building materials and influences heat and mass transfer in these materials. Materials with higher thermal conductivity diffuse more heat.

Mathematical Modeling of Multiple Capacitor Coupled Substations (CCS) Impact on Transmission Lines and Approaches for Ferroresonance Suppression

Rural electrification remains a critical challenge in achieving equitable access to electricity, a cornerstone for poverty alleviation, economic growth, and improved living standards. Capacitor Coupled Substations (CCS) offer a promising solution for delivering cost-effective electricity to these underserved areas. However, the integration of multiple CCS units along a transmission network introduces complex interactions that can significantly impact voltage, current, and power flow. This study presents a detailed mathematical model to analyze the effects of varying distances and configurations of multiple CCS units on a transmission network, with a focus on voltage stability, power quality, and reactive power fluctuations. Furthermore, the research addresses the phenomenon of ferroresonance, a critical issue in networks with multiple CCS units, by developing and validating suppression strategies to ensure stable operation. Through simulation and practical testing, the study provides insights into optimizing CCS deployment, ultimately contributing to more reliable and efficient rural electrification solutions.

Thermo-mechanical coupled particle model for rock

A thermo-mechanical coupled particle model for simulation of thermally-induced rock damage based on the particle simulation method was proposed.The simulation results of three verification examples,for which the analytical solutions are available,demonstrate the correctness and usefulness of the thermo-mechanical coupled particle model.This model is applied to simulating an application example with two cases：one is temperature-independent elastic modulus and strength,while the other is temperature-dependent elastic modulus and strength.The related simulation results demonstrate that microscopic crack initiation and propagation process with consideration of temperature-independent and temperature-dependent elastic modulus and strength are different and therefore,the corresponding macroscopic failure patterns of rock are also different.On the contrary,considering the temperature-dependent elastic modulus and strength has no or little effect on the heating conduction behavior.Numerical results,which are obtained by using the proposed model with temperature-dependent elastic modulus and strength,agree well with the experimental results.This also reveals that the rock subjected to heating experiences much more cracking than the rock subjected to cooling.

Prediction of Sedimentary Microfacies Distribution by Coupling Stochastic Modeling Method in Oil and Gas Energy Resource Exploitation

In view of the problem that a single modeling method cannot predict the distribution of microfacies, a new idea of coupling modeling method to comprehensively predict the distribution of sedimentary microfacies was proposed, breaking the tradition that different sedimentary microfacies used the same modeling method in the past. Because different sedimentary microfacies have different distribution characteristics and geometric shapes, it is more accurate to select different simulation methods for prediction. In this paper, the coupling modeling method was to establish the distribution of sedimentary microfacies with simple geometry through the point indicating process simulation, and then predict the microfacies with complex spatial distribution through the sequential indicator simulation method. Taking the DC block of Bohai basin as an example, a high-precision reservoir sedimentary microfacies model was established by the above coupling modeling method, and the model verification results showed that the sedimentary microfacies model had a high consistency with the underground. The coupling microfacies modeling method had higher accuracy and reliability than the traditional modeling method, which provided a new idea for the prediction of sedimentary microfacies.

Study on damage-stress loss coupling model of rock and prestressed anchor cable in dry-wet environment

The loss of anchoring force is one of the problems to be solved urgently.The anchorage loss is a key factor causing the failure of anchoring engineering,so it is crucial to study the time-dependent variation of anchoring force.Alternating dry and wet(D-W)conditions have a significant effect on deformation of rock.The anchoring system is composed of anchoring components and rock mass,and thus rock deformation has a significant impact on the loss of anchoring force.Quantifying rock deformation under the effects of D-W cycles is a prerequisite to understanding the factors that influence loss of anchoring force in anchor bolts.In this study,we designed an anchoring device that enabled real-time monitoring of the variation in strain during D-W periods and rock testing.Nuclear magnetic resonance(NMR)measurements showed that under D-W conditions,the increment in porosity was smaller for prestressed rock than unstressed rock.The trends of prestress loss and strain variation are consistent,which can be divided into three characteristic intervals:rapid attenuation stage,slow attenuation stage and relatively stable stage.At the same stress level,the rate of stress loss and strain for the soaking specimen was the highest,while that of the dried specimen was the lowest.In the same D-W cycling conditions,the greater the prestress,the smaller the strain loss rate of the rock,especially under soaking conditions.The characteristics of pore structure and physical mechanical parameters indicated that prestress could effectively suppress damage caused by erosion related to D-W cycles.The study reveals the fluctuation behavior of rock strain and prestress loss under D-W conditions,providing a reference for effectively controlling anchoring loss and ideas for inventing new anchoring components.

Multi-dimensional Simulation of Phase Change by a 0D-2D Model Coupling via Stefan Condition

Considering phase changes associated with a high-temperature molten material cooled down from the outside,this work presents an improvement of the modelling and the numerical simulation of such processes for an application pertaining to the safety of light water nuclear reactors.Postulating a core meltdown accident,the behaviour of the core melt(aka corium)into a steel vessel is of tremendous importance when evaluating the vessel integrity.Evaluating correctly the heat fluxes requires the numerical simulation of the interaction between the liquid material and its solid counterpart which forms during the solidification process,but also may melt back.To simulate this configuration,encoun-tered in various industrial applications,one considers a bi-phase model constituted by a liquid phase in contact and interaction with its solid phase.The liquid phase may solidify in presence of low energetic source,while the solid phase may melt due to an intense heat flux from the high-temperature liquid.In the frame of the in-house legacy code,several simplifying assumptions(0D multi-layer discretization,instantaneous heat transfer via a quadratic temperature profile in solids)are made for the modelling of such phase changes.In the present work,these shortcomings are illustrated and further overcome by solving a 2D heat conduction model in the solid by a mixed Raviart-Thomas finite element method coupled to the liquid phase due to heat and mass exchanges through Stefan condition.The liquid phase is modeled with a 0D multi-layer approach.The 0D-liquid and 2D-solid mod-els are coupled by a Stefan like phase change interface model.Several sanity checks are performed to assess the validity of the approach on 1D and 2D academical configurations for which exact or reference solutions are available.Then more advanced situations(genu-ine multi-dimensional phase changes and an"industrial-like scenario")are simulated to verify the appropriate behavior of the obtained coupled simulation scheme.