In this paper,we consider the truncated multiplicity finite range set problem of meromorphic functions on some complex disc.By using the value distribution theory of meromorphic functions,we establish a second main th...In this paper,we consider the truncated multiplicity finite range set problem of meromorphic functions on some complex disc.By using the value distribution theory of meromorphic functions,we establish a second main theorem for meromorphic functions with finite growth index which share meromorphic functions(may not be small functions).As its application,we also extend the result of a finite range set with truncated multiplicity.展开更多
This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy ...This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of ...For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.展开更多
The stability analysis of a finite Stokes layer is of practical importance in flow control. In the present work, the instantaneous stability of a finite Stokes layer with layer interactions is studied via a linear sta...The stability analysis of a finite Stokes layer is of practical importance in flow control. In the present work, the instantaneous stability of a finite Stokes layer with layer interactions is studied via a linear stability analysis of the frozen phases of the base flow. The oscillations of two plates can have different velocity amplitudes, initial phases, and frequencies. The effects of the Stokes-layer interactions on the stability when two plates oscillate synchronously are analyzed. The growth rates of two most unstable modes when δ < 0.12 are almost equal, and δ = δ*/h*, where δ*and h*are the Stokes-layer thickness and the half height of the channel, respectively. However, their vorticities are different. The vorticity of the most unstable mode is symmetric, while the other is asymmetric. The Stokes-layer interactions have a destabilizing effect on the most unstable mode when δ < 0.68, and have a stabilizing effect when δ > 0.68. However, the interactions always have a stabilizing effect on the other unstable mode. It is explained that one of the two unstable modes has much higher dissipation than the other one when the Stokes-layer interactions are strong. We also find that the stability of the Stokes layer is closely related to the inflectional points of the base-flow velocity profile. The effects of inconsistent velocity-amplitude, initial phase, and frequency of the oscillations on the stability are analyzed. The energy of the most unstable eigenvector is mainly distributed near the plate of higher velocity amplitude or higher oscillation frequency. The effects of the initial phase difference are complicated because the base-flow velocity is extremely sensitive to the initial phase.展开更多
Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyp...Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyperelastic-magnetic coupling problems is proposed.Since the solution space consists of two-phase domains,two separate networks are constructed to independently predict the solution for each phase region.In addition,a conscious point allocation strategy is incorporated to enhance the prediction precision of the PINN in regions characterized by sharp gradients.With the developed framework,the magnetic fields and deformation fields of magnetorheological elastomers(MREs)are solved under the control of hyperelastic-magnetic coupling equations.Illustrative examples are provided and contrasted with the reference results to validate the predictive accuracy of the proposed framework.Moreover,the advantages of the proposed framework in solving hyperelastic-magnetic coupling problems are validated,particularly in handling small data sets,as well as its ability in swiftly and precisely forecasting magnetostrictive motion.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and ...This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.展开更多
High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric fini...High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric finite element method.Firstly,the physical field is approximated by uniform B-spline interpolation,while geometry is represented by non-uniform rational B-spline interpolation.By introducing a transformation matrix,elements of types C^(0)and C^(1)are constructed in the isogeometric finite element method.Subsequently,the corresponding calculation formats for one-dimensional bars,beams,and two-dimensional linear elasticity in the isogeometric finite element method are derived through variational principles and parameter mapping.The proposed method combines element construction techniques of the finite element method with geometric construction techniques of isogeometric analysis,eliminating the need for mesh generation and maintaining flexibility in element construction.Two elements with interpolation characteristics are constructed in the method so that boundary conditions and connections between elements can be processed like the finite element method.Finally,the test results of several examples show that:(1)Under the same degree and element node numbers,the constructed elements are almost consistent with the results obtained by traditional finite element method;(2)For bar problems with large local field variations and beam problems with variable cross-sections,high-degree and multi-nodes elements constructed can achieve high computational accuracy with fewer degrees of freedom than finite element method;(3)The computational efficiency of isogeometric finite element method is higher than finite element method under similar degrees of freedom,while as degrees of freedom increase,the computational efficiency between the two is similar.展开更多
The diameter distribution function(DDF)is a crucial tool for accurately predicting stand carbon storage(CS).The current key issue,however,is how to construct a high-precision DDF based on stand factors,site quality,an...The diameter distribution function(DDF)is a crucial tool for accurately predicting stand carbon storage(CS).The current key issue,however,is how to construct a high-precision DDF based on stand factors,site quality,and aridity index to predict stand CS in multi-species mixed forests with complex structures.This study used data from70 survey plots for mixed broadleaf Populus davidiana and Betula platyphylla forests in the Mulan Rangeland State Forest,Hebei Province,China,to construct the DDF based on maximum likelihood estimation and finite mixture model(FMM).Ordinary least squares(OLS),linear seemingly unrelated regression(LSUR),and back propagation neural network(BPNN)were used to investigate the influences of stand factors,site quality,and aridity index on the shape and scale parameters of DDF and predicted stand CS of mixed broadleaf forests.The results showed that FMM accurately described the stand-level diameter distribution of the mixed P.davidiana and B.platyphylla forests;whereas the Weibull function constructed by MLE was more accurate in describing species-level diameter distribution.The combined variable of quadratic mean diameter(Dq),stand basal area(BA),and site quality improved the accuracy of the shape parameter models of FMM;the combined variable of Dq,BA,and De Martonne aridity index improved the accuracy of the scale parameter models.Compared to OLS and LSUR,the BPNN had higher accuracy in the re-parameterization process of FMM.OLS,LSUR,and BPNN overestimated the CS of P.davidiana but underestimated the CS of B.platyphylla in the large diameter classes(DBH≥18 cm).BPNN accurately estimated stand-and species-level CS,but it was more suitable for estimating stand-level CS compared to species-level CS,thereby providing a scientific basis for the optimization of stand structure and assessment of carbon sequestration capacity in mixed broadleaf forests.展开更多
Thermal quantities,including the the entropy density and gluon spectrum,of quark matter within a box that is finite in the longitudinal direction are calculated using a bag model.Under the assumption of entropy conser...Thermal quantities,including the the entropy density and gluon spectrum,of quark matter within a box that is finite in the longitudinal direction are calculated using a bag model.Under the assumption of entropy conservation,the corresponding gluon dissociation rate of J/ψis studied.It reaches a maximum at a certain longitudinal size𝐿m,below which the suppression is weak even if the temperature becomes higher than that without the finite size effect,and above which the dissociation rate approaches to the thermodynamic limit gradually with increasing longitudinal size of the fireball.展开更多
Semiconductor devices are often operated at elevated temperatures that are well above zero Kelvin,which is the temperature in most first-principles density functional calculations.Computational approaches to com-putin...Semiconductor devices are often operated at elevated temperatures that are well above zero Kelvin,which is the temperature in most first-principles density functional calculations.Computational approaches to com-puting and understanding the properties of semiconductors at finite temperatures are thus in critical demand.In this review,we discuss the recent progress in computationally assessing the electronic and phononic band structures of semiconductors at finite temperatures.As an emerging semiconductor with particularly strong temperature-induced renormalization of the electronic and phononic band structures,halide perovskites are used as a representative example to demonstrate how computational advances may help to understand the band struc-tures at elevated temperatures.Finally,we briefly illustrate the remaining computational challenges and outlook promising research directions that may help to guide future research in this field.展开更多
An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection...An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.展开更多
Dear Editor,This letter is concerned with the role of recurrent neural networks(RNNs)on the controller design for a class of nonlinear systems.Inspired by the architectures of RNNs,the system states are stacked accord...Dear Editor,This letter is concerned with the role of recurrent neural networks(RNNs)on the controller design for a class of nonlinear systems.Inspired by the architectures of RNNs,the system states are stacked according to the dynamic along with time while the controller is represented as the neural network output.To build the bridge between RNNs and finite-time controller,a novel activation function is imposed on RNNs to drive the convergence of states at finite-time and propel the overall control process smoother.Rigorous stability proof is briefly provided for the convergence of the proposed finite-time controller.At last,a numerical simulation example is presented to illustrate the efficiency of the proposed strategy.Neural networks can be classified as static(feedforward)and dynamic(recurrent)nets[1].The former nets do not perform well in dealing with training data and using any information of the local data structure[2].In contrast to the feedforward neural networks,RNNs are constituted by high dimensional hidden states with dynamics.展开更多
Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale pr...Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolso...In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.展开更多
In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication q...In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication quality.Tubular hydrogenerators are considered the optimal generator for exploiting low-head,high-flow hydro resources,and they have seen increasingly widespread application in China's power systems recent years.However,owing to the compact and constrained internal space of such generators,their internal magnetic-field harmonics are pronounced.Therefore,accurate calculation of their THD and THF is crucial during the analysis and design stages to ensure the quality of power communication.Especially in the electromagnetic field finite element modeling analysis of such generators,the type and order of the finite element meshes may have a significant impact on the THD and THF calculation results,which warrants in-depth research.To address this,this study takes a real 34 MW large tubular hydrogenerator as an example,and establishes its electromagnetic field finite element model under no-load conditions.Two types of meshes,five mesh densities,and two mesh orders are analyzed to reveal the effect of electromagnetic field finite element mesh types and orders on the calculation results of THD and THF for such generators.展开更多
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 t...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.展开更多
基金Supported by National Natural Science Foundation of China(12061041)Jiangxi Provincial Natural Science Foundation(20232BAB201003).
文摘In this paper,we consider the truncated multiplicity finite range set problem of meromorphic functions on some complex disc.By using the value distribution theory of meromorphic functions,we establish a second main theorem for meromorphic functions with finite growth index which share meromorphic functions(may not be small functions).As its application,we also extend the result of a finite range set with truncated multiplicity.
文摘This paper presents a mathematical model consisting of conservation and balance laws (CBL) of classical continuum mechanics (CCM) and ordered rate constitutive theories in Lagrangian description derived using entropy inequality and the representation theorem for thermoviscoelastic solids (TVES) with rheology. The CBL and the constitutive theories take into account finite deformation and finite strain deformation physics and are based on contravariant deviatoric second Piola-Kirchhoff stress tensor and its work conjugate covariant Green’s strain tensor and their material derivatives of up to order m and n respectively. All published works on nonlinear dynamics of TVES with rheology are mostly based on phenomenological mathematical models. In rare instances, some aspects of CBL are used but are incorrectly altered to obtain mass, stiffness and damping matrices using space-time decoupled approaches. In the work presented in this paper, we show that this is not possible using CBL of CCM for TVES with rheology. Thus, the mathematical models used currently in the published works are not the correct description of the physics of nonlinear dynamics of TVES with rheology. The mathematical model used in the present work is strictly based on the CBL of CCM and is thermodynamically and mathematically consistent and the space-time coupled finite element methodology used in this work is unconditionally stable and provides solutions with desired accuracy and is ideally suited for nonlinear dynamics of TVES with memory. The work in this paper is the first presentation of a mathematical model strictly based on CBL of CCM and the solution of the mathematical model is obtained using unconditionally stable space-time coupled computational methodology that provides control over the errors in the evolution. Both space-time coupled and space-time decoupled finite element formulations are considered for obtaining solutions of the IVPs described by the mathematical model and are presented in the paper. Factors or the physics influencing dynamic response and dynamic bifurcation for TVES with rheology are identified and are also demonstrated through model problem studies. A simple model problem consisting of a rod (1D) of TVES material with memory fixed at one end and subjected to harmonic excitation at the other end is considered to study nonlinear dynamics of TVES with rheology, frequency response as well as dynamic bifurcation phenomenon.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金supported by National Natural Science Foundation of China(11771257)the Shandong Provincial Natural Science Foundation of China(ZR2023YQ002,ZR2023MA007,ZR2021MA004)。
文摘For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.
基金Project supported by the National Natural Science Foundation of China (No. 11402211)。
文摘The stability analysis of a finite Stokes layer is of practical importance in flow control. In the present work, the instantaneous stability of a finite Stokes layer with layer interactions is studied via a linear stability analysis of the frozen phases of the base flow. The oscillations of two plates can have different velocity amplitudes, initial phases, and frequencies. The effects of the Stokes-layer interactions on the stability when two plates oscillate synchronously are analyzed. The growth rates of two most unstable modes when δ < 0.12 are almost equal, and δ = δ*/h*, where δ*and h*are the Stokes-layer thickness and the half height of the channel, respectively. However, their vorticities are different. The vorticity of the most unstable mode is symmetric, while the other is asymmetric. The Stokes-layer interactions have a destabilizing effect on the most unstable mode when δ < 0.68, and have a stabilizing effect when δ > 0.68. However, the interactions always have a stabilizing effect on the other unstable mode. It is explained that one of the two unstable modes has much higher dissipation than the other one when the Stokes-layer interactions are strong. We also find that the stability of the Stokes layer is closely related to the inflectional points of the base-flow velocity profile. The effects of inconsistent velocity-amplitude, initial phase, and frequency of the oscillations on the stability are analyzed. The energy of the most unstable eigenvector is mainly distributed near the plate of higher velocity amplitude or higher oscillation frequency. The effects of the initial phase difference are complicated because the base-flow velocity is extremely sensitive to the initial phase.
基金supported by the National Natural Science Foundation of China(Nos.12072105 and 11932006)。
文摘Recently,numerous studies have demonstrated that the physics-informed neural network(PINN)can effectively and accurately resolve hyperelastic finite deformation problems.In this paper,a PINN framework for tackling hyperelastic-magnetic coupling problems is proposed.Since the solution space consists of two-phase domains,two separate networks are constructed to independently predict the solution for each phase region.In addition,a conscious point allocation strategy is incorporated to enhance the prediction precision of the PINN in regions characterized by sharp gradients.With the developed framework,the magnetic fields and deformation fields of magnetorheological elastomers(MREs)are solved under the control of hyperelastic-magnetic coupling equations.Illustrative examples are provided and contrasted with the reference results to validate the predictive accuracy of the proposed framework.Moreover,the advantages of the proposed framework in solving hyperelastic-magnetic coupling problems are validated,particularly in handling small data sets,as well as its ability in swiftly and precisely forecasting magnetostrictive motion.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金Project supported by the National Natural Science Foundation of China(Nos.12072209,U21A2043012192211)+1 种基金the Natural Science Foundation of Hebei Province of China(No.A2020210009)the S&T Program of Hebei Province of China(No.225676162GH)。
文摘This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.
基金funded by the Zhejiang Province Science and Technology Plan Project under grant number 2023C01069the Hebei Provincial Program on Key Basic Research Project under grant number 23311808Dthe Wenzhou Major Science and Technology Innovation Project of China under grant number ZG2022004。
文摘High-performance finite element research has always been a major focus of finite element method studies.This article introduces isogeometric analysis into the finite element method and proposes a new isogeometric finite element method.Firstly,the physical field is approximated by uniform B-spline interpolation,while geometry is represented by non-uniform rational B-spline interpolation.By introducing a transformation matrix,elements of types C^(0)and C^(1)are constructed in the isogeometric finite element method.Subsequently,the corresponding calculation formats for one-dimensional bars,beams,and two-dimensional linear elasticity in the isogeometric finite element method are derived through variational principles and parameter mapping.The proposed method combines element construction techniques of the finite element method with geometric construction techniques of isogeometric analysis,eliminating the need for mesh generation and maintaining flexibility in element construction.Two elements with interpolation characteristics are constructed in the method so that boundary conditions and connections between elements can be processed like the finite element method.Finally,the test results of several examples show that:(1)Under the same degree and element node numbers,the constructed elements are almost consistent with the results obtained by traditional finite element method;(2)For bar problems with large local field variations and beam problems with variable cross-sections,high-degree and multi-nodes elements constructed can achieve high computational accuracy with fewer degrees of freedom than finite element method;(3)The computational efficiency of isogeometric finite element method is higher than finite element method under similar degrees of freedom,while as degrees of freedom increase,the computational efficiency between the two is similar.
基金funded by the National Key Research and Development Program of China(No.2022YFD2200503-02)。
文摘The diameter distribution function(DDF)is a crucial tool for accurately predicting stand carbon storage(CS).The current key issue,however,is how to construct a high-precision DDF based on stand factors,site quality,and aridity index to predict stand CS in multi-species mixed forests with complex structures.This study used data from70 survey plots for mixed broadleaf Populus davidiana and Betula platyphylla forests in the Mulan Rangeland State Forest,Hebei Province,China,to construct the DDF based on maximum likelihood estimation and finite mixture model(FMM).Ordinary least squares(OLS),linear seemingly unrelated regression(LSUR),and back propagation neural network(BPNN)were used to investigate the influences of stand factors,site quality,and aridity index on the shape and scale parameters of DDF and predicted stand CS of mixed broadleaf forests.The results showed that FMM accurately described the stand-level diameter distribution of the mixed P.davidiana and B.platyphylla forests;whereas the Weibull function constructed by MLE was more accurate in describing species-level diameter distribution.The combined variable of quadratic mean diameter(Dq),stand basal area(BA),and site quality improved the accuracy of the shape parameter models of FMM;the combined variable of Dq,BA,and De Martonne aridity index improved the accuracy of the scale parameter models.Compared to OLS and LSUR,the BPNN had higher accuracy in the re-parameterization process of FMM.OLS,LSUR,and BPNN overestimated the CS of P.davidiana but underestimated the CS of B.platyphylla in the large diameter classes(DBH≥18 cm).BPNN accurately estimated stand-and species-level CS,but it was more suitable for estimating stand-level CS compared to species-level CS,thereby providing a scientific basis for the optimization of stand structure and assessment of carbon sequestration capacity in mixed broadleaf forests.
基金supported by the National Natural Science Foundation of China(Grant No.12175165)。
文摘Thermal quantities,including the the entropy density and gluon spectrum,of quark matter within a box that is finite in the longitudinal direction are calculated using a bag model.Under the assumption of entropy conservation,the corresponding gluon dissociation rate of J/ψis studied.It reaches a maximum at a certain longitudinal size𝐿m,below which the suppression is weak even if the temperature becomes higher than that without the finite size effect,and above which the dissociation rate approaches to the thermodynamic limit gradually with increasing longitudinal size of the fireball.
基金supported by the National Natural Science Foundation of China(Grant Nos.11991060,52172136,12088101,12074029,and U2230402).
文摘Semiconductor devices are often operated at elevated temperatures that are well above zero Kelvin,which is the temperature in most first-principles density functional calculations.Computational approaches to com-puting and understanding the properties of semiconductors at finite temperatures are thus in critical demand.In this review,we discuss the recent progress in computationally assessing the electronic and phononic band structures of semiconductors at finite temperatures.As an emerging semiconductor with particularly strong temperature-induced renormalization of the electronic and phononic band structures,halide perovskites are used as a representative example to demonstrate how computational advances may help to understand the band struc-tures at elevated temperatures.Finally,we briefly illustrate the remaining computational challenges and outlook promising research directions that may help to guide future research in this field.
基金The National Key Research and Development Program of China under contract Nos 2022YFC3104804,2021YFC3101501,and 2017YFC1404103the National Programme on Global Change and Air-Sea Interaction of China under contract No.GASI-IPOVAI-04the National Natural Science Foundation of China under contract Nos 41876014,41606039,and 11801402.
文摘An anisotropic diffusion filter can be used to model a flow-dependent background error covariance matrix,which can be achieved by solving the advection-diffusion equation.Because of the directionality of the advection term,the discrete method needs to be chosen very carefully.The finite analytic method is an alternative scheme to solve the advection-diffusion equation.As a combination of analytical and numerical methods,it not only has high calculation accuracy but also holds the characteristic of the auto upwind.To demonstrate its ability,the one-dimensional steady and unsteady advection-diffusion equation numerical examples are respectively solved by the finite analytic method.The more widely used upwind difference method is used as a control approach.The result indicates that the finite analytic method has higher accuracy than the upwind difference method.For the two-dimensional case,the finite analytic method still has a better performance.In the three-dimensional variational assimilation experiment,the finite analytic method can effectively improve analysis field accuracy,and its effect is significantly better than the upwind difference and the central difference method.Moreover,it is still a more effective solution method in the strong flow region where the advective-diffusion filter performs most prominently.
文摘Dear Editor,This letter is concerned with the role of recurrent neural networks(RNNs)on the controller design for a class of nonlinear systems.Inspired by the architectures of RNNs,the system states are stacked according to the dynamic along with time while the controller is represented as the neural network output.To build the bridge between RNNs and finite-time controller,a novel activation function is imposed on RNNs to drive the convergence of states at finite-time and propel the overall control process smoother.Rigorous stability proof is briefly provided for the convergence of the proposed finite-time controller.At last,a numerical simulation example is presented to illustrate the efficiency of the proposed strategy.Neural networks can be classified as static(feedforward)and dynamic(recurrent)nets[1].The former nets do not perform well in dealing with training data and using any information of the local data structure[2].In contrast to the feedforward neural networks,RNNs are constituted by high dimensional hidden states with dynamics.
基金Supported by Science Center for Gas Turbine Project of China (Grant No.P2022-B-IV-014-001)Frontier Leading Technology Basic Research Special Project of Jiangsu Province of China (Grant No.BK20212007)the BIT Research and Innovation Promoting Project of China (Grant No.2022YCXZ019)。
文摘Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.
基金supported by the Key Laboratory of Road Construction Technology and Equipment(Chang’an University,No.300102253502)the Natural Science Foundation of Shandong Province of China(GrantNo.ZR2022YQ06)the Development Plan of Youth Innovation Team in Colleges and Universities of Shandong Province(Grant No.2022KJ140).
文摘In this study,we propose an efficient numerical framework to attain the solution of the extended Fisher-Kolmogorov(EFK)problem.The temporal derivative in the EFK equation is approximated by utilizing the Crank-Nicolson scheme.Following temporal discretization,the generalized finite difference method(GFDM)with supplementary nodes is utilized to address the nonlinear boundary value problems at each time node.These supplementary nodes are distributed along the boundary to match the number of boundary nodes.By incorporating supplementary nodes,the resulting nonlinear algebraic equations can effectively satisfy the governing equation and boundary conditions of the EFK equation.To demonstrate the efficacy of our approach,we present three numerical examples showcasing its performance in solving this nonlinear problem.
基金sponsored by the National Natural Science Foundation,Youth Foundation of China,Grant/Award Number:51607146Sichuan Natural Sciences Fund,Grant/Award Number:2023NSFSC0295。
文摘In generator design field,waveform total harmonic distortion(THD)and telephone harmonic factor(THF)are parameters commonly used to measure the impact of generator no-load voltage harmonics on the power communication quality.Tubular hydrogenerators are considered the optimal generator for exploiting low-head,high-flow hydro resources,and they have seen increasingly widespread application in China's power systems recent years.However,owing to the compact and constrained internal space of such generators,their internal magnetic-field harmonics are pronounced.Therefore,accurate calculation of their THD and THF is crucial during the analysis and design stages to ensure the quality of power communication.Especially in the electromagnetic field finite element modeling analysis of such generators,the type and order of the finite element meshes may have a significant impact on the THD and THF calculation results,which warrants in-depth research.To address this,this study takes a real 34 MW large tubular hydrogenerator as an example,and establishes its electromagnetic field finite element model under no-load conditions.Two types of meshes,five mesh densities,and two mesh orders are analyzed to reveal the effect of electromagnetic field finite element mesh types and orders on the calculation results of THD and THF for such generators.
基金This work was financially supported by the Key Science and Technology Project of Longmen Laboratory(No.LMYLKT-001)Innovation and Entrepreneurship Training Program for College Students of Henan Province(No.202310464050)。
文摘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.