Computer Simulation Research Articles

Numerical prediction model for long- and short-term concentration of indoor volatile organic compounds from building materials

ABSTRACT Emission models of volatile organic compounds (VOCs) from individual indoor building materials have been developed and validated. However, multiple indoor building materials release VOCs simultaneously, and neither single building material nor multiple building material emission models can predict the entire release cycle of VOCs accurately. This study established a long- and short-term numerical prediction model for indoor VOC concentration. The model includes an attenuation coefficient θ. To describe the decay rate of the total VOC content, which is mainly influenced by time, and by designing experiments and testing in environmental warehouses under different seasonal conditions, the value of θ was first obtained. Then, after successfully plotting the emission curve of indoor pollutant concentration over time through numerical solution and using θ, the VOC content was corrected for various seasonal conditions. On the basis of this model, an exposure dose integration algorithm was proposed to evaluate the environmental health risks, as an application of this model. In comparison with previous research results and experimental data, this model has better predictive performance.

The CARMENES search for exoplanets around M dwarfs. Revisiting the GJ 581 multi-planetary system with new Doppler measurements from CARMENES, HARPS, and HIRES

GJ\,581 is a nearby M dwarf known to host a packed multiple planet system composed of two super-Earths and a Neptune-mass planet. We present new orbital analyses of the GJ\,581 system, utilizing recent radial velocity (RV) data obtained from the CARMENES spectrograph combined with newly reprocessed archival data from the HARPS and HIRES spectrographs. Our aim was to analyze the post-discovery spectroscopic data of GJ\,581, which were obtained with CARMENES. In addition, we used publicly available HIRES and HARPS spectroscopic data to seek evidence of the known and disputed exoplanets in this system. We aimed to investigate the stellar activity of GJ\,581 and update the planetary system's orbital parameters using state-of-the-art numerical models and techniques. We performed a periodogram analysis of the available precise CARMENES, HIRES, and HARPS RVs and of stellar activity indicators. We conducted detailed orbital analyses by testing various orbital configurations consistent with the RV data. We studied the posterior probability distribution of the parameters fit to the data and we explored the long-term stability and overall orbital dynamics of the GJ\,581 system. We refined the orbital parameters of the GJ\,581 system using the most precise and complete set of Doppler data available. Consistent with the existing literature, our analysis confirms that the system is unequivocally composed of only three planets detectable in the present data, dismissing the putative planet GJ\,581\,d as an artifact of stellar activity. Our N-body fit reveals that the system's inclination is $i = $\,deg, which implies that the planets could be up to 30<!PCT!> more massive than their previously reported minimum masses. Furthermore, we report that the GJ\,581 system exhibits long-term stability, as indicated by the posterior probability distribution, characterized by secular dynamical interactions without the involvement of mean motion resonances.

Inverse interaction of shear between steel-stirrups and externally bonded carbon fiber-reinforced polymers in shear-strengthened reinforced concrete beams: Analytical and numerical models

Shear failure in reinforced concrete (RC) beams have always been a serious concern due to its brittle fracture mode. In addition, many questions are raised about the accuracy of current design guidelines for predicting the shear resistance contribution of externally bonded carbon fiber-reinforced polymers (EB-CFRPs) to the ultimate shear strength of strengthened RC beams, particularly with regard to the inverse interaction of the shear contribution between steel-stirrups and EB-CFRPs. The main objective of the present study is to implement experimental and numerical tests to develop analytical and numerical models for RC beams strengthened in shear using EB-CFRPs. Emphasis is placed on the negative inverse interaction between the steel-stirrups and the EB-CFRPs as the ratio of EB-CFRP-to-steel-stirrups increases with increasing the CFRP rigidity (CFRP thickness and CFRP width). The inverse interaction is not included in the design models proposed by most current guidelines, although it has a considerable effect on shear resistance as predicted by the guidelines. In fact, the shear contribution associated with EB-CFRP decreases as the ratio of EB-CFRP-to-steel-stirrups increases. Therefore, proposing reliable effective strains by including these parameters improve the calculated shear contribution of EB-CFRPs. First, an analytical model is proposed for CFRP effective strain considering the inverse interaction between EB-CFRPs and steel-stirrups. Afterward, a validation of the proposed model with experimental data is done by conducting a parametric study of the increasing trends with respect to the ratio of EB-CFRP-to-steel-stirrups. A numerical finite-element model for the reduction factor and the corresponding effective strain based on the inverse interactions between EB-CFRPs and steel-stirrups is also proposed, and the results are compared with various current guidelines. The results are presented in terms of shear crack patterns, load-midspan deflections, shear stresses, strain responses along the fibers, maximum strain profiles for all the CFRPs and specimens, applied shear forces and strains for all the steel-stirrups and EB-CFRPs, and interactions between steel-stirrups and EB-CFRPs based on their maximum strain contributions at the maximum shear forces and the maximum strain they experience after shear failure.

Real-time position and pose prediction for a self-propelled undulatory swimmer in 3D space with artificial lateral line system

This study aims to investigate the feasibility of using an artificial lateral line (ALL) system for predicting the real-time position and pose of an undulating swimmer with Carangiform swimming patterns. We established a 3D computational fluid dynamics simulation to replicate the swimming dynamics of a freely swimming mackerel under various motion parameters, calculating the corresponding pressure fields. Using the simulated lateral line data, we trained an artificial neural network to predict the centroid coordinates and orientation of the swimmer. A comprehensive analysis was further conducted to explore the impact of sensor quantity, distribution, noise amplitude and sampling intervals of the ALL array on predicting performance. Additionally, to quantitatively assess the reliability of the localization network, we trained another neural network to evaluate error magnitudes for different input signals. These findings provide valuable insights for guiding future research on mutual sensing and schooling in underwater robotic fish.

Aortic Valve Stenosis Causes Accumulation of Extracellular Hemoglobin and Systemic Endothelial Dysfunction.

Whether aortic valve stenosis (AS) can adversely affect systemic endothelial function independently of standard modifiable cardiovascular risk factors is unknown. We therefore investigated endothelial and cardiac function in an experimental model of AS mice devoid of standard modifiable cardiovascular risk factors and human cohorts with AS scheduled for transcatheter aortic valve replacement. Endothelial function was determined by flow-mediated dilation using ultrasound. Extracellular hemoglobin (eHb) concentrations and NO consumption were determined in blood plasma of mice and humans by ELISA and chemiluminescence. This was complemented by measurements of aortic blood flow using 4-dimensional flow acquisition by magnetic resonance imaging and computational fluid dynamics simulations. The effects of plasma and red blood cell (RBC) suspensions on vascular function were determined in transfer experiments in a murine vasorelaxation bioassay system. In mice, the induction of AS caused systemic endothelial dysfunction. In the presence of normal systolic left ventricular function and mild hypertrophy, the increase in the transvalvular gradient was associated with elevated eryptosis, increased eHb and plasma NO consumption; eHb sequestration by haptoglobin restored endothelial function. Because the aortic valve orifice area in patients with AS decreased, postvalvular mechanical stress in the central ascending aorta increased. This was associated with elevated eHb, circulating RBC-derived microvesicles, eryptotic cells, lower haptoglobin levels without clinically relevant anemia, and consecutive endothelial dysfunction. Transfer experiments demonstrated that reduction of eHb by treatment with haptoglobin or elimination of fluid dynamic stress by transcatheter aortic valve replacement restored endothelial function. In patients with AS and subclinical RBC fragmentation, the remaining circulating RBCs before and after transcatheter aortic valve replacement exhibited intact membrane function, deformability, and resistance to osmotic and hypoxic stress. AS increases postvalvular swirling blood flow in the central ascending aorta, triggering RBC fragmentation with the accumulation of hemoglobin in the plasma. This increases NO consumption in blood, thereby limiting vascular NO bioavailability. Thus, AS itself promotes systemic endothelial dysfunction independent of other established risk factors. Transcatheter aortic valve replacement is capable of limiting NO scavenging and rescuing endothelial function by realigning postvalvular blood flow to near physiological patterns. URL: https://www.clinicaltrials.gov; Unique identifier: NCT05603520. URL: https://www.clinicaltrials.gov; Unique identifier: NCT01805739.

Deconvoluting the Individual Effects of Nanoparticle Proximity and Size in Thermocatalysis.

Nanoparticle (NP) size and proximity are two physical descriptors applicable to practically all NP-supported catalysts. However, with conventional catalyst design, independent variation of these descriptors to investigate their individual effects on thermocatalysis remains challenging. Using a raspberry-colloid-templating approach, we synthesized a well-defined catalyst series comprising Pd12Au88 alloy NPs of three distinct sizes and at two different interparticle distances. We show that NP size and interparticle distance independently control activity and selectivity, respectively, in the hydrogenation of benzaldehyde to benzyl alcohol and toluene. Surface-sensitive spectroscopic analysis indicates that the surfaces of smaller NPs expose a greater fraction of reactive Pd dimers, compared to inactive Pd single atoms, thereby increasing intrinsic catalytic activity. Computational simulations reveal how a larger interparticle distance improves catalytic selectivity by diminishing the local benzyl alcohol concentration profile between NPs, thus suppressing its readsorption and consequently, undesired formation of toluene. Accordingly, benzyl alcohol yield is maximized using catalysts with smaller NPs separated by larger interparticle distances, overcoming activity-selectivity trade-offs. This work exemplifies the high suitability of the modular raspberry-colloid-templating method as a model catalyst platform to isolate individual descriptors and establish clear structure-property relationships, thereby bridging the materials gap between surface science and technical catalysts.

Numerical Investigation of Funicular Liquid Bridges between Three Spherical Grains in a Bidisperse Particulate System.

Appropriate capillary effects are beneficial for controlling the wet powder performance and agglomerate formation. As water content rises, the funicular regime supplants the pendular regime as the predominant state in wet granular media. The displacement of grains leading to the stretching of funicular liquid bridges until rupture is an interesting and common phenomenon. Utilizing Surface Evolver software (an energy minimization approach), this work develops an efficient and accurate numerical model to describe liquid interactions among three spherical grains. The effects of liquid volume, contact angle, grain size ratio, grain-pair gap, and separation distance on the capillary forces and rupture distances are investigated. Notably, we present a modified closed-form equation for predicting the rupture distance of funicular bridges between three grains, which reflects the coupled effects of the contact angle, grain size, and liquid volume on rupture distance. This present study provides insights for incorporating capillary effects into mechanical models relying on microassembly composed of several grains in bidisperse particulate systems. Additionally, the numerical findings confirm some findings regarding the splitting of funicular bridges.

AeoLiS: Numerical modelling of coastal dunes and aeolian landform development for real-world applications

The formation and evolution of coastal dunes result from a complex interplay of eco-morphodynamic processes. State-of-the-art models can simulate aeolian transports and morphological dune evolution under certain conditions. However, a model combining these processes for coastal engineering applications was not yet available. This study aims to develop a predictive tool for dune development to inform coastal management decisions and interventions. The aeolian sediment transport model AeoLiS is extended with functionalities that allow for simulations of coastal landforms. The added functionalities include the effect of topographic steering on wind shear, avalanching of steep slopes and vegetation processes in the form of growth and wind shear reduction. The model is validated by simulating four distinct coastal landforms; barchan-, parabolic-, embryo dunes and blowouts. Simulations, based on real-world conditions, replicate the landform formation, migration rates and seasonal variability.

Reliability analysis of shallow foundation using the important-region-based method

ABSTRACT To improve the computational efficiency and accuracy is a constant pursuit for evaluating the stability of geotechnical structures built in spatially variable soils. This study proposes an important-region-based Karhunen – Loève expansion and first-order reliability method (KL-FORM-IR method) to perform reliability analysis for shallow foundation resting on spatially variable soil. The KL-FORM-IR method further improves the computational efficiency by ignoring the variability of the soil parameters outside the important region. An identification method for locating the important region of shallow foundation is first proposed, and the effect of the scale factor (K f), autocorrelation distance of soil properties, and foundation width (B) on the identification of important regions are investigated. For the purpose of validation, the proposed KL-FORM-IR method compared with the KL-FORM and KL-MCSM by performing reliability analyses based on the same numerical model. Conclusively, the KL-FORM-IR method can improve the computational efficiency and it contributes a new approach to the reliability analysis of shallow foundation resting on a spatially variable soil.

Multihazard resilience and economic loss evaluation method for cable-stayed bridges under the combined effects of scour and earthquakes

Cable-stayed bridges (CSBs) commonly experience scour-induced foundation issues during their service life, impacting both structural seismic resilience and economic losses. Nonetheless, these aspects remain inadequately addressed in existing seismic codes. This paper presents a comprehensive probabilistic multihazard fragility, resilience, and economic loss evaluation method for CSBs influenced by scour and earthquakes. Numerical models with various scour depths considering epistemic uncertainties in loads, structural properties, and soil conditions are first established using the Latin-hypercube sampling (LHS) method. Incremental dynamic analyses (IDAs) are then conducted on these models in cases with near-fault and far-field bidirectional ground motions, yielding probabilistic multihazard demand models and component- and system-level multihazard fragility surfaces for diverse limit states. Subsequently, resilience and economic losses are comparatively evaluated by integrating damage probabilities with corresponding uncertain functional recovery functions. Eventually, a rapid bridge traffic capacity evaluation method based on the critical traffic capacity surface concept is devised to effectively assess CSB open patterns under various peak ground accelerations (PGAs), scour depths, and repair times. The numerical results show apparent shifts in the nonmonotonic seismic fragility, resilience, and economic loss of the case CSB due to scour effects, with the rates of change reaching 31.4 %, 34 %, and 24.8 %, respectively, associated with a critical scour depth of 12 m at IM = 0.8 g. The developed analysis framework provides a reference for the life-cycle design and maintenance of CSBs affected by multihazard scour and earthquakes.

Impact of High-Resolution Land Cover on Simulation of a Warm-Sector Torrential Rainfall Event in Guangzhou

This study on the warm-sector heavy rainfall event in Guangzhou on 7 May 2017, examined the effects and mechanisms of incorporating 30 m high-resolution land surface data into its numerical simulation. The updated 1km numerical model, integrating 30 m high-resolution land surface data, successfully captured the initiation, back-building, and organized development of warm-sector convections in Huadu and Zengcheng District. The analysis revealed that the high spatial resolution of the surface data led to a reduced urban area footprint (urban −6.31%), increased vegetation cover (forest 11.63%, croplands 1%), and enhanced surface runoff (water 2.77%) compared with a model’s default land cover (900 m). These changes mitigated the urban heat island (UHI) effect within the metropolitan area and decreased the surface sensible heat flux. This reduction contributed to a pronounced temperature gradient between Huadu Mountain and the urban area. Additionally, a stronger high-pressure recirculation and sea–land breezes facilitated the transport of warm and moist air from the sea inland, creating a humid corridor along the sea–land interface. The consistent influx of warm and moist air near the mountain front, where strong temperature gradients were present, forcibly triggered warm-sector convection, intensifying its organization. This study highlighted the critical role of high-resolution land surface data in the accurate numerical simulation of warm-sector heavy rainfall.

Revealing short-term dynamics of tropical cyclone wind speeds from satellite synthetic aperture radar

Both unresolved physics in numerical models and limited theoretical understanding of the small-scale diffusion processes occurring near the ocean surface hamper predictability of tropical cyclone (TC) wind changes. An analytical model is here developed to diagnose the short-term evolution of the TC wind profile. An effective frictional parameter is introduced to control the unknown diffusion effects. When this frictional parameter is adjusted to match the TC intensity change, solutions are consistent with observed high-resolution ocean surface wind speeds from satellite synthetic aperture radar (SAR). The initial high-resolution estimate of the near-core wind structure is then found to strongly modulate the wind profile evolution. The frictional parameter can, unfortunately, not efficiently be calibrated using outer-core wind speed changes. Low-resolution observations or standard numerical weather predictions may thus not be directly used to reinterpret and anticipate short-term TC wind changes. The expected accumulation of orbiting SAR sensors as well as improved measurements of the ocean-atmosphere boundary layer characteristics shall then become essential to more precisely monitor TC dynamics.

Synthesis and Molecular Dynamic Simulation of Novel Cationic and Non-cationic Pyrimidine Derivatives as Potential G-quadruplex-ligands.

Drug resistance has been a problem in cancer chemotherapy, which often causes shortterm effectiveness. Further, the literature indicates that telomere G-quadruplex could be a promising anti-cancer target. We synthesized and characterized two new pyrimidine derivatives as ligands for G-quadruplex DNA. The interaction of novel non-cationic and cationic pyrimidine derivatives (3a, b) with G-quadruplex DNA (1k8p and 3qsc) was explored by circular dichroism (CD) and ultraviolet-visible spectroscopy and polyacrylamide gel electrophoresis (PAGE) methods. The antiproliferative activity of desired compounds was evaluated by the MTT assay. Apoptosis induction was assessed by Propidium iodide (P.I.) staining and flow cytometry. Computational molecular modeling (CMM) and molecular dynamics simulation (MD) were studied on the complexes of 1k8p and 3qsc with the compounds. The van der Waals, electrostatic, polar solvation, solventaccessible surface area (SASA), and binding energies were calculated and analyzed. The experimental results confirmed that both compounds 3a and 3b interacted with 1k8p and 3qsc and exerted cytotoxic and proapoptotic effects on cancer cells. The number of hydrogen bonds and the RMSD values increased in the presence of the ligands, indicating stronger binding and suggesting increased structural dynamics. The electrostatic contribution to binding energy was higher for the cationic pyrimidine 3b, indicating more negative binding energies. Both experimental and MD results confirmed that 3b was more prone to form a complex with DNA G-quadruplex (1k8p and 3qsc), inhibit cell growth, and induce apoptosis, compared to the non-cationic pyrimidine 3a.

A plate containing a central rectangular hole made of PMMA resin and reinforced with graphene particles: low-velocity impact analysis

Purpose Due to the important and extensive use of graphene in the engineering structures, the purpose of this paper is to investigate the low-velocity impact (LVI) response of a graphene-reinforced plate with a rectangular central hole. Design/methodology/approach In this paper, Halpin-Tsai theory is used to model the mechanical properties of the plate made of polymethyl methacrylate (PMMA) resin and graphene particles. First, the plate displacement field is presented, and then the strains and stresses are obtained. The motion equations are extracted using energy equations, Ritz method and generalized Lagrange equations. The verification of theoretical formulation shall be carried out using the ABAQUS finite element software suite. Findings The effects of graphene volume fraction, approaching the impact point to the corner of the rectangular plate, and different boundary conditions of the plate are studied for LVI of impactor with spherical tip on the rectangular central hole reinforced by graphene. Important responses of the impact, such as the contact force between the impactor and the plate as well as the displacement of the plate at the impact place, are investigated and analyzed in this research. Originality/value Considering the wide application of graphene in the engineering structures, the simulation of the LVI on the graphene-reinforced plate made of PMMA resin with a rectangular central hole is carried out in this paper using numerical and theoretical modeling.

Influence of Mg and Cu on precipitation behaviors and mechanical properties of Al–Si alloys

Al–8Si-0.5 Mg, Al–8Si–2Cu, and Al–8Si–2Cu-0.5 Mg alloys have been designed in this study, aiming to reveal the influence of Mg and/or Cu additions on microstructures and mechanical properties of Al–Si alloys. Transmission electron microscopy (TEM), differential scanning calorimetry (DSC) and the Kampman-Wagner Numerical (KWN) model have been used to investigate age-hardening behaviors as well as their relation to mechanical properties. Microstructure analysis suggest that β'', θ′ and QP precipitates have formed in peak-aged Al–8Si-0.5 Mg, Al–8Si–2Cu and Al–8Si–2Cu-0.5 Mg alloys, respectively. Notably, the Al–8Si–2Cu-0.5 Mg alloy exhibits faster precipitation kinetics, larger age-hardening response and remarkable thermal stability compared to the other two alloys. Theoretical calculations indicate that the nucleation of QP precipitates is easier than β'' and θ′ precipitates, but the slower diffusivity of Cu limits the growth rate of QP precipitates during the growth stage. Furthermore, the contribution of β'', θ′ and QP precipitates to yield strength of peak-aged Al–Si alloys has been quantified, which shows a significant variation. β'' precipitates in the Al–8Si-0.5 Mg alloy exhibit the highest strengthening effect (∼186.4 MPa), followed by QP precipitates in Al–8Si–2Cu-0.5 Mg (∼160.2 MPa) and θ′ precipitates in the Al–8Si–2Cu alloy (∼13.4 MPa). This work provides valuable insights into the influence of Mg and/or Cu additions on the aging response and mechanical properties of Al–Si alloys, offering a reference for further enhancing the performance of Al–Si alloys.

Computer Simulation of Rural Landscape Design Based on Remote Sensing Image Technology

Introduction: The field of rural landscape design deals with the design and recovery of rural areas and landscape in a way that it can support natural biodiversity in addition to human needs in sustainable ways as well as maintaining its cultural character. It employs the use of plants, landforms and some water without interference with the environment in order to come up with useful and beautiful spaces. Recognizing the importance of preserving the environment while at the same time investing in developmental projects, it is centralized and focuses on the entire eco-system with the aim of enhancing the lives of the people. Aim: This study aims to develop theoretical aspect of an innovative computer simulation model for designing rural landscapes by applying the technology of remote sensing image. Research methodology: We suggest a new Starling Murmuration search-driven Adaptive YOLOv7 algorithm to identify and categorize several rural buildings and setting types. For the image data, we collected abundant data from several environments using UAV devices to train our proposed model. It is not surprising that our proposed model combined the use of three dimensional (3D) geographic information system (GIS) virtual imaging design model in the simulation of the rural landscape designs. Our recommended model is then extended using SM optimization to improve object detection with YOLOv7. By repeated adjustments of the network parameters in a somewhat similar fashion like flocking, we managed to enhance both accuracy and efficiency. This framework exploits crowdsourcing for delimiting rural buildings and landscapes with high-fidelity. Findings and Conclusion: We implemented our recommended model in Python software. During the phase of evaluation, we evaluate the efficacy of our recommended SM-AYOLOv7 model across a variety of parameters such as precision (91.72%), recall (92.34%), Intersection over Union (IoU) (90.23%), and f1 score (93.64%). Our experimental results precisely indicate that our approach outperforms traditional approaches. We demonstrate significant increases in accuracy and adaptability, especially when adjusting to dynamic configurations.

Studies on the properties and molecular dynamics simulations of nonionic surfactants based on succinic acid derivatives

AbstractNonionic surfactants have proven useful in various applications such as wastewater treatment, enhanced oil recovery, dyeing, and cosmetics. Novel nonionic surfactants such as PEMP, BEMP, HEMP and BEEP were synthesized based on succinic acid derivatives and using L‐isoleucine as the linking group and four polyether alcohols as the hydrophilic group. First, the structures of the four nonionic surfactants were studied using Fourier transform infrared (FTIR) spectroscopy and nuclear magnetic resonance (NMR) spectra. Then the critical micelle concentration (CMC) and surface tension at CMC (γCMC) of the four nonionic surfactants in aqueous solution were tested. γ‐lg c curves determined the relationships of their surface properties as BEEP > PEMP > BEMP > HEMP. In order to expand the range of applications for nonionic surfactants, we evaluated the salt‐resistant properties of four such surfactants. Our findings demonstrated that this class of surfactants indeed has superior salt‐resistant properties. Molecular dynamics (MD) simulations were used to study how surfactant molecules aggregate in the interfacial film. The study investigated the trend of solvent accessible surface area (SASA) over time. Results showed that the surfactant molecules interacted well with solvent molecules in the equilibrium state. This study investigates the performance differences among four types of surfactants using the electrostatic potential (ESP) distribution of their molecules. The study employed both experimental and computational simulations to provide a more comprehensive understanding of surfactant properties. The results offer insights into the theoretical research and application extension of this class of surfactants.

Time history response prediction of stochastic SDOF structures through Kriging-NARX modeling and adaptive real-time hybrid simulation

Emulating time history response is a critical benchmark for evaluating the seismic resilience of structures and assessing their reliability in the face of uncertainties. Traditional computational simulations, while commonly favored, often fall short in accurately describing the behavior of intricate structural components. Real-time hybrid simulation (RTHS) offers an effective and efficient experimental technique for replicating the responses of deterministic systems in a cyber-physical manner but encounters challenges when dealing with stochastic systems. This study proposes an innovative application of RTHS for time history response prediction of nonlinear stochastic structures. A cross-validation (CV)-Voronoi sampling strategy is applied to sequentially select uncertain parameter samples for RTHS. Nonlinear autoregressive with exogenous input (NARX) model is identified and updated from RTHS results in a sequential manner. Coefficients of final NARX model are surrogated using the Kriging meta-models to account for given uncertainties. Laboratory tests on single-degree-of-freedom (SDOF) structures with self-centering viscous dampers (SC-VD) were conducted to demonstrate the effectiveness of the proposed approach. Predicted time history responses from the final Kriging-NARX model are further compared and validated. It is demonstrated that the proposed approach provides an effective and efficient method for time history response prediction of stochastic systems especially when accurate numerical models are not available.

Mechanical behavior and cement sheath damage characteristics of casing pullout based on peridynamics

In petroleum exploration and production, the recovery and replacement of casing are crucial. However, difficulties often arise during casing pullout operations due to the strong bonding strength between the casing and the cement sheath. Currently, there are few researches on casing pullout process, which leads to the fact that casing pullout operations often relies on experience and lacks theoretical support. To comprehend the mechanical behavior, interface shear failure mechanisms of casing pullout process, and to propose rational solutions for practical engineering applications, this study establishes a Peridynamics (PD) numerical model for casing pullout to conduct research and analysis on its process. Based on the characteristics of the cement sheath, a damage model for the cement sheath is proposed in this paper. Additionally, a kernel function is introduced to reflect the variation of long-range forces with the distance between material points, aiming to characterize the changes in the stiffness of bond constant within the horizon of material points. The mechanical responses of the cement sheath and casing at different phases of the casing pullout process are analyzed and characterized. The research results indicate that an increase in the casing-cement sheath thickness ratio promotes the damage at the casing-cement sheath interface, reducing the time required for casing pullout but inhibiting the expansion of cement sheath cracks. As the thickness ratio increases from 0.1 to 0.18, the degree of cement sheath damage decreases from 6.904% to 6.453%. The increase in Young's modulus of the cement sheath is negatively correlated with the time required for casing pullout, positively correlated with the interface failure rate, and it does not exhibit a linear relationship with the degree of cement sheath damage. Within the range of Young's modulus from 4 to 20 GPa, the maximum difference in cement sheath damage degree is 3.715%. In engineering practice, the Young's modulus of well cement and the thickness ratio of casing-cement sheath can be appropriately increased to enhance pullout efficiency.

Three-Dimensional Discrete Element Analysis of Bearing Characteristics of Concrete–Cored Sand–Gravel Pile Composite Foundation under Cyclic Dynamic Load

Concrete-cored sand–gravel piles are a kind of composite pile formed by wrapping a concrete-cored pile with a sand–gravel shell, which has the advantages of both a rigid pile and bulk-material pile. The bearing characteristics of the concrete-cored sand–gravel pile composite foundation were investigated by establishing a three-dimensional discrete element numerical model for a cyclic dynamic loading test. The results show that the vertical stress of the core pile body fluctuates greatly at the beginning of loading, and the fluctuation amplitude decreases with the depth, and gradually tends to be stable in the middle and late stages, and the vertical-stress distribution is relatively uniform. The radial stress in the upper part of the core pile body fluctuates greatly, the fluctuation in the lower part is small, and the radial stress in each part of the core pile body gradually tends to be stable in the late-loading period. The radial stress factor of the core pile body reaches the stable speed with the foundation depth decreasing; the fluctuation amplitude of the pile-soil stress ratio decreases with the foundation depth and gradually tends to be stable with the increase in loading. The results of this study can provide a reference for the design and construction of a core sand pile composite foundation.