Optimizing wind energy harvesting performance remains a significant challenge.Machine learning(ML)offers a promising approach for addressing this challenge.This study proposes an ML-based approach using the radial bas...Optimizing wind energy harvesting performance remains a significant challenge.Machine learning(ML)offers a promising approach for addressing this challenge.This study proposes an ML-based approach using the radial basis function neural network(RBFNN)and differential evolution(DE)to predict and optimize the structural parameters(the diameter of the spherical bluff body D,the total spring stiffness k,and the length of the piezoelectric cantilever beam L)of the wind energy harvester(WEH).The RBFNN model is trained with theoretical data and validated with wind tunnel experimental results,achieving the coefficient-of-determination scores R2of 97.8%and 90.3%for predicting the average output power Pavgand aero-electro-mechanical efficiencyηaem,respectively.The DE algorithm is used to identify the optimal parameter combinations for wind speeds U ranging from 2.5 m/s to 6.5 m/s.The maximum Pavgis achieved when D=57.5 mm,k=28.8 N/m,L=112.1 mm,and U=4.6 m/s,while the maximumηaemis achieved when D=52.7 mm,k=29.2 N/m,L=89.2 mm,and U=4.7 m/s.Compared with that of the non-optimized structure,the WEH performance is improved by 28.6%in P_(avg)and 19.1%inη_(aem).展开更多
1 Legal requirements The author(s)guarantee(s)that the manuscript will not be published elsewhere in any language without the consent of the copyright holders,that the rights of third parties will not be violated,and ...1 Legal requirements The author(s)guarantee(s)that the manuscript will not be published elsewhere in any language without the consent of the copyright holders,that the rights of third parties will not be violated,and that the publisher will not be held legally responsible should there be any claims for compensation.展开更多
The narrow attenuation bands of traditional marine structures have long been a challenge in mitigating water waves.In this paper,a metastructure(MS)composed of floating periodic pontoons is proposed for broadband wate...The narrow attenuation bands of traditional marine structures have long been a challenge in mitigating water waves.In this paper,a metastructure(MS)composed of floating periodic pontoons is proposed for broadband water wave attenuation.The interaction of surface gravity waves with the MS is investigated using linear wave theory.The potential solutions of water waves by the MS with a finite array are developed by using the eigenfunction expansion matching method(EEMM),and the band structure of the MS is calculated by the transfer matrix method(TMM),in which the evanescent modes of waves are considered.The solution is verified against the existing numerical result for a special case.Based on the present solution,the association between Bragg resonance reflection and Bloch bandgaps is examined,the effects of pontoon geometry are analyzed,and the comparison between floating MS and bottom-mounted periodic structures is conducted.A computational fluid dynamics(CFD)model is further developed to assess the structures in practical fluid environments,and the floating MS presents excellent wave attenuation performance.The study presented here may provide a promising solution for protecting the coast and offshore structures.展开更多
In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascina...In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascinating.For this purpose,the effects of the wind speed and random excitation level are investigated with the stochastic averaging method(SAM)based on the energy envelope.The results of the analytical prediction are verified with the Monte-Carlo method(MCM).The numerical simulation shows that the introduction of wind can reduce the critical excitation level for triggering an inter-well jump and make a bi-stable energy harvester(BEH)realize the performance enhancement for a weak base excitation.However,as the strength of the wind increases to a particular level,the influence of the random base excitation on the dynamic responses is weakened,and the system exhibits a periodic galloping response.A comparison between a BEH and a linear energy harvester(LEH)indicates that the BEH demonstrates inferior performance for high-speed wind.Relevant experiments are conducted to investigate the validity of the theoretical prediction and numerical simulation.The experimental findings also show that strong random excitation is favorable for the BEH in the range of low wind speeds.However,as the speed of the incoming wind is up to a particular level,the disadvantage of the BEH becomes clear and evident.展开更多
Wetting of a liquid droplet on another liquid substrate is governed by the well-known Neumann equations.The present work aims to develop an explicit modified version of the Neumann equations for axisymmetric wetting o...Wetting of a liquid droplet on another liquid substrate is governed by the well-known Neumann equations.The present work aims to develop an explicit modified version of the Neumann equations for axisymmetric wetting of a liquid droplet on a highly stretched elastic membrane of non-zero bending rigidity.An explicit modified form of the Neumann equations is derived to determine the two contact angles,which is reduced to Young's equation for a liquid droplet on a rigid membrane or to the Neumann equations for a liquid droplet on another liquid substrate.Further implications of the modified Neumann equations are examined by comparison with some previous results reported in the recent literature,particularly considering the ranges of material and geometrical parameters of the liquid droplet-membrane system which have not been recently addressed in the literature.展开更多
Wrinkles in flat graded elastic layers have been recently described as a timevarying Hamiltonian system by the energy method.Cylindrical core/shell structures can also undergo surface instabilities under the external ...Wrinkles in flat graded elastic layers have been recently described as a timevarying Hamiltonian system by the energy method.Cylindrical core/shell structures can also undergo surface instabilities under the external pressure.In this study,we show that by treating the radial direction as a pseudo-time variable,the graded core/shell system with radially decaying elastic properties can also be described within the symplectic framework.In combination with the shell buckling equation,the present paper addresses the surface wrinkling of graded core/shell structures subjected to the uniform external pressure by solving a series of ordinary differential equations with varying coefficients.Three representative gradient distributions are showcased,and the predicted critical pressure and critical wave number are verified by finite element simulations.The symplectic framework provides an efficient and accurate approach to understand the surface instability and morphological evolution in curved biological tissues and engineered structures.展开更多
The warping may become an important factor for the precise transverse vibrations of curved beams.Thus,the first aim of this study is to specify the structural design parameters where the influence of cross-sectional w...The warping may become an important factor for the precise transverse vibrations of curved beams.Thus,the first aim of this study is to specify the structural design parameters where the influence of cross-sectional warping becomes great and the first-order shear deformation theory lacks the precision necessary.The outof-plane vibrations of the first-order shear deformation theory are compared with the warping-included vibrations as the curvature and/or thickness increase for symmetric and asymmetric transversely-functionally graded(TFG)curved beams.The second aim is to determine the influence of design parameters on the vibrations.The circular/exact elliptical beams are formed via curved mixed finite elements(MFEs)based on the exact curvature and length.The stress-free conditions are satisfied on three-dimensional(3D)constitutive equations.The variation of functionally graded(FG)material constituents is considered based on the power-law dependence.The cross-sectional warping deformations are defined over a displacement-type FE formulation.The warping-included MFEs(W-MFEs)provide satisfactory 3D structural characteristics with smaller degrees of freedom(DOFs)compared with the brick FEs.The Newmark method is used for the forced vibrations.展开更多
A nonlinear vibration isolation system is promising to provide a high-efficient broadband isolation performance.In this paper,a generalized vibration isolation system is established with nonlinear stiffness,nonlinear ...A nonlinear vibration isolation system is promising to provide a high-efficient broadband isolation performance.In this paper,a generalized vibration isolation system is established with nonlinear stiffness,nonlinear viscous damping,and Bouc-Wen(BW)hysteretic damping.An approximate analytical analysis is performed based on a harmonic balance method(HBM)and an alternating frequency/time(AFT)domain technique.To evaluate the damping effect,a generalized equivalent damping ratio is defined with the stiffness-varying characteristics.A comprehensive comparison of different kinds of damping is made through numerical simulations.It is found that the damping ratio of the linear damping is related to the stiffness-varying characteristics while the damping ratios of two kinds of nonlinear damping are related to the responding amplitudes.The linear damping,hysteretic damping,and nonlinear viscous damping are suitable for the small-amplitude,medium-amplitude,and large-amplitude conditions,respectively.The hysteretic damping has an extra advantage of broadband isolation.展开更多
The problem of cavity stability widely exists in deep underground engineering and energy exploitation.First,the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak stren...The problem of cavity stability widely exists in deep underground engineering and energy exploitation.First,the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak strength drop model considering the rock’s characteristics of constant modulus and double moduli.Then,the orthogonal non-associative flow rule is used to establish the displacement of the surrounding rock under constant modulus and double moduli,respectively,considering the stiffness degradation and dilatancy effects in the plastic region and assuming that the elastic strain in the plastic region satisfies the elastic constitutive relationship.Finally,the evolution of the displacement in the surrounding rock is analyzed under the effects of the double modulus characteristics,the strength drop,the stiffness degradation,and the dilatancy.The results show that the displacement solutions of the surrounding rock under constant modulus and double moduli have a unified expression.The coefficients of the expression are related to the stress field of the original rock,the elastic constant of the surrounding rock,the strength parameters,and the dilatancy angle.The strength drop,the stiffness degradation,and the dilatancy effects all have effects on the displacement.The effects can be characterized by quantitative relationships.展开更多
Fines migration induced by injection of low-salinity water(LSW) into porous media can lead to severe pore plugging and consequent permeability reduction. The deepbed filtration(DBF) theory is used to model the aforeme...Fines migration induced by injection of low-salinity water(LSW) into porous media can lead to severe pore plugging and consequent permeability reduction. The deepbed filtration(DBF) theory is used to model the aforementioned phenomenon, which allows us to predict the effluent concentration history and the distribution profile of entrapped particles. However, the previous models fail to consider the movement of the waterflood front. In this study, we derive a stochastic model for fines migration during LSW flooding, in which the Rankine-Hugoniot condition is used to calculate the concentration of detached particles behind and ahead of the moving water front. A downscaling procedure is developed to determine the evolution of pore-size distribution from the exact solution of a large-scale equation system. To validate the proposed model,the obtained exact solutions are used to treat the laboratory data of LSW flooding in artificial soil-packed columns. The tuning results show that the proposed model yields a considerably higher value of the coefficient of determination, compared with the previous models, indicating that the new model can successfully capture the effect of the moving water front on fines migration and precisely match the effluent history of the detached particles.展开更多
The direct numerical simulation (DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the...The direct numerical simulation (DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the turbulent drag reduction are investigated. The simulation is carried out at a friction Reynolds number of 180, which is based on the averaged friction velocity at the interface between the porous medium and the clear fluid domain. The depth of the porous layer ranges from 0.9 to 54 viscous units. The permeability in the spanwise direction is set to be lower than the other directions in the present simulation. The maximum drag reduction obtained is about 15.3% which occurs for a depth of 9 viscous units. The increasing of drag is addressed when the depth of the porous layer is more than 25 wall units. The thinner porous layer restricts the spanwise extension of the streamwise vortices which suppresses the bursting events near the wall. However, for the thicker porous layer, the wall-normal fluctuations are enhanced due to the weakening of the wall-blocking effect which can trigger strong turbulent structures near the wall.展开更多
The two-dimensional (2D) Eshelby tensors are discussed. Based upon the complex variable method, an integrity basis of ten isotropic invariants of the 2D Eshelby tensors is obtained. Since an integrity basis is always ...The two-dimensional (2D) Eshelby tensors are discussed. Based upon the complex variable method, an integrity basis of ten isotropic invariants of the 2D Eshelby tensors is obtained. Since an integrity basis is always a polynomial functional basis, these ten isotropic invariants are further proven to form an irreducible polynomial functional basis of the 2D Eshelby tensors.展开更多
For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. ...For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.展开更多
Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a...Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.展开更多
This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy pro...This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy projection （EEP） technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element （FE） model is reached. This conceptual dimension-by- dimension （D-by-D） discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional （2D） and three-dimensional （3D） problems can be obtained over the domain. Numerical examples of 3D Poisson＇s equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.展开更多
A splashing crown is commonly observed when a high-speed drop impacts a liquid film. The influence of the liquid viscosity on the crown＇s evolution is not yet clear. We review several existing theories of this proble...A splashing crown is commonly observed when a high-speed drop impacts a liquid film. The influence of the liquid viscosity on the crown＇s evolution is not yet clear. We review several existing theories of this problem, and carry out a series of numerical simulations. We find that a three-segment model can describe the crown＇s motion. In the very early stage when the crown is barely visible, the influence of viscosity is small. Later, a shallow water approach used in most existing models is applicable as long as the initial conditions are formulated properly. They depend on viscous dissipation in the intermediate period. Preliminary estimation based on a dissipation function is proposed to characterize the influence of viscosity in this problem.展开更多
Under the 3：1 internal resonance condition, the steady-state periodic response of the forced vibration of a traveling viscoelastic beam is studied. The viscoelastic behaviors of the traveling beam are described by th...Under the 3：1 internal resonance condition, the steady-state periodic response of the forced vibration of a traveling viscoelastic beam is studied. The viscoelastic behaviors of the traveling beam are described by the standard linear solid model, and the material time derivative is adopted in the viscoelastic constitutive relation. The direct multi-scale method is used to derive the relationships between the excitation frequency and the response amplitudes. For the first time, the real modal functions are employed to analytically investigate the periodic response of the axially traveling beam. The unde- termined coefficient method is used to approximately establish the real modal functions. The approximate analytical results are confirmed by the Galerkin truncation. Numerical examples are presented to highlight the effects of the viscoelastic behaviors on the steady-state periodic responses. To illustrate the effect of the internal resonance, the energy transfer between the internal resonance modes and the saturation-like phenomena in the steady-state responses is presented.展开更多
Nonlinear interactions of the two-dimensional （2D） second mode with oblique modes are studied numerically in a Mach 6.0 fiat-plate boundary layer, focusing on its selective enhancement effect on amplification of dif...Nonlinear interactions of the two-dimensional （2D） second mode with oblique modes are studied numerically in a Mach 6.0 fiat-plate boundary layer, focusing on its selective enhancement effect on amplification of different oblique waves. Evolution of oblique modes with various frequencies and spanwise wavenumbers in the presence of 2D second mode is simulated successively, using a modified parabolized stability equation （PSE） method, which is able to simulate interaction of two modes with different frequen- cies efficiently. Numerical results show that oblique modes in a broad band of frequencies and spanwise wavenumbers can be enhanced by the finite amplitude 2D second mode instability wave. The enhancement effect is accomplished by interaction of the 2D second mode, the oblique mode, and a forced mode with difference frequency. Two types of oblique modes are found to be more amplified, i.e., oblique modes with frequency close to that of the 2D second mode and low-frequency first mode oblique waves. Each of them may correspond to one type of transition routes found in transition experiments. The spanwise wavenumber of the oblique wave preferred by the nonlinear interaction is also determined by numerical simulations.展开更多
Transition prediction is of great importance for the design of long distance flying vehicles. It starts from the problem of receptivity, i.e., how external disturbances trigger instability waves in the boundary layer....Transition prediction is of great importance for the design of long distance flying vehicles. It starts from the problem of receptivity, i.e., how external disturbances trigger instability waves in the boundary layer. For super/hypersonic boundary layers, the external disturbances first interact with the shock ahead of the flying vehicles before entering the boundary layer. Since direct numerical simulation （DNS） is the only available tool for its comprehensive and detailed investigation, an important problem arises whether the numerical scheme, especially the shock-capturing method, can faithfully reproduce the interaction of the external disturbances with the shock, which is so far unknown. This paper is aimed to provide the answer. The interaction of weak disturbances with an oblique shock is investigated, which has a known theoretical solution. Numerical simulation using the shock-capturing method is conducted, and results are compared with those given by theoretical analysis, which shows that the adopted numerical method can faithfully reproduce the interaction of weak external disturbances with the shock.展开更多
1 Legal requirements The author（s） guarantee（s） that the manuscript will not be published elsewhere in any language without the consent of the copyright holders, that the rights of third parties will not be violat...1 Legal requirements The author（s） guarantee（s） that the manuscript will not be published elsewhere in any language without the consent of the copyright holders, that the rights of third parties will not be violated, and that the publisher will not be held legally responsible should there be any claims for compensation. Authors wishing to include figures or text passages that have already been published elsewhere are required to obtain permission from the copyright holder（s） and to include evidence that such permission has been granted when submitting their papers. Any material received without such evidence will be assumed to originate from the authors. Manuscripts at proof stage must be accompanied by a copyright statement.展开更多
基金Project supported by the National Key R&D Program of China(No.2021YFF0501001)the National Natural Science Foundation of China(Nos.52308315,51922046,and 52192661)+3 种基金the Research Funds of Huazhong University of Science and Technology(No.2023JCYJ014)the China Postdoctoral Science Foundation(No.2023M731206)the Research Funds of China Railway Siyuan Survey and Design Group Co.Ltd.(Nos.KY2023014S,KY2023126S,2021K085,2020K006,and 2020K172)the Autonomous Innovation Fund of Hubei Province of China(No.5003242027)。
文摘Optimizing wind energy harvesting performance remains a significant challenge.Machine learning(ML)offers a promising approach for addressing this challenge.This study proposes an ML-based approach using the radial basis function neural network(RBFNN)and differential evolution(DE)to predict and optimize the structural parameters(the diameter of the spherical bluff body D,the total spring stiffness k,and the length of the piezoelectric cantilever beam L)of the wind energy harvester(WEH).The RBFNN model is trained with theoretical data and validated with wind tunnel experimental results,achieving the coefficient-of-determination scores R2of 97.8%and 90.3%for predicting the average output power Pavgand aero-electro-mechanical efficiencyηaem,respectively.The DE algorithm is used to identify the optimal parameter combinations for wind speeds U ranging from 2.5 m/s to 6.5 m/s.The maximum Pavgis achieved when D=57.5 mm,k=28.8 N/m,L=112.1 mm,and U=4.6 m/s,while the maximumηaemis achieved when D=52.7 mm,k=29.2 N/m,L=89.2 mm,and U=4.7 m/s.Compared with that of the non-optimized structure,the WEH performance is improved by 28.6%in P_(avg)and 19.1%inη_(aem).
文摘1 Legal requirements The author(s)guarantee(s)that the manuscript will not be published elsewhere in any language without the consent of the copyright holders,that the rights of third parties will not be violated,and that the publisher will not be held legally responsible should there be any claims for compensation.
基金Project supported by the National Natural Science Foundation of China(No.12272128)the Natural Science Foundation of Jiangsu Province of China(No.BK20243019)。
文摘The narrow attenuation bands of traditional marine structures have long been a challenge in mitigating water waves.In this paper,a metastructure(MS)composed of floating periodic pontoons is proposed for broadband water wave attenuation.The interaction of surface gravity waves with the MS is investigated using linear wave theory.The potential solutions of water waves by the MS with a finite array are developed by using the eigenfunction expansion matching method(EEMM),and the band structure of the MS is calculated by the transfer matrix method(TMM),in which the evanescent modes of waves are considered.The solution is verified against the existing numerical result for a special case.Based on the present solution,the association between Bragg resonance reflection and Bloch bandgaps is examined,the effects of pontoon geometry are analyzed,and the comparison between floating MS and bottom-mounted periodic structures is conducted.A computational fluid dynamics(CFD)model is further developed to assess the structures in practical fluid environments,and the floating MS presents excellent wave attenuation performance.The study presented here may provide a promising solution for protecting the coast and offshore structures.
基金Project supported by the National Natural Science Foundation of China(Nos.12272355,1202520411902294)+1 种基金the Opening Foundation of Shanxi Provincial Key Laboratory for Advanced Manufacturing Technology of China(No.XJZZ202304)the Shanxi Provincial Graduate Innovation Project of China(No.2023KY629)。
文摘In the practical environment,it is very common for the simultaneous occurrence of base excitation and crosswind.Scavenging the combined energy of vibration and wind with a single energy harvesting structure is fascinating.For this purpose,the effects of the wind speed and random excitation level are investigated with the stochastic averaging method(SAM)based on the energy envelope.The results of the analytical prediction are verified with the Monte-Carlo method(MCM).The numerical simulation shows that the introduction of wind can reduce the critical excitation level for triggering an inter-well jump and make a bi-stable energy harvester(BEH)realize the performance enhancement for a weak base excitation.However,as the strength of the wind increases to a particular level,the influence of the random base excitation on the dynamic responses is weakened,and the system exhibits a periodic galloping response.A comparison between a BEH and a linear energy harvester(LEH)indicates that the BEH demonstrates inferior performance for high-speed wind.Relevant experiments are conducted to investigate the validity of the theoretical prediction and numerical simulation.The experimental findings also show that strong random excitation is favorable for the BEH in the range of low wind speeds.However,as the speed of the incoming wind is up to a particular level,the disadvantage of the BEH becomes clear and evident.
基金Project supported by the Natural Science&Engineering Research Council(NSERC)of Canada(No.NSERC-RGPIN204992)。
文摘Wetting of a liquid droplet on another liquid substrate is governed by the well-known Neumann equations.The present work aims to develop an explicit modified version of the Neumann equations for axisymmetric wetting of a liquid droplet on a highly stretched elastic membrane of non-zero bending rigidity.An explicit modified form of the Neumann equations is derived to determine the two contact angles,which is reduced to Young's equation for a liquid droplet on a rigid membrane or to the Neumann equations for a liquid droplet on another liquid substrate.Further implications of the modified Neumann equations are examined by comparison with some previous results reported in the recent literature,particularly considering the ranges of material and geometrical parameters of the liquid droplet-membrane system which have not been recently addressed in the literature.
基金Project supported by the National Natural Science Foundation of China(No.11972259)。
文摘Wrinkles in flat graded elastic layers have been recently described as a timevarying Hamiltonian system by the energy method.Cylindrical core/shell structures can also undergo surface instabilities under the external pressure.In this study,we show that by treating the radial direction as a pseudo-time variable,the graded core/shell system with radially decaying elastic properties can also be described within the symplectic framework.In combination with the shell buckling equation,the present paper addresses the surface wrinkling of graded core/shell structures subjected to the uniform external pressure by solving a series of ordinary differential equations with varying coefficients.Three representative gradient distributions are showcased,and the predicted critical pressure and critical wave number are verified by finite element simulations.The symplectic framework provides an efficient and accurate approach to understand the surface instability and morphological evolution in curved biological tissues and engineered structures.
基金Project supported by the Scientific and Technological Research Council of Turkey(TUBITAK)via 2209-A Programme。
文摘The warping may become an important factor for the precise transverse vibrations of curved beams.Thus,the first aim of this study is to specify the structural design parameters where the influence of cross-sectional warping becomes great and the first-order shear deformation theory lacks the precision necessary.The outof-plane vibrations of the first-order shear deformation theory are compared with the warping-included vibrations as the curvature and/or thickness increase for symmetric and asymmetric transversely-functionally graded(TFG)curved beams.The second aim is to determine the influence of design parameters on the vibrations.The circular/exact elliptical beams are formed via curved mixed finite elements(MFEs)based on the exact curvature and length.The stress-free conditions are satisfied on three-dimensional(3D)constitutive equations.The variation of functionally graded(FG)material constituents is considered based on the power-law dependence.The cross-sectional warping deformations are defined over a displacement-type FE formulation.The warping-included MFEs(W-MFEs)provide satisfactory 3D structural characteristics with smaller degrees of freedom(DOFs)compared with the brick FEs.The Newmark method is used for the forced vibrations.
基金Project supported by the National Natural Science Foundation of China(No.11902097)the China Postdoctoral Science Foundation(No.2019M661266)。
文摘A nonlinear vibration isolation system is promising to provide a high-efficient broadband isolation performance.In this paper,a generalized vibration isolation system is established with nonlinear stiffness,nonlinear viscous damping,and Bouc-Wen(BW)hysteretic damping.An approximate analytical analysis is performed based on a harmonic balance method(HBM)and an alternating frequency/time(AFT)domain technique.To evaluate the damping effect,a generalized equivalent damping ratio is defined with the stiffness-varying characteristics.A comprehensive comparison of different kinds of damping is made through numerical simulations.It is found that the damping ratio of the linear damping is related to the stiffness-varying characteristics while the damping ratios of two kinds of nonlinear damping are related to the responding amplitudes.The linear damping,hysteretic damping,and nonlinear viscous damping are suitable for the small-amplitude,medium-amplitude,and large-amplitude conditions,respectively.The hysteretic damping has an extra advantage of broadband isolation.
基金Project supported by the National Natural Science Foundation of China and Shandong Province Joint Program(No.U1806209)the National Natural Science Foundation of China(Nos.51774196 and 51774194)and Shandong University of Science and Technology(SDUST)Research Fund(No.2019TDJH101)。
文摘The problem of cavity stability widely exists in deep underground engineering and energy exploitation.First,the stress field of the surrounding rock under the uniform stress field is deduced based on a post-peak strength drop model considering the rock’s characteristics of constant modulus and double moduli.Then,the orthogonal non-associative flow rule is used to establish the displacement of the surrounding rock under constant modulus and double moduli,respectively,considering the stiffness degradation and dilatancy effects in the plastic region and assuming that the elastic strain in the plastic region satisfies the elastic constitutive relationship.Finally,the evolution of the displacement in the surrounding rock is analyzed under the effects of the double modulus characteristics,the strength drop,the stiffness degradation,and the dilatancy.The results show that the displacement solutions of the surrounding rock under constant modulus and double moduli have a unified expression.The coefficients of the expression are related to the stress field of the original rock,the elastic constant of the surrounding rock,the strength parameters,and the dilatancy angle.The strength drop,the stiffness degradation,and the dilatancy effects all have effects on the displacement.The effects can be characterized by quantitative relationships.
基金the National Natural Science Foundation of China(Nos.51804316,51734010,and U1762211)the National Science and Technology Major Project of China(No.2017ZX05009)the Science Foundation of China University of Petroleum,Beijing(No.2462017YJRC037)。
文摘Fines migration induced by injection of low-salinity water(LSW) into porous media can lead to severe pore plugging and consequent permeability reduction. The deepbed filtration(DBF) theory is used to model the aforementioned phenomenon, which allows us to predict the effluent concentration history and the distribution profile of entrapped particles. However, the previous models fail to consider the movement of the waterflood front. In this study, we derive a stochastic model for fines migration during LSW flooding, in which the Rankine-Hugoniot condition is used to calculate the concentration of detached particles behind and ahead of the moving water front. A downscaling procedure is developed to determine the evolution of pore-size distribution from the exact solution of a large-scale equation system. To validate the proposed model,the obtained exact solutions are used to treat the laboratory data of LSW flooding in artificial soil-packed columns. The tuning results show that the proposed model yields a considerably higher value of the coefficient of determination, compared with the previous models, indicating that the new model can successfully capture the effect of the moving water front on fines migration and precisely match the effluent history of the detached particles.
基金Project supported by the National Natural Science Foundation of China(Nos.11572183,91852111,and 11825204)the Program of Shanghai Municipal Education Commission(No.2019-01-07-00-09-E00018)
文摘The direct numerical simulation (DNS) is carried out for the incompressible viscous turbulent flows over an anisotropic porous wall. Effects of the anisotropic porous wall on turbulence modifications as well as on the turbulent drag reduction are investigated. The simulation is carried out at a friction Reynolds number of 180, which is based on the averaged friction velocity at the interface between the porous medium and the clear fluid domain. The depth of the porous layer ranges from 0.9 to 54 viscous units. The permeability in the spanwise direction is set to be lower than the other directions in the present simulation. The maximum drag reduction obtained is about 15.3% which occurs for a depth of 9 viscous units. The increasing of drag is addressed when the depth of the porous layer is more than 25 wall units. The thinner porous layer restricts the spanwise extension of the streamwise vortices which suppresses the bursting events near the wall. However, for the thicker porous layer, the wall-normal fluctuations are enhanced due to the weakening of the wall-blocking effect which can trigger strong turbulent structures near the wall.
基金Project supported by the National Natural Science Foundation of China(Nos.11271221,11771244,11571178,and 11771405)
文摘The two-dimensional (2D) Eshelby tensors are discussed. Based upon the complex variable method, an integrity basis of ten isotropic invariants of the 2D Eshelby tensors is obtained. Since an integrity basis is always a polynomial functional basis, these ten isotropic invariants are further proven to form an irreducible polynomial functional basis of the 2D Eshelby tensors.
基金Project supported by the National Natural Science Foundation of China(Nos.11702029,11771054,U1730118,91852207,and 11801036)the China Postdoctoral Science Foundation(No.2016M600967)
文摘For solid-fluid interaction, one of the phase-density equations in diffuse interface models is degenerated to a "0 = 0" equation when the volume fraction of a certain phase takes the value of zero or unity. This is because the conservative variables in phasedensity equations include volume fractions. The degeneracy can be avoided by adding an artificial quantity of another material into the pure phase. However, nonphysical waves,such as shear waves in fluids, are introduced by the artificial treatment. In this paper,a transport diffuse interface model, which is able to treat zero/unity volume fractions, is presented for solid-fluid interaction. In the proposed model, a new formulation for phase densities is derived, which is unrelated to volume fractions. Consequently, the new model is able to handle zero/unity volume fractions, and nonphysical waves caused by artificial volume fractions are prevented. One-dimensional and two-dimensional numerical tests demonstrate that more accurate results can be obtained by the proposed model.
基金Project supported by the National Natural Science Foundation of China(Nos.11572183 and 11272198)
文摘Interaction between turbulence and particles is investigated in a channel flow. The fluid motion is calculated using direct numerical simulation(DNS) with a lattice Boltzmann(LB) method, and particles are tracked in a Lagrangian framework through the action of force imposed by the fluid. The particle diameter is smaller than the Kolmogorov length scale, and the point force is used to represent the feedback force of particles on the turbulence. The effects of particles on the turbulence and skin friction coefficient are examined with different particle inertias and mass loadings. Inertial particles suppress intensities of the spanwise and wall-normal components of velocity, and the Reynolds shear stress. It is also found that, relative to the reference particle-free flow,the overall mean skin-friction coefficient is reduced by particles. Changes of near wall turbulent structures such as longer and more regular streamwise low-speed streaks and less ejections and sweeps are the manifestation of drag reduction.
基金supported by the National Natural Science Foundation of China(Nos.51378293 and 51078199)
文摘This paper presents a strategy for computation of super-convergent solutions of multi-dimensional problems in the finite element method （FEM） by recursive application of the one-dimensional （1D） element energy projection （EEP） technique. The main idea is to conceptually treat multi-dimensional problems as generalized 1D problems, based on which the concepts of generalized 1D FEM and its consequent EEP formulae have been developed in a unified manner. Equipped with these concepts, multi-dimensional problems can be recursively discretized in one dimension at each step, until a fully discretized standard finite element （FE） model is reached. This conceptual dimension-by- dimension （D-by-D） discretization procedure is entirely equivalent to a full FE discretization. As a reverse D-by-D recovery procedure, by using the unified EEP formulae together with proper extraction of the generalized nodal solutions, super-convergent displacements and first derivatives for two-dimensional （2D） and three-dimensional （3D） problems can be obtained over the domain. Numerical examples of 3D Poisson＇s equation and elasticity problem are given to verify the feasibility and effectiveness of the proposed strategy.
基金Project supported by the National Natural Science Foundation of China(Nos.11672310 and 11372326)the National Basic Research Program of China(No.2014CB04680202)
文摘A splashing crown is commonly observed when a high-speed drop impacts a liquid film. The influence of the liquid viscosity on the crown＇s evolution is not yet clear. We review several existing theories of this problem, and carry out a series of numerical simulations. We find that a three-segment model can describe the crown＇s motion. In the very early stage when the crown is barely visible, the influence of viscosity is small. Later, a shallow water approach used in most existing models is applicable as long as the initial conditions are formulated properly. They depend on viscous dissipation in the intermediate period. Preliminary estimation based on a dissipation function is proposed to characterize the influence of viscosity in this problem.
基金Project supported by the State Key Program of the National Natural Science Foundation of China(No.11232009)the National Natural Science Foundation of China(Nos.11372171 and 11422214)
文摘Under the 3：1 internal resonance condition, the steady-state periodic response of the forced vibration of a traveling viscoelastic beam is studied. The viscoelastic behaviors of the traveling beam are described by the standard linear solid model, and the material time derivative is adopted in the viscoelastic constitutive relation. The direct multi-scale method is used to derive the relationships between the excitation frequency and the response amplitudes. For the first time, the real modal functions are employed to analytically investigate the periodic response of the axially traveling beam. The unde- termined coefficient method is used to approximately establish the real modal functions. The approximate analytical results are confirmed by the Galerkin truncation. Numerical examples are presented to highlight the effects of the viscoelastic behaviors on the steady-state periodic responses. To illustrate the effect of the internal resonance, the energy transfer between the internal resonance modes and the saturation-like phenomena in the steady-state responses is presented.
基金Project supported by the National Natural Science Foundation of China(No.11332007)
文摘Nonlinear interactions of the two-dimensional （2D） second mode with oblique modes are studied numerically in a Mach 6.0 fiat-plate boundary layer, focusing on its selective enhancement effect on amplification of different oblique waves. Evolution of oblique modes with various frequencies and spanwise wavenumbers in the presence of 2D second mode is simulated successively, using a modified parabolized stability equation （PSE） method, which is able to simulate interaction of two modes with different frequen- cies efficiently. Numerical results show that oblique modes in a broad band of frequencies and spanwise wavenumbers can be enhanced by the finite amplitude 2D second mode instability wave. The enhancement effect is accomplished by interaction of the 2D second mode, the oblique mode, and a forced mode with difference frequency. Two types of oblique modes are found to be more amplified, i.e., oblique modes with frequency close to that of the 2D second mode and low-frequency first mode oblique waves. Each of them may correspond to one type of transition routes found in transition experiments. The spanwise wavenumber of the oblique wave preferred by the nonlinear interaction is also determined by numerical simulations.
基金supported by the National Natural Science Foundation of China(Nos.11472188 and11332007)the National Key Research and Development Program of China(No.2016YFA0401200)
文摘Transition prediction is of great importance for the design of long distance flying vehicles. It starts from the problem of receptivity, i.e., how external disturbances trigger instability waves in the boundary layer. For super/hypersonic boundary layers, the external disturbances first interact with the shock ahead of the flying vehicles before entering the boundary layer. Since direct numerical simulation （DNS） is the only available tool for its comprehensive and detailed investigation, an important problem arises whether the numerical scheme, especially the shock-capturing method, can faithfully reproduce the interaction of the external disturbances with the shock, which is so far unknown. This paper is aimed to provide the answer. The interaction of weak disturbances with an oblique shock is investigated, which has a known theoretical solution. Numerical simulation using the shock-capturing method is conducted, and results are compared with those given by theoretical analysis, which shows that the adopted numerical method can faithfully reproduce the interaction of weak external disturbances with the shock.
文摘1 Legal requirements The author（s） guarantee（s） that the manuscript will not be published elsewhere in any language without the consent of the copyright holders, that the rights of third parties will not be violated, and that the publisher will not be held legally responsible should there be any claims for compensation. Authors wishing to include figures or text passages that have already been published elsewhere are required to obtain permission from the copyright holder（s） and to include evidence that such permission has been granted when submitting their papers. Any material received without such evidence will be assumed to originate from the authors. Manuscripts at proof stage must be accompanied by a copyright statement.
