This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The me...This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The method leverages an efficient collocation technique to satisfy traction and electric charge equilibrium on the crack surface,aligning internal nodes with piezoelectric governing equations without needing integration or variational principles.It combines the strengths of the strong form collocation and finite element methods.The J-integral is derived analytically using the equivalent domain integral method,employing Green's formula and Gauss's divergence theorem to transform line integrals into area integrals for solving two-dimensional piezoelectric material problems.The accuracy of the method is validated through comparison with three typical examples,and it offers fracture prevention strategies for engineering piezoelectric structures under different electrical loading patterns.展开更多
In this paper,the convergence of the split-step theta method for stochastic differential equations is analyzed using stochastic C-stability and stochastic B-consistency.The fact that the numerical scheme,which is both...In this paper,the convergence of the split-step theta method for stochastic differential equations is analyzed using stochastic C-stability and stochastic B-consistency.The fact that the numerical scheme,which is both stochastically C-stable and stochastically B-consistent,is convergent has been proved in a previous paper.In order to analyze the convergence of the split-step theta method(θ∈[1/2,1]),the stochastic C-stability and stochastic B-consistency under the condition of global monotonicity have been researched,and the rate of convergence 1/2 has been explored in this paper.It can be seen that the convergence does not require the drift function should satisfy the linear growth condition whenθ=1/2 Furthermore,the rate of the convergence of the split-step scheme for stochastic differential equations with additive noise has been researched and found to be 1.Finally,an example is given to illustrate the convergence with the theoretical results.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a t...Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a three-impulse contingency return trajectory scheme is presented by combining the Lambert transfer and maneuver at the special point.Secondly,a calculation model of three-impulse contingency return trajectories is established.Then,fast calculation methods are proposed by adopting the high-order Taylor expansion of differential algebra in the twobody trajectory dynamics model and perturbed trajectory dynamics model.Finally,the performance of the proposed methods is verified by numerical simulation.The results indicate that the fast calculation method of two-body trajectory has higher calculation efficiency compared to the semi-analytical calculation method under a certain accuracy condition.Due to its high efficiency,the characteristics of the three-impulse contingency return trajectories under different contingency scenarios are further analyzed expeditiously.These findings can be used for the design of contingency return trajectories in future manned lunar landing missions.展开更多
This paper investigates the active traveling wave vibration control of an elastic supported rotating porous aluminium conical shell(CS)under impact loading.Piezoelectric smart materials in the form of micro fiber comp...This paper investigates the active traveling wave vibration control of an elastic supported rotating porous aluminium conical shell(CS)under impact loading.Piezoelectric smart materials in the form of micro fiber composites(MFCs)are used as actuators and sensors.To this end,a metal pore truncated CS with MFCs attached to its surface is considered.Adding artificial virtual springs at two edges of the truncated CS achieves various elastic supported boundaries by changing the spring stiffness.Based on the first-order shear deformation theory(FSDT),minimum energy principle,and artificial virtual spring technology,the theoretical formulations considering the electromechanical coupling are derived.The comparison of the natural frequency of the present results with the natural frequencies reported in previous literature evaluates the accuracy of the present approach.To study the vibration control,the integral quadrature method in conjunction with the differential quadrature approximation in the length direction is used to discretize the partial differential dynamical system to form a set of ordinary differential equations.With the aid of the velocity negative feedback method,both the time history and the input control voltage on the actuator are demonstrated to present the effects of velocity feedback gain,pore distribution type,semi-vertex angle,impact loading,and rotational angular velocity on the traveling wave vibration control.展开更多
This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeabi...This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.展开更多
The hydroelastic behavior of a moored oil storage vessel subjected to arbitrary time-dependent external loads,which include wind,waves,and currents with different incident directions,is investigated with the time-doma...The hydroelastic behavior of a moored oil storage vessel subjected to arbitrary time-dependent external loads,which include wind,waves,and currents with different incident directions,is investigated with the time-domain modal expansion method.First,the water boundary integral equations on the body surface of a quarter model,which can be obtained via the free-surface Green’s function method,are established.Then,the time-dependent elastic deflection of the moored oil storage vessel is expressed by a superposition of modal functions and corresponding modal amplitudes,and a Galerkin scheme is applied to derive the linear system of equations for the modal amplitudes.The second-order linear differential equations for modal amplitudes are solved via the fourth-order Runge−Kutta method.The present model is validated against existing frequency domain results for a truncated cylinder and a VLFS.Numerical calculations for the moored oil storage vessel are then conducted to obtain the time series of various modal amplitudes and elastic displacements of the measurement points and the corresponding spectra with different incident directions.展开更多
This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredepend...This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.展开更多
This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node...This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.展开更多
This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the reg...This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’...To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.展开更多
In this article, a modified version of the Differential Transform Method (DTM) is employed to examine soliton pulse propagation in a weakly non-local parabolic law medium and wave propagation in optical fibers. This s...In this article, a modified version of the Differential Transform Method (DTM) is employed to examine soliton pulse propagation in a weakly non-local parabolic law medium and wave propagation in optical fibers. This semi-analytic method has the advantage of overcoming the obstacle of the hardest nonlinear terms and is used to explain the origin of the bright and dark soliton solutions through the Schrödinger equation in its non-local form and the Radhakrishnan-Kundu-Laksmannan (RKL) equation. Numerical examples demonstrate the effectiveness of this method.展开更多
The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and ...The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.展开更多
The consistency of the cell has a significant impact on battery capacity,endurance,overall performance,safety,and service life extension.However,it is challenging to identify cells with high consistency and no loss of...The consistency of the cell has a significant impact on battery capacity,endurance,overall performance,safety,and service life extension.However,it is challenging to identify cells with high consistency and no loss of battery energy.This paper presents a cell screening algorithm that integrates genetic and numerical differentiation techniques.Initially,a mathematical model for battery consistency is established,and a multi-step charging strategy is proposed to satisfy the demands of fast charging technology.Subsequently,the genetic algorithm simulates biological evolution to efficiently search for superior cell combinations within a short time while evaluating capacity,voltage consistency,and charge/discharge efficiency.Finally,through experimental validation and comparative analysis with similar algorithms,our proposed method demonstrates notable advantages in terms of both search efficiency and performance.展开更多
A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of part...A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of partial differential equations defined in a 2D or 3D geological model. The new approach refers to a “strong version” against the “weak version” of the subspace spectral method based on the variational principle or Galerkin’s weighting scheme. We incorporate local nonlinear transformations and global spline interpolations in a curved coordinate system and make the discrete grid exactly matches geometry of the model so that it is achieved to convert the global domain into subdomains and apply Chebyshev points to locally sampling physical quantities and globally computing the spatial derivatives. This new approach not only remains exponential convergence of the standard spectral method in subdomains, but also yields a sparse assembled matrix when applied for the global domain simulations. We conducted 2D and 3D synthetic experiments and compared accuracies of the numerical differentiations with traditional finite difference approaches. The results show that as the points of differentiation vector are larger than five, the subdomain Chebyshev spectral method significantly improve the accuracies of the finite difference approaches.展开更多
This study focuses on determining the second-order irregular wave loads in the time domain without using the Inverse Fast Fourier Transform(IFFT).Considering the substantial displacement effects that Floating Offshore...This study focuses on determining the second-order irregular wave loads in the time domain without using the Inverse Fast Fourier Transform(IFFT).Considering the substantial displacement effects that Floating Offshore Wind Turbine(FOWT)support structures undergo when subjected to wave loads,the time-domain wave method is more suitable,while the frequency-domain method requiring IFFT cannot be used for moving bodies.Nonetheless,the computational challenges posed by the considerable computer time requirements of the time-domain wave method remain a significant obstacle.Thus,the paper incorporates various numerical schemes,including parallel computing and extrapolation of wave forces during specific time steps to improve overall efficiency.Despite the effectiveness of these schemes,the computational difficulties associated with the time-domain wave method persist.This study then proposes an innovative approach utilizing different randomnumbers in distinct segments,significantly reducing the computation of second-order wave loads.This random number interpolation ensures a smooth curve transition between two segments,emphasizingminimizing errors near the end of the first segment.Numerical analyses demonstrate substantial decreases in total computer time for FOWT structural analyses while maintaining consistent steel design results.The proposed method is uncomplicated,requiring only a simple subprogram modification in a conventional wave load computer program.展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
This study focuses on numerically investigating thermal behavior within a differentially heated cavity filled with nanofluid with and without obstacles.Numerical comparison with previous studies proves the consistency...This study focuses on numerically investigating thermal behavior within a differentially heated cavity filled with nanofluid with and without obstacles.Numerical comparison with previous studies proves the consistency and efficacy of the lattice Boltzmann method associated with a single relaxation time and its possibility of studying the nanofluid and heat transfer with high accuracy.Key parameters,including nanoparticle type and concentration,Rayleigh number,fluid basis,and obstacle position and dimension,were examined to identify optimal conditions for enhancing heat transfer quality.Principal findings indicated that increasing the Rayleigh number boosts buoyancy forces and alters vortex structure,improving the heat transfer efficiency across all nanofluid configu-rations.Moreover,nanoparticles with higher thermal conductivity,particularly Cu nanoparticles,exhibit slight improvements in heat transfer quality compared to Al2O3 nanoparticles,while higher nanoparticle concentrations generally lead to enhanced heat transfer effectiveness.Water-Cu nanofluids also demonstrate superior heat transfer performance over ethylene glycol-Cu nanofluids.Furthermore,the presence of obstacles at cavity extremities hampers overall heat transfer,whereas those positioned centrally augment heat exchange rates.This research offers valuable insights into optimizing convective heat transfer in nanofluid-filled cavities crucial for various engineering applications.展开更多
The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect ...The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.展开更多
基金Financial support of this work by the Technology Development program of China(Grant No.2022204B003)National Natural Science Foundation of China(12272083 and 12172078)the Fundamental Research Funds for the Central Universities(DUT24YJ136)is gratefully acknowledged.
文摘This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The method leverages an efficient collocation technique to satisfy traction and electric charge equilibrium on the crack surface,aligning internal nodes with piezoelectric governing equations without needing integration or variational principles.It combines the strengths of the strong form collocation and finite element methods.The J-integral is derived analytically using the equivalent domain integral method,employing Green's formula and Gauss's divergence theorem to transform line integrals into area integrals for solving two-dimensional piezoelectric material problems.The accuracy of the method is validated through comparison with three typical examples,and it offers fracture prevention strategies for engineering piezoelectric structures under different electrical loading patterns.
基金Supported by the National Natural Science Foundation of China (Grant No. 12301521)the Natural Science Foundation of Shanxi Province (Grant No. 20210302124081)。
文摘In this paper,the convergence of the split-step theta method for stochastic differential equations is analyzed using stochastic C-stability and stochastic B-consistency.The fact that the numerical scheme,which is both stochastically C-stable and stochastically B-consistent,is convergent has been proved in a previous paper.In order to analyze the convergence of the split-step theta method(θ∈[1/2,1]),the stochastic C-stability and stochastic B-consistency under the condition of global monotonicity have been researched,and the rate of convergence 1/2 has been explored in this paper.It can be seen that the convergence does not require the drift function should satisfy the linear growth condition whenθ=1/2 Furthermore,the rate of the convergence of the split-step scheme for stochastic differential equations with additive noise has been researched and found to be 1.Finally,an example is given to illustrate the convergence with the theoretical results.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金co-supported by the National Natural Science Foundation of China(No.12072365)the Technology Innovation Team of Manned Space Engineering,China。
文摘Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a three-impulse contingency return trajectory scheme is presented by combining the Lambert transfer and maneuver at the special point.Secondly,a calculation model of three-impulse contingency return trajectories is established.Then,fast calculation methods are proposed by adopting the high-order Taylor expansion of differential algebra in the twobody trajectory dynamics model and perturbed trajectory dynamics model.Finally,the performance of the proposed methods is verified by numerical simulation.The results indicate that the fast calculation method of two-body trajectory has higher calculation efficiency compared to the semi-analytical calculation method under a certain accuracy condition.Due to its high efficiency,the characteristics of the three-impulse contingency return trajectories under different contingency scenarios are further analyzed expeditiously.These findings can be used for the design of contingency return trajectories in future manned lunar landing missions.
基金Supported by the National Natural Science Foundation of China(Nos.12272056 and 11832002)。
文摘This paper investigates the active traveling wave vibration control of an elastic supported rotating porous aluminium conical shell(CS)under impact loading.Piezoelectric smart materials in the form of micro fiber composites(MFCs)are used as actuators and sensors.To this end,a metal pore truncated CS with MFCs attached to its surface is considered.Adding artificial virtual springs at two edges of the truncated CS achieves various elastic supported boundaries by changing the spring stiffness.Based on the first-order shear deformation theory(FSDT),minimum energy principle,and artificial virtual spring technology,the theoretical formulations considering the electromechanical coupling are derived.The comparison of the natural frequency of the present results with the natural frequencies reported in previous literature evaluates the accuracy of the present approach.To study the vibration control,the integral quadrature method in conjunction with the differential quadrature approximation in the length direction is used to discretize the partial differential dynamical system to form a set of ordinary differential equations.With the aid of the velocity negative feedback method,both the time history and the input control voltage on the actuator are demonstrated to present the effects of velocity feedback gain,pore distribution type,semi-vertex angle,impact loading,and rotational angular velocity on the traveling wave vibration control.
文摘This study presents a method for the inverse analysis of fluid flow problems.The focus is put on accurately determining boundary conditions and characterizing the physical properties of granular media,such as permeability,and fluid components,like viscosity.The primary aim is to deduce either constant pressure head or pressure profiles,given the known velocity field at a steady-state flow through a conduit containing obstacles,including walls,spheres,and grains.The lattice Boltzmann method(LBM)combined with automatic differentiation(AD)(AD-LBM)is employed,with the help of the GPU-capable Taichi programming language.A lightweight tape is used to generate gradients for the entire LBM simulation,enabling end-to-end backpropagation.Our AD-LBM approach accurately estimates the boundary conditions for complex flow paths in porous media,leading to observed steady-state velocity fields and deriving macro-scale permeability and fluid viscosity.The method demonstrates significant advantages in terms of prediction accuracy and computational efficiency,making it a powerful tool for solving inverse fluid flow problems in various applications.
基金financially supported by the Department of Natural Resources of Guangdong Province(Grant No.[2024]31)the National Natural Science Foundation of China(Grant No.52071145)+1 种基金the Natural Science Foundation of Guangdong Province,China(Grant No.2022B1515020071)the Fundamental Research Funds for the Central Universities(Grant No.2023ZYGXZR029).
文摘The hydroelastic behavior of a moored oil storage vessel subjected to arbitrary time-dependent external loads,which include wind,waves,and currents with different incident directions,is investigated with the time-domain modal expansion method.First,the water boundary integral equations on the body surface of a quarter model,which can be obtained via the free-surface Green’s function method,are established.Then,the time-dependent elastic deflection of the moored oil storage vessel is expressed by a superposition of modal functions and corresponding modal amplitudes,and a Galerkin scheme is applied to derive the linear system of equations for the modal amplitudes.The second-order linear differential equations for modal amplitudes are solved via the fourth-order Runge−Kutta method.The present model is validated against existing frequency domain results for a truncated cylinder and a VLFS.Numerical calculations for the moored oil storage vessel are then conducted to obtain the time series of various modal amplitudes and elastic displacements of the measurement points and the corresponding spectra with different incident directions.
基金supported by a grant from the National Science and Technology Council of the Republic of China(Grant Number:MOST 112-2221-E-006-048-MY2).
文摘This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.
基金Anhui Provincial Natural Science Foundation(2308085QD124)Anhui Province University Natural Science Research Project(GrantNo.2023AH050918)The University Outstanding Youth Talent Support Program of Anhui Province.
文摘This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.
基金supported by the National Natural Science Foundation of China (Grant No. 41475021)
文摘This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
文摘To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.
文摘In this article, a modified version of the Differential Transform Method (DTM) is employed to examine soliton pulse propagation in a weakly non-local parabolic law medium and wave propagation in optical fibers. This semi-analytic method has the advantage of overcoming the obstacle of the hardest nonlinear terms and is used to explain the origin of the bright and dark soliton solutions through the Schrödinger equation in its non-local form and the Radhakrishnan-Kundu-Laksmannan (RKL) equation. Numerical examples demonstrate the effectiveness of this method.
文摘The numerical approach for finding the solution of fractional order systems of boundary value problems (BPVs) is derived in this paper. The implementation of the weighted residuals such as Galerkin, Least Square, and Collocation methods are included for solving fractional order differential equations, which is broadened to acquire the approximate solutions of fractional order systems with differentiable polynomials, namely Legendre polynomials, as basis functions. The algorithm of the residual formulations of matrix form can be coded efficiently. The interpretation of Caputo fractional derivatives is employed here. We have demonstrated these methods numerically through a few examples of linear and nonlinear BVPs. The results in absolute errors show that the present method efficiently finds the numerical solutions of fractional order systems of differential equations.
文摘The consistency of the cell has a significant impact on battery capacity,endurance,overall performance,safety,and service life extension.However,it is challenging to identify cells with high consistency and no loss of battery energy.This paper presents a cell screening algorithm that integrates genetic and numerical differentiation techniques.Initially,a mathematical model for battery consistency is established,and a multi-step charging strategy is proposed to satisfy the demands of fast charging technology.Subsequently,the genetic algorithm simulates biological evolution to efficiently search for superior cell combinations within a short time while evaluating capacity,voltage consistency,and charge/discharge efficiency.Finally,through experimental validation and comparative analysis with similar algorithms,our proposed method demonstrates notable advantages in terms of both search efficiency and performance.
文摘A new numerical approach, called the “subdomain Chebyshev spectral method” is presented for calculation of the spatial derivatives in a curved coordinate system, which may be employed for numerical solutions of partial differential equations defined in a 2D or 3D geological model. The new approach refers to a “strong version” against the “weak version” of the subspace spectral method based on the variational principle or Galerkin’s weighting scheme. We incorporate local nonlinear transformations and global spline interpolations in a curved coordinate system and make the discrete grid exactly matches geometry of the model so that it is achieved to convert the global domain into subdomains and apply Chebyshev points to locally sampling physical quantities and globally computing the spatial derivatives. This new approach not only remains exponential convergence of the standard spectral method in subdomains, but also yields a sparse assembled matrix when applied for the global domain simulations. We conducted 2D and 3D synthetic experiments and compared accuracies of the numerical differentiations with traditional finite difference approaches. The results show that as the points of differentiation vector are larger than five, the subdomain Chebyshev spectral method significantly improve the accuracies of the finite difference approaches.
基金funded by National Science and Technology Council,grant number NSTC 113-2223-E-006-014.
文摘This study focuses on determining the second-order irregular wave loads in the time domain without using the Inverse Fast Fourier Transform(IFFT).Considering the substantial displacement effects that Floating Offshore Wind Turbine(FOWT)support structures undergo when subjected to wave loads,the time-domain wave method is more suitable,while the frequency-domain method requiring IFFT cannot be used for moving bodies.Nonetheless,the computational challenges posed by the considerable computer time requirements of the time-domain wave method remain a significant obstacle.Thus,the paper incorporates various numerical schemes,including parallel computing and extrapolation of wave forces during specific time steps to improve overall efficiency.Despite the effectiveness of these schemes,the computational difficulties associated with the time-domain wave method persist.This study then proposes an innovative approach utilizing different randomnumbers in distinct segments,significantly reducing the computation of second-order wave loads.This random number interpolation ensures a smooth curve transition between two segments,emphasizingminimizing errors near the end of the first segment.Numerical analyses demonstrate substantial decreases in total computer time for FOWT structural analyses while maintaining consistent steel design results.The proposed method is uncomplicated,requiring only a simple subprogram modification in a conventional wave load computer program.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.
文摘This study focuses on numerically investigating thermal behavior within a differentially heated cavity filled with nanofluid with and without obstacles.Numerical comparison with previous studies proves the consistency and efficacy of the lattice Boltzmann method associated with a single relaxation time and its possibility of studying the nanofluid and heat transfer with high accuracy.Key parameters,including nanoparticle type and concentration,Rayleigh number,fluid basis,and obstacle position and dimension,were examined to identify optimal conditions for enhancing heat transfer quality.Principal findings indicated that increasing the Rayleigh number boosts buoyancy forces and alters vortex structure,improving the heat transfer efficiency across all nanofluid configu-rations.Moreover,nanoparticles with higher thermal conductivity,particularly Cu nanoparticles,exhibit slight improvements in heat transfer quality compared to Al2O3 nanoparticles,while higher nanoparticle concentrations generally lead to enhanced heat transfer effectiveness.Water-Cu nanofluids also demonstrate superior heat transfer performance over ethylene glycol-Cu nanofluids.Furthermore,the presence of obstacles at cavity extremities hampers overall heat transfer,whereas those positioned centrally augment heat exchange rates.This research offers valuable insights into optimizing convective heat transfer in nanofluid-filled cavities crucial for various engineering applications.
文摘The static behavior of piezoelectric circular spherical shallow shells under both electrical and mechanical loads is studied by using the differential quadrature element method (DQEM). Geometrical nonlinearity effect is considered. Detailed formulations and procedures are given for the first time. Several examples are analyzed and accurate results are obtained by the DQEM. Based on the results in this paper, one may conclude that the DQEM is a useful tool for obtaining solutions of structural elements. It can be seen that the shell shape may be theore tically controlled and snap through may occur when the applied voltage reaches a critical value even without mechanical load for certain geometric configurations.
