Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

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.

Experimental Study of Heat Transfer in an Insulated Local Heated fromBelow and Comparison with Simulation by Lattice Boltzmann Method

In this paper,experimental and numerical studies of heat transfer in a test local of side H=0.8 m heated from below are presented and compared.All the walls,the rest of the floor and the ceiling are made from plywood and polystyrene in sandwich form(3 mmplywood-3 cm polystyrene-3 mmplywood)just on one of the vertical walls contained a glazed door(2 H/3×0.15 m).This local is heated during two heating cycles by a square plate of iron the width L=0.6 H,which represents the heat source,its temperature Th is controlled.The plate is heated for two cycles by an adjustable set-point heat source placed just down the center of it.For each cycle,the heat source is switched“on”for 6 h and switched“off”for 6 h.The outdoor air temperature is kept constant at a low temperature Tc<Th.All measurements are carried out with k-type thermocouples and with flux meters.Results will be qualitatively presented for two cycles of heating in terms of temperatures and heat flux densitiesϕfor various positions of the test local.The temperature evolution of the center and the profile of the temperature along the vertical centerline are compared by two dimensions simulation using the lattice Boltzmann method.The comparison shows a good agreement with a difference that does not exceed±1℃.

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.

A novel adaptive harmonic balance method with an asymptotic harmonic selection

The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The adaptive harmonic balance(AHB)method is an improved HBM method.This paper presents a modified AHB method with the asymptotic harmonic selection(AHS)procedure.This new harmonic selection procedure selects harmonics from the frequency spectra of nonlinear terms instead of estimating the contribution of each harmonic to the whole nonlinear response,by which the additional calculation is avoided.A modified continuation method is proposed to deal with the variable size of nonlinear algebraic equations at different values of path parameters,and then all solution branches of the amplitude-frequency response are obtained.Numerical experiments are carried out to verify the performance of the AHB-AHS method.Five typical nonlinear dynamic equations with different types of nonlinearities and excitations are chosen as the illustrative examples.Compared with the classical HBM and Runge-Kutta methods,the proposed AHB-AHS method is of higher accuracy and better convergence.The AHB-AHS method proposed in this paper has the potential to investigate the nonlinear vibrations of complex high-dimensional nonlinear systems.

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.

Harmonic Balance Methods:A Review and Recent Developments

The harmonic balance(HB)method is one of the most commonly used methods for solving periodic solutions of both weakly and strongly nonlinear dynamical systems.However,it is confined to low-order approximations due to complex symbolic operations.Many variants have been developed to improve the HB method,among which the time domain HB-like methods are regarded as crucial improvements because of their fast computation and simple derivation.So far,there are two problems remaining to be addressed.i)A dozen of different versions of HB-like methods,in frequency domain or time domain or in hybrid,have been developed;unfortunately,misclassification pervades among them due to the unclear borderlines of different methods.ii)The time domain HB-like methods suffer from non-physical solutions,which have been shown to be caused by aliasing(mixture of the high-order into the low-order harmonics).Although a series of dealiasing techniques have been developed over the past two decades,the mechanism of aliasing and the final solution to dealiasing are still not well known to the academic community.This paper aims to provide a comprehensive review of the development of HB-like methods and enunciate their principal differences.In particular,the time domain methods are emphasized with the famous aliasing phenomenon clearly addressed.

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.

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.

Amplitude and Period Effect on Heat Transfer in an Enclosure with Sinusoidal Heating from Below Using Lattice Boltzmann Method

This work presents a simulation of the phenomena of natural convection in an enclosure with a variable heating regime by the lattice Boltzmann method(LBM).We consider a square enclosure of side H filled with air(Pr=0.71)and heated from below,with a hot portion of length L=0.8 H,by imposing a sinusoidal temperature.The unheated segments of the bottom wall are treated as adiabatic,and one of the vertical walls features a cold region,while the remaining walls remain adiabatic.The outcomes of the two-dimensional(2D)problem are depicted through isotherms,streamlines,the temperature evolution within the enclosure,and the Nusselt number.These visualizations span various amplitude values“a”in the interval[0.2,0.8],and of the period T0 for Ra=107.The amplitude and period effect on the results is evaluated and discussed.The amplitude of the temperature at the heart of the enclosure increases with the increase in amplitude.This also increases with the period(T0)of the imposed temperature,something that is not observable on the global Nusselt number.

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.

Regulating metal-acid double site balance on mesoporous SiO_(2)-Al_(2)O_(3) composite oxide for supercritical n-decane cracking

The balance between metal and acid sites directly affects the preparation of high-performance cracking catalysts with high heat sink and low coking.Nevertheless,how to control acid-metal sites balance and its relationship with cracking performance are reported scarcely.In this work,a series of Pt/Al_(2)O_(3)-SiO_(2) dual sites catalysts with different metal to acid active sites ratio(C_(M)/C_(SA))were constructed by ethanolassisted impregnation method and the impact on n-decane cracking under supercritical conditions was systematically and deeply investigated.The results showed that the conversion and carbon deposition increased gradually with varied C_(M)/C_(SA)and reached the balance at C_(M)/C_(SA)of 0.13.The proper ratio C_(M)/C_(SA)(0.13)can balance the deep dehydrogenation coking over metal active sites and high heat sink of cracking over acid active sites,the chemical heat sink reaches amazing 1.75 MJ/kg and carbon deposition is only22.03 mg/cm^(2) at 750℃.Meanwhile,the few metal sites at low C_(M)/C_(SA)and the few strong acid sites at high C_(M)/C_(SA)are the main factors limiting the cracking activity.Low C_(M)/C_(SA)limit the activation of C-H bond and deep dehydrogenation of coking precursor,resulting in relative low cracking activity and carbon deposition,while high C_(M)/C_(SA)limit the activation of C-C bond and increase the deep dehydrogenation.In this contribution,design and construction of metal-acid dual sites can not only provide the technical solution for the preparation of high heat sink and low coking cracking catalyst,but also deepen the understanding of the cracking path of hydrocarbon fuel.

Machine learning models for the density and heat capacity of ionic liquid-water binary mixtures

Ionic liquids(ILs),because of the advantages of low volatility,good thermal stability,high gas solubility and easy recovery,can be regarded as the green substitute for traditional solvent.However,the high viscosity and synthesis cost limits their application,the hybrid solvent which combining ILs together with others especially water can solve this problem.Compared with the pure IL systems,the study of the ILs-H_(2)O binary system is rare,and the experimental data of corresponding thermodynamic properties(such as density,heat capacity,etc.)are less.Moreover,it is also difficult to obtain all the data through experiments.Therefore,this work establishes a predicted model on ILs-water binary systems based on the group contribution(GC)method.Three different machine learning algorithms(ANN,XGBoost,LightBGM)are applied to fit the density and heat capacity of ILs-water binary systems.And then the three models are compared by two index of MAE and R^(2).The results show that the ANN-GC model has the best prediction effect on the density and heat capacity of ionic liquid-water mixed system.Furthermore,the Shapley additive explanations(SHAP)method is harnessed to scrutinize the significance of each structure and parameter within the ANN-GC model in relation to prediction outcomes.The results reveal that system components(XIL)within the ILs-H_(2)O binary system exert the most substantial influence on density,while for the heat capacity,the substituents on the cation exhibit the greatest impact.This study not only introduces a robust prediction model for the density and heat capacity properties of IL-H_(2)O binary mixtures but also provides insight into the influence of mixture features on its density and heat capacity.

Harmonic balance simulation of the influence of component uniformity and reliability on the performance of a Josephson traveling wave parametric amplifier

A Josephson traveling wave parametric amplifier(JTWPA),which is a quantum-limited amplifier with high gain and large bandwidth,is the core device of large-scale measurement and control systems for quantum computing.A typical JTWPA consists of thousands of Josephson junctions connected in series to form a transmission line and hundreds of shunt LC resonators periodically loaded along the line for phase matching.Because the variation of these capacitors and inductors can be detrimental to their high-frequency characteristics,the fabrication of a JTWPA typically necessitates precise processing equipment.To guide the fabrication process and further improve the design for manufacturability,it is necessary to understand how each electronic component affects the amplifier.In this paper,we use the harmonic balance method to conduct a comprehensive study on the impact of nonuniformity and fabrication yield of the electronic components on the performance of a JTWPA.The results provide insightful and scientific guidance for device design and fabrication processes.

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.

A combined method using Lattice Boltzmann Method(LBM)and Finite Volume Method(FVM)to simulate geothermal reservoirs in Enhanced Geothermal System(EGS)

With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.

Numerical Predictions of Laminar Forced Convection Heat Transfer with and without Buoyancy Effects from an Isothermal Horizontal Flat Plate to Supercritical Nitrogen

Numerical predictions are made for Laminar Forced convection heat transfer with and without buoyancy effects for Supercritical Nitrogen flowing over an isothermal horizontal flat plate with a heated surface facing downwards.Computations are performed by varying the value ofΔT from5 to 30 K and P_(∞)/P_(cr)ratio from1.1 to 1.5.Variation of all the thermophysical properties of supercritical Nitrogen is considered.The wall temperatures are chosen in such a way that two values of Tw are less than T∗(T*is the temperature at which the fluid has a maximum value of Cp for the given pressure),one value equal to T∗and two values greater than T∗.Three different values of U∞are used to obtain Re∞range of 3.6×10_(4)to 4.74×10^(5)for forced convection without buoyancy effects and Gr_(∞)/Re^(2)_(∞)range of 0.011 to 3.107 for the case where buoyancy effects are predominant.Six different forms of correlations are proposed based on numerical predictions and are compared with actual numerical predictions.It has been found that in all six forms of correlations,the maximum deviations are found to occur in those cases where the pseudocritical temperature TT∗lies between the wall temperature and bulk fluid temperature.