To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
A new higher-order shear deformation theory based on global-local superposition technique is developed. The theory satisfies the free surface conditions and the geometric and stress continuity conditions at interfaces...A new higher-order shear deformation theory based on global-local superposition technique is developed. The theory satisfies the free surface conditions and the geometric and stress continuity conditions at interfaces. The global displacement components are of the Reddy theory and local components are of the internal first to third-order terms in each layer. A two-node beam element based on this theory is proposed. The solutions are compared with 3D-elasticity solutions. Numerical results show that present beam element has higher computational efficiency and higher accuracy.展开更多
With the development of aeronautic and astronautic techniques, radiation becomes much more significant while the structure is exposed to the higher and higher temperature. Most of the current finite element software p...With the development of aeronautic and astronautic techniques, radiation becomes much more significant while the structure is exposed to the higher and higher temperature. Most of the current finite element software packages treat it using the net-radiation method or absorbed radiation method based on the assumption of isothermal surface with uniform radiation heat flux, which brings the conflict between the precision and the quantity of grids. Using integral method to compute the variable radiation heat flux in higher-order finite element, the precision can be improved greatly while using the same quantity of grids, because it is more consistent with the distribution of real temperature. In this paper, the integral is only processed on the same integral points as those used for solving the finite element equations, so it may be of high efficiency. In an academic testing model, the result is contrast to which get in ANSYS, proving the high precision of the method. Then an actual sandwich panel used in the thermal protection system is analyzed with the method, and the error is comparatively low to the analytical answer while the computation being of high efficiency.展开更多
This article briefly reviews the topic of complex network synchronization,with its graph-theoretic criterion,showing that the homogeneous and symmetrical network structures are essential for optimal synchronization.Fu...This article briefly reviews the topic of complex network synchronization,with its graph-theoretic criterion,showing that the homogeneous and symmetrical network structures are essential for optimal synchronization.Furthermore,it briefly reviews the notion of higher-order network topologies and shows their promising potential in application to evaluating the optimality of network synchronizability.展开更多
The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical mo...The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.展开更多
In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introduc...In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.展开更多
This paper deals with extensions of higher-order optimality conditions for scalar optimization to multiobjective optimization.A type of directional derivatives for a multiobjective function is proposed,and with this n...This paper deals with extensions of higher-order optimality conditions for scalar optimization to multiobjective optimization.A type of directional derivatives for a multiobjective function is proposed,and with this notion characterizations of strict local minima of order k for a multiobjective optimization problem with a nonempty set constraint are established,generalizing the corresponding scalar case obtained by Studniarski[3].Also necessary not sufficient and sufficient not necessary optimality conditions for this minima are derived based on our directional derivatives,which are generalizations of some existing scalar results and equivalent to some existing multiobjective ones.Many examples are given to illustrate them there.展开更多
Higher-order displacement-based finite element methods are useful for simulating bending problems and potentially addressing mesh-locking associated with nearly-incompressible elasticity,yet are computationally expens...Higher-order displacement-based finite element methods are useful for simulating bending problems and potentially addressing mesh-locking associated with nearly-incompressible elasticity,yet are computationally expensive.To address the computational expense,the paper presents a matrix-free,displacement-based,higher-order,hexahedral finite element implementation of compressible and nearly-compressible(ν→0.5)linear isotropic elasticity at small strain with p-multigrid preconditioning.The cost,solve time,and scalability of the implementation with respect to strain energy error are investigated for polynomial order p=1,2,3,4 for compressible elasticity,and p=2,3,4 for nearly-incompressible elasticity,on different number of CPU cores for a tube bending problem.In the context of this matrix-free implementation,higher-order polynomials(p=3,4)generally are faster in achieving better accuracy in the solution than lower-order polynomials(p=1,2).However,for a beam bending simulation with stress concentration(singularity),it is demonstrated that higher-order finite elements do not improve the spatial order of convergence,even though accuracy is improved.展开更多
Potential field due to line sources residing on slender heterogeneities is involved in various areas,such as heat conduction,potential flow,and electrostatics.Often dipolar line sources are either prescribed or induce...Potential field due to line sources residing on slender heterogeneities is involved in various areas,such as heat conduction,potential flow,and electrostatics.Often dipolar line sources are either prescribed or induced due to close interaction with other objects.Its calculation requires a higher-order scheme to take into account the dipolar effect as well as net source effect.In the present work,we apply such a higher-order line element method to analyze the potential field with cylindrical slender heterogeneities.In a benchmark example of two parallel rods,we compare the line element solution with the boundary element solution to show the accuracy as a function in terms of rods distance.Furthermore,we use more complicated examples to demonstrate the capability of the line element technique.展开更多
Information spreading has been investigated for many years,but the mechanism of why the information explosively catches on overnight is still under debate.This explosive spreading phenomenon was usually considered dri...Information spreading has been investigated for many years,but the mechanism of why the information explosively catches on overnight is still under debate.This explosive spreading phenomenon was usually considered driven separately by social reinforcement or higher-order interactions.However,due to the limitations of empirical data and theoretical analysis,how the higher-order network structure affects the explosive information spreading under the role of social reinforcement has not been fully explored.In this work,we propose an information-spreading model by considering the social reinforcement in real and synthetic higher-order networks,describable as hypergraphs.Depending on the average group size(hyperedge cardinality)and node membership(hyperdegree),we observe two different spreading behaviors:(i)The spreading progress is not sensitive to social reinforcement,resulting in the information localized in a small part of nodes;(ii)a strong social reinforcement will promote the large-scale spread of information and induce an explosive transition.Moreover,a large average group size and membership would be beneficial to the appearance of the explosive transition.Further,we display that the heterogeneity of the node membership and group size distributions benefit the information spreading.Finally,we extend the group-based approximate master equations to verify the simulation results.Our findings may help us to comprehend the rapidly information-spreading phenomenon in modern society.展开更多
AIM:To compare the visual outcomes and corneal higherorder aberrations(HOAs)of patients with high or low myopic astigmatism after small incision lenticule extraction(SMILE).METHODS:A total of 157 eyes of 157 patients ...AIM:To compare the visual outcomes and corneal higherorder aberrations(HOAs)of patients with high or low myopic astigmatism after small incision lenticule extraction(SMILE).METHODS:A total of 157 eyes of 157 patients who underwent SMILE were included in this retrospective,nonrandomized,comparative study.All the eyes which were with the rule astigmatism were divided into high astigmatism group(HAG;astigmatism≤-2.00 D,73 eyes)and low astigmatism group(LAG;astigmatism≥-1.00 D,84 eyes).Visual and refractive examinations were performed,HOAs of the anterior surface,posterior surface,and total cornea of the eyes were evaluated preoperatively and 6mo postoperatively.RESULTS:At the postoperative 6-month follow-up,uncorrected distance visual acuity of 20/20 or better was achieved in 97%and 100%eyes in HAG and LAG respectively and 74%and 100%eyes were within-0.50 D.Vector analysis revealed no significant differences in the correction index(P=0.066),angle of error(P=0.091)or flattening index(P=0.987)between two groups.The magnitude of error was-0.37±0.31 D in HAG and-0.04±0.19 D in LAG(P<0.001).Index of success(IOS)was 0.22±0.09 in the HAG and 0.50±0.46 in the LAG(P<0.001).HOAs of most anterior,posterior and total cornea significantly increased after SMILE,especially the spherical aberration and coma.For HAG,the SMILE procedure induced significantly higher anterior,posterior and total cornea horizontal coma and total corneal total HOAs compared with LAG(P<0.001)and these surgically induced HOAs predominantly originated from the anterior surface of the cornea.CONCLUSION:SMILE surgery induces more HOAs and a mild under-correction of astigmatism in eyes with high astigmatism.The increment in HOAs after SMILE is related to preoperative astigmatism.展开更多
Higher-order modes of the neutron diffusion/transport equation can be used to study the temporal behavior of nuclear reactors and can be applied in modal analysis, transient analysis, and online monitoring of the reac...Higher-order modes of the neutron diffusion/transport equation can be used to study the temporal behavior of nuclear reactors and can be applied in modal analysis, transient analysis, and online monitoring of the reactor core. Both the deterministic method and the Monte Carlo(MC) method can be used to solve the higher-order modes. However, MC method, compared to the deterministic method, faces challenges in terms of computational efficiency and α mode calculation stability, whereas the deterministic method encounters issues arising from homogenization-related geometric and energy spectra adaptation.Based on the higher-order mode diffusion calculation code HARMONY, we developed a new higher-order mode calculation code, HARMONY2.0, which retains the functionality of computing λ and α higher-order modes from HARMONY1.0, but enhances the ability to treat complex geometries and arbitrary energy spectra using the MC-deterministic hybrid two-step strategy. In HARMONY2.0, the mesh homogenized multigroup constants were obtained using OpenMC in the first step,and higher-order modes were then calculated with the mesh homogenized core diffusion model using the implicitly restarted Arnoldi method(IRAM), which was also adopted in the HARMONY1.0 code. In addition, to improve the calculation efficiency, particularly in large higher-order modes, event-driven parallelization/domain decomposition methods are embedded in the HARMONY2.0 code to accelerate the inner iteration of λ∕α mode using OpenMP. Furthermore, the higher-order modes of complex geometric models, such as Hoogenboom and ATR reactors for λ mode and the MUSE-4 experiment facility for the prompt α mode, were computed using diffusion theory.展开更多
This paper focuses on the investigation of a hyperbolic Kirchhoff equation with nonlinear damping and higher-order dissipation terms.Initially,the existence and uniqueness of local weak solutions are rigorously establ...This paper focuses on the investigation of a hyperbolic Kirchhoff equation with nonlinear damping and higher-order dissipation terms.Initially,the existence and uniqueness of local weak solutions are rigorously established.Next,within the framework of potential well theory,the classification of solution behaviors,including blow-up and global existence,is systematically analyzed according to the relationships among the exponents of nonlinear source terms.Finally,explicit bounds for the blow-up time and decay estimates for global solutions are presented.展开更多
The higher-order Kaup-Newell equation is examined by applying the Fokas unified method on the half-line.We demonstrate that the solution can be expressed in relation to the resolution of the Riemann-Hilbert problem.Th...The higher-order Kaup-Newell equation is examined by applying the Fokas unified method on the half-line.We demonstrate that the solution can be expressed in relation to the resolution of the Riemann-Hilbert problem.The jump matrix for this problem is derived from the spectral matrix,which is calculated based on both the initial conditions and the boundary conditions.The jump matrix is explicitly dependent and expressed through the spectral functions,which are derived from the initial and boundary information,respectively.These spectral functions are interdependent and adhere to a so-called global relationship.展开更多
In this paper,the exact boundary controllability of the higher-order KdVtype equation on torus is studied.That is,given the initial and final states in the appropriate space,by adding the appropriate control function ...In this paper,the exact boundary controllability of the higher-order KdVtype equation on torus is studied.That is,given the initial and final states in the appropriate space,by adding the appropriate control function on the boundary,the solution of the system can transition from the initial state to the specified final value.Firstly,we establish the observability inequality for the higher-order KdV-type equation by Ingham inequality.Then,based on the observability inequality,Hilbert uniqueness method and a integral identity we obtain the exact boundary controllability of the higher-order KdV-type equation.展开更多
We are concerned with a Camassa-Holm type equation with higher-order nonlinearity including some integrable peakon models such as the Camassa-Holm equation,the Degasperis-Procesi equation,and the Novikov equation.We s...We are concerned with a Camassa-Holm type equation with higher-order nonlinearity including some integrable peakon models such as the Camassa-Holm equation,the Degasperis-Procesi equation,and the Novikov equation.We show that all the horizontal symmetric waves for this equation must be traveling waves.This extends the previous results for the Camassa-Holm and Novikov equations.展开更多
In recent years,there has been a surge of interest in higher-order topological phases(HOTPs)across various disciplines within the field of physics.These unique phases are characterized by their ability to harbor topol...In recent years,there has been a surge of interest in higher-order topological phases(HOTPs)across various disciplines within the field of physics.These unique phases are characterized by their ability to harbor topological protected boundary states at lower-dimensional boundaries,a distinguishing feature that sets them apart from conventional topological phases and is attributed to the higher-order bulk-boundary correspondence.Two-dimensional(2D)twisted systems offer an optimal platform for investigating HOTPs,owing to their strong controllability and experimental feasibility.Here,we provide a comprehensive overview of the latest research advancements on HOTPs in 2D twisted multilayer systems.We will mainly review the HOTPs in electronic,magnonic,acoustic,photonic and mechanical twisted systems,and finally provide a perspective of this topic.展开更多
The exploration of topological phases remains a cutting-edge research frontier,driven by their promising potential for next-generation electronic and quantum technologies.In this work,we employ first-principles calcul...The exploration of topological phases remains a cutting-edge research frontier,driven by their promising potential for next-generation electronic and quantum technologies.In this work,we employ first-principles calculations and tightbinding modeling to systematically investigate the topological properties of freestanding two-dimensional(2D)honeycomb Bi,HgTe,and Al_(2)O_(3)(0001)-supported HgTe.Remarkably,all three systems exhibit coexistence of intrinsic first-and higher-order topological insulator states,induced by spin-orbit coupling(SOC).These states manifest as topologically protected gapless edge states in one-dimensional(1D)nanoribbons and symmetry-related corner states in zero-dimensional(0D)nanoflakes.Furthermore,fractional electron charges may accumulate at the corners of armchair-edged nanoflakes.Among these materials,HgTe/Al_(2)O_(3)(0001)is particularly promising due to its experimentally feasible atomic configuration and low-energy corner states.Our findings highlight the importance of exploring higher-order topological phases in quantum spin Hall insulators and pave the way for new possibilities in device applications.展开更多
Recent engineering applications increasingly adopt smart materials,whose mechanical responses are sensitive to magnetic and electric fields.In this context,new and computationally efficient modeling strategies are ess...Recent engineering applications increasingly adopt smart materials,whose mechanical responses are sensitive to magnetic and electric fields.In this context,new and computationally efficient modeling strategies are essential to predict the multiphysic behavior of advanced structures accurately.Therefore,the manuscript presents a higher-order formulation for the static analysis of laminated anisotropic magneto-electro-elastic doubly-curved shell structures.The fundamental relations account for the full coupling between the electric field,magnetic field,and mechanical elasticity.The configuration variables are expanded along the thickness direction using a generalized formulation based on the Equivalent Layer-Wise approach.Higher-order polynomials are selected,allowing for the assessment of prescribed values of the configuration variables at the top and bottom sides of solids.In addition,an effective strategy is provided for modeling general surface distributions of mechanical pressures and electromagnetic external fluxes.The model is based on a continuum-based formulation which employs an analytical homogenization of the multifield material properties,based on Mori&Tanaka approach,of a magneto-electro-elastic composite material obtained from a piezoelectric and a piezomagnetic phase,with coupled magneto-electro-elastic effects.A semi-analytical Navier solution is applied to the fundamental equations,and an efficient post-processing equilibrium-based procedure is here used,based on the numerical assessment with the Generalized Differential Quadrature(GDQ)method,to recover the response of three-dimensional shells.The formulation is validated through various examples,investigating the multifield response of panels of different curvatures and lamination schemes.An efficient homogenization procedure,based on the Mori&Tanaka approach,is employed to obtain the three-dimensional constitutive relation of magneto-electro-elastic materials.Each model is validated against three-dimensional finite-element simulations,as developed in commercial codes.Furthermore,the full coupling effect between the electric and magnetic response is evaluated via a parametric investigation,with useful insights for design purposes of many engineering applications.The paper,thus,provides a formulation for the magneto-electro-elastic analysis of laminated structures,with a high computational efficiency,since it provides results with three-dimensional capabilities with a two-dimensional formulation.The adoption of higher-order theories,indeed,allows us to efficiently predict not only the mechanical response of the structure as happens in existing literature,but also the through-the-thickness distribution of electric and magnetic variables.A novel higher-order theory has been proposed in this work for the magneto-electro-elastic analysis of laminated shell structures with varying curvatures.This theory employs a generalized method to model the distribution of the displacement field components,electrostatic,and magneto-static potential,accounting for higher-order polynomials.The thickness functions have been defined to prescribe the arbitrary values of configuration variables at the top and bottom surfaces,even though the model is ESL-based.The fundamental governing equations have been derived in curvilinear principal coordinates,considering all coupling effects among different physical phenomena,including piezoelectric,piezomagnetic,and magneto-electric effects.A homogenization algorithm based on a Mori&Tanaka approach has been adopted to obtain the equivalent magneto-electro-mechanical properties of a two-phase transversely isotropic composite.In addition,an effective method has been adopted involving the external loads in terms of surface tractions,as well as the electric and magnetic fluxes.In the post-processing stage,a GDQ-based procedure provides the actual 3D response of a doubly-curved solid.The model has been validated through significant numerical examples,showing that the results of this semi-analytical theory align well with those obtained from 3D numerical models from commercial codes.In particular,the accuracy of the model has been verified for lamination schemes with soft layers and various curvatures under different loading conditions.Moreover,this formulation has been used to predict the effect of combined electric and magnetic loads on the mechanical response of panels with different curvatures and lamination schemes.As a consequence,this theory can be applied in engineering applications where the combined effect of electric and magnetic loads is crucial,thus facilitating their study and design.An existing limitation of this study is that the solution is that it is derived only for structures with uniform curvature,cross-ply lamination scheme,and simply supported boundary conditions.Furthermore,it requires that each lamina within the stacking sequence exhibits magneto-electro-elastic behavior.Therefore,at the present stage,it cannot be used for multifield analysis of classical composite structures with magneto-electric patches.A further enhancement of the research work could be the derivation of a solution employing a numerical technique,to overcome the limitations of the Navier method.In this way,the same theory may be adopted to predict the multifield response of structures with variable curvatures and thickness,as well as anisotropic materials and more complicated boundary conditions.Acknowledgement:The authors are grateful to the Department of Innovation Engineering of Univer-sity of Salento for the support.展开更多
Under investigation is the n-component nonlinear Schrödinger equation with higher-order effects,which describes the ultrashort pulses in the birefringent fiber.Based on the Lax pair,the eigenfunction and generali...Under investigation is the n-component nonlinear Schrödinger equation with higher-order effects,which describes the ultrashort pulses in the birefringent fiber.Based on the Lax pair,the eigenfunction and generalized Darboux transformation are derived.Next,we construct several novel higher-order localized waves and classified them into three categories:(i)higher-order rogue waves interacting with bright/antidark breathers,(ii)higher-order breather fission/fusion,(iii)higherorder breather interacting with soliton.Moreover,we explore the effects of parameters on the structure,collision process and energy distribution of localized waves and these characteristics are significantly different from previous ones.Finally,the dynamical properties of these solutions are discussed in detail.展开更多
基金Supported by the National Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金The project supported by the National Natural Science Foundation of China(10172023)
文摘A new higher-order shear deformation theory based on global-local superposition technique is developed. The theory satisfies the free surface conditions and the geometric and stress continuity conditions at interfaces. The global displacement components are of the Reddy theory and local components are of the internal first to third-order terms in each layer. A two-node beam element based on this theory is proposed. The solutions are compared with 3D-elasticity solutions. Numerical results show that present beam element has higher computational efficiency and higher accuracy.
文摘With the development of aeronautic and astronautic techniques, radiation becomes much more significant while the structure is exposed to the higher and higher temperature. Most of the current finite element software packages treat it using the net-radiation method or absorbed radiation method based on the assumption of isothermal surface with uniform radiation heat flux, which brings the conflict between the precision and the quantity of grids. Using integral method to compute the variable radiation heat flux in higher-order finite element, the precision can be improved greatly while using the same quantity of grids, because it is more consistent with the distribution of real temperature. In this paper, the integral is only processed on the same integral points as those used for solving the finite element equations, so it may be of high efficiency. In an academic testing model, the result is contrast to which get in ANSYS, proving the high precision of the method. Then an actual sandwich panel used in the thermal protection system is analyzed with the method, and the error is comparatively low to the analytical answer while the computation being of high efficiency.
基金Hong Kong Research Grants Council under the GRF(9043664).
文摘This article briefly reviews the topic of complex network synchronization,with its graph-theoretic criterion,showing that the homogeneous and symmetrical network structures are essential for optimal synchronization.Furthermore,it briefly reviews the notion of higher-order network topologies and shows their promising potential in application to evaluating the optimality of network synchronizability.
基金supported by the National Natural Science Foundation of China(Grant Nos.52271278 and 52111530137)the Natural Science Foundation of Jiangsu Province(Grant No.SBK2022020579)the Newton Advanced Fellowships by the Royal Society(Grant No.NAF\R1\180304).
文摘The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.
基金Supported by the National Natural Science Foundation of China(12061048)NSF of Jiangxi Province(20232BAB201026,20232BAB201018)。
文摘In this paper,a new technique is introduced to construct higher-order iterative methods for solving nonlinear systems.The order of convergence of some iterative methods can be improved by three at the cost of introducing only one additional evaluation of the function in each step.Furthermore,some new efficient methods with a higher-order of convergence are obtained by using only a single matrix inversion in each iteration.Analyses of convergence properties and computational efficiency of these new methods are made and testified by several numerical problems.By comparison,the new schemes are more efficient than the corresponding existing ones,particularly for large problem sizes.
文摘This paper deals with extensions of higher-order optimality conditions for scalar optimization to multiobjective optimization.A type of directional derivatives for a multiobjective function is proposed,and with this notion characterizations of strict local minima of order k for a multiobjective optimization problem with a nonempty set constraint are established,generalizing the corresponding scalar case obtained by Studniarski[3].Also necessary not sufficient and sufficient not necessary optimality conditions for this minima are derived based on our directional derivatives,which are generalizations of some existing scalar results and equivalent to some existing multiobjective ones.Many examples are given to illustrate them there.
基金The research relied on computational resources[29]provided by the University of Colorado Boulder Research Computing Group,which is supported by the National1302 CMES,2021,vol.129,no.3 Science Foundation(Awards ACI-1532235 and ACI-1532236)University of Colorado Boulder,and Colorado State University.
文摘Higher-order displacement-based finite element methods are useful for simulating bending problems and potentially addressing mesh-locking associated with nearly-incompressible elasticity,yet are computationally expensive.To address the computational expense,the paper presents a matrix-free,displacement-based,higher-order,hexahedral finite element implementation of compressible and nearly-compressible(ν→0.5)linear isotropic elasticity at small strain with p-multigrid preconditioning.The cost,solve time,and scalability of the implementation with respect to strain energy error are investigated for polynomial order p=1,2,3,4 for compressible elasticity,and p=2,3,4 for nearly-incompressible elasticity,on different number of CPU cores for a tube bending problem.In the context of this matrix-free implementation,higher-order polynomials(p=3,4)generally are faster in achieving better accuracy in the solution than lower-order polynomials(p=1,2).However,for a beam bending simulation with stress concentration(singularity),it is demonstrated that higher-order finite elements do not improve the spatial order of convergence,even though accuracy is improved.
文摘Potential field due to line sources residing on slender heterogeneities is involved in various areas,such as heat conduction,potential flow,and electrostatics.Often dipolar line sources are either prescribed or induced due to close interaction with other objects.Its calculation requires a higher-order scheme to take into account the dipolar effect as well as net source effect.In the present work,we apply such a higher-order line element method to analyze the potential field with cylindrical slender heterogeneities.In a benchmark example of two parallel rods,we compare the line element solution with the boundary element solution to show the accuracy as a function in terms of rods distance.Furthermore,we use more complicated examples to demonstrate the capability of the line element technique.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12305043 and 12165016)the Natural Science Foundation of Jiangsu Province(Grant No.BK20220511)+1 种基金the Project of Undergraduate Scientific Research(Grant No.22A684)the support from the Jiangsu Specially-Appointed Professor Program。
文摘Information spreading has been investigated for many years,but the mechanism of why the information explosively catches on overnight is still under debate.This explosive spreading phenomenon was usually considered driven separately by social reinforcement or higher-order interactions.However,due to the limitations of empirical data and theoretical analysis,how the higher-order network structure affects the explosive information spreading under the role of social reinforcement has not been fully explored.In this work,we propose an information-spreading model by considering the social reinforcement in real and synthetic higher-order networks,describable as hypergraphs.Depending on the average group size(hyperedge cardinality)and node membership(hyperdegree),we observe two different spreading behaviors:(i)The spreading progress is not sensitive to social reinforcement,resulting in the information localized in a small part of nodes;(ii)a strong social reinforcement will promote the large-scale spread of information and induce an explosive transition.Moreover,a large average group size and membership would be beneficial to the appearance of the explosive transition.Further,we display that the heterogeneity of the node membership and group size distributions benefit the information spreading.Finally,we extend the group-based approximate master equations to verify the simulation results.Our findings may help us to comprehend the rapidly information-spreading phenomenon in modern society.
文摘AIM:To compare the visual outcomes and corneal higherorder aberrations(HOAs)of patients with high or low myopic astigmatism after small incision lenticule extraction(SMILE).METHODS:A total of 157 eyes of 157 patients who underwent SMILE were included in this retrospective,nonrandomized,comparative study.All the eyes which were with the rule astigmatism were divided into high astigmatism group(HAG;astigmatism≤-2.00 D,73 eyes)and low astigmatism group(LAG;astigmatism≥-1.00 D,84 eyes).Visual and refractive examinations were performed,HOAs of the anterior surface,posterior surface,and total cornea of the eyes were evaluated preoperatively and 6mo postoperatively.RESULTS:At the postoperative 6-month follow-up,uncorrected distance visual acuity of 20/20 or better was achieved in 97%and 100%eyes in HAG and LAG respectively and 74%and 100%eyes were within-0.50 D.Vector analysis revealed no significant differences in the correction index(P=0.066),angle of error(P=0.091)or flattening index(P=0.987)between two groups.The magnitude of error was-0.37±0.31 D in HAG and-0.04±0.19 D in LAG(P<0.001).Index of success(IOS)was 0.22±0.09 in the HAG and 0.50±0.46 in the LAG(P<0.001).HOAs of most anterior,posterior and total cornea significantly increased after SMILE,especially the spherical aberration and coma.For HAG,the SMILE procedure induced significantly higher anterior,posterior and total cornea horizontal coma and total corneal total HOAs compared with LAG(P<0.001)and these surgically induced HOAs predominantly originated from the anterior surface of the cornea.CONCLUSION:SMILE surgery induces more HOAs and a mild under-correction of astigmatism in eyes with high astigmatism.The increment in HOAs after SMILE is related to preoperative astigmatism.
基金supported by the National Natural Science Foundation of China(No.U2267207)Science and Technology on Reactor System Design Technology Laboratory(No.KFKT-05-FWHTWU-2023004).
文摘Higher-order modes of the neutron diffusion/transport equation can be used to study the temporal behavior of nuclear reactors and can be applied in modal analysis, transient analysis, and online monitoring of the reactor core. Both the deterministic method and the Monte Carlo(MC) method can be used to solve the higher-order modes. However, MC method, compared to the deterministic method, faces challenges in terms of computational efficiency and α mode calculation stability, whereas the deterministic method encounters issues arising from homogenization-related geometric and energy spectra adaptation.Based on the higher-order mode diffusion calculation code HARMONY, we developed a new higher-order mode calculation code, HARMONY2.0, which retains the functionality of computing λ and α higher-order modes from HARMONY1.0, but enhances the ability to treat complex geometries and arbitrary energy spectra using the MC-deterministic hybrid two-step strategy. In HARMONY2.0, the mesh homogenized multigroup constants were obtained using OpenMC in the first step,and higher-order modes were then calculated with the mesh homogenized core diffusion model using the implicitly restarted Arnoldi method(IRAM), which was also adopted in the HARMONY1.0 code. In addition, to improve the calculation efficiency, particularly in large higher-order modes, event-driven parallelization/domain decomposition methods are embedded in the HARMONY2.0 code to accelerate the inner iteration of λ∕α mode using OpenMP. Furthermore, the higher-order modes of complex geometric models, such as Hoogenboom and ATR reactors for λ mode and the MUSE-4 experiment facility for the prompt α mode, were computed using diffusion theory.
基金Supported by Shandong Provincial Natural Science Foundation of China(Grant No.ZR2021MA003)。
文摘This paper focuses on the investigation of a hyperbolic Kirchhoff equation with nonlinear damping and higher-order dissipation terms.Initially,the existence and uniqueness of local weak solutions are rigorously established.Next,within the framework of potential well theory,the classification of solution behaviors,including blow-up and global existence,is systematically analyzed according to the relationships among the exponents of nonlinear source terms.Finally,explicit bounds for the blow-up time and decay estimates for global solutions are presented.
文摘The higher-order Kaup-Newell equation is examined by applying the Fokas unified method on the half-line.We demonstrate that the solution can be expressed in relation to the resolution of the Riemann-Hilbert problem.The jump matrix for this problem is derived from the spectral matrix,which is calculated based on both the initial conditions and the boundary conditions.The jump matrix is explicitly dependent and expressed through the spectral functions,which are derived from the initial and boundary information,respectively.These spectral functions are interdependent and adhere to a so-called global relationship.
文摘In this paper,the exact boundary controllability of the higher-order KdVtype equation on torus is studied.That is,given the initial and final states in the appropriate space,by adding the appropriate control function on the boundary,the solution of the system can transition from the initial state to the specified final value.Firstly,we establish the observability inequality for the higher-order KdV-type equation by Ingham inequality.Then,based on the observability inequality,Hilbert uniqueness method and a integral identity we obtain the exact boundary controllability of the higher-order KdV-type equation.
基金partially supported by the National Natural Science Foundation of China(Grant No.12201417)the Project funded by the China Postdoctoral Science Foundation(Grant No.2023M733173)partially supported by the National Natural Science Foundation of China(Grant No.12375006)。
文摘We are concerned with a Camassa-Holm type equation with higher-order nonlinearity including some integrable peakon models such as the Camassa-Holm equation,the Degasperis-Procesi equation,and the Novikov equation.We show that all the horizontal symmetric waves for this equation must be traveling waves.This extends the previous results for the Camassa-Holm and Novikov equations.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12304539,12074108,12474151,12347101)the Natural Science Foundation of Chongqing(Grant No.CSTB2022NSCQ-MSX0568)Beijing National Laboratory for Condensed Matter Physics(Grant No.2024BNLCMPKF025)。
文摘In recent years,there has been a surge of interest in higher-order topological phases(HOTPs)across various disciplines within the field of physics.These unique phases are characterized by their ability to harbor topological protected boundary states at lower-dimensional boundaries,a distinguishing feature that sets them apart from conventional topological phases and is attributed to the higher-order bulk-boundary correspondence.Two-dimensional(2D)twisted systems offer an optimal platform for investigating HOTPs,owing to their strong controllability and experimental feasibility.Here,we provide a comprehensive overview of the latest research advancements on HOTPs in 2D twisted multilayer systems.We will mainly review the HOTPs in electronic,magnonic,acoustic,photonic and mechanical twisted systems,and finally provide a perspective of this topic.
基金supported by the Program for Science and Technology Innovation Team in Zhejiang Province,China(Grant No.2021R01004)the Six Talent Peaks Project of Jiangsu Province,China(Grant No.2019-XCL-081)the Startup Funding of Ningbo University and Yongjiang Recruitment Project(Grant No.432200942).
文摘The exploration of topological phases remains a cutting-edge research frontier,driven by their promising potential for next-generation electronic and quantum technologies.In this work,we employ first-principles calculations and tightbinding modeling to systematically investigate the topological properties of freestanding two-dimensional(2D)honeycomb Bi,HgTe,and Al_(2)O_(3)(0001)-supported HgTe.Remarkably,all three systems exhibit coexistence of intrinsic first-and higher-order topological insulator states,induced by spin-orbit coupling(SOC).These states manifest as topologically protected gapless edge states in one-dimensional(1D)nanoribbons and symmetry-related corner states in zero-dimensional(0D)nanoflakes.Furthermore,fractional electron charges may accumulate at the corners of armchair-edged nanoflakes.Among these materials,HgTe/Al_(2)O_(3)(0001)is particularly promising due to its experimentally feasible atomic configuration and low-energy corner states.Our findings highlight the importance of exploring higher-order topological phases in quantum spin Hall insulators and pave the way for new possibilities in device applications.
基金funded by the Project PNRR M4C2—Innovation Grant DIRECT:Digital twIns foR EmergenCy supporT—CUP F83C22000740001.
文摘Recent engineering applications increasingly adopt smart materials,whose mechanical responses are sensitive to magnetic and electric fields.In this context,new and computationally efficient modeling strategies are essential to predict the multiphysic behavior of advanced structures accurately.Therefore,the manuscript presents a higher-order formulation for the static analysis of laminated anisotropic magneto-electro-elastic doubly-curved shell structures.The fundamental relations account for the full coupling between the electric field,magnetic field,and mechanical elasticity.The configuration variables are expanded along the thickness direction using a generalized formulation based on the Equivalent Layer-Wise approach.Higher-order polynomials are selected,allowing for the assessment of prescribed values of the configuration variables at the top and bottom sides of solids.In addition,an effective strategy is provided for modeling general surface distributions of mechanical pressures and electromagnetic external fluxes.The model is based on a continuum-based formulation which employs an analytical homogenization of the multifield material properties,based on Mori&Tanaka approach,of a magneto-electro-elastic composite material obtained from a piezoelectric and a piezomagnetic phase,with coupled magneto-electro-elastic effects.A semi-analytical Navier solution is applied to the fundamental equations,and an efficient post-processing equilibrium-based procedure is here used,based on the numerical assessment with the Generalized Differential Quadrature(GDQ)method,to recover the response of three-dimensional shells.The formulation is validated through various examples,investigating the multifield response of panels of different curvatures and lamination schemes.An efficient homogenization procedure,based on the Mori&Tanaka approach,is employed to obtain the three-dimensional constitutive relation of magneto-electro-elastic materials.Each model is validated against three-dimensional finite-element simulations,as developed in commercial codes.Furthermore,the full coupling effect between the electric and magnetic response is evaluated via a parametric investigation,with useful insights for design purposes of many engineering applications.The paper,thus,provides a formulation for the magneto-electro-elastic analysis of laminated structures,with a high computational efficiency,since it provides results with three-dimensional capabilities with a two-dimensional formulation.The adoption of higher-order theories,indeed,allows us to efficiently predict not only the mechanical response of the structure as happens in existing literature,but also the through-the-thickness distribution of electric and magnetic variables.A novel higher-order theory has been proposed in this work for the magneto-electro-elastic analysis of laminated shell structures with varying curvatures.This theory employs a generalized method to model the distribution of the displacement field components,electrostatic,and magneto-static potential,accounting for higher-order polynomials.The thickness functions have been defined to prescribe the arbitrary values of configuration variables at the top and bottom surfaces,even though the model is ESL-based.The fundamental governing equations have been derived in curvilinear principal coordinates,considering all coupling effects among different physical phenomena,including piezoelectric,piezomagnetic,and magneto-electric effects.A homogenization algorithm based on a Mori&Tanaka approach has been adopted to obtain the equivalent magneto-electro-mechanical properties of a two-phase transversely isotropic composite.In addition,an effective method has been adopted involving the external loads in terms of surface tractions,as well as the electric and magnetic fluxes.In the post-processing stage,a GDQ-based procedure provides the actual 3D response of a doubly-curved solid.The model has been validated through significant numerical examples,showing that the results of this semi-analytical theory align well with those obtained from 3D numerical models from commercial codes.In particular,the accuracy of the model has been verified for lamination schemes with soft layers and various curvatures under different loading conditions.Moreover,this formulation has been used to predict the effect of combined electric and magnetic loads on the mechanical response of panels with different curvatures and lamination schemes.As a consequence,this theory can be applied in engineering applications where the combined effect of electric and magnetic loads is crucial,thus facilitating their study and design.An existing limitation of this study is that the solution is that it is derived only for structures with uniform curvature,cross-ply lamination scheme,and simply supported boundary conditions.Furthermore,it requires that each lamina within the stacking sequence exhibits magneto-electro-elastic behavior.Therefore,at the present stage,it cannot be used for multifield analysis of classical composite structures with magneto-electric patches.A further enhancement of the research work could be the derivation of a solution employing a numerical technique,to overcome the limitations of the Navier method.In this way,the same theory may be adopted to predict the multifield response of structures with variable curvatures and thickness,as well as anisotropic materials and more complicated boundary conditions.Acknowledgement:The authors are grateful to the Department of Innovation Engineering of Univer-sity of Salento for the support.
基金Project supported by the National Natural Science Foundation of China(Grant No.12271096)the Natural Science Foundation of Fujian Province(Grant No.2021J01302)。
文摘Under investigation is the n-component nonlinear Schrödinger equation with higher-order effects,which describes the ultrashort pulses in the birefringent fiber.Based on the Lax pair,the eigenfunction and generalized Darboux transformation are derived.Next,we construct several novel higher-order localized waves and classified them into three categories:(i)higher-order rogue waves interacting with bright/antidark breathers,(ii)higher-order breather fission/fusion,(iii)higherorder breather interacting with soliton.Moreover,we explore the effects of parameters on the structure,collision process and energy distribution of localized waves and these characteristics are significantly different from previous ones.Finally,the dynamical properties of these solutions are discussed in detail.
