Galerkin-based quasi-smooth manifold element(QSME)method for anisotropic heat conduction problems in composites with complex geometry

The accurate and efficient analysis of anisotropic heat conduction problems in complex composites is crucial for structural design and performance evaluation. Traditional numerical methods, such as the finite element method(FEM), often face a trade-off between calculation accuracy and efficiency. In this paper, we propose a quasi-smooth manifold element(QSME) method to address this challenge, and provide the accurate and efficient analysis of two-dimensional(2D) anisotropic heat conduction problems in composites with complex geometry. The QSME approach achieves high calculation precision by a high-order local approximation that ensures the first-order derivative continuity.The results demonstrate that the QSME method is robust and stable, offering both high accuracy and efficiency in the heat conduction analysis. With the same degrees of freedom(DOFs), the QSME method can achieve at least an order of magnitude higher calculation accuracy than the traditional FEM. Additionally, under the same level of calculation error, the QSME method requires 10 times fewer DOFs than the traditional FEM. The versatility of the proposed QSME method extends beyond anisotropic heat conduction problems in complex composites. The proposed QSME method can also be applied to other problems, including fluid flows, mechanical analyses, and other multi-field coupled problems, providing accurate and efficient numerical simulations.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

Synthesis strategies of covalent organic frameworks: An overview from nonconventional heating methods and reaction media

Covalent organic frameworks(COFs), as an emerging class of porous crystalline materials constructed by covalent links between the building monomers, have gained tremendous attention. Over the past 15 years, COFs have made rapid progress and substantial development in the chemistry and materials fields. However, the synthesis of COFs has been dominated by solvothermal methods for a long time and it usually involves high temperature, high pressure and toxic organic solvents, which created many challenges for environmental considerations. Recently,the exploration of new approaches for facile fabrication of COFs has aroused extensive interest. Hence, in this review, we comprehensively describe the synthetic strategies of COFs from the aspects of nonconventional heating methods and reaction media. In addition, the advantages,limitations and properties of the preparation methods are compared. Finally, we outline the main challenges and development prospects of the synthesis of COFs in the future and propose some possible solutions.

A particle-resolved heat-particle-fluid coupling model by DEM-IMB-LBM

Multifield coupling is frequently encountered and also an active area of research in geotechnical engineering.In this work,a particle-resolved direct numerical simulation(PR-DNS)technique is extended to simulate particle-fluid interaction problems involving heat transfer at the grain level.In this extended technique,an immersed moving boundary(IMB)scheme is used to couple the discrete element method(DEM)and lattice Boltzmann method(LBM),while a recently proposed Dirichlet-type thermal boundary condition is also adapted to account for heat transfer between fluid phase and solid particles.The resulting DEM-IBM-LBM model is robust to simulate moving curved boundaries with constant temperature in thermal flows.To facilitate the understanding and implementation of this coupled model for non-isothermal problems,a complete list is given for the conversion of relevant physical variables to lattice units.Then,benchmark tests,including a single-particle sedimentation and a two-particle drafting-kissing-tumbling(DKT)simulation with heat transfer,are carried out to validate the accuracy of our coupled technique.To further investigate the role of heat transfer in particle-laden flows,two multiple-particle problems with heat transfer are performed.Numerical examples demonstrate that the proposed coupling model is a promising high-resolution approach for simulating the heat-particle-fluid coupling at the grain level.

A Novel Localized Meshless Method for Solving Transient Heat Conduction Problems in Complicated Domains

This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.

Numerical manifold method for thermo-mechanical coupling simulation of fractured rock mass

As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.

Research on Operation Optimization of Heating System Based on Electric Storage Coupled Solar Energy and Air Source Heat Pump

For heating systems based on electricity storage coupled with solar energy and an air source heat pump(ECSA),choosing the appropriate combination of heat sources according to local conditions is the key to improving economic efficiency.In this paper,four cities in three climatic regions in China were selected,namely Nanjing in the hot summer and cold winter region,Tianjin in the cold region,Shenyang and Harbin in the severe cold winter region.The levelized cost of heat(LCOH)was used as the economic evaluation index,and the energy consumption and emissions of different pollutants were analyzed.TRNSYS software was used to simulate and analyze the system performance.The Hooke-Jeeves optimization algorithm and GenOpt software were used to optimize the system parameters.The results showed that ECSA systemhad an excellent operation effect in cold region and hot summer and cold winter region.Compared with ECS system,the systemenergy consumption,and the emission of different pollutants of ECSA system can be reduced by a maximum of 1.37 times.In cold region,the initial investment in an air source heat pump is higher due to the lower ambient temperature,resulting in an increase in the LOCH value of ECSA system.After the LOCH value of ECSA system in each region was optimized,the heating cost of the system was reduced,but also resulted in an increase in energy consumption and the emission of different pollutant gases.

Clinical Study on Application of Xinrun Tongluo Method Based on the Theory of Collateral Diseases in Treating Androgenic Baldness of Blood Heat Wind Dryness Syndrome

[Objectives]To explore the intervention effect of the representative formula of Xinrun Tongluo method,Liangxue Xiaofeng Powder,on the incidence of androgenic alopecia in the syndrome of blood heat and wind dryness.[Methods]A total of 72 patients with androgenic alopecia in Suzhou TCM Hospital Affiliated to Nanjing University of Chinese Medicine from October,2022 to June,2023 were randomly divided into a control group(36 cases,treated with Western medicine)and a treatment group(36 cases,treated with Chinese herbal formula+Western medicine).The short-term and long-term efficacy of the two groups of patients was compared.[Results]The hair microscopic signs and short-term and long-term efficacy of the treatment group were significantly better than those of the control group before and after treatment,with a statistically significant difference(P<0.05).[Conclusions]The representative formula of Xinrun Tongluo method is Liangxue Xiaofeng Powder,which has better clinical efficacy as an auxiliary Western medicine in the treatment of androgenic alopecia patients with blood heat and wind dryness syndrome,and is worthy of further promotion and application in clinical practice.

CALCULATING METHOD OF AERODYNAMIC HEATING FOR HYPERSONIC AIRCRAFTS

A new calculating method of aerodynamic heating for unsteady hypersonic aircrafts with complex configuration is presented.This method,which considers the effects of high temperature chemical non-equilibrium and the heat transfer process in thermal protection structure,is based on the combination of the inviscid outerflow solution and the engineering method,where the Euler solver provides the flow parameters on boundary layer edge for engineering method in aerodynamic heating calculation.A high efficient interpolation technique,which can be applied to the fast computation of longtime aerodynamic heating for hypersonic aircraft,is developed for flying trajectory.In this paper,three hypersonic test cases are calculated,and the heat flux and temperature distribution of thermo-protection system are shown.The numerical results show the high efficiency of the developed method and the validation of thermal characteristics analysis on hypersonic aerodynamic heating.

Chemical Oxygen Demand of Seawater Determined with a Microwave Heating Method

This paper investigates a microwave heating method for the determination of chemical oxygen demand (COD) in seawater. The influences of microwave-power, heating time and standard substances on the results are studied. Using the proposed method, we analyzed the glucose standard solution, the coefficient of variation being less than 2%. Compared with the traditional electric stove heating method, the results of F-test and T-test showed that there was no significant difference between the two methods, but the microwave method had slightly higher precision and reproducibility than the electric stove method. With the microwave heating method, several seawater samples from Jiaozhou Bay and the South Yellow Sea were also analyzed. The recovery was between 97.5% and 104.3%. This new method has the advantages of shortening the heating time, improving the working efficiency and having simple operation and therefore can be used to analyze the COD in seawater.

Numerical Study by Imposing the Finite Difference Method for Unsteady Casson Fluid Flow with Heat Flux

This article presents an investigation into the flow and heat transfer characteristics of an impermeable stretching sheet subjected to Magnetohydrodynamic Casson fluid. The study considers the influence of slip velocity, thermal radiation conditions, and heat flux. The investigation is conducted employing a robust numerical method that accounts for the impact of thermal radiation. This category of fluid is apt for characterizing the movement of blood within an industrial artery, where the flow can be regulated by a material designed to manage it. The resolution of the ensuing system of ordinary differential equations (ODEs), representing the described problem, is accomplished through the application of the finite difference method. The examination of flow and heat transfer characteristics, including aspects such as unsteadiness, radiation parameter, slip velocity, Casson parameter, and Prandtl number, is explored and visually presented through tables and graphs to illustrate their impact. On the stretching sheet, calculations, and descriptions of the local skin-friction coefficient and the local Nusselt number are conducted. In conclusion, the findings indicate that the proposed method serves as a straightforward and efficient tool for exploring the solutions of fluid models of this kind.

The Regional Surface Heating Field over the Heterogeneous Landscape of the Tibetan Plateau Using MODIS and In Situ Data

In this study, a parameterization scheme based on Moderate Resolution Imaging Spectroradiometer （MODIS） data and in-situ data was tested for deriving the regional surface heating field over a heterogeneous landscape. As a case study, the methodology was applied to the whole Tibetan Plateau （TP） area. Four images of MODIS data （i.e., 30 January 2007, 15 April 2007, 1 August 2007, and 25 October 2007） were used in this study for comparison among winter, spring, summer, and autumn. The results were validated using the observations measured at the stations of the Tibetan Observation and Research Platform （TORP）. The results show the following： （1） The derived surface heating field for the TP area was in good accord with the land-surface status, showing a wide range of values due to the strong contrast of surface features in the area. （2） The derived surface heating field for the TP was very close to the field measurements （observations）. The APD （absolute percent difference） between the derived results and the field observations was 〈10%. （3） The mean surface heating field over the TP increased from January to April to August, and decreased in October. Therefore, the reasonable regional distribution of the surface heating field over a heterogeneous landscape can be obtained using this methodology. The limitations and further improvement of this method are also discussed.

A distributed measurement method for in-situ soil moisture content by using carbon-fiber heated cable

Moisture content is a fundamental physical index that quantifies soil property and is closely associatedwith the hydrological, ecological and engineering behaviors of soil. To measure in-situ soil moisturecontents, a distributed measurement system for in-situ soil moisture content （SM-DTS） is introduced.The system is based on carbon-fiber heated cable （CFHC） technology that has been developed to enhancethe measuring accuracy of in-situ soil moisture content. Using CFHC technique, a temperature characteristicvalue （Tt） can be defined from temperatureetime curves. A relationship among Tt, soil thermalimpedance coefficient and soil moisture content is then established in laboratory. The feasibility of theSM-DTS technology to provide distributed measurements of in-situ soil moisture content is verifiedthrough field tests. The research reported herein indicates that the proposed SM-DTS is capable ofmeasuring in-situ soil moisture content over long distances and large areas.

Study on the Pyrolytic Carbon Generated by the Electric Heating CVD Method

A new method of fabricating C/C composite materials, namely electric heating CVD method, was used, which electrified the carbon fiber directly by using the conductivity of itself. Acetylene was used as the carbon source with nitrogen as dilution gas, and the pyrolytic carbon started to deposit on the carbon fiber surface when the deposition temperature was reached. The morphology of pyrolytic carbon was characterized by SEM, and the surface properties of carbon fibers before and after CVD were characterized by Raman spectroscopy. The experimental results show that the electric heating method is a novel method to fabricate C/C composite materials, which can form a dense C/C composite material in a short time. The order degree and the average crystallite size of the carbon fiber surface were decreased after the experiment.

Lunar In Situ Large-Scale Construction:Quantitative Evaluation of Regolith Solidification Techniques

Lunar habitat construction is crucial for successful lunar exploration missions.Due to the limitations of transportation conditions,extensive global research has been conducted on lunar in situ material processing techniques in recent years.The aim of this paper is to provide a comprehensive review,precise classification,and quantitative evaluation of these approaches,focusing specifically on four main approaches:reaction solidification(RS),sintering/melting(SM),bonding solidification(BS),and confinement formation(CF).Eight key indicators have been identified for the construction of low-cost and highperformance systems to assess the feasibility of these methods:in situ material ratio,curing temperature,curing time,implementation conditions,compressive strength,tensile strength,curing dimensions,and environmental adaptability.The scoring thresholds are determined by comparing the construction requirements with the actual capabilities.Among the evaluated methods,regolith bagging has emerged as a promising option due to its high in situ material ratio,low time requirement,lack of hightemperature requirements,and minimal shortcomings,with only the compressive strength falling below the neutral score.The compressive strength still maintains a value of 2–3 MPa.The proposed construction scheme utilizing regolith bags offers numerous advantages,including rapid and large-scale construction,ensured tensile strength,and reduced reliance on equipment and energy.In this study,guidelines for evaluating regolith solidification techniques are provided,and directions for improvement are offered.The proposed lunar habitat design based on regolith bags is a practical reference for future research.

Analytical investigations of in situ stress inversion from borehole breakout geometries

This study aims to investigate the feasibility of deriving in situ horizontal stresses from the breakout width and depth using the analytical method.Twenty-three breakout data with different borehole sizes were collected and three failure criteria were studied.Based on the Kirsch equations,relatively accurate major horizontal stress(sH)estimations from known minor horizontal stress(sh)were achieved with percentage errors ranging from 0.33%to 44.08%using the breakout width.The Mogi-Coulomb failure criterion(average error:13.1%)outperformed modified Wiebols-Cook(average error:19.09%)and modified Lade(average error:18.09%)failure criteria.However,none of the tested constitutive models could yield reasonable sh predictions from known sH using the same approach due to the analytical expression of the redistributed stress and the nature of the constitutive models.In consideration of this issue,the horizontal stress ratio(sH/sh)is suggested as an alternative input,which could estimate both sH and sh with the same level of accuracy.Moreover,the estimation accuracies for both large-scale and laboratory-scale breakouts are comparable,suggesting the applicability of this approach across different breakout sizes.For breakout depth,conformal mapping and complex variable method were used to calculate the stress concentration around the breakout tip,allowing the expression of redistributed stresses using binomials composed of sH and sh.Nevertheless,analysis of the breakout depth stabilisation mechanism indicates that additional parameters are required to utilise normalised breakout depth for stress estimation compared to breakout width.These parameters are challenging to obtain,especially under field conditions,meaning utilising normalised breakout depth analytically in practical applications faces significant challenges and remains infeasible at this stage.Nonetheless,the normalised breakout depth should still be considered a critical input for any empirical and statistical stress estimation method given its significant correlation with horizontal stresses.The outcome of this paper is expected to contribute valuable insights into the breakout stabilisation mechanisms and estimation of in situ stress magnitudes based on borehole breakout geometries.

Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

Impact of a Magnetic Dipole on Heat Transfer in Non-Conducting Magnetic Fluid Flow over a Stretching Cylinder

The thermal behavior of an electrically non-conducting magnetic liquid flowing over a stretching cylinder under the influence of a magnetic dipole is considered.The governing nonlinear differential equations are solved numerically using a finite element approach,which is properly validated through comparison with earlier results available in the literature.The results for the velocity and temperature fields are provided for different values of the Reynolds number,ferromagnetic response number,Prandtl number,and viscous dissipation parameter.The influence of some physical parameters on skin friction and heat transfer on the walls of the cylinder is also investigated.The applicability of this research to heat control in electronic devices is discussed to a certain extent.

Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

Catalytic CVD Growth of Carbon Nanotubes by Electric Heating Method

Carbon nanotubes（CNTs） were synthesized by the electric heating catalytic chemical deposition method（CCVD） using acetylene（C2H2） as the carbon source and nitrogen（N2） as carrier gas,and nickel catalyst was loaded by electroplating.The electric heating method,as a new method,electrifies the carbon fiber directly by using its conductivity.The morphology and structure of CNTs were characterized by SEM and TEM,and the surface properties of carbon fibers before and after the growth of CNT were characterized by Raman spectroscopy.The experimental results show that the electric heating method is a new method to produce CNT,and can grow a large number of CNTs in a short time,the crystallization degree and surface average crystallite size of carbon fiber increased after the growth of CNT on it.In addition,electroplating loading catalyst can also be used as an ideal loading way,which can control the number,shape,and distribution of nickel particles by controlling the plating time.