In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results i...In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results in increased leakage current,decreased breakdown voltage,and lower nonlinearity,ultimately compromising their protective performance.To investigate the evolution in electrical properties during DC aging,this work developed a finite element model based on Voronoi networks and conducted accelerated aging tests on commercial varistors.Throughout the aging process,current-voltage characteristics and Schottky barrier parameters were measured and analyzed.The results indicate that when subjected to constant voltage,current flows through regions with larger grain sizes,forming discharge channels.As aging progresses,the current focus increases on these channels,leading to a decline in the varistor’s overall performance.Furthermore,analysis of the Schottky barrier parameters shows that the changes in electrical performance during aging are non-monotonic.These findings offer theoretical support for understanding the aging mechanisms and condition assessment of modern stable ZnO varistors.展开更多
This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit...This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit scheme are obtained by discretizing the equivalent integral form of the equation.For the explicit scheme with nodal values,the algebraic system can be solved by the Thomas method.For the implicit scheme with both nodal values and their derivatives,the system can be implemented by a prediction-correction procedure,where in the correction stage,an implicit formula for recovering the nodal derivatives is introduced.Taking two point boundary value problem as an example,we prove that both the explicit and implicit schemes are convergent with fourth order accuracy with respect to some standard discrete norms using the energy method.Two numerical examples demonstrate the correctness and effectiveness of the schemes,as well as the indispensability of using non-uniform meshes.展开更多
This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationa...This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationally prohibitive for time-sensitive applications,including patient-specific aneurysm hemodynamics where rapid turnaround is valuable.The proposed pipeline learns to reconstruct HR velocity-magnitude fields fromlow-resolution(LR)FEM solutions generated under the same governing equations and boundary conditions.It consistsof three modules:(i)offline pre-training of a residual network on representative vascular geometries;(ii)lightweightfine-tuning to adapt the pretrained model to geometric variability,including patient-specific aneurysm morphologies;and(iii)an unstructured-to-structured sampling strategy with region-of-interest upsampling that concentrates resolution in flow-critical zones(e.g.,the aneurysm sac)rather than the full domain.This targeted reconstruction substantiallyreduces inference and post-processing cost while preserving key HR flow features.Experiments on cerebral aneurysmmodels show that HR velocity-magnitude fields can be recovered with accuracy comparable to direct HR simulationsat less than 1%of the direct HR simulation cost per analysis(LR simulation and SR inference),while adaptation to newgeometries requires only lightweight fine-tuning with limited target-specific HR data.While clinical endpoints andadditional variables(e.g.,pressure or wall-based metrics)are left for future work,the results indicate that the proposedsurrogate SR approach can streamline FEM-based CFD workflows toward near real-time hemodynamic analysis acrossmorphologically similar vascular models.展开更多
In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is de...In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.展开更多
This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens ...This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens manufactured using extrusion-based 3D printing.Through comprehensive testing,including cyclic compression at strain rates ranging from 0.12 to 120 mm/min(0%-15%strain)and creep/relaxation experiments(10%-30%strain),the lumped parameters were independently determined using both analytical and numerical solutions of the models’differential equations,followed by cross-verification in additional experiments.Numerical solutions for creep and relaxation problems were obtained using finite element analysis,with the three-parameter Mooney-Rivlin model and Prony series employed to simulate elastic and viscous stress components,respectively.Energy dissipation per cycle was quantified during cyclic compression tests.The results demonstrate that all three models adequately describe material behavior within the 0%-15%strain range across various strain rates.Comparative analysis revealed the Burgers model’s superior performance in characterizing creep and stress relaxation at low strain levels.While Zener and Burgers model parameters from uniaxial compression showed limited applicability for energy dissipation calculations,the generalized Maxwell model effectively captured viscoelastic properties across different strain rates.Notably,parameters derived from creep tests provided a more universal assessment of dissipative properties due to optimization based on characteristic curve regions.Both parameter sets described polyurethane’s elastic-hysteretic behavior with approximately 20%error,proving significantly more accurate than the linear strain-time dependence hypothesis.Finite element analysis(FEA)complemented numerical modeling by demonstrating that while the generalized Maxwell model effectively describes initial rapid stress-strain changes,FEA provides superior characterization of steady-state processes.This computational approach yields more physically representative results compared to simplified analytical solutions,despite certain limitations in transient analysis.展开更多
Motivated by the recent experimental discovery of superconductivity in rhombohedral tetralayer graphene,we investigate the pairing mechanism arising from the density–density interactions within the random-phase appro...Motivated by the recent experimental discovery of superconductivity in rhombohedral tetralayer graphene,we investigate the pairing mechanism arising from the density–density interactions within the random-phase approximation.This approach successfully highlights the dominance of the chiral p-wave pairing between electrons with the same spin and valley index at low densities,while also predicting the superconducting range in agreement with experimental findings.Furthermore,we examine the characteristics of distinct superconducting regions:SC1 and SC2 exhibit chiral finite-momentum superconductivity with pronounced phase fluctuations,whereas SC4 displays zero-momentum spin-singlet superconductivity.展开更多
This paper presents a nonlinear micropolar nonclassical continuum theory (MPNCCT) for finite deformation, finite strain deformation physics of thermosviscoelastic solid medium with memory (polymeric micropolar solids)...This paper presents a nonlinear micropolar nonclassical continuum theory (MPNCCT) for finite deformation, finite strain deformation physics of thermosviscoelastic solid medium with memory (polymeric micropolar solids) based on classical rotations cΘand their rates. Contravariant second Piola-Kirchhoff stress and moment tensors, in conjunction with finite deformation measures derived by the authors in recent paper, are utilized in deriving the conservation and balance laws and the constitutive theories based on conjugate pairs in entropy inequality and the representation theorem. This nonlinear MPNCCT for TVES with rheology: 1) incorporates nonlinear ordered rate dissipation mechanism based on Green’s strain rates up to order n;2) also incorporates an additional ordered rate dissipation mechanism due to microconstituents, the viscosity of the medium and the rates of the symmetric part of the rotation gradient (of cΘ) tensor up to order n, referred to as micropolar dissipation or micropolar viscous dissipation mechanism;3) incorporates the primary mechanism of memory or rheology due to long chain molecules of the polymer and the viscosity of the medium by using the contravaraint second Piola-Kirchhoff stress tensor and its rates up to order m, resulting in a relaxation spectrum;4) incorporates second mechanism of memory or rheology due to nonclassical physics, interaction of microconstituents with the viscous medium and long chain molecules by considering rates of the contravariant second Piola-Kirchhoff moment tensor up to order m, resulting in relaxation of second Piola-Kirchhoff moment tensor. This results in another relaxation spectrum for the second Piola-Kirchhoff moment tensor due to microconstituents, referred to as micropolar relaxation spectrum consisting of micropolar relaxation time constants of the material. This nonlinear MPNCCT for TVES with memory is thermodynamically and mathematically consistent, and the mathematical model consisting of conservation and balance laws and the constitutive theories has closure and naturally reduces to linear MPNCCT based on infinitesimal deformation assumption. BMM is the essential balance law for all MPNCCT and is used in the present work as well. In the absence of this balance law, a valid thermodynamically and mathematically consistent nonlinear MPNCCT is not possible. The nonlinear MPNCCT based on rotations (cΘ+αΘ) and αΘ(ignoring cΘ) is not considered due to the fact that even the linear MPNCCT based on these rotations is invalid and is thermodynamically and mathematically inconsistent MPNCCT.展开更多
Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and ...Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.展开更多
A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FE...A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.展开更多
In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solut...In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.展开更多
A theory based on the superposition principle is developed to uncover the basic physics of wave behavior in a finite grating of N unit cells.The theory reveals that bound states in the continuum(BICs)of infinite quali...A theory based on the superposition principle is developed to uncover the basic physics of wave behavior in a finite grating of N unit cells.The theory reveals that bound states in the continuum(BICs)of infinite quality factor(Q-factor)can be supported by such a grating when perfect reflection is introduced at its boundaries.If geometrical perturbations are introduced into the structure,the dark BICs transform into bright quasi-BICs with finite Q-factor,maintaining spectral characteristics nearly identical to those of quasi-BICs supported by infinite gratings.When the boundaries are replaced with high-reflectivity metallic mirrors,the Q-factor of the resonant mode is reduced to be finite;however,it can be much larger than that in the corresponding nanostructure with open boundaries and can be tuned over a large range by varying the number of unit cells or boundary conditions.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
1-way multihead quantum finite state automata (1QFA(k)) can be thought of modified version of 1-way quantum finite state automata (1QFA) and k-letter quantum finite state automata (k-letter QFA) respectively. It has b...1-way multihead quantum finite state automata (1QFA(k)) can be thought of modified version of 1-way quantum finite state automata (1QFA) and k-letter quantum finite state automata (k-letter QFA) respectively. It has been shown by Moore and Crutchfield as well as Konadacs and Watrous that 1QFA can’t accept all regular language. In this paper, we show different language recognizing capabilities of our model 1-way multihead QFAs. New results presented in this paper are the following ones: 1) We show that newly introduced 1-way 2-head quantum finite state automaton (1QFA(2)) structure can accept all unary regular languages. 2) A language which can’t be accepted by 1-way deterministic 2-head finite state automaton (1DFA((2)) can be accepted by 1QFA(2) with bounded error. 3) 1QFA(2) is more powerful than 1-way reversible 2-head finite state automaton (1RMFA(2)) with respect to recognition of language.展开更多
Purpose–The brake pipe system was an essential braking component of the railway freight trains,but the existing E-type sealing rings had problems such as insufficient low-temperature resistance,poor heat stability an...Purpose–The brake pipe system was an essential braking component of the railway freight trains,but the existing E-type sealing rings had problems such as insufficient low-temperature resistance,poor heat stability and short service life.To address these issues,low-phenyl silicone rubber was prepared and tested,and the finite element analysis and experimental studies on the sealing performance of its sealing rings were carried out.Design/methodology/approach–The low-temperature resistance and thermal stability of the prepared lowphenyl silicone rubber were studied using low-temperature tensile testing,differential scanning calorimetry,dynamic thermomechanical analysis and thermogravimetric analysis.The sealing performance of the lowphenyl silicone rubber sealing ring was studied by using finite element analysis software abaqus and experiments.Findings–The prepared low-phenyl silicone rubber sealing ring possessed excellent low-temperature resistance and thermal stability.According to the finite element analysis results,the finish of the flange sealing surface and groove outer edge should be ensured,and extrusion damage should be avoided.The sealing rings were more susceptible to damage in high compression ratio and/or low-temperature environments.When the sealing effect was ensured,a small compression ratio should be selected,and rubbers with hardness and elasticity less affected by temperature should be selected.The prepared low-phenyl silicone rubber sealing ring had zero leakage at both room temperature(RT)and�508C.Originality/value–The innovation of this study is that it provides valuable data and experience for the future development of the sealing rings used in the brake pipe flange joints of the railway freight cars in China.展开更多
Total hip arthroplasty for adults with sequelae from childhood hip disorders poses significant challenges due to altered anatomy.The paper published by Oommen et al reviews the essential management strategies for thes...Total hip arthroplasty for adults with sequelae from childhood hip disorders poses significant challenges due to altered anatomy.The paper published by Oommen et al reviews the essential management strategies for these complex cases.This article explores the integration of finite element analysis(FEA)to enhance surgical precision and outcomes.FEA provides detailed biomechanical insights,aiding in preoperative planning,implant design,and surgical technique optimization.By simulating implant configurations and assessing bone quality,FEA helps in customizing implants and evaluating surgical techniques like subtrochanteric shortening osteotomy.Advanced imaging techniques,such as 3D printing,virtual reality,and augmented reality,further enhance total hip arthroplasty precision.Future research should focus on validating FEA models,developing patient-specific simulations,and promoting multidisciplinary collaboration.Integrating FEA and advanced technologies in total hip arthroplasty can improve functional outcomes,reduce complications,and enhance quality of life for patients with childhood hip disorder sequelae.展开更多
Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres wi...Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.展开更多
The phase volume fraction has an important role in the match of the strength and plasticity of dual phase steel.The different bainite contents(18–53 vol.%)in polygonal ferrite and bainite(PF+B)dual phase steel were o...The phase volume fraction has an important role in the match of the strength and plasticity of dual phase steel.The different bainite contents(18–53 vol.%)in polygonal ferrite and bainite(PF+B)dual phase steel were obtained by controlling the relaxation finish temperature during the rolling process.The effect of bainite volume fraction on the tensile deformability was systematically investigated via experiments and crystal plasticity finite element model(CPFEM)simulation.The experimental results showed that the steel showed optimal strain hardenability and strength–plasticity matching when the bainite reached 35%.The 3D-CPFEM models with the same grain size and texture characters were established to clarify the influence of stress/strain distribution on PF+B dual phase steel with different bainite contents.The simulation results indicated that an appropriate increase in the bainite content(18%–35%)did not affect the interphase strain difference,but increased the stress distribution in both phases,as a result of enhancing the coordinated deformability of two phases and improving the strength–plasticity matching.When the bainite content increased to 53%,the stress/strain difference between the two phases was greatly increased,and plastic damage between the two phases was caused by the reduction of the coordinated deformability.展开更多
This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
This paper mainly studies the well-posedness of steady incompressible impinging jet flow problem through a 3D axisymmetric finitely long nozzle.This problem originates from the physical phenomena encountered in practi...This paper mainly studies the well-posedness of steady incompressible impinging jet flow problem through a 3D axisymmetric finitely long nozzle.This problem originates from the physical phenomena encountered in practical engineering fields,such as in short take-off and vertical landing(STOVL)aircraft.Nowadays many intricate phenomena associated with impinging jet flows remain inadequately elucidated,which limits the ability to optimize aircraft design.Given a boundary condition in the inlet,the impinging jet problem is transformed into a Bernoulli-type free boundary problem according to the stream function.Then the variational method is used to study the corresponding variational problem with one parameter,thereby the wellposedness is established.The main conclusion is as follows.For a 3D axisymmetric finitely long nozzle and an infinitely long vertical wall,given an axial velocity in the inlet of nozzle,there exists a unique smooth incom‑pressible impinging jet flow such that the free boundary initiates smoothly at the endpoint of the nozzle and extends to infinity along the vertical wall at far fields.The key point is to investigate the regularity of the corner where the nozzle and the vertical axis intersect.展开更多
文摘In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results in increased leakage current,decreased breakdown voltage,and lower nonlinearity,ultimately compromising their protective performance.To investigate the evolution in electrical properties during DC aging,this work developed a finite element model based on Voronoi networks and conducted accelerated aging tests on commercial varistors.Throughout the aging process,current-voltage characteristics and Schottky barrier parameters were measured and analyzed.The results indicate that when subjected to constant voltage,current flows through regions with larger grain sizes,forming discharge channels.As aging progresses,the current focus increases on these channels,leading to a decline in the varistor’s overall performance.Furthermore,analysis of the Schottky barrier parameters shows that the changes in electrical performance during aging are non-monotonic.These findings offer theoretical support for understanding the aging mechanisms and condition assessment of modern stable ZnO varistors.
文摘This paper studies high order compact finite volume methods on non-uniform meshes for one-dimensional elliptic and parabolic differential equations with the Robin boundary conditions.An explicit scheme and an implicit scheme are obtained by discretizing the equivalent integral form of the equation.For the explicit scheme with nodal values,the algebraic system can be solved by the Thomas method.For the implicit scheme with both nodal values and their derivatives,the system can be implemented by a prediction-correction procedure,where in the correction stage,an implicit formula for recovering the nodal derivatives is introduced.Taking two point boundary value problem as an example,we prove that both the explicit and implicit schemes are convergent with fourth order accuracy with respect to some standard discrete norms using the energy method.Two numerical examples demonstrate the correctness and effectiveness of the schemes,as well as the indispensability of using non-uniform meshes.
文摘This study develops a surrogate super-resolution(SR)framework that accelerates finite element method(FEM)-based computational fluid dynamics(CFD)using deep learning.High-resolution(HR)FEM-based CFDremains computationally prohibitive for time-sensitive applications,including patient-specific aneurysm hemodynamics where rapid turnaround is valuable.The proposed pipeline learns to reconstruct HR velocity-magnitude fields fromlow-resolution(LR)FEM solutions generated under the same governing equations and boundary conditions.It consistsof three modules:(i)offline pre-training of a residual network on representative vascular geometries;(ii)lightweightfine-tuning to adapt the pretrained model to geometric variability,including patient-specific aneurysm morphologies;and(iii)an unstructured-to-structured sampling strategy with region-of-interest upsampling that concentrates resolution in flow-critical zones(e.g.,the aneurysm sac)rather than the full domain.This targeted reconstruction substantiallyreduces inference and post-processing cost while preserving key HR flow features.Experiments on cerebral aneurysmmodels show that HR velocity-magnitude fields can be recovered with accuracy comparable to direct HR simulationsat less than 1%of the direct HR simulation cost per analysis(LR simulation and SR inference),while adaptation to newgeometries requires only lightweight fine-tuning with limited target-specific HR data.While clinical endpoints andadditional variables(e.g.,pressure or wall-based metrics)are left for future work,the results indicate that the proposedsurrogate SR approach can streamline FEM-based CFD workflows toward near real-time hemodynamic analysis acrossmorphologically similar vascular models.
基金The Natural Science Foundation of Xinjiang Uygur Autonomous Region of China“RBF-Hermite difference scheme for the time-fractional kdv-Burgers equation”(2024D01C43)。
文摘In this paper,we present a finite volume trigonometric weighted essentially non-oscillatory(TWENO)scheme to solve nonlinear degenerate parabolic equations that may exhibit non-smooth solutions.The present method is developed using the trigonometric scheme,which is based on zero,first,and second moments,and the direct discontinuous Galerkin(DDG)flux is used to discretize the diffusion term.Moreover,the DDG method directly applies the weak form of the parabolic equation to each computational cell,which can better capture the characteristics of the solution,especially the discontinuous solution.Meanwhile,the third-order TVD-Runge-Kutta method is applied for temporal discretization.Finally,the effectiveness and stability of the method constructed in this paper are evaluated through numerical tests.
文摘This study presents and verifies a hybrid methodology for reliable determination of parameters in structural rheological models(Zener,Burgers,and Maxwell)describing the viscoelastic behavior of polyurethane specimens manufactured using extrusion-based 3D printing.Through comprehensive testing,including cyclic compression at strain rates ranging from 0.12 to 120 mm/min(0%-15%strain)and creep/relaxation experiments(10%-30%strain),the lumped parameters were independently determined using both analytical and numerical solutions of the models’differential equations,followed by cross-verification in additional experiments.Numerical solutions for creep and relaxation problems were obtained using finite element analysis,with the three-parameter Mooney-Rivlin model and Prony series employed to simulate elastic and viscous stress components,respectively.Energy dissipation per cycle was quantified during cyclic compression tests.The results demonstrate that all three models adequately describe material behavior within the 0%-15%strain range across various strain rates.Comparative analysis revealed the Burgers model’s superior performance in characterizing creep and stress relaxation at low strain levels.While Zener and Burgers model parameters from uniaxial compression showed limited applicability for energy dissipation calculations,the generalized Maxwell model effectively captured viscoelastic properties across different strain rates.Notably,parameters derived from creep tests provided a more universal assessment of dissipative properties due to optimization based on characteristic curve regions.Both parameter sets described polyurethane’s elastic-hysteretic behavior with approximately 20%error,proving significantly more accurate than the linear strain-time dependence hypothesis.Finite element analysis(FEA)complemented numerical modeling by demonstrating that while the generalized Maxwell model effectively describes initial rapid stress-strain changes,FEA provides superior characterization of steady-state processes.This computational approach yields more physically representative results compared to simplified analytical solutions,despite certain limitations in transient analysis.
基金supported by the National Natural Science Foundation of China (Grant Nos.12447125,12234016,and 12174317)the New Cornerstone Science Foundation。
文摘Motivated by the recent experimental discovery of superconductivity in rhombohedral tetralayer graphene,we investigate the pairing mechanism arising from the density–density interactions within the random-phase approximation.This approach successfully highlights the dominance of the chiral p-wave pairing between electrons with the same spin and valley index at low densities,while also predicting the superconducting range in agreement with experimental findings.Furthermore,we examine the characteristics of distinct superconducting regions:SC1 and SC2 exhibit chiral finite-momentum superconductivity with pronounced phase fluctuations,whereas SC4 displays zero-momentum spin-singlet superconductivity.
文摘This paper presents a nonlinear micropolar nonclassical continuum theory (MPNCCT) for finite deformation, finite strain deformation physics of thermosviscoelastic solid medium with memory (polymeric micropolar solids) based on classical rotations cΘand their rates. Contravariant second Piola-Kirchhoff stress and moment tensors, in conjunction with finite deformation measures derived by the authors in recent paper, are utilized in deriving the conservation and balance laws and the constitutive theories based on conjugate pairs in entropy inequality and the representation theorem. This nonlinear MPNCCT for TVES with rheology: 1) incorporates nonlinear ordered rate dissipation mechanism based on Green’s strain rates up to order n;2) also incorporates an additional ordered rate dissipation mechanism due to microconstituents, the viscosity of the medium and the rates of the symmetric part of the rotation gradient (of cΘ) tensor up to order n, referred to as micropolar dissipation or micropolar viscous dissipation mechanism;3) incorporates the primary mechanism of memory or rheology due to long chain molecules of the polymer and the viscosity of the medium by using the contravaraint second Piola-Kirchhoff stress tensor and its rates up to order m, resulting in a relaxation spectrum;4) incorporates second mechanism of memory or rheology due to nonclassical physics, interaction of microconstituents with the viscous medium and long chain molecules by considering rates of the contravariant second Piola-Kirchhoff moment tensor up to order m, resulting in relaxation of second Piola-Kirchhoff moment tensor. This results in another relaxation spectrum for the second Piola-Kirchhoff moment tensor due to microconstituents, referred to as micropolar relaxation spectrum consisting of micropolar relaxation time constants of the material. This nonlinear MPNCCT for TVES with memory is thermodynamically and mathematically consistent, and the mathematical model consisting of conservation and balance laws and the constitutive theories has closure and naturally reduces to linear MPNCCT based on infinitesimal deformation assumption. BMM is the essential balance law for all MPNCCT and is used in the present work as well. In the absence of this balance law, a valid thermodynamically and mathematically consistent nonlinear MPNCCT is not possible. The nonlinear MPNCCT based on rotations (cΘ+αΘ) and αΘ(ignoring cΘ) is not considered due to the fact that even the linear MPNCCT based on these rotations is invalid and is thermodynamically and mathematically inconsistent MPNCCT.
基金supported by the National Natural Science Foundation of China(No.42372287 and No.U24A20178)the Fundamental Research Funds for the Central Universities CHD(No.2024SHEEAR002)+3 种基金the Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shaanxi Province(No.2020024)the China Postdoctoral Science Foundation(GZC20232955,2024M753472,and 2024MD763937)the Science-Technology Foundation for Young Scientists of Gansu Province,China(No.24JRRA097)the Study of biodiversity survey and limiting factor analysis of Yinkentala(2023ZL01).
文摘Accurately simulating water flow movement in vadose zone is crucial for effective water resources assessment.Richards'equation,which describes the movement of water flow in the vadose zone,is highly nonlinear and challenging to solve.Existing numerical methods often face issues such as numerical dispersion,oscillation,and mass non-conservation when spatial and temporal discretization conditions are not appropriately configured.To address these problems and achieve accurate and stable numerical solutions,a finite analytic method based on water content-based Richards'equation(FAM-W)is proposed.The performance of the FAM-W is compared with analytical solutions,Finite Difference Method(FDM),and Finite Analytic Method based on the pressure Head-based Richards'equation(FAM-H).Compared to analytical solution and other numerical methods(FDM and FAM-H),FAM-W demonstrates superior accuracy and efficiency in controlling mass balance errors,regardless of spatial step sizes.This study introduces a novel approach for modelling water flow in the vadose zone,offering significant benefits for water resources management.
基金the National Natural Science Foundation of China(No.11672238)the 111 Project(No.BP0719007)the Shaanxi Province Natural Science Foundation(No.2020JZ-06)for the financial support.
文摘A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.
基金supported by the National Natural Science Foundation of China(No.12002290)。
文摘In the field of discretization-based meshfree/meshless methods,the improvements in the higher-order consistency,stability,and computational efficiency are of great concerns in computational science and numerical solutions to partial differential equations.Various alternative numerical methods of the finite particle method(FPM)frame have been extended from mathematical theories to numerical applications separately.As a comprehensive numerical scheme,this study suggests a unified resolved program for numerically investigating their accuracy,stability,consistency,computational efficiency,and practical applicability in industrial engineering contexts.The high-order finite particle method(HFPM)and corrected methods based on the multivariate Taylor series expansion are constructed and analyzed to investigate the whole applicability in different benchmarks of computational fluid dynamics.Specifically,four benchmarks are designed purposefully from statical exact solutions to multifaceted hydrodynamic tests,which possess different numerical performances on the particle consistency,numerical discretized forms,particle distributions,and transient time evolutional stabilities.This study offers a numerical reference for the current unified resolved program.
基金supported by the National Natural Science Foundation of China(Grant Nos.11874270 and 12174228)the Shenzhen Basic Research Special Project(Grant No.JCYJ20240813141606009)。
文摘A theory based on the superposition principle is developed to uncover the basic physics of wave behavior in a finite grating of N unit cells.The theory reveals that bound states in the continuum(BICs)of infinite quality factor(Q-factor)can be supported by such a grating when perfect reflection is introduced at its boundaries.If geometrical perturbations are introduced into the structure,the dark BICs transform into bright quasi-BICs with finite Q-factor,maintaining spectral characteristics nearly identical to those of quasi-BICs supported by infinite gratings.When the boundaries are replaced with high-reflectivity metallic mirrors,the Q-factor of the resonant mode is reduced to be finite;however,it can be much larger than that in the corresponding nanostructure with open boundaries and can be tuned over a large range by varying the number of unit cells or boundary conditions.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
文摘1-way multihead quantum finite state automata (1QFA(k)) can be thought of modified version of 1-way quantum finite state automata (1QFA) and k-letter quantum finite state automata (k-letter QFA) respectively. It has been shown by Moore and Crutchfield as well as Konadacs and Watrous that 1QFA can’t accept all regular language. In this paper, we show different language recognizing capabilities of our model 1-way multihead QFAs. New results presented in this paper are the following ones: 1) We show that newly introduced 1-way 2-head quantum finite state automaton (1QFA(2)) structure can accept all unary regular languages. 2) A language which can’t be accepted by 1-way deterministic 2-head finite state automaton (1DFA((2)) can be accepted by 1QFA(2) with bounded error. 3) 1QFA(2) is more powerful than 1-way reversible 2-head finite state automaton (1RMFA(2)) with respect to recognition of language.
基金supported by the Science and Technology Research and Development Plan of the China State Railway Group Company Limited(No.Q2023J012).
文摘Purpose–The brake pipe system was an essential braking component of the railway freight trains,but the existing E-type sealing rings had problems such as insufficient low-temperature resistance,poor heat stability and short service life.To address these issues,low-phenyl silicone rubber was prepared and tested,and the finite element analysis and experimental studies on the sealing performance of its sealing rings were carried out.Design/methodology/approach–The low-temperature resistance and thermal stability of the prepared lowphenyl silicone rubber were studied using low-temperature tensile testing,differential scanning calorimetry,dynamic thermomechanical analysis and thermogravimetric analysis.The sealing performance of the lowphenyl silicone rubber sealing ring was studied by using finite element analysis software abaqus and experiments.Findings–The prepared low-phenyl silicone rubber sealing ring possessed excellent low-temperature resistance and thermal stability.According to the finite element analysis results,the finish of the flange sealing surface and groove outer edge should be ensured,and extrusion damage should be avoided.The sealing rings were more susceptible to damage in high compression ratio and/or low-temperature environments.When the sealing effect was ensured,a small compression ratio should be selected,and rubbers with hardness and elasticity less affected by temperature should be selected.The prepared low-phenyl silicone rubber sealing ring had zero leakage at both room temperature(RT)and�508C.Originality/value–The innovation of this study is that it provides valuable data and experience for the future development of the sealing rings used in the brake pipe flange joints of the railway freight cars in China.
文摘Total hip arthroplasty for adults with sequelae from childhood hip disorders poses significant challenges due to altered anatomy.The paper published by Oommen et al reviews the essential management strategies for these complex cases.This article explores the integration of finite element analysis(FEA)to enhance surgical precision and outcomes.FEA provides detailed biomechanical insights,aiding in preoperative planning,implant design,and surgical technique optimization.By simulating implant configurations and assessing bone quality,FEA helps in customizing implants and evaluating surgical techniques like subtrochanteric shortening osteotomy.Advanced imaging techniques,such as 3D printing,virtual reality,and augmented reality,further enhance total hip arthroplasty precision.Future research should focus on validating FEA models,developing patient-specific simulations,and promoting multidisciplinary collaboration.Integrating FEA and advanced technologies in total hip arthroplasty can improve functional outcomes,reduce complications,and enhance quality of life for patients with childhood hip disorder sequelae.
基金the financial support provided by Tianfu Yongxing Laboratory Organized Research Project Funding(No.2023CXXM01)the ARC linkage program(No.LP200100420).
文摘Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.
基金supported by the Project of Liaoning Marine Economic Development(Development of high strength pipeline steel for submarine oil and gas transmission)State Key Laboratory of Metal Material for Marine Equipment and Application Funding(No.SKLMEA-K202205).
文摘The phase volume fraction has an important role in the match of the strength and plasticity of dual phase steel.The different bainite contents(18–53 vol.%)in polygonal ferrite and bainite(PF+B)dual phase steel were obtained by controlling the relaxation finish temperature during the rolling process.The effect of bainite volume fraction on the tensile deformability was systematically investigated via experiments and crystal plasticity finite element model(CPFEM)simulation.The experimental results showed that the steel showed optimal strain hardenability and strength–plasticity matching when the bainite reached 35%.The 3D-CPFEM models with the same grain size and texture characters were established to clarify the influence of stress/strain distribution on PF+B dual phase steel with different bainite contents.The simulation results indicated that an appropriate increase in the bainite content(18%–35%)did not affect the interphase strain difference,but increased the stress distribution in both phases,as a result of enhancing the coordinated deformability of two phases and improving the strength–plasticity matching.When the bainite content increased to 53%,the stress/strain difference between the two phases was greatly increased,and plastic damage between the two phases was caused by the reduction of the coordinated deformability.
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
文摘This paper mainly studies the well-posedness of steady incompressible impinging jet flow problem through a 3D axisymmetric finitely long nozzle.This problem originates from the physical phenomena encountered in practical engineering fields,such as in short take-off and vertical landing(STOVL)aircraft.Nowadays many intricate phenomena associated with impinging jet flows remain inadequately elucidated,which limits the ability to optimize aircraft design.Given a boundary condition in the inlet,the impinging jet problem is transformed into a Bernoulli-type free boundary problem according to the stream function.Then the variational method is used to study the corresponding variational problem with one parameter,thereby the wellposedness is established.The main conclusion is as follows.For a 3D axisymmetric finitely long nozzle and an infinitely long vertical wall,given an axial velocity in the inlet of nozzle,there exists a unique smooth incom‑pressible impinging jet flow such that the free boundary initiates smoothly at the endpoint of the nozzle and extends to infinity along the vertical wall at far fields.The key point is to investigate the regularity of the corner where the nozzle and the vertical axis intersect.
