Numerical Results Research Articles

Mesoscopic numerical simulation of chloride diffusion behavior in cracked recycled aggregate concrete

The cracking of recycled aggregate concrete (RAC) is well known to promotes the chloride diffusion, accelerates the corrosion of reinforcement embedded in RAC. To reveal the mechanism of chloride diffusion in RAC under cracking, a multiphase mesoscopic model for chloride diffusion in RAC was proposed. It should be noted that RAC is regarded as eight-phase composite materials consisting of coarse aggregate, reinforcement, new and old mortar, new and old interface transition zones (ITZ), cracks, and damage zones. The effects of the width and depth of cracks and damage zones on chloride diffusion behavior in RAC after cracking were further investigated. The numerical simulation results show that the damage zones accelerate the chloride diffusion and exacerbates the accumulation effect of chloride at the crack tip. Compared to the crack depth, the crack width of RAC has a small effect on chloride diffusion behavior, especially, the crack width is less than 50 µm. More importantly, the chloride diffusion streamline generated by numerical simulation reveals the mechanism of cracks promoting chloride diffusion. The research in this paper provides new insights into the durability design of RAC by revealing the diffusion behavior of chloride ions in RAC.

Efficient quantum state estimation with low-rank matrix completion

This paper introduces a novel and efficient technique for quantum state estimation, coined as low-rank matrix-completion quantum state tomography for characterizing pure quantum states, as it requires only non-entangling bases and 2n+1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$2n + 1$\\end{document} local Pauli operators. This significantly reduces the complexity of the process and increases the accuracy of the state estimation, as it eliminates the need for the entangling bases, which are experimentally difficult to implement on quantum devices. The required minimal post-processing, improved accuracy and efficacy of this matrix-completion-based method make it an ideal benchmarking tool for investigating the properties of quantum systems, enabling researchers to verify the accuracy of quantum devices, characterize their performance, and explore the underlying physics of quantum phenomena. Our numerical results demonstrate that this method outperforms contemporary techniques in its ability to accurately reconstruct multi-qubit quantum states on real quantum devices, making it an invaluable contribution to the field of quantum state characterization and an essential step toward the reliable deployment of intermediate- and large-scale quantum devices.

Resource Optimization for Quantum Dynamics with Tensor Networks: Quantum and Classical Algorithms.

The exponential scaling of the quantum degrees of freedom with the size of the system is one of the biggest challenges in computational chemistry and particularly in quantum dynamics. We present a tensor network approach for the time-evolution of the nuclear degrees of freedom of multiconfigurational chemical systems at a reduced storage and computational complexity. We also present quantum algorithms for the resultant dynamics. To preserve the compression advantage achieved via tensor network decompositions, we present an adaptive algorithm for the regularization of nonphysical bond dimensions, preventing the potentially exponential growth of these with time. While applicable to any quantum dynamical problem, our method is particularly valuable for dynamical simulations of nuclear chemical systems. Our algorithm is demonstrated using ab initio potentials obtained for a symmetric hydrogen-bonded system, namely, the protonated 2,2'-bipyridine, and compared to exact diagonalization numerical results.

Fractional-Order Boosted Hybrid Young’s Double-Slit Experimental Optimizer for Truss Topology Engineering Optimization

Inspired by classical experiments that uncovered the inherent properties of light waves, Young’s Double-Slit Experiment (YDSE) optimization algorithm represents a physics-driven meta-heuristic method. Its unique search mechanism and scalability have attracted much attention. However, when facing complex or high-dimensional problems, the YDSE optimizer, although striking a good balance between global and local searches, does not converge as fast as it should and is prone to fall into local optimums, thus limiting its application scope. A fractional-order boosted hybrid YDSE, called FYDSE, is proposed in this article. FYDSE employs a multi-strategy mechanism to jointly address the YDSE problems and enhance its ability to solve complex problems. First, a fractional-order strategy is introduced into the dark edge position update of FYDSE to ensure more efficient use of the search potential of a single neighborhood space while reducing the possibility of trapping in a local best. Second, piecewise chaotic mapping is constructed at the initial stage of the population to obtain better-distributed initial solutions and increase the convergence rate to the optimal position. Moreover, the low exploration space is extended by using a dynamic opposition strategy, which improves the probability of acquisition of a globally optimal solution. Finally, by introducing the vertical operator, FYDSE can better balance global exploration and local exploitation and explore new unknown areas. The numerical results show that FYDSE outperforms YDSE in 11 (91.6%) of cec2022 sets. In addition, FYDSE performs best in 8 (66.6%) among all algorithms. Compared with the 11 methods, FYDSE obtains the optimal best and average weights for the 20-bar, 24-bar, and 72-bar truss problems, which proves its efficient optimization capability for difficult optimization cases.

Financing and operational strategies for supply chains with yield uncertainty and capital constraints

Facing yield uncertainty is a common phenomenon that corporations encounter in their production and operations processes, especially small and medium-sized enterprises (SMEs). It is a challenge for SMEs that also have capital constraints, and thus financing from supply chain members or commercial banks is the normal way to improve the efficiency of supply chain operations. We consider a supply chain finance system consisting of an SME, a core enterprise (CE), and a bank, in which the SME is capital-constrained and faces productivity yield uncertainty. Therefore, the SME may use two financing sources to support his production: financing or advance payment from the CE. In the presence of bankruptcy risks for the SME, we model their strategic interaction as a Stackelberg game with the CE as the leader. We design the framework for the order penalty contract and derive the optimal strategies that include the planned production quantity decision of the SME and the pricing decisions of the CE under the two financing sources. We find that from the SME’s perspective, he will always prefer advance payments from the CE to bank financing. Furthermore, we propose a framework for partial credit guarantees. When it is difficult for banks to evaluate the credit history of SMEs, CEs provide guarantees to banks for financing SMEs. The analytical methods and numerical results provide decision-making methods and suggestions for production planning of the SME and wholesale prices formulated by CE.

A smart building energy management incorporating clustering-based tariffs in the presence of domestic solar energy, battery, and electric vehicle

Smart buildings play a crucial role in optimizing energy management within the power network. As end-users of the power network, they have the ability to not only reduce economic costs for householders but also modify the technical indices of the power network. To promote efficient device management in smart homes (SH), demand response programs are recommended for consumers. This research investigates the application of clustering-based electricity pricing strategy aimed at effectively managing the energy devices of a residential smart home. The utilized method categorizes the electricity tariff into five rates according to the clustering of the real-time pricing program. Ward’s clustering method is utilized to cluster and determine new electricity tariffs. The primary goal of the energy management program is to minimize the building’s energy cost, which is accomplished through the utilization of the multi-verse optimizer. The smart home consists of essential and manageable appliances, a photovoltaic panel (PV), a sodium-sulfur (NaS) battery, and an electric vehicle (EV). The initial parameters of the PV and EV are modeled stochastically by their probability distribution functions and calculated using the Latin hypercube sampling algorithm. The smart building’s performance is assessed by taking into account various demand response programs. The numerical results present that the application of the clustering-based management method has resulted in a significant reduction of 23–43 % in the electricity cost of smart homes. Additionally, the smart home exhibits a more linear consumption pattern when considering the electricity tariffs based on the clustering approach.

The effect of an auxiliary spring on the decrease of contact resistance and contact temperature for power automotive connector

Power automotive connectors used to connect two electrical systems in an electric vehicle should resist the rising of electric power or electric current. The latter can cause a rise in contact temperature by the Joule effect and can damage a connector improperly designed. The purpose of this work is to study experimentally and numerically the influence of the addition of a second spring on the ability of a power connector to transmit a high current. For this reason, an experimental test has been held on a connector without a second spring to measure contact resistance and contact temperature; then, a finite element model FEM of this connector was developed to compare and validate the numerical results with the experimental results. In addition, another FE model with a second spring was used to minimize contact resistance and contact temperature. Moreover, a theoretical model has shown its effectiveness in contact temperature and contact area calculations; it provides good agreement with experimental and numerical results. Numerical results show that the effect of the addition of a 2nd spring on the decrease of the contact temperature and contact resistance is significant. In addition, the thermo-electrical analysis by finite element shows that the current limit supported by a power connector with an auxiliary spring will equal twice the current supported by a power connector with one spring. FE mechanical analysis shows that the maximum mechanical stress in the connector with an auxiliary spring remains lower than the 2nd spring material yield stress. Finally, experimental and numerical results are in good agreement.

Theoretical investigation of optoelectronic properties in PbS/CdS core/shell spherical quantum dots under the effect of the electric field intensity, hydrogenic impurity and geometric parameters

In this paper, we theoretically investigated the combined effects of the external electric field (EF) strength and the geometric parameters on the linear, nonlinear and total dielectric functions as well as the effective dielectric function (DF) coefficients of PbS/CdS spherical core/shell quantum dots (CSQD) in the presence of the hydrogenic impurity located at the center. The subband energy eigenvalues and their corresponding wave functions are obtained by solving the time-independent Schrödinger equation and using the variational method (VM) in the framework of the effective mass approximation (EMA). The linear, nonlinear and total DF as well as the effective DF were discussed and evaluated under the influence of EF intensity, geometric parameters, the change of the number of quantum dots (QDs) per unit volume and optical intensity I based on the compact density matrix (CDM) approach. The obtained results show that the height peaks of the linear, nonlinear and Total DF as well as the effective DF coefficients increase and their resonant peak moves towards lower energies as the EF and geometric factors increase in both cases with and without considering the hydrogenic impurity effect. Furthermore, our findings show that the impact of optical intensity and the number of QDs per unit volume does not change the resonance peaks of imaginary parts of DF and effective DF but decreases their magnitudes. As a result, we believe that numerical results will present important developments and provide great contributions in designing new optoelectronic devices related to CSQD hetero-nanostructure.

Study on the Multiphase Flow Behavior in Jet Pump Drainage and Natural Gas Hydrate Production Wells with Combined Depressurization and Thermal Stimulation Method

Natural gas hydrate (NGH) trials have been performed successfully with different development methods and gas recovery drainage technologies. Multiphase flow in a wellbore and the drainage of natural gas hydrate are two important parts for its whole extraction process. Additionally, the choice of the drainage method is linked to the development method, making the drainage of NGH more complex. Jet pump drainage is usable for NGH production wells with the combined depressurization and thermal stimulation method. The objective of this study is to shed more light on the multiphase flow behavior in jet pump drainage and NGH production wells and put forward suggestions for adjusting heat injection parameters. The mechanism of jet pump drainage recovery technology for NGH wells was analyzed and its applicability to NGH development by the combined depressurization and thermal stimulation method was demonstrated. In addition, multiphase flow models of tubing and annulus were established, respectively, for the phenomenon of the countercurrent flow of heat exchange in the process of jet pump drainage and gas production, and the corresponding multiphase flow laws were derived. On the basis of these studies, sensitivity analysis and the optimization of thermal stimulation parameters were conducted. It is demonstrated that jet pump drainage gas recovery technology is feasible for the development of onshore NGH with the combined depressurization and thermal stimulation method. The laws of multiphase flow in the tubing and annulus of jet pump drainage and NGH production wells were disclosed in this study. Numerical simulation results show that the temperature and pressure profiles along the wellbore of jet pump drainage and NGH production wells during the drainage recovery process are affected by injection conditions. Increasing injection rate and injection temperature can both improve the effect of heat injection and reduce the hydrate reformation risk in the bottom of the annulus. This study offers a theoretical basis and technical support for production optimization and hydrate prevention and control in the wellbore of jet pump drainage and NGH production wells.

Meshing theory and contact analysis of double enveloping hourglass worm drive with planar generatrix

The meshing and limit equations of worm drive usually have strong nonlinearities such as multiple solutions, solution nonexistence and equation singularity. Meanwhile, the tooth surfaces of worm drive are in nonconforming line contact, which often requires mesh refinement of contact region for the loaded contact analysis. These two challenges cause that the modeling of worm drive heavily relies on manual adjustment and the loaded contact analysis of worm drive is still rare especially when edge contact and assembly error are concerned. Focusing on the double enveloping hourglass worm (DEHW) drive with planar generatrix, this work presents procedures to solve meshing and limit equations with global convergence. The instantaneous contact line, meshing limit line, curvature interference limit line and tooth surface grid discretization are adaptively generated, without manual trial or adjustment. On the basis of adaptive mesh refinement of tooth surface, the mortar virtual element method is adopted for loaded contact analysis of DEHW drive with edge contact and center distance error. Under sliding friction, the discrete system of governing equalities and inequalities is solved by semi-smooth Newton algorithm after constraint condensation. Numerical results for meshing theory and loaded contact analysis of DEHW drive with planar generatrix are discussed.

Numerical modeling and validation of dike-induced water flow dynamics using OpenFOAM

ABSTRACT The basic purpose of the present research is to numerically elucidate the flow around a single dike in four different geometric configurations (CN: no cylinders, C5: 5 cylinders, C10: 10 cylinders, C15: 15 cylinders), under conditions of a constant flow rate and subcritical flow. This involved substituting the impermeable dike with varying numbers of piles to validate the observed experimental flow patterns, particularly the conditions leading to reverse flow formation. OpenFOAM, an open-source algorithm, was utilized for the simulation, incorporating the volume of fluid (VOF) method. Accurately depicting the reverse flow formation with minimal numerical diffusion, a strategy involving multizone meshing and the application of mesh refinement zones was utilized. These refinements were particularly focused on the section between 1.8 m and 2.2 m within the numerical domain. The water surface profile in front section of the dike and the averaged velocity at five different locations, which were then compared with the numerical results. The numerically predicted streamwise velocity showed a percentage error ranging between 0.5% and 4.5%. The numerical results from this study are pivotal in dike design, aiming to mitigate accelerated wear during flood events and to protect both the dike head and the adjacent bank from high energy flow failures.

Multi-layered solid element for seismic analysis of rubber bearings considering stress continuity based on finite particle method

Rubber bearings are crucial components in seismic isolation systems to mitigate the damaging effects of earthquakes on structures. Rubber bearings generally consist of alternate layers of steel and rubber bonded together. Simplified two-point spring models omit this layered construction and have limitations in adaptability and programmability, while refined multiple solid models lead to substantial computational costs. This study presents an effective and efficient approach to simulate rubber bearings by using multi-layered solid element based on the finite particle method (FPM). The traditional “Zigzag” shape function is extended into the 3D situation to reflect the serrated deformation within the layered construction. Based on the simplified assumption of stress continuity, an iterative approach is presented for the deformation gradient to achieve the force continuity between layers in dynamic nonlinear scenarios. The multi-layered solid models for both natural rubber bearings (NRBs) and lead rubber bearings (LRBs) are proposed based on this element. The errors of the hysteresis curves of rubber bearings are below 3.1 % between the numerical results obtained by using the multi-layered solid model and the experimental results. The computational speedup ratio of the proposed multi-layered solid model is verified to be 19 relative to the refined multiple solid model. The proposed method is applied to a base-isolated frame structure to demonstrate its effectiveness in structural seismic analysis.

A length-free phase field model for epoxy adhesive materials based on modified Drucker–Prager criterion

Characterizing the cracking behavior of epoxy adhesive material is of great importance for composite structures. Recently, the phase field model has emerged as the preferred method for predicting the crack propagation. This paper presents a length-free phase field model considering the characteristics of epoxy adhesive materials. The proposed model is formulated within a framework that decomposes the phase field fracture driving force, ensuring that it adequately considers the tension/compression asymmetry of the epoxy adhesive materials. The effects of hydrostatic pressure, deviational stress, and friction are incorporated into the model by employing a modified Drucker-Prager failure criterion. The failure criterion is included in the phase field model using the effective stress invariants. Furthermore, the internal length scale, previously considered to lack real physical meaning, is reformulated based on material strength. This model addresses the challenge of measuring the internal length of epoxy materials. It serves as a valuable reference for applying phase field models to epoxy materials. The accuracy of the model is validated through an L-shaped plate simulation, and an attempt using the model to predict the crack development in thick epoxy adhesive joints is undertaken. Several representative samples are simulated, including three double cantilever beams with different configurations. The numerical results demonstrate a good correlation with experimental data, highlighting the potential of the model in simulating the crack propagation in epoxy adhesive joints.

Experimental and Numerical Investigation of Impact Behaviour of Anchored CFRP-to-Concrete Bonding Interface

In the last 20 years, Carbon Fiber Reinforced Polymers (CFRP) have become a common alternative to strengthening and retrofitting. The performance of the strengthened specimens is closely related to the bond slip between the CFRP strips and the adhesion surface. Therefore, researchers have begun to concentrate on the various anchoring methods to delay the debonding of the CFRP from the surface of strengthened specimens. When the literature review of the adhesion behavior between the anchored CFRP strips and the concrete surface is examined, it has been shown that the investigations are mostly done under the effect of static or reversible cyclic loads. In the literature review, there is not any research related to the adhesion behavior between anchored CFRP strips and the concrete surface under the effect of sudden dynamic impulsive impact loading. Therefore, the variables examined in the experimental study are the anchor type used in the CFRP strips, CFRP strip adhesion length, and the number of anchors placed on the strip. As a result of the experiments, the acceleration-time, displacement–time curves of the beams, and the strain–time curves along the CFRP strips were interpreted. In addition, simulation models of the test specimens were created with ABAQUS software, and non-linear finite element analyzes were performed. Finally, by comparing numerical results with experimental ones, the effects of variables such as bond length, applied anchor types and numbers on the impact behavior of the beam and the strain distributions along the CFRP strips were interpreted.

Numerical Analysis of the Cylindrical Shell Pipe with Preformed Holes Subjected to a Compressive Load Using Non-Uniform Rational B-Splines and T-Splines for an Isogeometric Analysis Approach

In this paper, we implement the finite detail technique primarily based on T-Splines for approximating solutions to the linear elasticity equations in the connected and bounded Lipschitz domain. Both theoretical and numerical analyses of the Dirichlet and Neumann boundary problems are presented. The Reissner–Mindlin (RM) hypothesis is considered for the investigation of the mechanical performance of a 3D cylindrical shell pipe without and with preformed hole problems under concentrated and compression loading in the linear elastic behavior for trimmed and untrimmed surfaces in structural engineering problems. Bézier extraction from T-Splines is integrated for an isogeometric analysis (IGA) approach. The numerical results obtained, particularly for the displacement and von Mises stress, are compared with and validated against the literature results, particularly with those for Non-Uniform Rational B-Spline (NURBS) IGA and the finite element method (FEM) Abaqus methods. The obtained results show that the computation time of the IGA based on the T-Spline method is shorter than that of the IGA NURBS and FEM Abaqus/CAE (computer-aided engineering) methods. Furthermore, the highlighted results confirm that the IGA approach based on the T-Spline method shows more success than numerical reference methods. We observed that the NURBS IGA method is very limited for studying trimmed surfaces. The T-Spline method shows its power and capability in computing trimmed and untrimmed surfaces.

Optimizing thermal efficiency: Advancements in flat plate heat exchanger performance through baffle integration

Compact heat exchangers have been gaining popularity in many industrial applications. Various types of passive turbulising structures such as corrugations, protrusions and ribs are introduced in the flow path to increase the effective heat transfer area and the level of turbulence in the flow path. This study investigates the impact of introduction of baffles on the performance of PHEs in terms of flow characteristics, pressure drop, and heat transfer. Two distinct types of baffle structures, namely wedge and aerofoil configurations, were introduced at varying numbers - 1, 3, and 5. An extensive experimentation is conducted for a FPHE in a thermal power plant of 500 × 2 MW and various flow and thermal parameters are measured. Computational Fluid Dynamics is utilized in this study to find the optimum baffle configuration. A detailed validation study is executed to obtain the correct computational algorithm, that is the right mesh count, optimum turbulence model, and precise numerical algorithm by comparing the numerical results with the available experimental results. Wedge- type baffles create increased turbulence and pressure drop, while aerofoil-type baffles minimize stagnation and exhibit lower pressure drop. Both baffle configurations lead to a substantial increase in heat transfer, with the 5-wedge-baffle setup showing the highest up to a 55% enhancement of Nusselt number. The Performance Evaluation Criterion (PEC) of wedge and aerofoil type is about 1.24 to 1.3 and 1.22 to 1.24 to that of conventional one respectively.

Electroosmotic flow modulation and dispersion of uncharged solutes in soft nanochannel.

We perform a systematic study on the modulation of electroosmotic flow (EOF), tuning the selectivity using electrolyte ions and hydrodynamic dispersion of the solute band across the soft nanochannel. The supporting walls of the channel are considered to be hydrophobic and bear non-zero surface charge. For such a channel, the inner side of the supporting rigid walls of the channel are coated with a soft polyelectrolyte layer (PEL). The inhomogeneous distribution of monomers and accompanying volume charge within the PEL is modelled via soft-step function. The dielectric permittivity of the PEL and electrolyte solution are in general different, which in turn leads to the ion partitioning effect. The impact of ion steric effects due to finite sized ions is further accounted through the modified ion activity coefficient. To model the EOF modulation considering the combined impact of the ion steric and ion partitioning effects as well as inhomogeneous distribution of monomers across the PEL, we adopt the modified Poisson-Boltzmann equation as the governing equation for electrostatic potential. The Debye-Bueche model is adopted to study the flow field across the PEL and the Stokes equation governs the EOF outside the PEL. In order to study the impact of the modulated EOF field on the dispersion of uncharged solution, we adopt three different models, i.e., a general 2D convective-diffusion model as well as cross-sectional averaged dispersion models due to Gill and late-time Taylor and Aris. Going beyond the widely employed Debye-Hückel approximation and uniform distribution of the monomer as well as accompanying volume charge, we find the results for the electric double layer (EDL) potential, EOF field and averaged throughput, by tuning the ion selectivity, etc., which is sufficient to analyze the transport of ionized liquid across the channel. The numerical results are supplemented with analytical results for the EDL potential as well as the EOF field under various limiting situations. Besides, we have further shown the impact of the modulated EOF field on the solute dispersion process. We have presented results that highlight the impact of parameters related to EOF field modulation, on solute dispersion governed by a convective-diffusive process, as well as obtaining the results for an effective dispersion coefficient. The dispersion models under the modulated EOF field adopted in the present study can thus be applied to study the dispersion process in engineered microdevices.

Numerical solution of fuzzy stochastic Volterra integral equations with constant delay

This paper considers the nonlinear fuzzy stochastic Volterra integral equations with constant delay, which are general and include many fuzzy stochastic integral and differential equations discussed in literature. Since Doob's martingale inequality is no longer applicable to such equations, a new maximum inequality is obtained. Combining with the Picard approximation method, the existence and uniqueness of solutions to nonlinear fuzzy stochastic Volterra integral equations with constant delay are given. Moreover we prove that the solution behaves stably in the case of small changes of initial values, kernels and nonlinearities. We further develop a Euler-Maruyama (EM) scheme and prove the strong convergence of the scheme. It is shown that the strong convergence order of the EM method is 0.5 under Lipschitz condition. Moreover, the strong superconvergence order is 1 if further, the kernel h(t,s) of the stochastic term satisfies h(t,t)=0. Numerical examples demonstrated that the numerical results are consistent with the theoretical research conclusions. Furthermore, the application model of the fuzzy stochastic Volterra integral equation with constant delay in population dynamics is considered, and the exact solution of the numerical example is given in explicit form.

A well-balanced finite volume solver for the 2D shallow water magnetohydrodynamic equations with topography

In this paper, a second-order finite volume Non-Homogeneous Riemann Solver is used to obtain an approximate solution for the two-dimensional shallow water magnetohydrodynamic (SWMHD) equations considering non-flat bottom topography. We investigate the stability of a perturbed steady state, as well as the stability of energy in these equations after a perturbation of a steady state using a dispersive analysis. To address the elliptic constraint ∇⋅hB=0, the GLM (Generalized Lagrange Multiplier) method designed specifically for finite volume schemes, is used. The proposed solver is implemented on unstructured meshes and verifies the exact conservation property. Several numerical results are presented to validate the high accuracy of our schemes, the well-balanced, and the ability to resolve smooth and discontinuous solutions. The developed finite volume Non-Homogeneous Riemann Solver and the GLM method offer a reliable approach for solving the SWMHD equations, preserving numerical and physical equilibrium, and ensuring stability in the presence of perturbations.

Emergoenvironmental analysis and sustainable flexible peaking of a novel solar-assisted solid oxide fuel cell cogeneration system

The development of green and sustainable conversion technologies is imperative to alleviate the current severe energy and environmental crises and to achieve the goal of carbon neutrality. This project presents a novel cogeneration system that integrates a solid oxide fuel cell afterburning system, an energy recovery system, and a parallel heating system equipped with multistage utilization, solar-assisted, and flexible peaking technologies. Following the establishment and validation of the system, power source analysis, parameter sensitivity studies, three-objective optimization, and internal exergy flow investigation are carried out. Subsequently, the positive impacts of multistage utilization, solar-assisted, and flexible peaking technologies, as well as the thermoeconomic and emergoenvironmental operating characteristics of the system, are discussed. The numerical results demonstrate that the emergoenvironmental factor and emergy environmental impact difference of the system are 70.15 % and 82.52 %, and the energy efficiency is improved by 1.41 %∼19.09 %, respectively, compared with similar solid oxide fuel cell systems, confirming its favorable thermodynamic performance and emergoenvironmental impact. The system achieves daily heating savings of 3302 kWh and life cycle economic benefits of $19,310,524, with a repayment of the investment cost in 4.17 years, reflecting considerable economic benefits. This project also demonstrates exergy flow, module defects, and improvement mechanisms from the perspective of emergoenvironmental, providing a valuable reference for the promotion of green energy.