Multiphysics Processes Research Articles

A New Performance-Based Test for Assessing Chloride-Induced Reinforcement Corrosion Resistance of Geopolymer Mortars

The widespread adoption of geopolymer concretes in the industry has been slow, mainly due to concerns over their long-term performance and durability. One of the main causes of concrete structures’ deterioration is chloride-induced corrosion of the reinforcement. The reinforcement corrosion process in concrete is composed of two main stages: the initiation phase, which is the amount of time required for chloride ions to reach the reinforcement, and the propagation phase, which is the active phase of corrosion. The inherent complexities associated with the properties of precursors and type of activators, and with the multi-physics processes, in which different transfer mechanisms (moisture, chloride, oxygen, and charge transfer) are involved and interact with each other, have been a major obstacle to predicting the durability of reinforced alkali-activated concretes in chloride environments. Alternatively, the durability of alkali-activated concretes can be assessed through testing. However, the performance-based tests that are currently available, such as the rapid chloride permeability test, the migration test or the bulk diffusion test, are only focusing on the initiation phase of the corrosion process. As a result, existing testing protocols do not capture every aspect of the material performance, which could potentially lead to misleading conclusions, particularly when involving an electrical potential to reduce the testing time. In this paper, a new performance-based test is proposed for assessing the performance of alkali-activated concretes in chloride environments, accounting for both the initiation and propagation phases of the corrosion process. The test is designed to be simple and to be completed within a reasonable time without involving any electrical potential.

Design and experimental study of solar-driven biomass gasification based on direct irradiation solar thermochemical reactor

Solar-driven biomass steam gasification is a promising technology to produce H2-rich syngas, meanwhile achieve flexible storage of renewable energy. However, the solar gasification application also faces numerous challenges, due to its involved complex chemical reaction characteristics and the inherent the multi-physics energy conversion processes. This work designs a novel direct irradiation solar gasification reactor, with the improved optical acceptance structure and reasonable test module, and the wheat straw particle is selected as experimental sample. Employing a 9 kWe high flux solar simulator, the maximum temperature of reaction bed in this prototype reactor reaches to 1260 °C, with the average solar flux of 1171.3 kW/m2. With adjustable radiation flux as experiment factor, under the highest total radiation of 3.35 kW, the molar fraction of H2 in the produced syngas is 47.1 % with the energy upgrade factor of 1.15, and the cellulose component completes decomposition. Increasing the mass flow rate of steam reduces the syngas output while raising the H2/CO ratio of the syngas. In addition, smaller biomass particle size facilitates the reaction, thereby improving the conversion efficiency of direct irradiation gasification. The investigation results provide a meaningful reference for the efficient utilization of abundant solar and biomass energy.

Using a multiphysics coupling-oriented flood modelling approach to assess urban flooding under various regulation scenarios combined with rainstorms and tidal effects

Climate change and rapid urbanisation have intensified urban flooding. Flood models constructed using hydrological and hydrodynamic approaches are crucial tools for simulating and regulating urban flooding processes. Urban flooding arises from the interplay of various physical processes in urban water systems, influenced by heterogeneous urban surfaces, underground structures, and intricate river networks influencing surface flow, which must be considered in flood models. For coastal cities with intricate drainage systems, the disasters of combined rainstorms and tidal flooding can be mitigated through river scheduling and real-time control of hydraulic structures. However, a lack of research on the effectiveness of flood regulations limits the application of flood models. This study proposes a flood modelling approach oriented towards the coupling of multiphysics processes, achieving a comprehensive simulation and regulation of the flood process. A surface flood control model was created by coupling hydrological and hydrodynamic processes using the Storm Water Management Model and the orthogonal storage cell model, enabled by data exchange synchronisation technology. Subsequently, the model’s accuracy was validated using observed rainstorm events, identifying drainage system overflows and land surface inundation. Finally, the developed model was applied to analyse the enhancing effect of tides on rainstorm-induced flooding and the impact of various control scenarios. This study offers theoretical direction for large-scale urban flood modelling and model development, offering significant references for revisiting urban planning and formulating flood prevention strategies under the joint impact of rainstorms and tides.

TC-SRNet: Reconstruction and prediction of typhoon high-resolution turbulent fields based on meteorological numerical forecast scale wind fields and deep learning method

Meteorological numerical forecast models can provide a more accurate typhoon wind fields compared to engineering typhoon models due to their incorporate atmospheric multiphysical processes. However, the latest advancements in supercomputing power indicate that the current finest level of real-time weather forecasting typically operates on grid scales ranging from 1 to 4 km, while the convergence of typhoon intensity and turbulent field characteristics occurs at scales as fine as 62–185 m. Therefore, the primary scientific inquiry lies in determining how to achieve high-precision turbulent wind fields while considering the realistic atmospheric multiphysical processes, that means establish a “bridge” of typhoon wind field between kilometer-level and hundred-meter scales. This study investigates the super-resolution reconstruction of wind fields across different horizontal grid scales, utilizing a benchmark wind field at a 62 m horizontal grid scale (ground truth), which is based on a hybrid down-sampling skip connection (DSC)/multi-scale (MS) model. The research findings demonstrate that compared to traditional interpolation methods, the DSC/MS method significantly improves reconstruction accuracy, albeit with some residual high-frequency energy dissipation issues. Additionally, the DSC/MS method currently exhibits better reconstruction performance for 62 m scale wind fields based on kilometer-scale and smaller horizontal grid scales (1 km, 555 m, 185 m), with improved reconstruction as grid scale decreases. However, significant errors are observed in reconstructing fine turbulent fields at 62 m scale based on wind fields at 1.67 km horizontal grid scale. The findings presented in the present study can provide real and high-precision turbulent wind fields for structural wind engineering and wind energy assessment studies, thereby holding significant scientific and engineering application value.

Multiphysical testing of strength development of cemented paste backfill containing superplasticizer

Current laboratory procedures for curing and testing the mechanical strength of cemented paste backfill (CPB) do not take into account the complex Multiphysics (thermal, T; hydraulic, H; mechanical, M; chemical, C) processes that CPB structures are subjected to in the field. This oversight can lead to unreliable measurements and unsafe designs. In this study, a multiphysical curing and testing procedures for CPB with superplasticizer (CPB-PES) has been developed to evaluate its strength development under THMC curing conditions close to those encountered in the field. The obtained results demonstrated that the strength development of CPB-PES is greatly affected by the THMC factors and their interactions. The contributions of each one of the investigated THMC factors on the strength are not equally similar and greatly depend on the interaction between these factors and curing time. CPB samples with 0.125 % PES that underwent drainage during THMC curing showed a strength increase of up to 809 % after 28 days of curing, compared to the control samples. The strength of CPB-PES samples cured under THMC can be up to 57 % higher than that of samples cured under THC conditions. The results indicate a significant interaction between thermal (T; elevated field curing temperature) and chemical (C; superplasticizer and cement hydration) factors, between chemical (C) and mechanical (M; field curing stress) factors, as well as between thermal and mechanical factors. The influence of the mechanical factor on strength development was observed to be less pronounced compared to the impact of chemical and thermal factor, and is reduced at elevated curing temperatures. The findings underscore the critical importance of accounting for field-relevant THMC factors and their interactions in the determination of the CPB-PES strength development, which is essential for the design of safer and more economical CPB structures, ultimately enhancing mine productivity and safety.

Gaussian smoothed particle hydrodynamics: A high-order meshfree particle method

Smoothed particle hydrodynamics (SPH) has attracted significant attention in recent decades, and exhibits special advantages in modeling complex flows with multiphysics processes and complex phenomena. Its accuracy depends heavily on the distribution of particles, and will generally be lower if the particles are distributed non-uniformly. A high-order SPH scheme is proposed in the present work for simulating both compressible and incompressible flows. It uses a Gaussian quadrature rule to perform the particle approximation of SPH by introducing Gaussian nodes. Unfortunately, the Gaussian nodes hardly overlap with SPH particles due to the Lagrangian feature, and thus we use a high-order interpolation method to obtain the corresponding physical quantities at the Gaussian nodes. The accuracy and robustness of the proposed Gaussian SPH are demonstrated by several numerical tests, including the Sod problem, Poiseuille flow, Couette flow, cavity flow, Taylor–Green vortex and dam break flow, and a convergence analysis is also conducted to evaluate the effects of particle resolution and distribution for reconstructing a given function. The simulation results for each test case are in good agreements with the available analytical, experimental or numerical results, showing that the proposed Gaussian SPH method is accurate and reliable but expensive for simulating compressible and incompressible flow problems.

Mathematical Modeling of Multiphysics Processes in Photoelectrochemical Water Oxidation at Near-Neutral pH

Mathematical Modeling of Multiphysics Processes in Photoelectrochemical Water Oxidation at Near-Neutral pH

State of the Art of Coupled Thermo–hydro-Mechanical–Chemical Modelling for Frozen Soils

AbstractNumerous studies have investigated the coupled multi-field processes in frozen soils, focusing on the variation in frozen soils and addressing the influences of climate change, hydrological processes, and ecosystems in cold regions. The investigation of coupled multi-physics field processes in frozen soils has emerged as a prominent research area, leading to significant advancements in coupling models and simulation solvers. However, substantial differences remain among various coupled models due to the insufficient observations and in-depth understanding of multi-field coupling processes. Therefore, this study comprehensively reviews the latest research process on multi-field models and numerical simulation methods, including thermo-hydraulic (TH) coupling, thermo-mechanical (TM) coupling, hydro-mechanical (HM) coupling, thermo–hydro-mechanical (THM) coupling, thermo–hydro-chemical (THC) coupling and thermo–hydro-mechanical–chemical (THMC) coupling. Furthermore, the primary simulation methods are summarised, including the continuum mechanics method, discrete or discontinuous mechanics method, and simulators specifically designed for heat and mass transfer modelling. Finally, this study outlines critical findings and proposes future research directions on multi-physical field modelling of frozen soils. This study provides the theoretical basis for in-depth mechanism analyses and practical engineering applications, contributing to the advancement of understanding and management of frozen soils.

Modeling and simulation of laser-ignited Al-PTFE reactive material in vacuum

Despite extensive investigations elucidating the initial response behavior of aluminized polytetrafluoroethylene (Al-PTFE) under atmospheric pressure, research on their response in lean-oxygen environments pertinent to propulsion and explosion scenarios is still lacking. This paper theoretically analyzes the multi-physical processes of laser-ignited Al-PTFE, identifying the initial response, main chemical reaction paths, and gas-phase product flow mechanisms. A dynamic reaction model of Al-PTFE considering material moving front, chemical reactions, gas-phase product dynamics, and variations in thermophysical properties, was developed. The robust coupling scheme in the model addresses nonlinear equations and fluid-solid coupling boundary conditions. Validation was achieved by comparing the ablation depth and plume number density field calculated by the model with the experimental results in the literature. The model reveals a positive correlation between laser fluence and both temperature growth rate and ablation depth of Al-PTFE, with the effects mainly observed in the full width at half maximum (FWHM) of the pulsed laser due to the thermal decomposition of the PTFE matrix. Simulations of the equilibrium constants and molar fractions of 19 components in the Al-PTFE ablation plume within the temperature range of 2000 K to 30,000 K at 0.001 Pa indicate that the equilibrium constants for dissociative reactions exceed those for ionization reactions, primarily in the early plume formation stages. Analysis of kinetic parameters and ambient pressure effects shows that during laser ablation, the maximum temperature of the plume initially appears near the ablation surface, shifting to the front, where the maximum velocity is always located. An elliptical shock wave forms due to adiabatic expansion of the plume, creating a cavity inside the plume. Lower ambient pressure increases plume expansion velocity, temperature, and size.

Multiphysics machine learning framework for on-demand multi-functional nano pattern design by light-controlled capillary force lithography

Nature finds ways to realize multi-functional surfaces by modulating nano-scale patterns on their surfaces, enjoying transparent, bactericidal, and/or anti-fogging features. Therein height distributions of nanopatterns play a key role. Recent advancements in nanotechnologies can reach that ability via chemical, mechanical, or optical fabrications. However, they require laborious complex procedures, prohibiting fast mass manufacturing. This paper presents a computational framework to help design multi-functional nano patterns by light. The framework behaves as a surrogate model for the inverse design of nano distributions. The framework’s hybrid (i.e., human and artificial) intelligence-based approach helps learn plausible rules of multi-physics processes behind the UV-controlled nano patterning and enriches training data sets. Then the framework’s inverse machine learning (ML) model can describe the required UV doses for the target heights of liquid in nano templates. Thereby, the framework can realize multiple functionalities including the desired nano-scale color, frictions, and bactericidal properties. Feasibility test results demonstrate the promising capability of the framework to realize the desired height distributions that can potentially enable multi-functional nano-scale surface properties. This computational framework will serve as a multi-physics surrogate model to help accelerate fast fabrications of nanopatterns with light and ML.

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

A Coupling Physics Model for Real-Time 4D Simulation of Cardiac Electromechanics

Cardiac simulators can assist in the diagnosis of heart disease and enhance human understanding of this leading cause of mortality. The coupling of multiphysics, such as electrophysiology and active–passive mechanics, in the simulation of the heart poses challenges in utilizing existing methodologies for real-time applications. The low efficiency of physically-based simulation is mostly caused by the need for electrical-stress conduction to use tiny time steps in order to prevent numerical instability. Additionally, the mechanical simulation experiences sluggish convergence when dealing with significant deformation and stiffness, and there are also concerns regarding volume inversion. We provide a coupling physics model that transforms the active–passive dynamics into multiphysics solving constraints, aiming at boosting the real-time efficiency of the cardiac electromechanical simulation. The multiphysics processes are initially divided into two levels: cell-level electrical stimulation and organ-level electrical-stress diffusion/conduction. This separation is achieved by employing operator splitting in combination with the quasi-steady-state method, which simplifies the system equations. Next, utilizing spatial discretization, we employ the matrix-free conjugate gradient approach to solve the electromechanical model, therefore improving the efficiency of the simulation. The experimental results illustrate that our simulation model is capable of replicating intricate cardiac physiological phenomena, including 3D spiral waves and rhythmic contractions. Moreover, our model achieves a significant advancement in real-time computation while maintaining a comparable level of accuracy to current methods. This improvement is advantageous for interactive medical applications.

Study on the fracture behavior in clayey geomaterials under moisture diffusion by phase field modeling

Clayey geomaterials are prone to cracking under changing moisture conditions, which substantially compromises their engineering properties. When moisture intrusion and boundary constraints act in tandem on geomaterials, the relevance of multi-physics field processes to complex cracking phenomena remains underexplored. To provide deeper insight into the effect of moisture transport on the cracking behavior of clayey geomaterials, a theoretical model for moisture-induced fracture is proposed. The model combines moisture diffusion, volume deformation dependent on moisture conditions, and phase field fracture, and is solved using the finite element method with a user-defined element subroutine. The developed subroutine is applied to simulate soil desiccation experiments and mudstone water absorption tests. The results reveal that moisture gradients induce non-uniform expansion or shrinkage within the clayey geomaterials. The neighboring geomaterial particles mutually constrain these deformations, resulting in an internal stress field that serves as the primary cause of cracking. Stability and sensitivity analyses validate the reliability of the numerical framework. The proposed model can effectively simulate crack initiation, extension, and convergence during humidification or desiccation processes in both two and three dimensions. This study charts a feasible path for numerical simulation and mechanical mechanism exploration of cracking behaviors in clayey geomaterials during water migration.

A conceptual model for evaluating the stability of high-altitude ice-rich slopes through coupled thermo-hydro-mechanical simulation

The collapse of glacial, permafrost, and ice-rich moraine slopes in high-altitude mountainous areas not only threatens downstream residents and infrastructure but also displaces ice mass to lower and warmer elevations, accelerating glacial ablation to some extent. Despite being considered rare, ice-rich slopes collapse more frequently than commonly thought due to climate change. For example, two adjacent mountain glaciers in the Aru Range of the Tibetan Plateau, characterized by large volumes (68 × 106 m3 and 83 × 106 m3, respectively) and low surface slope angles (12.3° and 12.9°, respectively), collapsed surprisingly in 2016. While the mechanisms behind these collapses have garnered broad attention, the ability to quantitatively assess the instability of ice-rich slopes remains limited due to the complex interplay of multi-physical processes. Taking the Aru glacier collapses as reference cases, this paper presents a conceptual model, implemented through coupled thermo-hydro-mechanical simulation, to evaluate the stability of high-altitude ice-rich slopes due to climate change, rainfall and ice ablation. Results indicate that the methodology captures the effects of temperature change, rainfall and meltwater on the instability events well, demonstrating promising potential in evaluating potential collapse zones of ice-rich slopes similar to the Aru glaciers. Furthermore, the role of climate change in the well-known Aru events is demonstrated using a state-of-the-art global climate reanalysis dataset. Findings reveal that the increase in liquid water infiltrating the Aru glaciers since 2010 was a critical factor leading to the instability events.

Two-phase flow thermo-hydro-mechanical modeling for a water flooding field case

Simulation of subsurface energy system involves multi-physical processes such as thermal, hydraulical, and mechanical (THM) processes, and requires a so-called THM coupled modeling approach. THM coupled modeling is commonly performed in geothermal energy production. However, for hydrocarbon extraction, we need to consider multiphase flow additionally. In this paper, we describe a three-dimensional numerical model of non-isothermal two-phase flow in the deformable porous medium by integrating governing equations of two-phase mixture in the porous media flow in the reservoir. To account for inter-woven impacts in subsurface conditions, we introduced a temperature-dependent fluid viscosity and a fluid density along with a strain-dependent reservoir permeability. Subsequently, we performed numerical experiments of a ten-year water flooding process employing the open-source parallelized code, OpenGeoSys. We considered different well patterns with colder water injection in realistic scenarios. Our results demonstrate that our model can simulate complex interactions of temperature, pore pressure, subsurface stress and water saturation simultaneously to evaluate the recovery performance. High temperature can promote fluid flow while cold water injection under non-isothermal conditions causes the normal stress reduction by significant thermal stress. Under different well patterns the displacement efficiency will be changed by the relative location between injection and production wells. This finding has provided the important reference for fluid flow and induced stress evolution during hydrocarbon exploitation under the environment of large reservoir depth and high temperature.

Analytical solution to non-isothermal pore-pressure evolution in hydrating minefill

Anomalous rises in minefill pressure during deposition shutdown have been increasingly reported in recent field investigations, raising growing concerns over contemporary underground tailings disposal. While constrained backfill expansion during hydration-heat release has been invoked to rationalize such perplexing observations, a benchmark for understanding the complex minefill response to thermal loading is still lacking. A novel analytical procedure is thus proposed in this study to characterize the non-isothermal evolution of pore pressure in cemented backfill hydrating in undrained adiabatic settings. By constructing a suitable variable system and selecting a proper auxiliary parameter for calculation, the proposed strategy based on canonical thermo-poroelasticity could enable straightforward implementation and has shown excellent compliance with computational and experimental results. The analytical results suggest that because higher temperature could stimulate cement hydration while enhancing simultaneously fluid expansivity, the pore pressure would decline initially with the rising curing temperature and then revert to increase after a critical threshold. Moreover, our investigation has shown that increasing binder usage is not always conducive to pore-pressure dissipation, as the payoff can be also hinged strongly upon the initial thermal condition. These findings could better elucidate the intriguing control of coupled multiphysics processes on backfill behaviour, and can also serve benchmarking purposes for sophisticated computational procedures.

Flexible and rigorous numerical modelling of multiphysics processes in fractured porous media using PorePy

Multiphysics processes in fractured porous media is a research field of importance for several subsurface applications and has received considerable attention over the last decade. The dynamics are characterized by strong couplings between processes as well as interaction between the processes and the structure of the fractured medium itself. The rich range of behaviour calls for explorative mathematical modelling, such as experimentation with constitutive laws and novel coupling concepts between physical processes. Moreover, efficient simulations of the strong couplings between multiphysics processes and geological structures require the development of tailored numerical methods.We present a modelling framework and its implementation in the open-source simulation toolbox PorePy, which is designed for rapid prototyping of multiphysics processes in fractured porous media. PorePy uses a mixed-dimensional representation of the fracture geometry and generally applies fully implicit couplings between processes. The code design follows the paradigms of modularity and differentiable programming, which together allow for extreme flexibility in experimentation with governing equations with minimal changes to the code base. The code integrity is supported by a multilevel testing framework ensuring the reliability of the code.We present our modelling framework within a context of thermo-poroelasticity in deformable fractured porous media, illustrating the close relation between the governing equations and the source code. We furthermore discuss the design of the testing framework and present simulations showcasing the extendibility of PorePy, as well as the type of results that can be produced by mixed-dimensional simulation tools.

Activated corrosion product contamination assessments of DEMO WCLL breeding blanket primary heat transport system

In water-cooled fusion reactors, the assessment of the primary system contamination is essential for waste management, machine availability, occupational radiation exposure, and radiological hazard determination. The primary cooling water is not only directly activated by the intense neutron field but is a contamination vector for a significant variety of gamma emitters with short to long decay half-lives. Corrosion products can be activated in those regions under neutron flux of the primary circuit and then released in the cooling water.In the EU-DEMO fusion power plant equipped with the Water-Cooled Lithium Lead Breeding Blanket (WCLL-BB) concept, the primary coolant undergoes intense neutron fields in the first wall and the breeding zone regions of the blanket. Activated Corrosion Products (ACPs) are then formed, released into the water, transported in the cooling loop and finally deposited onto the ex -vessel surfaces of the Primary Heat Transport System (PHTS), where working personnel are susceptible to being radiologically exposed.This work addresses the complete assessment of ACPs in the WCLL-BB PHTS of EU-DEMO. The simultaneous and multi-physical processes behind the ACP formation are tackled using the OSCAR-Fusion code, a comprehensive tool developed by the CEA (France) to assess contamination in fusion nuclear reactors. The whole system is modelled with zero-dimensional nodes with assigned geometrical, thermal-hydraulics, material and chemical parameters. Activation reaction rates integrated over the whole spectrum and calculated with MCNP are given to those regions exposed to the neutron flux. Results are provided in terms of mass and activity inventories of ACPs as deposit and inner oxide layers of components (pipes, heat exchangers, pumps...), ions in solution, particles in suspension, and filters and resins trapping.Mobilizable inventories such as ions, particles and deposits are important source terms in accidental scenario evolutions, while the whole activity inventory constitutes the main long-term gamma emitting source for dose rate maps determination in the tokamak building rooms housing the main PHTS equipment.

Development and verification of tritium transport code based on RELAP5/MOD3.3 with generic model towards COOL blanket

Tritium evaluation is a critical issue in fusion reactors. For the tritium self-sufficiency and safety concern, it is crucial to develop tritium analysis techniques to study its behavior in various conditions including steady states and dynamic states. The supercritical CO2 cOoled Lithium-Lead (COOL) blanket is an advanced blanket concept for the Chinese Fusion Engineering Testing Reactor (CFETR) and one of its major functions is to realize tritium self-sufficiency. Therefore, it is necessary to conduct comprehensive analyses for design and optimization from the perspective of tritium transport. In this work, the tritium transport module is integrated inside RELAP5/MOD3.3. Then the modified RELAP5 facilitates transient analyses to examine coupled thermal hydraulics and tritium transport responses. The code includes 1-D tritium permeation, tritium transport in fluid and tritium extraction, which are associated with thermodynamic parameters to cover the multi-physics processes of the COOL blanket. The modified RELAP5 is verified by the code-to-code benchmark. The results demonstrate the successful integration of the tritium transport capabilities in RELAP5. Furthermore, the code is applied to the heat exchanger and simplified PbLi loop to obtain tritium concentration to highlight its capability in both generic modeling and dynamic analysis.

Regularized coupling multiscale method for thermomechanical coupled problems

The coupling effects in multiphysics processes are often neglected in designing multiscale methods. The coupling may be described by a non-positive definite operator, which in turn brings significant challenges in multiscale simulations. In the paper, we develop a regularized coupling multiscale method based on the generalized multiscale finite element method (GMsFEM) to solve coupled thermomechanical problems, and it is referred to as the coupling generalized multiscale finite element method (CGMsFEM). The method consists of defining the coupling multiscale basis functions through local regularized coupling spectral problems in each coarse-grid block, which can be implemented by a novel design of two relaxation parameters. Compared to the standard GMsFEM, the proposed method can not only accurately capture the multiscale coupling correlation effects of multiphysics problems but also greatly improve computational efficiency with fewer multiscale basis functions. In addition, the convergence analysis is also established, and the optimal error estimates are derived, where the upper bound of errors is independent of the magnitude of the relaxation coefficient. Several numerical examples for periodic, random microstructure, and random material coefficients are presented to validate the theoretical analysis. The numerical results show that the CGMsFEM shows better robustness and efficiency than uncoupled GMsFEM.