Simulation Research Research Articles

A new evaluation method for sand control performance of slotted screen

Evaluation method of slotted screen performance is a critical concern for sand control. Despite its significance, comprehensive investigations into quantitative study of sand retaining during the complex flow process need further in-depth study, particularly in numerical simulation research model of the slotted screen. This paper adopts the two-way coupling calculation method of computational fluid dynamics (CFD) and discrete element method (DEM) for the particle movement and fluid flow at the slotted screen and establishes a numerical model of the sand particle migration process through the slotted screen. The sand particle deposition and plugging dynamics on the surface of the sand retaining medium and the micro sand plugging mechanism and blockage pattern are studied. The fractal dimension of formation sand is introduced to characterize different distribution of formation sand. The effects of fluid viscosity, width of the slotted screen, sand production concentration, fluid velocity and particle size distribution on sand production, sand passing rate and particle size distribution of formation sand passing through the slot and plugging behavior are studied. The sand plugging process of slotted screen can be divided into three stages: the beginning of blockage, the intensification of blockage, and the balance of blockage. Under simulated conditions, the initial sand plugging rate is 87%, and as the blockage intensifies, the sand plugging rate gradually increases to 96.5%. Based on the analysis of numerical simulation results, the index of total seepage resistance damage amplitude and velocity, general sand retention ratio and passed-sand size were defined to construct the integrated evaluation method of slotted screen media, which involves the performance of sand retention, anti-plugging and flow ability. A brand new set of ideas on quantifying the integrated evaluation method were introduced to calculate the individual and integrated performance indicators of the slotted screen based on simulated data. The quantitative evaluation of the sand control performance of the slotted screen provides guidance and reference for the parameter design of the slotted screen in the target oil region and the accuracy optimization of the sand control medium.

Research on Protective Effect of Lead Buffer Device under Explosion Load

In modern society, the frequency of explosion accidents is gradually increasing, which brings a great threat to the safety of people's lives and property and social stability. Explosive loads are extremely destructive, with huge energy and powerful impact force released in an instant, capable of causing serious damage to buildings and infrastructure. Therefore, it is of great practical significance to study explosion protection measures in depth and improve the damage resistance of buildings and infrastructure under explosion loads. The traditional underground dispersed layer can reduce the impact of blast load on the structure to a certain extent, but with the continuous improvement of protection requirements, its limitations are gradually revealed. For example, the protection of underground dispersed layers may be less than ideal in some complex explosion situations, and their performance is susceptible to environmental factors such as salinity, consolidation, precipitation and ice, resulting in reduced protection effectiveness. In order to find a more effective solution for explosion protection, the idea of replacing the underground dispersion layer with lead dampers was proposed. As a new type of protective device, lead damper has unique advantages. Lead has good plastic deformation ability, and can absorb a large amount of energy through its own deformation under the action of explosion load, so as to effectively reduce the impact of explosion on the structure. In addition, lead dampers can be customized to meet specific engineering needs, allowing them to better adapt to different explosion scenarios and structural types. The purpose of this study is to further investigate the protective effect of lead dampers instead of underground dispersed layers under blast loads. Through a combination of theoretical analysis, numerical simulation and experimental research, the performance of the lead damper is comprehensively evaluated, and its response mechanism under blast load is analyzed, as well as the advantages and disadvantages compared with the underground dispersed layer. At the same time, the influence of the design parameters of the lead damper on the protection effect will be studied, so as to provide scientific basis and guidance for its application in practical engineering. Through this study, it is expected to provide a new and more effective means of protection for explosion protection engineering, and improve the safety and reliability of buildings and infrastructure under explosion loads. This will not only help reduce the losses caused by explosion accidents and ensure the safety of people's lives and property, but also contribute to the stability and sustainable development of society.

An Effective Exponential Ratio and Product Type Estimator in Post-Stratification

In this paper we present a unique estimator for the population parameter under post-stratification. By segmenting the population into homogeneous subgroups, post-stratification is a widely used strategy in survey sampling that increases the accuracy of estimators. In the post-stratification scenario, this paper addressed the issue of estimating the mean of a finite population. For post- stratification, better separate ratio and product exponential type estimators are proposed. This study's primary goal is to evaluate our suggested estimators' performance against that of current estimators. We perform a thorough simulation research to assess the precision and effectiveness of our estimator. The mean squared errors and biases of the proposed estimators are obtained to the first degree of approximation. Theoretical and practical researches have demonstrated that the proposed estimators are more efficient than other estimators that were taken into consideration.

Physical simulation of downburst winds for civil structures: A review

Complementary to field measurement and computational simulation, experimental generation of downburst winds in laboratories is an effective approach to enhance the understanding of the wind characteristics and the induced loads on civil structures. This study reviews the literature in physical simulation of downburst winds, while capturing historic background, highlighting recent advances, and illuminating future trends. A total of 23 testing facilities, classified into three categories of standalone simulators (12 facilities), modified boundary layer wind tunnels (five facilities), and multi-fan wind tunnels (six facilities), are discussed with their key features and representative studies. This paper serves as a valuable reference for future research in physical simulation of downburst winds.

Numerical Simulation and Experimental Research on Effect of Flow Regimes for Refrigeration Performance of Wave Rotor

Abstract Wave rotors exchange energy through the interaction of fluids with different pressures to produce a cooling effect without the application of mechanical parts, which makes them suitable for the complex flow field environments. But the flow situation inside the wave rotor is complex, and research on the impact of various flow regimes for the refrigeration efficiency is currently incomplete. This study deals with the numerical simulation and the experimental verification of a four-port fixed rotor gas wave refrigerator. Various typical flow regimes within the wave rotor were obtained through simulation to analyze their impact on the refrigeration efficiency. To further demonstrate the effectiveness of the simulation results, experimental verification was conducted. The simulation results indicate that separation bubbles and vortices will be respectively generated on the upper and lower walls of the wave rotor tube during the gradual opening stage. Vortex affects the airflow flow and causes detachment of separation bubbles as the airflow moves, which results in the loss of pressure energy and kinetic energy of the incident gas. As a result, the loss caused by the gradual opening process leads to a decrease in the self-expansion ability of high-pressure gas during the gradual closing stage, thereby affecting the refrigeration effect. Therefore, reducing the influence of vortices plays an important role in improving the cooling efficiency of wave rotors. Finally, this study analyzes the main factors that affect the cooling effect within the wave rotor through flow regime and provides new insights for improving refrigeration efficiency.

Study on the formation and characteristics of unstable sheet cavitation on hydrofoils

Hydrofoils, as fundamental components of hydraulic machinery, directly influence the performance of such machinery. This study conducted an analysis of the cavitation characteristics of hydrofoils at a + 4° angle of attack under various cavitation numbers using numerical simulation and experimental research methods. The focus of the research was to explore the phenomenon of unstable sheet cavitation and its causes, as well as to reveal the characteristics of the re-entrant jet. The large eddy simulation method was employed to calculate the cavitation morphology under three different cavitation numbers. The method is highly consistent with the experimental results in simulating the small-scale detachment at the tail of unstable sheet cavitation and the large-scale shedding of cloud cavitation. The study found that the detachment of unstable sheet cavitation is closely related to the re-entrant jet, which exhibits transient and abrupt characteristics during the unstable sheet cavitation phase. Furthermore, by applying FFT processing to the distribution of maximum reverse velocity and the spatiotemporal changes of Ux on characteristic lines, eigenfrequency of the detachment of unstable sheet cavitation were identified. The research results indicate that cavitation mainly show as sheet cavitation when the cavitation closure point does not exceed the zero-slope point. Beyond this point, it transitions to cloud cavitation. This study provides new insights into the cavitation phenomenon of hydrofoils and offers quantitative research on the phenomenon of unstable sheet cavitation.

Research Progress in Spatiotemporal Dynamic Simulation of LUCC

Land Use and Land Cover Change (LUCC) represents the interaction between human societies and the natural environment. Studies of LUCC simulation allow for the analysis of Land Use and Land Cover (LULC) patterns in a given region. Moreover, these studies enable the simulation of complex future LUCC scenarios by integrating multiple factors. Such studies can provide effective means for optimizing and making decisions about the future patterns of a region. This review conducted a literature search on geographic models and simulations in the Web of Science database. From the literature, we summarized the basic steps of spatiotemporal dynamic simulation of LUCC. The focus was on the current major models, analyzing their characteristics and limitations, and discussing their expanded applications in land use. This review reveals that current research still faces challenges such as data uncertainty, necessitating the advancement of more diverse data and new technologies. Future research can enhance the precision and applicability of studies by improving models and methods, integrating big data and multi-scale data, and employing multi-model coupling and various algorithmic experiments for comparison. This would support the advancement of land use spatiotemporal dynamic simulation research to higher levels.

Modeling and simulation study of end around taxiway operation of Shanghai Hongqiao Airport

Abstract With the rapid expansion of the aviation industry, an increasing number of Close Spaced Parallel Runway (CSPR) airports are either planning or constructing End Around Taxiways (EAT) to alleviate field operation pressures and enhance safety. Taking Shanghai Hongqiao Airport’s typical CSPR EAT configuration as a case study, this research integrates the airport’s current operational status with the anticipated requirements for future structural renovations and increased flight volumes. Various operational scenarios are established, and simulation research on optimising EAT operations is conducted in advance. The simulation study proceeds as follows: first, an AirTOP simulation model is constructed based on Hongqiao Airport’s actual operational construction. Subsequently, leveraging existing operational scenarios, five simulation scenarios are devised by activating EATs at the departure and approach ends of the eastern zone. The merits and drawbacks of these scenarios are thoroughly analysed. The findings indicate that, with escalating flight volumes, the utilisation of EAT for larger aircraft can curtail their holding duration by nearly 8 min, consequently reducing overall arrival holding duration by 6 min. Departures from gates proximate to T1 experience a 3-min reduction in holding duration through the adoption of EAT at the approach end. Despite an increase in taxi distance due to a higher proportion of aircraft taxiing around, the overall taxi time is diminished. Activating EATs at the departure and approach ends of the eastern zone effectively mitigates the adverse effects of heightened flight volumes on field operational efficiency.

A review of resistance–capacitance thermal network model in urban building energy simulations

With the increasing population and urbanization promotion, energy consumption and carbon emissions have increased, and concern for inefficient energy use and deteriorating urban environments is growing. The prediction of building energy consumption and environmental evaluation, especially at large scales, are considered to be a major challenge confronting the research community. An outstanding strategy for mitigating energy consumption and carbon emissions resides in the field of energy modeling. As a simplified building energy modeling model, resistance–capacitance (RC) network model has the applications in fast predicting building energy consumption. Notably, in recent years, there has been an evident absence of thorough review endeavors related to RC model for urban building energy loads and climate. This review systematically explores the application of low-order models and reduction methods for urban or regional energy simulation starting from the typical RC model for building. In summation, the challenges associated with employing this model for urban building energy loads and climate can be succinctly summarized as follows: the absence of a unified platform for RC modeling; insufficient readily applicable urban datasets; the need for modeling tools that facilitate cross-platform analyses in conjunction with existing Geographic Information Systems (GIS) and urban thermal environment simulation research platforms, and the essential requirement for seamless integration with other complementary modules.

Large-scale molecular dynamics simulation and aggregate behavior research on asphalt

In this comprehensive study, a molecular dynamics (MD) model consisting of approximately 150,000 atoms was developed within the Strategic Highway Research Program (SHRP) framework to investigate the behavior of AAA asphalt, focusing on asphaltene aggregate dynamics behavior over a simulation time span of 0.282 microseconds. The findings reveal that the model achieves a stable state in terms of both energy and microstructure without significant changes or the formation of densely clustered molecules, despite the extensive simulation period. Notably, the patterns of asphaltene aggregation were found to significantly influence the MD system energy, with variations ranging from loosely organized structures due to dispersed asphaltenes to more structured, energy-efficient clusters formed by encapsulated asphaltene pentamers. The size of asphaltene aggregates emerged as a pivotal factor, affecting their encapsulation and leading to distinct formation patterns. Furthermore, the research highlighted that increasing shear rates prompt the depolymerization of asphaltene aggregates, shifting from large-scale structures to smaller ones and adapting dynamically to the flow conditions by aligning parallel to the shear direction. This study underscores the critical impact of asphaltene aggregate behavior and shear rates on the structural integrity and distribution within the asphalt model, offering profound insights into the MD of asphalt materials.

Simulating cyberattacks with extended Petri nets

Cybersecurity is an urgent concern. Cybersecurity simulation is an important part of the response to it. This article describes a research program consisting of several interconnected cybersecurity simulation research projects. Cyberattacks are modeled using Petri nets extended with features designed for modeling cyberattacks, including representations of the attacker’s and defender’s strategies, their actions, and their actions’ cost. A database of known attack patterns is automatically processed to generate cyberattack component models, one for each attack pattern. The models are verified and validated using multiple application-relevant methods that consider both Petri nets’ theoretical properties and cyberattacks’ practical characteristics. Because the source attack pattern database is attacker-centric, the cyberattack component models are enhanced to include defender actions and responses, as well as representations of normal user activities on the computer system being attacked. Cyberattack component models stored in a repository are selected and composed into complete models of target computer systems. Metadata associated with each model guides the selection and composition. The cyberattack models are executed to simulate cyberattacks. Multiple simulation iterations are used to train reinforcement learning algorithms that automatically learn improved attacker or defender strategies.

Hydrological Connectivity Response of Typical Soil and Water Conservation Measures Based on SIMulated Water Erosion Model: A Case Study of Tongshuang Watershed in the Black Soil Region of Northeast China

The black soil region of Northeast China is the largest commercial grain production base in China, accounting for about 25% of the total in China. In this region, the water erosion is prominent, which seriously threatens China’s food security. It is of great significance to effectively identify the erosion-prone points for the prevention and control of soil erosion on the slope of the black soil region in Northeast China. This article takes the Tongshuang small watershed (Heilongjiang Province in China) as an example, which is dominated by hilly landforms with mainly black soil and terraces planted with corn and soybeans. Based on the 2.5 cm resolution Digital Elevation Model (DEM) reconstructed by unmanned aerial vehicles (UAVs), we explore the optimal resolution for hydrological simulation research on sloping farmland in the black soil region of Northeast China and explore the critical water depth at which erosion damage occurs in ridges on this basis. The results show that the following: (1) Compared with the 2 m resolution DEM, the interpretation accuracy of field roads, wasteland, damaged points, ridges and cultivated land at the 0.2 m resolution is increased by 4.55–27.94%, which is the best resolution in the study region. (2) When the water depth is between 0.335 and 0.359 m, there is a potential erosion risk of ridges. When the average water depth per unit length is between 0.0040 and 0.0045, the ridge is in the critical range for its breaking, and when the average water depth per unit length is less than the critical range, ridge erosion damage occurs. (3) When local erosion damage occurs, the connectivity will change abruptly, and the remarkable change in the index of connectivity (IC) can provide a reference for predicting erosion damage.

Numerical Simulation and Experimental Research of Spatial Volume Distribution and Velocity Changes of Abrasive Suspension Jet Based on LES + VOF + DEM

The machining accuracy and cutting ability of abrasive suspension jet (ASJ) are affected by the spatial volume distribution and velocity changes of jet flow and abrasives in the free jet area outside the nozzle. In this paper, large eddy simulation (LES) flow model, volume of fluid (VOF) multiphase flow model and discrete element model (DEM) were used to simulate the flow field characteristics of ASJ. The spatial volume distribution and velocity changes of jet flow and abrasives under 10[Formula: see text]MPa jet pressure were studied and the model validation experiments were carried out. In the comparison between the numerical simulation results and the experimental results, it was found that the average error of the percentage of central jet flow is 3.4%, the percentage of central abrasive particles is 3.0% and the velocity of abrasive particles is 4.2%. The physical structure of the numerical model is fixed (so it does not support abrasive injection jet (AIJ)), and the material physical parameters, initial conditions and boundary conditions are editable. Therefore, it can be used to predict the degree of water divergence, abrasives divergence and abrasives velocity attenuation of ASJ with different jet pressure and abrasive types. This numerical model can provide guidance for the parameter setting of ASJ machining process in order to obtain better machining accuracy and cutting ability.

Study on organic carbon in medium-low mature shale oil reservoirs rich in type II kerogen: Oxidation behavior, thermal effect, and kinetics

In-situ combustion technology (ISC) was a potential technology for efficient development of medium–low maturity shale oil reservoirs. Clarifying the exothermic behavior of organic carbon oxidation will be beneficial for better control of ISC processes. In this work, the organic carbon composition, type II kerogen and residual oil, were separated and extracted from the shale samples of typical medium–low maturity shale oil reservoirs in the Ordos Basin, China. The synchronous thermal analyzer and mass spectrometry (STA-MS), Fourier transform infrared spectroscopy (FTIR), and ROCK-EVAL VI were used to study the products, heat release, and functional group changes of type II kerogen, residual oil, and shale during the oxidation process. The results indicated that the oxidation heat release of shale mainly came from the oxidation reaction between its residual oil and type II kerogen, but the oxidation heat release range of shale was relatively lagging compared to the single residual oil or type II kerogen affected by rock debris. Secondly, under atmospheric pressure, the type II kerogen experiences obviously higher heat release than residual oil and heavy oil, indicating its superior potential of in-situ heat release. In addition, the major productions that CO, H2O, SO2, and CO2 were quantitatively characterized, with CO2 having the highest production of 882.25 mg/g. Furthermore, combined with the molecular weight, an oxidation reaction equation of type II kerogen was constructed. Finally, the oxidation kinetics parameters of type II kerogen, residual oil, and shale were determined using Friedman and Ozawa-Flynn-Wall (OFW) models. It was found that the activation energy of shale was higher than II kerogen and residual in low temperature oxidation (LTO) and fuel deposition (FD) stages. This may be due to the joint reaction between kerogen and residual oil, thus exacerbating the challenge of shale in-situ ignition. However, due to the catalytic effect of rock debris, the activation energy of shale was lower than LTO in high temperature oxidation (HTO) stage, indicating that once fuel deposition was completed, organic carbon will undergo stable combustion. This study has theoretical guidance for the numerical simulation research of ISC development technology of maturity-low shale oil reservoirs.

Effective Patient Scheduling For Public Health Events Based On Deep Reinforcement Learning

Amidst the current worldwide pandemic of the novel coronavirus, there are still issues about the efficient allocation of emergency supplies and shortcomings in recovery protocols. To enhance individuals' health and overall welfare, it is imperative to enhance and refine the emergency patient scheduling framework for health emergencies while upholding the idea of a society with an inherent future for humanity. The paper proposed an Efficient Patient Scheduling for Public Health Emergencies (EPS-PHE) using Multiagent Reinforcement Learning (MRL). Originally intended to solve the PHE and improve EPS, the technology was conceived as a distributed multiagent system. MRL has proposed a decentralized incentive system evaluation technique for the EPS-PHE input. An incentive type called a Direct Incentive (DI) is used to show the relationships between individuals or details about two closely connected events. The quality of data transmission in EPS-PHE or the efficacy of the emergency network is measured using a Collective Incentive (CI). Simulation research has shown the value of the Deep Learning (DL) - MRL approach in EPS-PHE framework optimization.

Numerical simulation and experimental study on nonlinear ultrasonic characterization of graphite size in nodular cast iron

Graphite morphology has a significant influence on mechanical properties such as tensile strength and hardness in nodular cast iron (NCI). Compared to the detection of graphite nodularity, there is relatively less ultrasonic simulation and detection experiment research on the graphite size. Therefore, it is of great importance to study rapid nondestructive testing methods for the internal graphite size of NCI. This study establishes a simulation model of nonlinear ultrasonic penetration longitudinal wave detection for graphite size. The acoustic nonlinearity parameter (ANP) is used to characterize the graphite size. Nonlinear ultrasonic detection experiments on the internal graphite size are conducted to verify the reliability of the nonlinear ultrasonic simulation model. The simulation and experimental results show that the ANP of penetrating longitudinal waves decreases with the increase of average graphite diameter. The relationship between ANP and internal graphite size is established. Moreover, the experimental results verified the accuracy of the numerical model. The decrease in ANP may be related to the decrease in the number of grain boundaries. Therefore, the nonlinear ultrasonic technique is an effective method for characterizing the internal graphite size. The establishment of a nonlinear ultrasonic simulation model for graphite size in NCI lays the research foundation for nonlinear ultrasonic microstructure characterization experiments.

Simulation of ambient temperature and humidity distribution in eight-span greenhouse under different wind conditions and corn height

The greenhouse environment suitable for crop growth usually depends on reasonable ventilation, and outdoor air flow state and indoor crop height are important factors affecting the ventilation effect. This study conducted simulation and experimental research on corn cultivation in an eight-span plastic greenhouse, analyzing the effects of outdoor wind speed, wind direction, and indoor crop height variations on the uniformity of temperature and relative humidity distribution as well as the air flow dynamics within the greenhouse. The study shows that as outdoor wind speed increases from 0.5 m/s to 2.5 m/s, for every 1 m/s increment, temperatures decrease by 2.2 °C and 0.9 °C, respectively. The increase in outdoor wind speed resulted in uniformity of greenhouse temperature and humidity along the vertical height, but exacerbated the degree of inhomogeneity in the east-west direction. The change in outdoor wind direction leads to differences in airflow direction and ventilation volume. On the south side of the greenhouse, while temperature is lower, the relative humidity is higher, with maximum differences in temperature and humidity between the north and south sides being 0.6 °C and 1.5 %, respectively, when the northeast wind blows. The east wind has the best cooling effect. The uneven degree of temperature and humidity distribution in the east and west direction of the greenhouse is positively correlated with crop height. This study can provide valuable insights for adjusting and optimizing natural ventilation strategies in large multi-span greenhouses.

Design and conductivity optimization of La0.3Sr0.7TiO3/La0.8Sr0.2MnO3 bi-layer structures as tubular segmented-in-series solid oxide fuel cell interconnect☆

This study focused on meeting the stringent stability requirements of tubular segmented-in-series solid oxide fuel cells (SOFCs) in reducing and oxidizing atmospheres. To address this challenge, a bi-layer perovskite ceramic interconnect was designed by controlling the oxygen partial pressure, because of the strong correlation between the conductivity of strontium-doped lanthanum titanate (LST) and the oxygen partial pressure. The LST powder was prepared using solid-phase and sol-gel methods, and their influence on particle size and sintering behavior was compared. LST/lanthanum strontium manganite (LSM) bi-layer ceramic interconnects with varying thicknesses were fabricated through screen printing and co-sintering. The results demonstrate favorable interfacial bonding and excellent chemical compatibility between the ceramic layers. The conductivity of the bi-layer interconnect exhibits a temperature-dependent behavior, peaking at 550 °C. Simulation calculations and research findings validate that the conductivity of the bi-layer interconnect is determined by the thickness of the LSM layer and the oxygen partial pressure at the interconnect interface. Optimal conductivity is achieved with a bi-layer interconnect consisting of approximately 15 μm of LST and 4 μm of LSM. This can be attributed to the efficient regulation of oxygen partial pressure at the interface, effectively mitigating LSM decomposition caused by low oxygen partial pressure and the subsequent reduction in conductivity. These results provide valuable fundamental data and methodology for the development of high-performance interconnects for tubular segmented-in-series SOFCs.

Machine learning and numerical simulation research on specific energy consumption for gradated coarse particle two-phase flow in inclined pipes

In deep-sea mining engineering, accurately predicting the energy required per unit length of pipeline to transport a unit mass of solids (dimensionless specific energy consumption, DSEC) is crucial for ensuring energy conservation and efficiency in the project. Based on our previous work, we utilized the machine learning (ML) and the computational fluid dynamics (CFD)–discrete element method (DEM) method to study the transport characteristics and flow field variations of gradated coarse particles in inclined pipes (gradated particles refer to solid particles mixed in specific size and quantity ratios). First, we collect 1185 sets of data from 13 experimental literature, and after analyzing and processing them, an ensemble model based on four other ML models is developed. Both for pure substance particles (PS) and mixed particles (MP), the prediction accuracy of this ensemble model is relatively higher (PSs are spherical particles with uniform size and density, and MPs are particles with different shapes, sizes, and densities). Then, the CFD-DEM process and the operating conditions include low flow velocity with low volume concentration (2 m/s and 2.5%), low flow velocity with high volume concentration (2 m/s and 7.5%), and high flow velocity with low volume concentration (4 m/s and 2.5%). Under conditions of low flow velocity and low concentrations, as well as high flow velocity and low concentrations, the DSEC hardly changes with the variation of the pipe inclination angle. Under low flow velocity and high-concentration conditions, as the pipe gradually becomes vertical, the value of DSEC gradually increases.

Mechanical characteristics and constitutive model of cemented tailings backfill under temperature-time effects

The increase in original rock temperature during the deep metal mine backfill mining process significantly impacts the mechanical properties of cemented tailings backfill (CTB). Uniaxial compression experiments were conducted under various curing temperatures and times, analyzing their mechanical characteristics based on strain energy. The experimental results indicate that the mechanical properties of CTB are significantly affected by temperature. The uniaxial compressive strength (UCS) follows a quadratic trend, initially increasing and then decreasing as the curing temperature changes. At the curing temperature of 40 °C, the hydration reaction of CTB develops fully, the hydration products have strong compactness, and the energy absorbed when failure occurs is the largest. The CTB is divided into a solid framework and pore sections to develop a constitutive damage model using a combined macro and micro experimental approach. Polynomial functions represent the correlation between the damage model parameters and curing time and temperature, resulting in a unified temperature-time damage constitutive model. Using the secondary development platform in FLAC3D finite difference software, the damage constitutive model of the CTB was further developed and integrated into the software, enabling both its development and application. This customized model provides a more comprehensive representation of the overall deformation and failure mechanisms of the CTB during compression. The study can predict the strength of CTB under different curing conditions, and provide theoretical reference for the simulation research of filling mining in mines.